Corrigendum: Characterizing Antimicrobial Resistant Escherichia coli and Associated Risk Factors in a Cross-Sectional Study of Pig Farms in Great Britain
- 1Department of Bacteriology, Animal and Plant Health Agency, Weybridge, United Kingdom
- 2National Institute for Health Research, Health Protection Research Unit, University of Oxford in Partnership with Public Health England (PHE), Oxford, United Kingdom
- 3Department of Epidemiological Sciences, Animal and Plant Health Agency, Weybridge, United Kingdom
- 4Modernising Medical Microbiology Consortium, Nuffield Department of Medicine, John Radcliffe Hospital, University of Oxford, Oxford, United Kingdom
Combatting antimicrobial resistant (AMR) using a One-Health approach is essential as various bacteria, including Escherichia coli, a common bacteria, are becoming increasingly resistant and livestock may be a reservoir. The AMR gene content of 492 E. coli, isolated from 56 pig farms across Great Britain in 2014–2015, and purified on antibiotic selective and non-selective plates, was determined using whole genome sequencing (WGS). The E. coli were phylogenetically diverse harboring a variety of AMR profiles with widespread resistance to “old” antibiotics; isolates harbored up to seven plasmid Inc-types. None showed concurrent resistance to third-generation cephalosporins, fluoroquinolones and clinically relevant aminoglycosides, although ∼3% harbored AMR genes to both the former two. Transferable resistance to carbapenem and colistin were absent, and six of 117 E. coli STs belonged to major types associated with human disease. Prevalence of genotypically MDR E. coli, gathered from non-selective media was 35% and that of extended-spectrum-beta-lactamase E. coli was low (∼2% from non-selective). Approximately 72.6% of E. coli from ciprofloxacin plates and only 8.5% from the other plates harbored fluoroquinolone resistance due to topoisomerase mutations; the majority were MDR. In fact, multivariable analysis confirmed E. coli purified from CIP enrichment plates were more likely to be MDR, and suggested MDR isolates were also more probable from farms with high antibiotic usage, specialist finisher farms, and farms emptying their manure pits only after each batch. Additionally, farms from the South East were more likely to have MDR E. coli, whereas farms in Yorkshire and the Humber were less likely. Future investigations will determine whether suggested improvements such as better biosecurity or lower antimicrobial use decreases MDR E. coli on pig farms. Although this study focuses on pig farms, we believe the methodology and findings can be applied more widely to help livestock farmers in the United Kingdom and elsewhere to tackle AMR.
Introduction
Escherichia coli is a key species associated with AMR, including multidrug resistance (resistance to 3 or more antimicrobial classes; MDR) (World Health Organisation [WHO], 2014), threatening delivery of effective healthcare and challenging basic procedures used in modern human and veterinary medicine (O’Neill, 2016).
In Europe, since 2006, non-therapeutic use of antimicrobials has been limited in livestock (Cogliani et al., 2011). However, concerns remain that antibiotic use in food animals for therapeutics increases the risk of selection and emergence of AMR bacteria; and animals are a reservoir (Argudin et al., 2017). Recent livestock studies in Great Britain (GB), suggest that E. coli, chosen as an indicator organism for monitoring AMR in commensal bacteria in the gut flora of livestock, may carry resistance to high priority critically important antibiotics (HP-CIAs) such as cefotaxime or colistin (Anjum et al., 2011; Randall et al., 2014; Card et al., 2016; Duggett et al., 2017; Kirchner et al., 2017), that may spread horizontally to pathogens. Harmonized monitoring of AMR across Europe has indicated resistance of E. coli to HP-CIAs varies from country to country in humans, animals, and food (European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017); but typically only phenotypic analysis of AMR is performed on isolate(s) from a single representative animal on farm, recovered from selective and/or non-selective media. Recently, the fecal resistome of pigs and poultry from nine European countries was assessed using metagenomics, to define the on-farm AMR load (Munk et al., 2018). The study reported a higher AMR load in pigs than poultry, however, the organisms harboring the AMR genes was not defined, possibly due to limits of Illumina short read sequencing. Also, mobile genetic elements such as plasmids that play a key role in AMR dispersal could not be classified further.
The aim of this work was to characterize, using WGS, the AMR genes present in 492 E. coli from multiple pigs on farm, collected in a cross-sectional study, across GB. WGS also enabled phylogenetic diversity of the host E. coli to be explored, as well as putative AMR plasmids, which may play a key role in the transferring of AMR genes between E. coli. Presence of MDR isolates was analyzed with on-farm risk factors to identify features that affect their selection, as well as those that may help control their spread.
Materials and Methods
Farm Recruitment and Sampling
Fifty-six farms in GB were recruited, from two different production types (38 Farrow-to-Finish and 18 Finisher-only farms). Farms were geographically diverse, and loosely representative of the finisher-only pig industry in GB; no farms were selected from the North West, Wales or Scotland (Supplementary Table S1) (Pig Health and Welfare Council, 2017). The farms selected were not random but an opportunistic selection of 53 farms participating in a Salmonella study, and these included 19 farms with low prevalence of Salmonella (Martelli et al., 2017; Smith et al., 2017, 2018). Three additional farms were recruited solely for this study. All farmers provided written consent to allow their pigs to be sampled and their laboratory results used for research purposes. Sampling of farms took place at 12 different abattoirs, between March 2014 and October 2015. A comprehensive questionnaire, completed by farmers, collected information on farming practices that included hygiene and disinfection of houses, biosecurity, as well as antibiotic usage. Ethical approval was not sought as sampling from carcases is deemed outside of the Animal (Scientific Procedures) Act 1986.
Sampling and Isolate Collection
At abattoirs, cecal content from 10 randomly selected healthy finishing pigs per herd/farm, were collected. For each herd, 0.5 g of pig cecal content from each pig was suspended into 22.5 ml of 0.1 M PBS (pH7.2) and the pooled pig cecal samples diluted up to 10–5 in PBS.
Aliquots of 100 μl were plated on to the following agar; Brilliance UTI agar (Oxoid, Basingstoke, United Kingdom) plates containing either 1 mg/L cefotaxime, 1 mg/L ciprofloxacin, and without antibiotics. Six morphologically distinct colonies were selected from each plate, however, if less than six colonies were isolated on either 1 mg/L cefotaxime or 1 mg/L ciprofloxacin plates, additional colonies from antibiotic-free plates were selected, to reach 18 isolates per farm. In addition, Brilliance carbapenem-resistant Enterobacteriaceae (CRE) Agar (Oxoid, Basingstoke, United Kingdom), was used. Purified colonies were subsequently stored at −80°C in MicroBank beads (Pro-Lab Diagnostics, Neston, Cheshire, United Kingdom). The presumptive purified E. coli were identified to species level using Matrix Assisted Laser Desorption Ionization-Time of Flight mass spectrometry (MALDI-ToF) (Bruker, Coventry, United Kingdom) or 16S rRNA sequencing (Edwards et al., 2012). Antimicrobial susceptibility testing was performed using the BSAC agar dilution method on the four E. coli isolated from CRE agar to determine susceptibility to three carbapenem antibiotics (Doripenem, Imipenem, and Meropenem), for all E. coli isolated from ciprofloxacin and cefotaxime antibiotic plates (Andrews, 2001). The results were interpreted using the European Committee on Antimicrobial Susceptibility Testing (EUCAST) epidemiological cut-off (ECOFF) and clinical breakpoint values.
Whole Genome Sequence Analysis
DNA was extracted and Illumina HiSeq 4000 System used to perform WGS as described previously (Stoesser et al., 2013) on the 503 isolates identified as E. coli by MALDI-ToF. Kraken was used to confirm MALDI-ToF speciation results (Wood and Salzberg, 2014), which identified 11 isolates as E. fergusonii, a species that has been reported in livestock (Wragg et al., 2009), and excluded from further analysis. The presence of acquired AMR genes in the WGS of isolates was determined by mapping unassembled reads using APHA SeqFinder pipeline (Anjum et al., 2016). The criteria for determining gene presence using APHA SeqFinder pipeline was 100% gene mapping to the reference, and allowing between 1 and 10 non-synonymous SNP. The MIC results and the correlation between genotype and phenotype was >98% for E. coli isolates; these results have been presented elsewhere (Stubberfield et al., 2019), and due to the high correlation only the AMR genotype data is presented in this paper. Genotypic multidrug resistance was defined according to the European Food Safety Authority (EFSA) European surveillance reports (European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017) as an isolate harboring genes from the following 3 or more resistance classes: Extended spectrum β-lactamase (ESBL), Ampicillin (AMP), Tetracycline (TET), Gentamicin - clinically relevant aminoglycoside (GEN), Azithromycin (AZM), Chloramphenicol (CHL), Trimethoprim (TMP), Sulphonamide (SUL), and Fluoroquinolones (FQN).
Sequenced genomes were assembled using SPAdes 3.7.0 11 (Bankevich et al., 2012) and annotated using PROKKA 1.11.12 (Seemann, 2014). The Abricate1 tool was used to determine which AMR and plasmid replicon genes were co-located on a single contig. BRIG was used to compare plasmids (Alikhan et al., 2011). The Multilocus Sequence Type (MLST) of E. coli isolates (Wirth et al., 2006) was determined using either SRTS2 (Inouye et al., 2014) or DTU pipeline2 (Larsen et al., 2012). SNIPPY3 was used to generate a whole genome single nucleotide polymorphisms (SNP) alignment produced from WGS data of 492 E. coli isolates using the default settings. A maximum-likelihood tree under the General Time Reversible model of nucleotide substitution with among-site rate heterogeneity model (GTR-G) and 100 bootstrap replicates was inferred from the whole genome of 272,385 SNP alignment using E. coli MG1655 (Accession number: U00096.2) as reference in RAxML-NG4. The tree was visualized and annotated in EvolView (He et al., 2016). The E. coli raw sequence data generated and analyzed in this work are available in the European Nucleotide Archive (ENA) under study accession number PRJEB26317.
Statistical Analysis
The risk factor dataset contained 72 farm-level variables describing basic farm demographics, feeding and watering information including additive practices, and general farm management, which was provided to farmers for completion prior to cecal content sample collection and detection of the AMR gene content of isolates. Antimicrobial usage data were collected, however, due to difference in the detail of information provided by the farms, these responses were categorized subjectively into high, medium and low usage by the authors based upon answers to the total number of daily animal doses of antimicrobials given to pigs and the total cost of antimicrobials in a 12-month period. No standardized criteria existed to categorize the farms, but in general low usage farms were defined by an absence of prophylactic use of antibiotics and a rough cost of antibiotics under £0.5 per pig. Medium usage farms were defined by prophylactic use for short periods of time and a rough cost of antibiotics up to £2 per pig. High usage farms were generally either using in-feed or in-water treatment for multiple pig stages or using a combination of prophylactic/in-feed or in-water treatment and large quantities of injectables (e.g., routine injections to sows), with a rough cost of up to £3 per pig. The dataset also included a binary indicator of the MDR isolates gathered from WGS data analysis, for antimicrobials considered by EFSA (European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017).
A simple descriptive assessment was completed of each variable, excluding variables with little or no data and creating new categorical variables, where appropriate. Following this, a univariable analysis utilizing a mixed-effects logistic regression model, accounting for the non-independence of multiple isolates per farm, was carried out to determine variables for inclusion (p-value < 0.25) in the multivariable model. Finally, to assess the effect of including variables at the multivariable level, an iterative, forward step-wise approach was adopted with variables that improved the model fit [assessed via p-value and Akaike Information Criterion (AIC) (Akaike, 1974)], were selected at each step of the model. The categorical variable, related to antimicrobial usage on each farm over a 12-month period, was classed as a priori and automatically retained in the model as was a variable accounting for the use of the four different agar plates used in the study. A p-value of less than 0.05 was considered to indicate a statistically significant difference. Potential confounding between variables were investigated by monitoring whether the addition of variables at each step inflated the Odds Ratio of the variables retained in the model. All statistical analyses were completed in Stata 12 (StataCorp, 2011).
Results
E. coli AMR Genes and MDR
Four hundred and ninety-two E. coli were recovered from 56 pig farms; 51.6% were recovered on non-selective plates and the remainder on antibiotic selective plates (Table 1). The APHA SeqFinder pipeline (Anjum et al., 2016) identified 62 different AMR gene variants in 84.6% (416 of 492) of the E. coli isolated which were from 54 farms (Table 1). For the remaining two farms, which were both farrow-to-finish farms, the E. coli did not harbor any AMR genes. For the remaining farms, E. coli harbored between 1 and 15 AMR genes (Figure 1), and as expected, isolates from the antibiotic-free plate generally had less AMR genes (average = 4.9) than those from antibiotic containing plates (average = 7), indicating co-selection of multiple resistances. The number of AMR genes present in E. coli from each farm was variable with the mean number of AMR genes per isolate per farm being 5.1 (Figure 2); there was no association between farms and the presence of specific AMR genes. The most common AMR genes were the tetracycline resistance genes [tet(A)], followed by β-lactamases (blaTEM–1b) and the streptomycin resistance genes (strAB) (Table 1). Only 9.3% of isolates harbored genes conferring resistance to clinically relevant aminoglycosides (gentamicin, amikacin, netilmicin, and tobramycin), and no 16S rRNA methyltransferase enzyme (16S RMTase) genes were identified (Magiorakos et al., 2012; Public Health England, 2017). Fourteen farms had E. coli isolates harboring two blaTEM variants, with isolates from farm MSG54, harboring three different variants.
Table 1. Summary of antimicrobial resistance genes from different classes identified in Escherichia coli isolates.
Figure 1. Correlation of AMR genes harbored by isolates with the selection plate. The number of AMR genes detected by WGS in 492 isolates from 56 pigs farms. The colored bars indicate the number of isolates selected on agar supplemented with: antibiotic-free plate (green); 1 mg/L ciprofloxacin (blue); 1 mg/L cefotaxime (orange); or from Brilliance carbapenem -resistant Enterobacteriaceae Agar plates (CRE; pink).
Figure 2. Relationship between numbers of AMR genes, farm region, and production type. Box plot representing the interquartile range or midspread of AMR genes found in isolates from a farm (numbered as MSG01–MSG35 and MSG37–MSG57). Farms were ordered according to their region and production type (farrow-to-finisher = blue and finisher-only = red). Number to left of bars indicate the number of isolates from farm. Center lines show the median; box limits indicate the 25th and 75th percentiles as determined by R software. Whiskers extend 1.5 times the interquartile range from the 25th and 75th percentiles, with outliers represented by dots.
National (Public Health England, 2017) and European (European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017) surveillance programs monitor AMR in commensal and pathogenic bacteria, mainly at the phenotypic level. To facilitate comparison between these datasets with our AMR genotypes, we queried genes present in the AMR classes included in these programs (Supplementary Table S2). Using this criterion, 35% of E. coli recovered on non-selective plates, from 40 farms, were genotypically MDR (Supplementary Figure S1). A higher proportion of isolates (87%) recovered on antibiotic selective plates, from 49 farms, were genotypically MDR. Sixty-four MDR genotypes were found; the most common genotype (ampicillin, tetracycline, sulphonamide, trimethoprim, and fluoroquinolone) was present in isolates from 12 farms in both production systems (Supplementary Table S3). Forty-eight MDR profiles (221 isolates) included resistance classes reported in E. coli from human blood infections such as ESBL, fluoroquinolone or clinically relevant aminoglycosides (Public Health England, 2017); however, no isolate harbored all 3 classes. Fifteen isolates harboring ESBL and fluoroquinolone resistance were detected on five farms.
We considered the linkage of AMR genes and the plasmid incompatibility types (Inc-types) by observing the co-location of these genes on the same contigs, using the de novo assembled genomes of isolates. Most isolates harbored between 1 and 7 (average 3) Inc-types, with 90 Inc-types and AMR gene combinations being present (Supplementary Table S4). The most common contigs harboring an AMR gene and plasmid Inc-types ranged in size from ∼4 to >100 Kb (Table 2); several of these contigs showed similarity to plasmids of diverse origins, mostly reported from Salmonella and E. coli. The number of farms and isolates in which these putative plasmid sequences were identified varied, with some (e.g., IncQ1, IncF, and IncI) present in isolates with different STs and farms. The putative IncQ1 plasmid harboring strAB and sul2 was the most common and highly conserved, found in 48 isolates, but it only showed ∼50% identity to the closest references available in databases (Supplementary Figure S2).
Resistance to HP-CIAs in E. coli
Only 28 E. coli isolates (18 from cefotaxime plates) harbored genes, which mediate resistance to extended-spectrum cephalosporins (ESC) (Table 1); three farms had E. coli harboring two ESC resistance gene variants (blaCTX–M–1 and blaCTX–M–15, blaCMY–2 and blaCTX–M–15, or blaCMY–2 and blaSHV–12). ESC resistance genes were found co-located on contigs with plasmid replicons, with some showing high homology to published plasmids (Table 2). For example isolates from farm MSG17 and MSG43 harbored blaCTX–M–1 on a plasmid with high identity to the p369 IncI plasmid (Table 2 and Supplementary Figure S2). An additional 16 E. coli, isolated from cefotaxime plates, harbored mutations in the chromosomal AmpC promoter region associated with increased β-lactamase expression (Caroff et al., 2000).
Chromosomal mutations in quinolone resistance determining regions (QRDR) of gyrA and/or parC (Yoshida et al., 1991) were identified in 72.6% of E. coli (from 44 farms) recovered from ciprofloxacin plates (Table 1), with MICs ranging from 2 to 64 mg/L; 53.1% of these isolates were clinically resistant to ciprofloxacin (>4 mg/L). Only 8.5% of E. coli (from 18 farms) recovered from the other agar plates harbored mutations in QRDR, with MICs ranging between 2 and 128 mg/L; 48% of these were clinically resistant to ciprofloxacin (Table 1). Of the QRDR mutations identified, S83L/D87N double substitution (n = 80) was the most common for GyrA, and S80I substitution (n = 80) most common for ParC; 73 of these isolates harbored double-serine mutation (Fuzi et al., 2017), and the majority were recovered on ciprofloxacin plates (Supplementary Table S5). In addition, three plasmid mediated quinolone resistant (PMQR) variants (Redgrave et al., 2014) were detected in 67 E. coli (13.6%) with the most common gene being qnrS1 (Table 1). Majority of E. coli (n = 56) harboring PMQR genes did not harbor QRDR mutations; the MICs levels ranged between 0.06 and 8 mg/L, with only four isolate showing clinical resistance (>1 mg/L). Nineteen isolates had a contig encoding both qnrS1 and IncI (Table 2), which showed ∼79% sequence identity to the pSTM2 plasmid found in Salmonella Typhimurium (Supplementary Figure S2), which may be present in these isolates.
Four E. coli from two farms, recovered on carbapenem selection media, did not harbor transmissible carbapenem resistance genes and were phenotypically susceptible to carbapenems (Supplementary Table S6); these isolates were highly resistant to cefotaxime which may have contributed to the growth on CRE plates (MIC > 32 mg/L). Plasmid mediated colistin resistance genes were not detected.
Phylogenetic Diversity
A whole genome SNP-based maximum-likelihood phylogenetic tree was constructed and separated isolates into two distinct clades (1 and 2) with the majority of isolates (n = 394) clustering within clade 1 (Figure 3). There was little evidence of clustering by phylogeography, production type or farm, with the only exception being farm MSG31, where all isolates clustered within S7.
Figure 3. Phylogenetic tree of Escherichia coli isolates. A SNP-based maximum-likelihood phylogenetic tree was constructed in RaxML-NG using 272,385 single nucleotide polymorphisms (SNP) present in the core genome. Clades and subclades are labeled on the relevant branches. Isolate names are colored according to farm production type (blue, farrow-to-finish; pink, finisher-only) and nodes at the end of the leaves are colored according to the six regions in England from where the E. coli isolates originate. The circles at the edge indicate the following information (working from inner circle to outer circle): (1) the top 8 MLST clonal complexes identified; (2) MLST ST744; (3) QRDR mutations in gyrA; (4) QRDR mutations in parC; (5) presence of plasmid quinolone resistance (qnr genes); (6) ciprofloxacin susceptibility using EUCAST ECOFF; (7) ciprofloxacin susceptibility using EUCAST clinical breakpoints; (8) presence of ESBL gene variants; and (9) the numbers EFSA AMR classes harboring different resistance genes in each isolate (genes conferring resistant to the same antimicrobial were grouped within the same class).
Escherichia coli isolates belonged to 117 multilocus sequence types (MLST) which included 19 clonal complexes and 17 new ST variants. Nearly half the isolates belonged to seven STs, but there were no apparent association between AMR profiles and ST (Figure 3 and Supplementary Table S7). On average, farms had E. coli from five STs, with the exception of farm MSG17, which harbored nine STs. Fifty-eight percent of isolates, from different geographic regions, were located in clade 2; the majority were ST744, which has one SNP variation from ST10. The remainder clustered in three sub-clades within clade 1 and majority of isolates (n = 15) belonged to ST162 (ST469 complex). Sixty-five percent of the STs (n = 76) identified have previously been observed in human E. coli (Supplementary Table S7), however, due to lack of information in the MLST database it is not clear whether these are of clinical origin. Nevertheless, 89 isolates belonged to seven STs commonly associated with extra-intestinal pathogenic E. coli (ExPEC) from United Kingdom bacteraemia and UTI in man (Gibreel et al., 2012; Manges and Johnson, 2012; Horner et al., 2014; Kallonen et al., 2017; Day et al., 2019), which included ST10, ST117, and ST131 (Supplementary Table S7). Ten ST131 E. coli, from four farms, clustered within clade S7. However, unlike MDR human ST131 isolates (Nicolas-Chanoine et al., 2014) they were phenotypically susceptible at clinical breakpoints to ciprofloxacin (Figure 3); although seven ST131-fimH22 isolates from 3 farms harbored a single GyrA SNP (S83L), which led to resistances above the ECOFF and these isolates were MDR but not to any other HP-CIA; three ST131-fimH298 isolates, all from the same farm, had no detectable AMR genes. Forty-seven of the 79 isolates (60%) from the remaining six STs showed a MDR genotype which included resistance to HP-CIAs (15 harbored resistance to gentamicin, 5 to ESBLs and 31 to fluoroquinolone, although none harbored resistance to all three classes to HP-CIAs).
Clustering was noted based on mutations in gyrA and parC. Half the isolates within clade 2 (n = 42) harbored double gyrA (S83L/D87N or LN) and a parC QRDR mutations, as did a number of isolates primarily from sub-clades 2 and 3 within clade 1 (7.6%, n = 30). Seventy-one isolates harboring both LN and parC mutations were clinically ciprofloxacin resistant; 94% (n = 67) of these were genotypically MDR (Figure 3).
On-Farm Risk Factors
Multivariable risk factor analysis was performed at sample level rather than farm level using information gathered in a detailed questionnaire on farming practices, including hygiene and disinfection of houses, biosecurity, as well as antibiotic usage; the type of agar MDR isolates were selected on was also included. There were 29 potential risk factors (p < 0.25) identified at univariable screening but the final multivariable logistic regression model retained only six risk factors associated with farms harboring genotypically MDR E. coli (Table 3). The retained variables included those identified as a priori (antimicrobial usage and ciprofloxacin enrichment plate) and four additional variables selected by forward stepwise selection that were significantly (p < 0.05) associated with the outcome and provided the best model fit (indicated by a lower AIC). Not surprisingly, the greater the antibiotic usage on farms, the greater the risk of MDR E. coli. Farms with a medium usage of antimicrobial products had almost three times the odds of MDR than farms with low usage (p = 0.050), whereas farms with high usage of antimicrobial products had over five times the odds of MDR (p < 0.001). Isolates recovered from ciprofloxacin plates had a higher probability of being MDR compared to isolates from the cefotaxime and non-selective plates (both p < 0.001); the carbapenem-resistant Enterobacteriaceae plate had too few records to produce a credible result.
Table 3. Multivariable risk factor analysis, assessing for the association between genotypic multi-drug resistance and variables of interest (n = 492, of which 296 were multi-drug resistant).
Specialist finisher farms were at a higher risk of MDR than breeder-finisher farms (p < 0.001), while those farms that emptied their manure pit after every batch were at greater risk than those that emptied it weekly (p = 0.023). Isolates from farms that on average, never cleaned and disinfected finisher pens between batches of pigs, were less likely to have MDR detected (p = 0.012). This answer was only provided by breeder-finisher farms and might have been due to continuous production; variables related to this factor did not appear to effect the selection of this factor, indicating that they were unlikely to fully explain this outcome. The geographical region of the farm was also selected as a significant factor, with farms in the South East being at lower risk of MDR (p = 0.012) and farms in Yorkshire and the Humber region at higher risk (p = 0.038) than those in the South West.
Discussion
This study aimed to characterize the AMR genes, using WGS, present in 492 E. coli from pig farms in GB. The most common AMR genes harbored by these phylogenetically diverse E. coli of different STs were to “old” antimicrobials, including aminoglycosides (which were not within the clinically important group), tetracyclines, ampicillin, sulphonamides, and trimethoprim (Hopkins et al., 2007; Card et al., 2014, 2015; Kirchner et al., 2014; Ulstad et al., 2016). None of the porcine E. coli concurrently harbored resistance genes to ESBL, fluoroquinolone, and clinically relevant aminoglycosides, but 15 isolates from five farms carried AMR genes to the former two antimicrobials. This is in contrast to the AST data reported by Public Health England (PHE) for humans between 2012 and 2016 where ∼5% of E. coli isolated from blood and cerebrospinal fluid showed resistance to third-generation cephalosporin, fluoroquinolones, and aminoglycosides, specifically gentamicin and/or tobramycin, which nevertheless have previously been reported from livestock (Szmolka et al., 2012). Also, the 16S RMTase genes reported from human clinical isolates (Public Health England, 2017), were absent. Importantly, several E. coli STs that have been isolated from humans were present, including seven collected from bacteraemia or feces of people admitted to UK hospitals, indicating some overlap between these compartments; many of these isolates were genotypically MDR and included resistance to HP-CIAs. Of the 10 ST131 isolates, 7 harbored the fimH22 fimbriae, progenitor of fimH30 associated with ST131 human pathogens (Ben Zakour et al., 2016). A previous study in retail meat also found that H22 strains were common and suggested the potential importance of E. coli ST131-fimH22 as a foodborne pathogen (Liu et al., 2018). While this study only examined E. coli, the repertoire of AMR genes is likely to be greater due to other bacteria present in the pig intestine; for example Moraxella species harboring mcr (AbuOun et al., 2017) and three methicillin resistant Staphylococcus aureus (Sharma et al., 2016) has been reported from pig samples collected in this study. However, a thorough survey of the diverse bacterial population within the pig intestine was beyond the scope of this work.
In this study, 35% of E. coli isolates recovered from non-selective plates were genotypically MDR. This is lower than the average reported across Europe (38.1%) and United Kingdom (51%) from non-selective media from pigs during 2015 (European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017). Differences in the number of farms sampled, methods of collection and isolation, and numbers of isolates tested could account for the variations; this study included 18 isolates/farm compared to only one isolate/farm included in EU surveillance. In addition, we used WGS to determine presence of AMR genes, which is highly predictive of AMR presence and could be more sensitive than phenotyping (Ellington et al., 2017; Stubberfield et al., 2019).
Carbapenem usage is not permitted in food-producing animals in the United Kingdom (Grace et al., 2015; Broadfoot et al., 2016), therefore absence of plasmid mediated carbapenem resistance genes, a serious public health threat was reassuring (Logan and Weinstein, 2017). The sales/usage in animals of 3rd generation cephalosporins is low (Broadfoot et al., 2016), so the low percentage of farms (16%) positive for ESBL-harboring E. coli was expected; although it was lower than the figures reported from EU ESBL E. coli monitoring (EU prevalence = 31.9%, UK prevalence = 21.7%) (Randall et al., 2014; European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC], 2017). No isolates harboring a blaCTX–M group 1 enzyme belonged to ST73, ST131, or ST95, STs associated with human infections (Gibreel et al., 2012; Horner et al., 2014; Kallonen et al., 2017).
The highest level of resistance to HP-CIAs was seen in fluoroquinolones and multivariable analysis indicated significantly higher likelihood of E. coli harboring MDR being isolated from ciprofloxacin compared to cefotaxime or non-selective plates. Thus, ciprofloxacin plates enabled detection not only of ciprofloxacin resistant E. coli, most often with QRDR mutations, but additional AMR genes. Furthermore, these MDR E. coli were less likely to be detected on other plates because they probably represented a small percentage of the total E. coli, as previously shown for mcr-harboring E. coli isolated from pig farms (Randall et al., 2018). QRDR mutations in Salmonella can confer protection against other antimicrobials (Webber et al., 2013) and the presence of double-serine QRDR mutations (GyrA83Ser and ParC80Ser) confers a selective advantage in lineages of several bacterial pathogens (Fuzi et al., 2017). We speculate these mutations may also aid on-farm persistence of QRDR E. coli enabling them to gain MDR status by acquiring more AMR plasmids. It could then lead to increased antibiotic use due to possible treatment failures, as supported by multivariable analysis, which showed increased risk of MDR E. coli being present on farms with high antibiotic usage. Farmers should limit or cease in-feed or in-water antibiotic treatment of pigs at herd level in order to reduce the risk of MDR. More targeted treatment of individual pigs would be preferable, although this would only be suitable where any reduction in use did not impact upon the welfare of the pigs.
Additionally, improvements to biosecurity or more frequent emptying of the manure pit on farm may lead to reduced risk of MDR E. coli being present and hence less need for treatment. The frequent emptying of the pit below the flooring in pig housing has been recommended previously as a factor associated with reducing the prevalence of Salmonella on pig farms (Beloeil et al., 2004). Never using disinfectants on average when cleaning a pen was found to be protective, which may infer disinfectants are being diluted inappropriately. It is known that bacteria can become resistant to disinfection when exposed to sublethal concentrations (Zou et al., 2014). Finisher farms were shown to be at higher risk of MDR, which may reflect the differences in management (such as the number of incoming movements and sources used) and pig types present, when compared to breeder-finisher farms. The region of the farm was also shown to be associated with MDR, which may reflect the density of pig farm production in those significant areas, as Yorkshire and Humber is a known high farm-density area, whereas the South East is sparsely populated. Although in this study each farm was sampled in more depth than many other cross-sectional national study, including for the EFSA harmonized monitoring of AMR. These results should nevertheless be treated with some caution as only 10 pigs were sampled per farm, and only from healthy slaughtered pigs, so it may not have been fully representative of these farms and appropriate for all farm types, such as those with only breeding pig stock.
As plasmids play a role in the mobilization of AMR genes between bacteria (Frost et al., 2005), putative plasmid contigs were used as surrogates for detecting circulating AMR plasmids. Isolates harbored between one and seven Inc-types and the most commonly circulating AMR plasmids belonged to the broad host range IncQ family (Loftie-Eaton and Rawlings, 2012), followed by IncI plasmids, which is frequently associated with AMR dissemination (Rozwandowicz et al., 2018). Future studies performed with hybrid short and long-read sequencing of isolates will enable a more complete resolution of plasmid genomes and other mobile genetic elements (MGEs) to accurately define the regions flanking AMR genes, and their transmission. However, high homology of plasmid genomes from this study with others present in the databases illustrate plasmids are part of a global network, cycling and disseminating AMR genes, as already made apparent by the mcr-plasmid phenomenon (Wyrsch et al., 2016; Duggett et al., 2017, 2018; Figueiredo et al., 2019).
Conclusion
We used genomics to detect MDR E. coli present on pig farms and combined it with multivariable analysis to identify factors affecting their selection, and possible control measures to help mitigate their transmission. Future studies can define common circulating plasmid genomes and other MGEs harboring AMR genes more fully and verify whether measures recommended to control their dissemination were successful.
Limitations
The farmer questionnaires were not completed fully by all farmers. Many of the questions were in a free format answer, which made it difficult to standardize and compare the results between farms, for example, it was difficult to gage levels of antibiotics used by farmers. Classification of antimicrobial usage was subjective and may have led to some misclassification bias. For the risk factor analysis, the sample size was large enough to detect large associations (Odds Ratio >6.5) with 95% confidence and 80% power. However, as the number of pigs sampled on each farm was relatively low, and the number of E. coli isolates per farm in the analysis differed, the results may not be fully representative of the population on that farm. A random effect was utilized in the risk factor model to account for the differing levels of clustering of isolates from farms and their non-independence.
Data Availability Statement
The datasets generated for this study can be found in the SRA PRJEB26317.
Ethics Statement
APHA undertakes research using animals under the Animal (Scientific Procedures) Act 1986 (ASPA) which includes having it’s own Animal Welfare and Ethics Board (AWERB). The AWERB looks at all aspects of the Science Division’s use of animals. As we had approval from the owners of the animals and the animals sampled were dead, this was outside ASPA (which only concerns live animals) and so it was considered that ethical approval from the committee was not required. Written informed consent was received from all participating farmers to allow their pigs to be sampled and their laboratory results and questionnaire data to be used for research purposes.
Author Contributions
MFA and RS conceived and directed the project and designed the work. MFA and DC funding acquisition. MA, ES, and EJS performed bacteriology laboratory work. MA, ES, EJS, and JN-G performed the gene and WGS analysis. MA, HO’C, EJS, and MFA interpreted the results. MA, RS, and MFA wrote the manuscript.
Funding
This work was funded by the Veterinary Medicines Directorate (VMD) to MFA under VMD grant numbers VM0518 and VM0533. Sequencing of the Gram-negative isolates was funded by the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Healthcare Associated Infections and Antimicrobial Resistance at Oxford University in partnership with Public Health England (PHE) Grant ID: HPRU-2012-10041.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank the farm owners and abattoirs who participated in the study. We thank the Food standards Agency and British Pig Executive for the collection of cecal contents at the abattoirs. OZ0344 DEFRA project supplied many of the farms and farm questionnaire information. We are very thankful to Dr. Rodrick Card for helpful discussions. We thank past APHA colleagues for contribution to sample collection and microbiological analysis: Mr. Jack Garretty, Miss Luisa Dormer, Mrs. Elizabeth Marier. We thank Dr. Nick Duggett and Mr. Jason Bramwell for assistance in aspects of the bioinformatics analysis. We would like to thank the sequencing centre at the Wellcome Trust Centre for Human Genetics and the Modernising Medical Microbiology (MMM) bioinformatics team in Oxford, United Kingdom.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00861/full#supplementary-material
Abbreviations
AMR, antimicrobial resistance; CRE, Brilliance carbapenem-resistant Enterobacteriaceae agar; HP-CIAs, highest priority critically important antibiotics; MDR, multidrug resistance (resistance to 3 or more antimicrobial classes); WGS, whole genome sequencing.
Footnotes
- ^ https://github.com/tseemann/abricate
- ^ https://cge.cbs.dtu.dk/services/MLST/
- ^ https://github.com/tseemann/snippy
- ^ https://github.com/amkozlov/raxml-ng
References
AbuOun, M., Stubberfield, E. J., Duggett, N. A., Kirchner, M., Dormer, L., Nunez-Garcia, J., et al. (2017). mcr-1 and mcr-2 variant genes identified in Moraxella species isolated from pigs in Great Britain from 2014 to 2015. J. Antimicrob. Chemother. 72, 2745–2749. doi: 10.1093/jac/dkx286
Akaike, H. (1974). A new look at the statistical model identification. IEEE Trans. Autom. Control AC 19, 716–723. doi: 10.1109/TAC.1974.1100705
Alikhan, N. F., Petty, N. K., Ben Zakour, N. L., and Beatson, S. A. (2011). BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12:402. doi: 10.1186/1471-2164-12-402
Andrews, J. M. (2001). Determination of minimum inhibitory concentrations. J. Antimicrob. Chemother. 48(Suppl. 1), 5–16.
Anjum, M. F., Choudhary, S., Morrison, V., Snow, L. C., Mafura, M., Slickers, P., et al. (2011). Identifying antimicrobial resistance genes of human clinical relevance within Salmonella isolated from food animals in Great Britain. J. Antimicrob. Chemother. 66, 550–559. doi: 10.1093/jac/dkq498
Anjum, M. F., Duggett, N. A., AbuOun, M., Randall, L., Nunez-Garcia, J., Ellis, R. J., et al. (2016). Colistin resistance in Salmonella and Escherichia coli isolates from a pig farm in Great Britain. J. Antimicrob. Chemother. 71, 2306–2313. doi: 10.1093/jac/dkw149
Argudin, M. A., Deplano, A., Meghraoui, A., Dodemont, M., Heinrichs, A., Denis, O., et al. (2017). Bacteria from animals as a pool of antimicrobial resistance genes. Antibiotics 6:12. doi: 10.3390/antibiotics6020012
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Beloeil, P. A., Fravalo, P., Fablet, C., Jolly, J. P., Eveno, E., Hascoet, Y., et al. (2004). Risk factors for Salmonella enterica subsp. enterica shedding by market-age pigs in French farrow-to-finish herds. Prev. Vet. Med. 63, 103–120. doi: 10.1016/j.prevetmed.2004.01.010
Ben Zakour, N. L., Alsheikh-Hussain, A. S., Ashcroft, M. M., Khanh Nhu, N. T., Roberts, L. W., Stanton-Cook, M., et al. (2016). Sequential acquisition of virulence and fluoroquinolone resistance has shaped the evolution of Escherichia coli ST131. mBio 7:e00347-16. doi: 10.1128/mBio.00347-316
Broadfoot, F., Harris, C., Brown, S., Healey, K., Grace, K., Reeves, H., et al. (2016). “UK VARSS 2015. UK veterinary antibiotic resistance and sales surveillance report”, in Veterinary Medicines Directorate & Animal and Plant Health Agency, (Addlestone, UK: Veterinary Medicines Directorate). Available online at: https://www.gov.uk/government/publications/veterinary-antimicrobial-resistance-and-sales-surveillance-2015
Card, R., Vaughan, K., Bagnall, M., Spiropoulos, J., Cooley, W., Strickland, T., et al. (2016). Virulence characterisation of Salmonella enterica isolates of differing antimicrobial resistance recovered from UK livestock and imported meat samples. Front. Microbiol. 7:640. doi: 10.3389/fmicb.2016.00640
Card, R. M., Mafura, M., Hunt, T., Kirchner, M., Weile, J., Rashid, M. U., et al. (2015). Impact of ciprofloxacin and clindamycin administration on gram-negative bacteria isolated from healthy volunteers and characterization of the resistance genes they harbor. Antimicrob. Agents Chemother. 59, 4410–4416. doi: 10.1128/AAC.00068-15
Card, R. M., Warburton, P. J., MacLaren, N., Mullany, P., Allan, E., and Anjum, M. F. (2014). Application of microarray and functional-based screening methods for the detection of antimicrobial resistance genes in the microbiomes of healthy humans. PLoS One 9:e86428. doi: 10.1371/journal.pone.0086428
Caroff, N., Espaze, E., Gautreau, D., Richet, H., and Reynaud, A. (2000). Analysis of the effects of -42 and -32 ampC promoter mutations in clinical isolates of Escherichia coli hyperproducing ampC. J. Antimicrob. Chemother. 45, 783–788. doi: 10.1093/jac/45.6.783
Cogliani, C., Goossens, H., and Greko, C. (2011). Restricting antimicrobial use in food animals: lessons from Europe. Microbe Magazine 6, 274–279. doi: 10.1128/microbe.6.274.1
Day, M. J., Hopkins, K. L., Wareham, D. W., Toleman, M. A., Elviss, N., Randall, L., et al. (2019). Extended-spectrum beta-lactamase-producing Escherichia coli in human-derived and foodchain-derived samples from England, Wales, and Scotland: an epidemiological surveillance and typing study. Lancet Infect. Dis. 19, 1325–1335. doi: 10.1016/S1473-3099(19)30273-7
Duggett, N. A., Randall, L. P., Horton, R. A., Lemma, F., Kirchner, M., Nunez-Garcia, J., et al. (2018). Molecular epidemiology of isolates with multiple mcr plasmids from a pig farm in Great Britain: the effects of colistin withdrawal in the short and long term. J. Antimicrob. Chemother. 73, 3025–3033. doi: 10.1093/jac/dky292
Duggett, N. A., Sayers, E., AbuOun, M., Ellis, R. J., Nunez-Garcia, J., Randall, L., et al. (2017). Occurrence and characterization of mcr-1-harbouring Escherichia coli isolated from pigs in Great Britain from 2013 to 2015. J. Antimicrob. Chemother. 72, 691–695. doi: 10.1093/jac/dkw477
Edwards, K. J., Logan, J. M., Langham, S., Swift, C., and Gharbia, S. E. (2012). Utility of real-time amplification of selected 16S rRNA gene sequences as a tool for detection and identification of microbial signatures directly from clinical samples. J. Med. Microbiol. 61(Pt 5), 645–652. doi: 10.1099/jmm.0.041764-0
Ellington, M. J., Ekelund, O., Aarestrup, F. M., Canton, R., Doumith, M., Giske, C., et al. (2017). The role of whole genome sequencing in antimicrobial susceptibility testing of bacteria: report from the EUCAST Subcommittee. Clin. Microbiol. Infect 23, 2–22. doi: 10.1016/j.cmi.2016.11.012
European Food Safety Authority [EFSA], and European Centre for Disease Prevention and Control [ECDC] (2017). The European Union Summary Report On. (Antimicrobial)Resistance in Zoonotic and Indicator Bacteria from Humans, Animals and Food in 2015. Solna: European Centre for Disease Prevention and Control.
Figueiredo, R., Card, R. M., Nunez-Garcia, J., Mendonca, N., da Silva, G. J., and Anjum, M. F. (2019). Multidrug-resistant Salmonella enterica isolated from food animal and foodstuff may also be less susceptible to heavy metals. Foodborne Pathog. Dis. 16, 166–172. doi: 10.1089/fpd.2017.2418
Frost, L. S., Leplae, R., Summers, A. O., and Toussaint, A. (2005). Mobile genetic elements: the agents of open source evolution. Nat. Rev. Microbiol. 3, 722–732. doi: 10.1038/nrmicro1235
Fuzi, M., Szabo, D., and Csercsik, R. (2017). Double-serine fluoroquinolone resistance mutations advance major international clones and lineages of various multi-drug resistant bacteria. Front. Microbiol. 8:2261. doi: 10.3389/fmicb.2017.02261
Gibreel, T. M., Dodgson, A. R., Cheesbrough, J., Fox, A. J., Bolton, F. J., and Upton, M. (2012). Population structure, virulence potential and antibiotic susceptibility of uropathogenic Escherichia coli from Northwest England. J. Antimicrob. Chemother. 67, 346–356. doi: 10.1093/jac/dkr451
Grace, K., Harris, C., Reeves, H., Healey, K., Coyne, L., and Teale, C. (2015). “UK VARSS 2014. UK veterinary antibiotic resistance and sales surveillance report”, in Veterinary Medicines Directorate & Animal and Plant Health Agency, (Addlestone, UK: Veterinary Medicines Directorate). Available online at: https://www.gov.uk/government/publications/veterinary-antimicrobial-resistance-and-sales-surveillance-2014
He, Z., Zhang, H., Gao, S., Lercher, M. J., Chen, W. H., and Hu, S. (2016). Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees. Nucleic Acids Res. 44, W236–W241. doi: 10.1093/nar/gkw370
Hopkins, K. L., Batchelor, M. J., Anjum, M., Davies, R. H., and Threlfall, E. J. (2007). Comparison of antimicrobial resistance genes in nontyphoidal salmonellae of serotypes enteritidis, hadar, and virchow from humans and food-producing animals in England and wales. Microb. Drug Resist. 13, 281–288. doi: 10.1089/mdr.2007.779
Horner, C., Fawley, W., Morris, K., Parnell, P., Denton, M., and Wilcox, M. (2014). Escherichia coli bacteraemia: 2 years of prospective regional surveillance (2010-12). J. Antimicrob. Chemother. 69, 91–100. doi: 10.1093/jac/dkt333
Inouye, M., Dashnow, H., Raven, L. A., Schultz, M. B., Pope, B. J., Tomita, T., et al. (2014). SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 6:90. doi: 10.1186/s13073-014-0090-6
Kallonen, T., Brodrick, H. J., Harris, S. R., Corander, J., Brown, N. M., Martin, V., et al. (2017). Systematic longitudinal survey of invasive Escherichia coli in England demonstrates a stable population structure only transiently disturbed by the emergence of ST131. Genome Res. 27, 1437–1449. doi: 10.1101/gr.216606.116
Kirchner, M., Lemma, F., Randall, L., and Anjum, M. F. (2017). Loop-mediated isothermal amplification for extended spectrum beta-lactamase gene detection in poultry carcase. Vet. Rec. 181:119. doi: 10.1136/vr.104150
Kirchner, M., Mafura, M., Hunt, T., Abu-Oun, M., Nunez-Garcia, J., Hu, Y., et al. (2014). Antimicrobial resistance characteristics and fitness of Gram-negative fecal bacteria from volunteers treated with minocycline or amoxicillin. Front. Microbiol. 5:722. doi: 10.3389/fmicb.2014.00722
Larsen, M. V., Cosentino, S., Rasmussen, S., Friis, C., Hasman, H., Marvig, R. L., et al. (2012). Multilocus sequence typing of total-genome-sequenced bacteria. J. Clin. Microbiol. 50, 1355–1361. doi: 10.1128/JCM.06094-11
Liu, C. M., Stegger, M., Aziz, M., Johnson, T. J., Waits, K., Nordstrom, L., et al. (2018). Escherichia coli ST131-H22 as a foodborne uropathogen. mBio 9:e00470-18. doi: 10.1128/mBio.00470-18
Loftie-Eaton, W., and Rawlings, D. E. (2012). Diversity, biology and evolution of IncQ-family plasmids. Plasmid 67, 15–34. doi: 10.1016/j.plasmid.2011.10.001
Logan, L. K., and Weinstein, R. A. (2017). The epidemiology of carbapenem-resistant Enterobacteriaceae: the impact and evolution of a global menace. J. Infect. Dis. 215(Suppl._1), S28–S36. doi: 10.1093/infdis/jiw282
Magiorakos, A. P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x
Manges, A. R., and Johnson, J. R. (2012). Food-borne origins of Escherichia coli causing extraintestinal infections. Clin. Infect. Dis. 55, 712–719. doi: 10.1093/cid/cis502
Martelli, F., Lambert, M., Butt, P., Cheney, T., Tatone, F. A., Callaby, R., et al. (2017). Evaluation of an enhanced cleaning and disinfection protocol in Salmonella contaminated pig holdings in the United Kingdom. PLoS One 12:e0178897. doi: 10.1371/journal.pone.0178897
Munk, P., Knudsen, B. E., Lukjacenko, O., Duarte, A. S. R., Van Gompel, L., Luiken, R. E. C., et al. (2018). Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries. Nat. Microbiol. 3, 898–908. doi: 10.1038/s41564-018-0241-4
Nicolas-Chanoine, M. H., Bertrand, X., and Madec, J. Y. (2014). Escherichia coli ST131, an intriguing clonal group. Clin. Microbiol. Rev. 27, 543–574. doi: 10.1128/CMR.00125-13
O’Neill, J. (2016). Tackling Drug-Resistant Infections Globally: Final Report and Recommendations. The Review on Antimicrobial Resistance. London: HM Government and the Wellcome Trust. Available online at: https://amr-review.org/
Pig Health and Welfare Council (2017). Pig Health and Welfare Council Biennial Report 2015-2016. Dhaka: Pig Health and Welfare Council.
Public Health England (2017). English Surveillance Programme for Antimicrobial Utilisation and Resistance (ESPAUR) Report 2017. London: Public Health England.
Randall, L. P., Horton, R. A., Lemma, F., Martelli, F., Duggett, N. A. D., Smith, R. P., et al. (2018). Longitudinal study on the occurrence in pigs of colistin-resistant Escherichia coli carrying mcr-1 following the cessation of use of colistin. J. Appl. Microbiol. 125, 596–608. doi: 10.1111/jam.13907
Randall, L. P., Lemma, F., Rogers, J. P., Cheney, T. E., Powell, L. F., and Teale, C. J. (2014). Prevalence of extended-spectrum-beta-lactamase-producing Escherichia coli from pigs at slaughter in the UK in 2013. J Antimicrob Chemother 69, 2947–2950. doi: 10.1093/jac/dku258
Redgrave, L. S., Sutton, S. B., Webber, M. A., and Piddock, L. J. (2014). Fluoroquinolone resistance: mechanisms, impact on bacteria, and role in evolutionary success. Trends Microbiol. 22, 438–445. doi: 10.1016/j.tim.2014.04.007
Rozwandowicz, M., Brouwer, M. S. M., Fischer, J., Wagenaar, J. A., Gonzalez-Zorn, B., Guerra, B., et al. (2018). Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae. J. Antimicrob. Chemother. 73, 1121–1137. doi: 10.1093/jac/dkx488
Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153
Sharma, M., Nunez-Garcia, J., Kearns, A. M., Doumith, M., Butaye, P. R., Argudin, M. A., et al. (2016). Livestock-associated methicillin resistant Staphylococcus aureus (LA-MRSA) clonal complex (CC) 398 isolated from UK animals belong to european lineages. Front. Microbiol. 7:1741. doi: 10.3389/fmicb.2016.01741
Smith, R. P., Andres, V., Dormer, L., Gosling, R., Oastler, C., and Davies, R. H. (2017). Study of the impact on Salmonella of moving outdoor pigs to fresh land. Epidemiol. Infect. 145, 1983–1992. doi: 10.1017/S0950268817000930
Smith, R. P., Andres, V., Martelli, F., Gosling, B., Marco-Jimenez, F., Vaughan, K., et al. (2018). Maternal vaccination as a Salmonella Typhimurium reduction strategy on pig farms. J. Appl. Microbiol. 124, 274–285. doi: 10.1111/jam.13609
Stoesser, N., Batty, E. M., Eyre, D. W., Morgan, M., Wyllie, D. H., Del Ojo Elias, C., et al. (2013). Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data. J. Antimicrob. Chemother. 68, 2234–2244. doi: 10.1093/jac/dkt180
Stubberfield, E., AbuOun, M., Sayers, E., O’Connor, H. M., Card, R. M., and Anjum, M. F. (2019). Use of whole genome sequencing of commensal Escherichia coli in pigs for antimicrobial resistance surveillance, United Kingdom, 2018. Eurosurveillance 24:1900136. doi: 10.2807/1560-7917.ES.2019.24.50.1900136
Szmolka, A., Anjum, M. F., La Ragione, R. M., Kaszanyitzky, E. J., and Nagy, B. (2012). Microarray based comparative genotyping of gentamicin resistant Escherichia coli strains from food animals and humans. Vet. Microbiol. 156, 110–118. doi: 10.1016/j.vetmic.2011.09.030
Ulstad, C. R., Solheim, M., Berg, S., Lindbaek, M., Dahle, U. R., and Wester, A. L. (2016). Carriage of ESBL/AmpC-producing or ciprofloxacin non-susceptible Escherichia coli and Klebsiella spp. in healthy people in Norway. Antimicrob. Resist. Infect. Control. 5:57. doi: 10.1186/s13756-016-0156-x
Webber, M. A., Ricci, V., Whitehead, R., Patel, M., Fookes, M., Ivens, A., et al. (2013). Clinically relevant mutant DNA gyrase alters supercoiling, changes the transcriptome, and confers multidrug resistance. mBio 4:e00273-13. doi: 10.1128/mBio.00273-13
Wirth, T., Falush, D., Lan, R., Colles, F., Mensa, P., Wieler, L. H., et al. (2006). Sex and virulence in Escherichia coli: an evolutionary perspective. Mol. Microbiol. 60, 1136–1151. doi: 10.1111/j.1365-2958.2006.05172.x
Wood, D. E., and Salzberg, S. L. (2014). Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 15:R46. doi: 10.1016/j.mimet.2016.01.011
World Health Organisation [WHO] (2014). World Health Organization Antimicrobial Resistance: Global Report on Surveillance 2014. Geneva: WHO.
Wragg, P., La Ragione, R. M., Best, A., Reichel, R., Anjum, M. F., Mafura, M., et al. (2009). Characterisation of Escherichia fergusonii isolates from farm animals using an Escherichia coli virulence gene array and tissue culture adherence assays. Res. Vet. Sci. 86, 27–35. doi: 10.1016/j.rvsc.2008.05.014
Wyrsch, E. R., Roy Chowdhury, P., Chapman, T. A., Charles, I. G., Hammond, J. M., and Djordjevic, S. P. (2016). Genomic microbial epidemiology is needed to comprehend the global problem of antibiotic resistance and to improve pathogen diagnosis. Front. Microbiol. 7:843. doi: 10.3389/fmicb.2016.00843
Yoshida, H., Bogaki, M., Nakamura, M., Yamanaka, L. M., and Nakamura, S. (1991). Quinolone resistance-determining region in the DNA gyrase gyrB gene of Escherichia coli. Antimicrob. Agents Chemother. 35, 1647–1650. doi: 10.1128/aac.35.8.1647
Keywords: antimicrobial resistance, multidrug resistance, plasmids, epidemiology, risk factor analysis, pig farms, Escherichia coli, Great Britain
Citation: AbuOun M, O’Connor HM, Stubberfield EJ, Nunez-Garcia J, Sayers E, Crook DW, Smith RP and Anjum MF (2020) Characterizing Antimicrobial Resistant Escherichia coli and Associated Risk Factors in a Cross-Sectional Study of Pig Farms in Great Britain. Front. Microbiol. 11:861. doi: 10.3389/fmicb.2020.00861
Received: 15 May 2019; Accepted: 09 April 2020;
Published: 25 May 2020.
Edited by:
Ziad Daoud, University of Balamand, LebanonReviewed by:
Abasiofiok Mark Ibekwe, United States Department of Agriculture (USDA), United StatesLee W. Riley, University of California, Berkeley, United States
Torahiko Okubo, Hokkaido University, Japan
Miklos Fuzi, Semmelweis University, Hungary
Copyright © 2020 AbuOun, O’Connor, Stubberfield, Nunez-Garcia, Sayers, Crook, Smith and Anjum. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Muna F. Anjum, muna.anjum@apha.gov.uk