Hostname: page-component-7c8c6479df-24hb2 Total loading time: 0 Render date: 2024-03-28T15:08:48.078Z Has data issue: false hasContentIssue false

Spatial and temporal patterns in antimicrobial resistance of Salmonella Typhimurium in cattle in England and Wales

Published online by Cambridge University Press:  03 January 2012

R. COX*
Affiliation:
National Centre for Zoonosis Research, Leahurst, University of Liverpool, Wirral, UK Centre for Veterinary Epidemiological Research, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Prince Edward Island, Canada
T. SU
Affiliation:
Department of Mathematics and Statistics, Lancaster University, Lancaster, UK
H. CLOUGH
Affiliation:
National Centre for Zoonosis Research, Leahurst, University of Liverpool, Wirral, UK
M. J. WOODWARD
Affiliation:
Department of Bacteriology, Animal Health and Veterinary Laboratories Agency, Addlestone, Surrey, UK
C. SHERLOCK
Affiliation:
National Centre for Zoonosis Research, Leahurst, University of Liverpool, Wirral, UK Department of Mathematics and Statistics, Lancaster University, Lancaster, UK
*
*Author for correspondence: Dr R. Cox, Centre for Veterinary Epidemiological Research, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Prince Edward Island, Canada. (Email: rucox@upei.ca)
Rights & Permissions [Opens in a new window]

Summary

Salmonella is the second most commonly reported human foodborne pathogen in England and Wales, and antimicrobial-resistant strains of Salmonella are an increasing problem in both human and veterinary medicine. In this work we used a generalized linear spatial model to estimate the spatial and temporal patterns of antimicrobial resistance in Salmonella Typhimurium in England and Wales. Of the antimicrobials considered we found a common peak in the probability that an S. Typhimurium incident will show resistance to a given antimicrobial in late spring and in mid to late autumn; however, for one of the antimicrobials (streptomycin) there was a sharp drop, over the last 18 months of the period of investigation, in the probability of resistance. We also found a higher probability of resistance in North Wales which is consistent across the antimicrobials considered. This information contributes to our understanding of the epidemiology of antimicrobial resistance in Salmonella.

Type
Original Papers
Copyright
Copyright © Cambridge University Press 2012

INTRODUCTION

Salmonella is the second most commonly reported human foodborne pathogen in England and Wales, with more than 10 000 laboratory-confirmed cases in 2009 [1]. The majority of these cases involve foodborne transmission, mostly of animal origin [2]. Antimicrobial-resistant (AMR) serotypes of Salmonella are globally widespread [Reference Threlfall3] and are an increasing problem in human and veterinary medicine [Reference Threlfall4]; moreover, multidrug resistance (MDR; i.e. resistance to ⩾4 antimicrobials) has commonly been recorded [Reference Threlfall3, 5]. Despite legislation to control antimicrobial use, prevalence of AMR Salmonella isolates has increased in developed and developing countries in recent years [Reference Threlfall3, Reference Meakins6]. AMR serotypes can be detrimental to animal health and productivity [5, Reference Miller, McNamara and Singer7] and the resulting disease can be severe. Clinical salmonellosis in cattle, for example, can cause acute diarrhoea, abortion, decreased milk production and high mortality [Reference Smith and Smith8]. It results in economic loss to herd owners and impacts on future trading opportunities [Reference Miller, McNamara and Singer7]. AMR serotypes in livestock can cause human infections through the food chain [Reference Threlfall3, Reference Molbak9] and therefore also have an adverse impact on human health and welfare [Reference Frenzen10, 11], e.g. through failed treatment or prolonged illness [Reference Molbak9], increased hospitalization [Reference Varma12], or increased mortality [Reference Helms13].

The majority of salmonellosis incidents are caused by relatively few serotypes. Serotypes S. Typhimurium and S. Enteritidis account for 60–80% of all human salmonellosis [1]. S. Typhimurium is a particularly epidemic serotype with a ubiquitous host range that is commonly responsible for clinical disease in both livestock and humans [5]. In cattle, for example, it accounts for more than 7% of incidents and is the second most common serotype after S. Dublin. (Although S. Dublin causes more than 70% of outbreaks it is rarely associated with human foodborne infection [5].) Antimicrobial resistance is a major problem in S. Typhimurium compared to other Salmonella serotypes [Reference Prescott, Prescott, Baggot and Walker14]. In England and Wales at least 70% of isolates from livestock are resistant to one or more antimicrobials [5]; in humans the proportion is greater than 80% [Reference Threlfall4]. Multidrug-resistant S. Typhimurium was first identified in the UK in the 1960s and a number of different phage types have caused serious epidemics since then [Reference Threlfall3], e.g. multiresistant definitive phage type (DT) 104 [Reference Threlfall15].

