Next Article in Journal
Antibiotic Resistance in Shiga Toxigenic Escherichia coli Isolates from Surface Waters and Sediments in a Mixed Use Urban Agricultural Landscape
Next Article in Special Issue
Assessment and Antibiotic Resistance Profiling in Vibrio Species Isolated from Wild Birds Captured in Danube Delta Biosphere Reserve, Romania
Previous Article in Journal
Characterization of a Novel Variant of the Quinolone-Resistance Gene qnrB (qnrB89) Carried by a Multi-Drug Resistant Citrobacter gillenii Strain Isolated from Farmed Salmon in Chile
Previous Article in Special Issue
Insects, Rodents, and Pets as Reservoirs, Vectors, and Sentinels of Antimicrobial Resistance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Escherichia coli and Related Enterobacteriaceae and Examination of Their Phenotypic Antimicrobial Resistance Patterns: A Pilot Study at A Wildlife–Livestock Interface in Lusaka, Zambia

by
Emmanuel Kabali
1,
Girja Shanker Pandey
2,
Musso Munyeme
2,
Penjaninge Kapila
2,
Andrew Nalishuwa Mukubesa
2,
Joseph Ndebe
2,
John Bwalya Muma
2,
Charles Mubita
2,
Walter Muleya
3,
Elizabeth Muligisa Muonga
2,4,
Shuya Mitoma
5,
Bernard Mudenda Hang’ombe
6,
Anuwat Wiratsudakul
7,
Mai Thi Ngan
5,8,
Eslam Elhanafy
9,
Hala El Daous
5,9,
Nguyen Thi Huyen
10,
Wataru Yamazaki
11,
Tamaki Okabayashi
12,13,
Maiku Abe
14,
Junzo Norimine
12,13 and
Satoshi Sekiguchi
12,13,*
add Show full author list remove Hide full author list
1
Quality Assurance Unit, Director General’s Office, Zambia Medicines Regulatory Authority, Lusaka 10101, Zambia
2
Department of Disease Control, School of Veterinary Medicine, University of Zambia, Lusaka 10101, Zambia
3
Department of Biomedical Sciences, School of Veterinary Medicine, University of Zambia, Lusaka 10101, Zambia
4
Department of Environmental Health, School of Health Sciences, Eden University, Lusaka 10101, Zambia
5
Graduate School of Medicine and Veterinary Medicine, University of Miyazaki, Miyazaki 889-1692, Japan
6
Department of Para-Clinical Studies, School of Veterinary Medicine, University of Zambia, Lusaka 10101, Zambia
7
Department of Clinical Sciences and Public Health, Faculty of Veterinary Science, Mahidol University, Nakhon Pathom 73170, Thailand
8
Faculty of Veterinary Medicine, Vietnam National University of Agriculture, Hanoi 100000, Vietnam
9
Faculty of Veterinary Medicine, Benha University, Moshtohor, Toukh, Qalyubia 13736, Egypt
10
National Institute of Veterinary Research, Hanoi 100000, Vietnam
11
Center for Southeast Asian Studies, Kyoto University, Kyoto 606-8501, Japan
12
Department of Veterinary Science, Faculty of Agriculture, University of Miyazaki, Miyazaki 889-2192, Japan
13
Centre for Animal Disease Control, University of Miyazaki, Miyazaki 889-2192, Japan
14
Education and Research Center for Mathematical and Data Science, Hokkaido University, Hokkaido 060-0812, Japan
*
Author to whom correspondence should be addressed.
Submission received: 15 January 2021 / Revised: 22 February 2021 / Accepted: 22 February 2021 / Published: 26 February 2021
(This article belongs to the Special Issue Antimicrobial Resistance and the Environment: One Health Approach)

Abstract

:
A cross-sectional study was used to identify and assess prevalence and phenotypic antimicrobial resistance (AMR) profiles of Escherichia coli and other enterobacteria isolated from healthy wildlife and livestock cohabiting at a 10,000 acres game ranch near Lusaka, Zambia. Purposive sampling was used to select wildlife and livestock based on similarities in behavior, grazing habits and close interactions with humans. Isolates (n = 66) from fecal samples collected between April and August 2018 (n = 84) were examined following modified protocols for bacteria isolation, biochemical identification, molecular detection, phylogenetic analysis, and antimicrobial susceptibility testing by disc diffusion method. Data were analyzed using R software, Genetyx ver.12 and Mega 6. Using Applied Profile Index 20E kit for biochemical identification, polymerase chain reaction assay and sequencing, sixty-six isolates were identified to species level, of which Escherichia coli (72.7%, 48/66), E. fergusonii (1.5%, 1/66), Shigella sonnei (22.7%, 14/66), Sh. flexinerri (1.5%, 1/66) and Enterobacteriaceae bacterium (1.5%, 1/66), and their relationships were illustrated in a phylogenetic tree. Phenotypic antimicrobial resistance or intermediate sensitivity expression to at least one antimicrobial agent was detected in 89.6% of the E. coli, and 73.3% of the Shigella isolates. The E. coli isolates exhibited the highest resistance rates to ampicillin (27%), ceftazidime (14.3%), cefotaxime (9.5%), and kanamycin (9.5%). Multidrug resistance (MDR) was detected in 18.8% of E. coli isolates while only 13.3% Shigella isolates showed MDR. The MDR was detected among isolates from impala and ostrich (wild animals in which no antimicrobial treatment was used), and in isolates from cattle, pigs, and goats (domesticated animals). This study indicates the possible transmission of drug-resistant microorganisms between animals cohabiting at the wildlife–livestock interface. It emphasizes the need for further investigation of the role of wildlife in the development and transmission of AMR, which is an issue of global concern.

1. Introduction

Antimicrobial resistance (AMR) occurs when bacteria, viruses, fungi and parasites change over time and no longer respond to antimicrobial medicines (including antibiotics, antivirals, antifungals and antiparasitics) making infections harder to treat and increasing the risk of disease spread, severe illness and death [1]. According to WHO, AMR is occurring everywhere in the world, compromising humans’ ability to treat infectious diseases, and undermining many other advances in health and medicine [2]. Globally, it is estimated that antimicrobial resistance (AMR) will be responsible for 10 million deaths per year by 2050 if steps are not taken immediately to combat and prevent the spread of AMR. In addition, the cumulative economic cost of AMR, estimated to be 100 trillion USD, is more than one and a half times the annual world gross domestic product as of 2015 [2,3].
Therefore, the emerging threat of AMR to public health continues to attract significant attention worldwide [2,3,4]. In response to this threat, a global initiative led by the World Health Organization (WHO), the Food and Agricultural Organization and the World Organisation for Animal Health (OIE), encompassing a global action plan, has been implemented to combat AMR [2]. Various factors contribute to the emergence of AMR. These include: limited knowledge of the risk of AMR, inappropriate use of antimicrobial agents in humans, animals and plants, including overprescribing and dispensing [3,5], poor drug quality and lack of drug quality control and monitoring/testing facilities [3,6], lack of adherence to good distribution practices and inadequate regulation of distribution channels/facilities [6,7], patients not completing the entire prescribed treatment course, and poor/inadequate healthcare facilities for provision of clinical services to humans, animals and plants [3].
Specifically, limited knowledge of the risk of AMR at various levels within global and local ecosystems will continue to hamper effective action for the control of AMR. For instance, recent studies on AMR have shown that commensal microorganisms belonging to families such as Enterobacteriaceae are becoming increasingly important in the transmission of resistance genes [1,5,8]. Indicator organisms such as Escherichia coli are receiving increased research attention based on observation that most commensal microorganisms display inherent resistance to specific antimicrobial agents [1,8]. The importance of these organisms to AMR risks is reinforced by recent findings showing that the use of broad-spectrum antibiotic treatments tends to select for resistant microorganisms within the human microbiome [9,10]. With the ability to transmit AMR genes horizontally via mobile genetic elements, commensal bacteria may play a significant role in the development and spread of AMR among other microorganisms [5,11]. This becomes a concern when transmission occurs from commensal to pathogenic species, and/or when commensal microorganisms acquire virulence determinants through mutations, making them pathogenic. The AMR transmission dynamics become even more significant when it implicates immunocompromised individuals, due to chronic and lifelong conditions such as HIV/AIDS, diabetes and cancer, requiring them to depend on long-term treatments [1].
Relatively little is known about the prevalence and development of AMR in bacteria from wild animals or the environment, although there is an increasing number of reports on multidrug-resistant (MDR) organisms in natural environments [12]. An understanding of the occurrence of AMR in environmental pathogens and commensals is important, especially when there are close interactions among wildlife, livestock and humans. This is even more relevant given our increasing knowledge of horizontal transmission of resistance genes via mobile elements among microorganisms [5,11], and the likelihood of preservation of microorganisms in wild animals and the environment [13]. The increasing interactions among humans, livestock and wildlife increase the public health significance of microorganisms that can be transmitted from wildlife and the environment to humans and livestock [13,14,15,16,17,18,19]. Therefore, the purpose of the current study was to investigate the prevalence and types of AMR amongst bacteria from apparently healthy wildlife and livestock cohabiting at a game ranch in Lusaka, Zambia.

