 Research article
 Open Access
 Published:
Invasion and transmission of Salmonella Kentucky in an adult dairy herd using approximate Bayesian computation
BMC Veterinary Research volume 9, Article number: 245 (2013)
Abstract
Background
An outbreak of Salmonella Kentucky followed by a high level of sustained endemic prevalence was recently observed in a US adult dairy herd enrolled in a longitudinal study involving intensive fecal sampling. To understand the invasion ability and transmission dynamics of Salmonella Kentucky in dairy cattle, accurate estimation of the key epidemiological parameters from longitudinal field data is necessary. The approximate Bayesian computation technique was applied for estimating the transmission rate (β), the recovery rate (γ) and shape (n) parameters of the gamma distribution for the infectious (shedding) period, and the basic reproduction ratio (R_{0}), given a susceptibleinfectiousrecoveredsusceptible (SIRS) compartment model with a gamma distribution for the infectious period.
Results
The results report that the mean transmission rate (β) is 0.417 month^{1} (median: 0.417, 95% credible interval [0.406, 0.429]), the average infectious period (γ^{1}) is 7.95 months (median: 7.95, 95% credible interval [7.70, 8.22]), the mean shape parameter (n) of the gamma distribution for the infectious period is 242 (median: 182, 95% credible interval [16, 482]), and the mean basic reproduction ratio (R_{0}) is 2.91 (median: 2.91, 95% credible interval [2.83, 3.00]).
Conclusions
This study shows that Salmonella Kentucky in this herd was of mild infectiousness and had a long infectious period, which together provide an explanation for the observed prevalence pattern after invasion. The transmission rate and the recovery rate parameters are inferred with better accuracy than the shape parameter, therefore these two parameters are more sensitive to the model and the observed data. The estimated shape parameter (n) has large variability with a minimal value greater than one, indicating that the infectious period of Salmonella Kentucky in dairy cattle does not follow the conventionally assumed exponential distribution.
Background
Salmonella is one of the major causes of foodborne gastroenteritis worldwide and poses a considerable threat to public health. In the United States (US) alone there are approximately 1.4 million illnesses, 16,000 hospitalizations, and 400 deaths annually [1, 2]. Humans generally acquire salmonellosis through consumption of contaminated food or contact with infected animals or a contaminated environment [3, 4]. Emergence of multidrugresistant Salmonella in human infections is particularly serious due to increased morbidity and mortality [5].
More than 2,500 serotypes of Salmonella have been identified and significant variability has been found in virulence, infectious dose, and host. Most human Salmonella infections are caused by relatively few Salmonella subtypes (S. enterica Enteritidis, Typhimurium, and NewPort) [6], but all Salmonella serotypes are considered potentially pathogenic [7]. Farm animals are recognized as important reservoirs for Salmonella and other foodborne pathogens [8]. Many Salmonella serotypes have been found in samples from dairy animals and their environment, some of which have also been isolated in human cases [9]. Reduction of Salmonella prevalence in farm animals and bacterial loads in the contaminated environment is important to decrease the risk of zoonotic Salmonella infection [8].
Aiming to ensure a safe food supply through identifying pathogen transmission pathways and subsequent best management practices in dairy farms, the Regional Dairy Quality Management Alliance (RDQMA) and the Agriculture Research Service (ARS) of the USDA established a longitudinal observational study on three commercial dairy farms in the Northeastern US [10], one of which was located in Pennsylvania. Fecal samples were intensively collected from all adult animals in this herd beginning in the spring of 2004. These samples were tested for Salmonella and a number of other foodborne pathogens [10, 11]. After an initial outbreak of Salmonella Cerro, a subsequent outbreak of Salmonella Kentucky with a high level of prevalence and longterm endemic infection in this farm was found [12]. Drugresistant Salmonella Kentucky (ST198) from human cases has recently been identified from a study of National Salmonella Surveillance Systems from France, England and Wales, Denmark, and the United States [13]. The identification of Salmonella Kentucky is also a common occurrence in poultry [14].
Understanding the transmission dynamics that underlie observed shedding patterns from longitudinal field data is essential for the effective design of Salmonella prevention and intervention. As the transmission dynamics of Salmonella spp. are determined by complex interactions among host, pathogen, and environment, mathematical modeling approaches have been applied to provide insights in the understanding of transmission [15–19]. However, a common problem in mathematical models in epidemiology is how to accurately and reliably estimate the nonobservable model parameters such as the transmission rate, given the available longitudinal field data. To solve this problem, a number of Bayesian inference approaches for infectious disease transmission models have been developed [20–25]. Posterior distributions of parameters can generally be computed using an explicit likelihood function given parameter prior distributions with the help of Markov Chain Monte Carlo (MCMC) methods [20–22], or using a likelihoodfree approach, the approximate Bayesian computation (ABC) technique [23] with a newly proposed efficient sequential Monte Carlo algorithm [24, 25].
The objective of this study was to infer from longitudinal field data the key epidemiological parameters that are important to understand the invasion ability and transmission dynamics in an outbreak and subsequent endemicity of Salmonella Kentucky on a dairy farm. Specifically, we estimated the transmission rate (β), the rate (γ) and shape (n) parameters of the gamma distribution for the infectious (shedding) period, and the basic reproduction ratio (R_{0}) using the approximate Bayesian computation technique.
Methods
Longitudinal field data
The dairy herd (socalled “Farm B”) in Pennsylvania consisted of approximately 100110 adult cows housed in a free stall barn [10–12]. Calves from this herd were transferred to an offsite rearing center at 6 months of age and were returned to the herd as replacement animals within 1 to 2 months prior to their first calving. Intensive fecal samples were collected for all adult cows in the herd with a sampling interval of 6 to 8 weeks during the study period. Methods for isolation and serotyping of Salmonella were previously described [10–12].
The longitudinal data used in this study are composed of observed withinherd prevalence (proportion) of animals shedding Salmonella Kentucky for a total of 14 time points from January 2006 to December 2007. As the sensitivity of the culture test for Salmonella is generally estimated to be imperfect and relatively low, we used a correction to account for likely falsenegative culture results. This rule states that a negative test found between two immediate (neighboring) positive tests in an individual cow was assumed to be a falsenegative test and this negative test was corrected to be positive (+ + →+++). The specificity of the culture test for Salmonella was assumed to be one, so there were no falsepositive tests.
The SIRS model
Multiple episodes of shedding of Salmonella Kentucky were observed in individual cows’ test result profiles. Due to the assumed perfect test specificity and imperfect test sensitivity (see above), at least two or more consecutive negative tests had to appear between these testpositive episodes to define a second or higher shedding (infection) period. On the basis of the shedding pattern of individual cows, we developed a susceptibleinfectiousrecoveredsusceptible (SIRS) transmission model with the following assumptions:

