Next Article in Journal
Artificial Intelligence to Get Insights of Multi-Drug Resistance Risk Factors during the First 48 Hours from ICU Admission
Next Article in Special Issue
The Food Anti-Microbials β-Phenylethylamine (-HCl) and Ethyl Acetoacetate Do Not Change during the Heating Process
Previous Article in Journal
Identification of Escherichia coli and Related Enterobacteriaceae and Examination of Their Phenotypic Antimicrobial Resistance Patterns: A Pilot Study at A Wildlife–Livestock Interface in Lusaka, Zambia
Previous Article in Special Issue
Elicitation of Stress-Induced Phenolic Metabolites for Antimicrobial Applications against Foodborne Human Bacterial Pathogens
 
 
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 3
10.3390/antibiotics10030237
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

Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape

by
Yvonne Ma
1,†,
Jessica Chen
1,†,
Karen Fong
1,
Stephanie Nadya
1,
Kevin Allen
1,
Chad Laing
2,
Kim Ziebell
3,
Ed Topp
4,
Laura M. Carroll
5,
Martin Wiedmann
5,
Pascal Delaquis
6 and
Siyun Wang
1,*
1
Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
2
National Centre for Animal Diseases, Canadian Food Inspection Agency, Lethbridge, AB T1J 3Z4, Canada
3
Laboratory for Foodborne Zoonoses, Public Health Agency of Canada, Guelph, ON N1G 3W4, Canada
4
London Research and Development Centre, Agriculture and Agri-Food Canada, London, ON N5V 4T3, Canada
5
Department of Food Science, Cornell University, Ithaca, NY 14853, USA
6
Summerland Research and Development Centre, Agriculture and Agri-Food Canada, Summerland, BC V0H 1Z0, Canada
*
Author to whom correspondence should be addressed.
†
Authors contributed equally to the work.
Antibiotics 2021, 10(3), 237; https://doi.org/10.3390/antibiotics10030237
Submission received: 21 January 2021 / Revised: 16 February 2021 / Accepted: 19 February 2021 / Published: 26 February 2021
(This article belongs to the Special Issue Development of Novel Anti-microbials to Reduce Bacterial Contamination of Food)
Browse Figures
Versions Notes

Abstract

:
Antibiotic resistance (AR) phenotypes and acquired resistance determinants (ARDs) detected by in silico analysis of genome sequences were examined in 55 Shiga toxin-producing Escherichia coli (STEC) isolates representing diverse serotypes recovered from surfaces waters and sediments in a mixed use urban/agricultural landscape in British Columbia, Canada. The isolates displayed decreased susceptibility to florfenicol (65.5%), chloramphenicol (7.3%), tetracycline (52.7%), ampicillin (49.1%), streptomycin (34.5%), kanamycin (20.0%), gentamycin (10.9%), amikacin (1.8%), amoxicillin/clavulanic acid (21.8%), ceftiofur (18.2%), ceftriaxone (3.6%), trimethoprim-sulfamethoxazole (12.7%), and cefoxitin (3.6%). All surface water and sediment isolates were susceptible to ciprofloxacin, nalidixic acid, ertapenem, imipenem and meropenem. Eight isolates (14.6%) were multidrug resistant. ARDs conferring resistance to phenicols (floR), trimethoprim (dfrA), sulfonamides (sul1/2), tetracyclines (tetA/B), and aminoglycosides (aadA and aph) were detected. Additionally, narrow-spectrum β-lactamase blaTEM-1b and extended-spectrum AmpC β-lactamase (cephalosporinase) blaCMY-2 were detected in the genomes, as were replicons from plasmid incompatibility groups IncFII, IncB/O/K/Z, IncQ1, IncX1, IncY and Col156. A comparison with surveillance data revealed that AR phenotypes and ARDs were comparable to those reported in generic E. coli from food animals. Aquatic environments in the region are potential reservoirs for the maintenance and transmission of antibiotic resistant STEC, associated ARDs and their plasmids.

1. Introduction

The Shiga toxigenic Escherichia coli (STEC) pathotype is characterized by the production of one or more cytotoxic Shiga-like toxins, small proteins encoded by a temperate lambdoid phage on the main chromosome. Infection with STEC causes diseases of variable severity, ranging from comparatively mild to acute hemorrhagic diarrhea, and potentially fatal thrombotic thrombocytopenic purpura or hemolytic uremic syndrome (HUS). Foodborne outbreaks have been reported with increasing frequency since two epidemiologically significant disease clusters were definitively associated with the consumption of red meat in the US in 1983 [1], and STEC are now recognized as significant water and foodborne threats to global public health [2,3,4]. A recent food-borne source attribution analysis performed using disease and surveillance data from several countries showed that foods from a broad range of categories are implicated in outbreaks [5]. Moreover, the analysis confirmed that plant-based foods, primarily fresh fruits and vegetables, are now significant vehicles for the foodborne transmission of STEC, notably in the Americas.
The advisability of antibiotic therapies for the management of STEC infections is the subject of enduring debate due to reports of antibiotic-induced toxin release from dead cells or by the induction of Shiga toxin-converting prophage [6]. For example, the fluoroquinolone antibiotic ciprofloxacin was shown to cause the induction of bacteriophage and Shiga toxin production in both in vivo and in vitro studies in a mouse model [7]. However, definitive evidence of causative association with clinical symptoms of human STEC disease remains elusive. A review of available data from clinical studies led Agger et al. [8] to conclude that it remains unclear whether the administration of protein and cell wall synthesis inhibitors influences the duration of diarrhea or the development of HUS. Moreover, a mixture of ciprofloxacin and the carbapenem meropenem was successfully used to alleviate the symptoms of HUS in patients infected with E. coli serotype O104:H4 during a major foodborne outbreak in Germany [9]. Nevertheless, antibiotic therapies are presently contraindicated by clinicians and public health authorities in some jurisdictions, primarily in the developed world. There is little doubt, however, that antibiotics are often administered to patients infected with STEC where diagnosis is delayed or masked by infections caused by other agents of gastrointestinal disease. Moreover, their use persists in many countries, including Canada, for the control and treatment of livestock diseases caused by STEC, notably edema disease in weaned piglets and watery diarrhea in calves [10].
The persistent therapeutic use of antibiotics in human and veterinary medicine or for growth promotion in animal production has contributed to the emergence and proliferation of antibiotic resistant bacteria (ARB), including zoonotic pathogenic species [11]. The global emergence of acquired resistance and concomitant decline in efficacy are anticipated to lessen the therapeutic value of present-day antibiotics used in human and veterinary medicine, with potentially dire social and economic consequences [12]. A report by the World Health Organization highlights the importance of resistance in E. coli due to the prevalence and wide spectrum of infections caused by variable pathotypes of the species [13]. Early clinical human STEC stains (primarily of serotype O157:H7) were generally susceptible to antibiotics, but there is evidence that antibiotic resistance (AR) is increasingly common in the pathotype [14,15,16,17,18,19]. Some isolates are reported to be multidrug resistant (MDR) according to the proposed standardized terminology for acquired resistance in Magiorakos et al. [20], whereby so-designated bacteria are not susceptible to more than one agent in more than three antimicrobial classes specified for the family Enterobacteriaceae. E. coli share many enzymatic and non-enzymatic resistance mechanisms with other species of Gram-negative bacteria, including intrinsic mechanisms involving inactivating enzymes or efflux pumps or mechanisms acquired by the horizontal transfer of antibiotic resistance determinants (ARDs) on mobile genetic elements leading, to the impaired absorption, modification of the antibiotic target by mutation or post-translational modification, or inactivation by hydrolysis or other chemical means [21]. Horizontal transfer is favoured by integrons, site-specific recombination systems consisting of two conserved segments flanking a central region containing “cassettes” of functional genes [22]. Integrons of relevance in the development of antibiotic resistance in enteric bacterial pathogens are primarily borne on self-replicating plasmids [23]. Most of the 27 recognized main incompatibility (Inc) groups of plasmids in Enterobacteriaceae have been associated with resistance to antibiotics in E. coli, including STEC [24,25,26].
Nucleic acid sequencing and in silico analysis are widely used to characterize ARB and associated ARDs in species recovered from disparate environments [27,28]. The current knowledge about AR and ARDs in STEC is largely based on the analysis of human clinical specimens or isolates obtained through surveillance systems that primarily target livestock species and meat products. In contrast, little is known about resident or transient STEC populations in natural or anthropogenically impacted terrestrial and aquatic ecosystems. The latter can serve as secondary reservoirs for transmission to humans by indirect routes, notably the consumption of food crops grown in environments susceptible to pollution from urban or agricultural sources [29]. The assessment of potential exposure risks and the development of mitigation measures require knowledge about the abundance and distribution of ARB and ARDs in such environments. In the present work, we report on AR phenotypes and associated ARDs in STEC recovered from surface waters in the Fraser Valley of British Columbia, Canada, an urban–agricultural region that supports intensive livestock farming and food crop production under irrigation.

