Skip to main content

Comparison of antimicrobial resistant Escherichia coli isolated from Irish commercial pig farms with and without zinc oxide and antimicrobial usage

Abstract

Background

The prophylactic use of antimicrobials and zinc oxide (ZnO) in pig production was prohibited by the European Union in 2022 due to potential associations between antimicrobial and heavy metal usage with antimicrobial resistance (AMR) and concerns regarding environmental pollution. However, the effects of their usage on the bacterial AMR profiles on commercial pig farms are still not fully understood and previous studies examining the effect of ZnO have reported contrasting findings. The objective of this study was to examine the effects of antimicrobial and ZnO usage on AMR on commercial pig farms. Faecal and environmental samples were taken on 10 Irish commercial farms, of which 5 farms regularly used ZnO and antimicrobials (amoxicillin or sulphadiazine-trimethoprim) for the prevention of disease. The other 5 farms did not use ZnO or any other form of prophylaxis. Escherichia coli numbers were quantified from all samples using non-supplemented and supplemented Tryptone Bile X-glucuronide agar.

Results

In total 351 isolates were phenotypically analysed, and the genomes of 44 AmpC/ESBL-producing E. coli isolates from 4 farms were characterised using whole-genome sequencing. Phenotypic analysis suggested higher numbers of multi-drug resistant (MDR) E. coli isolates on farms using prophylaxis. Furthermore, farms using prophylaxis were associated with higher numbers of isolates resistant to apramycin, trimethoprim, tetracycline, streptomycin, and chloramphenicol, while resistance to ciprofloxacin was more associated with farms not using any prophylaxis. Thirty-four of the 44 AmpC/ESBL-producing E. coli strains harboured the blaCTX-M-1 resistance gene and were multi drug resistant (MDR). Moreover, network analysis of plasmids and analysis of integrons showed that antimicrobial and biocide resistance genes were frequently co-located on mobile genetic elements, indicating the possibility for co-selection during antimicrobial or biocide usage as a contributor to AMR occurrence and persistence on farms.

Conclusions

The results of this study showed evidence that antimicrobial and ZnO treatment of pigs post-weaning can favour the selection and development of AMR and MDR E. coli. Co-location of resistance genes on mobile genetic elements was observed. This study demonstrated the usefulness of phenotypic and genotypic detection of antimicrobial resistance by combining sequencing and microbiological methods.

Background

Antimicrobial resistance (AMR) is a global threat to human, animal, and environmental health. According to the One Health principles which were originally formulated in 2004 at the “One World, One Health” symposium in New York, transmission of disease and spread of AMR can only be prevented by acknowledging interconnections between humans, animals, and the environment [1]. Consequently, great efforts have been taken on a global scale to reduce the development and spread of AMR, especially in livestock production, where antimicrobials are extensively used, and have been associated with increased AMR in several studies [2]. Moreover, potential linkages of resistance to heavy metals and AMR due to cross-resistance, co-regulation, and co-occurrence of resistance genes on mobile genetic elements (MGEs) may select for persistence of AMR, even in the absence of antimicrobial-induced selective pressure on pig farms [3, 4].

Escherichia coli is a ubiquitous commensal for both humans and food-producing animals; however non-harmful strains co-exist with pathotypes which can cause gastrointestinal and urinary tract infections [5, 6]. In pig production, enterotoxigenic E. coli (ETEC) is one of the causative agents associated with post-weaning diarrhoea, which can result in ecological costs and impact negatively on farm profitability and animal health [7]. Therefore, antimicrobials and zinc oxide (ZnO) are commonly used for prophylactic treatment of piglets to prevent proliferation of ETEC during weaning [8]. However, in the European Union the prophylactic use of antimicrobials and therapeutic use of ZnO (2000–4000 ppm) has been banned due to possible associations with AMR development and spread, but also to reduce environmental pollution with heavy metals and antimicrobial residues via the use of manure in agriculture [9, 10]. Nevertheless, the impact of ZnO on AMR in livestock production is not yet fully understood.

The rise of MDR pathogenic E. coli in livestock is a major concern for animal, human, and environmental health [11]. In addition, MDR E. coli strains, or those carrying MGEs with resistance to last resort antimicrobials (e.g., mcr-1 plasmid-encoded strains), constitute a risk to human health by their potential transmission to humans, mainly through the food chain [12]. Due to its relevance and the straightforward low-cost monitoring methods for E. coli, it is often used as a microbial indicator of AMR in farm studies [13].

Prophylactic use of antibiotics and/or ZnO at weaning has favoured the selection of MDR resistant strains, which usually carry MGEs. For instance, two studies by Bednorz, et al. [14] and Ciesinski, et al. [15] which examined E. coli isolates from pigs, observed associations between zinc supplementation and an increased number of MDR E. coli. As discussed by these two studies, one possible explanation for their findings could be the co-selection of E. coli strains carrying multiple ARGs and heavy metal resistance genes (MRGs), which are able to withstand zinc-induced selective pressure. Bednorz, et al. [14] reported an overall high diversity of E. coli in the porcine gut but indicated that this diversity can reduce during zinc administration, leaving primarily MDR E. coli strains [14]. Furthermore, ZnO also affects the host, inducing changes in intestinal morphology, immune system, hepatic proteome, mucins, and methylation, thus host-mediated selective pressures in response to ZnO treatment may additionally affect E. coli diversity in the porcine gut [16,17,18,19].

A number of studies have analysed the molecular mechanisms underlying the observed increase in AMR. In a study by Vahjen, et al. [20], increases in copy numbers of the tet(A) and sul1 antimicrobial resistance genes (ARGs) were observed in response to high doses of dietary ZnO. These ARGs confer resistance to tetracyclines and sulphonamides respectively, suggesting that resistance to tetracycline or sulphonamides could potentially be increased with ZnO usage. This may occur because for instance, some efflux pumps that confer resistance to zinc can also confer cross-resistance to tetracycline and vice versa [20, 21].

Despite the information provided by these studies, most of them were small scale studies, performed under controlled conditions, and may not have captured the complexity of the commercial farm environment and other factors in pig production, such as the background resistance level due to the frequent use of antimicrobials or heavy metals on farms [22]. In addition, the new legislation which bans the prophylactic use of antibiotics and heavy metals to treat E. coli infections in pigs, creates a new window of opportunity to evaluate not only their link to AMR development, but also the potential impact of their removal on AMR [23, 24]. Consequently, the aim of this study was to investigate AMR E. coli diversity as a biomarker for monitoring the impact of the systematic use of antimicrobials and/or ZnO, or their removal, on AMR development on commercial pig farms.

Results

E. coli enumeration

In the first part of this study, E. coli colony counts of faecal and swab samples were enumerated, using non-supplemented and antimicrobial supplemented Tryptone Bile X-glucuronide agar (TBX), a chromogenic agar for the detection and enumeration of E. coli. Overall no differences between E. coli concentrations on AB-ZnO farms and AB-ZnO free farms were observed, independent of the of the TBX media used. On non-supplemented TBX, E. coli counts reached mean values of 7.42–8.08 log10 CFU/g and 7.35–8.21 log10 CFU/g in faecal samples for AB-ZnO and AB-ZnO free farms, respectively (Fig. 1). Interestingly, differences between AB-ZnO and AB-ZnO free farms were observed when comparing E. coli concentrations derived from apramycin-supplemented plates. However, these differences did not reach significance (p > 0.05). At the day of weaning, faecal apramycin resistant E. coli numbers on AB-ZnO and AB-ZnO free farms were almost identical, with mean values of 4.43 log10 CFU/g and 4.22 log10 CFU/g, respectively. However, while the mean value of apramycin-resistant E. coli numbers decreased by 0.67 log10 CFU/g in AB-ZnO free farms after 2 weeks, on AB-ZnO farms there was an increase of 2.08 log10 CFU/g. Although differences between farms were not significant, apramycin-resistant E. coli numbers were slightly higher on AB-ZnO farms than AB-ZnO free farms.

Fig. 1
figure 1

Enumeration of E. coli counts originating from faecal samples. E. coli counts from faecal samples from AB-ZnO using and AB-ZnO free farms using non-supplemented, apramycin- supplemented, ciprofloxacin-supplemented, and cefotaxime-supplemented TBX agar. Means and standard deviations are indicated in red. E. coli counts between AB-ZnO using and AB-ZnO free farms were compared per sample using the Mann-Witney U test (no significant differences were observed). FW0 = faecal sample at day of weaning; FW1 = faecal sample 1 week after weaning; FW2 = faecal sample 2 week after weaning

For ciprofloxacin- and cefotaxime-resistant E. coli numbers no clear observations could be made, as these resistance types were infrequently observed on farms.

Phenotypic antimicrobial resistance determination

After assessing phenotypic resistance profiles of 361 E. coli isolates which were taken from supplemented TBX agar, 351 isolates remained for further comparison (Table 1).

Table 1 Comparison of antimicrobial resistance patterns of Escherichia coli isolated on supplemented TBX agar plates

When analysing the phenotypic resistance profiles of the 351 isolates, it was observed that certain resistance patterns occurred more frequently, independently of selective pressure, which was mediated by culturing isolates on supplemented TBX (Fig. 2). Moreover, it was observed that specific resistance profiles occurred more frequently on AB-ZnO or AB-ZnO free farms (Table 1). Based on the analysis of isolates derived from the apramycin-supplemented TBX, resistance to, trimethoprim (AB-ZnO free: 63.41%; AB-ZnO: 86.61%; p < 0.001), tetracycline (AB-ZnO free: 63.51%; AB-ZnO: 90.55%; p < 0.001), streptomycin (AB-ZnO free: 69.51%; AB-ZnO: 92.13%; p < 0.001), chloramphenicol (AB-ZnO free: 19.51%; AB-ZnO: 48.03%; p < 0.001), and multi-drug resistance (AB-ZnO free: 63.29%; AB-ZnO: 94.49%; p < 0.001) was observed more often in isolates originating from AB-ZnO using farms. Similar results were observed for isolates derived from cefotaxime-supplemented TBX, showing more frequently resistance to apramycin (AB-ZnO free: 35.29%; AB-ZnO: 97.22%; p < 0.001), tetracycline (AB-ZnO free: 47.06%; AB-ZnO: 100.00%; p < 0.001), streptomycin (AB-ZnO free: 35.29%; AB-ZnO: 91.67%; p < 0.001), chloramphenicol (AB-ZnO free: 0.00%; AB-ZnO: 58.33%; p < 0.001), and multi-drug resistance (AB-ZnO free: 64.71%; AB-ZnO: 100.00%; p < 0.001) on AB-ZnO farms. Isolates derived from ciprofloxacin-supplemented TBX, however showed only more frequent resistance to chloramphenicol (AB-ZnO free: 32.61%; AB-ZnO: 72.09%; p < 0.001) on AB-ZnO farms. Conversely, resistance to ciprofloxacin was higher in isolates originating from AB-ZnO free farms based on the analysis of isolates derived from apramycin-supplemented TBX (AB-ZnO free: 28.05%; AB-ZnO: 16.54%; p = 0.046) and cefotaxime-supplemented TBX (AB-ZnO free: 64.71%; AB-ZnO: 25.00%; p = 0.005) (Table 1).

Fig. 2
figure 2

Phenotypic resistance patterns of 351 Escherichia coli isolates. Isolates were obtained from a apramycin-supplemented, b ciprofloxacin-supplemented, and c cefotaxime-supplemented TBX

Associations between AMR profiles were seen when visualizing phenotypic resistance patterns (Fig. 2). Specific patterns of co-occurrence of AMR were frequently observed. For instance, E. coli isolates showing phenotypic resistance to apramycin were also likely to be resistant to tetracycline, streptomycin, trimethoprim, and chloramphenicol (Fig. 2a, c).

Interestingly, the opposite was observed for E. coli isolates derived from the ciprofloxacin- supplemented plates. These isolates were less likely to be resistant to apramycin, but similar to the other isolates they commonly displayed phenotypic resistance to tetracycline, streptomycin, trimethoprim, and chloramphenicol, in addition to nalidixic acid resistance (Fig. 2b).

Whole genome sequencing analysis of AmpC/ESBL E. coli isolates

In total, 64 AmpC β-lactamase-producing and ESBL E. coli isolates were confirmed by MALDI-TOF MS, of which 39 isolates derived from Farm8 (AB-ZnO farm), 21 isolates from Farm10 (AB-ZnO free farm) and 2 isolates from Farm2 (AB-ZnO free farm) and Farm6 (AB-ZnO farm), respectively. Interestingly, all of these isolates originated from farms that either used amoxicillin trihydrate together with ZnO or no treatment. No AmpC β-lactamase-producing and ESBL E. coli isolates were obtained from farms with sulphadiazine-trimethoprim together with ZnO treatment. Of the 64 isolates 44 (3 AmpC-positive and 41 ESBL positive) were selected for WGS, based on differences in their origin (farm, sample type, and selective plate), E. coli strain (identified via MALDI-TOF MS) and phenotypically expressed AMR.

Information on strain origin, serotypes, and sequence types, as well as ESBL genes present in all sequenced E. coli isolates are summarized in Additional file 1. The majority of β-lactamase resistance genes carried were the TEM-type gene blaTEM-1B_1, CTX-M-type gene blaCTX-M-1_1, and the OXA-type gene blaOXA-1_1.

No AmpC or ESBL resistance genes were detected for the phenotypically confirmed AmpC β-lactamase-producing isolates KB213, KB635, and KB643 which may be due to the set thresholds for minimum identity and coverage when screening for ARGs against the ResFinder database. The results of screening against the PointFinder database however, revealed mutations in the AmpC promotor associated with an AmpC-positive phenotype.

Apart from this, further analysis of ARGs showed that ARGs such as tet(A) (confers resistance to tetracycline), mdf(A) (confers multi-drug resistance), and aac(3)-IV (confers resistance to apramycin) were commonly detected in the E. coli isolates (Fig. 3). In total, twelve different MRGs were identified conferring resistance to zinc. All isolates carried at least one zinc resistance gene.

Fig. 3
figure 3

Identified resistance genes in AmpC / ESBL E. coli isolates. Genotypic and phenotypic results were additionally compared (see heatmap on the right). Sample abbreviations: FW0 faecal sample at day of weaning, FW1 faecal sample 1 week after weaning, FW2 faecal sample 2 week after weaning, SWD environmental swab sample of cleaned drinkers and feeders before contact with piglets, SWW environmental swab sample of cleaned walls and floors before contact with piglets

Interestingly, the colistin resistance gene mcr-1.1 was detected in five isolates from an AB-ZnO free farm. Phenotypic expression of colistin resistance was confirmed for all five strains (≥ 16 µg/mL) using Vitek2® and the AST-GN96 susceptibility test card, which tested for polymyxin B. In most cases, phenotypic and genotypic AMR profiles of E. coli isolates matched, except for resistance to quinolones which often could only be confirmed phenotypically.

AB-ZnO and AB-ZnO free farms clustered in different branches in the Euclidean distance analyses performed based on ARGs in the samples (Fig. 3). These clusters however, seemed to be less impacted by the sample type when compared to farm and ZnO and antimicrobial usage. Although distinct clusters were formed, no difference in the frequency of ARGs between AB-ZnO farms and AB-ZnO free farms was observed.