(a)
First lactation animals entering the herd as replacement animals were assumed to be susceptible (S) because no positive tests were found in any heifers leaving for the offsite facility or returning to the herd.

(b)
Testpositive animals (shedding Salmonella Kentucky) were assumed to be in the infectious state (I).

(c)
The infectious (shedding) period was assumed to follow a gamma distribution f(tn, γ) = (n γ)^{n} t ^{n1} exp (nγt)/Г(n)) with two parameters, the rate (γ) and shape (n), reflecting that time since infection is important [15]. When the shape parameter (n) is one, the gamma distribution reduces to an exponential distribution.

(d)
Susceptible animals (S) remained in their susceptible state until they became infectious (I) at rate β I/ N(force of infection), where the transmission rate is denoted by β.

(e)
The time period covering two or more consecutive negative tests in individual animals between the two neighboring positive tests was relatively long, approximately 3 months or more. Within that time period, animals were assumed to build their immunity, becoming recovered (R). Eventually during that period, the recovered animals (R) lost their immunity and became susceptible (S) again at rate ϕ.

(f)
Susceptible animals were assumed not to be distinguished based on the presence or absence of prior exposure to Salmonella (i.e. no immunological memory).

(g)
Direct cowcow transmission was considered to represent the common fecaloral transmission route of Salmonella Kentucky in the herd.

(h)
Herd size was assumed to be a constant (108 cows in this study based on the average herd size across all 14 sampling time points), reflecting a constant size of the milk producing herd.