2. Results

2.1. Antibiotic Resistance Phenotypes of STEC Isolates Recovered from Surface Water and Sediments

The proportion of 55 STEC isolates resistant, of intermediate resistance or susceptible to eighteen antibiotics is shown in Figure 1 and individual resistance phenotypes are provided in Table 1. All isolates were susceptible to the quinolone nalidixic acid (NAL), the fluoroquinolone ciprofloxacin (CIP) and carbapenems imipenem (IPM), ertapenem (ETP) and meropenem (MEM). The percentage of isolates resistant or of intermediate resistance to the other classes/groups of antibiotics examined were (in decreasing order of abundance): phenicols florfenicol (FFC), 65.5% and chloramphenicol (CHL), 7.3%; tetracycline (TET), 52.7%; penicillin ampicillin (AMP), 49.1%; aminoglycosides streptomycin (STR) 34.5%, kanamycin (BCN) 20.0%, gentamicin (GEN) 10.9%, and amikacin (AMK) 1.8%; penicillin + β-lactam inhibitor amoxicillin/clavulanic acid (AMC), 21.8%; third generation cephalosporins ceftiofur (TIO), 18.2% and ceftriaxone (CRO), 3.6%; folate pathway inhibitor trimethoprim-sulfamethoxazole (SXT), 12.7%; and the cephamycin cefoxitin (FOX), 3.6%. Resistance at concentrations examined in this work was limited to 10 of the antibiotics including FFC (16.4% of isolates), STR (23.6%), TET (18.2%), AMP (20.0%), SXT (12.7%), BCN (7.3%), AMC (7.3%), CHL (3.6%), TIO (3.6%) and FOX (1.8%). According to the proposed standardized terminology for acquired resistance profiles specified for the family Enterobacteriaceae in Magiorakos et al. [20], 10 isolates (18.2%) were pansusceptible, 37 (67.2%) were of intermediate resistance and eight (14.6%) were MDR. Four of the MDR isolates were from serogroup O111 and the most complex phenotype was found in isolate 385-O111:NM, which was resistant to nine antibiotics (AMPR AMCR CHLR FFCR FOXR STRR SXTR TETR TIOR) and showed intermediate resistance to two additional aminoglycosides (CROI GENI) and one third generation cephalosporin (BCNI).

2.2. Genomic Analysis of the STEC Isolates

A phylogeny inferred from the analysis of whole genome sequences was indicative of close relationships between isolates within analogous serogroups, as illustrated by the clustering of all O26:H11, O111:H8/NM and O157:H7 isolates from water/sediment and clinical sources (Figure 1). Exceptions were noted for some pairs (O109:H5/O117:H7, O5:NM/O165:H25, O116:H25/O188:H25, O91:H21/O76:H19) that clustered despite differences in serological assignment and are occasionally isolated from human, animal or environmental sources, but reports of human infections caused by strains from these serogroups are infrequent [30,31]. Few isolates have been studied in detail and little is known about their phylogenetic relationship to other serogroups. The sequences of some O gene clusters are known to differ by only a few nucleotides [32]. Hence, clustering of the isolates with different serovars by whole genome SNPs may reflect a true close evolutionary relatedness of isolates representing these different serovars.
A presence–absence matrix displayed next to the tree in Figure 2 illustrates the array of ARDs detected in the genomes of water/sediment and clinical STEC isolates; the prevalence and distribution of specific genes are provided in Table 2. ARDs conferring resistance to phenicols (catA, floR) and fosfomycin (fosA) were found infrequently. ARDs that encode resistance to sulfonamides (sul1/2), trimethoprim (dfrA) tetracyclines (tetA/B) and aminoglycosides (aac, aadA, and aph) were more common, although the aac aminoglycoside-3-acetyltransferase gene was only detected in a single isolate from water. Narrow-spectrum β-lactamase blaTEM-1b and extended-spectrum AmpC β-lactamase (cephalosporinase) blaCMY-2 -encoding gene sequences were also present in isolates from both sources, while the extended-spectrum β-lactamase-encoding blaCTX-M and the chromosomal gyrA83 mutation associated with quinolone/fluoroquinolone resistance were restricted to one and three clinical isolates, respectively. In addition to these acquired mechanisms of β-lactam resistance, three environmental isolates displayed mutations in the ampC promoter region. The chloramphenicol resistance effector catA was detected in one isolate from water. While macrolide resistance was not measured in the present study, the presence of the mphA ARD in two environmental isolates was noted in light of reports that suggest the prevalence of the gene is increasing in E. coli [33]. Overall, categorical analysis using Fisher’s exact and Bonferroni test indicated that there were no significant differences (p > 0.05) in the frequency of ARDs in isolates from either source (Table 2).
The molecular basis of resistance was examined in more detail in the eight MDR isolates. ARDs with a known association to resistance to specific antibiotics are displayed adjacent to the resistance profiles in Table 3. Concordance was observed between phenotype determined in vitro and resistance predicted by in silico analysis in all isolates except 373-O165:H25, which was resistant to FFC despite our inability to detect known phenicol ARDs. The intermediate resistance in most isolates could be attributed to the presence of ARDs that impart reduced susceptibility to specific classes of antibiotics. However, none were detected in isolates with intermediate resistance to phenicols, in three isolates with intermediate resistance to the penicillins or penicillin + β-lactam inhibitors AMP and AMC, and one with intermediate resistance to the cephalosporin TIO. The plasmid replication sequences detected in each MDR isolate are also shown in Table 3. Replicons homologous to sequences in Inc incompatibility groups were detected in all the isolates: IncFII(pRSB107) was detected in three isolates, IncFIB replicons were detected in three isolates, IncFII(pHN7A8), IncFII, and IncQ1 were found in two isolates, and lastly IncX1 and IncY were each detected in one isolate. Replicons homologous to IncB/O/K/Z were common and were detected in six of eight isolates. ColE1 were identified in three isolates and Col156, a small colicin-encoding ColE-like plasmid which is often referred to as ColE(pIGJC156), was detected in two MDR isolates.

3. Discussion

