A new compartmental model of Mycobacterium avium subsp. paratuberculosis infection dynamics in cattle

Investigating effective testing strategies for the control of Johne's disease in western Canadian cow-calf herds using an agent-based simulation model

Johne's disease is an insidious infectious disease of ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease can have important implications for animal welfare and risks causing economic losses in affected herds due to reduced productivity, premature culling and replacement, and veterinary costs. Despite the limited accuracy of diagnostic tools, testing and culling is the primary option for controlling Johne's disease in beef herds. However, evidence to inform specific test and cull strategies is lacking. In this study, a stochastic, continuous-time agent-based model was developed to investigate Johne's disease and potential control options in a typical western Canadian cow-calf herd. The objective of this study was to compare different testing and culling scenarios that included varying the testing method and frequency as well as the number and risk profile of animals targeted for testing using the model. The relative effectiveness of each testing scenario was determined by the simulated prevalence of cattle shedding MAP after a 10-year testing period. A second objective was to compare the direct testing costs of each scenario to identify least-cost options that are the most effective at reducing within-herd disease prevalence. Whole herd testing with individual PCR at frequencies of 6 or 12 months were the most effective options for reducing disease prevalence. Scenarios that were also effective at reducing prevalence but with the lowest total testing costs included testing the whole herd with individual PCR every 24 months and testing the whole herd with pooled PCR every 12 months. The most effective method with the lowest annual testing cost per unit of prevalence reduction was individual PCR on the whole herd every 24 months. Individual PCR testing only cows that had not already been tested 4 times also ranked well when considering both final estimated prevalence at 10 years and cost per unit of gain. A more in-depth economic analysis is needed to compare the cost of testing to the cost of disease, taking into account costs of culling, replacements and impacts on calf crops, and to determine if testing is an economically attractive option for commercial cow-calf operations.

Modeling Paratuberculosis Transmission in a Small Dairy Herd Typical of Slovenia Suggests That Different Models Should Be Used to Study Disease Spread in Herds of Different Sizes

This study aimed to investigate the possible dynamics of paratuberculosis or Johne’s disease in a typical Slovenian dairy herd of about 17 cows. Paratuberculosis is a worldwide endemic disease of cattle caused by Mycobacterium avium subsp. paratuberculosis (MAP) and is associated with significant economic losses. We developed a stochastic compartmental model with two pathways of disease progression, infections of adult cows and infections of young animals through horizontal and vertical transmission, and transmission through animal movements. The average proportions of subclinically and clinically infected cows were 4% and 0.47%, respectively. The prevalence within the herd, which included latently infected animals, averaged 7.13% and ranged from 0% to 70.59%. Under the given circumstances, the results showed a relatively high rate of spontaneous elimination (0.22 per herd per year) of the disease and a high rate of reinfection (0.18 per herd per year) facilitated by active animal trade. To our knowledge, this stochastic compartmental model is the first to be developed specifically to represent a small dairy herd and could apply to other countries with a similar structure of dairy farms. The results suggest that different models should be used to study MAP spread in herds of various sizes.

Characterizing infectious disease progression through discrete states using hidden Markov models

Infectious disease management relies on accurate characterization of disease progression so that transmission can be prevented. Slowly progressing infectious diseases can be difficult to characterize because of a latency period between the time an individual is infected and when they show clinical signs of disease. The introduction of Mycobacterium avium ssp. paratuberculosis (MAP), the cause of Johne’s disease, onto a dairy farm could be undetected by farmers for years before any animal shows clinical signs of disease. In this time period infected animals may shed thousands of colony forming units. Parameterizing trajectories through disease states from infection to clinical disease can help farmers to develop control programs based on targeting individual disease state, potentially reducing both transmission and production losses due to disease. We suspect that there are two distinct progression pathways; one where animals progress to a high-shedding disease state, and another where animals maintain a low-level of shedding without clinical disease. We fit continuous-time hidden Markov models to multi-year longitudinal fecal sampling data from three US dairy farms, and estimated model parameters using a modified Baum-Welch expectation maximization algorithm. Using posterior decoding, we observed two distinct shedding patterns: cows that had observations associated with a high-shedding disease state, and cows that did not. This model framework can be employed prospectively to determine which cows are likely to progress to clinical disease and may be applied to characterize disease progression of other slowly progressing infectious diseases.