(i)
To maintain a constant herd size, we assumed the replacement rate (μ) of the milking herd was the same as the removal rate (μ) of animals from the herd.
There are a total of 5 parameters in the SIRS model described in Figure 1. Three of these are unknown and to be estimated: the transmission rate (β) and the rate (γ) and shape (n) parameters describing the gamma distribution for the infectious period. The replacement rate (μ, or the removal rate) of animals was assumed to be constant and was calculated from the animal movement data to be 0.03 (month)^{1}, representing the average duration of survival of cows in the milking herd at approximately 2.8 years. The rate of immunity loss (ϕ = 0.33 (month)^{1}) was also assumed constant and calculated to be around 3 months because it needed to be equal to or less than the time period for at least two consecutive negative tests (covering the recovered (R) state and its next susceptible (S) state).
The system of ordinary differential equations describing the transmission dynamic model (SIRS) with a gamma distribution for the infectious period in Figure 1 is given in Additional file 1: Appendix A.
Method of estimating parameters
We applied the approximate Bayesian computation (ABC) technique [23] to infer the unknown epidemiological parameters (β, γ, n and R_{0}) given the longitudinal data and the SIRS model. Uniform (flat) prior distributions were assumed for the 3 unknown parameters, β ∈ [0.01,2], γ ∈ [0.01,1], and n ∈ [1, 500].
Unknown parameters (β, γ, n) were sampled from their prior distributions and these sampled values were used to numerically solve the system of ordinary differential equations (the SIRS model in Additional file 1: Appendix A). The sum of squared errors between the fitted and observed prevalence was calculated. If the sum of squared errors was less than a desired tolerance value, then the sampled parameter values were accepted. However this rejection algorithm was not effective and was not able to put into in practice due to high computational demands. An efficient algorithm recently developed for ABC using the sequential Monte Carlo method was implemented for parameter estimate and model selection for nonlinear dynamic systems [24, 25]. In this study we used this efficient algorithm for parameter estimation.
To examine the invasion ability, we estimated the basic reproduction ratio (R_{0}), which represents the number of secondary cases caused by the introduction of a primary index case into a fully susceptible population during its whole infectious period. The basic reproduction ratio (R_{0}) of the SIRS model in this study is a function of the other three parameters (β, γ, n)and was obtained using the nextgeneration matrix [26]:
Posterior predictive check and crossvalidation
Using the estimates from the posterior distributions of parameters (β, γ, n), stochastic simulations based on the direct Gillespie algorithm for the SIRS model were performed [27]. The expected prevalence and its 2.5% and 97.5% quantiles were compared to the observed prevalence. As we had only one dataset of the observed longitudinal prevalence (consisting of 14 data points), crossvalidation was performed to check the behavior of the SIRS model. A total of 14 datasets for crossvalidation were formed by removing each data point once from the full observed prevalence data set.
Impact of the rate of loss immunity
As the rate at which immunity wanes (ϕ) was uncertain and partly based on the sampling interval, we repeated the ABC analyses with a slower rate, 0.25 (month)^{1}, to investigate the effects of varying the rate of immunity loss (ϕ) on the estimation of the transmission rate (β), the rate (γ) and shape (n) parameters for the gamma distribution for the infectious period, and the basic reproduction ratio (R_{0}).
Results
Posterior distributions of parameters (β, γ, n)
Posterior distributions of the transmission rate (β), and the rate (γ) and shape (n) parameters of the gamma distribution describing the infectious period are shown in the top row of Figure 2. All distributions were unimodal. For the transmission rate (β) and the rate of recovery (γ) the ranges of these two posterior distributions were relatively narrow. The transmission rate (β) showed a mean of 0.417 month^{1} with a 95% credible interval of 0.406 to 0.429. The average recovery rate (γ) was 0.126 month^{1} with a 95% credible interval of 0.1216 to 0.1298 month^{1}. The shape parameter (n) had large variability with a median of 182 (mean: 242) with a 95% credible interval of 16 to 482.
The bottom row of Figure 2 shows scatter plots of paired parameters. Each scatter plot illustrates three time snapshots from the start of simulation to the middle and finally to the stable state – the process of how the sequential Monte Carlo algorithm in ABC leading to the convergence posterior distribution.
The infectious period and the basic reproduction ratio (R_{0})
The distributions of the infectious period and the basic reproduction ratio (R_{0}) using the estimates from the posterior distributions of parameters (β, γ, n) are shown in Figure 3a and 3b, respectively. The average infectious period had a mean of 7.95 months (1/γ = 1/ 0.1258) with a 95% credible interval of 7.70 (1/γ = 1/ 0.1298) to 8.22 (1/γ = 1/ 0.1216) months. The R_{0} had a mean of 2.91 with a 95% credible interval of 2.83 to 3.00.
Posterior predictive check
Comparison between the observed prevalence and the expected (fitted) prevalence with the predicted 2.5% and 97.5% quantile is demonstrated in Figure 4. The expected prevalence exhibited agreed well with the observed prevalence before Salmonella Kentucky reached endemic prevalence. All of the observed prevalence data points were within the range from the predicted 2.5% to 97.5% quantile; two observed prevalence data points, the 12^{th} and 14^{th}, were quite close to the low 2.5% and high 97.5% quantiles, respectively.
Crossvalidation
Boxplots of the transmission rate (β), the recovery rate (γ), the infectious period, and the basic reproduction ratio (R_{0}) are shown for 14 subdatasets in Figure 5. The transmission rates (β) mostly ranged from 0.39 to 0.46 month^{1}, and the recovery rate parameters (γ) from 0.112 to 0.138 month^{1}. The infectious period was mostly between 7 and 9 months, and the R_{0} value was between 2.70 and 3.14 for almost all datasets.
Impact of varying the rate of immunity
We studied the impacts of varying the rate of immunity (ϕ) on the transmission rate (β), the rate (γ) and shape (n) parameters of the gamma distribution for the infectious period, and the basic reproduction ratio (R_{0}). When the rate of immunity (ϕ) was changed from 0.33 (month)^{1} to 0.25 (month)^{1} (i.e., the average duration in the recovery state (R) thereby increased from 3 months to 4 months), the estimated transmission rate (β) parameter (mean: 0.417 → 0.428 month^{1}; 95% credible interval [0.1216, 0.1298] → [0.1137, 0.1249] month^{1}), and the rate of recovery (γ) parameter (mean: 0.1258 → 0.1195 month^{1}; 95% credible interval [0.1216, 0.1298] → [0.1137, 0.1249] month^{1}), and the basic reproduction ratio (R_{0}) (mean: 2.9 → 3.09; 95% credible interval [2.83, 3.00] → [2.91, 3.31]) changed only slightly, but the shape parameter (n) had a relatively large change and was well inferred (median: 182 → 11; 95% credible interval [16, 482] → [4, 32]). Posterior distributions of these parameters are given in Figure B1 in Additional file 1: Appendix B.
Discussion
The mean estimate of the basic reproduction ratio (R_{0}) for Salmonella Kentucky in adult dairy cows from our data is 2.91, which indicates a relatively mild infectiousness (one primary shedding cow can, on average, infect 3 susceptible cows). The relatively modest R_{0} value also indicates that preventative efforts to reduce the reproduction rate to values below 1 may be reasonable and even socalled ‘leaky’ vaccines may prove to be sufficiently efficacious to provide herd immunity [17].
The mean duration of the infectious (or shedding) period is 8 months. This estimated duration is long compared to other Salmonella strains causing clinical signs in dairy cattle, such as Salmonella Typhimurium and Salmonella Dublin [18]. The mild infectiousness and long infectious period together explained the observed dynamic pattern of Salmonella Kentucky in Figure 4; the prevalence of animals shedding Salmonella Kentucky gradually increased during the epidemic phase and then a relatively stable longterm endemic infection was established in the herd.
The posterior distribution of the shape parameter (n = 182, 95% credible interval [16, 482]) of the gamma distribution for the infectious period did not include n = 1; therefore, the conventional assumption that the infectious period is described by an exponential distribution was found to be inappropriate and a more realistic gamma distribution for the infectious period was favored. This indicates that time since infection is important in the transmission dynamics of Salmonella Kentucky in adult dairy cattle. Consequently, strategies for prevention and intervention could be affected by the infection time of animals in the herd [28]. Although large variability of the shape parameter is observed with an assumption of a 3month period of immunity loss, it can be significantly reduced if a 4month period of immunity loss is assumed (Additional file 1: Appendix B).
We expect that for many more persistent infections, the assumption of an exponential rate of disappearance from the infectious state will be incorrect. Implementation of nonexponential distributions in ODEbased parameter estimation is not straightforward. The ability to use distributions other than exponential, with the gamma distribution being an attractive alternative, may be one of the key benefits of using ABC for parameter estimation, as no explicit likelihood function needed to be defined. This was especially helpful in this study, where we used a gamma distribution for the infectious period. However, the lack of an explicit likelihood function requires highly demanding computational efforts. Therefore, the implementation of an efficient algorithm in ABC as shown here became imperative.
Posterior distributions of the transmission rate and the rate of recovery (β, γ) obtained from ABC (Figure 2) were relatively narrow. This was also described in a previous study [24], partly due to the use of a deterministic model (a system of differential equations) in the simulation of ABC. If a stochastic SIRS model (in the formulation of either stochastic differential equations or continuoustime Markov Chain) was used in ABC, the credible intervals would be wider due to the addition of stochasticity into the SIRS model. The transmission rate (β) and the recovery rate (γ) parameters had significantly better accuracy than the shape parameter (n) as shown in the top row of Figure 2, which indicated that these two parameters were more sensitive to the model and data than the shape parameter (n) [24].
We did not apply ABC to the stochastic SIRS model because the estimated values of the transmission rate (β), the rate of recovery (γ) and the shape parameter (n) from the deterministic SIRS model were sufficient to capture the observed prevalence pattern (Figure 4) using a stochastic SIRS model implementing the direct Gillespie algorithm. We also did not perform elaborate model selection in this study because the SIRS model appears to be capable of explaining the observed transmission dynamics of Salmonella Kentucky. As mentioned before, the observed intermittent shedding may be explained by the relatively poor sensitivity of culture methods and we therefore corrected for the occasional assumed falsenegative result. Further research to distinguish between assuming true intermittent shedding [29] and a nonperfect test sensitivity in continuous shedding may be necessary.
The posterior predictive check shown in Figure 4 indicated that the estimated transmission rate (β) and the rate (γ) and shape (n) parameters in Figure 2 were reasonable, as stochastic simulations for the SIRS model with these parameter estimates from their posterior distributions were able to capture the observed dynamic (prevalence) pattern. The cross validation (Figure 5) suggested that the estimate of the basic reproduction ratio (R_{0}) was fairly consistent. Although the shape parameter had large variability (the top row of Figure 2c), it did not have a substantial effect on the basic reproduction ratio. In other words, neither the basic reproduction ratio (R_{0}) nor the dynamic pattern were sensitive to changes in the shape parameter (n).
When varying the period of immunity from 3 months to 4 months, no significant changes were found in the transmission rate and the rate of the gamma distribution for the infectious period. However, serological data that are able to distinguish the recovered state (R) and the susceptible (S) state would help remove the uncertainty in the rate of immunity loss and increase the accuracy of parameter estimation, especially for the estimation of the shape parameter (n).
Conclusions
We developed a susceptibleinfectiousrecoveredsusceptible (SIRS) model to describe the transmission of Salmonella Kentucky in an adult dairy herd. The important epidemiological parameters of the SIRS model were estimated from a longitudinal data set using the approximate Bayesian computation method. This study shows that Salmonella Kentucky has a mild invasion ability (R_{0 =} 2.91, 95% credible interval [2.83, 3.00]) and has a long average infectious period (7.95 months, 95% credible interval [7.70, 8.22]) in dairy cattle. These findings together provide an explanation for the observed prevalence pattern after invasion. The transmission rate and the recovery rate parameters are inferred with better accuracy than the shape parameter, therefore these two parameters are more sensitive to the model and the observed data. The estimated shape parameter (n) has large variability with a minimal value greater than one, indicating that the infectious period of Salmonella Kentucky in dairy cattle does not follow the conventionally assumed exponential distribution.
References
 1.