Data on the prevalence of antibiotic resistance in STEC from either clinical or environmental sources in the province of British Columbia is scarce. Allen et al. [37] reported that 10 E. coli O157:H7 and 15 non- serotype O157 human clinical isolates were susceptible to the second-generation fluoroquinolone ciprofloxacin (CIP) and the carbapenem imipenem (IPM), antimicrobial drugs that are presently listed as Category I based on their very high importance in human medicine [38]. Maal-Bared et al. examined the resistance to six antibiotics in 27 E. coli O157:H7 isolates recovered from water, sediment and biofilms in a British Columbia watershed (Elk Creek) heavily impacted by agriculture and found that all were susceptible to CIP [39]. In the present study, 55 STEC isolates representing different serotypes and recovered from surface water and sediments in four watersheds located in the same region were also susceptible to CIP and the carbapenems IPM, ertapenem (ETP) and meropenem (MEM). Directed studies performed in different countries [40,41,42,43,44,45] and passive surveillance of laboratory-confirmed O157 STEC in the US [46] indicate that resistance to fluoroquinolones and carbapenems is historically rare in STEC from clinical, animal or food sources. The provincial public health authority (British Columbia Centre for Disease Control, BCCDC) compiles data and analyzes antibiotic resistance trends in clinical human E. coli isolates of all pathotypes captured by surveillance programs. A recent report was examined to assess the scope of resistance to Category I antibiotics in strains that cause human infections for comparison with rates measured in the environmental STEC isolates. In 2014, 23.8% of clinical E. coli isolates were resistant to CIP [47]. CIP is the third most widely prescribed antibiotic in Canada, primarily as a second-line drug for the treatment of persistent human infections (principally of the urinary tract). The response to drug pressure in community and hospital settings is likely contributing to increased resistance to this antibiotic in some human E. coli pathotypes [48,49]. While CIP is no longer licensed for veterinary purposes, use of other fluoroquinolones (eg., danofloxacin, enrofloxacin) persists, primarily for the control of respiratory diseases in cattle and pigs [50]. However, the susceptibility of environmental isolates suggests that the on-farm use of fluoroquinolones is not presently a significant driver of resistance in STEC. The incidence of carbapenamase-producing clinical strains also remains very low in the province [47]. Hence, the prudent use of “last-line agent” carbapenems in human medicine and their prohibition for the treatment, disease prevention and/or growth promotion of farm animals appears to be forestalling resistance to these antibiotics in E. coli from all pathotypes, including STEC.
The susceptibility of environmental STEC to fluoroquinolones and carbapenems was reassuring in light of alarming rates of resistance reported in some jurisdictions. However, resistance or intermediate resistance to other Category I (penicillin-β-lactamase inhibitor combinations, third generation cephalosporins), high importance Category II (aminoglycosides, penicillins, cephamycins, folate pathway inhibitors) and to several classes of Category III antibiotics was recurrent in comparatively small collections of STEC isolates from the region. The comparison with provincial surveillance data indicated that rates of resistance to penicillins (AMP, 49.8%, 49.1%), aminoglycosides (GEN, 9.1%, 10.9%), folate pathway inhibitors (SXT, 24.0%, 12.7%) and third generation cephalosporins (cefotaxime, 11.2%, CRO, 3.6%) are similar in human clinical E. coli and environmental STEC isolates [47]. Provincial surveillance does not presently provide information about resistance to tetracyclines or phenicols in clinical E. coli. Both are administered to food-producing animals in Canada to prevent disease and for growth promotion. The Canadian Antimicrobial Resistance Surveillance System (CARSS) captures data on resistance to a range of antibiotics, including tetracyclines and phenicols, in generic E. coli isolated during production, slaughter and in retail meats derived from poultry, swine and cattle [51]. Because the surface waters sampled in this work are situated in a region of intensive livestock production, resistance rates in the environmental STEC were also examined against surveillance data derived from analysis of E. coli from meat animals and their products. Data from 2014 indicate that 55.3% of generic E. coli were resistant to TET, 32.2% to AMP, 29.7% to STR, 12.0% to STX and 8.1% to CRO [51]. Reciprocal resistance patterns in clinical [37] and environmental STEC described in previous work [39] and the present study suggests that common selection pressures, modes of transmission and mechanisms are likely driving resistance to antibiotics from these classes in the species. Unfortunately, none of the active surveillance systems capture data on phenicol resistance. CHL is prohibited from use in food producing animals but FFC, a fluorinated derivative of CHL, has a long history of use for the treatment of bovine respiratory diseases and to control pen disease in the local fish farming industry [52]. While 92.7% of STEC from water and sediments were susceptible to CHL, 16.4% were resistant and 49.1% were of intermediate resistance to FFC. Surprisingly, the overall prevalence of ARDs associated with resistance to CHL/FFC (floR, 1.20%) or CHL (catA, 0.60%) was very low. Moreover, none were detected in the MDR isolate 373-O165:H25 which was resistant to FFC. Consequently, the extent and molecular basis of resistance to phenicols in the STEC isolates could not be explained on the basis of the analyses performed in the present work and are clearly deserving of further investigation.
In contrast, ARDs associated with resistance to tetracyclines (tetA/B), penicillins (blaTEM), aminoglycosides (aac, aph) and folate pathway inhibitors (dfrA, sul1/2) were detected in the sequences of environmental and clinical STEC at frequencies that were not significantly different (p > 0.05). Several studies have shown that sul1/2, dfr, aac, aph, tetA/B and blaTEM ARDs are widely distributed in E. coli from aquatic environments impacted by human activity [53,54,55,56]. The β-lactamase encoding blaTEM-1B was detected at the highest frequency in the present work, although two MDR isolates lacking the ARD were of intermediate resistance to AMP, possibly due to the presence of undetected or unknown efflux pumps which can lessen susceptibility to antibiotics [57,58]. The AmpC β-lactamase encoding blaCMY-2 ARD associated with resistance to third-generation cephalosporins is reported to be widespread in E. coli [59] but was only detected in one MDR isolate (385-O111-NM). Interestingly, this isolate displayed the most complex antibiophenotype and was resistant to both TIO and FFC. The floR gene is often situated on chromosomal DNA, although carriage of both floR and blaCMY-2 on an IncA/C plasmid has been reported [60,61]. IncA/C plamids were not detected in the isolates, but other incompatibility types previously shown to carry blaCMY-2 genes in E. coli from human and animal sources were found, including IncB/O and IncFII [24,62].
In conclusion, the analysis of STEC isolates from surface waters and sediments in a region of intensive agricultural activity revealed AR phenotypes and the presence of ARDs comparable to those reported in generic E. coli from food animals. This finding suggests that common selective pressures are contributing to the development of AR in commensal and pathogenic variants of the species, and that aquatic environments can serve as reservoirs for the maintenance of resistant STEC, associated ARDs and their plasmids. Hence, the use of surface waters for the irrigation of food crops could introduce the risk of exposure, particularly in commodities that are commonly consumed raw.

4. Materials and Methods

4.1. Bacterial Strains and Culture Conditions

STEC isolates (n = 55) from 21 different serotypes isolated from surface waters and sediments from four watersheds in a mixed use urban–agricultural landscape located in the Fraser Valley of British Columbia, Canada, by Nadya et al. [29] and Falardeau et al. [63] were examined in this study (Table S1). The isolates were recovered from sampling sites located upstream or downstream of water sources used for the irrigation of various crops, including berry fruits and market vegetables. Samplings were carried out monthly over two years. All isolates were confirmed to be STEC using a hydrophobic grid membrane immunoblot method that targets both Stx1 and Stx2, and the presence of the associated genes was confirmed by multiplex PCR [29,63]. E. coli strain ATCC® 25922™ was used as a positive control in antibiotic resistance assays (Cedarlane Labs, Burlington, ON, Canada). All strains were stored at −80 °C in tryptic soy broth (TSB; Becton, Dickson and Company, Sparks, MD, USA) with 20% glycerol. Working cultures were prepared by application of thawed cultures to the surface of tryptic Soy agar (TSA; Becton, Dickson and Company), followed by incubation at 37 °C for 24 h. Working cultures were held at 4 °C for a maximum of 30 days.

4.2. Antibiotic Susceptibility Testing

Susceptibility to antibiotics was determined by the Kirby–Bauer disc diffusion assay [64]. Briefly, 70 µL aliquots of 18–24 h cultures grown in Mueller–Hinton broth (Becton, Dickson and Company) at 37 °C were mixed with 7 mL of 0.75% molten agar (44 °C) and added to the surface of petri plates containing Mueller–Hinton agar (Becton, Dickson and Company). The plates were incubated at an ambient temperature for 10 min and antibiotic discs were applied to the surface of the agar. The zones of inhibition were measured to the nearest millimeter after incubation at 37 °C for 24 h. The assignment of the isolates to the “susceptible”, “intermediate resistance” and “resistant” categories was based on the diameter of the inhibition zone according to the Clinical and Laboratory Standards Institute (CLSI) guidelines [65,66,67]. E. coli ATCC® 25922™ was used as the reference strain to verify assay performance. Eighteen antibiotics from several classes or groups of antimicrobials used in animal and/or human medicine were tested, including: aminoglycosides (amikacin (AMK; 20/10 μg), gentamicin (GEN; 10 μg), kanamycin (BCN; 30 μg), streptomycin (STR; 10 μg)); phenicols (chloramphenicol (CHL; 30 μg) and florfenicol (FFC; 30 µg)); quinolones (nalidixic acid (NAL; 30 μg)); second generation fluoroquinolones (ciprofloxacin (CIP; 5 μg)); tetracyclines (tetracycline HCl (TET; 30 μg)); folate pathway inhibitors (trimethoprim-sulfamethoxazole (SXT)); penicillin + ß-lactamase inhibitors (amoxicillin/clavulanic acid (AMC; 30 μg)); penicillins (ampicillin (AMP; 10 μg)); cephamycins (cefoxitin (FOX; 30 μg)); third generation cephalosporins (ceftiofur (TIO; 30 μg), ceftriaxone (CRO; 30 µg)); and carbapenems (ertapenem (ETP; 10 µg), imipenem (IPM; 10 μg) and meropenem (MEM; 10 µg)).

4.3. Genomic Analyses