Modelling spread and surveillance of Mycobacterium avium subsp. paratuberculosis in the Swedish cattle trade network

To monitor a state of disease freedom and to ensure a timely detection of new introductions of disease, surveillance programmes need be evaluated prior to implementation. We present a strategy to evaluate surveillance of Mycobacterium avium subsp. paratuberculosis (MAP) using simulated testing of bulk milk in an infectious disease spread model. MAP is a globally distributed, chronic infectious disease with substantial animal health impact. Designing surveillance for this disease poses specific challenges because methods for surveillance evaluation have focused on estimating surveillance system sensitivity and probability of freedom from disease and do not account for spread of disease or complex and changing population structure over long periods. The aims of the study were to 1. define a model that describes the spread of MAP within and between Swedish herds; 2. define a method for simulation of imperfect diagnostic testing in this framework; 3. to compare surveillance strategies to support surveillance design choices. The results illustrate how this approach can be used to identify differences between the probability of detecting disease in the population based on choices of the number of herds sampled and the use of risk-based or random selection of these herds. The approach was also used to assess surveillance to detect introduction of disease and to detect a very low prevalence endemic state. The use of bulk milk sampling was determined to be an effective method to detect MAP in the population with as few as 500 herds tested per year if the herd-level prevalence was 0.2 %. However, detection of point introductions in the population was unlikely in the 13-year simulation period even if as many as 2000 herds were tested per year. Interestingly, the use of a risk-based selection strategy was found to be a disadvantage to detect MAP given the modelled disease dynamics.

Implications of Cattle Trade for the Spread and Control of Infectious Diseases in Slovenia

The objectives of this study were to gain insight into the structure of the cattle trade network in Slovenia and to evaluate the potential for infectious disease spread through movements. The study considered cattle movements between different types of premises that occurred between August 1, 2011 and July 31, 2016 with the exclusion of the movements to the end nodes (e.g., slaughterhouses). In the first part, we performed a static network analysis on monthly and yearly snapshots of the network. These time scales reflect our interest in slowly spreading pathogens; namely Mycobacterium avium subsp. paratuberculosis (MAP), which causes paratuberculosis, a worldwide economically important disease. The results showed consistency in the network measures over time; nevertheless, it was evident that year to year contacts between premises were changing. The importance of individual premises for the network connectedness was highly heterogeneous and the most influential premises in the network were collection centers, mountain pastures, and pastures. Compared to random node removal, targeted removal informed by ranking based on local network measures from previous years was substantially more effective in network disassociation. Inclusion of the latest movement data improved the results. In the second part, we simulated disease spread using a Susceptible-Infectious (SI) model on the temporal network. The SI model was based on the empirically estimated true prevalence of paratuberculosis in Slovenia and four scenarios for probabilities of transmission. Different probabilities were realized by the generation of new networks with the corresponding proportion of contacts which were randomly selected from the original network. These diluted networks served as substrates for simulation of MAP spread. The probability of transmission had a significant influence on the velocity of disease spread through the network. The peaks in daily incidence rates of infected herds were observed at the end of the grazing period. Our results suggest that network analysis may provide support in the optimization of paratuberculosis surveillance and intervention in Slovenia. The approach of simulating disease spread on a diluted network may also be used to model other transmission pathways between herds.

Reviewing age-structured epidemiological models of cattle diseases tailored to support management decisions: Guidance for the future