To further analyse which resistance and virulence genes were carried by plasmids, a network analysis was conducted. Based on the network analysis, IncI1_1_Alpha and IncFIB (AP001918) were the main carriers of blaCTX-M-1 and blaTEM-1B respectively (Fig. 4). Additionally, aac(3)-IV was frequently co-located together with blaCTX-M-1 on contigs identified as IncI1_1_Alpha. Similarly, tet(A) was frequently found on IncI1_1_Alpha and IncFIC(FII), but rarely together with other resistance genes. Overall, ARGs conferring resistance to aminoglycosides were predominantly found on plasmids. Other ARGs conferring resistance to trimethoprim, sulfonamide, tetracycline, and polymyxin were also identified, often together with ARGs of another resistance classes. The colistin resistance gene mcr-1.1 was only found on the IncX4 plasmid in five isolates, which did not carry any other virulence or resistance genes. Only two MRGs were found on plasmids, namely terC (confers resistance to tellurium) on IncH2A_1 and merA (confers resistance to mercury) on IncFII_1. Virulence factors on plasmids included an aerobactin-producing iucABCD-iutA operon located on the IncFIA_1 plasmid which was carried by a single isolate from Farm2 (KB213), the senB encoding for the Shigella enterotoxin 2 found in three isolates from Farm10 (KB1012, KB1015, and KB1038) carried by the Col156 plasmid, and the gene paa, encoding for adherence factors was only found in an isolate from Farm8 (KB835). All virulence factors identified were not co-located with other virulence factors or resistance genes. Virulence factors were only found in isolates originating from faecal samples.

Fig. 4
figure 4

Network analysis of identified plasmids, resistance genes, and virulence genes of AmpC/ESBL E. coli isolates. Each connection between plasmids and genes indicates that the gene was found on the plasmid. Connections between genes indicate their co-occurrence on the same plasmid per isolate. Connections between plasmids indicate their co-occurrence in the same isolate. Node placement is based on the Fruchterman-Reingold algorithm

In total, 30 unique integrons carrying ARGs were identified. Of these integrons, 21 were identified in isolates from Farm8, eight in isolates from Farm10, and one in the isolate from Farm2. Twenty-two integrons were found in isolates from faecal samples, while eight were found in isolates from environmental swab samples. Considering that three isolates which originated from faeces on Farm8 carried two AMR integrons, it can be summarized that of 19 of 31 isolates (61.29%) from faecal samples and 8 of 13 isolates (61.54%) from environmental samples carried at least one AMR integron. Overall, the ARGs aadA (confers aminoglycoside resistance), dfrA (confers trimethoprim resistance), cmlA (confers chloramphenicol resistance), and qacF (confers resistance to quaternary ammonium compounds) were predominant on all integrons. Additionally, the ARG ant(3’’)-Ia (confers aminoglycoside resistance) was also found on two integrons.

Eight of the AMR integrons carried more than three resistance genes, with the aminoglycoside ARGs aadA1 and aadA2, the chloramphenicol ARG cmlA1, and the trimethoprim ARG dfrA12 being predominant (Fig. 5).

Fig. 5
figure 5

Integrons with at least three ARGs, MRGs, or BRGs found in all samples. Potential integron classes are written in brackets for each gene in accordance with entries in the INTEGRALL database. Sample abbreviations: FW0 faecal sample at day of weaning, SWD environmental swab sample of cleaned drinkers and feeders in empty pens, SWW environmental swab sample of cleaned walls and floors in empty pens

All sequenced AmpC / ESBL E. coli isolates were compared using a pangenome analysis approach (Fig. 6). The pangenomic analysis resulted in 9,901 gene clusters with 218,672 genes, in the 44 E. coli genomes, which were identified using NCBI blastp. The core genome of all isolates consisted of 3,229 gene clusters with 147,907 genes. Additionally, Euclidean and Ward clustering resulted in several clusters independent of sample type. An apparent association of clusters with phylotypes of the E. coli isolates was observed. Phylotype A (30 isolates; 23 from AB-ZnO farms & 7 from AB-ZnO free farms) was predominantly identified, followed by B1 (11 isolates; 4 from AB-ZnO farms & 7 from AB-ZnO free farms) and E (3 isolates; 1 from a AB-ZnO farm & 2 from AB-ZnO free farms).

Fig. 6
figure 6

Pangenome analysis for 44 AmpC / ESBL E. coli isolates. Gene clusters are indicated as opaque coloured elements and ordered based on their presence and absence, using the Euclidean distance and Ward clustering (see also inner dendrogram). Samples are hierarchically clustered based on the frequency of gene clusters (see also outer dendrogram). Sample abbreviations: FW0 faecal sample at day of weaning, FW1 faecal sample 1 week after weaning, FW2 faecal sample 2 week after weaning, SWD environmental swab sample of cleaned drinkers and feeders in empty pens, SWW environmental swab sample of cleaned walls and floors in empty pens

Discussion

In this study, E. coli isolates derived from 10 Irish pig farms, 5 farms that regularly used either amoxicillin or sulphadiazine-trimethoprim, together with ZnO, for 2-weeks as prophylactic treatment during weaning, and 5 farms that did not use any prophylaxis during weaning, were compared for the presence of AMR E. coli.

E. coli counts were first compared between farms by culturing samples on non-supplemented and antimicrobial-supplemented TBX agar. E. coli counts obtained from non-supplemented TBX were almost identical between AB-ZnO farms and AB-ZnO free farms. Furthermore, no decrease in total E. coli counts in faecal samples was observed 2 weeks post-weaning. Similar results were observed in a study of Ciesinski, et al. [15], who also did not report any differences in E. coli counts in faeces of ZnO-treated and non-treated pigs. When comparing E. coli counts obtained from apramycin-containing TBX, numerical differences between faecal samples originating from AB-ZnO farms and AB-ZnO free farms were observed from 1 week post-weaning. However, these differences were not statistically different, possibly due to the limited number of farms and samples taken. Therefore, larger number of farms and samples may be necessary in future studies to further investigate these differences. As discussed by Fairbrother, et al. [8] and Broom, et al. [25], administration of ZnO-supplemented feed may not result in an overall reduction of E. coli, but may rather result in a shift of the diverse E. coli community to more resistant strains, as reported by Ghazisaeedi, et al. [26] and Bednorz, et al. [14]. Comparison of AMR profiles of E. coli isolates confirmed the observations made by Ghazisaeedi, et al. [26] and Bednorz, et al. [14], showing a higher frequency of E. coli isolates resistant to apramycin, trimethoprim, tetracycline, streptomycin, chloramphenicol, and with multi-drug resistance on AB-ZnO farms. Interestingly, resistance to ciprofloxacin was more frequently observed for isolates originating from AB-ZnO free farms, although the reasons for this were not clear.

When analysing phenotypic resistance profiles of the 351 isolates, it was observed that certain resistance patterns occurred more frequently, independent of the selective pressure mediated by culturing isolates on supplemented TBX. Commonly observed resistance patterns across all supplemented agars included phenotypic resistance to tetracycline, streptomycin, trimethoprim, and chloramphenicol. Similar resistance patterns were observed by Ciesinski, et al. [15], who investigated the effects of zinc oxide (2103 mg/kg diet) use on MDR E. coli in the porcine gut. In their study, they described three frequently observed resistance patterns: (1) Ampicillin-Streptomycin-Tetracycline; (2) Streptomycin-Sulfamethoxazole/Trimethoprim-Tetracycline; and (3) Ampicillin-Streptomycin-Sulfamethoxazole/Trimethoprim-Tetracycline. Although they used different antimicrobials for phenotypic AMR determination, the observed resistance patterns were similar.

These results may indicate a selective adaptive response of E. coli to treatment on farms using sulphadiazine-trimethoprim and ZnO, but also non-selective adaption possibly due to cross-resistance or co-location of ARGs on MGEs, as previously reported by a study of Rovira, et al. [27]. To further investigate possible co-selection events, 44 AmpC/ESBL E. coli isolates were sequenced using a whole-genome sequencing approach.

Analysis of resistance profiles of AmpC / ESBL E. coli isolates showed that the ARGs aac(3)-IV (N-acetyltransferase AAC(3)-IV) and tet(A) (efflux pump) were two of the most frequently observed ARGs in all AmpC / ESBL E. coli isolates. Both of these ARGs were previously reported to be frequently found on plasmids [28, 29]. This was also confirmed by the network analysis in the present study, showing the frequent occurrence of tet(A) on IncI1_1_Alpha and IncFIC(FII). The ARG aac(3)-IV was also frequently found on IncI1_1_Alpha and was commonly co-located with blaCTX-M-1 on the same plasmids, carried by the majority of the AmpC / ESBL E. coli isolates. The occurrence of IncI1_1_Alpha-blaCTX-M-1 plasmids in livestock production has been regularly reported in Europe and is of great concern for public health [30,31,32]. These plasmids are highly transmissible between Enterobacteriaceae, including Salmonella enterica serovar Typhimurium, and have been reported in food-producing animals but also in clinical settings, with evidence that this ARG may disseminate via E. coli contamination in the food chain [33,34,35,36]. A study by Hiley, et al. [36] reported the occurrence of the plasmid IncI1 in S. Typhimurium isolates associated with a food poisoning outbreak in 2015. Similarly, blaCTX-M-1 has been reported in a Klebsiella pneumoniae outbreak in Spain within an intensive care unit by Mena and colleagues [37]. However, to the best of the author’s knowledge, no study yet reported a clinical disease outbreak caused by Enterobacteriaceae carrying IncI1_1_Alpha-blaCTX-M-1.

On Farm10 five mcr-1.1-positive ESBL E. coli isolates were obtained which showed phenotypic resistance to polymyxin B. Network analysis revealed that mcr-1.1 was located on the plasmid IncX4, which has been reported in previous clinical and environmental studies [38,39,40,41]. To the best of the authors’ knowledge this is the first study reporting mcr-1.1-positive E. coli isolates on Irish commercial pig farms. Co-occurrence of the IncX4-mcr-1.1 plasmid with the IncI1_1_Alpha plasmid in the same bacterial isolate was observed for two isolates.

Integron analysis also showed frequent co-occurrence of ARGs such as dfrA, aadA, cmlA, and qacF conferring trimethoprim, aminoglycoside, chloramphenicol, and resistance to quaternary ammonium compounds respectively. These ARGs were reported before on pig farms by Mencía-Ares, et al. [42], mainly associated with class 1 integrons and may partially explain the frequent occurrence of the previously discussed resistance patterns. The majority of these integrons were found in isolates from Farm8. Thus, the use of antimicrobials such as amoxicillin (β-lactam) and sulphadiazine-trimethoprim (sulphonamide and trimethoprim) may facilitate the co-selection of ARGs conferring resistance to these agents, with ARGs conferring resistance to aminoglycoside such as apramycin and streptomycin, chloramphenicol, quaternary ammonium compounds, and tetracyclines. Although the number of isolates may not allow clear conclusions to be drawn, co-selection and co-location of ARGs of different antimicrobial classes on MGEs may explain the increased phenotypic resistance observed to these antimicrobial classes, including multi-drug resistance. Furthermore, carriage of multiple ARG-plasmids may additionally allow the persistence of AMR isolates, such as mcr-1.1-positive E. coli strains, on farms due to co-selection and fitness advantages against other niche competitors during antimicrobial and ZnO usage.

Apart from that, it was observed that all sequenced E. coli strains carried at least one MRG conferring resistance to zinc. No zinc-resistance genes were found on MGEs and it remains unclear in the current study settings, if ZnO plays a crucial role in comparison to AMU in terms of AMR selection or the observation of the increased numbers of MDR E. coli observed on AB-ZnO farms. Ghazisaeedi, et al. [26] and Bednorz, et al. [14] reported increases in MDR E. coli associated with ZnO supplementation. While Bednorz, et al. [14] concluded in their study that increases in MDR E. coli could be due to co-selection under selective pressure and enhanced bacterial conjugation rates [43, 44], Ghazisaeedi, et al. [26] suggested that co-selection may not allow higher tolerance to zinc, but MDR isolates may be more resistant to environmental stresses such as zinc exposure per se.

Phenotypic differences between farms were also reflected in the genotypic results, although phenotypically observed resistance could not always be confirmed genotypically. Mutations in gyrA (encodes DNA gyrase) and parC (encodes DNA topoisomerase IV) often confer resistance to quinolones, as well as qnr or aac(6')-Ib-cr resistance genes. These mutations are not included in the ResFinder database, which would explain why resistance to quinolones was frequently confirmed phenotypically, but not genotypically in many isolates [45]. However, additional screening of these 13 isolates, showing phenotypic resistance to quinolones, against the PointFinder database confirmed mutations in gyrA (in 12 cases), gyrA and parC (in 4 cases), and gyrA, parC, and parE (in 2 cases). Moreover, the three isolates which were confirmed phenotypically as AmpC β-lactamase-producers, but for which no AmpC β-lactamase encoding resistance gene was found, could be further confirmed using the PointFinder database.

Distinct clusters of E. coli isolates from AB-ZnO and AB-ZnO free farms were observed in the ARG heatmap and the pangenomic analysis. Clusters in the ARG heatmap were more associated with differences in ARG carriage, while clusters in the pangenome analysis seemed to be more associated with E. coli phylotypes. The predominant phylotype observed on all pig farms was phylotype A, and phylotype A and B1 in the case of AB-ZnO free farms. These phylogenetic groups have been reported before in faecal samples of pigs and are known to represent mainly commensal, but also intestinal, pathogenic strains [46]. Additionally, three strains belonged to phylogenetic group E, of which E. coli O157:H7 is a well-known member [47, 48]. Two of these three strains originated from Farm10 and belonged to the serotype O167:H4. The other strain serotype could not be identified. The serogroup O167 has been reported before as enterotoxigenic or enteroinvasive E. coli. Virulence factors associated with these pathogroups were also found on MGEs, such as paa and senB, however not in the same isolates.

The results of the pangenome analysis showed that the core genome of all sequenced isolates consisted of 3,229 gene clusters, while the total pangenome consisted of 9,901 gene clusters. This is in accordance with results of Cummins, et al. [49], who reported that pangenomes vary between E. coli sequence types. Thus, they reported that E. coli core genome sizes can range from 3,066 to 4,058 clusters, with a pangenome size range from 6,825 to 27,634 gene clusters.

Future studies in this field should include in-depth analysis methods such as long-read sequencing or combined hybrid approaches, which with current costs, do not permit sequencing the number of isolates included in this study. As shown in this study, the information gained by using short-sequencing methods is sufficient to identify potential pathogens, their plasmids, and resistance genes, and demonstrates the benefits of combining culture-dependent and culture-independent analyses to provide a more complete picture.

Conclusions

In conclusion, the results of this study showed prophylactic use of antimicrobials and ZnO for 14 days post-weaning on farms favoured the selection and development of AMR and MDR E. coli. Frequently observed AMR patterns and higher prevalence of AMRs on AB-ZnO farms could be partially explained by co-location of ARGs on MGEs such as integrons or plasmids, which are potentially more frequently transferred under ZnO/antimicrobial exposure. The study also identified the first mcr-1-encoding E. coli isolated in Irish pig production. This study demonstrated the usefulness of combining microbiological and sequencing methods to evaluate the impact of antimicrobial and heavy metal removal on AMR on pig farms.

Methods

Farm sampling

Samples were collected from 10 pig farms in Ireland between September 2020 and April 2021. Five farms regularly used prophylactic treatments with antimicrobials and ZnO for the first 2 weeks post-weaning (AB-ZnO farms) for more than 3 years and five farms did not use any such medication for the 3 years prior to sample collection (AB-ZnO free farms). Of the AB-ZnO farms Farm4, Farm6, and Farm8 used amoxicillin trihydrate (Stabox [15 mg/kg], Virbac, France) and Farm1 and Farm3 used sulphadiazine-trimethoprim (Sulfoprim 15% [15 mg/kg], Univet Limited, Ireland), all in combination with ZnO (3000 ppm, i.e. 3 g ZnO/kg feed). A total of five different samples were taken on each farm, including two environmental and three faecal samples. Piglets on each farm were weaned at an age of approximately 28 days and allocated to nursery pens. Immediately before entry of piglets to the pens, environmental samples of cleaned walls, floors, feeders, and drinkers were taken using a sponge-stick soaked in 10 mL neutralizing buffer (3 M, Ireland). One sponge was used to swab an area of 1 m2 of cleaned walls and floors, and another sponge was used to swab cleaned drinkers and feeders. Freshly voided faecal droppings of piglets were collected on the day of weaning (day 0), as well as on days 7 and 14 post-weaning. Three faecal droppings per batch were pooled together in a sterile plastic container and mixed with a sterile micro-spatula. All samples were transported to the laboratory on ice in a polystyrene cooling box and processed within 24 h.

E. coli isolation

Faecal samples were mixed with a sterile micro-spatula before processing in the laboratory. From each sample, 4 g of faeces were transferred into a sterile 150 mL container and mixed with 36 mL maximum recovery diluent (MRD, Oxoid, United Kingdom) to obtain a 1:10 dilution. After vortexing, the suspended faecal samples were transferred into filter stomacher bags (Grade products, United Kingdom) and were homogenized for 60 s with a Star Blender LB 400 (VWR, Ireland).

For the swab samples, 20 mL sterile Tryptone Soy Broth (TSB, Oxoid, United Kingdom) was added directly to the swab bag and the sponge was squeezed several times. Swabs were incubated for 1 h at 37 °C for bacterial enrichment.

Faecal and environmental sample suspensions were serially diluted and plated in duplicate on non-supplemented and supplemented Tryptone Bile X-glucuronide agar (TBX, Merck KGaA, Germany) plates containing either no antibiotics or 16 mg/L apramycin sulfate salt (Sigma-Aldrich, Germany, according to Yates, et al. [50], 1 mg/L ciprofloxacin (Thermo Scientific, Ireland, or 4 mg/L cefotaxime sodium salt (Sigma-Aldrich, Germany), in accordance with the Clinical and Laboratory Standards Institute; CLSI M100-ED32:2022). Antimicrobials were dissolved in sterilized deionised water and filtered using sterile 0.22 µm syringe filters (GVS, USA), before addition to media. After samples were plated they were incubated for 24 h at 44 °C. It is important to mention that enumeration of colony counts after incubation was only conducted for faecal samples, due to the enrichment step for the environmental samples. Up to 5 colonies per sample were selected and transferred into cryogenic tubes (Technical Service Consultant Ltd., UK), which were stored at −80 °C.

Phenotypic antimicrobial resistance determination

In total, 361 isolates were analysed using the Kirby-Bauer disc diffusion susceptibility test [51]. Isolates were streaked out on Tryptone Soya Agar (TSA, Oxoid, United Kingdom) and incubated for 24 h at 37 °C. Following this, one colony of each isolate was transferred to 10 mL TSB and incubated for 16 h at 37 °C to prepare overnight cultures (ONCs). ONCs were serial diluted in MRD. Suspensions were spread on Mueller–Hinton Agar (MHA, Sigma-Aldrich, Germany) followed by the application of antibiotic discs on the plates. Nine antimicrobials were tested, including apramycin (15 µg), trimethoprim (5 µg), nalidixic acid (30 µg), tetracycline (30 µg), meropenem (10 µg), cefotaxime (5 µg), ciprofloxacin (5 µg), streptomycin (10 µg), and chloramphenicol (30 µg) (Oxoid, UK). The MHA plates were then incubated for 24 h at 37 °C. Susceptibility tests were conducted in duplicate and E. coli strain ATCC® 25922 was used as a reference. Inhibition zones were measured and evaluated in accordance with clinical breakpoints published by the CLSI and the European Committee on Antimicrobial Susceptibility Testing (EUCAST). The susceptibility test results of 10 of 361 isolates were not consistent with the results of the selective plating of the faecal and environmental samples and were excluded from further analysis, thus 351 isolates remained. Where a colistin resistance gene was detected via WGS analysis, phenotypic expression of polymyxin B resistance was assessed using the Vitek2® and the AST-GN96 susceptibility test card (BioMérieux, United Kingdom).

Identification and confirmation of AmpC/ESBL E. coli

Isolates showing phenotypic resistance to cefotaxime were screened for AmpC and extended-spectrum β-lactamases (ESBLs) using the MASTDISCS® AmpC ESBL detection set (Product code: D68C, Mast Group Ltd., UK).

Isolates that were confirmed as ESBL and AmpC β-lactamase producers were analysed via MALDI-TOF MS (Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry) at DAFM Laboratories Backweston (Department of Agriculture, Food and the Marine, Ireland) to confirm that they were E. coli. Briefly, AmpC and ESBL isolates were streaked out on TSA and incubated at 37 °C for 24 h. One colony of each isolate was transferred into 2 mL screw cap micro tubes and suspended in 300 µl ultrapure water (18.2 MΩxcm). Subsequently, 900 µl pure ethanol (Merck, Ireland) were added and vortexed. The suspensions were centrifuged at 13,000 rpm for 2 min, the supernatant removed, and the pellet allowed to dry for 5–10 min. Pellets were re-suspended in 25 µl of 70% formic acid (Merck, Ireland) before 25 µl acetonitrile (Sigma-Aldrich, Ireland) was added. The suspensions were centrifuged at the same settings as previously. For each isolate, 1 µl of supernatant was transferred onto different spots of a MALDI target plate. After the spots dried, 1 µl of a-Cyano-4-hydroxycinnamic acid matrix solution (HCCA, Bruker Daltronics GmbH, Germany) was added. After another drying step, isolates were analysed via MALDI-TOF using the Bruker MALDI Biotyper IVD (version 10) and the FlexControl software (version 3.4) (Bruker Daltronics GmbH, Germany). The US IVD Bacterial Test Standard (BTS, Bruker Daltronics GmbH, Germany) was used for validation and as a quality control for the MALDI-TOF analysis. Isolates were analysed using the ‘MBT_AutoX’ AutoXecute acquisition method with a mass-to-charge ratio (m/z) range of 1.06–20.14 kDA. Obtained protein spectra were aligned against the BDAL Bruker database. Logarithmic alignment scores and consistencies were calculated by the MBT Compass Explorer software (version 4.1; Build 100), ranging from 0.00 to 3.00 and A(+ + +)−C(−) respectively. Only results with a score greater than 2.00 and a consistency of A were accepted, resulting in 64 confirmed E. coli isolates, of which 44 isolates were selected based on their origin, taxonomy, and resistance profiles for further analysis via whole-genome sequencing.

Whole genome sequencing analysis of AmpC/ESBL E. coli isolates

DNA of confirmed AmpC / ESBL E. coli isolates was extracted using the DNeasy® UltraClean® Microbial kit (Qiagen, Germany). DNA extracts were quantified using the NanodropTM 1000 Spectrophotometer and the Qubit 4.0 Fluorometer (Invitrogen, ThermoFisher Scientific, United Kingdom). DNA libraries of extracts were prepared using the Illumina DNA prep kit in accordance with the manufacturer’s instructions, followed by 2 × 150 bp sequencing with P2 reagent on the Illumina Nextseq2000 platform. Obtained raw reads were cleaned using Trimmomatic (v0.38) [52] and the quality of trimmed raw reads was assessed using FastQC (v0.11.8) and MultiQC (v1.9) [53, 54]. The taxonomic assignment of E. coli strains was confirmed using Kraken 2 (v2.1.1) and the Kraken 2 PlusPF database (last update: 2021-May-17). Subsequently, reads were assembled to scaffolds using SPAdes (v3.15.3) and the –isolate tag [55]. QUAST and MultiQC were used to assess the quality of assembled scaffolds [53, 56]. Afterwards, scaffolds were examined for antimicrobial resistance genes and heavy metal resistance genes via ABRicate (v1.0.1; https://github.com/tseemann/abricate) using the ResFinder database (last update: 2021-Mar-27) with –minid 80 –mincov 75, and Bacmet-Scan (v1.0) with default settings and the BacMet2 database containing experimentally confirmed resistance genes [57]. Additionally, scaffolds were screened against the PointFinder database (last update: 2022-Apr-22) with a minimum threshold for identity of 70% and a minimum length of 60% via the CGE online platform (Center for Genomic Epidemiology; http://www.genomicepidemiology.org/) to account for resistance caused by point mutations [58]. E. coli serotypes were determined using ABRicate and the EcOH database (last update: 2021-Mar-27) with the same identity and coverage thresholds. Multi-locus sequence typing was conducted using mlst (v2.18.0, https://github.com/tseemann/mlst) [59]. To further examine for MGEs and their association to specific ARGs and MRGs, assembled scaffolds were separated into chromosomal and plasmidic scaffolds using Platon (v1.6) [60]. Plasmids based on plasmidic scaffolds were identified using PlasmidFinder (last update: 2021-Mar-27) and were analysed for ARGs, BRGs and MRGs, and virulence factors using ABRicate together with the ResFinder and VFDB database (last update: 2021-Mar-27) with the previous mentioned thresholds, as well as Bacmet-Scan (v1.0). Only plasmidic scaffolds which were identified by PlasmidFinder and carried resistance or virulence genes were included in the subsequent network analysis. Scaffolds obtained by SPAdes were screened for integrons using IntegronFinder (v2.0) [61]. Resistance genes on integrons were identified using ABRicate, BacMet2, and the mobilome analysis pipeline created by José F. Cobo-Díaz (https://github.com/JoseCoboDiaz/ARG-contig_mobilome_analysis). The pangenome of all 44 isolates was reconstructed using anvi’o (v7.1).

Data analysis and visualization

All graphs except for the pangenome analysis were visualized using R v4.0.2 (2020-06-22) and RStudio v2021.09.2 [62]. E. coli concentrations were visualized using the ggplot2 package [63] and were compared using the Mann-Witney U test via the rstatix package (v0.7.0) [64]. In cases where no colonies were obtained at the lowest dilution, a value of 0 log10(CFU/mL) was set, representing 1 cfu in an undiluted sample. Upset plots for visualizing phenotypic resistance patterns of E. coli isolates (Fig. 2) were done using the ComplexUpset package [65]. Statistical comparison of phenotypic resistance patterns was performed by conducting the Chi-square test via Minitab 17 (v17.1.0). Heatmaps showing the presence and absence of ARGs of AmpC/ESBL E. coli isolates were visualized using the pheatmap package [66]. Network analysis of plasmids, virulence factors, and resistance genes was conducted using the igraph package (v1.3.0) and the Fruchterman-Reingold algorithm for node placement [67, 68]. Integrons were visualized using the gggenes package (v0.4.1) [69]. The pangenome analysis of the 44 AmpC/ESBL E. coli isolates was performed using anvi’o (v7.1) in accordance with the workflow for microbial pangenomics (https://merenlab.org/2016/11/08/pangenomics-v2/) [70]. Briefly, assembled scaffold.fasta files were first converted into anvi’o contig databases using the ‘anvi-gen-contigs-database’ program. Afterwards, genes in scaffolds were identified by using Prodigal to screen for open reading frames, which were further annotated using the NCBI’s Clusters of Orthologous Groups database (‘anvi-run-ncbi-cogs’ program) [71, 72]. Genes were annotated against four provided HMM profiles of anvi’o using hidden Markov models (‘anvi-run-hmms’ program). The pangenome was computed by (1) calculating the amino acid sequence similarities and comparison across all genomes via NCBI blastp, (2) eliminating weak matches between amino acid sequences using minbit heuristics of 0.5 [73], and (3) identifying clusters using the MCL algorithm [74] (‘anvi-pan-genome’ program). The pangenome was visualized and edited using the anvi’o interactive interface and the program ‘anvi-display-pan’. Figures were arranged and modified using Affinity Designer (v1.8.5.703, Serif).

Availability of data and materials

Raw sequences from the whole-genome sequencing analysis of the AmpC / ESBL E. coli isolates are publicly accessible under BioProject: PRJNA901700.

References

  1. Destoumieux-Garzón D, Mavingui P, Boetsch G, Boissier J, Darriet F, Duboz P, et al. The one health concept: 10 years old and a long road ahead. Front Vet Sci. 2018;5:14.

    Article  PubMed  PubMed Central  Google Scholar 

  2. The Food and Agriculture Organization: drivers, dynamics and epidemiology of antimicrobial resistance in animal production. 2016. https://www.fao.org/3/i6209e/i6209e.pdf. Accessed 26 Jul 2022.

  3. Peng S, Zheng H, Herrero-Fresno A, Olsen JE, Dalsgaard A, Ding Z. Co-occurrence of antimicrobial and metal resistance genes in pig feces and agricultural fields fertilized with slurry. Sci Total Environ. 2021;792:148259.

    Article  CAS  PubMed  Google Scholar 

  4. Cheng G, Ning J, Ahmed S, Huang J, Ullah R, An B, et al. Selection and dissemination of antimicrobial resistance in Agri-food production. Antimicrob Resist Infect Control. 2019;8:158.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ramos S, Silva V, Dapkevicius MLE, Caniça M, Tejedor-Junco MT, Igrejas G, et al. Escherichia coli as commensal and pathogenic bacteria among food-producing animals: health implications of extended spectrum β-lactamase (ESBL) production. Animals. 2020;10:2239.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Russo TA, Johnson JR. Medical and economic impact of extraintestinal infections due to Escherichia coli: focus on an increasingly important endemic problem. Microbes Infect. 2003;5:449–56.

    Article  PubMed  Google Scholar 

  7. Luppi A. Swine enteric colibacillosis: diagnosis, therapy and antimicrobial resistance. Porc Health Manag. 2017;3:16.

    Article  Google Scholar 

  8. Fairbrother JM, Nadeau E, Gyles CL. Escherichia coli in postweaning diarrhea in pigs: an update on bacterial types, pathogenesis, and prevention strategies. Anim Health Res Rev. 2005;6:17–39.

    Article  CAS  PubMed  Google Scholar 

  9. Bonetti A, Tugnoli B, Piva A, Grilli E. Towards zero zinc oxide: feeding strategies to manage post-weaning diarrhea in piglets. Animals. 2021;11:642.

    Article  PubMed  PubMed Central  Google Scholar 

  10. European Commission: Veterinary medicines: new rules to promote animal health and fight antimicrobial resistance now apply. 2022. https://ec.europa.eu/commission/presscorner/detail/en/ip_22_663. Accessed 26 Jul 2022.

  11. European Food Safety Authority. The European Union summary report on antimicrobial resistance in zoonotic and indicator bacteria from humans, animals and food in 2019–2020. EFSA J. 2022;20:e07209.

    Article  CAS  PubMed Central  Google Scholar 

  12. Caruso G. Antibiotic resistance in Escherichia coli from farm livestock and related analytical methods: a review. J AOAC Int. 2018;101:916–22.

    Article  CAS  PubMed  Google Scholar 

  13. Anjum MF, Schmitt H, Börjesson S, Berendonk TU, Donner E, Stehling EG, et al. The potential of using E. coli as an indicator for the surveillance of antimicrobial resistance (AMR) in the environment. Curr Opin Microbiol. 2021;64:152–8.

    Article  CAS  PubMed  Google Scholar 

  14. Bednorz C, Oelgeschläger K, Kinnemann B, Hartmann S, Neumann K, Pieper R, et al. The broader context of antibiotic resistance: zinc feed supplementation of piglets increases the proportion of multi-resistant Escherichia coli in vivo. Int J Med Microbiol. 2013;303:396–403.

    Article  CAS  PubMed  Google Scholar 

  15. Ciesinski L, Guenther S, Pieper R, Kalisch M, Bednorz C, Wieler LH. High dietary zinc feeding promotes persistence of multi-resistant E. coli in the swine gut. PLoS One. 2018;13:e0191660.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Liu P, Pieper R, Rieger J, Vahjen W, Davin R, Plendl J, et al. Effect of dietary zinc oxide on morphological characteristics, mucin composition and gene expression in the colon of weaned piglets. PLoS ONE. 2014;9:e91091.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Liu P, Pieper R, Tedin L, Martin L, Meyer W, Rieger J, et al. Effect of dietary zinc oxide on jejunal morphological and immunological characteristics in weaned piglets. J Anim Sci. 2014;92:5009–18.

    Article  CAS  PubMed  Google Scholar 

  18. Karweina D, Kreuzer-Redmer S, Müller U, Franken T, Pieper R, Baron U, et al. The zinc concentration in the diet and the length of the feeding period affect the methylation status of the ZIP4 zinc transporter gene in piglets. PLoS ONE. 2015;10:e0143098.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bondzio A, Pieper R, Gabler C, Weise C, Schulze P, Zentek J, et al. Feeding low or pharmacological concentrations of zinc oxide changes the hepatic proteome profiles in weaned piglets. PLoS ONE. 2013;8:e81202.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Vahjen W, Pietruszyńska D, Starke IC, Zentek J. High dietary zinc supplementation increases the occurrence of tetracycline and sulfonamide resistance genes in the intestine of weaned pigs. Gut Pathog. 2015;7:23.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Silver S, Phung LT. Bacterial heavy metal resistance: new surprises. Annu Rev Microbiol. 1996;50:753–89.

    Article  CAS  PubMed  Google Scholar 

  22. Byrne N, O’Neill L, Dίaz JAC, Manzanilla EG, Vale AP, Leonard FC. Antimicrobial resistance in Escherichia coli isolated from on-farm and conventional hatching broiler farms in Ireland. Ir Vet J. 2022;75:7.

    Article  PubMed  PubMed Central  Google Scholar 

  23. European Commission: Regulation (EU) 2019/6 of the European Parliament and of the Council of 11 December 2018 on veterinary medicinal products and repealing Directive 2001/82/EC. 2019. https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32019R0006&from=EN%0Ahttps://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32019R0006&qid=1552299700950&from=EN. Accessed 13 Jan 2023.

  24. European Commission: Regulation (EU) 2019/4 of the European Parliament and of the Council of 11 December 2018 on the manufacture, placing on the market and use of medicated feed, amending Regulation (EC) No 183/2005 of the European Parliament and of the Council and repealing Council Directive 90/167/EEC. 2019. https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32019R0004&from=EN. Accessed 13 Jan 2023.

  25. Broom LJ, Miller HM, Kerr KG, Knapp JS. Effects of zinc oxide and Enterococcus faecium SF68 dietary supplementation on the performance, intestinal microbiota and immune status of weaned piglets. Res Vet Sci. 2006;80:45–54.

    Article  CAS  PubMed  Google Scholar 

  26. Ghazisaeedi F, Ciesinski L, Bednorz C, Johanns V, Pieper L, Tedin K, et al. Phenotypic zinc resistance does not correlate with antimicrobial multi-resistance in fecal E. coli isolates of piglets. Gut Pathog. 2020;12:4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Rovira P, McAllister T, Lakin SM, Cook SR, Doster E, Noyes NR, et al. Characterization of the microbial resistome in conventional and “raised without antibiotics” beef and dairy production systems. Front Microbiol. 2019;10:1980.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Roberts MC. Tetracycline resistance determinants: mechanisms of action, regulation of expression, genetic mobility, and distribution. FEMS Microbiol Rev. 1996;19:1–24.

    Article  CAS  PubMed  Google Scholar 

  29. Davies J, O’Connor S. Enzymatic modification of aminoglycoside antibiotics: 3-N-acetyltransferase with broad specificity that determines resistance to the novel aminoglycoside apramycin. Antimicrob Agents Chemother. 1978;14:69–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zurfluh K, Wang J, Klumpp J, Nüesch-Inderbinen M, Fanning S, Stephan R. Vertical transmission of highly similar blaCTX-M-1-harbouring IncI1 plasmids in Escherichia coli with different MLST types in the poultry production pyramid. Front Microbiol. 2014;5:519.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Dahmen S, Haenni M, Madec J-Y. IncI1/ST3 plasmids contribute to the dissemination of the blaCTX-M-1 gene in Escherichia coli from several animal species in France. J Antimicrob Chemother. 2012;67:3011–2.

    Article  CAS  PubMed  Google Scholar 

  32. Börjesson S, Bengtsson B, Jernberg C, Englund S. Spread of extended-spectrum beta-lactamase producing Escherichia coli isolates in Swedish broilers mediated by an incI plasmid carrying blaCTX-M-1. Acta Vet Scand. 2013;55:3.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Mo SS, Telke AA, Osei KO, Sekse C, Slettemeås JS, Urdahl AM, et al. blaCTX–M–1/IncI1-Iγ plasmids circulating in Escherichia coli from Norwegian broiler production are related, but distinguishable. Front Microbiol. 2020;11:333.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Leverstein-van Hall MA, Dierikx CM, Stuart JC, Voets GM, van den Munckhof MP, van Essen-Zandbergen A, et al. Dutch patients, retail chicken meat and poultry share the same ESBL genes, plasmids and strains. Clin Microbiol Infect. 2011;17:873–80.

    Article  CAS  PubMed  Google Scholar 

  35. Abraham S, Kirkwood RN, Laird T, Saputra S, Mitchell T, Singh M, et al. Dissemination and persistence of extended-spectrum cephalosporin-resistance encoding IncI1-blaCTXM-1 plasmid among Escherichia coli in pigs. ISME J. 2018;12:2352–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hiley L, Graham RMA, Jennison AV. Characterisation of IncI1 plasmids associated with change of phage type in isolates of Salmonella enterica serovar Typhimurium. BMC Microbiol. 2021;21:92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Mena A, Plasencia V, García L, Hidalgo O, Ayestarán JI, Alberti S, et al. Characterization of a large outbreak by CTX-M-1-producing Klebsiella pneumoniae and mechanisms leading to in vivo carbapenem resistance development. J Clin Microbiol. 2006;44:2831–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Hassan J, Eddine RZ, Mann D, Li S, Deng X, Saoud IP, et al. The mobile colistin resistance gene, mcr-1.1, is carried on IncX4 plasmids in multidrug resistant E. coli isolated from rainbow trout aquaculture. Microorganisms. 2020;8:1636.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Majewski P, Gutowska A, Smith DGE, Hauschild T, Majewska P, Hryszko T, et al. Plasmid mediated mcr-1.1 colistin-resistance in clinical extraintestinal Escherichia coli strains isolated in Poland. Front Microbiol. 2021;12:547020.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Furlan JPR, Lopes R, Ramos MS, dos Santos LDR, da Silva RR, Savazzi EA, et al. Colistin-resistant mcr-1-positive Escherichia coli ST1775-H137 co-harboring blaCTX-M-2 and blaCMY-2 recovered from an urban stream. Infect Genet Evol. 2021;96:105156.

    Article  CAS  PubMed  Google Scholar 

  41. Girardello R, Piroupo CM, Martins J, Maffucci MH, Cury AP, Franco MRG, et al. Genomic characterization of mcr-1.1-producing Escherichia coli recovered from human infections in São Paulo. Brazil Front Microbiol. 2021;12:663414.

    Article  PubMed  Google Scholar 

  42. Mencía-Ares O, Cabrera-Rubio R, Cobo-Díaz JF, Álvarez-Ordóñez A, Gómez-García M, Puente H, et al. Antimicrobial use and production system shape the fecal, environmental, and slurry resistomes of pig farms. Microbiome. 2020;8:164.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Ou JT. Effect of Zn2+ on bacterial conjugation: increase in ability of F- cells to form mating pairs. J Bacteriol. 1973;115:648–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Ekhlas D, Soro AB, Leonard FC, Manzanilla EG, Burgess CM. Examining the impact of zinc on horizontal gene transfer in Enterobacterales. Sci Rep. 2022;12:20503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Johnning A, Kristiansson E, Fick J, Weijdegård B, Larsson DGJ. Resistance mutations in gyrA and parC are common in Escherichia communities of both fluoroquinolone-polluted and uncontaminated aquatic environments. Front Microbiol. 2015;6:1355.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Carlos C, Pires MM, Stoppe NC, Hachich EM, Sato MIZ, Gomes TAT, et al. Escherichia coli phylogenetic group determination and its application in the identification of the major animal source of fecal contamination. BMC Microbiol. 2010;10:161.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  48. Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8:207–17.

    Article  CAS  PubMed  Google Scholar 

  49. Cummins EA, Hall RJ, Connor C, McInerney JO, McNally A. Pangenome evolution in Escherichia coli is sequence type, not phylogroup, specific. bioRxiv. 2022. https://0-doi-org.brum.beds.ac.uk/10.1101/2022.05.20.492775.

    Article  Google Scholar 

  50. Yates CM, Pearce MC, Woolhouse MEJ, Amyes SGB. High frequency transfer and horizontal spread of apramycin resistance in calf faecal Escherichia coli. J Antimicrob Chemother. 2004;54:534–7.

    Article  CAS  PubMed  Google Scholar 

  51. Bauer AW, Kirby WM, Sherris JC, Turck M. Antibiotic susceptibility testing by a standardized single disk method. Am J Clin Pathol. 1966;45:493–6.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed 14 Nov 2022.

  55. Prjibelski A, Antipov D, Meleshko D, Lapidus A, Korobeynikov A. Using SPAdes De Novo Assembler. Curr Protoc Bioinformatics. 2020;70:e102.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DGJ. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2013;42:D737–43.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Jolley KA, Maiden MCJ. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Schwengers O, Barth P, Falgenhauer L, Hain T, Chakraborty T, Goesmann A. Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores. Microb Genom. 2020;6:mgen000398.

    PubMed  PubMed Central  Google Scholar 

  61. Néron B, Littner E, Haudiquet M, Perrin A, Cury J, Rocha EPC. IntegronFinder 2.0: identification and analysis of integrons across bacteria, with a focus on antibiotic resistance in Klebsiella. Microorganisms. 2022;10:700.

    Article  PubMed  PubMed Central  Google Scholar 

  62. R Core Team. R: a language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria. 2020. https://www.R-project.org/. Accessed 14 Nov 2022.

  63. Wickham H. Package ‘ggplot2’: elegant graphics for data analysis. New York: Springer; 2016.

    Google Scholar 

  64. Kassambara A. rstatix: Pipe-friendly framework for basic statistical tests. 2021. https://CRAN.R-project.org/package=rstatix. Accessed 15 Nov 2022.

  65. Lex A, Gehlenborg N, Strobelt H, Vuillemot R, Pfister H. UpSet: visualization of intersecting sets. IEEE Trans Vis Comput Graph. 2014;20:1983–92.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Kolde R. pheatmap: pretty heatmaps. 2019. https://cran.r-project.org/web/packages/pheatmap/index.html. Accessed 15 Nov 2022.

  67. Csardi G, Nepusz T. The igraph software package for complex network research. Inter J Complex Syst. 2006;1695:1–9.

    Google Scholar 

  68. Fruchterman TMJ, Reingold EM. Graph drawing by force-directed placement. Softw Pract Exp. 1991;21:1129–64.

    Article  Google Scholar 

  69. Wilkins D. gggenes: draw gene arrow maps in 'ggplot2'. 2020. https://CRAN.R-project.org/package=gggenes. Accessed 15 Nov 2022.

  70. Eren AM, Kiefl E, Shaiber A, Veseli I, Miller SE, Schechter MS, et al. Community-led, integrated, reproducible multi-omics with anvi’o. Nat Microbiol. 2021;6:3–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.

    Article  Google Scholar 

  72. Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Benedict MN, Henriksen JR, Metcalf WW, Whitaker RJ, Price ND. ITEP: an integrated toolkit for exploration of microbial pan-genomes. BMC Genom. 2014;15:8.

    Article  Google Scholar 

  74. van Dongen S, Abreu-Goodger C. Using MCL to extract clusters from networks. In: van Helden J, Toussaint A, Thieffry D, editors. Bacterial molecular networks: methods and protocols. New York: Springer; 2012. p. 281–95.

    Chapter  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge Dr. Fiona Crispie and Dr. Gaston Allendez for their help in sequencing our AmpC / ESBL E. coli isolates. Furthermore, we would like to thank Elena Anedda, Teagasc and Gillian Madigan and the staff of the Department of Agriculture, Food and the Marine for their help in sample preparation and analysis via MALDI-TOF. Additional thanks goes to the technical staff of the UCD Veterinary Hospital Bacteriology Diagnostic Laboratory for phenotypic resistance testing of five mcr-1.1-positive isolates via Vitek2®. We would like to thank the pig farmers for providing us access to their farms and for their co-operation.

Funding

This study was funded by the Teagasc Walsh Scholarship scheme (ref 2018027). Héctor Argüello is granted by the Spanish Ministry of Education under the programme “Beatriz Galindo” (BEAGAL-18-106).

Author information

Authors and Affiliations

Authors

Contributions

DE conducted the sampling, sample processing, analysis, data interpretation, and wrote the manuscript. All authors contributed in defining the research questions, designing the study, revision and editing of the manuscript. JMOS and EGM conducted together with DE the sampling on farms. HA, EGM, and CMB contributed to the statistical analysis and data interpretation. Supervision was provided by EGM, FL and CMB. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Catherine M. Burgess.

Ethics declarations

Ethics approval and consent to participate

This study was licensed by Teagasc Animal Ethics Committee.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1.

Summary of analysed AmpC / ESBL E. coli isolates obtained in this study. Sample abbreviations: FW0 = faecal sample at day of weaning; FW1 = faecal sample 1 week after weaning; FW2 = faecal sample 2 week after weaning; SWD = environmental swab sample of cleaned drinkers and feeders before interaction with piglets; SWW = environmental swab sample of cleaned walls and floors before interaction with piglets.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ekhlas, D., Sanjuán, J.M.O., Manzanilla, E.G. et al. Comparison of antimicrobial resistant Escherichia coli isolated from Irish commercial pig farms with and without zinc oxide and antimicrobial usage. Gut Pathog 15, 8 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s13099-023-00534-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13099-023-00534-3

Keywords