The genomes of the STEC isolates from water/sediment were sequenced and assembled as described in Nadya et al. [31] and Fong et al. [68]. All sequences are available in NCBI’s Sequence Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra/ accessed on 5 January 2021) database under Bioproject Accessions PRJNA287560 and PRJNA649237 [69,70]. De novo assemblies were generated using shovill v. 1.0.94, using gsize 5000000 and mincov of 10% of the average genome coverage, as determined by the run_assembly_metrics.pl script in CG-pipeline (https://github.com/lskatz/CG-Pipeline accessed on 5 January 2021) [71]. An additional 112 assembled genomes from Canadian human clinical STEC isolates of variable serotypes from the same database were included in some of the analyses for comparative purposes. Tables S1 and S2 provide a list of accession numbers for each isolate, along with resistance genes, mutations and plasmids detected as described below. A single nucleotide polymorphism (SNP) phylogeny was generated using Parsnp v. 1.5.2 with E. coli O157:H7 EDL933 as a reference. The –x option to remove recombination regions was employed. An unrooted tree was visualized using Tree of Life [72]. Isolates of the same serotype with identical AMR pattern or lacking AMR determinants were omitted from the tree. Assemblies were queried for ARDs in the ResFinder database version from 28 October 2020, using staramr v. 0.4.0, with 90% nucleotide identity and 50% gene coverage cutoffs. Relevant mutations in the ampC promoter region, gyrA, gyrB, parC, and parE were detected by comparing to the PointFinder Database (https://bitbucket.org/genomicepidemiology/pointfinder_db/src/master/ accessed on 5 January 2021) using ARIBA v 2.12.0 where raw reads were available. For the clinical assemblies, the ampC, the ampC promoter, gyrA, gyrB, parC, and parE regions were extracted using a perl script (https://github.com/lskatz/lskScripts/blob/master/scripts/blastAndExtract.pl accessed on 5 January 2021), and subsequently aligned using MUSCLE in Geneious (version 6.0.6, Biomatters Ltd.). Relevant resistance conferring mutations were noted. Plasmid replicons were detected in each assembly using abricate 0.8.10 and a custom database based off the PlasmidFinder Enterobacteriacae database with a nucleotide identity cutoff of 90% and a coverage cutoff of 60% (database fasta available at https://github.com/StaPH-B/resistanceDetectionCDC/blob/master/plasmidDatabase.fasta accessed on 5 January 2021). Fisher’s exact test and Bonferroni correction to correct for multiple testing was calculated using R Studio (version 1.2.1335, http://www.rstudio.com/ accessed on 5 January 2021), and was used to determine if the prevalence of ARDs was non-randomly associated in water/sediment and clinical isolates [73].

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-6382/10/3/237/s1, Table S1: serotypes, IDs, accession numbers for genomic sequence reads and assemblies, resistance genes, resistance mutations and plasmids identified in Shiga toxigenic Escherichia coli isolates recovered from surface waters and sediments in the Lower Mainland of British Columbia between 2012–2016, Table S2: serotypes, IDs, accession numbers for genomic sequence reads and assemblies, resistance genes, resistance mutations and plasmids identified in Shiga toxigenic Escherichia coli clinical isolates used for comparison of antibiotic resistance gene frequencies.

Author Contributions

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

Funding

This research was funded by the Canadian Federal Genomics Research and Development Initiative (GRDI) and the National Sciences and Engineering Research Council of Canada [NSERC Discovery Grant RGPIN-2015-04871]. We acknowledge the contributions of the GRDI—Food and Water Safety Consortium (comprised of researchers from Agriculture and Agri-Food Canada, the Canadian Food Inspection Agency, Environment Canada, Health Canada, the National Research Council of Canada, and the Public Health Agency of Canada) for useful discussions and contributions.

Data Availability Statement

All data derived from the study are provided in the article and in supplementary material.

Acknowledgments

We gratefully acknowledge the help provided by Adam Olson of the PHAC National Microbiology Laboratory, Winnipeg, Manitoba for invaluable assistance with the sequencing, assembly, and annotation of STEC genomes.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Riley, L.W.; Remis, R.S.; Helgerson, S.D.; McGee, H.B.; Wells, J.G.; Davis, B.R.; Hebert, R.J.; Olcott, E.S.; Johnson, L.M.; Hargrett, N.T.; et al. hemorrhagic colitis associated with a rare Escherichia coli serotype. N. Engl. J. Med. 1983, 308, 681–685. [Google Scholar] [CrossRef] [PubMed]
  2. Tarr, P.I.; Gordon, C.A.; Chandler, W.L. Shiga-toxin-producing Escherichia coli and haemolytic uraemic syndrome. Lancet 2005, 365, 1073–1086. [Google Scholar] [CrossRef]
  3. Nathanson, S.; Kwon, T.; Elmaleh, M.; Charbit, M.; Launay, E.A.; Harambat, J.; Brun, M.; Ranchin, B.; Bandin, F.; Cloarec, S.; et al. Acute neurological involvement in diarrhea-associated hemolytic uremic syndrome. Clin. J. Am. Soc. Nephrol. 2010, 5, 1218–1228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Melton-Celsa, A.; Mohawk, K.; Teel, L.; O’Brien, A. Pathogenesis of Shiga-Toxin Producing Escherichia coli; Springer: Berlin/Heidelberg, Germany, 2011; Volume 357, pp. 67–103. [Google Scholar]
  5. Pires, S.M.; Majowicz, S.; Gill, A.; Devleesschauwer, B. Global and regional source attribution of Shiga toxin-producing Escherichia coli infections using analysis of outbreak surveillance data. Epidemiol. Infect. 2019, 147, e236. [Google Scholar] [CrossRef] [Green Version]
  6. Goldwater, P.N.; Bettelheim, K.A. Treatment of enterohemorrhagic Escherichia coli (EHEC) infection and hemolytic uremic syndrome (HUS). BMC Med. 2012, 10, 12. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, X.; McDaniel, A.D.; Wolf, L.E.; Keusch, G.T.; Waldor, M.K.; Acheson, D.W.K. Quinolone antibiotics induce shiga toxin–encoding bacteriophages, toxin production, and death in mice. J. Infect. Dis. 2000, 181, 664–670. [Google Scholar] [CrossRef]
  8. Agger, M.; Scheutz, F.; Villumsen, S.; Mølbak, K.; Petersen, A.M. Antibiotic treatment of verocytotoxin-producing Escherichia coli (VTEC) infection: A systematic review and a proposal. J. Antimicrob. Chemother. 2015, 70, 2440–2446. [Google Scholar] [CrossRef] [Green Version]
  9. Menne, J.; Nitschke, M.; Stingele, R.; Abu-Tair, M.; Beneke, J.; Bramstedt, J.; Bremer, J.P.; Brunkhorst, R.; Busch, V.; Dengler, R.; et al. Validation of treatment strategies for enterohaemorrhagic Escherichia coli O104:H4 induced haemolytic uraemic syndrome: Case-control study. BMJ 2012, 345, e4565. [Google Scholar] [CrossRef] [Green Version]
  10. Wieler, L.H.; Bauerfeind, R. STEC as a veterinary problem. Diagnostics and prophylaxis in animals. Methods Mol. Med. 2003, 73, 75–89. [Google Scholar] [PubMed]
  11. Premanandh, J.; Samara, B.S.; Mazen, A.N. Race against antimicrobial resistance requires coordinated action—An overview. Front. Microbiol. 2016, 6, 1536. [Google Scholar] [CrossRef]
  12. Laxminarayan, R.; Duse, A.; Wattal, C.; Zaidi, A.K.M.; Wertheim, H.F.L.; Sumpradit, N.; Vlieghe, E.; Hara, G.L.; Gould, I.M.; Goossens, H.; et al. Antibiotic resistance—The need for global solutions. Lancet Infect. Dis. 2013, 13, 1057–1098. [Google Scholar] [CrossRef] [Green Version]
  13. World Health Organization (WHO). Antimicrobial Resistance Global Report on Surveillance 2014. Available online: http://www.who.int/drugresistance/documents/surveillancereport/en (accessed on 18 January 2021).
  14. Kim, H.H.; Samadpour, M.; Grimm, L.; Clausen, C.R.; Besser, T.E.; Baylor, M.; Kobayashi, J.M.; Neill, M.A.; Schoenknecht, F.D.; Tarr, P.I. Characteristics of antibiotic-resistant Escherichia coli O157:H7 in Washington State, 1984–1991. J. Infect. Dis. 1994, 170, 1606–1609. [Google Scholar] [CrossRef]
  15. Zhao, S.; White, D.G.; Ge, B.; Ayers, S.; Friedman, S.; English, L.; Wagner, D.; Gaines, S.; Meng, J. Identification and characterization of integron-mediated antibiotic resistance among Shiga toxin-producing Escherichia coli isolates. Appl. Environ. Microbiol. 2001, 67, 1558–1564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Singh, R.; Schroeder, C.M.; Meng, J.; White, D.G.; McDermott, P.F.; Wagner, D.D.; Yang, H.; Simjee, S.; Debroy, C.; Walker, R.D.; et al. Identification of antimicrobial resistance and class 1 integrons in Shiga toxin-producing Escherichia coli recovered from humans and food animals. J. Antimicrob. Chemother. 2005, 56, 216–219. [Google Scholar] [CrossRef] [Green Version]
  17. Murinda, S.E.; Ebner, P.D.; Nguyen, L.T.; Mathew, A.G.; Oliver, S.P. Antimicrobial resistance and class 1 integrons in pathogenic Escherichia coli from dairy farms. Foodborne Pathog. Dis. 2005, 2, 348–352. [Google Scholar] [CrossRef] [PubMed]
  18. Guerra, B.; Junker, E.; Schroeter, A.; Helmuth, R.; Guth, B.E.C.; Beutin, L. Phenotypic and genotypic characterization of antimicrobial resistance in Escherichia coli O111 isolates. J. Antimicrob. Chemother. 2006, 57, 1210–1214. [Google Scholar] [CrossRef] [Green Version]
  19. Karmali, M.A.; Gannon, V.; Sargeant, J.M. Verocytotoxin-producing Escherichia coli (VTEC). Vet. Microbiol. 2010, 140, 360–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Magiorakos, A.-P.; Srinivasan, A.; Carey, R.B.; Carmeli, Y.; Falagas, M.E.; Giske, C.G.; Harbarth, S.; Hindler, J.F.; 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]
  21. Blair, J.M.A.; Webber, M.A.; Baylay, A.J.; Ogbolu, D.O.; Piddock, L.J.V. Molecular mechanisms of antibiotic resistance. Nat. Rev. Microbiol. 2015, 13, 42–51. [Google Scholar] [CrossRef] [PubMed]
  22. Stokes, H.W.; Hall, R.M. A novel family of potentially mobile DNA elements encoding site-specific gene-integration functions: Integrons. Mol. Microbiol. 1989, 3, 1669–1683. [Google Scholar] [CrossRef] [PubMed]
  23. Gillings, M.R. Integrons: Past, Present, and Future. Microbiol. Mol. Biol. Rev. 2014, 78, 257–277. [Google Scholar] [CrossRef] [Green Version]
  24. Carattoli, A. Resistance plasmid families in Enterobacteriaceae. Antimicrob. Agents Chemother. 2009, 53, 2227–2238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shintani, M.; Sanchez, Z.K.; Kimbara, K. Genomics of microbial plasmids: Classification and identification based on replication and transfer systems and host taxonomy. Front. Microbiol. 2015, 6, 242. [Google Scholar] [CrossRef] [PubMed]
  26. Losada, L.; Debroy, C.; Radune, D.; Kim, M.; Sanka, R.; Brinkac, L.M.; Kariyawasam, S.; Shelton, D.R.; Fratamico, P.M.; Kapur, V.; et al. Whole genome sequencing of diverse Shiga toxin-producing and non-producing Escherichia coli strains reveals a variety of virulence and novel antibiotic resistance plasmids. Plasmid 2016, 83, 8–11. [Google Scholar] [CrossRef]
  27. Köser, C.U.; Ellington, M.J.; Peacock, S.J. Whole-genome sequencing to control antimicrobial resistance. Trends Genet. 2014, 30, 401–407. [Google Scholar] [CrossRef] [Green Version]
  28. 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. [Google Scholar] [CrossRef] [Green Version]
  29. Nadya, S.; Delaquis, P.; Chen, J.; Allen, K.; Johnson, R.P.; Ziebell, K.; Laing, C.; Gannon, V.; Bach, S.; Topp, E. Phenotypic and genotypic characteristics of Shiga toxin-producing Escherichia coli isolated from surface waters and sediments in a Canadian urban-agricultural landscape. Front. Cell. Infect. Microbiol. 2016, 6, 36. [Google Scholar] [CrossRef] [PubMed]
  30. Hussein, H.S. Prevalence and pathogenicity of Shiga toxin-producing Escherichia coli in beef cattle and their products. J. Anim. Sci. 2007, 85, E63–E72. [Google Scholar] [CrossRef] [Green Version]
  31. Bai, X.; Hu, B.; Xu, Y.; Sun, H.; Zhao, A.; Ba, P.; Fu, S.; Fan, R.; Jin, Y.; Wang, H.; et al. Molecular and phylogenetic characterization of non-O157 Shiga toxin-producing Escherichia coli strains in China. Front. Cell. Infect. Microbiol. 2016, 6, 143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Liu, Y.; Fratamico, P.; Debroy, C.; Bumbaugh, A.C.; Allen, J.W. DNA Sequencing and identification of serogroup-specific genes in the Escherichia coli O118 O antigen gene cluster and demonstration of antigenic diversity but only minor variation in DNA sequence of the O antigen clusters of E. coli O118 and O151. Foodborne Pathog. Dis. 2008, 5, 449–457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wu, G.; Day, M.J.; Mafura, M.T.; Nunez-Garcia, J.; Fenner, J.J.; Sharma, M.; Van Essen-Zandbergen, A.; Rodríguez, I.; Dierikx, C.; Kadlec, K.; et al. Comparative analysis of ESBL-positive Escherichia coli isolates from animals and humans from the UK, The Netherlands and Germany. PLoS ONE 2013, 8, e75392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. 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] [Green Version]
  35. He, L.; Partridge, S.R.; Yang, X.; Hou, J.; Deng, Y.; Yao, Q.; Zeng, Z.; Chen, Z.; Liu, J. Complete nucleotide sequence of pHN7A8, an F33:A-:B- type epidemic plasmid carrying blaCTX-M-65, fosA3 and rmtB from China. J. Antimicrob. Chemother. 2013, 68, 46–50. [Google Scholar] [CrossRef] [Green Version]
  36. Johnson, T.J.; Bielak, E.M.; Fortini, D.; Hansen, L.H.; Hasman, H.; Debroy, C.; Nolan, L.K.; Carattoli, A. Expansion of the IncX plasmid family for improved identification and typing of novel plasmids in drug-resistant Enterobacteriaceae. Plasmid 2012, 68, 43–50. [Google Scholar] [CrossRef]
  37. Allen, K.J.; Laing, C.R.; Cancarevic, A.; Zhang, Y.; Mesak, L.R.; Xu, H.; Paccagnella, A.; Gannon, V.P.J.; Hoang, L. Characteristics of clinical Shiga toxin-producing Escherichia coli Isolated from British Columbia. BioMed Res. Int. 2013, 2013, 1–11. [Google Scholar] [CrossRef]
  38. Health Canada. Categorization of Antimicrobial Drugs Based on Importance in Human Medicine. Available online: https://www.canada.ca/en/health-canada/services/drugs-health-products/veterinary-drugs/antimicrobial-resistance/categorization-antimicrobial-drugs-based-importance-human-medicine.html (accessed on 18 January 2021).
  39. Maal-Bared, R.; Bartlett, K.H.; Bowie, W.R.; Hall, E.R. Phenotypic antibiotic resistance of Escherichia coli and E. coli O157 isolated from water, sediment and biofilms in an agricultural watershed in British Columbia. Sci. Total. Environ. 2013, 443, 315–323. [Google Scholar] [CrossRef] [PubMed]
  40. Meng, J.; Zhao, S.; Doyle, M.P.; Joseph, S.W. Antibiotic resistance of Escherichia coli O157:H7 and O157:NM isolated from animals, food, and humans. J. Food Prot. 1998, 61, 1511–1514. [Google Scholar] [CrossRef]
  41. Cabal, A.; Gómez-Barrero, S.; Porrero, C.; Bárcena, C.; López, G.; Cantón, R.; Gortazar, C.; Domínguez, L.; Alvarez, J. Assessment of virulence factors characteristic of human Escherichia coli pathotypes and antimicrobial resistance in O157:H7 and non-O157:H7 isolates from livestock in Spain. Appl. Environ. Microbiol. 2013, 79, 4170–4172. [Google Scholar] [CrossRef] [Green Version]
  42. Krüger, A.; Lucchesi, P.M.A.; Sanso, A.M.; Etcheverría, A.I.; Bustamante, A.V.; Burgán, J.; Fernández, L.; Fernandez, D.; Leotta, G.; Friedrich, A.W.; et al. Genetic characterization of Shiga toxin-producing Escherichia coli O26:H11 strains isolated from animal, food, and clinical samples. Front. Cell. Infect. Microbiol. 2015, 5, 74. [Google Scholar] [CrossRef] [Green Version]
  43. Amézquita-López, B.A.; Quiñones, B.; Soto-Beltrán, M.; Lee, B.G.; Yambao, J.C.; Lugo-Melchor, O.Y.; Chaidez, C. Antimicrobial resistance profiles of Shiga toxin-producing Escherichia coli O157 and Non-O157 recovered from domestic farm animals in rural communities in Northwestern Mexico. Antimicrob. Resist. Infect. Control. 2016, 5, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Day, M.; Doumith, M.; Jenkins, C.; Dallman, T.J.; Hopkins, K.L.; Elson, R.; Godbole, G.; Woodford, N. Antimicrobial resistance in Shiga toxin-producing Escherichia coli serogroups O157 and O26 isolated from human cases of diarrhoeal disease in England, 2015. J. Antimicrob. Chemother. 2017, 72, 145–152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Mora, A.; Blanco, J.E.; Blanco, M.; Alonso, M.P.; Dhabi, G.; Echeita, A.; González, E.A.; Bernárdez, M.I.; Blanco, J. Antimicrobial resistance of Shiga toxin (verotoxin)-producing Escherichia coli O157:H7 and non-O157 strains isolated from humans, cattle, sheep and food in Spain. Res. Microbiol. 2005, 156, 793–806. [Google Scholar] [CrossRef]
  46. Centers for Disease Control and Prevention (CDC). National Antimicrobial Resistance Monitoring System for Enteric Bacteria. Available online: https://www.cdc.gov/narms/index.html (accessed on 18 January 2021).
  47. British Columbia Center of Disease Control (BCCDC). Antimicrobial Resistance Trends in the Province of British Columbia. 2014. Available online: http://www.bccdc.ca/resource-gallery/Documents/AMR%202014%20Report-3NOV2015_FINAL.pdf (accessed on 18 January 2021).
  48. Redgrave, L.S.; Sutton, S.B.; Webber, M.A.; Piddock, L.J. Fluoroquinolone resistance: Mechanisms, impact on bacteria, and role in evolutionary success. Trends Microbiol. 2014, 22, 438–445. [Google Scholar] [CrossRef]
  49. Fasugba, O.; Gardner, A.; Mitchell, B.G.; Mnatzaganian, G. Ciprofloxacin resistance in community- and hospital-acquired Escherichia coli urinary tract infections: A systematic review and meta-analysis of observational studies. BMC Infect. Dis. 2015, 15, 545. [Google Scholar] [CrossRef] [Green Version]
  50. Public Health Agency of Canada (PHAC). Human Antimicrobial Use Report. 2014. Available online: https://www.canada.ca/en/public-health/services/publications/drugs-health-products/human-antimicrobial-use-report-2014.html (accessed on 1 February 2021).
  51. Public Health Agency of Canada (PHAC). Antimicrobial Resistance Surveillance System (CARSS) Report. 2016. Available online: https://www.canada.ca/en/public-health/services/publications/drugs-health-products/canadian-antimicrobial-resistance-surveillance-system-report-2016.html (accessed on 1 February 2021).
  52. Morrison, D.B.; Saksida, S. Trends in antimicrobial use in Marine Harvest Canada farmed salmon production in British Columbia (2003–2011). Can. Vet. J. 2013, 54, 1160–1163. [Google Scholar] [PubMed]
  53. Gomi, R.; Matsuda, T.; Matsumura, Y.; Yamamoto, M.; Tanaka, M.; Ichiyama, S.; Yoneda, M. Whole-genome analysis of antimicrobial-resistant and extraintestinal pathogenic Escherichia coli in river water. Appl. Environ. Microbiol. 2017, 83, e02703-16. [Google Scholar] [CrossRef] [Green Version]
  54. 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]
  55. Stange, C.; Sidhu, J.; Tiehm, A.; Toze, S. Antibiotic resistance and virulence genes in coliform water isolates. Int. J. Hyg. Environ. Health 2016, 219, 823–831. [Google Scholar] [CrossRef] [Green Version]
  56. 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]
  57. van Hoek, A.H.; Mevius, D.; Guerra, B.; Mullany, P.; Roberts, A.P.; Aarts, H.J. Acquired antibiotic resistance genes: An overview. Front. Microbiol. 2011, 2, 1–203. [Google Scholar]
  58. Poole, K. Efflux-mediated antimicrobial resistance. In Antibiotic Discovery and Development; Dougherty, T.J., Pucci, M.J., Eds.; Springer: Boston, MA, USA, 2012; pp. 349–395. [Google Scholar]
  59. Jacoby, G.A. AmpC β-lactamases. Clin. Microbiol. Rev. 2009, 22, 161–182. [Google Scholar] [CrossRef] [Green Version]
  60. Blickwede, M.; Schwarz, S. Molecular analysis of florfenicol-resistant Escherichia coli isolates from pigs. J. Antimicrob. Chemother. 2004, 53, 58–64. [Google Scholar] [CrossRef] [PubMed]
  61. Meunier, D.; Jouy, E.; Lazizzera, C.; Doublet, B.; Kobisch, M.; Cloeckaert, A.; Madec, J.-Y. Plasmid-borne florfenicol and ceftiofur resistance encoded by the floR and blaCMY-2 genes in Escherichia coli isolates from diseased cattle in France. J. Med. Microbiol. 2010, 59, 467–471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Bortolaia, V.; Hansen, K.H.; Nielsen, C.A.; Fritsche, T.R.; Guardabassi, L. High diversity of plasmids harbouring blaCMY-2 among clinical Escherichia coli isolates from humans and companion animals in the upper Midwestern USA. J. Antimicrob. Chemother. 2014, 69, 1492–1496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Falardeau, J.; Johnson, R.P.; Pagotto, F.; Wang, S. Occurrence, characterization, and potential predictors of verotoxigenic Escherichia coli, Listeria monocytogenes, and Salmonella in surface water used for produce irrigation in the Lower Mainland of British Columbia, Canada. PLoS ONE 2017, 12, e0185437. [Google Scholar] [CrossRef] [Green Version]
  64. Bauer, A.W.; Kirby, W.M.; Sherris, J.C.; Turck, M. Antibiotic susceptibility testing by a standardized single disk method. Tech. Bull. Regist. Med. Technol. 1966, 36, 49. [Google Scholar] [CrossRef]
  65. Clinical and Laboratory Standards Institute (CLSI). Performance Standards for Antimicrobial Susceptibility Testing: Twenty-First Informational Supplement; CLSI Document M100-S21; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2011; ISBN 1-56238-742-1. [Google Scholar]
  66. Clinical and Laboratory Standards Institute (CLSI). Performance Standards for Antimicrobial Disk and Dilution Susceptibility Tests for Bacteria Isolated From Animals, 5th ed.; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2018; ISBN 978-1-68440-009-6. [Google Scholar]
  67. Ahmed, M.O.; Clegg, P.D.; Williams, N.J.; Baptiste, K.E.; Bennett, M. Antimicrobial resistance in equine faecal Escherichia coli isolates from North West England. Ann. Clin. Microbiol. Antimicrob. 2010, 9, 12. [Google Scholar] [CrossRef] [Green Version]
  68. Fong, K.; Lu, Y.T.L.; Brenner, T.; Falardeau, J.; Wang, S. CRISPR and prophage diversity across Salmonella & verotoxin-producing Escherichia coli. BMC Genom. 2021. [Google Scholar]
  69. Leinonen, R.; Sugawara, H.; Shumway, M.; International Nucleotide Sequence Database Collaboration. The sequence read archive. Nucleic Acids Res. 2011, 39, D19–D21. [Google Scholar] [CrossRef] [Green Version]
  70. Barrett, T.; Clark, K.; Gevorgyan, R.; Gorelenkov, V.; Gribov, E.; Karsch-Mizrachi, I.; Kimelman, M.; Pruitt, K.D.; Resenchuk, S.; Tatusova, T.; et al. BioProject and BioSample databases at NCBI: Facilitating capture and organization of metadata. Nucleic Acids Res. 2012, 40, D57–D63. [Google Scholar] [CrossRef] [PubMed]
  71. Kislyuk, A.O.; Katz, L.S.; Agrawal, S.; Hagen, M.S.; Conley, A.B.; Jayaraman, P.; Nelakuditi, V.; Humphrey, J.C.; Sammons, S.A.; Govil, D.; et al. A computational genomics pipeline for prokaryotic sequencing projects. Bioinformatics 2010, 26, 1819–1826. [Google Scholar] [CrossRef] [Green Version]
  72. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256–W259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Walter, J.; Noll, I.; Feig, M.; Weiss, B.; Claus, H.; Werner, G.; Eckmanns, T.; Hermes, J.; Abu Sin, M. Decline in the proportion of methicillin resistance among Staphylococcus aureus isolates from non-invasive samples and in outpatient settings, and changes in the co-resistance profiles: An analysis of data collected within the Antimicrobial Resistance Surveillance Network, Germany 2010 to 2015. BMC Infect. Dis. 2017, 17, 169. [Google Scholar] [CrossRef] [Green Version]