Mead PS, Slutsker L, Dietz V, McCaig LF, Bresee JS, Shapiro C, Griffin PM, Tauxe RV: Foodrelated illness and death in the United States. Emerg Infect Dis. 1999, 5: 607625. 10.3201/eid0505.990502.
 2.
Voetsch AC, Van Gilder TJ, Angulo FJ, Farley MM, Shallow S, Marcus R, Cieslak PR, Deneen VC, Tauxe RV: FoodNet estimate of the burden of illness caused by nontyphoidal Salmonella infections in the United States. Clin Infect Dis: Offic Publ Infect Dis Soc Am. 2004, 38 (Suppl 3): S127S134.
 3.
Majowicz SE, Musto J, Scallan E, Angulo FJ, Kirk M, O'Brien SJ, Jones TF, Fazil A, Hoekstra RM: The global burden of nontyphoidal Salmonella gastroenteritis. Clin Infect Dis: Offic Publ Infect Dis Soc Am. 2010, 50: 882889. 10.1086/650733.
 4.
LP F, Wierup M: Salmonella contamination: a significant challenge to the global marketing of animal food products. Rev Sci Tech. 2006, 25: 541554.
 5.
Varma JK, Molbak K, Barrett TJ, Beebe JL, Jones TF, RabatskyEhr T, Smith KE, Vugia DJ, Chang HG, Angulo FJ: Antimicrobialresistant nontyphoidal Salmonella is associated with excess bloodstream infections and hospitalizations. J Infect Dis. 2005, 191: 554561. 10.1086/427263.
 6.