2. Results

2.1. Microbiological, Biochemical and Molecular Analysis for Identification of Enterobacteriaceae Fecal Isolates

A total of 84 fecal samples were collected, from which 66 presumptive E. coli isolates were obtained. From API biochemical tests of the 66 presumptive E. coli isolates, 59.1% (n = 39) were identified as E. coli Group 1, one isolate was identified as Burkholderia cepacia and one isolate was identified as Kluyvera species. The remaining 25 isolates (37.9%) could not be conclusively identified using available API 20E testing kits (Table 1).
Using PCR and cycle sequencing, the majority of the isolates were identified to strain level. Among isolates from impala, E. coli strains detected included E. coli FR–2, E. coli LD93–1, E. coli ECPF–16, and E. coli WP2–S18–ESBL–07. Among isolates from cattle, E. coli strains included E. coli O157 CFSAN 076619, E. coli 824422, and E. coli LD93–1. In ostrich, E. coli RHB04–C05 and E. coli YJ4 strains were detected, while in goats, E. coli UF–153 and E. coli WP5–W18–ESBL 11, and in pigs, E. coli PB7 DCRUST SK strains were detected. Shigella sonnei SE6–1 strain was detected at least in one of the samples from each animal species, while Shigella spp. NCCP–460 was detected in a sample from pig. Shigella flexinerri SD5 strain was detected in a sample from cattle.

2.2. Phylogenetic Analysis

Phylogenetic analysis of the16S rRNA gene sequences showed that all the isolates under study separated into two clusters namely A and B (Figure 1).
Cluster B was divided into minor clusters B1 and B2 with B2 further dividing into B2a and B2b. Within cluster A, only two wildlife sequences were present, these being E. coli originating from a buffalo and an ostrich. On the other hand, cluster B consisted of sequences from both wildlife and domestic animals as well as E. coli and Shigella reference sequences from Brazil and India. Most sequences in cluster B were E. coli with a few being Shigella. Many sequences in cluster B1 were E. coli obtained from pigs (n = 4) (Figure 1). Within cluster B2a, most sequences originated from cattle and these were two Shigella and eight E. coli sequences, respectively, while in cluster B2b, most sequences originated from impala (n = 5) with two being Shigella and the remaining three being E. coli. Overall, in cluster B, clustering based on type of organism or origin of the sample was not observed, instead sequences from different animals and bacterial organisms formed similar clusters, indicating the close similarity of the 16S rRNA gene in E. coli and Shigella from both wildlife and domestic animals (Figure 1).

2.3. Antimicrobial Susceptibility Testing

Antimicrobial susceptibility testing (AST) results indicated that Shigella isolates showed resistance or intermediate sensitivity to seven of the thirteen antimicrobials screened in this study, while the E. coli isolates showed resistance or intermediate sensitivity to eleven of the thirteen antimicrobials. Intermediate resistance to at least one antimicrobial agent was shown by 89.6% (n = 43) of the E. coli isolates and 73.3% (n = 11) of the Shigella isolates. Isolates obtained from impala showed the greatest diversity in antibiotic resistance, exhibiting resistance to eight different antimicrobial agents in total. In comparison, isolates from buffalo showed the least diversity, exhibiting resistance to none of the antimicrobial agents tested (Table 2). Of the E. coli isolates, 18.8% (n = 9) were MDR, while 1.3% (n = 2) of the Shigella isolates were MDR. The higher number of MDR isolates were recovered from impala (n = 3), while others were isolated from cattle (n = 2), ostrich (n = 2), pig (n = 2) and goat (n = 2) samples. Shigella MDR isolates were each recorded in ostrich and pig (Table 3). The MDR was observed in various combinations, including penicillins, cephalosporins, and quinolones combination (most commonly detected); penicillins, cephalosporins, and aminoglycosides combination; as well as penicillins, cephalosporins, and carbapenems combination.

3. Discussion

