- Research article
- Open Access
Spatio-temporal modelling of verotoxigenic Escherichia coli O157 in cattle in Sweden: exploring options for control
- Stefan Widgren^{1, 2}Email authorView ORCID ID profile,
- Stefan Engblom^{2},
- Ulf Emanuelson^{3} and
- Ann Lindberg^{1}
- Received: 22 January 2018
- Accepted: 20 July 2018
- Published: 2 August 2018
Abstract
A spatial data-driven stochastic model was developed to explore the spread of verotoxigenic Escherichia coli O157 (VTEC O157) by livestock movements and local transmission among neighbouring holdings in the complete Swedish cattle population. Livestock data were incorporated to model the time-varying contact network between holdings and population demographics. Furthermore, meteorological data with the average temperature at the geographical location of each holding was used to incorporate season. The model was fitted against observed data and extensive numerical experiments were conducted to investigate the model’s response to control strategies aimed at reducing shedding and susceptibility, as well as interventions informed by network measures. The results showed that including local spread and season improved agreement with prevalence studies. Also, control strategies aimed at reducing the average shedding rate were more efficient in reducing the VTEC O157 prevalence than strategies based on network measures. The methodology presented in this study could provide a basis for developing disease surveillance on regional and national scales, where observed data are combined with readily available high-resolution data in simulations to get an overview of potential disease spread in unobserved regions.
Introduction
Verotoxigenic Escherichia coli O157:H7 (VTEC O157) is a major causative agent of human enterohemorrhagic E. coli (EHEC) cases, with cattle being an important reservoir for the bacteria [1]. It has been suggested that a substantial proportion of the EHEC cases could be prevented if VTEC O157 shedding from cattle could be controlled [2, 3]. However, interventions aimed at reducing the VTEC O157 prevalence on a regional and national scale have so far not been successful. The distribution of VTEC O157 in the Swedish cattle population is not uniform in space and time; the prevalence is higher in southern Sweden [4, 5], and more cattle herds are infected during the autumn [6]. Risk-based surveillance is an approach to cost-efficiently identify infected premises by directing the sampling to those herds which are more likely to be infected [7]. Information regarding relevant risk factors for the pathogen on a geographical, herd, and individual level are required for planning such activities. A better understanding of the mechanisms behind the spatial and seasonal patterns of VTEC O157 infections in cattle could improve surveillance and prevention strategies.
Mathematical modelling and computer simulations provide means by which various assumptions and simplifications of pathogen spread in a complex system may be explored. A recent simulation study of VTEC O157 transmission involving the complete Swedish cattle population, showed that the spatial pattern in prevalence may be due to regional differences in livestock movements [8]. However, the model seemingly overestimated the prevalence of VTEC O157 in the northern regions. There are biologically plausible factors that could be incorporated in the model to improve the modelling accuracy and thus also our understanding. For example, it has been suggested that proximity to VTEC O157 infected farms is a risk factor for presence of VTEC O157 on farms because of local spread [6, 9, 10]. Hence heterogeneities in holding density might explain the different spatial pattern in prevalence. It has further been suggested that the seasonal variation in ambient temperature could have a considerable effect on the survival and growth of VTEC O157 in the environment and consequently on the indirect transmission dynamics [11, 12]. The cattle herd density is higher [13] and the duration of the warmer period longer [14] in the southern part of Sweden. These two parameters can be incorporated in the model using livestock and meteorological data.
It is commonly observed that VTEC O157 infection on farms cease without intervention [6, 15]. This implies that mitigation measures should rather strive at reducing the probability of introduction and decreasing the likelihood of spread, for example by increasing biosecurity, limit movement of infected cattle and reducing transmission from infectious animals to the environment by influencing shedding. Interventions aimed at reducing the impact of high-shedding individuals has been shown, in a modelling framework, to efficiently reduce the VTEC O157 prevalence [16]. However, to identify high-shedding individuals by faecal sampling is a challenge; the shedding pattern is highly variable, and it is possible that high-shedding events appear in many cattle for brief periods [17, 18]. Alternatively, mitigations could target farms that constitute a high risk for transmitting the infection. For example, using information from pig movement data to target control was reported to efficiently reduce the disease spread in simulations [19]. Simulations investigating interventions based on cattle movement data e.g. targeting highly connected herds, could further inform efficient control strategies for VTEC O157.
A mathematical model of infection dynamics in spatially segregated herds can be approached using a metapopulation model, i.e., a collection of subpopulations (herds) with their own infection dynamics [20]. To capture infection dynamics on a regional and national scale, the subpopulations must include interactions among them to model between-herd spread, for example, from livestock movements. Additionally, since each herd size is small, it is generally necessary to use stochastic models and run many iterations, for example, to account for the random event that an infection will go extinct [21]. Thus, performance of the simulator becomes critical when using these methods for national-scale disease spread modelling. There exist several recent approaches in computational epidemiology to implement such models, for example, KENDRICK [22], EMuLSion [23] and SimInf [24]. In the present study we used the SimInf framework which integrates infection dynamics in subpopulations as continuous-time Markov chains (CTMC) using the Gillespie stochastic simulation algorithm [25] and is designed to efficiently incorporate population demographics and livestock movement data.
The aims of this study were to use mathematical modelling and simulations to explore: (i) the importance of local spread of VTEC O157 among proximal holdings, (ii) the seasonal effect on the spatial distribution of VTEC O157 in Sweden, and (iii) the effectiveness of various interventions to reduce the prevalence of VTEC O157 in the Swedish cattle population.
Materials and methods
Input data
The complete Swedish livestock data for 37 221 holdings during the period 1 July 2005 to 31 December 2013 was incorporated in the simulations to drive the population demographics and the time-varying contact network from livestock movements. The data originated from the Swedish national cattle database managed by the Swedish Board of Agriculture. The dataset, described in detail in [8], contained 18 649 921 events reported on an individual animal level. Briefly, the livestock data included the following information about each cow: (i) the date and the holding for its birth, (ii) the date and the source and destination holding for any movements, and (iii) the date for slaughter or death. Moreover, all holdings (n = 37 221) were geopositioned using the available information with the highest spatial resolution for each holding. Exact coordinates were available for 31 188 holdings. For other farms, the 5-digit postal code was known for 4748 holdings and the locations were randomly sampled within the postal code area using the R package sp [26]. Similarly, for 1283 holdings the 3-digit postal area code was known (contains multiple postal codes) and each location was randomly sampled within the postal area. Finally, for two holdings the postal area code was missing. For these, other complementary information was used to determine the postal area to which they belonged. Random sampling of their coordinates was then performed in the postal areas. Each holding should correspond to a single geographical location, e.g. a farm building or a pasture, where animals are kept. However, some coordinates (n = 337) were identical between holdings. It was assumed that these coordinates correspond to e.g. a farm building and a nearby pasture managed by the same farmer. The distance d_{ik} was calculated for every holding i to all neighbours k. A minimum distance of 100 m was used to separate any two holdings to avoid division by zero in Equation (2) below.
Data from the Swedish Meteorological and Hydrological Institute (SMHI) was used to determine the duration of the seasons (winter, spring, summer, and autumn) depending on geographical location. The seasons were defined by SMHI per the average temperature for the reference period 1961–1990: (winter) below 0 °C, (spring) between 0 and 10 °C, (summer) above 10 °C and (autumn) between 0 and 10 °C [14]. Determined by the location, each holding had its own definition of when (day of the year) each season started. Note that this approach does not incorporate a year-to-year variation in when the seasons begin.
The dataset that was used to calibrate the model parameters (described below), was from a longitudinal herd-level observational study conducted over 38 months (2009–2012) in 126 cattle herds in four distinct regions in southern Sweden [6]. In that study, the VTEC O157 herd status was repeatedly determined, on average at 17 occasions 64 days apart. VTEC O157 was detected in 67 of the herds on 224 occasions. The VTEC O157 herd status was determined from environmental sampling among calves (younger than 120 days) and young stock (between 120 and 365 days).
To evaluate the simulated spatial distribution of infected holdings, it was visually compared with the observed distribution of a nation-wide monitoring of VTEC O157 in Swedish cattle that was carried out in abattoirs from autumn 2011 to autumn 2012 [27]. In total, the pathogen was detected in 73 out of 2376 randomly collected faecal samples. Moreover, the farm that each sampled animal originated from was geopositioned (Figure 3).
Pathogen transmission model
The VTEC O157 transmission was modelled using an extended variant of the stochastic compartment model SIS_{E} described elsewhere [8, 28]. The model contained the compartments susceptible (S) and infected (I) as well as the environmental compartment (E), contaminated with free living pathogens by infected and shedding animals. The susceptible and infected compartments were further divided into three age categories to match the age categories in the longitudinal herd-level observational study [6] used for calibration. Let index j denote the age category where: (1) calves—younger than 120 days, (2) young stock—between 120 and 365 days, and (3) adults—older than 365 days.
To incorporate seasonality in the infection dynamics, β could vary during the year. However, compared to the model in [8], where the year was evenly divided into four quarters, it was instead divided into the four seasons: spring, summer, autumn and winter, with β_{1}, β_{2}, β_{3}, and β_{4} reflecting the rate of decay in each season, respectively. In this way, the parameter β could indirectly capture the effect of varying seasons.
Model calibration
Parameters in a stochastic SIS _{ E} VTEC O157 model
Parameter | Description (unit) | Value |
---|---|---|
Α | Rate of contribution to the environmental infectious pressure per infected individual (units per day) | ^{a}1.00 × 10^{0} |
β _{1} | Decay of environmental infectious pressure during spring (per day) | 1.57 × 10^{−1} |
β _{2} | Decay of environmental infectious pressure during summer (per day) | ^{a}1.44 × 10^{−1} |
β _{3} | Decay of environmental infectious pressure during autumn (per day) | 1.50 × 10^{−1} |
β _{4} | Decay of environmental infectious pressure during winter (per day) | 1.57 × 10^{−1} |
υ _{1} | Indirect transmission rate of the environmental infectious pressure in calves (per animal per day) | 2.48 × 10^{−2} |
υ _{2} | Indirect transmission rate of the environmental infectious pressure in young stock (per animal per day) | 2.48 × 10^{−2} |
υ _{3} | Indirect transmission rate of the environmental infectious pressure in adults (per animal per day) | 1.37 × 10^{−2} |
γ | The recovery rate of infection (per day) | ^{a}1.00 × 10^{−1} |
D | The rate of local spread among proximal holdings (per day per m) | 0.11 × 10^{−5} |
The longitudinal herd-level observational study [6] was repeatedly replicated in simulations to determine \(G_{1} \left(\varvec{\theta}\right)\) and \(G_{2} \left(\varvec{\theta}\right)\), as follows. Let \(Y_{in}^{*}\) denote the nth observed VTEC O157 status (1-positive; 0-negative) in holding i at time t_{n}. Similarly, let \(Y_{in} \left(\varvec{\theta}\right)\) denote the simulated status, corresponding to \(Y_{in}^{*}\). To determine the status \(Y_{in} \left(\varvec{\theta}\right)\) that could have been found if simulated holdings had been sampled, the sampling strategy from the observed data was replicated as previously described [8]. Briefly, the sampling was simulated at each sample point as follows. First, pools (pool size = 3) were randomly created within each age category from the number of susceptible and infected individuals at the time for the sample point in the simulation. Given the proportion of infected individuals in a pool, it was randomly classified as positive or negative, with P (positive) equal to the estimated test sensitivity [29]. Finally, using the estimated pool prevalence, the simulated herd status was randomly classified as positive or negative given the sensitivity of the sampling protocol [30].
The Nelder–Mead algorithm [31] in R [32] was used to find the parameter combination \(\varvec{\theta}\) that minimised the objective function \(G\left(\varvec{\theta}\right)\) under the constraint that all parameters \(\varvec{\theta}\ge 0\). Exploratory analysis indicated that fitting a model with the parameters υ_{1}, υ_{2}, υ_{3}, β_{1}, β_{2}, β_{3}, β_{4}, and D, was not feasible i.e. the system was not identifiable. Therefore, to maintain model parsimony, υ_{1} and υ_{2} (calves and young stock indirect transmission rates) were fitted as a single parameter, and β_{2} (summer decay) was fixed at a decimal reduction rate (the time required at a given temperature to kill 90% of VTEC O157) of 16 days [33–35]. Moreover, the average recovery rate was fixed at γ = 0.1 [36]. Finally, the contamination rate α, was fixed at 1.0 per day, thus defining the unit of the environmental infectious pressure variable \({\varphi_i}(t)\). The disease spread simulations were performed with the SimInf package in R [8, 21, 24, 37].
The Nelder–Mead algorithm was started from two different parameter conditions and then restarted a third time from the parameters found for the minimum \(G\left(\varvec{\theta}\right)\). Each Nelder–Mead optimisation ran for 250 iterations. To quantify stochasticity in \(G\left(\varvec{\theta}\right)\), the coefficient of variation was calculated from the last 10 iterations of the Nelder–Mead.
Initialisation
Simulations started 1 July 2005 with a uniform geographical distribution with 25% of the holdings randomly sampled to be infected. To have an overall individual prevalence of 4% [5], the within-holding prevalence in initially infected holdings was set to 16% in all age categories. Moreover, the initial environmental infectious pressure \({\varphi_i}\) was set to zero in each holding. The model was simulated for a burn-in of 1570 days before the parameter estimation started in October 2009.
Evaluation of the calibrated model
Using the calibrated model parameters, the model was graphically evaluated based on the outcome of 1000 trajectories. First, the distributions of the simulated data \(\overline{a}_{q}\) and \(\overline{b}_{q}\) in each quarter were plotted in a boxplot and compared with the observed data \(a_{q}^{*}\) and \(b_{q}^{*}\). Secondly, distributions of weekly time-series of the herd-level prevalence and the individual-level prevalences in each age category were generated for the period 1 January 2008 to 31 December 2013. Finally, the spatial distribution of the proportion of trajectories where each holding had at least one infected animal was visualised for the following dates: 1 October 2011, 1 January 2012, 1 April 2012, and 1 July 2012. The spatial distribution was compared with the outcome from the nation-wide monitoring of VTEC O157 in Swedish cattle 2011–2012 [27].
Sensitivity analysis
Sensitivity analysis was performed to explore how perturbations of the six calibrated model parameters (υ_{1,2}, υ_{3}, β_{1}, β_{3}, β_{4}, D) would influence \(G\left(\varvec{\theta}\right)\). First, each parameter was perturbed at multiple values, keeping the others fixed at the calibrated value. For every perturbation, the average \(G\left(\varvec{\theta}\right)\) was estimated from 40 trajectories. Secondly, a global sensitivity analysis was conducted since there may be interactions among parameters, thus affecting \(G\left(\varvec{\theta}\right)\). The Extended Fourier Amplitude Sampling Test (eFAST) [38, 39] was used for the global sensitivity analysis. This measure gives an indication of the influence of each parameter and groups of parameters on \(G\left(\varvec{\theta}\right)\). Based on the exploratory perturbation analysis above, the parameters were constrained to the following respective range: υ_{1,2} (0.022–0.027), υ_{3} (0.012–0.015), β_{1,3,4} (0.13–0.2), and D (0–0.3). Additionally, a “dummy” (1–10) parameter that does not affect the simulation was added for statistical comparison. The eFAST analysis was performed using 65 samples per search curve and a resampling size of 3 and 40 replicates per parameter set. The spartan R package [40], was used to generate the 1365 parameter value sets and for analysing the response in \(G\left(\varvec{\theta}\right)\). The first-order (Si) and total-order (STi) sensitivity indexes were calculated for each parameter and plotted. A two-sample t-test (significance level 0.95) was estimated to indicate significance of each parameters sensitivity index, contrasted to the “dummy” parameter.
Implementation of control strategies
The effectiveness of control strategies was investigated from numerical experiments comparing a baseline, i.e. the outcome from simulations with the calibrated model parameters, with the outcome from simulations with an adjusted model that reflected the control scenario in question. One thousand trajectories were generated for each investigated control scenario, and weekly time series of the average herd-level prevalence were generated for the period 1 January 2008 to 31 December 2013. All control scenarios started 1 January 2009.
First, control scenarios aimed at reducing the between-herd transmission were investigated. One approach to limit connectivity between herds is simply to remove some livestock movements from the data. However, removing the movements will affect the demography. Therefore, infected animals about to be moved, were replaced with susceptible animals, leading to the same results for the risk of infection, but keeping the demography in adequacy with the observed data. Two network-based strategies were explored in the simulations to inform when to apply control. One strategy was to reduce the number of holdings from which an infection could be introduced. This was done by considering a threshold for applying control based on the in-degree (ID) [41] and the ingoing contact chain (ICC) [42] of the receiving holding. The ID of a holding is defined as the number of holdings that move animals directly to the holding for a defined period. The ICC metric additionally includes all holdings that have indirect movements i.e. all holdings that can be reached when tracing the sequence of movements over time. Similarly, the other strategy was to reduce the number of holdings to which an infection could spread. This was done by instead applying control based on the out-degree (OD) [41] and the outgoing contact chain (OCC) [43] of the sending holding. The OD of a holding is defined as the number of holdings that the holding moves animals to for a defined period. The OCC metric additionally includes all holdings that can be reach indirectly when tracing the sequence of movements over time. The ID, OD, ICC and OCC over 90 days were calculated weekly for each holding during the period 1 July 2005 to 31 December 2013 using the R package EpiContactTrace [44]. After inspecting the 80–97.5% percentiles of ID, OD, ICC and OCC summarised for the whole period and over all holdings, the following thresholds were considered for ID and OD: >1, >2, >3 and >4, and for ICC and OCC: >2, >4, >6, >8 and >10. In the simulations, the model simply changed the state of an infected individual to susceptible when it was moved, conditional on the weekly network measures, so that infection was not transferred to the destination holding [37].
Next, control scenarios aimed at reducing the animal-to-animal transmission within a herd were explored. A reduction of the between-animal transmission can be accomplished by making cattle less likely to become infected and excrete the pathogen. One approach to achieve this is to use vaccination [reviewed by 45, 46], which has been shown to inhibit VTEC O157 colonisation of the terminal recta and thereby leading to a significant reduction in shedding. Another approach is to adjust the microflora of the gastrointestinal tract to prevent colonisation of VTEC O157 by feeding probiotics [reviewed by 47]. Although the complex mechanisms involved in treatments, such as vaccination or probiotics, are not specifically included in the model, it was assumed that modifying the parameters α (rate of contamination to the environment) and υ (indirect transmission rate) would mimic the effect of such a treatment. The following parameter adjustments were considered to explore effects of reduced animal-to-animal transmission: (i) 10% reduction of the contamination rate α, (ii) 10% reduction of the indirect transmission rate υ, and (iii) both (i) and (ii) together.
Finally, although unrealistic to consider for control, scenarios were also generated to compare the outcome when blocking all between-holding transmission routes via: (i) livestock movements, (ii) the spatial coupling among proximal holdings, and (iii) both (i) and (ii).
Results
Model calibration
The model reached a minimum at \(G\left(\varvec{\theta}\right) = 16.75\), where \(G_{1} \left(\varvec{\theta}\right) = 6.42\) and \(G_{2} \left(\varvec{\theta}\right) = 10.33\) for the parameter values of \(\varvec{\theta}\) shown in Table 1. The coefficient of variation from 10 replicates of \(G\left(\varvec{\theta}\right)\) was estimated to 2.85%. The indirect transmission rate was higher for animals under 1 year of age compared to older animals in the calibrated model. Furthermore, the highest rate of the bacterial reduction per day was during the winter and spring seasons and the least reduction was during the summer.
Sensitivity analysis
Implementation of control strategies
Discussion
The present simulation study assumed that between-herd disease transmission of VTEC O157 occurs by moving infected animals and that there also exists spatial spread among proximal herds. In addition, it was assumed that season affects the survival of VTEC O157 and thus infection dynamics. The results show that the location of geographical regions with the most pronounced clustering of infected holdings was in agreement with previously reported prevalence studies [4, 5, 48]. The clusters at these locations are likely caused by a complex interaction between the season [14], the local cattle holding density in these regions [13] and that most of the animal movements are over short to medium distances (i.e. <100 km) [49]. The simulated holding-level prevalence was also consistent with the 8.9% reported in a nationwide study of Swedish dairy herds [48]. The results are also in agreement with the individual-level prevalence reported in another previous nationwide study [5], although they are not directly comparable since the age categories are different. In that study, the observed prevalence was 3.5% in young stock (12–24 months) and 1.5% in adults (>36 months). In comparison, the simulated prevalence was 4.6–7.3% in young stock (3–12 months) and 1.2–2.1% in adults (>12 months). These findings suggest that the proposed transmission model captures important aspects of the VTEC O157 spatio-temporal dynamics on a national scale.
Incorporating both season, defined as the average temperature at the location of the holding, and spatial spread among neighbouring holdings, improved the agreement with observations compared to the model proposed by Widgren et al. [8] which lacked these characteristics. The results show that VTEC O157 transmission occurs at a much lower rate in northern Sweden, and that if infection was introduced, it became extinct over time. This is consistent with previous observations, as the VTEC O157 bacteria has rarely been detected in cattle originating from northern Sweden in prevalence studies [4, 5]. The environmental infectious pressure of VTEC O157 had the slowest decay during the summer, which leads to higher levels of the pathogen during that season. This finding is in agreement with the hypothesis that increased summertime VTEC O157 colonisation results from increased seasonal oral exposure to the pathogen [12].
Among the investigated control strategies, the one that most efficiently reduced the holding-level prevalence, was the combination of reducing both the indirect transmission rate and the rate of contribution to the environmental infectious pressure, per infected individual. These findings are consistent with previous modelling studies, where reducing the level at which infected cattle shed the pathogen [16, 50, 51], and decreasing the indirect transmission rate [52] has been shown to efficiently reduce the VTEC O157 prevalence in cattle. In the present study, a spatial coupling among proximal herds was included in the model to capture between-herd transmission of VTEC O157 unrelated to cattle movements, for example, via fomites such as vehicles or equipment, or by personnel. Although it is unrealistic to completely remove the spatial coupling among proximal herds, the results show a substantial reduction in prevalence after removing this transmission route, which highlights the importance of good external biosecurity to prevent spread and introduction of pathogens.
The results also show that the prevalence was, in principle, unaffected by network interventions, targeted on the measures in-degree, out-degree and ingoing contact chain as well as outgoing contact chain. This is in agreement with work by Zhang and Woolhouse [51] who reported that reducing movement related transmission has, at best, a modest impact in reducing the steady-state prevalence of VTEC O157. One approach to explore interventions involving animal trade is to generate a new contact network that incorporates the intervention of interest and then study the disease transmission in the new network [53]. However, a drawback of that approach is that it requires rules for how to rewire trading partners in the network; and those rules could disrupt other inherent properties of the livestock data that are important for the transmission process, for example, the population demographics. In this study, an alternative approach was used, which allowed the same network data to be used in all simulations. Instead of removing some movements from the network data to reduce between-herd spread, all movements were processed, but any infected cattle were replaced with susceptible cattle before adding them to the destination herd. Even this approach was not sufficient to efficiently reduce VTEC O157 within the Swedish cattle population. This also indicates a higher contribution to the between-herd spread from local spatial transmission compared to livestock movements.
The model was unable to completely capture the quarter-to-quarter fluctuations in the longitudinal observational study [6] used for calibration. This outcome might be explained by how the seasons were defined in this study, where the average temperature for the reference period 1961–1990 was used to classify seasons and not the actual temperature data for each year. This approach was used since reference data was readily available [14]. Although daily weather data was available, implementing a model informed by this detailed information at the herd-level would be computationally costly. Because precipitation and relative humidity might also affect conditions for bacterial replication and decay and have also been associated with increased VTEC O157 shedding [54], we suggest future research to explore the effect on the infection dynamics after incorporating high resolution spatio-temporal meteorological data in simulations. This is also supported by the global sensitivity analysis which suggests the seasonal parameters have a significant effect on the simulation behaviour. Another limitation of the disease spread model presented in this paper was that herd type was not included. It has been reported from several field studies that the production type e.g. dairy or beef, influences the risk for presence of VTEC O157 [6, 55, 56]. However, herd type was not available in the data that were used for modelling the population demographic and the temporal network.
The sensitivity analysis showed that the first order effect of the spatial coupling was not statistically significant for the simulation behaviour. A similar finding was reported by Zhang et al. [57] from fitting stochastic models for spread of VTEC O157 infection among Scottish cattle farms. However, in the present study, the sensitivity analysis showed that non-linear effects between the spatial coupling and other model parameters were statistically significant, suggesting that the spatial coupling should be included in the proposed VTEC O157 spread model. One limitation in the analysis of the spatial coupling is that exact coordinates were only available for 83.8% of the herds in the register data. This highlights the importance of a complete and validated cattle register to conduct data-driven disease spread simulations to evaluate potential risk factors and explore control strategies.
Using stochastic disease spread models is helpful to increase our understanding of the complex infection dynamics among interconnected herds. However, stochastic models introduce several challenges for inference with respect to model specification, parameter calibration, algorithm complexity and computational time. In this study, the parameter calibration was formulated as an optimisation problem using the Nelder–Mead algorithm to determine point estimates. However, since the model was stochastic, many realisations of the model had to be simulated for every iteration of the Nelder–Mead algorithm, and this increased the computational time. This suggests that parameterisation might require more sophisticated approaches e.g. to use a Bayesian methodology with suitable prior information on the parameter values. Another advantage of using Bayesian methods would be to estimate confidence intervals for the parameter values. The results presented in this study are based on fitting a data-driven pathogen-transmission model to a dataset consisting of repeated environmental sampling in 126 holdings, all located in southern Sweden. Due to the relatively limited size of that dataset, careful judgements must be considered when generalising to a national scale. However, the fact that the proposed methodology seems to capture the observed prevalence and infection dynamics of VTEC O157 on a national and regional scale, warrants further studies on other hazards of importance to animal and public health, e.g. Salmonella or antimicrobial resistance, using similar methodology.
It has been suggested that a substantial proportion of the VTEC O157 human cases could be prevented by vaccinating cattle against the pathogen [3]. The results of this study showed that reducing both the indirect transmission rate and the rate of contribution to the environmental infectious pressure, per infected individual, efficiently decreased the prevalence of VTEC O157 in the Swedish cattle population. This supports that vaccination could be a viable option to control VTEC O157 in cattle and hence have a public health benefit.
Declarations
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors designed the simulation study. SW performed the simulations and drafted the manuscript. SE and SW designed the mathematical model and implemented the model in the SimInf R package. All authors participated in the epidemiological analyses and critically reviewing the paper. All authors read and approved the final manuscript.
Acknowledgements
The simulations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at UPPMAX.
Funding
This work was financially supported by the Swedish Research Council Formas (Stefan Widgren, Stefan Engblom, Ulf Emanuelson, Ann Lindberg), by the Swedish Research Council within the UPMARC Linnaeus center of Excellence (Stefan Engblom), and by the Swedish strategic research program eSSENCE (Stefan Widgren).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Hancock D, Besser T, Lejeune J, Davis M, Rice D (2001) The control of VTEC in the animal reservoir. Int J Food Microbiol 66:71–78View ArticlePubMedGoogle Scholar
- Naylor SW, Nart P, Sales J, Flockhart A, Gally DL, Low C (2007) Impact of the direct application of therapeutic agents to the terminal recta of experimentally colonized calves on Escherichia coli O157:H7 shedding. Appl Environ Microbiol 73:1493–1500View ArticlePubMedPubMed CentralGoogle Scholar
- Matthews L, Reeve R, Gally DL, Low JC, Woolhouse MEJ, McAteer SP, Locking ME, Chase-Topping ME, Haydon DT, Allison LJ, Hanson MF, Gunn GJ, Reid SWJ (2013) Predicting the public health benefit of vaccinating cattle against Escherichia coli O157. Proc Natl Acad Sci U S A 110:16265–16270View ArticlePubMedPubMed CentralGoogle Scholar
- Albihn A, Eriksson E, Wallen C, Aspán A (2003) Verotoxinogenic Escherichia coli (VTEC) O157:H7—a nationwide swedish survey of bovine faeces. Acta Vet Scand 44:1–10View ArticleGoogle Scholar
- Boqvist S, Aspan A, Eriksson E (2009) Prevalence of verotoxigenic Escherichia coli O157:H7 in fecal and ear samples from slaughtered cattle in Sweden. J Food Prot 72:1709–1712View ArticlePubMedGoogle Scholar
- Widgren S, Söderlund R, Eriksson E, Fasth C, Aspan A, Emanuelson U, Alenius S, Lindberg A (2015) Longitudinal observational study over 38 months of verotoxigenic Escherichia coli O157:H7 status in 126 cattle herds. Prev Vet Med 121:343–352View ArticlePubMedGoogle Scholar
- Stärk KDC, Regula G, Hernandez J, Knopf L, Fuchs K, Morris RS, Davies P (2006) Concepts for risk-based surveillance in the field of veterinary medicine and veterinary public health: review of current approaches. BMC Health Serv Res 6:20View ArticlePubMedPubMed CentralGoogle Scholar
- Widgren S, Engblom S, Bauer P, Frössling J, Emanuelson U, Lindberg A (2016) Data-driven network modeling of disease transmission using complete population movement data: spread of VTEC O157 in Swedish cattle. Vet Res 47:81View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen EM, Skov MN, Madsen JJ, Lodal J, Jespersen JB, Baggesen DL (2004) Verocytotoxin-producing Escherichia coli in wild birds and rodents in close proximity to farms. Appl Environ Microbiol 70:6944–6947View ArticlePubMedPubMed CentralGoogle Scholar
- Herbert LJ, Vali L, Hoyle DV, Innocent G, McKendrick IJ, Pearce MC, Mellor D, Porphyre T, Locking M, Allison L, Hanson L, Matthews L, Gunn GJ, Woolhouse MEJ, Chase-Topping ME (2014) E. coli O157 on Scottish cattle farms: evidence of local spread and persistence using repeat cross-sectional data. BMC Vet Res 10:95View ArticlePubMedPubMed CentralGoogle Scholar
- Gautam R, Bani-Yaghoub M, Neill WH, Dopfer D, Kaspar C, Ivanek R (2011) Modeling the effect of seasonal variation in ambient temperature on the transmission dynamics of a pathogen with a free-living stage: example of Escherichia coli O157:H7 in a dairy herd. Prev Vet Med 102:10–21View ArticlePubMedGoogle Scholar
- Sheng H, Shringi S, Baker KNK, Minnich SA, Hovde CJ, Besser TE (2016) Standardized Escherichia coli O157:H7 exposure studies in cattle provide evidence that bovine factors do not drive increased summertime colonization. Appl Environ Microbiol 82:964–971View ArticlePubMedPubMed CentralGoogle Scholar
- Widgren S, Frössling J (2010) Spatio-temporal evaluation of cattle trade in Sweden: description of a grid network visualization technique. Geospat Health 5:119–130View ArticlePubMedGoogle Scholar
- Anonymous (2016) Maps with the average date for the beginning of each season (spring, summer, autumn, winter) for the reference period 1961–1990. The Swedish Meteorological and Hydrological Institute (SMHI), Norrköping, Sweden. http://opendata-catalog.smhi.se/. Accessed 25 Aug 2016
- Robinson SE, Wright EJ, Hart CA, Bennett M, French NP (2004) Intermittent and persistent shedding of Escherichia coli O157 in cohorts of naturally infected calves. J Appl Microbiol 97:1045–1053View ArticlePubMedGoogle Scholar
- Matthews L, Low JC, Gally DL, Pearce MC, Mellor DJ, Heesterbeek JAP, Chase-Topping M, Naylor SW, Shaw DJ, Reid SWJ, Gunn GJ, Woolhouse MEJ (2006) Heterogeneous shedding of Escherichia coli O157 in cattle and its implications for control. Proc Natl Acad Sci U S A 103:547–552View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson SE, Brown PE, Wright EJ, Hart CA, French NP (2009) Quantifying within- and between-animal variation and uncertainty associated with counts of Escherichia coli O157 occurring in naturally infected cattle faeces. J R Soc Interface 6:169–177View ArticlePubMedGoogle Scholar
- Stanford K, Agopsowicz CA, McAllister TA (2012) Genetic diversity and antimicrobial resistance among isolates of Escherichia coli O157:H7 from feces and hides of super-shedders and low-shedding pen-mates in two commercial beef feedlots. BMC Vet Res 8:178View ArticlePubMedPubMed CentralGoogle Scholar
- Buttner K, Krieter J, Traulsen A, Traulsen I (2016) Epidemic spreading in an animal trade network—comparison of distance-based and network-based control measures. Transbound Emerg Dis 63:e122–e134View ArticlePubMedGoogle Scholar
- Grenfell B, Harwood J (1997) (Meta)population dynamics of infectious diseases. Trends Ecol Evol 12:395–399View ArticlePubMedGoogle Scholar
- Bartlett MS (1957) Measles periodicity and community size. J R Stat Soc Ser A 120:48–70View ArticleGoogle Scholar
- Bui TMA, Ziane M, Stinckwich S, Ho TV, Roche B, Papoulias N (2016) Separation of concerns in epidemiological modelling. In: 15th international conference on modularity, ACM, pp 196–200Google Scholar
- Picault S, Huang YL, Sicard V, Beaudeau F, Ezanno P (2017) A multi-level multi-agent simulation framework in animal epidemiology. In: International conference on practical applications of agents and multi-agent systems, LNCS, 10349, Springer, pp 209–221Google Scholar
- Widgren S, Bauer P, Engblom S (2016) SimInf: a framework for data-driven stochastic disease spread simulations. R package available on CRAN. https://CRAN.R-project.org/package=SimInf. Accessed 12 Oct 2016
- Gillespie DT (1977) Exact stochastic simulation of coupled chemical reactions. J Phys Chem 81:2340–2361View ArticleGoogle Scholar
- Pebesma EJ, Bivand RS (2005) Classes and methods for spatial data in R. R News 5:9–13Google Scholar
- Anonymous (2012) Surveillance of infectious diseases in animals and humans in Sweden 2012. National Veterinary Institute (SVA), Uppsala, Sweden. http://www.sva.se/globalassets/redesign2011/pdf/om_sva/publikationer/surveillance2012.pdf. Accessed 9 Sept 2017
- Bauer P, Engblom S, Widgren S (2016) Fast event-based epidemiological simulations on national scales. Int J High Perform Comput Appl 30:438–453View ArticleGoogle Scholar
- Arnold ME, Ellis-Iversen J, Cook AJ, Davies RH, McLaren IM, Kay AC, Pritchard GC (2008) Investigation into the effectiveness of pooled fecal samples for detection of verocytotoxin-producing Escherichia coli O157 in cattle. J Vet Diagn Invest 20:21–27View ArticlePubMedGoogle Scholar
- Widgren S, Eriksson E, Aspan A, Emanuelson U, Alenius S, Lindberg A (2013) Environmental sampling for evaluating verotoxigenic Escherichia coli O157:H7 status in dairy cattle herds. J Vet Diagn Invest 25:189–198View ArticlePubMedGoogle Scholar
- Nelder JA, Mead R (1965) A simplex method for function minimization. Comput J 7:308–313View ArticleGoogle Scholar
- R Core Team (2016) R: a language and environment for statistical computing. R Foundation for Statistical Computing, ViennaGoogle Scholar
- Fenlon DR, Ogden ID, Vinten A, Svoboda I (2000) The fate of Escherichia coli and E. coli O157 in cattle slurry after application to land. Symp Ser Soc Appl Microbiol 88:149S–156SView ArticleGoogle Scholar
- Ogden LD, Fenlon DR, Vinten AJ, Lewis D (2001) The fate of Escherichia coli O157 in soil and its potential to contaminate drinking water. Int J Food Microbiol 66:111–117View ArticlePubMedGoogle Scholar
- Strachan NJ, Fenlon DR, Ogden ID (2001) Modelling the vector pathway and infection of humans in an environmental outbreak of Escherichia coli O157. FEMS Microbiol Lett 203:69–73View ArticlePubMedGoogle Scholar
- Davis MA, Rice DH, Sheng H, Hancock DD, Besser TE, Cobbold R, Hovde CJ (2006) Comparison of cultures from rectoanal-junction mucosal swabs and feces for detection of Escherichia coli O157 in dairy heifers. Appl Environ Microbiol 72:3766–3770View ArticlePubMedPubMed CentralGoogle Scholar
- Widgren S, Bauer P, Eriksson R, Engblom S (2018) SimInf: an R package for data-driven stochastic disease spread simulations. ArXiv e-prints, arXiv:1605.01421 [q-bio.PE]:http://arxiv.org/abs/1605.01421
- Saltelli A, Bollardo R (1998) An alternative way to compute Fourier amplitude sensitivity test (FAST). Comput Stat Data Anal 26:445–460View ArticleGoogle Scholar
- Marino S, Hogue IB, Ray CJ, Kirschner DE (2008) A methodology for performing global uncertainty and sensitivity analysis in systems biology. J Theor Biol 254:178–196View ArticlePubMedPubMed CentralGoogle Scholar
- Alden K, Read M, Andrews P, Cosgrove J, Coles M, Timmis J (2018) spartan: Simulation Parameter Analysis R Toolkit ApplicatioN: ‘spartan’. R package avialable on CRAN. https://CRAN.R-project.org/package=spartan. Accessed 21 June 2018
- Wasserman S, Faust K (1994) Social network analysis: methods and applications, vol 8. Cambridge University Press, CambridgeView ArticleGoogle Scholar
- Nöremark M, Håkansson N, Sternberg-Lewerin S, Lindberg A, Jonsson A (2011) Network analysis of cattle and pig movements in Sweden: measures relevant for disease control and risk based surveillance. Prev Vet Med 99:78–90View ArticlePubMedGoogle Scholar
- Dubé C, Ribble C, Kelton D, McNab B (2008) Comparing network analysis measures to determine potential epidemic size of highly contagious exotic diseases in fragmented monthly networks of dairy cattle movements in Ontario, Canada. Transbound Emerg Dis 55:382–392View ArticlePubMedGoogle Scholar
- Nöremark M, Widgren S (2014) EpiContactTrace: an R-package for contact tracing during livestock disease outbreaks and for risk-based surveillance. BMC Vet Res 10:71View ArticlePubMedPubMed CentralGoogle Scholar
- Snedeker KG, Campbell M, Sargeant JM (2012) A systematic review of vaccinations to reduce the shedding of Escherichia coli O157 in the faeces of domestic ruminants. Zoonoses Public Health 59:126–138View ArticlePubMedGoogle Scholar
- Varela NP, Dick P, Wilson J (2013) Assessing the existing information on the efficacy of bovine vaccination against Escherichia coli O157:H7—a systematic review and meta-analysis. Zoonoses Public Health 60:253–268View ArticlePubMedGoogle Scholar
- LeJeune JT, Wetzel AN (2007) Preharvest control of O157 in cattle. J Anim Sci 85:E73–E80View ArticlePubMedGoogle Scholar
- Eriksson E, Aspan A, Gunnarsson A, Vågsholm I (2005) Prevalence of verotoxin-producing Escherichia coli (VTEC) 0157 in Swedish dairy herds. Epidemiol Infect 133:349–358View ArticlePubMedPubMed CentralGoogle Scholar
- Nöremark M, Håkansson N, Lindström T, Wennergren U, Lewerin SS (2009) Spatial and temporal investigations of reported movements, births and deaths of cattle and pigs in Sweden. Acta Vet Scand 51:37View ArticlePubMedPubMed CentralGoogle Scholar
- Turner J, Bowers RG, Begon M, Robinson SE, French NP (2006) A semi-stochastic model of the transmission of Escherichia coli O157 in a typical UK dairy herd: dynamics, sensitivity analysis and intervention/prevention strategies. J Theor Biol 241:806–822View ArticlePubMedGoogle Scholar
- Zhang X, Woolhouse MEJ (2011) Escherichia coli O157 infection on Scottish cattle farms: dynamics and control. J R Soc Interface 8:1051–1058View ArticlePubMedGoogle Scholar
- Turner J, Begon M, Bowers RG, French NP (2003) A model appropriate to the transmission of a human food-borne pathogen in a multigroup managed herd. Prev Vet Med 57:175–198View ArticlePubMedGoogle Scholar
- Gates CM, Woolhouse MEJ (2015) Controlling infectious disease through the targeted manipulation of contact network structure. Epidemics 12:11–19View ArticlePubMedPubMed CentralGoogle Scholar
- Williams KJ, Ward MP, Dhungyel OP, Hall EJS (2015) Risk factors for Escherichia coli O157 shedding and super-shedding by dairy heifers at pasture. Epidemiol Infect 143:1004–1015View ArticlePubMedGoogle Scholar
- Cobbold RN, Rice DH, Szymanski M, Call DR, Hancock DD (2004) Comparison of shiga-toxigenic Escherichia coli prevalences among dairy, feedlot, and cow-calf herds in Washington State. Appl Environ Microbiol 70:4375–4378View ArticlePubMedPubMed CentralGoogle Scholar
- Cobbaut K, Berkvens D, Houf K, De Deken R, De Zutter L (2009) Escherichia coli O157 prevalence in different cattle farm types and identification of potential risk factors. J Food Prot 72:1848–1853View ArticlePubMedGoogle Scholar
- Zhang XS, Chase-Topping ME, McKendrick IJ, Savill NJ, Woolhouse MEJ (2010) Spread of E. coli O157 infection among Scottish cattle farms: stochastic models and model selection. Epidemics 2:11–20View ArticlePubMedPubMed CentralGoogle Scholar