Previous research has highlighted that there is limited understanding about certain aspects of the development of AMR serotypes, including risk factors for carriage of resistant organisms [Reference Warnick16], seasonal prevalence [Reference Pangloli17] and the emergence of multidrug-resistant serotypes [11, Reference Lanzas18]. This lack of knowledge hinders attempts to develop effective targeted Salmonella control programmes in the UK and worldwide. Critical to this is a full understanding of the spatial and temporal development of AMR serotypes in livestock populations. Similar antimicrobial resistance patterns occurring at similar periods in time (‘temporal components’) or locations (‘spatial components’) may reflect the effect of explanatory variables which themselves are structured in time and/or space [Reference Diggle19] or might provide evidence of a contagious mechanism. Large-scale regional variations in infections may indicate large-scale risk factors, e.g. use of certain antimicrobials in certain regions or at certain time of the year, whereas small-scale patterns, e.g. clustering of farms with resistance to the same antimicrobial may indicate a local risk factor, e.g. environmental contamination [Reference Murray, Wray and Wray20]. While some explanatory variables may be recorded, others may currently be unmeasured or unknown. However, quantitative description of such spatial and temporal patterns will inform knowledge of the underlying epidemiology and biological processes [Reference Fenton21].

As a step towards the goal of better understanding the development of antimicrobial resistance, we investigated temporal and spatial patterns of antimicrobial resistance in S. Typhimurium isolates from one livestock sector: cattle in England and Wales. We focused on three particular antimicrobials: streptomycin, sulphonamide compounds, and chloramphenicol, as well as examining MDR. Specifically, we address the following questions:

  • Is there evidence of changes from year to year in the probability of antimicrobial resistance being observed in a given incident?

  • Is there evidence of a seasonal effect in the occurrence of antimicrobial resistance?

  • Is there evidence of geographical variation across England and Wales in the probability of resistant organisms being observed in a given incident?

  • What are the similarities and differences in the temporal and spatial resistance patterns for the different antimicrobials?

MATERIALS AND METHODS

Data were obtained from the Veterinary Laboratories Agency (VLA) ‘Farmfile’ database [1], which documents all livestock incidents of Salmonella in the UK. The database is one of the UK's most extensive livestock databases (for review see [Reference Gibbens22]) and includes details of each incident's location, date, Salmonella serotype and antimicrobial sensitivity. It is a passive surveillance system; samples are submitted to the VLA regional laboratories by veterinary practitioners. Salmonella is a notifiable disease in the UK and therefore all cases are reported. The database does include reports from statutory monitoring and surveillance; however, the majority of these are from poultry flocks where surveillance is common practice in contrast to the majority of reports from other species which are the result of examinations of clinical disease.

Data from 1 January 2003 to 31 December 2006 were used in this analysis. Data were restricted to this time period to strike a compromise between being able to look for seasonality and compatibility of data over time: in January 2003 the VLA established the current version of the ‘Farmfile’ database which integrated Salmonella recording and reporting [5]. Data were selected according to reason for submission to the VLA. We included any incident that was reported as a result of examinations performed to diagnose clinical disease and excluded any reports that resulted from statutory monitoring or surveillance activities. We focused our analysis on cattle, a species in which isolations must by law be reported and a known reservoir for antimicrobial resistance to S. Typhimurium [5]. An incident of Salmonella was defined as the first isolation and any subsequent isolations of the same serotype of a particular Salmonella, following diagnosis of a clinical case, from an animal, group of animals or their environment on a single premises, within 30 days [5].

Samples from each incident had been tested for their in vitro sensitivity to 16 antimicrobials using the British Society for Antimicrobial Chemotherapy (BSAC) disc diffusion technique (www.bsac.org.uk) on Oxoid ‘Isosensitest’ agar. The method is described in detail in a Defra publication [5]. The antimicrobials (with disc concentrations in parentheses) were ampicillin (10 μg), amoxicillin (30 μg), streptomycin (25 μg), sulphonamide compounds (comprising of 37% sulphadiazine; 37% sulphathiazole and 26% sulphamerazine) (300 μg), chloramphenicol (10 μg), sulphamethoxozole (25 μg), nalidixic acid (30 μg), tetracycline (10 μg), neomycin (10 μg), furazolidone (15 μg), ceftazidime (30 μg), amikacin (30 μg), gentamicin (10 μg), cefotaxime (30 μg), apramycin (15 μg) and ciprofloxacin (1 μg). Separate tests for extended spectrum β-lactamases were not routinely performed until after 2006 and so these data are not included. The choice of antimicrobials, which is reviewed periodically, is designed to comprise a core set which has been used in veterinary practice for many years, some of the more recently licensed antimicrobials and some of limited usage in animals in Great Britain which are used in other European countries. Analysis excluded any incidents where the samples were not tested because there was no approved sensitivity test. We refer to an incident of S. Typhimurium as a resistant case if the sample from this incident shows resistance to a certain antimicrobial; otherwise it is a susceptible case. We refer to an incident as multidrug resistant if a sample was resistant to ⩾4 antimicrobial agents in the panel of 16 [5].

Exploratory data analysis

Between January 2003 and December 2006 there were 294 incidents of S. Typhimurium in England and Wales on 256 farms, the locations of which are shown in Figure 1. Of these farms 226 experienced a single incident, 26 farms had two incidents, one farm had three incidents, two farms had four incidents, and one farm had five incidents. The largest distance from a farm to its nearest neighbour was 67·9 km. For data protection purposes in all spatial plots the coordinates of each individual farm incident have been randomized uniformly to a disk of radius 10 km centred on the true location.