In this study, several bacteria strains were detected in wildlife using both biochemical and molecular techniques. A few similar studies have previously been conducted to examine AMR in apparently healthy pastoralist cattle in an interface area in Zambia by Mubita et al. [20], and in a more recent study, Mubita et al. [21] also applied a combined approach similar to this study (molecular characterization and phenotypic detection of AMR) in examining AMR in Salmonella isolates at the wildlife–livestock interface in Zambia. This combined approach has, however, not been applied in examining AMR in E. coli isolates at the interface in Zambia.
In Impala, for instance, there has not been any reported cases in which E. coli FR–2, E. coli LD93–1, E. coli ECPF–16, and E. coli WP2–S18–ESBL–07 strains were detected. These E. coli strains may not necessarily be commensal, even though they were isolated from apparently healthy animals. Specifically, the isolation of extended beta-lactamase producing bacteria in animals in which antimicrobials may not have been used might be evidence relevant in elucidating the role of wildlife in the transmission and maintenance of antimicrobial resistance genes in the wild. The isolation of Shigella sonnei SE6–1 strain in all animal species sampled, domesticated and wild alike, is also a first report in Zambia at a wildlife–livestock interface. This is an important finding as Shigella species have been reported to be among the common pathogens found in animals and are responsible for causing diarrheal diseases in humans in some African countries [22].
The use of the combined approach made it possible to identify isolates to strain level. Difference in identity of some isolates at species level were noted, between phenotypic and molecular identification. It was interesting to note that on selective media, all the 66 isolates showed phenotypes characteristic of E. coli. In most resource-limited settings, these isolates would be presumed to be E. coli, and this may result in suboptimal handling of the otherwise more potent and pathogenic organisms. Using the API system, 37.9 % of the isolates could not be categorically identified. This was either because the organisms identified were not in the API data base, or the identity confidence was less than 90 %. The cytochrome oxidase test, and supplementary tests (reduction of nitrates to nitrites (NO2) and N2 gas (N2); motility (MOB); oxidation of glucose (OF-O); fermentation of glucose (OF-F)) were not conducted, as confirmatory identification was done using available PCR, and sequencing techniques. Two isolates identified by API system as Burkholderia cepacia and Kluyvera species, both were identified as strains of E. coli on genotypic testing. In the final instance, the genotypic identities were relied upon.
In addition, the clustering on phylogenetic analysis (Figure 4) implied that the examined sections of the 16S rRNA gene in both the E. coli and Shigella bacteria were very similar and illustrates how closely related E. coli and some Shigella species are, based on the highly conserved 16S rRNA gene. From a One Health point of view, the close relation between E. coli and certain Shigella species has been illustrated for some time. For instance, in 2000, Johnson [23] describes the role of Shigella flexinerri and Shigella sonnei in causing endemic shigellosis in developed countries, while Shigella dysenteriae was reported to be responsible for the same disease in developing countries. Our findings, however, show that Shigella sonnei and Shigella flexinerri isolated from apparently healthy animals have the potential to cause shigellosis in humans, especially considering that the common routes of transmission of these obligate pathogens are the fecal–oral route, with food, water, fomites, insects, and direct contact [23], all of which are possible at the wildlife–livestock–human interface. Further, the presence of bacteria that are so closely related, and exhibiting AMR may complicate accurate diagnosis and institution of effective treatment, especially in settings where both diagnostic and treatment options are limited. The detection of MDR in some of these isolates greatly increases the importance of their presence in apparently healthy animals, from the above mentioned One Health viewpoint.
In this study, evidence of AMR occurring at the wildlife–livestock interface was demonstrated. Of note were the observed higher levels of AMR amongst isolates from impala compared with domestic animals. These preliminary findings are uniquely important given that the wildlife–livestock–human interface is a fertile location for the exchange of disease-causing microorganisms, including pathogens of animal health and zoonotic importance [14,15,16,17,18,19,24,25]. Zoonotic diseases of global importance have previously been shown to be transmitted at this interface, including rabies, Ebola, anthrax, trypanosomiasis and tuberculosis [14,15,16,17,18,19,26,27,28,29]. This complex disease situation is compounded by the presence of AMR [1,8]. With the implication of transmission via commensal microorganisms such as E. coli, AMR genes in these species become a grave concern. In our view, the transmission of resistance at the wildlife–livestock interface complicates the fight against AMR in human pathogens, as it has done for the treatment of many important animal and zoonotic diseases, because wild animals are a stable reservoir of microbial pathogens. This is further reinforced by the fact that AMR is inherent in every treatable disease complex.
AMR amongst bacteria of wildlife origin has been detected in several countries but has not been widely reported in Zambia [21,30]. In our study, higher levels of resistance were observed amongst isolates from impala and ostrich. This is an interesting finding because antimicrobials are rarely used in wild animals. It was, therefore, hypothesized that resistance may have been acquired via either direct contact with livestock and humans or indirectly via the environment, as levels of resistance were higher in isolates from both wild (impala and ostrich) and domestic (cattle and goat) animals. Drinking water from dams and streams, which drain from communities around the game ranch, is another possible route for the acquisition of resistance. In this case, it would have been expected that the prevalence of AMR amongst isolates from buffalo samples would have been higher. In addition, our findings were contrary to our expectation that the prevalence and patterns of AMR in isolates from buffalos would be alike those observed in cattle as they have similar grazing patterns. It was thought that the solitary and aggressive nature of wild buffalos, which minimizes close contact with other wildlife, domestic animals and humas, may have been a factor in the lower levels of resistance detected in the current study. Despite this, the detection of resistance at the levels observed in wildlife-derived isolates in the current study is a cause for concern. Most microorganisms found in the environment and in wild animals are highly resilient and can survive harsh environmental conditions, allowing them to persist in the environment. Their preservation in the community can also be influenced by wildlife conservation policies and practices [13].
The prevalence of AMR among isolates from pigs was relatively low. However, the pigs are raised under an intensive production system, which lowers the risk of transmission of AMR microorganisms from other animals. In this case, resistance would likely be acquired through contact with human attendants and through inappropriate antimicrobial use. The low levels of resistance in isolates from goats were interesting and contrary to our expectations. At the study site, goats were used for tick control and had more frequent contact with other livestock, wildlife and humans, which would make them more likely to acquire AMR microorganisms.
Our findings raise concerns about some wildlife products. In most instances, strict food safety procedures are not followed in the preparation and trade of game meat [13]. For instance, poachers of wild animals for game meat and trophies often do not follow approved meat processing procedures [13]. Therefore, subsequent handlers and eventual consumers of the illicit products risk infection with the resistant microorganisms. This problem can be compounded in situations where food control systems are poor, and hygiene in food production and processing is compromised, with bacterial contamination occurring at high rates [31]. Trophy hunters and handlers may also have some risk of contracting these resistant bacteria through direct contact. Observance and enforcement of basic biosafety, biosecurity, and food hygiene and safety measures in handling products of wild animal origin could address most of these risks.
Our study showed that the highest levels of resistance were expressed against ampicillin in most E. coli isolates. Resistance to third generation cephalosporins was higher amongst E. coli isolates obtained from impala, cattle and pig samples, with the greatest number of resistant isolates originating from impala. This finding was deemed important, as the tested third generation cephalosporins were all relatively new members with enhanced activity against Gram-negative bacteria, including β-lactamase-producing strains [32,33]. This is concerning given the higher proportion of E. coli isolates, some of which may be pathogenic and could also contribute to the transmission of resistance to other bacteria. There were a number of isolates which showed resistance to penicillin and also to some cephalosporins, hence considered as extended-spectrum β-lactamase-producing, which is particularly important given that β-lactam antimicrobials are classified by the WHO as being amongst the highest priority of the critically important antimicrobials for human medicine, and by the OIE as critically important veterinary antimicrobial agents [34,35]. In most recent studies conducted in Zambia, resistance to these antimicrobials was common, with similar trends observed worldwide [29,36].
The prevalence of MDR was slightly higher among E. coli isolates (18.8 %) compared with the Shigella isolates (13.4%). The high levels of MDR among these isolates justify the need to further investigate the mode(s) of transmission of the corresponding resistance gene(s) and identify the elements responsible for the resistance at the molecular level. However, in the absence of requisite capacities to conduct molecular characterization, the results obtained in this study could be useful for estimating the prevalence of AMR in interface areas. Specifically, our results showed that although some isolates were resistant to up to six individual antimicrobials, none of the isolates were resistant to more than four classes of antimicrobials, an important indicator when considering availability of appropriate antimicrobial treatment options for use in both humans and animals. It is important, however, to note that the observed MDR profiles could rapidly change considering the variety of drug resistance phenotypes observed in this study and the ease of horizontal transfer of resistance determinants among bacteria [5,8,11].

4. Materials and Methods

4.1. Ethical Consideration

This study was designed as part of a larger study, for which ethical approval was sought and obtained from the ERES Converge Institutional Review Board in Lusaka, Zambia (Ref: no. 2018–Apr–006). Consent to conduct the study and collect samples from selected wild animals was also obtained from the game ranch management. The sex of the animals was not considered as a specific factor in the current study. The study was conducted according to relevant guidelines for scientific research using animal subjects and was strictly observational.

4.2. Study Design