Jones TF, Ingram LA, Cieslak PR, Vugia DJ, TobinD'Angelo M, Hurd S, Medus C, Cronquist A, Angulo FJ: Salmonellosis outcomes differ substantially by serotype. J Infect Dis. 2008, 198: 109114. 10.1086/588823.
 7.
CDC: Salmonella Surveillance: Annual Summary, 2004. Altanta, GA: U.S. Department of Health and Human Servicies, Centers for Disease Control and Prevention (CDC) 2005.
 8.
Stevens MP, Humphrey TJ, Maskell DJ: Molecular insights into farm animal and zoonotic Salmonella infections. Phil Trans Roy Soc Lond Ser B Biol Sci. 2009, 364: 27092723. 10.1098/rstb.2009.0094.
 9.
Callaway TR, Keen JE, Edrington TS, Baumgard LH, Spicer L, Fonda ES, Griswold KE, Overton TR, VanAmburgh ME, Anderson RC, et al: Fecal prevalence and diversity of Salmonella species in lactating dairy cattle in four states. J Dairy Sci. 2005, 88: 36033608. 10.3168/jds.S00220302(05)730459.
 10.
Pradhan AK, Van Kessel JS, Karns JS, Wolfgang DR, Hovingh E, Nelen KA, Smith JM, Whitlock RH, Fyock T, Ladely S, et al: Dynamics of endemic infectious diseases of animal and human importance on three dairy herds in the northeastern United States. J Dairy Sci. 2009, 92: 18111825. 10.3168/jds.20081486.
 11.