Fig. 1. Locations of the 256 farms which experienced at least one incident of S. Typhimurium between January 2003 and December 2006.

Table 1 presents the frequencies of S. Typhimurium incidents which were resistant or susceptible to each antimicrobial. None of the incidents showed evidence of resistance to amoxicillin, ceftazidime, amikacin, gentamicin, apramycin, or ciprofloxacin; resistance to cefotaxime, furazolidone, and neomycin was found in only one, one, and four of the incidents, respectively, over the 4-year period. These antimicrobials are therefore not included.

Table 1. Number of Salmonella Typhimurium incidents per year and frequency of antimicrobial resistance to each antimicrobial (only antimicrobials for which there was at least one incident of resistance are listed)

* Resistant to ⩾4 of the 16 antimicrobials tested.

We were interested in discerning patterns in the spatial or temporal variability of resistance and in comparing the patterns for antimicrobials between which a genetic link is suspected as well as in comparing compounds where no such relationship is believed to exist. For this reason our analysis focused on streptomycin, sulphonamide compounds, and chloramphenicol, to which 62%, 77%, and 70% of the cattle isolates demonstrated resistance, respectively.

These antimicrobials often occur in the characteristic multiresistant pentavalent-resistant pattern (resistance to tetracycline, ampicillin, chloramphenicol, streptomycin and sulphonamide compounds) due to chromosomal integration [Reference Davison23]. We could not analyse antimicrobials to which isolates were either almost entirely susceptible or almost entirely resistant. Therefore we could not include antimicrobials which are not part of the multiresistant pentavalent pattern, and which are unlinked in their mode of resistance, e.g. nalidixic acid. Independent spatiotemporal analyses was performed on the relative incidences of resistant and susceptible cases for each of the three antimicrobial datasets identified above and for MDR.

Figure 2 maps all of the incidents of S. Typhimurium indicating the presence or absence of bacteria susceptible to chloramphenicol. Superficially, at least, plots for streptomycin, sulphonamide compounds, and MDR appear very similar, with, e.g. some separation between groups of resistant samples and groups of non-resistant samples, especially in Wales and the Southwest. Exploratory fitting of Generalized Additive Models (e.g. [Reference Hastie and Tibshirani24]), showed a likely spatial pattern, motivating the need to allow for spatial correlation in the full statistical analysis.

Fig. 2. Incidents of S. Typhimurium in cattle farms in England and Wales that were resistant to chloramphenicol between 1 January 2003 and 31 December 2006. Incidents with bacteria resistant to chloramphenicol are indicated by an open symbol (○), while incidents susceptible to chloramphenicol are indicated by a cross (×).

Statistical framework

We focused on the probability that an observed outbreak of S. Typhimurium at time t and farm i is resistant to ‘A’, where ‘A’ denotes one of the antimicrobials of interest, or MDR. This quantity, which we denote p t A (t) is our main focus since we are interested in the pattern of resistance in observed incidents of S. Typhimurium. We were particularly interested in how the chance that a given incident is resistant to ‘A’ varies across England and Wales as well as over the 4-year period, and so we express this probability as a function of the spatial location of farm i, x i ; i.e. p A (x i ,t)=p t A (t). We studied each antimicrobial (and MDR) in turn, and drop the superscript to ease notation.

A Generalized Linear Spatial Model (GLSM; e.g. [Reference Diggle and Ribeiro25]) with a logit (i.e. log-odds) link for p(x i , t) is as follows,

Temporal trends are included through the deterministic covariate term f(t) as in a traditional Generalized Linear Model (GLM); however, the GLSM also allows for spatial variation via a spatially structured random-effects term S(x). The time-varying contribution to the log-odds that an incident is resistant to ‘A’ is assumed to be the same for all farms. Conversely, we assume that there is a spatially varying contribution to the log-odds that an incident is resistant to ‘A’, but that this does not change over the 4 years.

We model S as a Gaussian process with mean 0, and variance σ2. This model allows the log-odds that S. Typhimurium incidents at two different farms are resistant to a given antimicrobial to be correlated. For a Gaussian process, at any given farm location x, the value of S, is a realization from a Gaussian distribution, i.e. S(x)~N(0, σ2). Both σ2 and S(x) are unknown, and are estimated during the model-fitting process.

The correlation between farms is assumed to decay exponentially with distance; for farms at locations x and x′,

(1)

Correlation of this kind allows S(x) and S(x′) to be positively linked, with the link growing weaker as the distance between x and x′ increases. The scale over which the spatial dependence tapers off is also estimated during the model-fitting process.

The correlation between the surface at each farm and any arbitrary point can be calculated using equation (1); we may thus characterize the distribution of the surface at, e.g. a fine grid of points across the country in terms of the values at the observation points, and hence provide a point estimate of the value and a measure of the uncertainty at each point on the grid. For more details of this process, which is known as kriging, see e.g. [Reference Diggle and Ribeiro25]. In the current context these kriging predictions make the most sense at points where there are cattle farms; at other points they simply indicate the likely variation in risk due to location that would occur if a cattle farm were to be situated there.