Mechanistic simulation models are being increasingly used as tools to assist with animal health decision-making in the cattle sector. We reviewed scientific literature for studies reporting age-structured cattle management models in application to infectious diseases. Our emphasis was on papers dedicated to support decision making in the field. In this systematic review we considered 1290 manuscripts and identified 76 eligible studies. These are based on 52 individual models from 10 countries addressing 9 different pathogens. We provide an overview of these models and present in detail their theoretical foundations, design paradigms and incorporated processes. We propose a structure of the characteristics of cattle disease models using three main features: [1] biological processes, [2] farming-related processes and [3] pathogen-related processes. It would be of benefit if future cattle disease models were to follow this structure to facilitate science communication and to allow increased model transparency.

Mastitis risk effect on the economic consequences of paratuberculosis control in dairy cattle: A stochastic modeling study

The benefits and efficacy of control programs for herds infected with Mycobacterium avium subsp. paratuberculosis (MAP) have been investigated under various contexts. However, most previous research investigated paratuberculosis control programs in isolation, without modeling the potential association with other dairy diseases. This paper evaluated the benefits of MAP control programs when the herd is also affected by mastitis, a common disease causing the largest losses in dairy production. The effect of typically suggested MAP controls were estimated under the assumption that MAP infection increased the rate of clinical mastitis. We evaluated one hundred twenty three control strategies comprising various combinations of testing, culling, and hygiene, and found that the association of paratuberculosis with mastitis alters the ranking of specific MAP control programs, but only slightly alters the cost-benefit difference of particular MAP control components, as measured by the distribution of net present value of a representative U.S. dairy operation. In particular, although testing and culling for MAP resulted in a reduction in MAP incidence, that control led to lower net present value (NPV) per cow. When testing was used, ELISA was more economically beneficial than alternative testing regimes, especially if mastitis was explicitly modeled as more likely in MAP-infected animals, but ELISA testing was only significantly associated with higher NPV if mastitis was not included in the model at all. Additional hygiene was associated with a lower NPV per cow, although it lowered MAP prevalence. Overall, the addition of an increased risk of mastitis in MAP-infected animals did not change model recommendations as much as failing to consider.

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

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.

Control measures to prevent the increase of paratuberculosis prevalence in dairy cattle herds: an individual-based modelling approach

Paratuberculosis, a gastrointestinal disease caused by Mycobacterium avium subsp. paratuberculosis (Map), can lead to severe economic losses in dairy cattle farms. Current measures are aimed at controlling prevalence in infected herds, but are not fully effective. Our objective was to determine the most effective control measures to prevent an increase in adult prevalence in infected herds. We developed a new individual-based model coupling population and infection dynamics. Animals are characterized by their age (6 groups) and health state (6 states). The model accounted for all transmission routes and two control measures used in the field, namely reduced calf exposure to adult faeces and test-and-cull. We defined three herd statuses (low, moderate, and high) based on realistic prevalence ranges observed in French dairy cattle herds. We showed that the most relevant control measures depend on prevalence. Calf management and test-and-cull both were required to maximize the probability of stabilizing herd status. A reduced calf exposure was confirmed to be the most influential measure, followed by test frequency and the proportion of infected animals that were detected and culled. Culling of detected high shedders could be delayed for up to 3 months without impacting prevalence. Management of low prevalence herds is a priority since the probability of status stabilization is high after implementing prioritized measures. On the contrary, an increase in prevalence was particularly difficult to prevent in moderate prevalence herds, and was only feasible in high prevalence herds if the level of control was high.

Evaluating the impact of transmission mode, calibration level and farmer compliance in simulation models of paratuberculosis in dairy herds

Simulation models can predict the outcome of different strategies for the control and eradication of paratuberculosis (PTB) in dairy herds. Two main transmission modes have previously been used to simulate the spread of PTB: direct (contact between animals) and indirect (through the environment). In addition, previous models were calibrated to either low or high within-herd prevalence levels, which we refer to as normal and low hygiene levels, respectively. We simulated both direct and indirect transmission with the same model in both normal and low hygiene level scenarios. The effectiveness of a test-and-cull strategy was dependent on the calibration level of the simulation model, and eradication occurred less frequently with the more biologically plausible indirect transmission mode. The results were compared to within-herd prevalence records from 314 dairy herds. The prevalence in 50% of the herds varied less than 0.9% per year on average, and less than 4% in 90% of the herds. We therefore conclude that the normal-hygiene scenario best describes most dairy herds in Denmark. Finally, we simulated different levels of farmer compliance with a test-and-cull strategy and found that a 60% compliance level was not sufficient to reach eradication within 10 years.

An agent-based model evaluation of economic control strategies for paratuberculosis in a dairy herd

This paper uses an agent-based simulation model to estimate the costs associated with Mycobacterium avium ssp. paratuberculosis (MAP), or Johne's disease, in a milking herd, and to determine the net benefits of implementing various control strategies. The net present value (NPV) of a 1,000-cow milking herd is calculated over 20 yr, parametrized to a representative US commercial herd. The revenues of the herd are generated from sales of milk and culled animals. The costs include all variable and fixed costs necessary to operate a representative 1,000-cow milking herd. We estimate the NPV of the herd with no MAP infection, under an expected endemic infection distribution with no controls, and under an expected endemic infection distribution with various controls. The initial number of cows in a herd with an endemic MAP infection is distributed as 75% susceptible, 13% latent, 9% low MAP shedding, and 3% high MAP shedding. Control strategies include testing using ELISA and fecal culture tests and culling of cows that test positive, and culling based on observable milk production decrease. Results show that culling cows based on test results does not increase the herd's NPV and in most cases decreases NPV due to test costs as well as false positives and negatives with their associated costs (e.g., culling healthy cows and keeping infected cows). Culling consistently low producing cows when MAP is believed to be present in the herd produces higher NPV over the strategy of testing and culling MAP infected animals, and over the case of no MAP control.

Evaluation of fecal shedding and antibody response in dairy cattle infected with paratuberculosis using national surveillance data in Japan

Paratuberculosis or Johne's disease (JD), is a chronic infectious disease causing intractable diarrhea in cattle, which leads to less productivity, such as decreased milk yield, and lower daily weight gain. As a control measure against JD in cattle, national serological surveillance has been conducted in Japan since 1998. To conduct modeling studies that are useful to evaluate the effectiveness of control measures against JD, reliable parameter values, such as length of time from infection to the start of fecal shedding or antibody expression, are especially important. These parameters in the Japanese cattle population are assumed to be different from those in other countries with a higher prevalence of JD or in experimental infection settings; therefore, they must be estimated for the cattle population in Japan. Data from national surveillance conducted in Tokachi District, Hokkaido Prefecture, were used for this study. Using data from JD diagnostic tests for all cattle in Tokachi District between 1998 and 2014, all testing histories for infected animals were estimated as the number of tested cattle and positive cattle at each age of month for both fecal and antibody tests. A deterministic mathematical model for JD development, from infection to fecal shedding and antibody expression in infected cattle, was constructed to obtain the probability of testing positive when applied to both fecal and antibody tests at a given age. Likelihood was obtained from these estimated test results and best values for parameters were obtained using the Markov Chain Monte-Carlo method. Fifty-five percent of infected cattle were projected to have a transient shedding period, which was estimated to start 12 months after infection and last for 4 months. Persistent shedding was projected to occur in all infected cattle, and estimated to begin 7-84 months from infection. Following persistent shedding, antibody expression was estimated to start 7 months later. These values are useful for developing models to evaluate the status of JD infection and the effectiveness of control measures in the Japanese cattle population.

Which phenotypic traits of resistance should be improved in cattle to control paratuberculosis dynamics in a dairy herd: a modelling approach

Paratuberculosis is a worldwide disease causing production losses in dairy cattle herds. Variability of cattle response to exposure to Mycobacterium avium subsp. paratuberculosis (Map) has been highlighted. Such individual variability could influence Map spread at larger scale. Cattle resistance to paratuberculosis has been shown to be heritable, suggesting genetic selection could enhance disease control. Our objective was to identify which phenotypic traits characterising the individual course of infection influence Map spread in a dairy cattle herd. We used a stochastic mechanistic model. Resistance consisted in the ability to prevent infection and the ability to cope with infection. We assessed the effect of varying (alone and combined) fourteen phenotypic traits characterising the infection course. We calculated four model outputs 25 years after Map introduction in a naïve herd: cumulative incidence, infection persistence, and prevalence of infected and affected animals. A cluster analysis identified influential phenotypes of cattle resistance. An ANOVA quantified the contribution of traits to model output variance. Four phenotypic traits strongly influenced Map spread: the decay in susceptibility with age (the most effective), the quantity of Map shed in faeces by high shedders, the incubation period duration, and the required infectious dose. Interactions contributed up to 12% of output variance, highlighting the expected added-value of improving several traits simultaneously. Combinations of the four most influential traits decreased incidence to less than one newly infected animal per year in most scenarios. Future genetic selection should aim at improving simultaneously the most influential traits to reduce Map spread in cattle populations.

Use of an Individual-based Model to Control Transmission Pathways of Mycobacterium avium Subsp. paratuberculosis Infection in Cattle Herds

Johne’s disease (JD) is a chronic enteric disease in cattle caused by Mycobacterium avian subsp. paratuberculosis (MAP). Eradicating JD is a difficult task due to the long incubation period of MAP, inefficient diagnostic tests, and delayed clinical signs. Effective control strategies can help farmers to reduce prevalence, but those most acceptable to farmers combine specific information about lactation performance and testing results, which existing models do not provide. This paper presents an individual-based model of MAP infection dynamics and assesses the relative performance of the applied alternative control strategies. The base dairy herd model included the daily life events of a dairy cow and reflects several current dairy management processes. We then integrated MAP infection dynamics into the model. The model adopted four different test-based control strategies based on risk-based culling decisions and three hygiene scenarios. The model tracked the source of each infection and quantified the efficacy of each control strategy in reducing the risks of different transmission routes. The results suggest that risk-based culling can reduce prevalence compared with no control, but cannot eliminate the infection. Overall, this work provides not only a valuable tool to investigate MAP transmission dynamics but also offers adaptability to model similar infectious diseases.

Economic consequences of paratuberculosis control in dairy cattle: A stochastic modeling study

The cost of paratuberculosis to dairy herds, through decreased milk production, early culling, and poor reproductive performance, has been well-studied. The benefit of control programs, however, has been debated. A recent stochastic compartmental model for paratuberculosis transmission in US dairy herds was modified to predict herd net present value (NPV) over 25 years in herds of 100 and 1000 dairy cattle with endemic paratuberculosis at initial prevalence of 10% and 20%. Control programs were designed by combining 5 tests (none, fecal culture, ELISA, PCR, or calf testing), 3 test-related culling strategies (all test-positive, high-positive, or repeated positive), 2 test frequencies (annual and biannual), 3 hygiene levels (standard, moderate, or improved), and 2 cessation decisions (testing ceased after 5 negative whole-herd tests or testing continued). Stochastic dominance was determined for each herd scenario; no control program was fully dominant for maximizing herd NPV in any scenario. Use of the ELISA test was generally preferred in all scenarios, but no paratuberculosis control was highly preferred for the small herd with 10% initial prevalence and was frequently preferred in other herd scenarios. Based on their effect on paratuberculosis alone, hygiene improvements were not found to be as cost-effective as test-and-cull strategies in most circumstances. Global sensitivity analysis found that economic parameters, such as the price of milk, had more influence on NPV than control program-related parameters. We conclude that paratuberculosis control can be cost effective, and multiple control programs can be applied for equivalent economic results.