Article Text


Predictors and risk factors for the intestinal shedding of Escherichia coli O157 among working donkeys (Equus asinus) in Nigeria
  1. Jesse T. Jedial1,
  2. Aminu Shittu2,
  3. Faruk M. Tambuwal1,
  4. Mikail B. Abubakar1,
  5. Muhammed K. Garba3,
  6. Jacob P. Kwaga4 and
  7. Folorunso O. Fasina5
  1. 1Department of Pathology and Microbiology, Faculty of Veterinary Medicine, Usmanu Danfodiyo University, Sokoto, Nigeria
  2. 2Department of Public Health and Preventive Medicine, Faculty of Veterinary Medicine, Usmanu Danfodiyo University, Sokoto, Nigeria
  3. 3School of Medical Laboratory Sciences, College of Health Sciences, Usmanu Danfodiyo University, Sokoto, Nigeria
  4. 4Department of Public Health and Preventive Medicine, Faculty of Veterinary Medicine, Ahmadu Bello University, Zaria, Nigeria
  5. 5Department of Production Animal Studies, University of Pretoria, Pretoria, Onderstepoort, South Africa
  1. Correspondence to Dr Aminu Shittu; ameen_vet{at}


Objectives Escherichia coli are an important group of bacteria in the normal gastrointestinal system but can sometimes cause infections in domestic animals and man. Donkeys are routinely used as multipurpose animal but details of burdens of potentially infectious bacteria associated with it are limited. The prevalence and associations between intestinal shedding of E. coli O157 and animal characteristics and management factors were studied among 240 randomly selected working donkeys in north-western Nigeria.

Design Four local government areas, of Sokoto State in north-western Nigeria were recruited in this study. A multistage randomised cluster design was used to select subjects and donkey owners within selected zones. Confirmation of infection was based on bacterial culture, isolation and biochemical test for E. coli O157 from faecal samples.

Results Of the total bacteria isolated, 203 of the 329 (61.70 per cent) were E. coli, 76 of which was E. coli serotype O157. A multivariable logistic regression model was used to examine the relation between intestinal shedding of E. coli O157 and selected variables. The analysis yielded five potential predictors of shedding: soft faeces in donkeys, Akaza and Fari ecotypes of donkey were positive predictors while maize straw as feed and sampling during the cold dry period were negative predictors.

Conclusions This study concludes that controlling intestinal shedding of E. coli O157 among working donkeys in Nigeria is possible using the identified predictors in planning appropriate interventions to reduced human risk of infection.

Statistics from


Escherichia coli are a group of intestinal microflora in human beings and animals that are usually harmless (Greenland and others 2009). While many of the E. coli gut microflora are innocuous, their build-up due to gut stasis, disruptions of the intestinal activities or a sudden change of diet can enhance the accumulation of certain toxins associated with these organisms and cause disease conditions. Some strains of E. coli particularly produce powerful toxins that can cause intestinal or extraintestinal diseases (Kaper and others 2004). These strains include the verotoxin-producing E. coli (VTEC) also known as shiga toxin-producing E. coli (STEC) (Bettelheim and Beutin 2003). The STEC that causes haemorrhagic colitis and haemolytic uraemic syndrome is called enterohemorrhagic E. coli (Nataro and Kaper 1998), and it is recognised as an important foodborne pathogen (Gyles 2007). Even healthy animals can harbour human enteric pathogens, many of which have a low infectious dose (Bell and others 1994). Although cattle are the main reservoirs of human pathogenic VTEC, there is evidence that sheep, deer, dogs, poultry and goats can also carry the VTEC strains (Callaway and others 2006, Reinstein and others 2007, Greenland and others 2009). E. coli organisms are also classified into serogroups based on the heat-stable ‘O’ and heat-labile surface ‘K’ or flagellar ‘H’ antigens (Vosti and others 1964). The E. coli O157:H7 serogroup is the most important cause of severe foodborne illnesses in living organisms and severe infection can result in case fatality of up to 50 per cent. This serogroup among others carries the shigatoxin 1and 2 and the gene responsible for effacement (eae gene) (Gannon and others 1993).

The primary mode of transmission for E. coli is the faecal-oral route (contaminated food, milk and water), but there are other possible means of transmission, because animal fur, hair, skin and saliva often harbour faecal material with the infective organism (Keen and Elder 2002, Varma and others 2003, De Schrijver and others 2008).

Donkeys (Equus africanus asinus) belong to the family Equidae. They are especially useful in arid and semiarid locations in sub-Saharan Africa, because of their hardiness and ability to survive where oxen and most other animals do not thrive. In Nigeria, the population of working donkeys is estimated at over one million; they are classified according to the ecoclimatic conditions and most of them are distributed in the extreme northern region (Blench and others 1992, Pearson and others 1999). They are particularly useful for transport, pulling carts, farm tillage, threshing, fetching/carrying water, milling and other energy-intensive activities, especially in the rural and semiurban locations where road networks are unavailable to motorised vehicles. An adult donkey weighs between 90 kg and 210 kg depending on the breed. Donkeys depend on a low calorie fibre diet (e.g. straw, husks, hay, fresh fodder and dried grasses) for most of their lives (Pearson and others 1999).

Donkeys are known to harbour certain diseases, including trypanosomiasis, babesiosis, African horse sickness, the equine herpes virus, equine influenza, rabies, horse pox, mange and glanders, but their importance in transmitting some zoonotic diseases is not well known. Given the current effort to promote the use of donkeys in drier areas of the world and consequent potential for increased human-animal (donkey) interaction, it is vital to understand donkeys’ susceptibility and resistance to disease, and the zoonotic potential of donkey diseases (Pearson and others 1999). Previous study on animal-originated human VTEC O157 infection have shown that only cattle and sheep play more role than donkeys in the epidemiology of the disease (Pritchard and others 2009).

To the authors’ knowledge, no prior study has investigated predictors associated with intestinal shedding of E. coli O157 among working donkeys in Nigeria and few data exist elsewhere. In addition, there is little information on the diseases affecting donkeys based on data from the veterinary clinics in the country. The objective of this study therefore was to estimate the prevalence of E. coli O157 among working donkeys in parts of Nigeria and to identify animal characteristics and husbandry practices that could be potential predictors of intestinal shedding of E. coli O157.

Materials and methods

Study area and sampled population

Sokoto is a state in north-western Nigeria. Geographically, the state lies between the longitudes 4°0′E and 6°54′E and latitudes 12°N and 13°58N. The main occupation of the people in the state is arable farming and rearing livestock. Sokoto State has the second-largest livestock population in Nigeria, with an estimated 3 million cattle, 3 million sheep, 5 million goats, 4600 camels, 52,000 donkeys and hosts of local and exotic poultry species. The State consists of 23 local government areas (LGAs) and is broadly divided into two agricultural zones, namely the northern (comprising of 12 LGAs) and the western (11 LGAs) zones (Fig 1).

FIG 1:

Map of Sokoto State showing boundaries (national and international), agricultural zones and selected (local government areas) LGAs for the study

This study was conducted between May 2009 and April 2010 in Sokoto State previously described above. To select the study sample, simple balloting was used to pick LGAs in the stratified agricultural zones. Four LGAs were selected from each zone, namely Wurno, Isah, Illela and Tangaza from the northern zone, and Tambawal, Bodinga, Yabo and Sokoto North from the western zone. A multistage randomised cluster design sampling method was then used to select subjects and donkey owners within each of the selected LGAs. Sample collections were done principally to target market days for each of the selected LGAs.

Selection of participants, questionnaire administration, data retrieval and sampling

Four research assistants were recruited and trained in questionnaire administration, retrieval, assessment of scoring criteria and sample collection in a prestudy orientation. Severity of loss of body condition in the donkeys was assessed using the scorecard according to Pearson and Ouassat (2000) as follows: (1) Severe loss of body condition (very thin to less thin); (2) Mild loss of body condition (less than moderate to just moderate); and (3) No loss of body condition (less fat to very fat) and the present health conditions were scored based on the current observations/presenting conditions mentioned by the owners/handlers.

A pilot study was conducted using seven donkeys and their owners/handlers to test and validate the questionnaire and procedures for the data and samples required. For the main study, 30 donkey owners/handlers who had a minimum of one working donkey each were included from each selected LGAs, adding up to a total of 120 (30 owners/handlers×4 LGAs) individuals per zone. Since the two zones were included in the study, a total of 240 closed-ended structured questionnaires were administered to selected donkey owners/handlers on the spot during the visitations in each of the four selected LGAs.

The questionnaire was used to collect data comprising 21 variables with 55 options and it focused on management and husbandry practices relating to the working donkeys. A comprehensive list of variables that have been thought to predict E. coli O157 among working donkeys were analysed (Table 1), and the descriptive result on total number of donkeys positive for E. coli O157 serotype is shown in Table 1. Since an on-the-spot assessment was conducted, all the 240 donkey owners and handlers contacted responded to the interview (100 per cent return rate) and no respondent had more than 1 donkey. Faecal samples were collected from each donkey at the same time that the questionnaire was being administered, and these samples were taken directly from the rectum of each donkey (n=240), and were transported to the Microbiology Laboratory, Department of Medical Laboratory Sciences, College of Health Sciences, Usmanu Danfodiyo University Sokoto, Nigeria, immediately for processing.


Donkey ecotypes, sex, definition of potential predictors of intestinal shedding of Escherichia coli O157 and prevalence of Escherichia coli O157 among working donkeys in Nigeria

Microbiological and biochemical tests were conducted to characterise the bacterial organisms according to standard techniques (Klein and others 2002, Voetsch and others 2004, Kotgire 2012). Briefly, faeces were observed for consistency, presence of blood, mucous, worms and colour. Cellular exudates were checked for using methylene blue as described by Kotgire (2012). Wet mount were carried out using the hanging drop method. Prepared faecal samples were plated on MacConkey's agar and a selective medium (Sorbitol MacConkey's agar) to identify E. coli O157. Further screening for E. coli O157 was performing using an antigen-specific latex-agglutination test (Rosario and others 1991). Although few other bacterial organisms were picked alongside the E. coli O157, further identification was not carried out on them since they fell outside the objective of this study.

Statistical analysis

Statistical analysis was performed using Stata version 10.0 by StataCorp, Lakeway Drive, College Station, Texas, USA. An initial descriptive analysis of the proportions of faecal samples confirmed positive for E. coli O157 in donkeys, stratified by animal and management explanatory variables (Table 1), was conducted in order to identify variables that might be of value for further investigation. Logistic regression analysis was performed with intestinal shedding of E. coli O157 (yes=1/no=0) as the observed outcome (dependent variable) based on the microbiological and biochemical test results. Only variables that had unconditional associations with the outcome that were significant using a probability value of P<0.25 in the univariable analysis (Table 2) were included for further analysis using a multivariable logistic regression model (Dohoo and others 2009). The correlation among explanatory variables was first checked using a multicollinearity analysis and the observation of mean variance inflation factors (VIF) (Table 3). A preliminary model was developed using stepwise backward removal of variables. Differences between nested models were assessed using likelihood ratio tests. Variables were retained in the multivariable model if they improved model fit significantly (P≤0.05), using the Hosmer-Lemeshow goodness-of-fit (GOF) χ2 test. Results are reported as crude OR (ORc) in univariable analysis and adjusted OR (ORa) in multivariable analysis with 95% CIs.


Results of the univariable logistic regression analysis for associations between the intestinal shedding of Escherichia coli O157 and various animal and management factors for 240 samples


Final multivariable logistic regression model of predictors and risk factors associated with the intestinal shedding of Escherichia coli O157 among working donkeys in Nigeria for 240 faecal samples


Bacteriological analysis

A total of 235 faecal samples were positive for the presence of one or more bacterial organisms (97.92 per cent) by culture and specifically E. coli was identified in 202 samples representing 84.58 per cent of the total faecal samples (240). Only 31.67 per cent (76/240) of the faeces collected from the seven different ecotypes of working donkeys were positive for E. coli O157.

Results of univariable logistic regression analysis

No evidence of multicollinearity amongst any of the categorical variables was identified (mean VIF=1.23, Table 3). Of the 48 potential predictor variables tested in a χ2 univariable analysis, only 11 variables (with bold font P values, Table 2) showed potential association (P<0.25) with intestinal shedding of E coli O157 and these were considered for inclusion in the multivariable analysis (Table 2).

Using the stepwise backward elimination procedure, only five variables were retained in the final model (Table 3). The Fari (ORa=3.40, P=0.05, 95% CI 0.98 to 11.78) and Akaza (ORa=2.77, P=0.02, 95% CI 1.16 to 6.66) ecotypes of donkey were the most susceptible to the risk of intestinal shedding of E. coli O157 compared with other ecotypes in this study (Table 3). Similarly, soft faecal consistency is the most important predictor for intestinal shedding of E. coli O157 (ORa=4.16, P=0.002, 95% CI 1.70 to 10.19). Contrastingly, both the feeding of maize straw (ORa=0.52, P=0.05, 95% CI 0.27 to 0.99) and sampling during the cold dry period (harmattan) in Nigeria (ORa=0.53, P=0.05, 95% CI 0.28 to 1.00) were negative predictors for intestinal shedding of E. coli O157 (Table 3).

The logistic regression χ2 was 22.75 with 4 degrees of freedom and the Hosmer-Lemeshow GOF χ2 was 2.59; P=0.63, an indication of the good model fit to the data.