We used the package geoRglm in the R (www.r-project.org) environment, which employs Markov Chain Monte Carlo (MCMC) methods under the Bayesian paradigm. The Bayesian approach requires a description of prior beliefs about the model parameters. The choice for each prior is as follows.

  1. (1) Covariate parameters β were assumed to be independent a priori, and each was given a vague Gaussian prior, N(0, 10).

  2. (2) The range parameter φ was given a discrete prior with possible values of 1 km, 2 km, …, 200 km, and prior probability proportional to exp[−φ/20]. This favours low φ values (i.e. short-range spatial correlation) and forces the data to assert the existence of any real spatial correlation.

  3. (3) The variance σ2 was given a scaled inverse χ2 prior on 6 d.f. with scale value 1 (so that the expected value of 1/σ2 is 1). It is well known (e.g. [Reference Hastie and Tibshirani24, Reference Diggle and Ribeiro25]) that φ and σ2 can be very strongly correlated a posteriori and that this can make individual identification of φ and σ2 difficult. It is often only σ/φ½ which is well identified, but fortunately it is also this combination which is important in predicting the underlying surface S(x). To avoid large, unnecessary values for σ2 and poor mixing of the MCMC algorithm we truncated the prior, allowing only values of σ2⩽10.

RESULTS

Choice of temporal covariates to be used in GLSM

In a GLSM analysis the temporal effect is considered through a set of time-varying covariates. We considered a number of possible temporal covariate effects for each antimicrobial, and for MDR, via the simple logistic regression GLM. The temporal covariate effects were: a linear trend, a quadratic trend, continuous piecewise linear interpolation between knots at each year end, sine and cosine terms for the annual cycle and for the bi-annual cycle. These models were compared using Akaike's Information Criterion. For sulphonamide compounds, chloramphenicol and MDR the best model included terms for both annual and bi-annual cycles. For streptomycin the best model included these terms and a continuous piecewise-linear trend in each year (this is similar to allowing for piecewise-constant year effects but removes the artificial jump that occurs in such models between 31 December of one year and 1 January of the next).

A GLSM was fitted for each individual antimicrobial (and MDR), and it is desirable that the spatial fields which are estimated for each antimicrobial be directly comparable so that we are able to determine any similarities. We therefore require the same temporal covariates for each antimicrobial, and so choose the covariate model with both annual and bi-annual cycles, and with continuous piecewise-linear trends in each year. The covariate terms in our model for the log-odds of resistance are therefore

with a different set of βs for each antimicrobial (and MDR). Here S03, S04, S05, and S06 allow for the piecewise linear trend and are maximal at the knot points of 1 January 2003, 2004, 2005 and 2006, respectively.

Fitting the GLSMs

For each antimicrobial (and MDR) the MCMC algorithm was run for 20 million iterations, of which the first 0·5 million were treated as burn-in; to keep file sizes manageable, only one in every 800 iterations was actually stored. For all antimicrobials the chain mixed thoroughly for the temporal covariate parameters β2, and for the scale parameter φ. Mixing for the variance σ2 and for the more important quantity w2/φ was adequate for all four MCMC runs.

Table 2 shows the parameter estimates (posterior medians) together with a 95% credibility interval (CI) for each of the four model fits.

Table 2. Point estimates (median) and 95% credibility interval

Credibility intervals marked with an asterisk

(*) do not contain zero and therefore show either a high probability that the parameter is positive or a high probability that the parameter is negative.

In terms of temporal effects, the most important terms (95% CI does not include zero) for streptomycin appear to be a bi-annual cycle and a peak in probability of resistance in the winter of 2004/2005. For sulphonamide compounds the most important terms involve the annual cycle, although the CIs for the coefficients of the bi-annual cycle only just include zero. The bi-annual cycle is also important for chloramphenicol, and one of the annual cycle terms only just includes zero. Finally the temporal MDR signal also appears to be mainly composed of an annual and bi-annual cycle. All of these results are consistent with the findings from the earlier exploratory GLM fits.

To visualize the temporal signals for each of the four fits we chose an ‘average’ point where the spatial signal S=0. For each iteration of the (thinned) MCMC sample we then calculated the probability that an incident of S. Typhimurium would be resistant to antimicrobial ‘A’ (or MDR). Figure 3 shows the median predicted value together with the 2·5th and 97·5th percentiles for each of the four model fits. These values are from the posterior predictive distribution of the temporal signal; the CI is not directly analogous to a confidence interval.

Fig. 3. Temporal signal for the probability of resistance for each of the four antimicrobials. Graphs show the median together with the 2·5% and 97·5% quantiles of N≈25 000 Markov Chain Monte Carlo samples. The probabilities are calculated using the average value of the spatial signal, i.e. with S(x)=0.