A cross-sectional study was used to examine the AMR profiles of bacteria from selected wild animals and livestock at a game ranch located about 40 km north of Lusaka, Zambia’s capital city (Figure 2A,B). The game ranch is completely fenced off, covering an area of about 10,000 acres, and is a conservancy with multiple species of wildlife and livestock, with more than 72 wild animal species reported. In the current study, focus was mainly on antelope and bovine species for sampling, with candidate species including blue wildebeest (Connochaetes taurinus), buffalo (Syncerus caffer), bushbuck (Tragelaphus sylvaticus), eland (Taurotragus oryx), hartebeest (Alcelaphus buselaphus), impala (Aepyceros melampus), Kafue lechwe (Kobus leche kafuensis), puku (Kobus vardonii), roan antelope (Hippotragus equinus), sable antelope (Hippotragus niger), and the common Tsessebe (Damaliscus lunatus).
Probability assessment, based on population size and frequency of citing at selected locations, was used to determine which wildlife species were most likely to come into close contact with livestock and humans. Grazing patterns were also considered in sample selection, and seasonality was a factor for convenience of sample collection. Impala and buffalo, with estimated park populations of 1000 and 70, respectively, were selected for sampling because they were determined to have been more interactive, and/or showing similar grazing patterns to the livestock kept on the ranch (goats and cattle). Ostriches (Struthio camelus) were also selected as they had close contact with the four other selected species, while pigs were selected because they had close contact with humans but were reared in enclosed pig houses.

4.3. Sample Collection

Fresh fecal samples were collected from the target species using the identification method described by Chame [37], between April and August 2018. The research team followed a buffalo herd and identified feces that had been dropped within the previous 4 h. Five herds of impala were also followed, with fecal samples collected upon observation of defecation when the animals stopped. Samples were also collected from a flock of ostriches, and samples were routinely collected from livestock (cattle, goats and pigs). All samples were collected within a 10 km radius of the main office at the game ranch, located at latitude −15.155472 and longitude 28.498028 (Figure 2C). Samples were collected aseptically and placed in sterile sampling bottles before being assigned identification numbers and labelled. Samples were placed in ice cooled boxes and transported for processing at the Microbiology and Public Health laboratories, Department of Disease Control, School of Veterinary Medicine, University of Zambia. For all samples, the process of bacterial isolation and identification was initiated on the day of sample collection (Supplement A). Fecal samples were then stored at −20 °C for further use or reference.

4.4. Bacteria Isolation

Approximately 1 g of each fecal sample was suspended in 9 mL of buffered peptone water (Oxoid, Cheshire, England) and incubated at 37 °C overnight (18–24 h) for pre-enrichment. Selection for E. coli was conducted by inoculating 10 µL of the enriched suspension onto lactose MacConkey agar (Oxoid) plates using sterile wire loops and incubated at 37 °C overnight. Four presumptive E. coli isolates from lactose MacConkey agar were stained using the Gram staining method and further inoculated onto eosin-methylene blue (EMB) agar (Oxoid) and incubating aerobically at 37 °C for 24 h. Up to two presumptive E. coli colonies (demonstrating the characteristic metallic green sheen on EMB agar) from each sample were then inoculated onto non-selective blood agar (Oxoid). Atypical E. coli colony from each sample was also inoculated onto blood agar and incubated at 37 °C overnight. Isolates were subcultured on nutrient agar to ensure purity.

4.5. Biochemical Identification

Isolates were subjected to biochemical tests by using the API® 20 E Gram-negative Microbial Identification Kit (bioMérieux, Midrand, South Africa) for identification of Enterobacteriaceae to species level. E. coli strain ATCC 25922 was used as positive control. Isolates that were not conclusively identified to genus or species level were collectively referred to as “unidentified isolates”.

4.6. Molecular Identification

4.6.1. DNA Extraction and Polymerase Chain Reaction

Chromosomal DNA was extracted from purified colonies of every isolate, after it was cultured on nutrient agar. DNA extraction was done using the QIAamp DNA Mini Kit (Qiagen, South Africa). Identification was done using a combination of three primers, 16E1 (forward primer), 16E2 and 16E3 (both reverse primers), targeted at the 16S gene and designed to detect all pathogenic and non-pathogenic E. coli strains, Shigella sonnei and Shigella flexineri, applying a modified PCR procedure as described by Tsen et al. [38]. Each 20 µL PCR mixture contained 1 x PCR buffer (2.0 µL), 260 µM/L dNTP (1.6 µL), 1 µM/l of each primer, 0.5 U of Taq DNA polymerase enzyme (0.1 µL), milli q water (11.9 µL), and template DNA (from sample) (2.0 µL). For each PCR cycle, denaturation, annealing, and extension were carried out at 95.0 °C for 2 min, 60.0 °C for 30 s, and 72 °C for 2 min, respectively. Final extension was done at 72 °C for 10 min. A total of 35 PCR cycles were performed using the Veriti™ 96-Well Thermal Cycler (Applied Biosystems™, Thermo Fisher Scientific, Waltham, MA, USA). PCR products were examined for a fragment of about 584 bp visualized in ethidium bromide-stained 1.5% agarose gel under UV light. Identification of the product band was done by molecular weight maker of the Takara 100 bp DNA Ladder (Takara, Japan).
For isolates that could not be detected using the 16E1/16E2/16E3 primers, the more general primers (P3mod and P5) for identification of bacteria by targeting the 16S gene were used to identify the isolates, using a modified PCR procedure described by Tsen et al. [38]. Each 20 µL PCR mixture contained 1 x PCR buffer (2.0 µL), 260 µM/L dNTP (1.6 µL), 1 µM/L of each primer, 0.5 U of Taq DNA polymerase enzyme (0.1 µL), milli q water (12.3 µL), and template DNA (from sample) (2.0 µL). Each PCR cycle, denaturation, annealing, and extension were carried out at 95.0 °C for 1 min 30 s, 57.0 °C for 1 min 30 s, and 72 °C for 3 min, respectively. Final extension was done at 72 °C for 10 min. A total of 35 PCR cycles were performed using the Veriti™ 96-Well Thermal Cycler from Applied Biosystems™ (Thermo Fisher Scientific, Waltham, MA, USA). A 700 bp fragment was visualized in ethidium bromide-stained 1.5% agarose gel.

4.6.2. Cycle Sequencing

The PCR product obtained was purified using the DNA purification kit (Promega, USA). The purified PCR products were sequenced using Big dye terminator version 3.1 (Thermo Fischer) according to the manufacturer’s specifications. The sequence PCR products obtained were further purified using the ethanol purification method followed by the denaturation of the purified products and capillary electrophoresis using the ABI 3500 Genetic Analyzer (Applied Biosystems).

4.6.3. Sequence Analysis

