Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Dec 14;13(12):e0203177.
doi: 10.1371/journal.pone.0203177. eCollection 2018.

A data-driven individual-based model of infectious disease in livestock operation: A validation study for paratuberculosis

Affiliations

A data-driven individual-based model of infectious disease in livestock operation: A validation study for paratuberculosis

Mohammad A Al-Mamun et al. PLoS One. .

Abstract

Chronic livestock diseases cause large financial loss and affect animal health and welfare. Controlling these diseases mostly requires precise information on both individual animal and population dynamics to inform the farmer's decisions, but even successful control programmes do by no means assure elimination. Mathematical models provide opportunities to test different control and elimination options rather than implementing them in real herds, but these models require robust parameter estimation and validation. Fitting these models to data is a difficult task due to heterogeneities in livestock processes. In this paper, we develop an infectious disease modeling framework for a livestock disease (paratuberculosis) that is caused by Mycobacterium avium subsp. paratuberculosis (MAP). Infection with MAP leads to reduced milk production, pregnancy rates, and slaughter value and increased culling rates in cattle and causes significant economic losses to the dairy industry. These economic effects are particularly important motivations in the control and elimination of MAP. In this framework, an individual-based model (IBM) of a dairy herd was built and MAP infection dynamics was integrated. Once the model produced realistic dynamics of MAP infection, we implemented an evaluation method by fitting it to data from three dairy herds from the Northeast region of the US. The model fitting exercises used least-squares and parameter space searching methods to obtain the best-fitted values of selected parameters. The best set of parameters were used to model the effect of interventions. The results show that the presented model can complement real herd statistics where the intervention strategies suggest a reduction in MAP prevalence without elimination. Overall, this research not only provides a complete model for MAP infection dynamics in a dairy herd but also offers a method for estimating parameters by fitting IBM models.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. A flow diagram of animal movements among infection categories for the adult, calves, and heifers within the herd.
Each horizontal gray box classifies the animals according to their initial age group. The green and red boxes define the susceptible and infected states, respectively, for each animal in the three age categories. The probabilities of exit at each time point from susceptible to latent, latent to low shedding and low shedding to high shedding animals are s1, h1, and y1, respectively. Vertical transmission probabilities from latent, low shedding and high shedding animals are Vh, Vy1, and Vy2, respectively. Horizontal transmission probabilities to calves from low shedding and high shedding animals are Hy1 and Hy2, respectively. The probability an animal gets infected by the environment is βenvironment. Calf-to-calf and heifer-to-heifer transmission probabilities are Cinf and Yinf, respectively. Stochastic death/sale probabilities for adult, calves, and heifers are μa, μc, and μh, respectively. μ is the replacement animals coming from heifer compartment upon completion of two years.
Fig 2
Fig 2. The comparison of observed and model predicted milk yield distribution for 1% simulation using best-fit parameters for the milk yield.
In the box plot, the bottom and top end of the bars are minimum and maximum values respectively, the top of the box is the 75th percentile, the bottom of the box is the 25th percentile, and the horizontal line within the box is median; outliers are presented as a solid black circle and the density of the milk yield is presented by the width of the violin.
Fig 3
Fig 3. The fitting results of three in silico herds A (top), B (middle), and C (bottom) compared to the observed apparent prevalence for 7 years by biannual sampling.
The shaded region shows the 95% confidence interval of the best 1% simulation runs.
Fig 4
Fig 4. The apparent prevalence during the pre- and post-intervention period during the simulation of three in silico herds with two control strategies.
Control I: culling red animals immediately and control II: culling only red animal with a delay of 305 days in milk. The two control measures are simulated in separating runs of the three in silico herds.

References

    1. BROUGHAN JM, JUDGE J, ELY E, DELAHAY RJ, WILSON G, CLIFTON-HADLEY RS, et al. A review of risk factors for bovine tuberculosis infection in cattle in the UK and Ireland. Epidemiol Infect. Cambridge University Press; 2016;144: 2899–2926. 10.1017/S095026881600131X - DOI - PMC - PubMed
    1. Rathnaiah G, Zinniel DK, Bannantine JP, Stabel JR, Gröhn YT, Collins MT, et al. Pathogenesis, Molecular Genetics, and Genomics of Mycobacterium avium subsp. paratuberculosis, the Etiologic Agent of Johne’s Disease. Front Vet Sci. 2017;4: 187 10.3389/fvets.2017.00187 - DOI - PMC - PubMed
    1. Brooks-Pollock E, Roberts GO, Keeling MJ. A dynamic model of bovine tuberculosis spread and control in Great Britain. Nature. Nature Publishing Group; 2014;511: 228–231. 10.1038/nature13529 - DOI - PubMed
    1. NAHMS. National Animal Health Monitoring System (2007) Dairy 2007-Johne’s Disease on U.S. Dairies, 1991–2007,National Animal Health Monitoring Systems. Fort Collins, CO; 2007.
    1. Behr M, Collin D. Paratuberculosis:Organism, Disease, Control. CABI; 2010.

Publication types

MeSH terms