E. coli O157 serotypes were shed by a moderately large proportion of the donkeys (31.67 per cent). This organism may serve as a contaminant or it may have direct or indirect zoonotic implications. The authors isolated certain other bacteria from donkey faeces but these were not the main focus of the present study. Derlet and Carlson (2002) previously confirmed that donkeys are reservoirs of certain zoonotic organisms. Since these animals are used for work regularly and the handlers may eat between tasks, sometimes without thorough washing of their hands or observation of strict hygienic measures, they stand the risk of inadvertent exposure to these organisms from donkey. Interestingly, donkey faeces is sometimes used in rural communities to rub/coat the inner walls of mud buildings where human beings live, and this creates a strong potential for contamination and/or infection of those persons who perform this work and people who live in and touch these surfaces (Derlet and Carlson 2002, Pritchard and others 2009). In addition, there is a huge risk of water and environmental contamination since these animals are reared extensively.

Significantly, the predictors that were positively associated with intestinal shedding of E. coli O157 included soft faecal consistency, the Fari and Akaza ecotypes of donkey (Tables 3 and 4). Since these Fari and Akaza ecotypes of donkeys seem to have greater susceptibility compared with the other ecotypes, perhaps these ecotypes are not as adapted to the highly fibrous feeds that are fed to these animals in Nigeria or there are intrinsic genetic differences that are yet to be studied. However, there is no proof to support these postulations. It may be possible also that some unknown innate factors predispose these ecotypes to more infection by E. coli O157. Since these animals are fed with the largely fibrous diet which is known to increase intestinal transit time in animals with probable increased levels of intestinal E coli O157 shedding and associated colienteritis-colisepticaemia, there may be severe diarrhoea (Harvey and others 1973). As such, subclinical diarrhoea may be a consequence of the feed type and E. coli O157 build-up in the intestine rather than a direct predictor of the organism. It should be noted, however, that the shedding of E. coli in faeces does not always correlate with clinical disease since animals carry this organism regularly.


Correlation matrix for the set of explanatory variables in the analysis of potential predictors of intestinal shedding of Escherichia coli O157 among working donkeys in Nigeria

Finally, sampling during the harmattan-cold dry period (which is usually associated with a surge in animal and human diseases) reduced the risk of shedding of E. coli O157 by about half. The reason for this association cannot be immediately established but it is known that during the dry cold harmattan period, free-range donkeys have access to less water, often pass hard dry faeces which are further subjected to the desiccating effects of the harsh environment. It is instructive that faecal sampling for prevalence study in animal within the country may not reveal a true prevalence if collected only during the harmattan period (dry, dusty and cold north-easterly trade wind).

In the course of the interviews, the authors observed that all the donkeys sampled are kept under a traditional (extensive) system of management despite their enormous contributions to the livelihood of the rural dwellers. This husbandry practice is probably not much different from what exists in other parts of Nigeria and indeed Africa (exposure of these donkeys to inclement weather, diseases and nutritional deficiencies) (Research Inventory Management 1992). The generally held perception that donkeys are very hardy and will withstand adverse conditions, and are resistant to diseases without a drop in outputs may underlie this observed poor management system.

More male donkeys were observed and sampled because they are more frequently owned. The owners’ preference for male donkeys is linked to their perception of comparative greater strength as draught animals and probably because female donkeys do not work satisfactorily particularly when they are pregnant (more so, in their last trimester). In this analysis, the authors found out that female donkeys are more likely to shed E. coli O157 compared with male donkeys. Keeping of donkeys or with other animals feeding them together does not influence the shedding of E. coli O157.

This work has some limitations. While the E. coli serotypes were isolated, the isolates were not checked for shigatoxin principles and or the eae gene. The authors are aware that these factors are important in the virulence properties of E. coli O157 and its bacteria zoonotic principles. The authors acknowledge that the outcome of this analysis is based on the limited sampled population and small geographical spread. Future studies may benefit from using geographically more diverse samples from different parts of the West African subregion. The test systems may also present with a degree of non-specificities and poor sensitivities. It is in the authors’ opinion though that the multiple microbiological tests applied in this study will reduce errors/misclassifications due to this effect. Nevertheless, this study has confirmed the presence of some potential infectious and zoonotic pathogens and contaminants in working donkeys in Nigeria and has associated these pathogens and contaminants with particular predictors thus opening up opportunities for further research in this field.


The authors thank local owners of working donkeys for allowing them to sample their donkeys and for assisting them with the animal restraint.

  • Received July 14, 2014.
  • Revision received May 1, 2015.
  • Accepted May 5, 2015.


View Abstract


  • Twitter Follow Aminu Shittu at @ameenvet

  • Contributors JTJ, AS, FMT and FOF designed the study. JTJ, MBA, MKG and AS collected processed samples. AS and FOF conducted statistical analysis and interpretation of data. FMT, MBA, FOF and JPK coordinated the study. All authors contributed to drafting the manuscript.

  • Funding This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.