For sequence analysis, raw sequences were initially subjected to blast analysis on the NCBI website (http://blast.ncbi.nlm.nih.gov/Blast.cgi access on 16 January 2021) followed by sequence assembly using the ATGC plug-in incorporated in Genetyx ver. 12. (Genetyx Co., Tokyo, Japan). Previous published representative sequences of E. coli, Shigella sonnei, Shigella flexinerri, and Brucella were downloaded from the GenBank and together with the obtained sequences under study, were used to generate a multiple sequence alignment file using Clustal W1.6. The multiple sequence alignment file was further converted to a MEGA file format and utilized to generate a neighbor joining tree with 1000 bootstrap replicates as a confidence interval using MEGA ver. 6 [39]. All the sequences generated in this study have been deposited in the DNA Data Base of Japan with serial numbers LC599922 to LC599973.

4.7. Antimicrobial Susceptibility Testing (AST)

Antimicrobial susceptibility testing was conducted following standard operating procedure developed through the WHO Advisory Group on Integrated Surveillance of Antimicrobial Resistance (AGISAR) project in Zambia. The Clinical and Laboratory Standards Institute (CLSI) 2018 M100–S28 guidelines [40] were utilized for both testing method and interpretation of the AST. For each of the identified isolates and control strain E. coli 25922, 1–2 colonies were picked from the blood agar plates, inoculated onto nutrient agar (Oxoid) and incubated at 37 °C for 18 h. Following incubation, 1–2 colonies were picked from each of the nutrient agar plates and transferred into vials containing 4 mL of sterile normal saline (bioMérieux, Midrand, South Africa) to a McFarland turbidity of ~0.5. Each bacterial suspension was then inoculated onto two 100 mm Muller Hinton (Oxoid) agar plates (poured to a depth of 4 mm) to form a lawn.
Antimicrobial susceptibility test discs (BD BBL Sensi–Disc; BD Biosciences, Franklin Lakes, NJ, USA) were applied to the surfaces of the Muller Hinton plates in two groupings, as follows: Panel 1—ampicillin, 30 µg; ceftazidime, 30 µg; cefotaxime, 30 µg; kanamycin, 30 µg; imipenem, 10 µg; ceftriaxone, 30 µg; amoxicillin/clavulanic acid, 30 µg; Panel 2—nalidixic acid, 30 µg; amikacin, 30 µg; ciprofloxacin, 5 µg; trimethoprim/sulfamethoxazole, 25 µg; gentamicin, 10 µg; chloramphenicol, 30 µg. The antimicrobials were selected based on their availability on the Zambian market and/or their importance to public health globally.
Plates were then incubated at 37 °C for 18 h. Assessment of antimicrobial susceptibility was carried out manually using digital Vernier calipers by measuring the diameter (mm) of the zone of inhibition against a dark surface. AST results for control strain E. coli ATCC 25922 were compared against expected results to ensure the quality of the test. AST results were entered into an excel spreadsheet for interpretation using the Clinical and Laboratory Standards Institute 2018 guidelines [40].

4.8. Statistical Analysis

Statistical analysis of generated data was conducted for 2 parts, first for cycle sequence results, and second for antimicrobial susceptibility tests. Descriptive statistics were computed for all variables, while univariable analysis by cross-tabulation using Pearson’s chi-square test was carried out to test for variability among the isolates across animal species [41]. Further analysis of the variability in AMR profiles among isolates across animal species was also conducted. Isolates were considered MDR if they showed resistance to three or more classes of antimicrobials [42]. All data analyses were conducted using R software (version 3.6.0; R Foundation for Statistical Computing, Vienna, Austria).

4.9. Study Limitations

It is acknowledged that this pilot study had several limitations. First, the number of samples investigated was small; larger studies encompassing multiple wildlife–livestock–human interface areas and larger numbers of samples and isolates are recommended to further examine and confirm our findings. Second, while interesting AMR profiles were observed amongst the isolates, further molecular investigation should be conducted to characterize the presence and types of AMR genes responsible for the resistance expressed phenotypically. This information is essential for confirming the mode of resistance and similarities among resistance genes in isolates from wildlife and livestock. Third, the sample size was limited to a single wildlife–livestock interface, and therefore cannot generalize the findings to all interface areas in Zambia. However, these findings are useful for informing further surveillance activities in these areas that are not well studied in terms of AMR.

5. Conclusions

The role of the environment and wildlife in the maintenance and transmission of AMR is of growing concern. By conducting a cross-sectional survey of AMR amongst E. coli and closely related enterobacteria isolated from selected wild animals and livestock at a game ranch in Zambia, it was possible to identify the presence and estimate the prevalence of AMR within the sampled population. Our findings indicate that there is indeed a likelihood of transmission of AMR between livestock and wildlife cohabiting closely, although the transmission dynamics could not be explained based on these results. Our study also confirmed the urgent need to further investigate the role of wildlife and the shared environments in the transmission and maintenance of AMR, which is a One Health issue of global concern. This information is important for improving wildlife and livestock management, and for making decisions about the control of AMR in interface areas. Further investigation would also help ensure that the possibility of transmission of AMR genes from wildlife-associated bacteria to human microbiome through livestock or directly is minimized. Whole-genome sequencing could be useful in such investigation, and for resource-limited settings, the availability of online resources for bacterial typing offers rapid classification and source tracing, which is increasingly important in a globalized community. Such online resources offer rapid typing and phylogenetic relatedness linked to antibiotic resistance genes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2079-6382/10/3/238/s1, Supplement A: list of detected strains of bacteria isolates.

Author Contributions

E.K. collected samples, performed culture and sensitivity assays, analyzed the data and drafted the manuscript. W.M. drafted some sections of the manuscript. P.K. and E.M.M. collected samples and performed culture and sensitivity assays. S.M., C.M., A.W., M.T.N., H.E.D., N.T.H. and E.E. helped in study design and laboratory analysis. P.K., A.N.M., and E.K. conducted molecular analysis. E.K., J.N. (Joseph Ndebe) and W.M. carried out the bioinformatics analysis and W.M. generated the phylogenetic tree. C.M. and B.M.H. contributed to protocol validation, laboratory analysis and resource mobilization (reagents). M.M., W.Y., T.O. and J.N. (Junzo Norimine) contributed to the study design and the validation of the results. G.S.P. and S.S., J.B.M.; A.W.; M.A. conceptualized and supervised the study. S.S. also participated in sample collection and analyzed and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was part-funded by a grant from the Japanese Society for the Promotion of Science (JSPS) through the Ronpaku scholarship grants program (ID number R11712, JPY2017–JPY2019, both inclusive).

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and, ap-proved by ERES Converge Institutional Review Board in Lusaka, Zambia (Ref: no. 2018–Apr–006).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. Sequence data will be available on the DDJB website as soon as the applicable embargo period lapses.

Acknowledgments

We would like to acknowledge the support received from the Dean of the School of Veterinary Medicine, University of Zambia, and for availing the laboratory facilities for the study. We also thank the AMR technical teams at the Food and Agricultural Organization of the United Nations, OIE and WHO–Afro for use of protocols for AMR surveillance developed through their tripartite work. We also thank the Zambia AMR Coordinating Committee (AMRCC) for being supportive of this study. We are also grateful to the TAFS forum, for providing technical support and advice during the study design phase. We also thank Tamsin Sheen, PhD, from Edanz Group (www.edanzediting.com/ac) for editing a draft of this manuscript. Thanks also to Markus Lipp for final checks of manuscript technical soundness and compliance to ethical issues, prior to submission.

Conflicts of Interest

This work is a product of the authors research work and does not reflect the opinions of any employer or the organizations the authors are affiliated to. The authors declare no conflict of interest. The funders had no role in study design, in the collection, analysis or interpretation of data, in the writing of the report, or in the decision to submit the article for publication.

References

  1. WHO. Antimicrobial Resistance: Fact. Sheet. 2018. Available online: https://www.who.int/news-room/fact-sheets/detail/antimicrobial-resistance (accessed on 26 May 2019).
  2. WHO. World Health Assembly Addresses Antimicrobial Resistance, Immunization Gaps and Malnutrition. 2015. Available online: http://www.who.int/mediacentre/news/releases/2015/wha-25-may-2015/en/ (accessed on 12 February 2017).
  3. O’Neill, J. Review on Antimicrobial Resistance: Tackling Drug-Resistant Infections Globally. Tackling a Global Health Crisis: Initial Step. 2015. Available online: https://amr-review.org/sites/default/files/Report-52.15.pdf (accessed on 1 December 2018).
  4. Yang, Y.; Li, B.; Ju, F.; Zhang, T. Exploring variation of antibiotic resistance genes in activated sludge over a four-year period through a metagenomic approach. Environ. Sci. Technol. 2013, 47, 10197–10205. [Google Scholar] [CrossRef]
  5. Van Bijnen, E.M.E.; Paget, J.; de Lange-de Klerk, E.S.M.; den Heijer, C.D.J.; Versporten, A.; Stobberingh, E.E.; Goossens, H.; Schellevis, F.G. Antibiotic exposure and other risk factors for antimicrobial resistance in nasal commensal Staphylococcus aureus: An ecological study in 8 European countries. PLoS ONE 2015, 10, e0135094. [Google Scholar] [CrossRef] [Green Version]
  6. Pisano, E. Antimicrobial Resistance: What Does Medicine Quality Have to Do with it? For the Review on Antimicrobial Resistance. 2015. Available online: https://amr-review.org/sites/default/files/ElizabethPisaniMedicinesQualitypaper.pdf (accessed on 1 December 2018).
  7. Kabali, E.; Ward, K.; Bapoo, R. Review of the Effectiveness of Medicines Regulatory Systems in Zambia; Lambert Academic Publishing: Beau Bassin, Mauritius, 2018. [Google Scholar]
  8. FAO. Drivers, Dynamics and Epidemiology of Antimicrobial Resistance in Animal Production. 2016. Available online: http://www.fao.org/3/a-i6209e.pdf (accessed on 27 May 2019).
  9. Kho, Z.Y.; Lal, S.K. The human gut microbiome—A potential controller of wellness and disease. Front. Microbiol. 2018, 9, 1835. [Google Scholar] [CrossRef] [Green Version]
  10. Modi, S.R.; Collins, J.J.; Relman, D.A. Antibiotics and the gut microbiota. J. Clin. Investig. 2014, 124, 4212–4218. [Google Scholar] [CrossRef] [Green Version]
  11. Jaja, I.F.; Bhembe, N.L.; Green, E.; Oguttu, J.; Muchenje, V. Molecular characterization of antibiotic-resistant Salmonella enterica isolates recovered from meat in South Africa. Acta Tropica 2019, 190, 129–136. [Google Scholar] [CrossRef]
  12. Nikaido, H. Multidrug resistance in bacteria. Annu. Rev. Biochem. 2009, 78, 119–146. [Google Scholar] [CrossRef] [Green Version]
  13. Kabali, E.; Pandey, G.S.; Maboshe, M.; Munyeme, M. The wildlife-livestock-human interface, legislation and its impact on communities around Mosi-Oa-Tunya National Park in Zambia. J. Sustain. Dev. Africa 2015, 17, 125–145. [Google Scholar]
  14. Woods, R.; Reiss, A.; Cox-Witton, K.; Grillo, T.; Peters, A. The importance of wildlife disease monitoring as part of global surveillance for zoonotic diseases: The role of Australia. Trop Med. Inf. Dis. 2019, 4, 29. [Google Scholar] [CrossRef] [Green Version]
  15. UNEP Frontiers 2016 Report. Zoonoses: Blurred Lines of Emergent Disease and Ecosystem Health. 2016. Available online: http://uneplive.org/media/docs/early_warning/zoonoses.pdf (accessed on 26 May 2019).
  16. Wang, L.F.; Crameri, G. Emerging zoonotic viral diseases. Rev. Sci Tech. 2014, 33, 569–581. Available online: https://www.oie.int/doc/ged/D14089.PDF (accessed on 26 May 2019). [CrossRef]
  17. Kooriyama, T.; Okamoto, M.; Yoshida, T.; Nishida, T.; Tsubota, T.; Saito, A.; Tomonaga, M.; Matsuzawa, T.; Akari, T.; Nishimura, H.; et al. Epidemiological study of zoonoses derived from humans in captive chimpanzees. Primates 2013, 54, 89–98. [Google Scholar] [CrossRef]
  18. Alexander, K.A.; Lewis, B.L.; Marathe, M.; Eubank, S.; Blackburn, J.K. Modeling of wildlife-associated zoonoses: Applications and caveats. Vector Borne Zoonotic Dis. 2012, 12, 1005–1018. [Google Scholar] [CrossRef] [PubMed]
  19. OIE. Training Manual on Wildlife Diseases and Surveillance. 2010. Available online: https://www.oie.int/doc/ged/D12066.PDF (accessed on 27 May 2019).
  20. Mubita, C.; Syakalima, M.; Chisenga, C.; Munyeme, M.; Muma, B.; Chifumpa, G.; Hangombe, B.; Sinkala, P.; Simuunza, M.; Fukushi, H.; et al. Antimicrobial resistance of faecal Escherichia coli and Enterococci spp. from Apparently Health Pastoralist Cattle in the Interface Areas of the Kafue basin of Zambia-short communication. Veterinarski Arhiv 2008, 78, 179–185. [Google Scholar]
  21. Mubita, C.M.; Muma, B.J.; Nalubamba, K.; Pandey, G.S.; Samui, K.; Munyeme, M.; Masahiro, K.; Qui, Y.; Saasa, N.; Hang’ombe, B.M. Characterization of Salmonella isolated from domestic animals and wildlife from selected areas of Zambia. Sci. Afr. J. 2020. [Google Scholar] [CrossRef]
  22. Mshana, S.E.; Matee, M.; Rweyemamu, M. Antimicrobial resistance in human and animal pathogens in Zambia, Democratic Republic of Congo, Mozambique and Tanzania: An urgent need of a sustainable surveillance system. Ann. Clin. Microbiol. Antimicrob 2013, 12, 28. [Google Scholar] [CrossRef] [Green Version]
  23. Johnson, J.R. Shigella and Escherichia coli at the crossroads: Machiavellian masqueraders or taxonomic treachery? J. Med. Microbiol. 2000, 49, 583–585. [Google Scholar]
  24. Simulundu, E.; Lubaba, C.H.; Heerden, J.V.; Kajihara, M.; Matta, L.; Chambaro, H.M.; Sinkala, Y.; Munjita, S.M.; Munang’andu, H.M.; Nalubamba, K.S.; et al. The epidemiology of African swine fever in “nonendemic” regions of Zambia (1989–2015): Implications for disease prevention and control. Viruses 2017, 9, 236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Nalubamba, K.S.; David, S.; Munyeme, M.; Kamboyi, H.K.; Saasa, N.; Mkandawire, E.; Munang’andu, H.M. Theileria spp. in free ranging giraffes (Giraffa camelopardalis) in Zambia. J. Vet. Med. Res. 2015, 2, 1039. Available online: https://www.researchgate.net/publication/311107984 (accessed on 26 May 2019).
  26. Kamboyi, H.K.; Garine-Wichatitsky, M.D.; Hang’ombe, B.M.; Munyeme, M. Risk mapping and eco-anthropogenic assessment of anthrax in the upper Zambezi basin. Vet. Med. Sci. 2019, 5, 419–427. [Google Scholar] [CrossRef]
  27. Wiethoelter, A.K.; Beltrán-Alcrudo, D.; Kock, R.; Mor, S.M. Global trends in infectious diseases at the wildlife livestock interface. Proc. Natl. Acad. Sci. USA 2015, 112, 9662–9667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Gortázar, C.; Ferroglio, E.; Höfle, U.; Frölich, K.; Vicente, J. Diseases shared between wildlife and livestock: A European perspective. Eur, J. Wildl. Res. 2007, 53, 241–256. [Google Scholar] [CrossRef]
  29. Fischer, J.R.; Gerhold, R. Wildlife as a risk factor in animal health and zoonoses. Conf. OIE 2002, 273–280. Available online: https://pdfs.semanticscholar.org/1590/5c1e9c2e64edc7b0dd7422b99052ebfe74ba.pdf (accessed on 27 May 2019).
  30. Arnold, K.E.; Williams, N.J.; Bennett, M. ‘Disperse abroad in the land’: The role of wildlife in the dissemination of antimicrobial resistance. Biol. Lett. 2016, 12, pii:20160137. [Google Scholar] [CrossRef] [Green Version]
  31. Mpundu, P.; Mbewe, A.R.; Muma, J.B.; Zgambo, J.; Munyeme, M. Evaluation of bacterial contamination in dressed chickens in Lusaka abattoirs. Front. Public Health 2019, 7, 19. [Google Scholar] [CrossRef] [PubMed]
  32. Watkins, R.R.; Bonomo, R.A. β-Lactam Antibiotics. In Infectious Diseases, 4th ed.; Cohen, J., Powderly, W.G., Opal, S.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2017; Volume 2, pp. 1203–1216.e2. [Google Scholar] [CrossRef]
  33. Alcaide, B.; Aragoncillo, C.; Almendros, P. ‘Cephalosporins’ in Katritzky AR, Ramsden CA, Scrive EFV, Taylor RJK (eds.-in-chief) Comprehensive Heterocyclic Chemistry III. Elsevier Sci. 2008, 2, 111–171. [Google Scholar] [CrossRef]
  34. WHO Critically Important Antimicrobials for Human Medicine, 5th rev. 2016. Available online: http://who.int/foodsafety/publications/antimicrobials-fifth/en/ (accessed on 25 May 2019).
  35. OIE List of Antimicrobial Agents of Veterinary Importance. 2015. Available online: http://www.oie.int/fileadmin/Home/eng/Our_scientific_expertise/docs/pdf/Eng_OIE_List_antimicrobials_May2015.pdf (accessed on 25 May 2019).
  36. Kabwe, M.; Tembo, J.; Chilukutu, L.; Chilufya, M.; Ngulube, F.; Lukwesa, C.; Kapasa, M.; Enne, V.; Wexner, H.; Mwananyanda, L.; et al. Etiology, antibiotic resistance and risk factors for neonatal sepsis in a large referral center in Zambia. Pediatr Infect. Dis. J. 2016, 35, e191–e198. [Google Scholar] [CrossRef] [Green Version]
  37. Chame, M. Terrestrial mammal feces: A morphometric summary and description. Mem. Inst. Oswaldo Cruz 2003, 98, 71–94. [Google Scholar] [CrossRef]
  38. Tsen, H.Y.; Lin, C.K.; Chi, W.R. Development and use of 16S rRNA gene targeted PCR primers for the identification of Escherichia coli cells in water. J. Appl. Microbiol. 1998, 85, 554–560. [Google Scholar] [CrossRef]
  39. Tamura, K.; Stecher, G.; Paaterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  40. Clinical and Laboratory Standards Institute. Performance Standards for Antimicrobial Susceptibility Testin, 28th ed.; CLSI supplementM100; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2018. [Google Scholar]
  41. Kozak, G.K.; Boerlin, P.; Janecko, N.; Reid-Smith, R.J.; Jardine, C. Antimicrobial resistance in Escherichia coli isolates from swine and wild small mammals in the proximity of swine farms and in natural environments in Ontario, Canada. Appl. Environ. Microbiol. 2009, 75, 559–566. [Google Scholar] [CrossRef] [Green Version]
  42. Magiorakos, A.P.; Srinivasan, A.; Carey, R.B.; Carmeli, Y.; Falagas, M.E.; Giske, C.G.; Harbarth, S.; Hindler, J.F.; Kahlmeter, G.; Olsson-Liljequist, B.; et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: An international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 2012, 18, 268–281. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Phylogenetic tree. The 16S rRNA gene neighbor joining phylogenetic tree based on 500 bp nucleotide sequences from Escherichia coli and Shigella isolated from wildlife and domestic animals. The tree was constructed using mega 6 using 1000 bootstrap replicates as confidence interval. The color codes represent origin from the same animal species. The isolates under study were separated into two clusters namely A and B, with cluster A having only 2 isolates from wild animals (buffalo and ostrich) while cluster B was divided into minor clusters B1 and B2 with B2 further dividing into B2a and B2b and most isolates were in this cluster.
Figure 1. Phylogenetic tree. The 16S rRNA gene neighbor joining phylogenetic tree based on 500 bp nucleotide sequences from Escherichia coli and Shigella isolated from wildlife and domestic animals. The tree was constructed using mega 6 using 1000 bootstrap replicates as confidence interval. The color codes represent origin from the same animal species. The isolates under study were separated into two clusters namely A and B, with cluster A having only 2 isolates from wild animals (buffalo and ostrich) while cluster B was divided into minor clusters B1 and B2 with B2 further dividing into B2a and B2b and most isolates were in this cluster.
Antibiotics 10 00238 g001
Figure 2. The location of the study site. Zambia is a landlocked country located in the Central Southern Africa (A); The study site is a game ranch located about 40 km north of Lusaka, Zambia’s capital city (B); Fresh fecal samples were collected from sites within 10 km radius of the main office at the game ranch and GPS coordinates were recorded immediately after sample collection (C).
Figure 2. The location of the study site. Zambia is a landlocked country located in the Central Southern Africa (A); The study site is a game ranch located about 40 km north of Lusaka, Zambia’s capital city (B); Fresh fecal samples were collected from sites within 10 km radius of the main office at the game ranch and GPS coordinates were recorded immediately after sample collection (C).
Antibiotics 10 00238 g002
Table 1. Microbiological, biochemical, and molecular identification of isolates. Escherichia coli and closely related isolates from selected animal species as identified using an API 20E Gram-negative Microbial Identification Kit, and molecular identification methods (PCR and Sequencing) using 16E1/16E2/16E3 and P3mod/P5 Primers.
Table 1. Microbiological, biochemical, and molecular identification of isolates. Escherichia coli and closely related isolates from selected animal species as identified using an API 20E Gram-negative Microbial Identification Kit, and molecular identification methods (PCR and Sequencing) using 16E1/16E2/16E3 and P3mod/P5 Primers.
Animal SpeciesNumber of Samples ProcessedTotal Isolates
(% Yield)
Identified Bacteria by Species
Biochemical Identification Using API® 20E
(n = 66)
Molecular Detection Using 16E1/16E2/16E3 and P3mod/P5 Primers
(n = 65)
Escherichia coli Group 1
n (%)
Burkholderia cepacian
(%)
Kluyvera spp.
n (%)
Unidentified
n (%)
Escherichia coli spp.
n (%)
Escherichia fergusoniin
(%)
Shigella spp.
n (%)
Enterobacteriaceae bacteriumn
(%)
Impala1615 (93.8)6 (40)01 (6.7)8 (53.3)11 (78.6)03 (21.4)0
Buffalo138 (61.5)5 (62.5)003 (37.5)6 (75)1 (12.5)1 (12.5)0
Ostrich76 (85.7)5 (83.3)001 (16.7)5 (83.3)01 (16.7)0
Cattle2016 (80)11 (68.8)005 (31.3)13 (81.2)03 (18.8)0
Pig1815 (83.3)11 (73.3)004 (26.7)10 (66.7)04 (26.7)1 (6.7)
Goat106 (60)1 (16.7)1 (16.7)04 (66.7)3 (50)03 (50)0
Total8466 (78.6)39 (59.1)1 (1.5)1 (1.5)25 (37.9)48 (72.7)1 (1.5)15 (22.7)1 (1.5)
Sequence percentage homology
(%)
Minimum90.001st Quantile98.10
Mean98.29Median99.00
Maximum100.003rd Quantile99.30
Table 2. Relative frequencies of isolates exhibiting resistance (R) or intermediate susceptibility (I) to the tested antimicrobials. A dash (–) denotes neither resistance nor intermediate sensitivity was detected (n = 63).
Table 2. Relative frequencies of isolates exhibiting resistance (R) or intermediate susceptibility (I) to the tested antimicrobials. A dash (–) denotes neither resistance nor intermediate sensitivity was detected (n = 63).
Antimicrobial AgentAMR DetectedAll Isolates
N = 63
Impala
N = 15
Buffalo
N = 8
Ostrich
N = 6
Cattle
N = 16
Pig
N = 15
Goat
N = 6
E. coli
spp.
(n = 48)
n (%)
Shigella
spp.
(n = 15)
n (%)
E. coli
spp.
(n = 11)
n (%)
Shigella
spp.
(n = 3)
n (%)
E. coli
spp.
(n = 6)
n (%)
Shigella
spp.
(n = 1)
n (%)
E. coli spp.
(n = 5)
n (%)
Shigella
spp.
(n = 1)
n (%)
E. coli spp.
(n = 13)
n (%)
Shigella
spp.
(n = 3)
n (%)
E. coli
spp. (n = 10)
n (%)
Shigella
spp.
(n = 4)
n (%)
E. coli
spp. (n = 3)
n (%)
Shigella
spp.
(n = 3)
n (%)
AmpicillinR17 (27)3 (4.8)6 (42.9)---3 (50)1 (16.7)7 (43.8)-1 (6.7)1 (6.7)-1 (16.7)
I10 (15.9)3 (4.8)2 (14.3)--1 (12.5)--2 (12.5)1 (6.2)5 (33.3)-1 (16.7)1 (16.7)
Amoxicillin + Clavulanic AcidR5 (7.9)-4 (28.6)---1 (16.7)-------
I2 (3.2)-1 (7.1)---1 (16.7)-------
CeftazidimeR9 (14.3)1 (1.6)2 (14.3)---1 (16.7)-5 (31.2)--1 (6.7)1 (16.7)-
I4 (6.3)-------1(6.2)-2 (13.3)-1 (16.7)-
CefotaximeR6 (9.5)2 (3.2)2 (14.3)---1(16.7)1(16.7)3 (18.8)--1 (6.7)--
I9 (14.3)1 (1.6)5 (35.7)-----2 (12.5)-2(13.3)--1 (16.7)
CeftriaxoneR4 (6.3)-2 (14.3)---1 (16.7)-1 (6.2)-----
I7 (11.1)2 (3.2)--1 (12.5)--1 (16.7)4 (25)2 (12.5)1 (6.7)---
KanamycinR6 (9.5)3 (4.8)1 (7.1)---2 (33.3)1 (16.7)--1 (6.7)1 (6.7)2 (33.3)1 (16.7)
I20 (31.7)5 (7.9)5 (35.7)1 (7.1)3 (37.5)1 (12.5)--6 (37.5)1 (6.2)4 (26.7)2 (13.3)--
AmikacinR--------------
I--------------
GentamicinR--------------
I1 (1.6)-------1 (6.2)-----
ImipenemR1 (1.6)-1 (7.1)-----------
I2 (3.2)-1 (7.1)-----2 (12.5)-----
Nalidixic AcidR3 (4.8)1(4)1 (7.1)-----2 (12.5)-1 (6.7)---
I1 (1.6)-1 (7.1)-----------
CiprofloxacinR--------------
I1 (1.6)-------1 (6.2)-----
Trimethoprim + SulfamethoxazoleR--------------
I--------------
ChloramphenicolR2 (3.2)1 (1.6)-------1 (6.2)1 (6.7)-1 (16.7)-
I--------------
Table 3. Distribution of relative frequencies of isolates exhibiting resistance to various antimicrobials. The categories of antimicrobials included in this study are penicillins, cephalosporins, aminoglycosides, quinolones, phenicols, folate pathway inhibitors/sulfonamides, and carbapenems. An isolate was considered multidrug-resistant if it showed resistance to at least one antimicrobial agent each among three or more antimicrobial categories (n = 63).
Table 3. Distribution of relative frequencies of isolates exhibiting resistance to various antimicrobials. The categories of antimicrobials included in this study are penicillins, cephalosporins, aminoglycosides, quinolones, phenicols, folate pathway inhibitors/sulfonamides, and carbapenems. An isolate was considered multidrug-resistant if it showed resistance to at least one antimicrobial agent each among three or more antimicrobial categories (n = 63).
Number of Antimicrobial ClassesAll IsolatesImpalaBuffaloOstrichCattlePigGoat
E. coli
spp.
(n = 48)
n (%)
Shigella
spp.
(n = 15)
n (%)
E. coli
spp.
(n = 11)
n (%)
Shigella
spp.
(n = 3)
n (%)
E. coli
spp.
(n = 6)
n (%)
Shigella
spp.
(n = 1)
n (%)
E. coli
spp.
(n = 5)
n (%)
Shigella
spp.
(n = 1)
n (%)
E. coli
spp.
(n = 13)
n (%)
Shigella
spp.
(n = 3)
n (%)
E. coli
spp.
(n = 10)
n (%)
Shigella
spp.
(n = 4)
n (%)
E. coli
spp.
(n = 3)
n (%)
Shigella
spp.
(n = 3)
n (%)
14 (8.4)-2 (18.2)-------2 (20)---
210 (20.9)1 (6.7)1 (9.1)---2 (40)-6 (46.2)-1 (10)--1 (33.3)
38 (16.7)1 (6.7)3 (27.3)---1 (20)1 (100)2 (15.4)-1 (10)-1 (33.3)-
41 (2.1)1 (6.7)---------1 (25)1 (33.3)-
5--------------
6--------------
7--------------
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kabali, E.; Pandey, G.S.; Munyeme, M.; Kapila, P.; Mukubesa, A.N.; Ndebe, J.; Muma, J.B.; Mubita, C.; Muleya, W.; Muonga, E.M.; et al. Identification of Escherichia coli and Related Enterobacteriaceae and Examination of Their Phenotypic Antimicrobial Resistance Patterns: A Pilot Study at A Wildlife–Livestock Interface in Lusaka, Zambia. Antibiotics 2021, 10, 238. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10030238

AMA Style

Kabali E, Pandey GS, Munyeme M, Kapila P, Mukubesa AN, Ndebe J, Muma JB, Mubita C, Muleya W, Muonga EM, et al. Identification of Escherichia coli and Related Enterobacteriaceae and Examination of Their Phenotypic Antimicrobial Resistance Patterns: A Pilot Study at A Wildlife–Livestock Interface in Lusaka, Zambia. Antibiotics. 2021; 10(3):238. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10030238

Chicago/Turabian Style

Kabali, Emmanuel, Girja Shanker Pandey, Musso Munyeme, Penjaninge Kapila, Andrew Nalishuwa Mukubesa, Joseph Ndebe, John Bwalya Muma, Charles Mubita, Walter Muleya, Elizabeth Muligisa Muonga, and et al. 2021. "Identification of Escherichia coli and Related Enterobacteriaceae and Examination of Their Phenotypic Antimicrobial Resistance Patterns: A Pilot Study at A Wildlife–Livestock Interface in Lusaka, Zambia" Antibiotics 10, no. 3: 238. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics10030238

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop