ORIGINAL RESEARCH article
Monte Carlo Simulations Suggest Current Chlortetracycline Drug-Residue Based Withdrawal Periods Would Not Control Antimicrobial Resistance Dissemination from Feedlot to Slaughterhouse
- 1Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY, United States
- 2Department of Diagnostic Medicine/Pathobiology, Institute of Computational Comparative Medicine, College of Veterinary Medicine, Kansas State University, Manhattan, KS, United States
Antimicrobial use in beef cattle can increase antimicrobial resistance prevalence in their enteric bacteria, including potential pathogens such as Escherichia coli. These bacteria can contaminate animal products at slaughterhouses and cause food-borne illness, which can be difficult to treat if it is due to antimicrobial resistant bacteria. One potential intervention to reduce the dissemination of resistant bacteria from feedlot to consumer is to impose a withdrawal period after antimicrobial use, similar to the current withdrawal period designed to prevent drug residues in edible animal meat. We investigated tetracycline resistance in generic E. coli in the bovine large intestine during and after antimicrobial treatment by building a mathematical model of oral chlortetracycline pharmacokinetics-pharmacodynamics and E. coli population dynamics. We tracked three E. coli subpopulations (susceptible, intermediate, and resistant) during and after treatment with each of three United States chlortetracycline indications (liver abscess reduction, disease control, disease treatment). We compared the proportion of resistant E. coli before antimicrobial use to that at several time points after treatment and found a greater proportion of resistant enteric E. coli after the current withdrawal periods than prior to treatment. In order for the proportion of resistant E. coli in the median beef steer to return to the pre-treatment level, withdrawal periods of 15 days after liver abscess reduction dosing (70 mg daily), 31 days after disease control dosing (350 mg daily), and 36 days after disease treatment dosing (22 mg/kg bodyweight for 5 days) are required in this model. These antimicrobial resistance withdrawal periods would be substantially longer than the current U.S. withdrawals of 0–2 days or Canadian withdrawals of 5–10 days. One published field study found similar time periods necessary to reduce the proportion of resistant E. coli following chlortetracycline disease treatment to those suggested by this model, but additional carefully designed field studies are necessary to confirm the model results. This model is limited to biological processes within the cattle and does not include resistance selection in the feedlot environment or co-selection of chlortetracycline resistance following other antimicrobial use.
In many developed countries, the average person has very little or no contact with livestock, yet there is undoubtedly a risk of infectious agent transmission between livestock and human populations through the food supply chain and environmental intermediaries. Antimicrobial-resistant pathogens are capable of moving between animal and human populations via these pathways. Bacterial gene sequencing and multi-locus sequence typing suggest that whole bacterium transfer has occurred between at least poultry meat and humans (Overdevest et al., 2011; Kluytmans et al., 2013) and plasmid restriction-enzyme typing shows horizontal gene transfer of resistance elements between animal and human strains (Yan et al., 2004). Additionally, a recent metagenomic study on human feces found a higher prevalence of resistance genes to antimicrobials used in both animals and humans than resistance to human-limited antimicrobials (Forslund et al., 2013).
Consumers are therefore concerned about antimicrobial use in livestock and the risk of antimicrobial resistance spreading from livestock to humans by the consumption of animal products. Companies such as McDonalds, Tyson Foods, and Chipotle have responded with public commitments to limit or eliminate antimicrobials used in livestock in their supply chains (Chipotle, 2013; Tyson Foods, 2015; McDonalds, 2016). The Food and Drug Administration recently limited the use of medically-important antimicrobials in livestock for production purposes (i.e., growth promotion) and required veterinary oversight for antimicrobial use (US Food and Drug Administration, 2012, 2013). However, it is unrealistic to think that antimicrobials can be completely phased out of animal agriculture. Antimicrobials will continue to play a vital role in keeping livestock healthy and producing animal products efficiently and sustainably for a growing human population.
We must therefore devise antimicrobial stewardship policies to reduce the risk of antimicrobial resistance spread from livestock to humans. One possibility is to modify withdrawal periods required after antimicrobial use in livestock to account for resistance dynamics (Volkova et al., 2016). Antimicrobials used in livestock are already subject to meat and milk withdrawal times to prevent violative levels of drug residues from entering the human food supply. These withdrawal periods vary by country and by antimicrobial indication, varying from 0 to 10 days for oral chlortetracycline (P.R. Vademecum1; Hoffmann La-Roche Inc., 1997; Canadian Food Inspection Agency, 2014; Australian Pesticides and Veterinary Medicines Authority, 2015; Zoetis, 2016; EU Directorate General Health Consumers, 2017; New Zealand Ministry for Primary Industries, 2017; UK Veterinary Medicines Directorate, 2017). “Resistance withdrawals” could be used to decrease the prevalence of resistant bacteria and genes in an animal's enteric microbiome before it is sent for slaughter, where enteric bacteria can contaminate meat products (De Filippis et al., 2013).
In order to investigate the efficacy and length of resistance withdrawal periods that would reduce the prevalence of resistance in the livestock enteric microbiota, we built a model of the pharmacokinetics-pharmacodynamics of chlortetracycline (CTC) fed to beef steers, and the population dynamics of their enteric generic Escherichia coli. Chlortetracycline is a commonly used antimicrobial on beef feedlots (United States Department of Agriculture, 2013), and tetracyclines were identified as a highly important antimicrobial class for human medicine by the FDA (US Food Drug Administration, 2003), making CTC a suitable candidate for evaluating the impact of antimicrobial resistance withdrawal periods.
The model consists of 3 sub-models that were connected and parameterized for CTC and generic E. coli in beef cattle: a pharmacokinetic model for the concentrations of orally administered drugs in the gastrointestinal tract (Cazer et al., 2014), a bacterial population dynamics model, and a pharmacodynamic model (Ayscue et al., 2009; Volkova et al., 2012). The framework for connecting these models and using them to investigate resistance withdrawal periods has been described (Volkova et al., 2016). Each model is described briefly below. Model equations are listed in Table 1. Each parameter was considered to be either a constant or a random variable. Each simulation of the three connected models represents one realization of a treated animal, with the values of the random variables drawn from the animal population-level parameter distributions. The model was parameterized using data available from literature wherever possible. For a parameter with only two published estimates, a Uniform distribution was used with the estimates as the maximum and minimum. For a parameter with only one published estimate, a Uniform distribution with the maximum and minimum being ±25% of the estimate was assigned (Volkova et al., 2017). For a parameter with at least 3 published estimates, the parameter distribution was determined by fitting candidate distributions to the dataset of the published estimates of the parameter values, with PROC CAPABILITY in SAS® 9.4 software for Windows (SAS Institute Inc., Cary, NC, U.S.); the best-fit distribution for the parameter was selected using goodness-of-fit tests (Anderson-Darling and Kolmogorov-Smirnov) and visual examination of quantile-quantile plots. The parameter distributions used in the model simulations are listed in Table 2. (The data used to fit the parameter value distributions for the parameters with at least 3 published estimates are included in Supplementary Materials).
Equations of the pharmacokinetic model are given in Table 1A. The daily dose of CTC was fed to a 300 kg steer in equal parts throughout a 12-h day, during each of the treatment days (Table 1A, Equations 1–4). The drug moved through the steer's gastrointestinal tract (Table 1A, Equations 5–7) at the rate of concentrates, determined from studies using isotope-labeled feedstuffs (Shaver et al., 1986; Zebeli et al., 2007). In the small intestine, some CTC was absorbed into plasma and from there distributed into tissues (Bradley et al., 1982), excreted into urine and via bile back into the small intestine (Eisner and Wulf, 1963) (Table 1A, Equations 6, 8, 9). In all gastrointestinal, plasma, and tissue compartments, CTC degraded abiotically (Eisner and Wulf, 1963) into compounds with greatly reduced antimicrobial activity (Halling-Sørensen et al., 2002).
Escherichia coli Population Model
The model of E. coli growth considered only free-living, luminal E. coli in the large intestine because the amount of colon mucosal-associated E. coli is 1-2 magnitudes less than the luminal bacteria (Krause et al., 2003; Laven et al., 2003; Volkova et al., 2012). The previous model of an E. coli population divided the bacteria into two subpopulations, resistant and susceptible (Volkova et al., 2012). We expanded the model to include E. coli intermediately susceptible to tetracyclines (Clinical and Laboratory Standards Institute, 2015). We assumed that tetracycline resistance genes were carried on conjugative elements (transposons or plasmids) (Speer et al., 1992). Each subpopulation of bacteria can expand or contract due to logistic growth until the total population reaches a carrying capacity (Table 1B, Equation 11). The bacterial subpopulations are also affected by fitness costs of resistance elements, inflow and outflow of E. coli from the large intestine (Table 1B, Equations 12, 13), and horizontal plasmid/transposon transfer among resistant, intermediate, and susceptible strains (Table 1B, Equation 14) (Volkova et al., 2012).
Tetracyclines are bacteriostatic antimicrobials at physiologically achievable concentrations that inhibit protein synthesis and suppress the growth of Enterobacteriaceae (Norcia et al., 1999). The pharmacodynamic model was based on a sigmoid Emax model in which the pharmacodynamic effect is a reduction of the bacterial population growth rate (Table 1B, Equation 15) (Mouton and Vinks, 2005; Goutelle et al., 2008; Volkova et al., 2012; Ahmad et al., 2015b; Wen et al., 2016).
It has been shown that susceptible isolates of Enterobacteriaceae exhibit some growth suppression when exposed to tetracycline concentrations below their MIC, whereas resistant isolates are more robust and not significantly growth-suppressed until the concentration of tetracycline equals or exceeds their MIC (Gullberg et al., 2011; Ahmad et al., 2015b). Additionally, the drug concentration-effect relationship (captured by the pharmacodynamic parameter Hill coefficient) changes depending on the isolate MIC in E. coli and other Gram-negative pathogens (Ahmad et al., 2015b; Wen et al., 2016). Therefore, different Hill coefficients were applied for the susceptible, intermediate, and resistant E. coli subpopulations. For tetracycline in E. coli, EC50 is a function of MIC (Ahmad et al., 2015a,b). MIC is routinely measured under aerobic laboratory conditions but the lumen of the large intestine is anaerobic. This change in environment has been shown to affect the MIC of E. coli (DeMars et al., 2016) so we applied an anaerobic penalty to the standard aerobic MIC to reflect this difference.
In the United States, the chlortetracycline label indications for beef cattle are the control of active Anaplasma marginale infection and shipping fever (bacterial pneumonia) caused by Pasteurella species (350 mg per head daily; “ADC”), reduction of incidence of liver abscesses (70 mg per head daily; “ARLA”), and treatment of Escherichia coli enteritis or Pasteurella pneumonia (22 mg per kg bodyweight per day for 5 days [6,600 mg per head daily for a 300 kg steer]; “ADT”) (Zoetis Aureomycin 502; ChlorMax 503; Zoetis, 2016). Chlortetracycline was previously labeled for growth promotion and increased feed efficiency but those label claims were removed as of January 2017 in response to the FDA's Guidance for Industry #213 (US Food and Drug Administration, 2013). However, the dosage that was used for growth promotion (70 mg per head daily) is still labeled for the reduction of liver abscesses. In the model, CTC was fed at ARLA or ADC dosage for 28 days or at ADT dosage for 5 days. The treatment period for ADT was considered to last an additional 3 days because CTC concentrations in the large intestine remained increased for 3 days after ADT ended.
Some brand-name CTC products sold in the U.S. have no required withdrawal period for any of the dosages (Zoetis Aureomycin 502). Generic CTC products sold in the U.S. have a 0-day withdrawal following ARLA, 2-day withdrawal following ADC, and 1-day withdrawal following ADT (Zoetis ChlorMax 50). In the model, the 0-day withdrawal was interpreted to mean that cattle could ingest CTC until they were put on a truck for transportation to slaughter and we assumed that they were slaughtered 6 h later. The travel time to the slaughterhouse was assumed to be included in the longer withdrawal periods (1 or 2 days).
The pharmacokinetic parameters were the same as in previous models for CTC (Cazer et al., 2014; Volkova et al., 2017) and are listed in Table 2A. In summary, sufficient data for fitting a distribution to the parameter values from literature were available for only the CTC degradation rate. The best-fit model for the unweighted, mean CTC degradation rate at physiological pH and temperature (Eisner and Wulf, 1963; Carlson and Mabury, 2006; Arikan, 2008; Dolliver et al., 2008; Arikan et al., 2009) was a Beta distribution (Anderson-Darling P > 0.250, Chi-square P = 0.052) (Volkova et al., 2017). Only one estimate for the fraction of CTC excreted in bile was available (Eisner and Wulf, 1963); hence, the Uniform distribution of the parameter values was assigned the boundaries ±25% of the estimate. In the absence of cattle-specific data, the fraction of CTC adsorbed to digesta was parameterized with data from the bioavailability of tetracycline in rat feces (Bahl et al., 2004) and protein-bound CTC in dog serum (Pindell et al., 1959); these estimates were used as the boundaries for the Uniform distribution of the parameter values. Finally, a Uniform distribution for the volume of the contents in the large intestine was estimated using the data on the content weights at slaughter (Murray et al., 1977) and the density of fresh cattle feces (Volkova et al., 2017) (Table 2A). The parameters related to the distribution of CTC in central circulation and in tissues were kept constant in order to focus on the variability of CTC concentration in the large intestine. For this modeling exercise, we assumed that the cattle ate a grain-based diet and the CTC flow rates through gastrointestinal compartments were parameterized accordingly (Shaver et al., 1986; Zebeli et al., 2007). Only single estimates were available for the transit rates of concentrates through the small intestine and large intestine (Shaver et al., 1986); hence, each of these transit rates was assigned a Uniform distribution with the boundaries ±25% of the published estimate. Two estimates for the transit rate from the stomachs to small intestine were located (Shaver et al., 1986; Zebeli et al., 2007) but there was only a 0.0005 h−1 difference between the estimates; hence, a Uniform distribution with the boundaries ±25% of the average of the two estimates was assigned to the stomachs' transit rate. The steer's body weight was kept constant since growth over 40 days was previously found to have only minor impacts on CTC concentration in the large intestine in this model (Cazer et al., 2014).
Escherichia coli Population Model
Parameter distributions for the E. coli population dynamics parameters and parameters of CTC pharmacodynamics against E. coli are given in Table 2B. The net growth rate of E. coli, accounting for the anaerobic conditions in the cattle large intestine, was adopted from a previously published model (Volkova et al., 2012), in which the growth rate was estimated based on the data of laboratory experiments of E. coli growth in animal cecal contents (Freter et al., 1983b). Other literature (Durso et al., 2004) supports that the population growth rate of generic E. coli in anaerobic conditions is reduced compared to aerobic growth. We were able to identify only two estimates for the fitness cost of tetracycline resistance in E. coli (Nguyen et al., 1989; Ahmad et al., 2015b). In the first study, E. coli demonstrating phenotypic tetracycline resistance grew at a rate equal to susceptible E. coli in the absence of antimicrobial exposure, indicating no fitness cost of tetracycline resistance (Ahmad et al., 2015b). However, the second study found a fitness cost to E. coli from carrying plasmids, excluding the cost of tetracycline genes (1–2% decreased growth rate). Constitutive tetracycline-resistant mutants carry a relatively high fitness cost (2–4%), although inducible tetracycline resistance genes only confer a fitness cost of at most 0.3% (Nguyen et al., 1989). The data for E. coli agree with data for Salmonella; tetracycline resistance plasmids confer a small fitness cost to Salmonella (Gullberg et al., 2011). We assumed that susceptible E. coli in our model did not carry conjugative resistance elements and that intermediate and resistant E. coli carried plasmids or transposons with inducible tetracycline genes and therefore may experience a fitness cost in the absence of tetracycline. In the absence of additional evidence, we assumed that intermediate and resistant E. coli experienced the same fitness cost, which was estimated to have a Uniform distribution from 0 (Ahmad et al., 2015b) to a 3% reduction in the bacterial hourly growth rate (Nguyen et al., 1989; Gullberg et al., 2011).
In the implemented model, E. coli in the enteric population can grow until they reach the carrying capacity; thus, the population growth plus inflowing bacteria equals death plus outflowing bacteria. Estimates of the E. coli carrying capacity in the large intestine of feedlot cattle (log10CFU/g feces) were obtained from nine studies (n = 44 reported averages from 3,610 fecal or colon samples) (Callaway et al., 2002; Krause et al., 2003; Laven et al., 2003; Aslam et al., 2004; Branham, 2007; Lowrance et al., 2007; Alexander et al., 2009, 2010; Agga et al., 2016). Candidate distributions were fit to this data and a Weibull distribution was found to have the best fit (Kolmogorov-Smirnov P = 0.125, Anderson-Darling P = 0.041), and thus was used to model the inter-individual variability in the large intestine E. coli carrying capacity. The total E. coli population at time 0 was allowed to vary from 10 to 90% of the carrying capacity.
Plasmids or transposons carrying antimicrobial resistance genes can transfer between susceptible, intermediate, and resistant bacteria, affecting the fraction of E. coli within each subpopulation in the large intestine. The rate of conjugative transfer of tetracycline resistance among free-living E. coli in the bovine large intestine is not known. In vitro density-dependent plasmid transfer rates (mL/cell*hour) for E. coli and the plasmid R1 were extracted from four studies (n = 138 transfer experiments) (Levin et al., 1979; Freter et al., 1983b; Simonsen et al., 1990; Gordon, 1992) and converted to frequency-dependent rates (hour−1) by multiplying the density-dependent rate by the initial cell density. The best-fit distribution for the frequency-dependent plasmid transfer rate was a Gamma distribution (Kolmogorov-Smirnov P > 0.50, Anderson-Darling P = 0.176). In the absence of conjugative element transfer data specifically for tetracycline resistance in E. coli, we assumed that the rate of gene transfer from resistant or intermediate bacteria to recipients (susceptible or intermediate bacteria) would come from the plasmid-transfer population parameter distribution described above.
Ingestion of E. coli from the pen environment, feed, and water, as well as transit of resident small intestine populations, and shedding of mucosal-associated E. coli into the lumen provides an inflow of new bacteria into the large intestine digesta. The defecation of feces containing E. coli constitutes the outflow of E. coli from the large intestine (Ayscue et al., 2009). We kept the same range of inflow and outflow rates as our previous E. coli population model (Volkova et al., 2012) and divided the range between inflow and outflow via an iterative process in order to keep E. coli populations at or below carrying capacity and prevent the accumulated losses of E. coli (outflow and death) from being greater than accumulated gains (inflow and growth) over the time period of the model. This maintained a stable total E. coli population in the model.
We assumed that the majority of ingested E. coli came from fecal contamination of the pen environment, feed, and water. We also assumed that the resident small intestinal and mucosal-associated E. coli would have the same resistance levels as colon luminal E. coli. Therefore, the proportions of inflow bacteria that are resistant or intermediate were estimated from studies of E. coli antimicrobial resistance in feedlot cattle feces. These studies cover a wide range of management and sampling strategies so the 10th and 90th percentiles of tetracycline resistance prevalence across these studies were used as the minimum and maximum in a Uniform distribution (Wagner et al., 2002; Branham, 2007; Lowrance et al., 2007; Alexander et al., 2008; Carson et al., 2008; Platt et al., 2008; Rao et al., 2010; Benedict, 2011; Morley et al., 2011; Kanwar et al., 2013; McGowan, 2014) of the resistance in inflow bacteria. Non-overlapping Uniform distributions were used in order to keep the sum of the incoming resistant, intermediate and susceptible E. coli proportions at 1. The starting proportions of resistant, intermediate, and susceptible E. coli in the large intestine were drawn from the same Uniform distributions.
Data on the tetracycline EC50 of E. coli isolates from cattle is not available, so data for 50 E. coli swine isolates were used to determine the relationship between EC50 and MIC (Ahmad et al., 2015a). A linear regression model (deemed appropriate based on the data for swine E. coli; Ahmad et al., 2015b) for the value of log2(EC50) depending on the isolate log2(MIC) was fit in PROC REG in SAS. The resulting model (R2 = 0.99) is given in Equation 16 (Table 1B) and was used to predict the EC50 of cattle E. coli isolates from their MIC. The breakpoints (established by the Clinical and Laboratory Standards Institute) for E. coli interpretation in a veterinary infection as resistant, intermediate, and susceptible to tetracycline were used to define the MIC of each of the E. coli subpopulations; Uniform MIC distributions were used in the model. The upper limit of resistant MIC considered was 128 μg/mL. The anaerobic MIC penalty was assumed to range between the 1st percentile (−1.3) and 99th percentile (0) reported difference between anaerobic and aerobic tetracycline MIC for generic E. coli (DeMars et al., 2016). Therefore the lower-bound of the aerobic MIC distributions was decreased by 1.3, except for the susceptible MIC which could not be less than 0.
The Hill coefficients used in the pharmacodynamic model were based on the estimates for the same set of 50 E. coli swine isolates (Ahmad et al., 2015b). A best-fit model for the Hill coefficients could not be found so the interquartile ranges were used as the boundaries of the Uniform distributions of the Hill coefficients for resistant and susceptible E. coli. Only one intermediate isolate was examined (Ahmad et al., 2015b) so a ±25% interval around the Hill coefficient value for this isolate was used as the distribution boundaries for the intermediate E. coli Hill coefficient.
The differential equations of the model were implemented in MatLab® R2016b (MathWorks, Natick, MA, U.S.) using a time-step of 0.1 h. The Escherichia coli population was allowed to reach equilibrium of the resistant, intermediate, and susceptible proportions before CTC was fed. The model was run for a 90 day simulation period, including the initial 48 h required to reach the pre-treatment equilibrium (Day 0–2), the treatment period (Day 2–10 for ADT, Day 2–30 for ARLA and ADC), and post-treatment period (Day 11–90 for ADT, Day 31–90 for ARLA and ADC). A thousand simulations were run for each treatment scenario, representing 1,000 treated cattle; simulations were also run for a 1,000 cattle with no CTC treatment. Descriptive statistics of the E. coli subpopulations were calculated in MatLab using the IOSR statistics toolbox. The bacterial subpopulation proportion medians were tested for equality with Kruskal-Wallis tests followed by Dunn-Sidak post-hoc comparisons. All significance tests were two-sided with α = 0.05. The following time-periods were defined for further investigation: pre-treatment (hour 42–48), during treatment (the steady-state from hour 552 to 720 for ARLA and ADC; hour 48 to 240 for ADT), day 90 (last 24 h), 0-day withdrawal period (6 h after treatment ends), 1-day withdrawal period (24 h after treatment ends), and 2-day withdrawal period (48 h after treatment ends). Figures were created in MatLab using the IOSR statistics toolbox and SubPlot toolbox. The MatLab code for the model is provided in Supplementary Materials.
All three scenarios of CTC indications for beef cattle (antimicrobial liver abscess reduction—ARLA; disease control—ADC; disease treatment—ADT) were implemented in the model of CTC pharmacokinetics-pharmacodynamics and generic E. coli population dynamics in the bovine large intestine. The average CTC concentration in the large intestine during treatment was calculated for the median steer, 5th percentile steer, and 95th percentile steer. For ARLA, the median steer had an average large intestine CTC concentration of 0.22 μg/mL, the 5th percentile steer had a concentration of 0.06 μg/mL, and the 95th percentile steer had a concentration of 0.62 μg/mL. For steers given ADC, the averages were 1.17, 0.34, and 2.96 μg/mL for median, 5th percentile and 95th percentile, and for ADT the average CTC concentrations were 12.99, 3.70, and 34.12 μg/mL, respectively. This is consistent with the results from our previous deterministic pharmacokinetic model (Cazer et al., 2014).
The use of CTC for any indication resulted in an increase in the median proportion (Figures 1D,G,J) and absolute number (Figure S1) of resistant E. coli during treatment in 1,000 simulated cattle. The median proportion resistant was 38% before treatment and increased to 46% at the end of ARLA, 60% at the end of ADC, and 66% at the end of ADT. In terms of the number of resistant E. coli in the large intestine, there were 105.76 CFU/g before treatment, 105.82 CFU/g at the end of ARLA, 105.93 CFU/g at the end of ADC, and 105.89 CFU/g at the end of ADT. The median proportion of intermediate E. coli also increased during ARLA (10% at the end of treatment) and ADC (11%) but decreased during ADT (5%) relative to pre-treatment levels (8%) (Figures 1B,E,H,K). However ADT still greatly increased the median proportion (71% at the end of ADT) and amount (105.92 CFU/g) (Figure S1) of non-susceptible bacteria compared to before treatment (46%; 105.84 CFU/g), accompanying the drop in the median proportion of susceptible E. coli (54% before treatment and 29% at the end of ADT) (Figure 1L). The proportions and amounts of resistant, intermediate, and susceptible bacteria reached an approximate steady state during ARLA and ADC (Figures 1D,G; Figure S1), whereas a steady state was not reached during the shorter ADT and the proportion of resistance peaked, on average, on the second day after ADT treatment was discontinued (Figure 1J). In general, ARLA and ADC resulted in an upward shift of the median and upper limit (95%) of the proportion non-susceptible distribution but little change in the lower limit (5%) during antimicrobial treatment (Figures 1F,I). In contrast, during the ADT period and following 3 days the entire proportion resistant distribution shifted upwards (Figure 1J).
Figure 1. Simulation of Escherichia coli subpopulations (resistant [A, D, G, J]; intermediate [B, E, H, K]; susceptible [C, F, I, L]) in the large intestine of beef cattle in the presence (D–L) and absence (A–C) of oral chlortetracycline treatment. Proportions (Y-axis) are presented for 1,000 simulations of each treatment scenario. Shaded band is the 95% non-parametric confidence interval of the median, black dashed lines are 25 and 75% percentiles and purple dashed lines are 5 and 95% percentiles of the overall distribution. The red shaded band represents the median proportion in the absence of chlortetracycline (CTC) treatment. Blue shaded bands represent the median proportion in antimicrobial treatment scenarios for reduction of liver abscesses (ARLA; D–F), disease control (ADC; G–I) or disease treatment (ADT; J–L). The green and red vertical lines mark the beginning and end of CTC treatment, respectively.
The median proportion of resistant E. coli reversed to the pre-treatment median confidence interval 15 days after ARLA ended, 31 days after ADC ended, and 36 days after ADT ended (Figures 1B–D). Only 35% of ARLA cattle, 50% of ADC cattle and 61% of ADT cattle returned to within 10% of their individual pre-treatment resistant E. coli proportion by the end of the 90 day simulation period. This includes 18% of ARLA cattle, 9% of ADC cattle, and 6% of ADT cattle that were already within 10% of pre-treatment levels at the end of antimicrobial treatment. Among the individual cattle that did return to pre-treatment levels, the median time from the end of treatment to the resistance reversion was significantly different between the treatment groups (P < 0.001): <1 day (mean 6 days) following ARLA, 8 days (mean 11 days) following ADC, and 11 days (mean 15 days) following ADT. Since ADT duration is shorter, the 90 day simulation period captured 60 days post-ARLA and -ADC but it captured 83 days post-ADT. At 60 days post-treatment, 59% of ADT cattle had returned to within 10% of pre-treatment proportion resistant with a median return time of 11 days (mean 14 days) after the treatment.
In the absence of treatment, 38% of generic E. coli in the cattle large intestine were resistant, 8% intermediate, and 54% susceptible to CTC on average (Figure 2). There was no significant difference in the average proportions of resistant bacteria during ADC (58%) and ADT (59%) (P = 0.850, Figure 2A), but ADC had a significantly higher proportion of intermediate bacteria (12%) compared to ADT (6%) (P < 0.001) (Figure 2B). This was associated with ADC having a lower average susceptible proportion of E. coli (30%) than ADT (35%, Figure 2C) despite ADT having the highest maximum proportion resistant (72%). At the steady state during ARLA, 46% of E. coli were resistant, 10% intermediate, and 44% susceptible on average. The average proportion of resistant E. coli during any of the three CTC treatments was significantly different from that in the absence of treatment (38%) (P < 0.001, Figure 2A). By the end of the 90 day simulation period, the mean proportion resistant reversed and was not significantly different from that in the absence of treatment (P > 0.050, Figure 2A).
Figure 2. Distribution of proportions of tetracycline resistant (A), intermediate (B), and susceptible (C) Escherichia coli in the large intestine of beef cattle at four time periods. Day 2 is the mean proportion during the 6 h before chlortetracycline treatment starts. Treatment is the mean proportion from Day 2 to 10 for ADT and from Day 23 to 30 for ARLA, ADC, and in the absence of treatment. Max is the maximum proportion between Days 2 and 12 for ADT and between Days 2 and 35 for ARLA, ADC, and in the absence of treatment. Day 90 is the mean proportion during the 24 h of the last day of the simulation period. Shaded box extends from 25th to 75th percentile of the simulated proportion. Middle line is the median and plus-symbol is the mean. Whiskers extend to the minimum or maximum data point within 1.5 times the interquartile range from the first or third quartile, respectively. The results of 1,000 simulations of each treatment scenario are summarized in each boxplot.
Brand-name U.S. CTC feed products have a 0-day meat withdrawal period for beef cattle so they can be sent to slaughter immediately after the consumption of CTC. We assumed that 6 h elapse between leaving the feedlot and slaughter. The median cattle given ARLA had 46% resistant enteric E. coli at slaughter; the median ADC cattle had 60% resistant E. coli, and the median ADT cattle had 67% resistant E. coli (Figure 3). Generic U.S. CTC products have a 2-day withdrawal period following ADC and a 1-day withdrawal period following ADT. These withdrawals resulted in statistically significant different proportions of resistant E. coli at slaughter compared to the 0-day withdrawal: the median cattle had 58% resistant E. coli at 2 days after ADC (P = 0.012) and 71% resistant 1 day after ADT (P < 0.001) (Figure 3). The cattle treated with CTC had a wider distribution and greater maximum proportion resistant at the currently recommended withdrawal times compared to those in the absence of treatment (Figure 3). The lower bound of the resistance distribution at the withdrawal times was higher for ADT relative to the other two treatments (Figure 3).
Figure 3. Distribution of proportion tetracycline resistant Escherichia coli in the large intestine of beef cattle at the current withdrawal times. For “0 Day” withdrawal period (e.g., “No Withdrawal”) the proportion resistant averged over the 6 h after treatment ended is given. For “48 h Withdrawal” and “24 h Withdrawal” periods, the proportion resistant averaged over 48 and 24 h, respectively, after treatment ended is given. Shaded box extends from 25th to 75th percentile of the simulated proportion resistant. Middle line is the median and plus-symbol is the mean. The notch is the 95% non-parametric confidence interval of the median. Whiskers extend to the minimum or maximum data point within 1.5 times the interquartile range from the first or third quartile, respectively. The results of 1,000 simulations of each treatment scenario are summarized in each boxplot.
Before and after CTC treatments, the proportions of CTC resistant, intermediate, and susceptible enteric E. coli in the large intestine were sensitive to variability in the inflowing E. coli distribution of resistance (pj). During the treatments, the proportions were sensitive to variability in the bacterial population dynamics parameters and in the CTC pharmacokinetic-pharmacodynamic parameters (Figure 4, Spearman coefficients smaller than the Bonferroni corrected α = 0.0006 are shown). In general, higher CTC dosages made the proportions less sensitive to the distribution of resistance in the inflowing E. coli but more sensitive to the overall inflow rate (λin) (Figure 4). In addition to the inflow rate and the inflow resistance distribution, the proportions during ADT were sensitive to the starting proportions (startj) and the CTC minimum inhibitory concentration (MIC) for the resistant E. coli. Conversely, the proportions during ARLA and ADC were sensitive to the MIC for the susceptible E. coli. Variability in the following pharmacokinetic parameters had a significant impact on the E. coli proportions only during ARLA and ADC (Figures 4A,C): CTC abiotic degradation throughout the gastrointestinal tract and other organs (δ), CTC binding to the digesta in the large intestine (η), and volume of the digesta in the large intestine (Vli). Variability in the Hill coefficient (Hj), a pharmacodynamic parameter, had no significant impact on the E. coli proportions during treatments. Variability in the bacterial population dynamics parameters of the growth rate (r), resistance fitness cost (α), carrying capacity (Nmax), and the resistance horizontal transfer rate (β), plus the pharmacokinetic parameters of the CTC flow rates through the gastrointestinal compartments (γs, γupper_si, γrest_si, γli) and the fraction excreted in bile (Eb), had no significant effect on the proportions of resistant, intermediate, and susceptible E. coli during treatments.
Figure 4. Sensitivity of proportions of tetracycline resistant (A), intermediate (B), and susceptible (C) Escherichia coli in the large intestine of beef cattle during chlortetracycline treatment to the parameter values of the drug pharmacokinetics, pharmacodynamics, and the bacterial population dynamics. Only parameters that significantly correlate with the proportion (Spearman correlation coefficient P < 0.0006, Bonferroni correction) are included. Parameters are defined in Table 2. The time period “during chlortetracycline treatment” was from Day 2 to 10 for ADT and from Day 23 to 30 for ARLA, ADC, and in the absence of treatment. The correlations between the mean proportions during those time periods and parameter values were evaluated using the outputs of 1,000 simulations of each treatment scenario.
Despite modeling a long follow-up period after ending CTC treatment (60 days after ARLA and ADC, 83 days after ADT), a significant number of the simulated cattle had still not returned to their pre-treatment proportion of resistant E. coli after CTC treatment. The dosage of CTC appeared to be associated with the percent of cattle that returned to their pre-treatment levels, with higher dosages (ADT, ADC) having greater proportions of cattle returning to pre-treatment resistance than low-dosage ARLA. This suggests that the low-dosage long-term ARLA may more frequently shift the enteric E. coli equilibrium to a state that is stable in the absence of CTC. This is evident in the change in proportion resistant from the end of treatment to the end of the simulation period (day 90). ARLA cattle that never returned to pre-treatment levels (65% of ARLA cattle) had a median difference in proportion resistant of 1.7% from the end of ARLA to day 90, indicating that they had reached an equilibrium during treatment that was relatively stable in the absence of CTC. On the other hand, ARLA cattle that did return to pre-treatment levels (35% of ARLA cattle) had a larger median change in proportion resistant (6.4%) from the end of treatment to the end of the simulation period. It is possible that the high levels of resistance observed during ADT and ADC were unsustainable after treatment given the model parameters (particularly the resistant inflow proportion). Only 2.4% of ADC cattle and 4.9% of ADT cattle achieved an equilibrium during treatment that was relatively stable in the absence of CTC (<1.7% difference between end of treatment and day 90 proportion resistant). However, the ARLA cattle that returned to their pre-treatment levels did so much faster than ADC and ADT cattle, likely because there was a smaller change in resistance level from pre-treatment to the end of ARLA (Figure 1).
Additional investigation is needed to study the impact of treatment duration separately from dosage, including testing off-label ADT treatment durations. Pulse-therapy with ADT dosages (3 treatments, 1 day apart) has been implemented in a beef cattle field study; an approximate steady state of 80% resistant E. coli isolates was reached by the second treatment (double the pre-treatment proportion of 40%) (Platt et al., 2008). This is similar to the mean pre-treatment proportion resistant (38%) and the mean maximum proportion resistant after ADT (72%) in our study. At the end of the 17 day post-treatment follow-up period, the proportion of resistant isolates in that study was still 1.5 times greater than the pre-treatment proportion. If the decrease continued at the same rate, the proportion resistant would reverse to the pre-treatment level 34 days after the treatment ended, similar to the 36 days required in our model for the median proportion resistant to be within 10% of the pre-treatment median. Additionally, the average proportion resistant during ADT in our study (59%) fell within the 95% confidence interval of the average proportion resistant during ADT pulse-therapy (58.7–77.3%) (Platt et al., 2008), confirming that our model is projecting biologically plausible resistance dynamics.
Unfortunately it is difficult to compare our model results to the data from field trials on E. coli tetracycline resistance prevalence following CTC treatment due to the variable field study design and management conditions. Those include pooling of isolates across treatment groups and time (Alexander et al., 2008), combination antimicrobial therapy (Kanwar et al., 2013), and reporting of animal-level rather than isolate-level results (Sharma et al., 2008). In order to fully validate the resistance prevalence predictions of this model, field studies are required that sample individual animals, determine the prevalence of tetracycline resistant E. coli isolates before and at several points after CTC treatment, and do not alter diet or environmental conditions during the study period.
This model illustrated important differences in individual and population level antimicrobial resistance. Each simulated population of 1,000 cattle reached its pre-treatment population distribution of tetracycline resistance within the follow-up period (Figure 1). This return occurred approximately 2 weeks after ARLA and 4–5 weeks after ADC or ADT ended. Even though the resistance levels of many individual animals were still significantly altered by CTC treatment, the population or herd resistance was able to return to its initial distribution. This implies that there would be no net increase in resistance exported from the feedlot to slaughterhouses if a sufficient antimicrobial resistance withdrawal period was applied.
The resistance dynamics following different dosages should be considered in antimicrobial use policies. When Denmark stopped using growth promoting doses of antimicrobials (like ARLA) in weaner pigs, the amount of antimicrobials used for disease treatment (like ADT) increased due to higher disease incidence (Kjeldsen and Callesen, 2006). A similar change in the use of CTC may result in significant increases in resistance exported from farms to slaughterhouses, depending on disease incidence and CTC withdrawal period before slaughter. However, long-term ARLA or ADC may release significantly more resistant E. coli, cumulatively, into the feedlot environment compared to one 5-day course of ADT. This highlights the need for additional rigorous, longitudinal trials comparing CTC labeled uses and collecting data on E. coli in cattle, in the feedlot environment, and on beef products.
Based on the model outcomes, it is evident that the current U.S. withdrawal periods (designed to mitigate violative drug residues in the edible meat) for in-feed CTC (Zoetis Aureomycin 50; ChlorMax 50; Hoffmann La-Roche Inc., 1997) in beef cattle are insufficient to reduce the treatment-associated prevalence of resistant E. coli in the large intestine, which can contaminate processing plants and beef products at slaughter (Sheridan, 1998; Aslam et al., 2003). Currently, U.S. cattle can consume CTC until immediately prior to slaughter and hence the proportion of tetracycline resistant bacteria in their intestine is significantly greater than in untreated cattle (Figure 3). The 24-h withdrawal after ADT may further increase the proportion of resistant E. coli at slaughter because of the CTC intestinal transit times.
Apart from the United States, the only other country with a large beef cattle industry that permits the use of oral CTC is Canada. They allow CTC fed at 70 mg/head per day (ARLA dosage) to prevent foot rot in beef cattle with a pre-slaughter withdrawal period of 5 days (Canadian Food Inspection Agency, 2016). Canada also has an approved CTC-sulfamethazine combination product that is fed at 350 mg CTC per head per day (ADC dosage) with a 10 day pre-slaughter withdrawal period (Canadian Food Inspection Agency, 2016). The European Union, United Kingdom, Australia, New Zealand, Brazil, Argentina, and Mexico have no approved oral CTC products for beef cattle over 6 months old, although a few permit topical CTC sprays with a 0 day pre-slaughter withdrawal period and New Zealand has an approved intrauterine CTC product with a 7 day pre-slaughter withdrawal period (P.R. Vademecum1; Australian Pesticides and Veterinary Medicines Authority, 2015; EU Directorate General Health Consumers, 2017; New Zealand Ministry for Primary Industries, 2017; UK Veterinary Medicines Directorate, 2017). This modeling study demonstrated that a 15 day withdrawal after ARLA, a 31 day withdrawal after ADC, and a 36 day withdrawal after ADT (Figure 1) may be necessary to reduce the median prevalence of tetracycline resistance among enteric E. coli that could contaminate beef products. Such withdrawal periods would be longer than any current CTC withdrawal periods (which are currently used to control only the CTC residues in the animal edible tissues) in countries with large beef industries. Similar or longer withdrawal periods have been used in antimicrobial feeding studies without detriment to cattle health and welfare (Alexander et al., 2008; Sharma et al., 2008; Agga et al., 2016). However, the economic and welfare costs of resistance withdrawal periods will depend on the specific disease treated and management strategy.
The same field studies (Alexander et al., 2008; Agga et al., 2016) did not find the current withdrawal periods effective in decreasing tetracycline resistance levels. However, one study (Alexander et al., 2008) observed a large increase in tetracycline resistant fecal E. coli when switching from a silage to grain based diet and this effect overwhelmed the ability to detect any changes due to antimicrobial withdrawal. The second study (Agga et al., 2016) found a decrease in fecal E. coli tetracycline resistance during a 27-day withdrawal after ADT but then observed a steady increase in both control and treated groups for an additional 90 days. They concluded that resistance is selected for in cattle-occupied pens independently of antimicrobial use (Agga et al., 2016). Since our model includes only effects of the administered antimicrobials, it is unable to simulate the pattern they observed. The model could be expanded to include the effects of diet (Alexander et al., 2008; Volkova et al., 2016, 2017), environmental E. coli populations (Ayscue et al., 2009; Volkova et al., 2013), environmental contamination with antimicrobial drugs (Call et al., 2013), and animal feeding or environmental contamination with other antimicrobial compounds such as heavy metals (Berendonk et al., 2015).
The sensitivity of the outputs of our model to E. coli population dynamics parameters, particularly the inflow and outflow of E. coli in the large intestine, emphasizes the lack of knowledge about E. coli and antimicrobial resistance circulation in feedlot environments. It also suggests that interventions aimed at changing the turnover of E. coli in the large intestine and the incoming proportions of resistant bacteria may be highly effective at reducing antimicrobial resistance. For example, moving cattle to a previously unoccupied pen for the resistance withdrawal period may reduce the ingestion of resistant bacteria (Agga et al., 2016) and shift the enteric population to a lower-resistance equilibrium. Additional research is needed on the circulation of antimicrobial resistant bacteria and genes in the feedlot environment and how the environmental bacterial populations intersect with the enteric ones. Probiotics and some feed additives have been shown to alter E. coli O157 dynamics in beef cattle (Sargeant et al., 2007) and similar strategies could complement an antimicrobial withdrawal period. Altering pharmacokinetic parameters, such as CTC degradation rate in the gastrointestinal tract or its sorption to digesta, may also be an avenue for decreasing the antimicrobial exposure of enteric microbiomes.
The level of resistant bacteria during CTC treatment was not significantly associated with the fitness cost value or the plasmid transfer rate value in the global sensitivity analysis. These findings are consistent with mathematical models of tetracycline resistant E. coli in the swine intestine, which found that the effect of the plasmid transfer rate was inconsequential compared to the effect of the growth rate of resistant bacteria (Graesboll et al., 2014; Ahmad et al., 2015a). The same models also found that resistant bacteria could persist in the swine intestine with a fitness cost up to an 80% reduction in the growth rate (Graesboll et al., 2014), although this was much greater than the fitness cost they observed in vitro (Ahmad et al., 2015b). Consistent with our model of tetracycline resistance in E. coli in the cattle intestine, the models of tetracycline resistant E. coli in the swine intestine found that bacterial inflow and outflow had a stronger correlation with the level of resistance during treatment than fitness cost or plasmid transfer rates (Graesboll et al., 2014).
By necessity, mathematical models are a simplified representation of biologic systems and can omit confounding or integral parameters. Since our model was limited to modeling E. coli within the large intestine, we were unable to directly account for environmental factors that may impact resistance dissemination and persistence in beef feedlots. In addition, this model did not account for co-selection of CTC resistance as a result of administering other antimicrobials, which may be an important mechanism for resistance dissemination and persistence in food animals (Love et al., 2016). The model parameter value distributions were assigned based on the available published literature. However, specific estimates of the relevant parameters for E. coli and CTC in cattle are often lacking and the parameter value distributions built on sparse data from cattle and other species introduce an unknown amount of error into the model outputs. In addition, in vitro data were used when in vivo data were not available for parameterization. This should be considered when interpreting the model outputs. For example, two in vitro estimates of tetracycline fitness cost in E. coli were used to parameterize fitness cost in our model (Nguyen et al., 1989; Ahmad et al., 2015b), resulting in a narrow range of fitness cost values (Table 2B). We performed additional model simulations considering all biologically possible values of fitness cost by specifying Uniform (0, 0.99) for the fitness cost parameter (data not shown). In that case, the median withdrawal period for ADT decreases from 36 days to 12 days.
In conclusion, our model of CTC pharmacokinetics-pharmacodynamics and the enteric bacterial population dynamics in beef cattle demonstrated that a withdrawal period before slaughter can be effective in reducing the population distribution of tetracycline resistant E. coli to pre-treatment levels. However, the withdrawal periods necessary to mitigate this microbiological food safety risk are significantly longer than current CTC withdrawal periods designed to mitigate toxicological food safety risks. The dynamics of resistant E. coli during and after CTC treatment vary by CTC dosage.
CC developed the chlortetracycline pharmacokinetic model, adapted the other model components for this research, parameterized the model, edited the MatLab code, analyzed the model outputs, and drafted this manuscript. LD adapted model components for chlortetracycline pharmacodynamics and tetracycline-resistant Escherichia coli in cattle, wrote the MatLab code, and revised this manuscript. VV developed the pharmacodynamic and Escherichia coli population models, contributed to the pharmacokinetic model development and model parameterization, and revised this manuscript. YG contributed to the development of the models, model parameterization, model output analysis, and revised this manuscript.
This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2016-68003-24607 from the USDA National Institute of Food and Agriculture. CC and LD were supported by Cornell University College of Veterinary Medicine. VV was supported via the Kansas Bioscience Authority funding for the Institute of Computational Comparative Medicine, Kansas State University.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank Gail Ragan (Institute of Computational Comparative Medicine, Kansas State University) for her assistance in identifying the indications, dosages, and withdrawal periods for chlortetracycline in beef cattle in different regions of the world.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01753/full#supplementary-material
CTC, chlortetracycline; ARLA, antimicrobial reduction of liver abscesses; ADC, antimicrobial disease control; ADT, antimicrobial disease treatment; MIC, minimum inhibitory concentration.
2. ^Zoetis “Aureomycin 50 Granular A [package insert]”, in: Chlortetracycline Type A Medicated Article.). http://www.zoetisus.com/contact/pages/product_information/msds_pi/pi/aureomycin_50.pdf.
3. ^Zoetis “ChlorMax 50 [package insert]”, in: Chlortetracycline Type A Medicated Article Antibacterial.). http://www.zoetisus.com/contact/pages/product_information/msds_pi/pi/ChlorMax_50.pdf.
Agga, G. E., Schmidt, J. W., and Arthur, T. M. (2016). Effects of in-feed chlortetracycline prophylaxis in beef cattle on animal health and antimicrobial-resistant Escherichia coli. Appl. Environ. Microbiol. 82, 7197–7204. doi: 10.1128/AEM.01928-16
Ahmad, A., Graesboll, K., Christiansen, L. E., Toft, N., Matthews, L., and Nielsen, S. S. (2015a). Pharmacokinetic-pharmacodynamic model to evaluate intramuscular tetracycline treatment protocols to prevent antimicrobial resistance in pigs. Antimicrob. Agents Chemother. 59, 1634–1642. doi: 10.1128/aac.03919-14
Ahmad, A., Zachariasen, C., Christiansen, L. E., Graesboll, K., Toft, N., Matthews, L., et al. (2015b). Pharmacodynamic modelling of in vitro activity of tetracycline against a representative, naturally occurring population of porcine Escherichia coli. Acta Vet. Scand. 57:79. doi: 10.1186/s13028-015-0169-0
Alexander, T. W., Inglis, G. D., Yanke, L. J., Topp, E., Read, R. R., Reuter, T., et al. (2010). Farm-to-fork characterization of Escherichia coli associated with feedlot cattle with a known history of antimicrobial use. Int. J. Food Microbiol. 137, 40–48. doi: 10.1016/j.ijfoodmicro.2009.11.008
Alexander, T. W., Reuter, T., Sharma, R., Yanke, L. J., Topp, E., and McAllister, T. A. (2009). Longitudinal characterization of resistant Escherichia coli in fecal deposits from cattle fed subtherapeutic levels of antimicrobials. Appl. Environ. Microbiol. 75, 7125–7134. doi: 10.1128/aem.00944-09
Alexander, T. W., Yanke, L. J., Topp, E., Olson, M. E., Read, R. R., Morck, D. W., et al. (2008). Effect of subtherapeutic administration of antibiotics on the prevalence of antibiotic-resistant Escherichia coli in feedlot cattle. Appl. Environ. Microbiol. 74, 4405–4416. doi: 10.1128/aem.00489-08
Arikan, O. A. (2008). Degradation and metabolization of chlortetracycline during the anaerobic digestion of manure from medicated calves. J. Hazard. Mater. 158, 485–490. doi: 10.1016/j.jhazmat.2008.01.096
Arikan, O. A., Mulbry, W., and Rice, C. (2009). Management of antibiotic residues from agricultural sources: use of composting to reduce chlortetracycline residues in beef manure from treated animals. J. Hazard. Mater. 164, 483–489. doi: 10.1016/j.jhazmat.2008.08.019
Aslam, M., Greer, G. G., Nattress, F. M., Gill, C. O., and McMullen, L. M. (2004). Genetic diversity of Escherichia coli recovered from the oral cavity of beef cattle and their relatedness to faecal E. coli. Lett. Appl. Microbiol. 39, 523–527. doi: 10.1111/j.1472-765X.2004.01620.x
Aslam, M., Nattress, F., Greer, G., Yost, C., Gill, C., and McMullen, L. (2003). Origin of Contamination and Genetic Diversity of Escherichia coli in Beef Cattle. Appl. Environ. Microbiol. 69, 2794–2799. doi: 10.1128/aem.69.5.2794-2799.2003
Australian Pesticides and Veterinary Medicines Authority (2015). Public Chemical Registration Information System Search. Available online at: portal.apvma.gov.au/pubcris
Bahl, M. I., Hansen, L. H., Licht, T. R., and Sorensen, S. J. (2004). In vivo detection and quantification of Tetracycline by use of a whole-cell biosensor in the rat intestine. Antimicrob. Agents Chemother. 48, 1112–1117. doi: 10.1128/aac.48.4.1112-1117.2004
Benedict, K. M. (2011). Antimicrobial Resistance Surveillance in Feedlot Cattle. Ph.D. Dissertation, Colorado State University. Available online at: dspace.library.colostate.edu/bitstream/handle/10217/47323/Benedict_colostate_0053A_10312.pdf?sequence=1
Berendonk, T. U., Manaia, C. M., Merlin, C., Fatta-Kassinos, D., Cytryn, E., Walsh, F., et al. (2015). Tackling antibiotic resistance: the environmental framework. Nat. Rev. Microbiol. 13, 310–317. doi: 10.1038/nrmicro3439.
Bradley, B., Allen, E., Showalter, D., and Colaianne, J. (1982). Comparative pharmacokinetics of chlortetracycline in milk-fed versus conventionally fed calves. J. Vet. Pharmacol. Ther. 5, 267–278. doi: 10.1111/j.1365-2885.1982.tb00442.x
Branham, L. A. (2007). Antimicrobial Susceptibility of Generic Escherichia coli Following Administration of Subtherapeutic Antimicrobial Drugs to Feedlot Cattle. Ph.D. Dissertation, Texas Tech University. Available online at: http://ttu-ir.tdl.org/ttu-ir/bitstream/handle/2346/17388/Branham_Loree_Diss.pdf?sequence=1
Call, D. R., Matthews, L., Subbiah, M., and Liu, J. (2013). Do antibiotic residues in soils play a role in amplification and transmission of antibiotic resistant bacteria in cattle populations? Front. Microbiol. 4:193. doi: 10.3389/fmicb.2013.00193
Callaway, T. R., Anderson, R. C., Genovese, K. J., Poole, T. L., Anderson, T. J., Byrd, J. A., et al. (2002). Sodium chlorate supplementation reduces E. coli O157:H7 populations in cattle. J. An. Sci. 80:1683. doi: 10.2527/2002.8061683x
Canadian Food Inspection Agency (2014). “Chlortetracycline hydrochloride and sulfamethazine - MIB #49,” in Compendium of Medicating Ingredient Brochures. Available online at: http://www.inspection.gc.ca/animals/feeds/medicating-ingredients/mib/mib-49/eng/1331048641123/1331048701288
Canadian Food Inspection Agency (2016). “Chlortetracycline hydrochloride - MIB #34,” in Compendium of Medicating Ingredient Brochures. Available online at: http://www.inspection.gc.ca/animals/feeds/medicating-ingredients/mib/mib-34/eng/1330984471781/1330984546791
Carlson, J. C., and Mabury, S. A. (2006). Dissipation kinetics and mobility of chlortetracycline, tylosin, and monensin in an agricultural soil in Northumberland County, Ontario, Canada. Environ. Toxicol. Chemistry 25:1. doi: 10.1897/04-657r.1
Carson, C. A., Reid-Smith, R., Irwin, R. J., Martin, W. S., and McEwen, S. A. (2008). Antimicrobial resistance in generic fecal Escherichia coli from 29 beef farms in Ontario. Can. J. Veter. Res. 72, 119–128.
Cazer, C. L., Volkova, V. V., and Grohn, Y. T. (2014). Use of pharmacokinetic modeling to assess antimicrobial pressure on enteric bacteria of beef cattle fed chlortetracycline for growth promotion, disease control, or treatment. Foodborne Pathog. Dis. 11, 403–411. doi: 10.1089/fpd.2013.1677
Chipotle (2013). Chipotle Responsibly Raised Beef Standards Remain Unchanged. Available online at: ir.chipotle.com/phoenix.zhtml?c=194775&p=irol-newsArticle&ID=1847311
Clinical and Laboratory Standards Institute (2015). VET01S: Performance Standards for Antimicrobial Disk and Dilution Susceptibility Tests for Bacteria Isolated from Animals. Available online at: vet01s.edaptivedocs.info/dashboard.aspx
Daniels, J. B., Call, D. R., Hancock, D., Sischo, W. M., Baker, K., and Besser, T. E. (2009). Role of ceftiofur in selection and dissemination of blaCMY-2-mediated cephalosporin resistance in Salmonella enterica and commensal Escherichia coli isolates from cattle. Appl. Environ. Microbiol. 75, 3648–3655. doi: 10.1128/AEM.02435-08
De Filippis, F., La Storia, A., Villani, F., and Ercolini, D. (2013). Exploring the sources of bacterial spoilers in beefsteaks by culture-independent high-throughput sequencing. PLoS ONE 8:e70222. doi: 10.1371/journal.pone.0070222
DeMars, Z., Biswas, S., Amachawadi, R. G., Renter, D. G., and Volkova, V. V. (2016). Antimicrobial susceptibility of enteric gram negative facultative anaerobe bacilli in aerobic versus anaerobic conditions. PLoS ONE 11:e0155599. doi: 10.1371/journal.pone.0155599
Durso, L. M., Smith, D., and Hutkins, R. W. (2004). Measurements of fitness and competition in commensal Escherichia coli and E. coli O157: H7 strains. Appl. Environ. Microbiol. 70, 6466–6472. doi: 10.1128/AEM.70.11.6466-6472.2004
EU Directorate General Health and Consumers (2017). Community Register of Veterinary Medicinal Products. Available online at: ec.europa.eu/health/documents/community-register/html/vatc.htm
Forslund, K., Sunagawa, S., Kultima, J. R., Mende, D. R., Arumugam, M., Typas, A., et al. (2013). Country-specific antibiotic use practices impact the human gut resistome. Genome Res. 23, 1163–1169. doi: 10.1101/gr.155465.113
Freter, R., Brickner, H., Botney, M., Cleven, D., and Aranki, A. (1983a). Mechanisms that control bacterial populations in continuous-flow culture models of mouse large intestinal flora. Infect. Immun. 39, 676–685.
Goutelle, S., Maurin, M., Rougier, F., Barbaut, X., Bourguignon, L., Ducher, M., et al. (2008). The Hill equation: a review of its capabilities in pharmacological modelling. Fundam. Clin. Pharmacol. 22, 633–648. doi: 10.1111/j.1472-8206.2008.00633.x
Graesboll, K., Nielsen, S. S., Toft, N., and Christiansen, L. E. (2014). How fitness reduced, antimicrobial resistant bacteria survive and spread: a multiple pig-multiple bacterial strain model. PLoS ONE 9:e100458. doi: 10.1371/journal.pone.0100458
Gullberg, E., Cao, S., Berg, O. G., Ilbäck, C., Sandegren, L., Hughes, D., et al. (2011). Selection of resistant bacteria at very low antibiotic concentrations. PLoS Pathog. 7:e1002158. doi: 10.1371/journal.ppat.1002158
Halling-Sørensen, B., Sengeløv, G., and Tjørnelund, J. (2002). Toxicity of tetracyclines and tetracycline degradation products to environmentally relevant bacteria, including selected tetracycline-resistant bacteria. Arch. Environ. Contam. Toxicol. 42, 263–271. doi: 10.1007/s00244-001-0017-2
Hoffmann La-Roche Inc. (1997). NADA 048-761 Aureomycin® - Supplemental Approval. US Food and Drug Administration. Available online at: http://www.fda.gov/AnimalVeterinary/Products/ApprovedAnimalDrugProducts/FOIADrugSummaries/ucm064717.htm
Kanwar, N., Scott, H. M., Norby, B., Loneragan, G. H., Vinasco, J., McGowan, M., et al. (2013). Effects of ceftiofur and chlortetracycline treatment strategies on antimicrobial susceptibility and on tet(A), tet(B), and bla CMY-2 resistance genes among E. coli isolated from the feces of feedlot cattle. PLoS ONE 8:e80575. doi: 10.1371/journal.pone.0080575
Kjeldsen, N., and Callesen, J. (2006). “Terminated use of antimicrobial growth promoters in pig production in Denmark–effects on pig welfare and productivity,” in BOKU-Symposium Tierernährung ohne Antibiotische Leistungsförderer, eds D. Barug, D. de Jong, A.K. Kies, and M. W. A. Verstegen (Vienna: Universität für Bodenkultur Wien: The National Committee for Pig Production), 127–135.
Kluytmans, J. A., Overdevest, I. T., Willemsen, I., Kluytmans-van den Bergh, M. F., van der Zwaluw, K., Heck, M., et al. (2013). Extended-spectrum beta-lactamase-producing Escherichia coli from retail chicken meat and humans: comparison of strains, plasmids, resistance genes, and virulence factors. Clin. Infect. Dis. 56, 478–487. doi: 10.1093/cid/cis929
Krause, D. O., Smith, W. J., Conlan, L. L., Gough, J. M., Williamson, M. A., and McSweeney, C. S. (2003). Diet influences the ecology of lactic acid bacteria and Escherichia coli along the digestive tract of cattle: neural networks and 16S rDNA. Microbiology 149(Pt. 1), 57–65. doi: 10.1099/mic.0.25685-0
Laven, R. A., Ashmore, A., and Stewart, C. S. (2003). Escherichia coli in the rumen and colon of slaughter cattle, with particular reference to E. coli O157. Veter. J. 165, 78–83. doi: 10.1016/s1090-0233(02)00162-4
Love, W. J., Zawack, K. A., Booth, J. G., Grhn, Y. T., and Lanzas, C. (2016). Markov networks of collateral resistance: national antimicrobial resistance monitoring system surveillance results from Escherichia coli isolates, 2004-2012. PLoS Comput. Biol. 12:e1005160. doi: 10.1371/journal.pcbi.1005160
Lowrance, T. C., Loneragan, G. H., Kunze, D. J., Platt, T. M., Ives, S. E., Scott, H. M., et al. (2007). Changes in antimicrobial susceptibility in a population of Escherichia coli isolated from feedlot cattle administered ceftiofur crystalline-free acid. Am. J. Vet. Res. 68, 501–507. doi: 10.2460/ajvr.68.5.501
Martin, C., Philippeau, C., and Michalet-Doreau, B. (1999). Effect of wheat and corn variety on fiber digestion in beef steers fed high-grain diets. J. Anim. Sci. 77, 2269–2278. doi: 10.2527/1999.7782269x
McDonalds (2016). McDonald's USA Announces Big Changes to Its Food. Available online at: news.mcdonalds.com/US/releases/McDonald-s-USA-Announces-Big-Changes-to-its-Food
McGowan, M. T. (2014). The Epidemiology of Tetracycline and Ceftiofur Resistance in Commensal Escherichia coli. MS Dissertation, Kansas State University. Available online at: https://hdl.handle.net/2097/17641
Morley, P. S., Dargatz, D. A., Hyatt, D. R., Dewell, G. A., Patterson, J. G., Burgess, B. A., et al. (2011). Effects of restricted antimicrobial exposure on antimicrobial resistance in fecal Escherichia coli from feedlot cattle. Foodborne Pathog. Dis. 8, 87–98. doi: 10.1089/fpd.2010.0632
Mouton, J. W., and Vinks, A. A. (2005). Pharmacokinetic/pharmacodynamic modelling of antibacterials in vitro and in vivo using bacterial growth and kill kinetics: the minimum inhibitory concentration versus stationary concentration. Clin. Pharmacokinet. 44, 201–210. doi: 10.2165/00003088-200544020-00005
New Zealand Ministry for Primary Industries (2017). Agricultural Compound and Veterinary Medicines Register. Available online at: eatsafe.nzfsa.govt.nz/web/public/acvm-register
Nguyen, T., Phan, Q. G., Duong, L. P., Bertrand, K. P., and Lenski, R. E. (1989). Effects of carriage and expression of the Tn10 tetracycline-resistance operon on the fitness of Escherichia coli K12. Mol. Biol. Evol. 6, 213–225. doi: 10.1093/oxfordjournals.molbev.a040545
Norcia, L. J., Silvia, A. M., and Hayashi, S. F. (1999). Studies on time-kill kinetics of different classes of antibiotics against veterinary pathogenic bacteria including Pasteurella, Actinobacillus and Escherichia coli. J. Antibiot. 52, 52–60. doi: 10.7164/antibiotics.52.52
Overdevest, I., Willemsen, I., Rijnsburger, M., Eustace, A., Xu, L., Hawkey, P., et al. (2011). Extended-spectrum beta-lactamase genes of Escherichia coli in chicken meat and humans, The Netherlands. Emerg. Infect. Dis. 17, 1216–1222. doi: 10.3201/eid1707.110209
Platt, T. M., Lonergan, G. H., Scott, M., Norby, B., Thomson, D. U., Brown, M. S., et al. (2008). Antimicrobial susceptibility of enteric bacteria recovered from feedlot cattle administered chlortetracycline in feed. Am. J. Vet. Res. 69, 988–996. doi: 10.2460/ajvr.69.8.988
Rao, S., Van Donkersgoed, J., Bohaychuk, V., Besser, T., Song, X. M., Wagner, B., et al. (2010). Antimicrobial drug use and antimicrobial resistance in enteric bacteria among cattle from Alberta feedlots. Foodborne Pathog. Dis. 7, 449–457. doi: 10.1089/fpd.2009.0400
Reinbold, J. B., Coetzee, J. F., Gehring, R., Havel, J. A., Hollis, L. C., Olson, K. C., et al. (2010). Plasma pharmacokinetics of oral chlortetracycline in group fed, ruminating, Holstein steers in a feedlot setting. J. Vet. Pharmacol. Ther. 33, 76–83. doi: 10.1111/j.1365-2885.2009.1116.x
Sargeant, J. M., Amezcua, M. R., Rajic, A., and Waddell, L. (2007). Pre-harvest interventions to reduce the shedding of E. coli O157 in the faeces of weaned domestic ruminants: a systematic review. Zoonoses Public Health 54, 260–277. doi: 10.1111/j.1863-2378.2007.01059.x
Sharma, R., Munns, K., Alexander, T., Entz, T., Mirzaagha, P., Yanke, L. J., et al. (2008). Diversity and distribution of commensal fecal Escherichia coli bacteria in beef cattle administered selected subtherapeutic antimicrobials in a feedlot setting. Appl. Environ. Microbiol. 74, 6178–6186. doi: 10.1128/aem.00704-08
Shaver, R., Nytes, A., Satter, L., and Jorgensen, N. (1986). Influence of amount of feed intake and forage physical form on digestion and passage of prebloom alfalfa hay in dairy cows. J. Dairy Sci. 69, 1545–1559. doi: 10.3168/jds.S0022-0302(86)80571-9
Simonsen, L., Gordon, D. M., Stewart, F. M., and Levin, B. R. (1990). Estimating the rate of plasmid transfer: an end-point method. J. Gen. Microbiol. 136, 2319–2325. doi: 10.1099/00221287-136-11-2319
Speer, B. S., Shoemaker, N. B., and Salyers, A. A. (1992). Bacterial resistance to tetracycline: mechanisms, transfer, and clinical significance. Clin. Microbiol. Rev. 5, 387–399. doi: 10.1128/CMR.5.4.387
Tyson Foods (2015). Tyson Foods Strives to Eliminate Human Antibiotics from Broiler Chicken Flocks by 2017. Available online at: ir.tyson.com/investor-relations/news-releases/news-releases-details/2015/Tyson-Foods-Strives-to-Eliminate-Human-Antibiotics-From-Broiler-Chicken-Flocks-by-2017/default.aspx
US Food and Drug Administration (2003). Guidance for Industry #152: Evaluating the Safety of Antimicrobial New Animal Drugs with Regard to Their Microbiological Effects on Bacteria of Human Health Concern. Rockville, MD. Available online at: http://www.fda.gov/downloads/AnimalVeterinary/GuidanceComplianceEnforcement/GuidanceforIndustry/UCM052519.pdf
US Food and Drug Administration (2012). Guidance for Industry #209: The Judicious use of Medically Important Antimicrobial Drugs in Food-Producing Animals. Rockville, MD. Available online at: http://www.fda.gov/downloads/AnimalVeterinary/GuidanceComplianceEnforcement/GuidanceforIndustry/UCM216936.pdf
US Food and Drug Administration (2013). Guidance for Industry #213, New Animal Drugs and New Animal Drug Combination Products Administered in or on Medicated Feed or Drinking Water of Food-Producing Animals: Recommendations for Drug Sponsors for Voluntarily Aligning Product Use Conditions with GFI# 209. Rockville, MD. Available online at: http://www.fda.gov/downloads/AnimalVeterinary/GuidanceComplianceEnforcement/GuidanceforIndustry/UCM299624.pdf
United States Department of Agriculture (2013). Feedlot 2011 Part IV: Health and Health Management on US Feedlots with a Capacity of 1,000 or More Head. Fort Collins, CO. Available online at: http://www.aphis.usda.gov/animal_health/nahms/feedlot/downloads/feedlot2011/Feed11_dr_PartIV.pdf
UK Veterinary Medicines Directorate (2017). Product Information Database. Available online at: http://www.vmd.defra.gov.uk/ProductInformationDatabase/Default.aspx
Volkova, V. V., Cazer, C. L., and Grohn, Y. T. (2017). Models of antimicrobial pressure on intestinal bacteria of the treated host populations. Epidemiol. Infect. 145, 1–14. doi: 10.1017/S095026881700084X
Volkova, V. V., KuKanich, B., and Riviere, J. E. (2016). Exploring post-treatment reversion of antimicrobial resistance in enteric bacteria of food animals as a resistance mitigation strategy. Foodborne Pathog. Dis. 13, 610–617. doi: 10.1089/fpd.2016.2152
Volkova, V. V., Lanzas, C., Lu, Z., and Grohn, Y. T. (2012). Mathematical model of plasmid-mediated resistance to ceftiofur in commensal enteric Escherichia coli of cattle. PLoS ONE 7:e36738. doi: 10.1371/journal.pone.0036738
Volkova, V. V., Lu, Z., Lanzas, C., and Grohn, Y. T. (2013). Evaluating targets for control of plasmid-mediated antimicrobial resistance in enteric commensals of beef cattle: a modelling approach. Epidemiol. Infect. 141, 2294–2312. doi: 10.1017/S0950268812002993
Wagner, B. A., Dargatz, D. A., Salman, M. D., Morley, P. S., Wittum, T. E., and Keefe, T. J. (2002). Comparison of sampling techniques for measuring the antimicrobial susceptibility of enteric Escherichia colirecovered from feedlot cattle. Am. J. Vet. Res. 63, 1662–1670. doi: 10.2460/ajvr.2002.63.1662
Wen, X., Gehring, R., Stallbaumer, A., Riviere, J. E., and Volkova, V. V. (2016). Limitations of MIC as sole metric of pharmacodynamic response across the range of antimicrobial susceptibilities within a single bacterial species. Sci. Rep. 6:37907. doi: 10.1038/srep37907
Yan, J. J., Hong, C. Y., Ko, W. C., Chen, Y. J., Tsai, S. H., Chuang, C. L., et al. (2004). Dissemination of blaCMY-2 among Escherichia coli isolates from food animals, retail ground meats, and humans in Southern Taiwan. Antimicrobial Agents Chemother. 48, 1353–1356. doi: 10.1128/aac.48.4.1353-1356.2004
Zebeli, Q., Tafaj, M., Weber, I., Dijkstra, J., Steingass, H., and Drochner, W. (2007). Effects of varying dietary forage particle size in two concentrate levels on chewing activity, ruminal mat characteristics, and passage in dairy cows. J. Dairy Sci. 90, 1929–1942. doi: 10.3168/jds.2006-354
Zoetis (2016). N048-761: Chlortetracycline Beef Cattle Feed – AUREOa – LA>400 lb Type B Medicated Beef Cattle Feed. Blue Bird Feed Company. Food & Drug Association. Available online at: http://www.zoetisus.com/products/vfd/cattle-pdfs/23.revised_2016-6-2-aureomycin-_048-761__cattle-over-400-lb-la__type-b_.pdf
Keywords: beef cattle, antibiotic resistance, enteric bacteria, food-borne pathogens, mathematical modeling, population pharmacokinetics, pharmacodynamics
Citation: Cazer CL, Ducrot L, Volkova VV and Gröhn YT (2017) Monte Carlo Simulations Suggest Current Chlortetracycline Drug-Residue Based Withdrawal Periods Would Not Control Antimicrobial Resistance Dissemination from Feedlot to Slaughterhouse. Front. Microbiol. 8:1753. doi: 10.3389/fmicb.2017.01753
Received: 02 May 2017; Accepted: 29 August 2017;
Published: 20 September 2017.
Edited by:Ed Topp, Agriculture and Agri-Food Canada, Canada
Copyright © 2017 Cazer, Ducrot, Volkova and Gröhn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Casey L. Cazer, email@example.com
†Present Address: Lucas Ducrot, Ecole Centrale de Nantes, Nantes, France