Van Kessel JS, Karns JS, Wolfgang DR, Hovingh E, Schukken YH: Longitudinal study of a clonal, subclinical outbreak of Salmonella enterica subsp. enterica serovar Cerro in a U.S. dairy herd. Foodborne Pathogens Dis. 2007, 4: 449461. 10.1089/fpd.2007.0033.
 12.
Van Kessel JA, Karns JS, Wolfgang DR, Hovingh E, Schukken YH: Dynamics of Salmonella serotype shifts in an endemically infected dairy herd. Foodborne Pathogens Dis. 2012, 9: 319324. 10.1089/fpd.2011.1054.
 13.
Le Hello S, Hendriksen RS, Doublet B, Fisher I, Nielsen EM, Whichard JM, Bouchrif B, Fashae K, Granier SA, JourdanDa Silva N, et al: International spread of an epidemic population of Salmonella enterica serotype Kentucky ST198 resistant to ciprofloxacin. J Infect Dis. 2011, 204: 675684. 10.1093/infdis/jir409.
 14.
Fricke WF, McDermott PF, Mammel MK, Zhao S, Johnson TJ, Rasko DA, FedorkaCray PJ, Pedroso A, Whichard JM, Leclerc JE, et al: Antimicrobial resistanceconferring plasmids with similarity to virulence plasmids from avian pathogenic Escherichia coli strains in Salmonella enterica serovar Kentucky isolates from poultry. Appl Environ Microbiol. 2009, 75: 59635971. 10.1128/AEM.0078609.
 15.
Chapagain PP, van Kessel JS, Karns JS, Wolfgang DR, Hovingh E, Nelen KA, Schukken YH, Grohn YT: A mathematical model of the dynamics of Salmonella Cerro infection in a US dairy herd. Epidemiol Infect. 2008, 136: 263272.
 16.