Antibiotics 10 00237 g001 550
Figure 1. Proportion of 55 Shiga toxigenic Escherichia coli (STEC) isolates that were susceptible, of intermediate resistance or resistant to eighteen antibiotics: amikacin (AMK), amoxicillin/clavulanic acid (AMC), ampicillin (AMP), ceftiofur (TIO), cefoxitin (FOX), ceftriaxone (CRO), ciprofloxacin (CIP), gentamicin (GEN), kanamycin (BCN), nalidixic acid (NAL), streptomycin (STR), trimethoprim-sulfamethoxazole (SXT), chloramphenicol (CHL), tetracycline (TET), imipenem (IPM), ertapenem (ETP), meropenem (MEM) and florfenicol (FFC).
Figure 1. Proportion of 55 Shiga toxigenic Escherichia coli (STEC) isolates that were susceptible, of intermediate resistance or resistant to eighteen antibiotics: amikacin (AMK), amoxicillin/clavulanic acid (AMC), ampicillin (AMP), ceftiofur (TIO), cefoxitin (FOX), ceftriaxone (CRO), ciprofloxacin (CIP), gentamicin (GEN), kanamycin (BCN), nalidixic acid (NAL), streptomycin (STR), trimethoprim-sulfamethoxazole (SXT), chloramphenicol (CHL), tetracycline (TET), imipenem (IPM), ertapenem (ETP), meropenem (MEM) and florfenicol (FFC).
Antibiotics 10 00237 g001
Antibiotics 10 00237 g002 550
Figure 2. Maximum likelihood phylogeny constructed using core single nucleotide polymorphisms (SNPs) identified among the genomes of 55 STEC isolates recovered from water and sediments (isolate prefix FWS and LMFS) and 112 Canadian clinical STEC isolates and reference strain EDL933. Isolates of the same serotype with identical antimicrobial resistance (AMR) pattern or lacking AMR genes were omitted from the tree; tip colors denote different serotypes. Black squares in the matrix denote the presence of antibiotic resistance determinants common to the isolates. The scale corresponds to the number of substitutions per site. Resistance determinant functions: ampC: mutation in the promoter region of the chromosomal ampC gene; gyrA83 mutation: mutation in DNA gyrase codon 83; aac: aminoglycoside-3-acetyltransferase; aadA: aminoglycoside-3-adenylyltransferase; aph: aminoglycoside-3-phosphotransferase; blaCMY-2: AmpC β-lactamase; blaCTX-M: extended-spectrum β-lactamase; blaTEM: β-lactamase; catA: chloramphenicol acetyltransferase; dfrA: resistance to trimethoprim; floR: florfenicol/chloramphenicol exporter; fosA: fosfomycin thiol transferase; mphA: macrolide-2-phosphotransferase; sul1: dihydropteroate synthase; sul2: dihydropteroate synthase; tetA: tetracycline efflux protein; tetB: tetracycline efflux protein.
Figure 2. Maximum likelihood phylogeny constructed using core single nucleotide polymorphisms (SNPs) identified among the genomes of 55 STEC isolates recovered from water and sediments (isolate prefix FWS and LMFS) and 112 Canadian clinical STEC isolates and reference strain EDL933. Isolates of the same serotype with identical antimicrobial resistance (AMR) pattern or lacking AMR genes were omitted from the tree; tip colors denote different serotypes. Black squares in the matrix denote the presence of antibiotic resistance determinants common to the isolates. The scale corresponds to the number of substitutions per site. Resistance determinant functions: ampC: mutation in the promoter region of the chromosomal ampC gene; gyrA83 mutation: mutation in DNA gyrase codon 83; aac: aminoglycoside-3-acetyltransferase; aadA: aminoglycoside-3-adenylyltransferase; aph: aminoglycoside-3-phosphotransferase; blaCMY-2: AmpC β-lactamase; blaCTX-M: extended-spectrum β-lactamase; blaTEM: β-lactamase; catA: chloramphenicol acetyltransferase; dfrA: resistance to trimethoprim; floR: florfenicol/chloramphenicol exporter; fosA: fosfomycin thiol transferase; mphA: macrolide-2-phosphotransferase; sul1: dihydropteroate synthase; sul2: dihydropteroate synthase; tetA: tetracycline efflux protein; tetB: tetracycline efflux protein.
Antibiotics 10 00237 g002
Table
Table 1. Antibiotic resistance phenotypes of 55 STEC isolates from surface water and sediments in the Lower Mainland of British Columbia, Canada.
Table 1. Antibiotic resistance phenotypes of 55 STEC isolates from surface water and sediments in the Lower Mainland of British Columbia, Canada.
Isolate Number and SerotypeAntibiotic Resistance Phenotype 1Source
Pansusceptible
347-O174:H21NoneWater
005-O103:H2NoneWater
007-O103:H2NoneWater
008-O103:H2NoneWater
010-O109:H5NoneWater
033-O34:H32NoneWater
036-O34:H32NoneWater
039-O22:H8NoneWater
043-O153:NMNoneWater
047-O153:NMNoneWater
Intermediate resistant
292-O177:NMTETI FFCIWater
293-O157:NMTIOI GENI BCNI STRI TETI FFCIWater
294-O26:H11TETI FFCIWater
300-O26:H11FFCIWater
301-O26:H11FFCIWater
337-O8:H19TETI FFCIWater
338-O168:H8FFCR AMPI TIOI GENI BCNI STRI TETIWater
340-O116:H25AMCR AMPR AMKI TIOI FOXI CROI GENI BCNI STRI TETIWater
341-O5:NMAMCI AMPI TETI FFCIWater
343-O98:NMTIOI GENI STRI TETI FFCIWater
344-O5:NMFFCR AMCI AMPI TETIWater
345-O22:H8FFCR AMCI AMPI TETIWater
360-OR:H21FFCIWater
362-O5:NMAMCI FFCIWater
363-O8:H19TETI FFCIWater
365-O128:H2AMPI TETI FFCIWater
369-O98:NMSTRI CHLI TETI FFCIWater
370-O98:NMCHLI TETI FFCIWater
371-O157:H7FFCIWater
375-O157:H7FFCIWater
377-O103:H2FFCIWater
378-O103:H25AMPI TIOI FFCIWater
379-O26:H11FFCIWater
380-O165:NMFFCR AMPI TETIWater
381-O174:H8STRR AMPI TIOI GENI BCNIWater
383-O163:H19BCNI STRI TETI FFCIWater
384-O128:H2AMPI TETIWater
386-O8:H9AMPI FFCIWater
387-O103:H25FFCR AMPIWater
389-O103:H2AMPI FFCISediment
390-O103:H2AMPI FFCISediment
003-O116:H25AMXR, AMPRWater
012-O116:H25AMXR, AMPRWater
017-O76:H19STRR, TETRWater
021-O69:H11AMPR, STRR, SXTR, TETRWater
025-O69:H11AMPR, STRR, SXTR, TETIWater
029-O69:H11AMPR, STRR, SXTR, TETIWater
Multidrug resistant
296-O111:NMAMPR STRR SXTR TETR AMCI FFCIWater
298-O111:NMAMPR BCNR STRR TETR FFCR AMCI TIOIWater
299-O111:H8AMPR TIOR BCNR STRR TETR AMCI FFCIWater
356-O69:H11AMPR STRR SXTR TETR AMCI FFCIWater
373-O165:H25BCNR STRR TETR FFCR AMPI TIOIWater
374-O165:NMBCNR STRR TETR FFCR AMPIWater
385-O111:NMAMCR AMPR TIOR FOXR STRR SXTR CHLR TETR FFCR CROI GENI BCNIWater
391-O103:H2STRR SXTR CHLR TETR AMPI FFCISediment
1 R denotes resistance and I intermediate resistance according to the Clinical and Laboratory Standards Institute (CLSI) guidelines. Antibiotics: amikacin (AMK), amoxicillin/clavulanic acid (AMC), ampicillin (AMP), cefoxitin (FOX), ceftiofur (TIO), ceftriaxone (CRO), chloramphenicol (CHL), ciprofloxacin (CIP), ertapenem (ETP), gentamicin (GEN), imipenem (IPM), kanamycin (BCN), meropenem (MEM), nalidixic acid (NAL), streptomycin (STR), tetracycline (TET) and trimethoprim-sulfamethoxazole (SXT) and florfenicol (FFC).
Table
Table 2. Distribution of acquired resistance genes detected in whole genome sequences of 55 STEC isolates recovered from water and sediments in the Lower Mainland of British Columbia, Canada and 112 Canadian clinical STEC isolates. Probabilities that the prevalence of specific genes were nonrandomly distributed were calculated using Fisher’s exact test. Resistance determinant functions: ampC: mutation in the promoter region of the chromosomal ampC gene; gyrA83: mutation in DNA gyrase codon 83; aac: aminoglycoside 3-N-acetyltransferase; aadA: aminoglycoside-3-adenylyltransferase; aph: aminoglycoside-3-phosphotransferase; blaCMY-2: β-lactamase; blaCTX-M: extended spectrum β-lactamase; blaTEM: β-lactamase; catA: chloramphenicol acetyltransferase; dfrA: resistance to trimethoprim; floR: florfenicol/chloramphenicol exporter; fosA: fosfomycin thiol transferase; mphA: macrolide-2-phosphotransferase; sul1: dihydropteroate synthase; sul2: dihydropteroate synthase; tetA: tetracycline efflux protein; tetB: tetracycline efflux protein.
Table 2. Distribution of acquired resistance genes detected in whole genome sequences of 55 STEC isolates recovered from water and sediments in the Lower Mainland of British Columbia, Canada and 112 Canadian clinical STEC isolates. Probabilities that the prevalence of specific genes were nonrandomly distributed were calculated using Fisher’s exact test. Resistance determinant functions: ampC: mutation in the promoter region of the chromosomal ampC gene; gyrA83: mutation in DNA gyrase codon 83; aac: aminoglycoside 3-N-acetyltransferase; aadA: aminoglycoside-3-adenylyltransferase; aph: aminoglycoside-3-phosphotransferase; blaCMY-2: β-lactamase; blaCTX-M: extended spectrum β-lactamase; blaTEM: β-lactamase; catA: chloramphenicol acetyltransferase; dfrA: resistance to trimethoprim; floR: florfenicol/chloramphenicol exporter; fosA: fosfomycin thiol transferase; mphA: macrolide-2-phosphotransferase; sul1: dihydropteroate synthase; sul2: dihydropteroate synthase; tetA: tetracycline efflux protein; tetB: tetracycline efflux protein.
Prevalence (%) in All IsolatesWater/Sediment IsolatesClinical Isolatesp Value
Fisher’s Extract		Bonferroni Corrected
p Value
Gene (+)(−)(+)(−)
ampC1.8035201120.0340.585
gyrA831.8005531090.5521.000
aac0.6015401120.3291.000
aadA32.9935221100.3331.000
aph16.17124315970.1841.000
blaCMY-21.8015421101.0001.000
blaCTX-M0.6005511111.0001.000
blaTEM12.5774814981.0001.000
catA0.6015401120.3291.000
dfrA5.9974831090.0150.263
floR1.2015411110.5521.000
fosA2.4035211110.1051.000
mphA1.2025301120.1071.000
sul12.9935221100.3331.000
sul213.1794613990.4661.000
tetA8.98550101021.0001.000
tetB7.7884751070.0310.531
Table
Table 3. Acquired resistance determinants (ARDs) and plasmid replicon sequences detected in whole-genome sequences of eight multidrug resistant (MDR) STEC isolates from surface water and sediments.
Table 3. Acquired resistance determinants (ARDs) and plasmid replicon sequences detected in whole-genome sequences of eight multidrug resistant (MDR) STEC isolates from surface water and sediments.
IsolateAntibiotic Class/Group and Resistance ProfileARDs Detected in the Sequences and Putative Role in Resistance 1Plasmid Replicon Sequences 2,3
296-O111:NMAminoglycosides STRR
Folate pathway inhibitors SXTR
Penicillins, penicillin +β-lactam inhibitors AMPR AMCI
Phenicols FFCI
Tetracyclines TETR
NT 4		aadA2, aph(3″)-Ib, aph(6)-Id
drfA8, dfrA12, sul1, sul2
blaTEM-1B
-
tetA, tetB
mphA		ColE1 (2), IncB/O/K/Z,
IncFII(pHN7A8), IncX1
298-O111:NMAminoglycosides BCNR STRR
Cephalosporins TIOI
Penicillins, penicillin +β-lactam inhibitors
AMPR AMCI
Phenicols FFCR
Tetracyclines TETR
NT 4		aph(3′)-Ia, aph(3″)-Ib, aph(6)-Id
blaTEM-1B
blaTEM-1B
-
tetA
sul2		Col156, IncB/O/K/Z,
IncFII(pRSB107), IncQ1
299-O111:H8Aminoglycosides BCNR STRR
Cephalosprins TIOR
Penicillins, penicillin +β-lactam inhibitors
AMPR AMCI
Phenicols FFCI
Tetracyclines TETR
NT 4		aph(3′)-Ia, s aph(3″)-Ib, aph(6)-Id
blaTEM-1B
blaTEM-1B
-
tetA
sul2		Col156, IncB/O/K/Z,
IncFII(pRSB107), IncQ1
356-O69:H11Aminoglycosides STRR
Folate pathway inhibitors SXTR
Penicillins penicillin +β-lactam inhibitors
AMPR AMCI Phenicols FFCI
Tetracyclines TETR		aph(3″)-Ib, aph(6)-Id
dfrA8, sul2
blaTEM-1B
-
tetB		IncB/O/K/Z, IncFII(pHN7A8), IncFIB(AP001918)
373-O165:H25Aminoglycosides BCNR STRR
Cephalosporins TIOI
Penicillins AMPI
Phenicols FFCR
Tetracyclines TETR		aph(3′)-1a, aph(3″)-Ib, aph(6)-Id
-
-
-
tetB		IncFII, IncFIB(AP001918)
374-O165:NMAminoglycosides BCNR STRR
Phenicols FFCR
Penicillins AMPI
Tetracyclines TETR		aph(3′)-Ia, aph(3″)-Ib,aph(6)-Id
-
-
tetB		IncFII, IncFIB(AP001918),
385-O111:NMAminoglycosides STRR BCNI GENI
Cephalosporins TIOR CROI
Cephamycins FOXR
Penicillins, penicillin + β-lactam inhibitors
AMPR AMCR
Folate pathway inhibitors SXTR
Phenicols CHLR FFCR
Tetracyclines TETR
NT		aac(3)-VIa, aadA1, aadA2, aph(3″)-Ib,
aph(6)-Id
blaCMY-2
blaCMY-2
dfrA12, sul1, sul2
floR
tetA
mphA		ColEI, IncB/O/K/Z,
IncFII(pRSB107), IncY
391-O103:H2Aminoglycosides STRR
Folate pathway inhibitors SXTR
Phenicols CHLR FFCI
Penicillin AMPI
Tetracyclines TETR		aadA1, aph(3″)-Ib, aph(6)-Id
dfrA1, sul1, sul2
catA1
-
tetB		ColEI, IncB/O/K/Z
1 Resistance gene function: aac(3)-VIa: aminoglycoside 3-N-acetyltransferase; aadA1: aminoglycoside-3′-adenylyltransferase; aadA2: aminoglycoside-3-adenylyltransferase; aph(3′)-Ia: aminoglycoside-3- phosphotransferase; aph(3″)-1b: aminoglycoside-3-phosphotransferase; aph(6)-Id: aminoglycoside-3-phosphotransferase; blaTEM-1B: β-lactamase; blaCMY-2: AmpC β-lactamase-encoding β-lactamase; catA1: chloramphenicol acetyltransferase; dfrA1: resistance to trimethoprim; drfA8: resistance to trimethoprim; dfrA12: resistance to trimethoprim; floR: florfenicol/chloramphenicol exporter; mphA: macrolide-2-phosphotransferase, aph(3″)-Ib: streptomycin phosphotransferase; aph(6)-Id: streptomycin phosphotransferase; sul1: dihydropteroate synthase; sul2: dihydropteroate synthase; tetA: tetracycline efflux protein; tetB: tetracycline efflux protein. Some genes are listed multiple times where resistance to drug classes may be induced by the same gene. 2 Italicized parentheses denote the number of times the sequence was detected in the genome. 3 Predicted incompatibility group. Plasmid descriptions and putative role: ColEI: small plasmid, colicin-encoding. Col156: small plasmid, colicin-encoding [24,34]; IncB/O/K/Z: MDR and virulence genes [24]; IncFIB(AP001918): MDR and virulence genes [34]; IncFII(pRSB107): MDR and virulence genes [34]; IncFII(pHN7A8): MDR and virulence genes [35]; IncQ1: very broad host range plasmid, MDR genes [34]; IncX1: narrow host range conjugative plasmid, MDR and virulence genes [34,36]; IncY: Phage-like plasmid [34].4 NT: not tested.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, Y.; Chen, J.; Fong, K.; Nadya, S.; Allen, K.; Laing, C.; Ziebell, K.; Topp, E.; Carroll, L.M.; Wiedmann, M.; et al. Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape. Antibiotics 2021, 10, 237. https://doi.org/10.3390/antibiotics10030237

AMA Style

Ma Y, Chen J, Fong K, Nadya S, Allen K, Laing C, Ziebell K, Topp E, Carroll LM, Wiedmann M, et al. Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape. Antibiotics. 2021; 10(3):237. https://doi.org/10.3390/antibiotics10030237

Chicago/Turabian Style

Ma, Yvonne, Jessica Chen, Karen Fong, Stephanie Nadya, Kevin Allen, Chad Laing, Kim Ziebell, Ed Topp, Laura M. Carroll, Martin Wiedmann, and et al. 2021. "Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape" Antibiotics 10, no. 3: 237. https://doi.org/10.3390/antibiotics10030237

APA Style

Ma, Y., Chen, J., Fong, K., Nadya, S., Allen, K., Laing, C., Ziebell, K., Topp, E., Carroll, L. M., Wiedmann, M., Delaquis, P., & Wang, S. (2021). Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape. Antibiotics, 10(3), 237. https://doi.org/10.3390/antibiotics10030237

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