Genetic relatedness of multidrug resistant Escherichia coli isolated from humans, chickens and poultry environments
Antimicrobial Resistance & Infection Control volume 10, Article number: 58 (2021)
Inappropriate use of antimicrobial agents in animal production has led to the development of antimicrobial resistance (AMR) in foodborne pathogens. Transmission of AMR foodborne pathogens from reservoirs, particularly chickens to the human population does occur. Recently, we reported that occupational exposure was a risk factor for multidrug-resistant (MDR) Escherichia coli (E. coli) among poultry-workers. Here we determined the prevalence and genetic relatedness among MDR E. coli isolated from poultry-workers, chickens, and poultry environments in Abuja, Nigeria. This study was conducted to address the gaps identified by the Nigerian AMR situation analysis.
We conducted a cross-sectional study among poultry-workers, chickens, and poultry farm/live bird market (LBM) environments. The isolates were tested phenotypically for their antimicrobial susceptibility profiles, genotypically characterized using whole-genome sequencing (WGS) and in silico multilocus sequence types (MLST). We conducted a phylogenetic single nucleotide polymorphism (SNPs) analysis to determine relatedness and clonality among the isolates.
A total of 115 (26.8%) out of 429 samples were positive for E. coli. Of these, 110 isolates were viable for phenotypic and genotypic characterization. The selection comprised 47 (42.7%) isolates from poultry-workers, 36 (32.7%) from chickens, and 27 (24.5%) from poultry-farm or LBM environments. Overall, 101 (91.8%) of the isolates were MDR conferring resistance to at least three drug classes. High frequency of resistance was observed for tetracycline (n = 102; 92.7%), trimethoprim/sulfamethoxazole (n = 93; 84.5%), streptomycin (n = 87; 79.1%) and ampicillin (n = 88; 80%). Two plasmid-mediated colistin genes—mcr-1.1 harboured on IncX4 plasmids were detected in environmental isolates. The most prevalent sequence types (ST) were ST-155 (n = 8), ST-48 (n = 8) and ST-10 (n = 6). Two isolates of human and environmental sources with a SNPs difference of 6161 originating from the same farm shared a novel ST. The isolates had similar AMR genes and plasmid replicons.
MDR E.coli isolates were prevalent amongst poultry-workers, poultry, and the poultry farm/LBM environment. The emergence of MDR E. coli with novel ST in two isolates may be plasmid-mediated. Competent authorities should enforce AMR regulations to ensure prudent use of antimicrobials to limit the risk of transmission along the food chain.
Antimicrobial resistance (AMR), is one of the biggest threats to food safety and considered a One-Health issue with the potential of spreading to other countries since resistant pathogens do not recognize boundaries [1, 2]. Recently, we have shown the transmission of AMR E. coli among chickens, humans, and the poultry environment [3, 4]. Globally, antimicrobial agents are used in food animal production to ensure good health and productivity of the animals [5,6,7]. Multiple studies have shown that inappropriate use of these antimicrobial agents in food animal production particularly poultry has led to the development of AMR [8,9,10].
Commensal E. coli are known to be part of the normal flora of the gastrointestinal tracts of man and animals without causing any harm to their host [11, 12]. Several E. coli strains have been used as indicator organisms in various studies on AMR [11, 13]. Although commensal E. coli are harmless to the host, the bacteria can acquire resistance genes and act as a reservoir for the spread of multidrug resistance (MDR) to and from food to humans . The genetic structure of E. coli strains is usually influenced by several factors including the host and environment enabling the bacteria to acquire various AMR mechanisms [13, 14].
In September 2016, 193 member countries including Nigeria signed the United Nations General Assembly resolution to develop national action plans (NAP) on AMR . In November 2016, Nigeria established its AMR coordinating body at the Nigeria Center for Disease Control (NCDC), and in January 2017, a One-Health AMR Technical working group was inaugurated to conduct AMR situation analysis and develop Nigeria’s NAP . One of the data gaps identified from the AMR situation analysis was the paucity of AMR studies done in Nigeria across humans, food-producing animals, and the environment .
It has been documented that the continuous use of antimicrobial agents for therapeutic purposes against infections has led to the emergence of drug-resistant bacteria such as MDR E. coli . MDR bacteria have made it difficult to treat certain infections effectively with modern or conventional antimicrobial agents . AMR has resulted in treatment failure in human and animal populations, because of the emergence of MDR foodborne pathogens like E. coli arising from the abuse or misuse of antimicrobial agents . This scenario further deteriorates in Nigeria because of the increasing number of farmers who practice self-prescription as well as self-administration of antimicrobials to their animals [5, 20]. Poultry farmers have easy access to antimicrobials that are available over-the-counter without prescription  and evidence shows that farmers administer the antimicrobials repeatedly against non-responsive infections [20, 21]. These actions by the farmers further promote the emergence and spread of antimicrobial-resistant foodborne pathogens with serious implications on public health . Continuous administration of antimicrobial agents to chickens for prophylaxis, therapeutic, or growth promotion purposes increases the antibiotic selection pressure for resistance in the bacteria . Our recent publication demonstrates that occupational exposure over ten years to chickens on poultry farms or live bird markets (LBMs) was a risk factor for MDR E. coli among poultry workers in Abuja .
We hypothesized that chickens harbouring MDR E. coli as well as contaminated poultry farm or LBM environment can become potential sources for transmission of resistance genes to poultry workers exposed to chickens and the environment on poultry farms or markets. To better understand the association between MDR E. coli isolates recovered from humans, chickens and poultry environment, we investigated the genetic relatedness of MDR E. coli isolates from poultry-workers, chickens, and selected poultry farms/LBM environments in Abuja, Nigeria.
Our current study was part of a larger project conducted from December 2018 to February 2020 exploring MDR E. coli in humans, chickens, and the poultry farm/market environment. An aspect of this research exploring the risk factors for MDR E. coli among poultry workers has already been previously published .
Characterization of E. coli isolates
Of 429 samples collected in the course of the present study, 110 E. coli strains isolated from the stool of apparently healthy poultry workers, faecal samples obtained from chickens as well as from poultry litter and water obtained from farm and LBM environments were characterized. The sample collection procedures, isolation of E. coli from these samples, and antimicrobial susceptibility profiling of E. coli using the disk diffusion method have been described previously . Briefly, suspected dark pink E. coli colonies on MacConkey agar were streaked on Eosin Methylene Blue agar and incubated at 37 °C for 24 h. Isolates were confirmed as E. coli using Microbact GNB 24E (Oxoid, UK).
Genotypic Detection of E. coli isolates
Whole genome sequencing (WGS) of E. coli isolates
All isolates were subjected to WGS as previously described . Briefly, libraries for each E. coli isolates were prepared for WGS using a Nextera XT kit. We processed 0.3 ng/µL of DNA from each isolate using a Nextera XT DNA Sample Prep Kit (Illumina Inc., San Diego, CA), pooled together, and sequenced on an Illumina Miseq platform using a 2 × 250 paired-end approach (Illumina Inc., San Diego, CA). Raw sequencing reads were de-multiplexed and converted to fastq files using CLC Genomics workbench version 9.4 (Qiagen bioinformatics, Valencia, CA). The DNA sequences for each isolate were transferred to the National Center for Biotechnology Information (NCBI) database after which each isolate was assigned an accession number.
Bioinformatic analysis of WGS data
Antimicrobial susceptibility determinants of E. coli isolates
High-quality Illumina paired-end reads were assembled de novo into the draft genome sequence for each isolate using SPAdes assembler v.3.13.1 . In silico detection and typing of resistance genes was done using ResFinder 3.2, a Center for Genomics Epidemiology (CGE) bioinformatics tool (database version 2020–02-11), to determine the acquired AMR genes as well as assess chromosomal point mutations . For each isolate, we used between 95–100% identity to match individual genes to an annotated resistance gene. . In silico determination of the existing plasmid replicon types of each E. coli isolate was done by submitting the assembled genomes to PlasmidFinder 2.1, a CGE bioinformatics tool (database version 2020–04-02). The selected threshold for minimum percentage identity was 95% while the minimum coverage of the contig was set at 60% . The in silico plasmid MLST typing of replicons (IncHI2 and IncF) were determined by submitting the assembled genome to pMLST 2.0 (database version: 2020–04-20) bioinformatics tool on the CGE website .
Multilocus sequence typing (MLST) of MDR E. coli isolates
As previously described  in silico MLST-analyses of the E. coli isolates were determined using schemes demonstrated by Achtman which made use of allelic variation amongst seven housekeeping genes (adk, fumC, gyrB, icd, mdh, purA, and recA) to assign sequence types (STs) . We used whole-genome sequence data to generate the E. coli MLST assignment for each isolate that perfectly matched the alleles in the MLST database. MLST Finder 2.0, a CGE bioinformatics tool was used to assign STs to the isolates with 100% match against known MLST alleles while those without perfect matches were identified as unknown . Some isolates were assigned as a new type after matching with MLST alleles of unknown ST in the MLST database.
Determination of E. coli Phylogroups, SNPs calling and Phylogeny
Phylogroups of E. coli genomes were determined using an in silico Clermont typing method . The Clermont Typer web interface is hosted by CATIBioMed (IAME UMR 1137) and accessible at http://clermontyping.iame-research.center/.
Phylogenetic trees were constructed to determine the phylogenetic relatedness of the E. coli isolates using the technique known as SNP calling described by Kaas et al. . Briefly, the tool CSI Phylogeny, a CGE bioinformatics tool accessed online at https://cge.cbs.dtu.dk/services/CSIPhylogeny/ was used for SNP calling. The CSI phylogeny uses BWA to map raw reads to a reference sequence and uses Samtools for SNP calling. E. coli strain NCTC11129 was used as the reference strain for SNPs calling to identify variants present in the chromosome of each isolate. The selected thresholds used were: cut-offs for depth = 10x; SNP quality = 30; mapping quality = 25 and Z score = 1.96. The phylogenetic SNP-based maximum likelihood tree were annotated and visualized using the programs Figtree (http://tree.bio.ed.ac.uk/software/figtree/) and interactive Tree of Life tool—iTOL (http://itol.embl.de/itol.cgi). Pairwise SNP differences between genomes were computed to determine if isolates of different origins were related.
Antimicrobial susceptibility data were analyzed using Epi Info 7 software by computing frequencies and proportions. The 108 assembled E. coli genomes of the present study have been deposited by the Thakur Molecular Epidemiology Laboratory, NC State University (GenomeTrakr Project) in the NCBI database under the Bioproject ID number PRJNA293225. The remaining two isolates have accession obtained from the DNA Data Bank of Japan (DDBJ) as previously reported .
Antimicrobial susceptibility profile of E. coli isolates
A total of 110 E. coli strains were isolated from 122 human stool samples obtained from poultry workers on farms and LBMs; 111 faecal samples obtained from chickens on farms and LBMs; and 196 poultry litter and water samples obtained from farm and LBM environments. Of the 110 E. coli strains 42.7% (n = 47) were recovered from humans; 37.7% (n = 36) from chickens and 24.5% (n = 27) from poultry environment. High resistance rates were observed for tetracycline, trimethoprim/ sulfamethoxazole, streptomycin, ampicillin, nalidixic acid and gentamicin. On the contrary resistance to colistin, imipenem, ceftazidime, amoxicillin/clavulanic acid, cefuroxime, cefotaxime and ceftriaxone were quite low although colistin resistance rate of 11.8% in commensal E. coli is quite worrisome (Table 1).
Analysis of resistance profiles of the 110 isolates showed that a single isolate (0.9%) from a poultry farmer was susceptible to all antimicrobial drugs tested; 4 (3.6%) were resistant to only one antimicrobial drug, 4 (3.6%) were resistant to two antimicrobial drugs and interestingly 101 (91.8%) were MDR (resistant to three or more classes of antimicrobial drugs). The number of antimicrobials against which each isolate showed resistance was between one and thirteen. Surprisingly, a single isolate from a poultry farm was resistant to 13 out of 16 antimicrobials tested. The AMR phenotypes with AMP, CEP, CHL, CT, GEN, NAL, S, SXT, and TET profile had the highest frequency of 13.6% (n = 15). Figure 1 summarizes the multiple AMR patterns exhibited by the isolates.
Prevalence of MDR E. coli in humans, chickens and poultry farm/LBM environment
The overall prevalence of E. coli from all sources was 26.8% (n = 115), however, only 110 were further characterized due to viability as the remaining five isolates were mistakenly discarded. Of the 110 E. coli isolates, 91.8% (n = 101) were MDR E. coli. Of these MDR E. coli isolates 38.6% (n = 39), 34.7% (n = 35), and 26.7% (n = 27) were recovered from humans, chickens and poultry environment respectively (Fig. 2). Surprisingly, all the poultry environment isolates were MDR. Of the 101 MDR E. coli isolates 47.5% (n = 48) were MDR5 (resistant to more than 5 classes) and 38.6% (n = 39) were classified as XDR (resistant to 8 or more classes i.e. extensively drug-resistant isolates). Overall, 36.6% (n = 37) of the isolates originated from the LBMs while 63.4% (n = 64) originated from farms. Of the 39 XDR E. coli isolates 41% (n = 16), 33.3% (n = 13), and 25.6% (n = 10) were recovered from chickens, humans and the poultry environment respectively.
In silico AMR gene analysis of MDR E. coli isolates in humans, chickens and poultry environment
This study identified 57 different resistance determinants from 101 MDR E. coli isolates. Genes encoding resistance to aminoglycosides accounted for the majority with about 14 different determinants (aadA1, aadA2, aadA2b, aadA5, aadA16, armA, aac(3)-IIa, aac(3)-IId, aac(3)-Ib, aac(6)-Ib-cr, aph(3)-Ia, aph(3)-Ib, aph(6)-Id, ant(2)-Ia) detected. A high prevalence (70.3%) of aph(6)-Id, which is a plasmid-encoded gene, was also observed. About two-thirds of the isolates (67.3%) exhibited aph(3)-Ib gene, a metabolic enzyme that confers aminoglycoside resistance. The aac(3)-IId gene responsible for conferring gentamicin resistance was observed in 27.7% of the MDR E. coli isolates. We also detected aac(6)-Ib-cr gene, responsible for the reduction in ciprofloxacin activity in two MDR E. coli isolates. Six different variants of β-lactam resistance genes were detected (blaTEM-1, blaOXA-1, blaOXA-10, blaOXA-129, blaCTX-M-15, blaCTX-M-65) out of which blaCTX-M type was classical of the ESBL producing E. coli. Ten different fluoroquinolone resistance determinants were observed, an important antimicrobial on the WHO list, (qnrB1, qnrB19, qnrB52, qnrS1, qnrS2, qnrS3, qnrS7, qnrS11, qnrS13, aac(6)-Ib-cr) and associated with mutations in the gyrA, parC, and parE genes. We detected other important resistance determinants such as trimethoprim resistance (dfrA1, dfrA8, dfrA12, dfrA14, dfrA15, dfrA17, dfrA21, and dfrA27), macrolide resistance (mdfA, mphA, mefB, ermB, ereA, mphE and msrE), phenicol resistance (cmlA1, catA1, catA2, catB3, floR), rifampicin resistance (ARR-2 and ARR-3), sulphonamide resistance (sul1, sul2, sul3), tetracycline resistance (tetA, tetB, tetM) and plasmid-mediated colistin resistance gene (PMCR)—mcr-1.1.
Multi-locus sequence determination of MDR E. coli isolates
The 101 MDR E. coli isolates belonged to 66 different sequence types (ST), out of which one (1) was non-conclusive and eight (8) were new types. In the in silico MLST analysis of E. coli isolates, the following were observed to appear more than once: ST155 (7.9%; n = 8), ST48 (7.9%; n = 8), ST10 (5.9%; n = 6), ST1638 (4%; n = 4), ST398 (3%; n = 3), ST216 (3%; n = 3), ST226 (3%; n = 3), ST101 (2%; n = 2), ST117 (2%; n = 2), ST165 (2%; n = 2), ST206 (2%; n = 2), ST4663 (2%; n = 2), ST1286 (2%; n = 2), and ST1196 (2%; n = 2). The most prevalent STs are shown in Fig. 3.
Phylogroups of E. coli isolates from humans, chickens and poultry environment
A majority of the isolates belonged to phylogroup A (n = 61, 55.5%) followed by phylogroup B1 (n = 36, 32.7%) while the rest belonged to phylogroup G (n = 3, 2.7%); D (n = 2, 1.8%); E (n = 2, 1.8%); F (n = 2, 1.8%); B2 (n = 1, 0.9%); C (n = 1, 0.9%); clade I (n = 1, 0.9%) and clade IV (n = 1, 0.9%). Isolates with phylogroup A originated from workers (n = 36) and poultry environment (n = 13) while isolates recovered from chickens mostly belonged to phylogroup B1 (Fig. 4). Of the 36 E. coli isolates, belonging to phylogroup B1, 22.2% (n = 8); 50% (n = 18) and 27.8% (n = 10) were recovered from humans, chickens and the poultry environment respectively.
All isolates assigned ST10 (n = 6), ST218 (n = 3), ST398 (n = 3) and ST1638 (n = 4) belonged to phylogroup A. However, all but one isolate assigned ST48 (7/8) and ST226 (2/3) also belonged to phylogroup A while a majority with ST155 (7/8) and novel ST (5/8) belonged to phylogroup B1.
Plasmid replicon profiles of MDR E. coli isolates from humans, chickens and poultry environment
Forty (40) different plasmid replicon types were detected among 97 MDR E. coli isolates however, 4% (n = 4) did not harbour any plasmid replicons. The most prevalent plasmid replicons detected in descending order were p0111 (36.6%, n = 37); IncFIB(AP001918) (33.7%, n = 34); IncFII (18.8%, n = 19); ColpHAD28 (14.9%, n = 15); IncQ1 (13.9%, n = 14); IncFIB(K) (13.9%, n = 14); ColpVC (12.9%, n = 13); IncFIC(FII) (12.9%, n = 13); IncR (9.9%, n = 10); IncFII(pCoo) (9.9%, n = 10); IncY (9.9%, n = 10); IncX1 (8.9%, n = 9) and IncI1-I(gamma) (8.9%, n = 9). The plasmid replicons recovered from human isolates were more genetically diverse than those recovered from chickens and the poultry environment. Eighteen replicon types were common to isolates from all sources: p0111, IncFIB(AP001918), IncFII, ColpHAD28, IncQ1, IncFIB(K), ColpVC, IncFIC(FII), IncX1, IncFII(pCoo), IncI1-I (gamma), IncFII (29), IncFII(pHN7A8), IncFIA, Col156, IncHI2, IncHI2A and IncX4.
IncFIB(AP001918) was the most common among human isolates (n = 12) while p0111 was commonly detected in both chicken (n = 15) and poultry environment isolates (n = 14). Interestingly, IncFIB (pLF82), a phage plasmid was detected in one isolate recovered from the LBM environment. Eight different plasmids were observed to harbour AMR genes. The following AMR genes were carried on plasmid replicons: mcr-1.1 + IncX4 (n = 2); tetA + IncX1 (n = 1); qnrB19 + Col440I (n = 7); sul2 + IncQ1 (n = 5); aph(3)-Ib + IncQ1 (n = 1); blaTEM-1 + IncFIC(FII) (n = 1); mdf(A) + IncFIB (n = 1); qnrS13 + IncFII (n = 1) and aac(3)-IIa + IncHI1B (n = 1). The plasmid replicons harbouring the AMR genes was commonly detected in commensal E. coli isolates recovered from poultry workers, chickens and the poultry environment.
Determination of pMLST for IncHI2 and IncF plasmid replicons
In silico pMLST identification and typing of IncHI2 and IncF plasmid replicons, were based on the combination of the alleles identified for the genes. For the IncHI2 the assigned ST was ST4 for isolates (MA_251 and MA_252) originating from a poultry farmworker and poultry litter on the same farm. The pMLST analysis assigned the two IncF plasmids for isolates MA_251 and MA_252 with ST[F18:A-B1]. It is interesting to note that although the plasmid structures of the two isolates were so similar, there was no clonal relationship between them.
Phylogenetic analysis of E. coli isolates from humans, chickens and poultry environments
All isolates assigned a phylogenetic group and ST were used to construct phylogenetic trees to determine if the isolates were genetically related or very diverse. Three phylogenetic trees were constructed: one for all the isolates (Fig. 5), one focusing on isolates with novel STs (Fig. 6a) and one with isolates of different origins assigned the same ST (Fig. 6b).
Overall, 110 isolates used to construct a maximum likelihood phylogenetic tree showed that the isolates were genetically diverse. The isolates were grouped based on similarities among them. Whole-genome (wg) SNPs-based phylogenetic analysis showed that some isolates sharing the same ST and phylogroups were not clonally related. The strains were in completely different clades in the SNP tree, separated by strains belonging to other STs. Three isolates with ST-1638 recovered from human, chicken and poultry environment were clustered together on the same clade. Pairwise SNP differences between the genomes of the isolates showed that they were not clonally related (Fig. 6b).
Two isolates of human and environmental origin with SNPs difference of 6161 were not clonally related although the isolates shared a novel ST and belonged to phylogroup B1 (Fig. 6a). The two isolates originating from the same farm had similar AMR gene profile (qnrB19, qnrS1, mdfA, mefB, sul 2, sul 3, blaTEM1, tetA, tetM, floR); as well as plasmid replicons (p0111, IncFIC(FII), IncHI2A, IncHI2, Col(pHAD28).
To the best of our knowledge, this is the first study to investigate the prevalence of MDR E. coli in poultry workers, chickens, and the poultry farm/market environments in Nigeria.
The first objective of this study was to characterize E. coli from poultry workers, chickens, and poultry environments. The unhygienic LBM environment where these chickens are sold acts as a reservoir of antimicrobial-resistant bacteria and eventually poses a health risk to people working in such an environment. A similar study done in the Netherlands reported a lower prevalence of MDR E. coli in chickens (23%) and chicken farmers (22%) when compared with the present study where a prevalence of 34.7% and 38.6% was detected in chickens and poultry workers respectively . Access to antimicrobials is better regulated in the Netherlands as well as the implementation of antimicrobial stewardship programs when compared to Nigeria and could explain the differences observed in both studies. A related study conducted in Bangladesh among poultry and poultry environment reported a much higher prevalence of E. coli (82%) from chicken faecal samples when compared with the findings from this study, with a much lower prevalence of 32.2% . Two similar studies performed among chickens from poultry farms in northern Nigeria also reported a much higher E. coli prevalence of 67.7%  and 69.8%  from cloacal swabs obtained from chickens on the farm. A possible explanation for the difference between studies carried out in northern Nigeria and our study could be due to the sample types collected as our study isolated E. coli from freshly dropped chicken faecal samples as opposed to cloacal swabs. A study conducted in Pakistan, reported a slightly higher E. coli prevalence of 36% from the poultry farm environment when compared to 26.1% obtained from the poultry environment in the present study . Our study findings are consistent with the reports of a related study carried out in Egypt where E. coli prevalence of 26.8% was obtained from the poultry environment . The similarity observed between our study findings and that of the Egyptian study may be due to similarities in poultry farming practices.
Our study examined AMR in E. coli isolates from poultry farm workers and chicken sellers and compared them to resistance rates of E. coli isolates from chickens and poultry farm/market environment. The patterns of resistance were similar for human and chicken isolates. High resistance rates were observed in isolates recovered from humans, chickens, and poultry farm/market environments for tetracycline, trimethoprim/ sulfamethoxazole, ampicillin, and streptomycin. This is in agreement with the findings of a study conducted in southwest Nigeria, where high resistance rates of E. coli isolates to beta-lactams, tetracyclines, macrolides, and sulfonamides were reported . This finding is not surprising as these antimicrobials are easily accessible and commonly used in poultry production in Nigeria for therapeutic purposes especially in the absence of antimicrobial stewardship programs .
Our study revealed a very high proportion (91.8%) of MDR E. coli isolates from all the sources. Interestingly, 83% of human, 97% of chicken, and 100% of poultry environment isolates were MDR E. coli. A possible explanation for this very high level of resistance observed could be because of the lack of prudent use of antimicrobials and the required regulation to support it resulting in over-the-counter availability often without prescription as reported in many studies [16, 38,39,40]. The potential transmission of the drug-resistant strains between different hosts could also be responsible for this observation because E. coli is a known zoonotic bacteria .
The most common beta-lactamase gene observed in this study was the blaTEM-1 gene, which confers ampicillin resistance in E. coli isolates and is in agreement with a previous study that reported ampicillin-resistant E. coli isolates in food, humans, and healthy animals . Our study however, did not detect any genes encoding carbapenem-hydrolyzing enzymes in any of the isolates although phenotypic characterization showed that 10.9% of the isolates were carbapenem-resistant. This may possibly be as a result of borderline interpretation of breakpoint settings between resistance and susceptibility. The present study identified one of the most important AMR genes , being the PMCR gene—mcr 1.1 harboured on IncX4 plasmids in two isolates recovered from the poultry environment. Evidence shows that the IncX4 plasmids harbouring mcr-1 genes have been detected in human and animal E. coli isolates however, our study recovered these plasmids from the poultry environment . Another study conducted in China also detected PMCR genes—mcr 1 in E. coli isolates sourced from the aquatic environment  however, the mcr 3.1 gene was detected in a human Salmonella case in the US . This further buttresses that mcr gene has spread across multiple pathogens.
Our study highlights that poultry workers, chickens, and the poultry environments share identical plasmid replicons and this is consistent with the literature [46, 47]. The IncF plasmids reported as one of the epidemic plasmids were observed in humans, chickens, and the poultry environments to harbour different AMR genes; blaTEM-1, mdf(A) and qnrS13 in the present study and these are consistent with the literature . The IncQ1 plasmids were detected in isolates with ST48 recovered from chickens and poultry farm environments harbouring the sul2 genes that confer sulphonamide resistance and this is consistent with reports of other studies [43, 48]. The poor biosecurity measures, unhygienic practices in poultry farms and LBMs, and occupational exposure of over ten years are factors that predispose these humans to get infected with these drug-resistant bacterial strains .
To determine the genetic relatedness of the isolates, we analyzed by WGS, E. coli recovered from humans, chickens, and poultry environments. Our results revealed that these isolates showed very diverse genetic profiles. Common STs were assigned based on MLST including ST155, ST48, ST10, ST1638, and ST398 in isolates from humans, chickens, and poultry environments, although ST155 was mostly detected in isolates of poultry origin at the LBM. The most common ST detected among isolates recovered from the poultry farm environment was ST48. Previous studies have reported that E. coli with ST48 in phylogroup A has been detected in healthy volunteers, seafood, and water [49,50,51]. Our study detected ST48 in E. coli recovered from healthy people, chickens, and the poultry environment. E. coli strains with ST10 have previously been reported as being emerging and pathogenic as implicated in human infections although MDR strains with ST10 have also been detected in poultry and other animal sources . Our study detected MDR E. coli strains with ST10 in healthy individuals, poultry manure, and water. A possible explanation could be that this is becoming an emerging public health issue arising due to possible mutations in the bacteria.
The majority of E. coli isolates in this present study belonged to phylogroup A (55.5%) and phylogroup B1 (32.7%). Most human and poultry environment isolates belonged to phylogroup A while majority of the chicken isolates belonged to phylogroup B1. Our study findings are in agreement with the results of a similar study conducted in Pakistan that reported that phylogroups B1 and A were the most prevalent detected among human and animal E. coli isolates . Interestingly, a study carried out in south-west Nigeria reported that chicken E. coli isolates were evenly distributed into phylogroups A and B1 while phylogroup B1 was the most prevalent among human isolates . Previous studies also showed that E. coli isolates belonging to phylogroup B2 are usually the most virulent, hence MDR [54,55,56]. However, our study observed that majority of the isolates, which belonged to either phylogroups A, and B1 were MDR. This is consistent with findings from a similar study conducted in south-west Nigeria which reported that isolates belonging to phylogroups B1 and A were MDR . Our study findings are not surprising and consistent with the literature that most commensal E. coli belong to phylogenetic groups A and B1 [57, 58]. However, it is worrisome that these indicator bacteria have become MDR with a negative impact on public health since they could be transferred to more virulent strains or species thus causing disease.
The phylogenetic SNP tree rooted using NCTC11129 reference strain revealed that the isolates were genetically diverse among the identified STs. Two unrelated isolates of human and environmental origin belonging to phylogroup B1 and sharing a novel ST, had Col440I replicons harbouring the qnrB19 genes that confer quinolone resistance and consistent with the literature . In silico pMLST typing of the two isolates further confirmed that the isolates shared the same plasmids: IncHI2[ST-4] and IncF[ST-F18:A-:B1]. The two isolates although not clonally related, shared the same plasmids (Col440I) harbouring AMR genes (qnrB19) possibly due to horizontal gene transfer. Studies have shown that the IncF and IncHI2 plasmids mainly found in E. coli strains, are frequently detected in humans and animals serving as reservoirs for the spread of AMR genes and have been associated with MDR E. coli [43, 59]. This evidence supports our study results and explanation of a possibility of horizontal gene transfer of AMR genes harboured in the plasmids. Our study did not find evidence of the clonal spread of MDR E. coli at the human-animal-environment interface; however, our findings suggest that mobile genetic elements may have facilitated the horizontal transfer of MDR genes between the plasmids among commensal E. coli which could potentially mutate into real pathogens with serious public health implications .
MDR E.coli isolates were found to be prevalent amongst poultry-workers, chickens, and poultry farm/market environments in Abuja, Nigeria. The highest resistance rates among MDR E. coli isolates were observed to tetracycline, sulphonamides, penicillins, aminoglycosides, and quinolones which are classes of antimicrobials commonly used in poultry production for treating avian diseases in Abuja. ST-155, ST-48, and ST-1638 were the only STs detected in humans, chickens, and poultry farm/LBM environments in our study. Our findings showed the emergence and spread of MDR E. coli with novel-ST from a poultry farm environment to a poultry farmer, which may have resulted from horizontal transfer of AMR genes harboured in plasmids. Consequent upon these, healthcare and poultry-workers should be educated on the fact that people in proximity with poultry are a high-risk group for faecal carriage of MDR E. coli. Competent authorities should enforce AMR regulation to ensure prudent use of antimicrobials to limit the risk of transmission along the food chain and to poultry workers. Farmers should be discouraged from indiscriminate use of antimicrobials in poultry production and encouraged to adopt preventive measures by observing biosecurity as well as good management practices on their farms.
Availability of data and materials
The datasets used and analyzed during the present study are available upon reasonable request from the corresponding author. All data generated or analyzed during this study are also included in this published article [and Additional file 1]. Raw sequence data have been submitted to NCBI (https://www.ncbi.nlm.nih.gov/) under accession no PRJNA293225.
Center for Genomics Epidemiology
DNA Data Bank of Japan
- E. coli :
Interactive Tree of Life
Live bird market
Multilocus sequence typing
National Action Plan
National Center for Biotechnology Information
Nigeria Center for Disease Control
Plasmid multilocus sequence typing
Plasmid-mediated colistin resistance
Single nucleotide polymorphisms
Simple sequence repeat
Parmley J, Leung Z, Léger D, Finley R, Irwin R, Pintar K, et al. One Health and Food Safety—The Canadian Experience: A Holistic approach toward enteric bacteria pathogens and antimicrobial resistance surveillance. Natl Acad Press. 2012. Available online at https://www.ncbi.nlm.nih.gov/books/NBK114511/.
Manyi-Loh C, Mamphweli S, Meyer E, Okoh A. Antibiotic use in agriculture and its consequential resistance in environmental sources: potential public health implications. Molecules. 2018. https://doi.org/10.3390/molecules23040795.
Aworh MK, Kwaga J, Okolocha E, Mba N, Thakur S. Prevalence and risk factors for multi-drug resistant Escherichia coli among poultry workers in the Federal Capital Territory, Abuja, Nigeria. PLoS ONE. 2019;14:e0225379. https://doi.org/10.1371/journal.pone.0225379.
Aworh MK, Kwaga J, Okolocha E, Harden L, Hull D, Hendriksen RS, et al. Extended-spectrum ß-lactamase-producing Escherichia coli among humans, chickens and poultry environments in Abuja, Nigeria. One Health Outlook. 2020;2:8. https://doi.org/10.1186/s42522-020-00014-7.
Alhaji NB, Isola TO. Antimicrobial usage by pastoralists in food animals in North-central Nigeria: the associated socio-cultural drivers for antimicrobials misuse and public health implications. One Health. 2018;6:41–7.
Van TTH, Yidana Z, Smooker PM, Coloe PJ. Antibiotic use in food animals worldwide, with a focus on Africa: pluses and minuses. J Glob Antimicrob Resist. 2020;20:170–7.
Hao H, Cheng G, Iqbal Z, Ai X, Hussain HI, Huang L, et al. Benefits and risks of antimicrobial use in food-producing animals. Front Microbiol. 2014;5:288. https://doi.org/10.3389/fmicb.2014.00288.
Marshall BM, Levy SB. Food animals and antimicrobials: impacts on human health. Clin Microbiol Rev. 2011;24:718–33.
Zhuge X, Ji Y, Tang F, Sun Y, Jiang M, Hu W. Population structure and antimicrobial resistance traits of Avian-origin. Transbound Emerg Dis. 2017;66:1920–9.
Agyare C, Etsiapa Boamah V, Ngofi Zumbi C, Boateng Osei F. Antibiotic use in poultry production and its effects on bacterial resistance. In: Antimicrobial resistance—a global threat. IntechOpen; 2019. https://doi.org/10.5772/intechopen.79371.
Delmani F, Jaran AS, Tarazi YA, Masaadeh H, Zaki O, Irbid T. Characterization of ampicillin resistant gene (blaTEM-1) isolated from E. coli in Northern Jordan. Asian J Biomed Pharm Sci. 2017;7:11–5.
Sarowska J, Futoma-Koloch B, Jama-Kmiecik A, Frej-Madrzak M, Ksiazczyk M, Bugla-Ploskonska G, et al. Virulence factors, prevalence and potential transmission of extraintestinal pathogenic Escherichia coli isolated from different sources: recent reports. Gut Pathogens. 2019;11:10. https://doi.org/10.1186/s13099-019-0290-0.
Szmolka A, Nagy B. Multidrug resistant commensal Escherichia coli in animals and its impact for public health. Front Microbiol. 2013;4:258. https://doi.org/10.3389/fmicb.2013.00258.
Tenaillon O, Skurnik D, Picard B, Denamur E. The population genetics of commensal Escherichia coli. Nat Rev Microbiol. 2010;8:207–17. https://doi.org/10.1038/nrmicro2298.
WHO. At UN, global leaders commit to act on antimicrobial resistance. 2016. https://www.who.int/news-room/detail/21-09-2016-at-un-global-leaders-commit-to-act-on-antimicrobial-resistance. Accessed 8 May 2020.
Nigeria Centre for Disease Control. Antimicrobial use and resistance in Nigeria. 2017;:1–158. http://www.ncdc.gov.ng/themes/common/docs/protocols/56_1510840387.pdf.
Tanwar J, Das S, Fatima Z, Hameed S. Multidrug resistance: an emerging crisis. Interdiscip Perspect Infect Dis. 2014. https://doi.org/10.1155/2014/541340.
Frieri M, Kumar K, Boutin A. Antibiotic resistance. J Infect Public Health. 2017;10:369–78.
Manyi-Loh C, Mamphweli S, Meyer E, Okoh A. Antibiotic use in agriculture and its consequential resistance in environmental sources: Potential public health implications. Molecules. 2018;23:795.
Adesokan HK, Akanbi IO, Akanbi IM, Obaweda RA. Pattern of antimicrobial usage in livestock animals in south-western Nigeria: the need for alternative plans. Onderstepoort J Vet Res. 2015;82(1):816. https://doi.org/10.4102/ojvr.v82i1.816.
Moyer MW. The looming threat of factory superbugs. Sci Am. 2016;315:70–9.
Trung NV, Carrique-Mas JJ, Thi Hoa N, Mai HH, Tuyen HT, Campbell JI, et al. Prevalence and risk factors for carriage of antimicrobial-resistant Escherichia coli on household and small-scale chicken farms in the Mekong Delta of Vietnam. J Antimicrob Chemother. 2015;70:2144–52. https://doi.org/10.1093/jac/dkv053.
Miles TD, McLaughlin W, Brown PD. Antimicrobial resistance of Escherichia coli isolates from broiler chickens and humans. BMC Vet Res. 2006;2:7. https://doi.org/10.1186/1746-6148-2-7.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.
Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4. https://doi.org/10.1093/jac/dks261.
Carattoli A, Zankari E, Garciá-Fernández A, Larsen MV, Lund O, Villa L, et al. In Silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58:3895–903.
Wirth T, Falush D, Lan R, Colles F, Mensa P, Wieler LH, et al. Sex and virulence in Escherichia coli: an evolutionary perspective. Mol Microbiol. 2006;60:1136–51.
Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, et al. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50:1355–61.
Beghain J, Bridier-Nahmias A, Le NH, Denamur E, Clermont O. ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping. Microb Genom. 2018;4(7):e000192. https://doi.org/10.1099/mgen.0.000192.
Kaas RS, Leekitcharoenphon P, Aarestrup FM, Lund O. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS ONE. 2014;9:e104984. https://doi.org/10.1371/journal.pone.0104984.
van den Bogaard AE. Antibiotic resistance of faecal Escherichia coli in poultry, poultry farmers and poultry slaughterers. J Antimicrob Chemother. 2001;47:763–71.
Akond MA, Alam S, Hassan SMR, Shirin M. Antibiotic resistance of Escherichia coli Isolated from poultry and poultry environment of Bangladesh. Am J Environ Sci. 2009;5(1):47–52. https://doi.org/10.3844/ajessp.2009.47.52.
Kwoji ID, Musa JA, Daniel N, Mohzo DL, Bitrus AA, Ojo AA, et al. Extended-spectrum beta-lactamase-producing Escherichia coli in chickens from small-scale (backyard) poultry farms in Maiduguri, Nigeria. Int J One Health. 2019;5:26–30.
Geidam YA, Ambali AG, Onyeyili PA. Detection and antibiotic sensitivity pattern of avian pathogenic Escherichia coli strains among rural chickens in the arid region of north-eastern Nigeria. Vet World. 2012;5:325. https://doi.org/10.5455/vetworld.2012.325-329.
Ur Rahman S, Ahmad S, Khan I. Incidence of ESBL-producing-Escherichia coli in poultry farm environment and retail poultry meat. Pak Vet J. 2018;39:116–20. https://doi.org/10.29261/pakvetj/2018.091.
Enany ME, Algammal AM, Nasef SA, Abo-Eillil SAM, Bin-Jumah M, Taha AE, et al. The occurrence of the multidrug resistance (MDR) and the prevalence of virulence genes and QACs resistance genes in E. coli isolated from environmental and avian sources. AMB Express. 2019;9:192. https://doi.org/10.1186/s13568-019-0920-4.
Olowe OA, Adefioye OJ, Ajayeoba TA, Schiebel J, Weinreich J, Ali A, et al. Phylogenetic grouping and biofilm formation of multidrug resistant Escherichia coli isolates from humans, animals and food products in South-West Nigeria. Sci Afr. 2019;6:e00158. https://doi.org/10.1016/j.sciaf.2019.e00158.
Oluwasile B, Agbaje M, Ojo O, Dipeolu M. Antibiotic usage pattern in selected poultry farms in Ogun state. Sokoto J Vet Sci. 2014;12:45. https://doi.org/10.4314/sokjvs.v12i1.7.
Akinyandenu O, Akinyandenu A. Irrational use and non-prescription sale of antibiotics in Nigeria: a need for change. J Sci Innov Res. 2014;3(2):251–7.
Awosan KJ, Ibitoye PK, Abubakar AK. Knowledge, risk perception and practices related to antibiotic resistance among patent medicine vendors in Sokoto metropolis, Nigeria. Niger J Clin Pract. 2018;21(11):1476–83. https://doi.org/10.4103/njcp.njcp_69_18.
Briñas L, Zarazaga M, Sáenz Y, Ruiz-Larrea F, Torres C. β-lactamases in ampicillin-resistant Escherichia coli isolates from foods, humans, and healthy animals. Antimicrob Agents Chemother. 2002;46:3156–63.
WHO. WHO | WHO list of Critically Important Antimicrobials (CIA). 2019. http://www.who.int/foodborne_disease/resistance/cia/en/#.UiMEZ7zmSDA.mendeley.
Rozwandowicz M, Brouwer M, Fischer J, Wagenaar JA, Gonzalez-Zorn B, Guerra B, et al. Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae. J Antimicrob Chemother. 2018;73:1121–37. https://doi.org/10.1093/jac/dkx488.
Zhu L, Zhou Z, Liu Y, Lin Z, Shuai X, Xu L, et al. Comprehensive understanding of the plasmid-mediated colistin resistance gene mcr-1 in aquatic environments. Cite This Environ Sci Technol. 2020;54:1603–13. https://doi.org/10.1021/acs.est.9b05919.
Monte DF, Nelson V, Cerdeira L, Keelara S, Greene S, Griffin D, et al. Multidrug-and colistin-resistant Salmonella enterica 4,,12:i:-sequence type 34 carrying the mcr-3.1 gene on the IncHI2 plasmid recovered from a human. J Med Microbiol. 2019;68:986–90.
Nguyen VT, Jamrozy D, Matamoros S, Carrique-Mas JJ, Ho HM, Thai QH, et al. Limited contribution of non-intensive chicken farming to ESBL-producing Escherichia coli colonization in humans in Vietnam: an epidemiological and genomic analysis. J Antimicrob Chemother. 2019;74:561–70. https://doi.org/10.1093/jac/dky506.
Salinas L, Cárdenas P, Johnson TJ, Vasco K, Graham J, Trueba G. Diverse commensal escherichia coli clones and plasmids disseminate antimicrobial resistance genes in domestic animals and children in a semirural community in ecuador. mSphere. 2019;4(3):1–8. https://doi.org/10.1128/msphere.00316-19.
Clément M, Keller PM, Bernasconi OJ, Stirnimann G, Frey PM, Bloemberg GV, et al. First clinical case of in vivo acquisition of DHA-1 plasmid-mediated AmpC in a Salmonella enterica subsp. Enterica Isolate. Antimicrob Agents Chemother. 2019. https://doi.org/10.1128/AAC.00992-19.
Ben Sallem R, Ben Slama K, Estepa V, Jouini A, Gharsa H, Klibi N, et al. Prevalence and characterisation of extended-spectrum beta-lactamase (ESBL)-producing Escherichia coli isolates in healthy volunteers in Tunisia. Eur J Clin Microbiol Infect Dis. 2012;31:1511–6.
Said LB, Hamdaoui M, Jouini A, Boudabous A, Slama KB, Torres C, et al. First detection of CTX-M-1 in extended-spectrum b-lactamase-producing Escherichia coli in seafood from Tunisia. J Food Prot. 2017;80:1877–81.
Madec JY, Haenni M, Ponsin C, Kieffer N, Rion E, Gassilloud B. Sequence type 48 Escherichia coli carrying the blaCTX-M-1 IncI1/ST3 plasmid in drinking water in France. Antimicrob Agents Chemother. 2016;60:6430–2.
Manges AR, Johnson JR. Food-borne origins of escherichia coli causing extraintestinal infections. Clin Infect Dis. 2012;55:712–9.
Umair M, Mohsin M, Ali Q, Qamar MU, Raza S, Ali A, et al. Prevalence and genetic relatedness of extended spectrum-β-lactamase-producing Escherichia coli among humans, cattle, and poultry in Pakistan. Microb Drug Resist. 2019;25:1374–81.
Iranpour D, Hassanpour M, Ansari H, Tajbakhsh S, Khamisipour G, Najafi A. Phylogenetic groups of escherichia coli strains from patients with urinary tract infection in Iran based on the new Clermont phylotyping method. Biomed Res Int. 2015;2015:5–12.
Yang F, Zhang S, Shang X, Wang L, Li H, Wang X. Characteristics of quinolone-resistant Escherichia coli isolated from bovine mastitis in China. J Dairy Sci. 2018;101:6244–52.
Moreno E, Prats G, Sabaté M, Pérez T, Johnson JR, Andreu A. Quinolone, fluoroquinolone and trimethoprim/sulfamethoxazole resistance in relation to virulence determinants and phylogenetic background among uropathogenic Escherichia coli. J Antimicrob Chemother. 2006;57:204–11.
Clermont O, Bonacorsi S, Bingen E, Bonacorsi P. Rapid and Simple Determination of the Escherichia coli Phylogenetic Group Rapid and Simple Determination of the Escherichia coli Phylogenetic Group Downloaded from http://aem.asm.org/ on January 18 , 2013 by University of Pretoria : Academic Information S. Appl Environ Microbiol. 2000;66:4555–8.
Ramadan HH, Jackson CR, Taha SA, Moawad AA, Barrett JB, Woodley TA. Contribution of healthy chickens to antimicrobial-resistant Escherichia coli associated with human extraintestinal infections in Egypt. Vector-Borne Zoonotic Dis. 2018;18:408–16.
Yang Q-E, Sun J, Li L, Deng H, Liu B-T, Fang L-X, et al. IncF plasmid diversity in multi-drug resistant Escherichia coli strains from animals in China. Front Microbiol. 2015;6:964. https://doi.org/10.3389/fmicb.2015.00964.
The authors would like to appreciate the support of Dr. Chikwe Ihekweazu, Director General, Nigeria Center for Disease Control, Abuja for providing laboratory space towards this research at the National Reference Laboratory, NCDC. Special appreciation goes to Dr. Adebayo Adedeji, Mrs. Mba Nwando, Dr. Abiodun Egwuenu, Miss Eme Ekeng, Mr. Akinpelu Afolabi, Mr. Micheal Popoola and Mr. Chris Chukwu of the NCDC, Abuja, towards the success of this research. The authors appreciate the efforts of members of Thakur Molecular laboratory, North Carolina State University; Lyndy Harden, Erin Harrell, Nigatu Atlaw, Ayanna Glaize, Dawn Hall, Morgan Young, Luke Raymond, Bryson Staley and Steven Branz towards the success of the research. Special appreciation goes to Dr. Muhammad Shakir Balogun and Dr. Bola Biliaminu Lawal of Nigeria FELTP for their contributions to the success of this research. We also acknowledge the contributions of Judit Szarvas and Pimlapas Leekitcharoenphon at the National Food Institute, Denmark Technical University towards the success of this research.
The author(s) received funding from FDA GenomeTrakr program grant 1U18FD00678801 for whole-genome sequencing of the isolates.
Ethics approval and consent to participate
Ethics approval for this study was obtained from the FCT Health Research Ethics Committee in July 2018 (Approval Number: FHREC/2018/01/84/16–07-18). Permission was obtained from farm and LBM managers before the commencement of the study. We assured participants of confidentiality for the information obtained. An informed consent form clearly stating the study details and benefits was duly signed by participants who were willing to participate in addition to the farm and LBM managers before sample collection. Ethics approval for this study was also obtained from the Scientific and Ethical Committee of the Ahmadu Bello University Committee on Animal Use and Care (Approval Number: ABUCAUC/2020/35).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Aworh, M.K., Kwaga, J.K.P., Hendriksen, R.S. et al. Genetic relatedness of multidrug resistant Escherichia coli isolated from humans, chickens and poultry environments. Antimicrob Resist Infect Control 10, 58 (2021). https://doi.org/10.1186/s13756-021-00930-x
- Escherichia coli
- Antimicrobial resistance
- Multidrug resistance
- Genetic relatedness
- One health