All four of the profiles showed a peak in the probability of resistance in mid to late spring, with a smaller peak in late autumn. The profiles for sulphonamide compounds, chloramphenicol, and MDR are very similar, all showing a consistent median level throughout the study period. The median level for streptomycin mirrored that of the other antimicrobials (and MDR) for the first 18 months before rising to a slight peak and then dropping off sharply over the final 18 months. This is consistent with Table 1 which shows that the proportion of positive results for streptomycin decreased to 49% in 2006 compared to 66% in 2003.

It is clear from Table 2 that there is a great deal of uncertainty in both spatial and temporal parameters. This is unsurprising given the paucity of data (294 binary observations), and indicates that any observed patterns should be treated with some caution, as the CIs will be large.

Figure 4 shows the posterior median estimate of the kriging surface S, the spatial contribution to the log-odds of an incident being resistant to a given antimicrobial. Each MCMC iteration supplies realizations from the joint posterior distribution of σ2 and φ, and the value of the surface S at each data of the 294 data points. Values for S on a fine grid over England and Wales can then be estimated via ordinary kriging [Reference Diggle and Ribeiro25].

Fig. 4. Posterior median of the kriging surface S, the spatial contribution to the log-odds of an incident being resistant to a given antimicrobial, from 25 000 Markov Chain Monte Carlo samples.

Patterns for the three antimicrobials and MDR show similarities, with higher odds in much of Wales, Wiltshire, and much of Devon, and with a further small peak around Leicestershire. All but the chloramphenicol plot also show an increased probability of resistance in North West England.

DISCUSSION

We analysed incidents of S. Typhimurium in cattle in England and Wales between 2003 and 2006. We looked for patterns in the spatial and temporal variability in the risk of antimicrobial resistance by fitting a GLSM to each of the four datasets. We focused on three antimicrobials (streptomycin, sulphonamide compounds, chloramphenicol) and MDR.

All four of the temporal profiles showed peak probability of resistance in mid-late spring and a lesser peak in late autumn. The mean signals for chloramphenicol, sulphonamide compounds, and MDR varied little from year to year, whereas that for streptomycin showed a sharp drop over the last 18 months. We are unsure why a sharp drop occurred in streptomycin as there were no obvious changes in farm management or laboratory procedures during 2005. This drop is surprising since streptomycin is linked to sulphonamide compounds in mode of resistance. One possibility is that it could have resulted from temporal changes in the dominant serotypes, which vary in their level of resistance [Reference Pangloli17, Reference Michel26]. Indeed, variation in resistance prevalence has been related to the clonal spread of particular strains as a result of husbandry and animal movement factors as well as to the variation in selective pressure exerted by antimicrobial usage [5].

It is possible that the levels of resistance that we recorded may be linked to antimicrobial usage since usage gives rise to selection pressure for resistance [11]. Sulphonamides and streptomycin are commonly used for treatment or as a prophylactic for respiratory and gastrointestinal infections in cattle [Reference Bateman, Prescott, Baggot and Walker27] and the seasonal patterns that we observed may be associated with seasonal patterns of cattle management involving antimicrobial use.

We hypothesize a connection between the peaks in antimicrobial resistance in spring and autumn that we recorded and the strong seasonal patterns in cattle births. Cattle births are characterized by a large spring peak and a smaller autumn peak [Reference Robinson and Christley28]. Seasonal patterns occur in both dairy and beef cattle; beef cattle calving peaks in spring (with the largest monthly births in April) and troughs in winter (December), while dairy cattle calving peaks in autumn (August or September) and troughs around May [29]. At these times antimicrobials (e.g. sulphonamides) are administered (when calves are removed from the dam), to treat or prevent diarrhoea and pneumonia [Reference McEwen and Fedorka-Cray30]. It has previously been noted that the selection pressure for antimicrobial resistance in S. Typhimurium is highest in calves due to the extensive use of antimicrobials in calf rearing and the type of Salmonella infection [Reference Prescott, Prescott, Baggot and Walker14]. Outbreaks of salmonellosis in calves in Great Britain have often been caused by phage types that are characteristically multiresistant, e.g. definitive phage types (DTs) 29, 204 and 104 and have been linked to therapeutic antimicrobial use. They have also caused serious illness in the human population [Reference Prescott, Prescott, Baggot and Walker14].

Seasonal patterns of adult cattle management involving antimicrobial usage could also play a part in seasonality of antimicrobial resistance. Antimicrobial usage and therefore selection pressure, tends to be high during the non-lactating phase for dairy cows and during cattle movement, and both of these activities peak in spring and autumn. Antimicrobials are routinely administered to entire adult cattle herds to prevent mastitis during the non-lactating period [Reference McEwen and Fedorka-Cray30], which is commonly in late spring (2 months prior to calving). Prophylactic treatments are also typically used during transport; a high-risk period for infectious disease [Reference McEwen and Fedorka-Cray30], which peaks strongly in spring and autumn (the autumn peak being the larger) [Reference Robinson and Christley28]. Other contributors to the increase in risk of antimicrobial resistance during transport include the movement of carrier animals between herds and the assembly of susceptible animals in close confinement [Reference McEwen and Fedorka-Cray30].

