INTRODUCTION
Bacterial resistance to antimicrobial agents now ranks as one of the most important global health concerns (
1), particularly since the development of new antimicrobials has stalled over recent years (
2). Antimicrobials are extensively used in food animal production to prevent and treat animal diseases, as well as to promote growth (
3). Antimicrobials used in animal farming and human medicine are largely similar (
4). In Vietnam, extensive use of antimicrobials in farming has been documented, including compounds considered to be critical for human medicine, such as fluoroquinolones, macrolides, and aminoglycosides (
5,
6). As a result, high levels of antimicrobial resistance (AMR) in zoonotic bacteria, such as nontyphoid
Salmonella and
Campylobacter spp., have been reported both on farms and in meat products (
7–9).
The use of antimicrobials on farms raising domestic animals for food is suspected to play an important role in generating/disseminating AMR organisms (organisms that demonstrate AMR) and determinants through environmental pathways (
10). Humans may become infected with AMR bacteria through occupational exposure and ingestion of contaminated food and water (
11). In addition, antimicrobials used on farms have been reported to increase AMR in organisms isolated from areas exposed to farm waste (
12), as well as in the flora of wildlife trapped in their proximity (
13,
14).
Escherichia coli organisms are generally commensal bacteria of warm-blooded animals. They are often used as indicators for the presence of AMR in monitoring and surveillance programs, given the availability of simple and efficient isolation procedures (
15,
16).
To date, no studies have investigated the prevalence of AMR in E. coli strains isolated from the most common farmed and wildlife species in the Mekong Delta of Vietnam. The Mekong Delta is a region criss-crossed by a maze of canals and rivers, with rice as the predominant agricultural crop. In addition, the area has a high density of poultry and pigs farmed under conditions of low biosecurity and biocontainment. We performed a survey of AMR in commensal E. coli isolates from the three most commonly farmed species (chickens, ducks, and pigs) and small wild mammals trapped in the same district. The aims were (i) to describe AMR phenotypic profiles, (ii) to investigate risk factors for AMR (farms), and (iii) to investigate the level of similarity between profiles of E. coli from farmed animals and wildlife.
MATERIALS AND METHODS
Farm survey design.
The survey was conducted between February and May 2012 in the district of Cao Lanh (population, 193,184), located in Dong Thap Province. In the district of Cao Lanh, the median farm sizes (number of animals) are 22 (interquartile range [IQR], 16 to 45), 57 (IQR, 40 to 200), and 4 (IQR, 3 to 8) for chicken, duck, and pig farms, respectively (data provided by the Sub-Department of Animal Health in Dong Thap {lsqb]SDAH-DT]). The selected study area is representative of the Mekong Delta of Vietnam. The area has a tropical monsoon climate, with temperatures ranging from 25 to 28°C with slight annual variations, and a rainy season from May to October.
We aimed to investigate ∼450
E. coli isolates from 90 randomly selected farms for the three main farmed species (pigs, ducks, and chickens; 30 each, i.e., 5 colonies per farm). The purpose was to estimate the probability of a randomly selected
E. coli colony being resistant to a panel of antimicrobials, as well as being multiresistant. This sample size, based on a cluster study design, would theoretically allow the estimation of an AMR prevalence of 25%, with a relative error of 20% and 80% confidence, assuming a design effect of 4 (i.e., a medium level of farm level clustering) (
17). Farms were classified by size (small, medium, and large; total, 9 strata) to maximize the chance of capturing differences between production systems. The cutoffs for animal numbers that determined farm size were as follows: for chickens, 20 to 50 (small), 51 to 100 (medium), and >100 (large); for ducks, 50 to 200 (small), 201 to 1,000 (medium), and >1,000 (large); and for pigs, 5 to 20 (small), 21 to 50 (medium), and >50 (large).
Farm sampling and data collection.
From each farm, naturally pooled feces were collected from the area housing the target species using one pair of boot swabs. These were used to walk for at least 30 steps on areas where fresh droppings were visible. When it was not practicable to use boot swabs (i.e., farrowing crates or chicken houses built on stilts), visible fresh fecal material from at least 10 different points was collected using hand gauze swabs. Immediately after collection, the swab samples were placed in individually sealed containers and transferred to the laboratory at 4°C. Data on farm management practices, experience of disease in the flocks/herds, and other variables were collected using structured questionnaires aimed at the person with primary responsibility for animal husbandry. The questionnaire also enquired about any antimicrobial formulation used over the previous 90 days on the target species (see the supplemental material). Both interviews and fecal sample collection were conducted by trained veterinarians from the SDAH-DT.
Trapping and processing of small wild mammals.
Small wild mammals (rats and shrews) were trapped in different locations (pig and poultry farms, forests, and the perimeters of rice fields) in Cao Lanh District during March 2013 using live-cage metal “tomahawk” traps (∼0.04 cubic meters in volume each). A total of 30 lines (one per location) of 10 traps were set up, checked for animals trapped, and reset over 10 subsequent days (totaling 3,000 trap-days). Traps with a caught animal were collected and replaced by a new trap. The animals were humanely culled by overdose of an inhalant anesthetic (isoflurane) following American Veterinary Medical Association (AVMA) guidelines (
18). After measuring the main biological and morphometric parameters, ∼0.2 g of cecal contents extracted from each animal was mixed with 0.5 ml glycerol and transferred to the laboratory at 4°C. Species identification was carried out by examination of morphological characteristics following existing guidelines (
19). Species identification was confirmed for a representative number of mammals (
n = 22) by amplification of a conserved housekeeping gene (cytochrome oxidase I [COI]) using PCR followed by molecular sequencing of the product and blasting using the RatSEA barcoding tool (
http://www.ceropath.org/barcoding_tool/rodentsea).
Laboratory processing of samples.
All samples were processed in the laboratory within 24 h of collection. A fixed volume of saline solution (100 ml) was added to each swab sample and was then manually shaken. Swab samples (100 μl of the mixture described above) and small wild fecal samples (100 μg of fecal material) were further diluted 1:1,000 in saline solution, and 50 μl of this suspension was plated onto MacConkey agar and incubated at 37°C overnight. From each sample, up to five colonies showing typical E. coli morphology were confirmed using standard biochemical tests (motility, indole, lactose/glucose fermentation, methyl red, citrate, urease, hydrogen sulfide, and gas production).
Antimicrobial resistance of
E. coli isolates was determined by the disk diffusion method using the procedure of the Clinical and Laboratory Standards Institute (CLSI) (
20). The following antimicrobials were used at the concentrations shown: amoxicillin-clavulanic acid (30 μg), ceftazidime (30 μg), ampicillin (10 μg), chloramphenicol (30 μg), tetracycline (30 μg), trimethoprim-sulfamethoxazole (1.25/23.75 μg), gentamicin (10 μg), and ciprofloxacin (5 μg). Isolates were classified as susceptible, intermediate, or resistant according to interpretation of the zone diameter as recommended by CLSI. Potential production of extended-spectrum beta-lactamase (ESBL), as indicated by resistance to ceftazidime and by an inhibitory effect of clavulanic acid, was confirmed using a double disk diffusion test. A reference
E. coli strain (ATCC 25922) was used to verify the quality and accuracy of the testing procedures. Colonies that were intermediate-resistant based on the measured zone were also regarded as resistant. A multidrug-resistant (MDR) strain was defined as a strain resistant to at least three different classes of antimicrobials.
Statistical analyses.
We aimed to model the overall probability of one
E. coli colony from a given farm or from a small wild mammal being AMR. The farm
E. coli data were adjusted for the sampling frame by assigning a stratum-specific sampling weight (
Wi) to each observation unit (i.e., an isolate obtained from a farm) using the following equation:
, where
NT is the total number of pig and poultry farms in Cao Lanh (35,248) and
Ni is the total number of farms of each stratum sampled (
i = 1…9). This adjustment gives more weight to isolates from more common farm types (i.e., small chicken farms). Standard errors were corrected to take into account potential similarities of usage between farms in each stratum (
21). Survey prevalence calculations were carried out using R software (
http://www.r-project.org) with the package “survey.”
The potential statistical association between resistances against any combination of two antimicrobials of eight tested (number of comparisons, 28) was investigated using Pearson's correlation test and was visualized using a correlogram built with the R package “corrgram.”
Potential factors associated with MDR were investigated by hierarchical logistical regression, taking into account the survey sampling frame, with “farm” as the clustering variable, using the R package “lmer.” The variables investigated were main farmed species, farm size, farmer demographics, recent use of antimicrobials (90 days), and type of water and feed. For
E. coli isolates from small mammals, generalized estimating equation (GEE) models were used to calculate the adjusted prevalence of AMR for isolates from animals trapped on farms versus those from animals trapped in forests and rice fields. The species and site of collection (trapped on farms versus trapped at other sites) were modeled as independent variables, and the animal identity was taken as the clustering variable. Model building was carried out in a step-forward fashion from univariable models, and variables were deemed significant if their
P values were <0.05 (
22). Interaction between main significant variables was tested.
Each isolate was characterized by a profile defined by its unique combination of susceptibility/resistance to the 8 antimicrobials tested. The Shannon-Weaver entropy (H) was used to estimate the diversity of the AMR profiles (i.e., the combination of susceptibilities) observed among isolates from each farmed species (chickens, ducks, and pigs), as well as among isolates from each farm separately. H was defined as follows: , where pi is the proportion of isolates with the ith resistance profile and R is the total number of profiles. Higher values of H correspond to greater diversity in AMR profiles.
A minimum spanning tree was plotted to visualize the relationships between the most common AMR profiles using the “adegenet” package in R.
Permutational nonparametric multivariate analysis of variance of profiles using distance matrices (
23) was carried out to investigate the fraction of the total variation in AMR profiles among
E. coli isolates attributable to the host type (pig/chicken/duck/wild mammal trapped on a farm/wild mammal trapped at a nonfarm site). The variability of profiles attributed to the farm was calculated separately for each farmed species. In addition, the relative similarity between profiles of isolates from small wild mammals trapped on farms and at nonfarm sites compared with profiles of isolates from farmed species (chickens, ducks, and pigs) was assessed by analysis of variance. This was performed by fitting two separate models with isolates from each of the three target species: one including data from isolates from farmed species and small wild mammals trapped on farms and another one including isolates from farmed species and small wild mammals trapped on nonfarm sites. The variability was tested by fitting the variable “species” (farmed animal/small wild mammal). Analyses were carried out using the R package “vegan.”
Discriminant analysis of principal components (DAPC) (
24) was used to investigate overall differences between AMR profiles in
E. coli isolates from different sources. Because of the small number of variables compared to the number of observations, no information was discarded in the prior-dimension reduction step. This analysis was carried out using the “adegenet” R package, and the results were visualized using the “ggplot2” R package.
DISCUSSION
Among E. coli isolates from domestic animal species, we observed very high levels (∼80%) of ampicillin and tetracycline resistance, as well as moderate levels (∼25 to 50%) of trimethoprim-sulfamethoxazole, chloramphenicol, amoxicillin-clavulanic acid, and ciprofloxacin resistance. In contrast, the levels of gentamicin and ceftazidime resistance were low (∼8% and <1%, respectively).
Levels of AMR among
E. coli isolates from farms were considerably higher than AMR levels recently reported in
E. coli isolates from meat and farmed animals in developed countries: for example, monitoring data from the European Union (EU) (2012) indicate that 22 to 23% of isolates from broiler meat and pig farms were tetracycline resistant, 20 to 50% were ampicillin resistant (highest in broilers), 23 to 29% were sulfonamide resistant, 2 to 10% were chloramphenicol resistant, and 1.5 to 3.8% were gentamicin resistant (
25).
Pig isolates were much more resistant to ciprofloxacin (41.4%) and gentamicin (38.1%) than isolates from pigs in Europe (5.6% and 1.1%, respectively). In contrast, levels of ciprofloxacin and gentamicin resistance among chicken isolates (21% and 10.8%) were not too dissimilar to levels reported in Europe (29.1% and 3.2%, respectively) (
25). Although levels of AMR and MDR were very high among pig isolates, they were still marginally lower than those reported in Thailand (i.e., 96.2% tetracycline resistant, 91.6% ampicillin resistant, and 98% MDR) (
26).
We have limited quantitative data on the use of antimicrobials in different farmed species on Mekong Delta farms. Because of this, we cannot establish whether the differences in AMR against different antimicrobials or the higher prevalence of AMR in the pig species is the result of differences in antimicrobial consumption. Overall relatively high levels of tetracycline and ampicillin resistance have been reported in other studies and probably reflect the long history of use of both of these antimicrobials as both therapeutic and prophylactic agents (
16).
A previous survey of chicken farms in the nearby province of Tien Giang indicated that over an ∼90-day period, polypeptides (colistin), tetracyclines, and penicillins were the three most commonly used classes of antimicrobials (by 18.6%, 17.5%, and 11.3% of farms, respectively). The study indicated a total usage of 690.4 g (±203.6 g) of antimicrobial compound per chicken produced, which is very high compared with European standards (
6).
A (qualitative) survey in the Red River Delta of Vietnam indicated a higher frequency of antimicrobial usage among pig farms both as growth promoters and for disease prevention and chemotherapy than on chicken farms (
5). In addition, in the Mekong Delta and elsewhere in Vietnam, anecdotal field observations indicate that antimicrobials, such as tetracycline, tylosin, avilamcyin, enramycin, and colistin, are commonly used in feed formulations. Both tetracycline and colistin are important antimicrobials for human medicine. Unfortunately, due to logistical limitations,
E. coli isolates were not tested for colistin resistance in our study.
Neither chloramphenicol nor ciprofloxacin is an antimicrobial approved for use in terrestrial-animal farming in Vietnam (
27). The observed levels of chloramphenicol resistance may partly be explained by the common use of florfenicol, a fluorinated derivative that shows some cross-resistance with chloramphenicol (
28). The observed high levels of resistance against ciprofloxacin (a quinolone antimicrobial not licensed for veterinary use) is probably explained by cross-resistance with other quinolones licensed for veterinary use, such as enrofloxacin and norfloxacin (
29). In contrast, levels of resistance against third-generation cephalosporins (ceftazidime, 0.8% [95% CI, 0.3 to 1.4%]) were very low and comparable with those reported in the EU (cefotaxime) (1.1 to 3.0%) (
25).
A large fraction of farmers (∼90%) reported having used antimicrobials over the previous 3 months. Unfortunately our study did not succeed in gathering detailed information on antimicrobial usage to be able to explain the observed differences in the levels of AMR at the farm level. More precise estimates of antimicrobial usage practices are required to establish this. Obtaining this type of data requires a longitudinal study design and dedicated farmer training in record keeping.
The observed lower levels of AMR on farms run by older farmers is probably a reflection of more experience in farming, resulting in better disease control with fewer antimicrobials. Overall, a high percentage of the total variability (41 to 55%) in AMR profiles was explained by the term “farm,” suggesting that on-farm interventions (likely including antimicrobial usage) may explain the observed AMR profiles. The overall highest diversity was found among isolates from duck flocks. A possible explanation is the characteristic management system of duck flocks that involves daily grazing on communal rice fields, representing a greater opportunity for infection with foreign
E. coli strains (
30).
We found relatively high levels of AMR among isolates from wild mammals and an unequivocally higher prevalence of AMR among isolates from wild mammals trapped in farming locations, suggesting that AMR in wildlife may be driven partly by on-farm antimicrobial usage and/or AMR on farms. These results are consistent with findings from a Canadian study, where levels of AMR among
E. coli isolates from wild mammals trapped on pig farms were considerably higher than AMR levels among isolates from mammals trapped in natural areas (30.9% versus 5.4%, respectively) (
13). A further Canadian study comparing
E. coli isolates from wild mammals trapped in different environments reported the highest levels of AMR among wild mammals trapped on pig farms (48.0%) compared to isolates from wild mammals in residential areas (9.0%) (
31). All our study farms and trapping sites were no less than 5 km away from the closest urban center (the city of Cao Lanh) and 15 km from the closest sewage treatment plant. Therefore, the impact of human sewage from farming populations on the AMR patterns of wildlife flora cannot be separated from the impact of the farmed animals themselves.
AMR may develop in the gut flora of wild mammals when the animals consume antimicrobials present in feed and water on farms, but also indirectly as a result of exposure to farm waste (
32,
33). However, comparisons using pulsed-field profiles between
E. coli isolates from pigs and from rodents trapped on the same farms indicated that most
E. coli isolates from rodents were not of porcine origin. Their observed similarities in resistance profiles were attributed to horizontal transfer of the resistance genes via plasmids and integrons rather than transfer of the strains themselves (
14). Areas exposed to farm waste may represent “hot spots” of AMR by providing conditions that promote horizontal gene transmission between commensal, pathogenic, and environmental bacteria. Unlike most antimicrobial agents, which are easily degraded in nature, the genetic platforms containing resistance genes may be rather stable (
34). It is likely that the particular features of the Mekong Delta region, with a high density of animal farms close to waterways and low levels of biosecurity and biocontainment, as well as the common practice of discharging animal excreta directly into ponds or waterways (often used to grow fish) (
35), may indeed contribute to the dissemination and environmental persistence of resistant bacteria and resistance genes. A survey conducted in several Southeast Asian countries (including Vietnam) has documented the widespread presence of antimicrobials in the environment, probably as a result of these practices (
33).
The observed high prevalence of ampicillin resistance among isolates from trapped wild mammals is intriguing, although these results are consistent with studies elsewhere. The prevalence of ampicillin resistance among
E. coli isolates from bank voles and wood mice from forest sites in England was also very high (90.0%) (
36). Furthermore, in a Canadian study, the prevalence of ampicillin resistance was higher among wild mammals trapped in a natural environment than among those trapped on swine farms (8.0% versus 2.0%, respectively) (
31). Further molecular studies are needed to elucidate the molecular basis of this.
The relative ease of rodent/shrew trapping in the Mekong Delta suggests that routine trapping could be integrated into surveillance programs to monitor AMR in the environment, as suggested by some authors (
37).
Since AMR genes originate in environmental bacteria, there is a concern that antimicrobials used on farms and discharged into the environment may accelerate evolution toward resistance among environmental bacteria, inducing even more drastic resistance among both commensal and pathogenic bacteria (
34).
In summary, we report high levels of AMR among commensal E. coli isolates for a number of common antimicrobials, including some considered to be of critical importance, such as quinolones and aminoglycosides. Levels of resistance and MDR were highest among pig isolates. We also report an unequivocal association between AMR among isolates from small wild mammals and a farm trapping location. These results should help promote awareness of the direct links between AMR in food production and the environment. We propose that similar data should be collected over time to be used as part of integrated AMR surveillance systems to monitor the emergence of new resistance types, as well as for evaluation of interventions aimed at encouraging more sustainable farming practices, including less use of antimicrobials and better biocontainment on farms in the Mekong Delta of Vietnam.