Xiao Y, Bowers RG, Clancy D, French NP: Understanding the dynamics of Salmonella infections in dairy herds: a modelling approach. J Theor Biol. 2005, 233: 159175. 10.1016/j.jtbi.2004.09.015.
 17.
Lu Z, Grohn YT, Smith RL, Wolfgang DR, Van Kessel JA, Schukken YH: Assessing the potential impact of Salmonella vaccines in an endemically infected dairy herd. J Theor Biol. 2009, 259: 770784. 10.1016/j.jtbi.2009.04.028.
 18.
Van Schaik G, Klinkenberg D, Veling J, Stegeman A: Transmission of Salmonella in dairy herds quantified in the endemic situation. Vet Res. 2007, 38: 861869. 10.1051/vetres:2007036.
 19.
Lanzas C, Warnick LD, Ivanek R, Ayscue P, Nydam DV, Grohn YT: The risk and control of Salmonella outbreaks in calfraising operations: a mathematical modeling approach. Vet Res. 2008, 39: 6110.1051/vetres:2008038.
 20.
O'Neill PD: Introduction and snapshot review: relating infectious disease transmission models to data. Stat Med. 2010, 29: 20692077. 10.1002/sim.3968.
 21.
Gibson GJ, Renshaw E: Estimating parameters in stochastic compartmental models using Markov chain methods. IMA J Math Appl Med Biol. 1998, 15: 1940. 10.1093/imammb/15.1.19.
 22.
Streftaris G, Gibson GJ: Bayesian inference for stochastic epidemics in closed populations. Stat Model. 2004, 4: 6375. 10.1191/1471082X04st065oa.
 23.
Beaumont MA: Approximate Bayesian computation in evolution and ecology. Annu Rev Ecol Evol Syst. 2010, 41: 379406. 10.1146/annurevecolsys102209144621.
 24.
Toni T, Welch D, Strelkowa N, Ipsen A, Stumpf MP: Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems. J Roy Soc Interface Roy Soc. 2009, 6: 187202. 10.1098/rsif.2008.0172.
 25.
Liepe J, Barnes C, Cule E, Erguler K, Kirk P, Toni T, Stumpf MP: ABCSysBio–approximate Bayesian computation in Python with GPU support. Bioinformatics. 2010, 26: 17971799. 10.1093/bioinformatics/btq278.
 26.
van den Driessche P, Watmough J: Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002, 180: 2948. 10.1016/S00255564(02)001086.
 27.
Gillespie DT: Stochastic simulation of chemical kinetics. Annu Rev Phys Chem. 2007, 58: 3555. 10.1146/annurev.physchem.58.032806.104637.
 28.
Wearing HJ, Rohani P, Keeling MJ: Appropriate models for the management of infectious diseases. PLoS Med. 2005, 2: e17410.1371/journal.pmed.0020174.
 29.
Ivanek R, Osterberg J, Gautam R, Sternberg Lewerin S: Salmonella fecal shedding and immune responses are dose and serotype dependent in pigs. PLoS One. 2012, 7: e3466010.1371/journal.pone.0034660.
Acknowledgements
We thank Professor James Booth of the Department of Biological Statics and Computational Biology at Cornell University for his help in calibrating the tolerance value in the implementation of ABC using additional information from our data set. We also thank anonymous reviewers for their constructive and valuable comments and suggestions on our previous submission. We acknowledge the collaboration of the dairy producer who allowed us to collect all the data on his farm and the support of the field staff at Pennsylvania State University and the laboratory staff at the Agricultural Research Services.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no any competing interests.
Authors' contributions
ZL, RMM, YHS and YTG conceived the study. YHS, YTG, ZL, RMM, and RLS contributed to the model development. ZL implemented the model and performed parameter estimation. JSK, JASVK, and DRW collected the field data and analyzed the field prevalence. All authors contributed to the writing of the manuscript, and read and approved the final version of the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.
About this article
Cite this article
Lu, Z., Mitchell, R.M., Smith, R.L. et al. Invasion and transmission of Salmonella Kentucky in an adult dairy herd using approximate Bayesian computation. BMC Vet Res 9, 245 (2013). https://doi.org/10.1186/174661489245
Received:
Accepted:
Published:
Keywords
 Epidemiological modeling
 Approximate Bayesian computation
 Transmission dynamics
 Salmonella
 Dairy cattle