If the seasonal patterns described were the result of antimicrobial use, then this suggests a fairly rapid decrease in resistance at times when antimicrobial use in cattle is reduced. We found little evidence of rapid seasonal changes in antimicrobial resistance in the literature; however, quick response to the removal of antimicrobials is possible. There was a decrease in resistant bacteria in healthy animals and humans in the years immediately following the ban on animal antimicrobials as growth promoters in the European Union (EU) [Reference Wegener31]. In Enterococcus faecium isolates from chickens and pigs, for example, resistance prevalence declined within 1 year of removal of growth promoters in many European countries. Indeed declines were seen within the first 3 months for four different antimicrobials. It was suggested that in the absence of selection pressure, a susceptible population began to replace phenotypically resistant strains [Reference Bywater32].

Clearly antimicrobial consumption cannot explain the levels of resistance of chloramphenicol that we recorded because this substance has been banned from use in food-producing animals in the EU for many years. Co-resistance with other compounds is a likely explanation [Reference Bywater33].

A spatial signal was visible for each antimicrobial, with certain attributes, such as high probabilities of antimicrobial resistance in North Wales shared between all four datasets.

Figure 1 demonstrated that North Wales has a very high density of farms with incidents of S. Typhimurium and it is possible that these two facts are related. Indeed some of the highest densities of cattle farms are located in Wales and the west of England. Beef cattle tend to be concentrated in the east and southwest of Wales and South West England, while dairy farms tend to be in southwest Wales, central England and the west coast, extending further north than beef farms [34]. Detailed information about the location of all cattle holdings was not available during our study. However, more recently, the density of registered premises in Great Britain by species and the number of incidents of Salmonella and have been reported [35] and confirm that S. Typhimurium aggregates around the Welsh borders and South West England. Despite the consistency of the spatial patterns across our four datasets, the credibility intervals (not shown) are large and so the trends should be viewed as tentative at best.

The seasonal patterns of resistance were similar for all antimicrobials. This is unsurprising since all three antimicrobials often occur in the characteristic multiresistant pentavalent resistance pattern (resistance to tetracycline, ampicillin, chloramphenicol, streptomycin and sulphonamide compounds) because they are chromosomally integrated as a single genetic island [Reference Davison23]. In brief, resistance can be exchanged between different bacteria by mobile genetic elements including plasmids, transposons, integrons and bacteriophages that carry genes for antimicrobial resistance. In some cases these elements may accumulate or co-integrate in a single host bacterium to give MDR [Reference Baker-Austin36]. MDR can then be transferred through one coherent piece of DNA (plasmid) that encodes specific resistance genes [Reference Randall37]. The pentavalent pattern has often been recorded in Salmonella isolates from dairy farms in England and Wales [Reference Davison23] and is common in S. Typhimurium, particularly definitive phage type (DT)104 [Reference Randall37]. This multiresistant phage type has caused numerous infections in food animals and humans in the UK and worldwide since the 1990s [Reference Threlfall15]. Not all resistances are carried on mobile genetic elements, as described above. In future, if sufficient data become available, it would be useful to compare other antimicrobials that are not linked in their mode of resistance, e.g. nalidixic acid. This analysis might demonstrate a different antimicrobial resistance response to local risk factors such as changes in antimicrobial usage.

In this work we do not specify any risk factor other than temporal covariates while studying patterns of antimicrobial resistance in S. Typhimurium. However, we can suggest potential risk factors for antimicrobial resistance, which have previously been identified for Salmonella infection that may act locally or regionally. Local risk factors include feed suppliers [Reference Davis38], herd size or composition (large farms and high stocking densities facilitate Salmonella dissemination) [Reference Warnick16, Reference McEwen and Fedorka-Cray30, Reference Davison39] and environmental contamination, e.g. bedding and water [Reference Pangloli17]. An increase in temperature in summer can increase the range of potential sources of contamination [Reference Pangloli17]. Herd-associated methods to reduce disease include purchasing replacement stock from direct sources rather than dealers, quarantine of purchased cattle for 4 weeks, housing sick animals in isolation areas and preventing wild birds access to cattle feed stores [Reference Evans40]. On a larger scale, studies have demonstrated contagious processes through farm-to-farm transmission. Risk factors acting on a regional scale include movement of contaminated humans or equipment between farms [Reference McEwen and Fedorka-Cray30] the presence of suitable habitat for wildlife vectors [Reference Warnick16, Reference Wilson41] or environmental factors such as runoff from pastures and wastewater contaminating local water and spreading between farms [Reference Murray, Wray and Wray20]. Infection can also be influenced by spatiotemporal factors. Infection in dairy herds, for example, is not constant over time and farms are more likely to become Salmonella positive if there are more positive farms within 30 km [Reference Fenton21]. A more thorough (and necessarily more complex) analysis might try to gauge the extent of these effects by allowing for such factors explicitly. To begin with we suggest stratification of data, when sufficient is available, to assess the differences between dairy and beef herds. We also suggest incorporating trade relations (e.g. movement of cattle between farms) in the analysis. The vast majority of cattle movements occur over <50 km [Reference Robinson and Christley42]; incorporating network parameters that account for the rapidly changing structure of the livestock industry would be very informative. Ultimately this could lead to the development of targeted surveillance and prevention measures aimed at reducing the incidence of antimicrobial resistance.

ACKNOWLEDGEMENTS

The authors thank Robert Smeatham and Christina Papadopoulou at the VLA for help with the Farmfile database, and Nicola Williams of the National Centre for Zoonoses, University of Liverpool for helpful discussions. Part of this work was funded through North West Development Agency project no. N0003212.

DECLARATION OF INTEREST

None.

References

REFERENCES

1. Department for Environment Food and Rural Affairs. Salmonellosis (Salmonella species). In: Zoonoses Report UK 2009 . London, 2009, pp. 121.Google Scholar
2. Health Protection Agency. Salmonella (http://www.hpa.org.uk/Topics/InfectiousDiseases/InfectionsAZ/Salmonella/). Accessed June 2011.Google Scholar
3. Threlfall, EJ. Antimicrobial drug resistance in Salmonella: problems and perspectives in food- and water-borne infections. Microbiology Reviews 2002; 26: 141148.Google Scholar
4. Threlfall, EJ, et al. A comparison of antimicrobial susceptibilities in nontyphoidal Salmonellas from humans and food animals in England and Wales in 2000. Microbial Drug Resistance 2003; 9: 183189.Google Scholar
5. Department for Environment Food and Rural Affairs. Salmonella in livestock production in Great Britain. Weybridge: Veterinary Laboratories Agency, 2005, pp. 1168.Google Scholar
6. Meakins, A, et al. Antimicrobial drug resistance in human nontyphoidal Salmonella isolates in Europe 2000–2004: a report from the Enter-net international surveillance network. Microbial Drug Resistance 2008; 14: 3135.Google Scholar
7. Miller, GY, McNamara, PE, Singer, RS. Stakeholder position paper: Economist's perspectives on antibiotic use in animals. Preventive Veterinary Medicine 2006; 73: 163168.Google Scholar
8. Smith, BP. Salmonellosis in ruminants. In: Smith, BP, ed. Large Animal Internal Medicine , 4th edn. St Louis: Mosby, 2008, pp. 877880.Google Scholar
9. Molbak, K, et al. An outbreak of multidrug-resistant, quinolone-resistant Salmonella enterica serotype Typhimurium DT104. New England Journal of Medicine 1999; 341: 14201425.Google Scholar
10. Frenzen, PD, et al. Salmonella cost estimate updated using FoodNet data. Food Review 1999; 22: 1015.Google Scholar
11. World Health Organization. Joint FAO/OIE/WHO expert workshop on non-human antimicrobial usage and antimicrobial resistance: scientific assessment. Geneva, Switzerland, 2004, pp. 15.Google Scholar
12. Varma, JK, et al. Hospitalization and antimicrobial resistance in Salmonella outbreaks, 1984–2002. Emerging Infectious Diseases 2005; 11: 943946.Google Scholar
13. Helms, M, et al. Excess mortality associated with antimicrobial drug-resistant Salmonella Typhimurium. Emerging Infectious Diseases 2002; 8: 490495.Google Scholar
14. Prescott, JF. Antimicrobial drug resistance and its epidemiology. In: Prescott, JF, Baggot, JD, Walker, RD, eds. Antimicrobial Therapy in Veterinary Medicine . Ames, IA, USA: Iowa State University Press, 2000, pp. 2749.Google Scholar
15. Threlfall, EJ, et al. The emergence and spread of antibiotic resistance in food-borne bacteria. International Journal of Food Microbiology 2000; 62: 15.Google Scholar
16. Warnick, LD, et al. Risk factors for clinical salmonellosis in Virginia, USA cattle herds. Preventive Veterinary Medicine 2001; 49: 259275.Google Scholar
17. Pangloli, P, et al. Seasonal incidence and molecular characterization of Salmonella from dairy cows, calves, and farm environment. Foodborne Pathogens and Disease 2008; 5: 8796.Google Scholar
18. Lanzas, C, et al. Transmission dynamics of a multidrug-resistant Salmonella Typhimurium outbreak in a dairy farm. Foodborne Pathogens and Disease 2010; 7: 467474.Google Scholar
19. Diggle, PJ. Statistical Analysis of Spatial Point Patterns, 2nd edn. Oxford, UK: Hodder Arnold, 2003.Google Scholar
20. Murray, CJ. Environmental aspects of Salmonella. In: Wray, C, Wray, A, eds. Salmonella in Domestic Animals . Wallingford: CABI Publishing, 2000, pp. 265284.Google Scholar
21. Fenton, S, et al. Spatial and spatio-temporal analysis of Salmonella infection in dairy herds in England and Wales. Epidemiology and Infection 2009; 137: 847857.Google Scholar
22. Gibbens, JC, et al. Use of laboratory data to reduce the time taken to detect new diseases: VIDA to FarmFile. Veterinary Record 2008; 162: 771776.Google Scholar
23. Davison, HC, et al. Prevalence, incidence and geographical distribution of serovars of Salmonella on dairy farms in England and Wales. Veterinary Record 2005; 157: 703711.Google Scholar
24. Hastie, TJ, Tibshirani, RJ. Generalized Additive Models: Chapman and Hall/CRC (Monographs on Statistics and Applied Probability), 1987, pp. 352.Google Scholar
25. Diggle, PJ, Ribeiro, PJ Jr.. Model-based Geostatistics. New York: Springer, 2007.Google Scholar
26. Michel, P, et al. Regional, seasonal, and antimicrobial resistance distribution of Salmonella Typhimurium in Canada. Canadian Journal of Public Health 2006; 97: 470474.Google Scholar
27. Bateman, KG. Antimicrobial drug use in cattle. In: Prescott, JF, Baggot, JD, Walker, RD, eds. Antimicrobial Therapy in Veterinary Medicine . Ames, IA, USA: Iowa State University Press, 2000, pp. 576590.Google Scholar
28. Robinson, SE, Christley, RM. Identifying temporal variation in reported births, deaths and movements of cattle in Britain. BMC Veterinary Research 2006; 2.Google Scholar
29. Department for Environment Food and Rural Affairs. The cattle book 2008. Descriptive statistics of cattle numbers in Great Britain on 1 June 2008. Defra, 2008.Google Scholar
30. McEwen, SA, Fedorka-Cray, PJ. Antimicrobial use and resistance in animals. Clinical Infectious Diseases 2002; 34: S93–106.Google Scholar
31. Wegener, HC. Antibiotics in animal feed and their role in resistance development. Current Opinion in Microbiology 2003; 6: 439445.Google Scholar
32. Bywater, R, et al. The susceptibility to growth-promoting antibiotics of Entercoccus faecium isolates from pigs and chickens in Europe. Journal of Antimicrobial Chemotherapy 2005; 56: 11391143.Google Scholar
33. Bywater, R, et al. A European survey of antimicrobial susceptibility among zoonotic and commensal bacteria isolated from food-producing animals. Journal of Antimicrobial Chemotherapy 2004; 54: 744754.Google Scholar
34. Department for Environment Food and Rural Affairs. Land use and livestock (http://www.defra.gov.uk/evidence/statistics/foodfarm/landuselivestock/index.htm). Accessed 18 March 2010.Google Scholar
35. Department for Environment Food and Rural Affairs. Salmonella in livestock production in Great Britain. Weybridge: Veterinary Laboratories Agency, 2008, pp. 1168.Google Scholar
36. Baker-Austin, C, et al. Co-selection of antibiotic and metal resistance. Trends in Microbiology 2006; 14: 176182.Google Scholar
37. Randall, LP, et al. Antibiotic resistance genes, integrons and multiple antibiotic resistance in thirty-five serotypes of Salmonella enterica isolated from humans and animals in the UK. Journal of Antimicrobial Chemotherapy 2004; 53: 208216.Google Scholar
38. Davis, M, et al. Feedstuffs as a vehicle of cattle exposure to Escherichia coli O157:H7 and Salmonella enterica . Veterinary Microbiology 2003; 95: 199210.Google Scholar
39. Davison, HC, et al. Risk factors associated with the Salmonella status of dairy farms in England and Wales. Veterinary Record 2006; 159: 871880.Google Scholar
40. Evans, SJ. A case control study of multiple resistant Salmonella Typhimurium DT104 infection of cattle in Great Britain. Cattle Practice 1996; 4: 259266.Google Scholar
41. Wilson, SJ, et al. Nontyphoidal Salmonellae in United Kingdom Badgers: Prevalence and spatial distribution. Applied and Environmental Microbiology 2003; 69: 43124315.Google Scholar
42. Robinson, SE, Christley, RM. Exploring the role of auction markets in cattle movements within Great Britain. Preventive Veterinary Medicine 2007; 81: 2137.Google Scholar
Figure 0

Fig. 1. Locations of the 256 farms which experienced at least one incident of S. Typhimurium between January 2003 and December 2006.

Figure 1

Table 1. Number of Salmonella Typhimurium incidents per year and frequency of antimicrobial resistance to each antimicrobial (only antimicrobials for which there was at least one incident of resistance are listed)

Figure 2

Fig. 2. Incidents of S. Typhimurium in cattle farms in England and Wales that were resistant to chloramphenicol between 1 January 2003 and 31 December 2006. Incidents with bacteria resistant to chloramphenicol are indicated by an open symbol (○), while incidents susceptible to chloramphenicol are indicated by a cross (×).

Figure 3

Table 2. Point estimates (median) and 95% credibility interval

Figure 4

Fig. 3. Temporal signal for the probability of resistance for each of the four antimicrobials. Graphs show the median together with the 2·5% and 97·5% quantiles of N≈25 000 Markov Chain Monte Carlo samples. The probabilities are calculated using the average value of the spatial signal, i.e. with S(x)=0.

Figure 5

Fig. 4. Posterior median of the kriging surface S, the spatial contribution to the log-odds of an incident being resistant to a given antimicrobial, from 25 000 Markov Chain Monte Carlo samples.