- Research Article
- Open Access
Which phenotypic traits of resistance should be improved in cattle to control paratuberculosis dynamics in a dairy herd: a modelling approach
© The Author(s) 2017
Received: 24 February 2017
Accepted: 22 August 2017
Published: 10 October 2017
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.
Bovine paratuberculosis or Johne’s disease (JD) is a bacterial infection caused by Mycobacterium avium subsp. paratuberculosis (Map). It mainly affects domestic ruminants. Paratuberculosis has a worldwide distribution with a high prevalence, herd prevalence being around 50% in Europe . The progressive evolution of the infection leads to a chronic diarrhoea, an emaciation and death. This infection is responsible for significant weight losses, a decrease in milk production, an increase in mortality, and the early culling of infected animals, inducing economic losses . Infectious animals shed bacteria in their faeces, milk, and colostrum. Susceptible animals are infected by ingesting Map or in utero. Calves are known to be the age group most susceptible to infection .
Individual response to a given exposure to Map differs among animals. Within-herd prevalence is usually low, with 2.8–27% of infected animals [4, 5]. Field observations have reported substantial variation in individual response to Map exposure: among birth cohorts which are assumed to have been similarly exposed to Map, some are later shown to be infected/infectious, while others remain not infected/infectious. In addition, the following observations have been made following experimental infection of similar aged calves with similar infectious dose of Map: (1) a wide range of paratuberculosis lesion severity have subsequently been observed , (2) different quantities of Map are shed in their faeces , and (3) different antibody responses have been detected, suggesting a variable duration of the latency period (being the period between infection and later detection by direct or indirect tests) . The duration of the incubation period (which is defined as the period between infection and clinical signs) varied greatly between animals, ranging from 4 months to 15 years [8–12]. The amount of bacteria shed by infectious cattle is also highly variable, some being high shedders, while others are low shedders. Both intermittent and continuous shedding has been observed.
Individual resistance to paratuberculosis is assumed to be highly variable among, and expresses as different courses of infection. The phenotype of cattle resistance to paratuberculosis can be divided into (1) the ability to prevent infection and (2) the ability to cope with infection. This resistance in response to Map exposure involves different mechanisms and individual characters. Each of these characters will be denoted thereafter as phenotypic traits, a phenotype being defined by combined phenotypic traits. At the population scale, the distribution of phenotypic traits among individuals will influence the level of herd immunity, and therefore impair Map spread.
Strategies to control Map spread within dairy cattle herds usually consist in two main actions: hygiene improvement to reduce environmental and food contamination by Map, and a test-and-cull strategy to identify and remove infected animals. These control measures are not sufficient to control Map spread at herd and regional scales [13–15]. Vaccines against paratuberculosis have also been developed. Available vaccines decrease shedding of Map by infectious animals and decrease clinical signs of the disease [14, 16]. However, they do not prevent the infection of susceptible animals. In addition, most licensed vaccines show a cross reaction with tuberculosis diagnostic tests . Therefore, the use of vaccination is restricted in many countries.
The observed variability of the individual response to Map exposure could support the development of innovative control measures applied at population scale if the most resistant animals can be selected. Several studies demonstrated a heritability of resistance to paratuberculosis in cattle ranging from 0.01 to 0.23 [17–21]. Recent studies highlighted an association between genetic markers and the course of Map infection [22–26]. Other genome markers were associated with Map shedding in faeces, presence of Map in several tissues, and seropositivity, in animals from comparable herds regarding paratuberculosis infection and of the same age group. Therefore, these animals were assumed to have been exposed in a similar way [25, 27, 28]. This highlights the potential to select for cattle more resistant to paratuberculosis. However, there are still gaps of knowledge concerning the phenotypic traits of resistance that would be the most relevant to improvements in the control of Map spread at population scale.
Modelling is the most appropriate approach to investigate the dynamics of complex systems such as within-herd Map transmission. Observational and experimental studies are both difficult to implement and expensive regarding the long evolution of paratuberculosis. In addition, a modelling approach allows us to overpass the lack of knowledge on genetic resistance of cattle to paratuberculosis by assuming improved phenotypic traits as if they were already selected for. Simulations then provide information on how such modifications of phenotypic traits would influence Map spread. Only one recent study investigated the potential effectiveness of hypothetical genetic selection as a strategy to control paratuberculosis at herd scale . The authors assessed the effect of varying three phenotypic traits of resistance: (1) length of the susceptibility period, (2) level of the susceptibility to infection (expressed as the dose of Map required resulting in infection), and (3) duration of the latency period. Each tested phenotypic trait has been tested one-at-a-time and ranked by the time required to reach eradication. Modelling predictions showed that, when only genetic selection is implemented, eradication takes hundreds of years. However, this study did not investigate the potential progress in disease control when combining variations in several traits. In addition, other traits also could influence Map spread including intensity of shedding by infectious animals, in utero transmission, and progress of the infection course through different infection stages.
Our objective was to identify which phenotypic traits of resistance to paratuberculosis have the strongest influence on Map spread within a dairy cattle herd. The purpose was to identify ranges of phenotypic trait variations and trait combinations that limit Map spread in the herd. We assessed three categories of phenotypic traits both one-at-a-time and in combination, including: infection susceptibility, delays in the infection course, and shedding levels.
Materials and methods
Overall study design and model choice
A modelling approach was used to predict the effect of varying phenotypic traits of resistance to paratuberculosis on Map spread in a dairy cattle herd. We compared a situation where phenotypic traits were set at current observed levels with situations reached after a successful hypothetical genetic selection of more resistant animals in response to Map exposure. For each change of a trait, the resistance level was simulated as constant over time assuming that this average level had been reached in the population after a (not modelled) selection period. Several scenarios were simulated where one or several phenotypic traits were varied. The scenarios were compared regarding Map spread in the herd.
Several models have been published that represent Map spread within a dairy cattle herd (reviewed in , and more recently [29, 31–38]). We selected a stochastic compartmental model that offers an up-to-date description of Map spread within a dairy cattle herd. This model takes into account all of the major processes involved (according to the most recent literature) and allowed us to represent phenotypic traits of resistance corresponding to all of our hypotheses of interest. This model adequately combines demographical and infection dynamics, and accounts for herd structure, all these processes having been shown to highly influence Map spread . The chosen model is mechanistic: each step and mechanism of the infection course is represented by a model parameter. This allowed us to simulate changes in phenotypic traits of resistance by minimal changes in the model.
Main features of the model
The within-herd transmission model and the corresponding equations are fully described in Marcé et al.  and Beaunée et al. . And a detailed description of the model is presented in Additional file 1.
Susceptible animals can be infected through five transmission routes : (1) contact with bacteria present in the general environment of the farm contaminated by all shedders, (2) contact with bacteria present in the local environment of calves contaminated by shedding calves, (3) in utero transmission from infected cows to their foetus, and (4) ingestion of contaminated milk or (5) colostrum from infectious cows.
Phenotypic traits of resistance to paratuberculosis
In this study we assessed the effect of varying 14 phenotypic traits of resistances to paratuberculosis and combination of them on Map transmission in the herd. Each of the tested scenarios corresponds to a variation, or combinations in variation, in one or more phenotypic trait of resistance to paratuberculosis.
The phenotype of cattle resistance to paratuberculosis is classically divided into (1) resistance to infection defined as the ability to prevent infection when exposed to a given dose of Map, and (2) tolerance to infection defined as the ability to cope with infection when infected [40, 41]. On the one hand, animals are considered to be resistant (ability to prevent infection) if they show a decrease in susceptibility to infection if they are no longer susceptible at a younger age, if they need to be exposed to a higher dose of Map to be infected, or if they show a faster decrease in susceptibility with age than less resistant animals. On the other hand, animals are considered tolerant (ability to cope with infection once infected) if they show longer latency and incubation periods, and a lower shedding level when infectious than less tolerant animals. In addition, foetuses of the latter may have a lower chance to be infected in utero.
Parameters coding for the phenotypic traits of resistance: definition and values
Univariate simulations: [min–max]
Multivariate simulations: tested values
Susceptibility period duration
Decay in susceptibility with age (coefficient)
Required infectious dose to be infected
1.5 × 106
2 × 106
2.5 × 106
3 × 106
Duration of latent state
ν L + ν Is
Duration before high shedding and clinically affected state
v T with ν T + ν L = constant
Duration of transiently infectious state with constant duration before moderate shedding state
νT = 25 weeks
(νT + νL = 77)
ν Is with ν L + ν Is = constant
Duration of moderate shedding state with constant duration before high shedding or clinically affected state
νIs = 104 weeks
(νL + νIs = 156 weeks )
Factor of decrease of Map shed in milk by animals in health state X
Moderate shedding state (Is)
High shedding and clinically affected state (Ic)
Factor of decrease of Map shed in faeces by animals in health state X
Transient state (T)
Moderate shedding state (Is)
High shedding or clinically affected state (Ic)
Factor of decrease of probability of in utero transmission for cows in health state X
Latent and moderate shedding states (LIs)
High shedding or clinically affected state (Ic)
Based on the literature, we defined a realistic variation of resistance levels to simulate within observed values for the investigated traits. The reference value was the worst one. Changes were simulated from reference to the most favourable value observed value, Indeed, calves susceptibility can sharply reduce, and animals are no longer susceptible, as soon as their first week of life [53–56]. Some susceptible animals have been shown to need a dose of bacteria as high as 1012 to become infected . After a transient shedding period, infected animals can have a barely detectable level of shedding for about 4 years (208 weeks) [10, 11]. Infected animals can show clinical signs of the disease up to more than 9 years after infection (468 weeks). Concerning the probability of transmission of Map in utero from infected dam to its foetus and the quantities of bacteria shed through different routes, only partial information was available. Hence, we chose to test for extreme values by assuming that animals can stop shedding completely with no further in utero transmission of the infection. Nevertheless, it has been shown that high shedders and clinically affected animals can shed as few as 108 bacteria/kg of faeces, which corresponds to 1/100th of the reference value that we have assumed in our model [58, 59].
Initial conditions and model outputs
Map spread was initiated by the introduction of a moderate shedding cow into a fully naïve herd of 260 animals. We assumed that herd renewal is mainly driven by internal demographic processes (no further introduction), which is typical of western Europe farming systems. Map spread was predicted over 25 years. To obtain accurate outputs from the stochastic model, we ran 500 repetitions for each of the tested scenario. A scenario represented one phenotype of interest. Each phenotype was defined by a set of values of 14 parameters.
Contribution of the four most influential phenotypic traits to the model output variance
Prevalence of infected animals
Prevalence of affected animals
νL + νIs
h:νL + νIs
νL + νIs:φFecesIc
α:νL + νIs
h:α:νL + νIs
2 × 10−3
h:νL + νIs:φFecesIc
6 × 10−3
2 × 10−4
5 × 10−3
α:νL + νIs:φFecesIc
8 × 10−3
3 × 10−4
2 × 10−3
4 × 10−3
h:α:νL + νIs:φFecesIc
3 × 10−3
3 × 10−3
4 × 10−4
7 × 10−5
Simulation protocol and output analysis
First, we performed a univariate simulation study: each of the traits of interest was varied one-at-a-time, assuming they varied independently (Table 1). Second, we performed a multivariate simulation study: combinations of phenotypic traits were studied to test for a potential enhanced effect of simultaneously improving several phenotypic traits simultaneously. The R programming language  was used for data analyses. Results obtained in the univariate simulation study revealed that some parameters—when analysed one-at-a-time—did not influence model outputs. Instead of keeping numerous parameters or removing some of them expected not to be influential, we grouped in this second step non-influential parameters when they are untangled in the same trait or when involved in a given transmission route. This decrease in the number of considered parameters without losing information eased the interpretation of the multivariate simulation study results. Ranges of variation of phenotypic traits were represented by five possible values per trait (including the reference value) combined in the multivariate simulation study using a complete factorial design, leading to 390 625 scenarios (Table 1). Five levels of variation per trait appeared to be a good compromise between parameter space exploration and number of scenarios to investigate interactions. A complete factorial design was required to assess all interaction orders.
In order to identify the most effective combinations of variation in phenotypic traits to decrease Map spread, we used the cumulative incidence output as an indicator of a successful Map control at herd scale. For each combination, the cumulative incidence was plotted and visually described (Figure 6). In addition, we chose two thresholds to evidence the most effective combinations with emphasis on the ones with the lowest variations in parameters (Figures 6, 7): (1) 25 newly infected animals over the 25 years of simulation, interpreting such a level of one newly infected animal per year as an infection under control, and (2) half this threshold, i.e. 12 newly infected animals over the 25 years of simulation.
Eight of the traits investigated in the univariate simulation study did not influence Map spread dynamics. These traits were: a shorter transiently infectious state when assuming a constant duration before the moderate shedding state (v T with v T + v L = constant), a shorter latent state when assuming a constant duration before the high shedding and clinically affected state (v Is with v L + v Is = constant), a decrease in the quantity of Map shed in milk by moderate shedders (φMilk Is ) and high shedders and clinically affected animals (φMilk Ic ), a decrease in the quantity of Map shed by transiently infectious animals (φFeces T ), and by moderate shedders (φFeces Is ), and a decrease in the probability of in utero transmission by latent infectious animals and moderate shedders (φP LIs ), and by high shedders and clinically affected animals (φP Ic ).
We chose traits to be included in the multivariate simulation study in light of these results, noting that it was not possible to evaluate interactions among traits could have been evaluated with such a univariate analysis. Among traits highlighted as influential, we kept all except u that was redundant with h. Among other traits, we kept v Is (assuming v L + v Is = constant) and we grouped traits related to Map shedding in milk and colostrum (φMilk), to in utero transmission (φP), and to Map shedding in faeces by transiently infectious animals and moderate shedders (φFeces TIs ).
Combined variations of phenotypic traits of resistance contributed to decrease Map spread dynamics. Interactions among traits showed contributions to model output variance ranging from 0.007% to up to 12% (Table 2). The interaction between increased decay in susceptibility with age (h) and a lengthened incubation period (ν L + ν Is ) contributed to 12% of the variance of the prevalence of affected animals, and was also the most contributing interaction for other model outputs. In addition, h was involved in all of the contributing interactions, thus having both the highest principal and interaction effects.
Variations of 4 of 14 phenotypic traits strongly reduced Map spread within a dairy cattle herd: the decay in susceptibility with age, this being the most influential trait, the quantity of Map shed in faeces by high shedders and clinically affected animals, the duration of the incubation period, and the required infectious dose. Combining these phenotypic traits was the sole way to effectively control Map spread at the herd scale. Most tested combinations of these influential phenotypic traits allowed the cumulative incidence to be reduced to <25 newly infected animals over the 25 years of simulation, which was interpreted here as an infection under control. Interestingly, such a low level of cumulative incidence could not be reached when varying a single phenotypic trait.
The increase in the decay in susceptibility with age is largely related to a shorter susceptibility period. We also highlighted the required infectious dose as an influential phenotypic trait. Our results concerning these traits are in agreement with van Hulzen et al. , who in a theoretical study also identified that an earlier resistance acquisition would be crucial when it comes to control paratuberculosis using genetic selection. However, there is nowadays no available knowledge to implement a genetic selection on these traits. These traits are not easily measurable in field conditions.
A decrease in the quantity of Map shed in faeces by high shedders and clinically affected animals, which was also identified as an influential phenotypic trait, might be achieved thanks to genetic selection. Currently, it has been shown that genetic markers could be associated with the occurrence of shedding versus no shedding at all by animals in infected herds [23, 27, 47]. More precise knowledge is needed concerning our ability to select cattle that will shed less Map in faeces while in their last stage of infection.
While van Hulzen et al.  identified the increase in duration of the latency period as an effective phenotypic trait in controlling paratuberculosis through genetic selection, we highlighted that an increase in this latency period (this being the period between infection and the occurrence of a moderate detectable shedding) without delaying the start of the high shedding or clinically affected state did not influence Map spread dynamics in the herd. We have shown that it will be more interesting to lengthen the incubation period, as this delays the occurrence of the high shedding or clinically affected state.
Phenotypic traits identified as influencing Map spread dynamics at the herd scale also are related to control measures currently implemented in infected herds in the field . Therefore, a valuable interaction can be expected between routine control plans and innovative control through genetic selection.
The variation of several other traits did not influence Map spread dynamics: decrease in duration of transiently infectious state with a constant duration before moderate shedding state, decrease in quantity of Map shed in milk and colostrum irrespective of the animal infection state, decrease in quantity of Map shed in faeces by transiently infectious animals and moderate shedders, and decrease in probability of in utero transmission irrespective of animal infection state. A decrease in duration of moderate shedding state (from 104 to 60 weeks) did not influence Map spread dynamics. The range of variation modelled for this trait was lower than for other traits due to limitations inherent to the compartmental model. Nevertheless, as no effect was evidenced with a reduction of one-third of that duration, we assumed that this trait was not highly influential over the simulated range.
The four traits identified as influential are well described in the literature therefore we can assume that their tested ranges of variation were realistic. We assumed extreme ranges of variation for traits for which information was missing. The other traits assessed were not influential even with such extreme, non-realistic, variations. Using a different set of variation in the investigated traits is not expected to change our conclusions concerning which traits influence Map spread within dairy cattle herds.
As our objective was to assess Map spread in herds in which phenotypic traits would have been improved, we did not account for the long time needed  to reach such targeted levels of phenotypic traits by a potential genetic selection. On the one hand, recent studies identified several genetic markers associated with resistance (reviewed in [48, 49]), but genes and mechanisms responsible for the tested phenotypic traits are still unknown. Further genetic studies of resistance of cattle to paratuberculosis are required, especially to identify genes and mechanisms involved in these relevant phenotypic traits to allow potential future selection of more resistant cattle. In addition, diagnostic tests currently available in the field do not allow identifying animals having the phenotypic traits identified here. Concerning future genetic selection, tests more sensitive during the early stage of the infection would be needed to distinguish infected animals from others and to better quantify the individual duration of incubation periods. On the other hand, there is a risk of a negative association between phenotypic traits of resistance to paratuberculosis and other traits of economic importance. For example, it has been shown that genetic markers associated with susceptibility to paratuberculosis could be associated to lactation persistence .
Our model represents a typical western Europe farming system for dairy cattle herds. Demographic processes have been shown to highly influence the disease dynamics , and therefore, different farming systems could change the influence of the studied phenotypic traits on Map spread dynamics in the herd. We assumed a closed herd without introduction of animals from other herds. Animal exchanges between herds could reintroduce Map in free herds and thus influence disease dynamics. However, it is not expected to modify our conclusions as regards the identification of crucial phenotypic traits to better manage infected herds. Indeed, a single Map introduction can lead to infection persistence in 40% of the cases under current situation as regards phenotypic traits , with a huge cumulative incidence reached after 25 years if no control is applied. Animal movements are not expected to modify significantly this finding. However, the occurrence of animal movements might increase the cumulative incidence under controlled situations with improved traits.
This study highlighted four phenotypic traits of resistance of cattle to paratuberculosis influencing Map spread within a dairy herd: decay in susceptibility with age, quantity of Map shed in faeces by high shedders and clinically affected animals, duration of the incubation period, and required infectious dose. A combination of these traits strongly contributes to limit Map spread. Further genetic study should aim at better identifying cattle genes involved in these traits in order to allow their potential selection.
The authors declare that they have no competing interests.
Conceived and designed the study: RBR, CF, PE, and RG. Adapted the model: GC, RBR, GB, and PE. Performed the simulations: RBR. Performed statistical analyses of simulated data: RBR and PE. Interpreted the results and wrote the paper: RBR, CF, and PE. All authors read and approved the final manuscript.
The authors acknowledge the funding from the INRA metaprogram “Integrated Management of Animal Health” for the “GISA-PICSAR” project, from the Regional Council of Pays de la Loire, from APIS-GENE and GDS France for the “PARADIGM” project, with additional support (computational resources and expertise) from The French Research Agency, program investments for the future, Grant #ANR-10-BINF-07 (MIHMES), and the European fund for the regional development of Pays de la Loire. The funders had no role in the production, analysis, interpretation and writing of the data. Alain Joly (Animal Health Services of Brittany), and Elisabeta Vergu and Sandie Arnoux (INRA) are thanked for their helpful comments on model analyses and output interpretation. Pr. Simon MORE is thanked for his help in improving the language editing of the manuscript.
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.
- Nielsen SS, Toft N (2009) A review of prevalences of paratuberculosis in farmed animals in Europe. Prev Vet Med 88:1–14View ArticlePubMedGoogle Scholar
- Whittington RJ, Windsor PA (2009) In utero infection of cattle with Mycobacterium avium subsp. paratuberculosis: a critical review and meta-analysis. Vet J 179:60–69View ArticlePubMedGoogle Scholar
- Mortier R, Orsel K, Barkema HW, Atkins G, De Buck J (2011) Age and dose dependent susceptibility to Mycobacterium avium subsp. paratuberculosis in dairy cattle. WCDS Adv Dairy Technol 23:375Google Scholar
- Good M, Clegg T, Sheridan H, Yearsely D, O’Brien T, Egan J, Mullowey P (2009) Prevalence and distribution of paratuberculosis (Johne’s disease) in cattle herds in Ireland. Ir Vet J 62:597–606View ArticlePubMedPubMed CentralGoogle Scholar
- Raizman EA, Wells SJ, Muñoz-Zanzi CA, Tavornpanich S (2011) Estimated within-herd prevalence (WHP) of Mycobacterium avium subsp. paratuberculosis in a sample of Minnesota dairy herds using bacterial culture of pooled fecal samples. Can J Vet Res 75:112–116PubMedPubMed CentralGoogle Scholar
- Mortier RA, Barkema HW, Orsel K, Wolf R, De Buck J (2014) Shedding patterns of dairy calves experimentally infected with Mycobacterium avium subspecies paratuberculosis. Vet Res 45:71View ArticlePubMedPubMed CentralGoogle Scholar
- Mortier RA, Barkema HW, De Buck J (2015) Susceptibility to and diagnosis of Mycobacterium avium subspecies paratuberculosis infection in dairy calves: a review. Prev Vet Med 121:189–198View ArticlePubMedGoogle Scholar
- Stewart DJ, Vaughan JA, Stiles PL, Noske PJ, Tizard ML, Prowse SJ, Michalski WP, Butler KL, Jones SL (2007) A long-term bacteriological and immunological study in Holstein-Friesian cattle experimentally infected with Mycobacterium avium subsp. paratuberculosis and necropsy culture results for Holstein-Friesian cattle, Merino sheep and Angora goats. Vet Microbiol 122:83–96View ArticlePubMedGoogle Scholar
- van Roermund HJW, Bakker D, Willemsen PTJ, de Jong MC (2007) Horizontal transmission of Mycobacterium avium subsp. paratuberculosis in cattle in an experimental setting: calves can transmit the infection to other calves. Vet Microbiol 122:270–279View ArticlePubMedGoogle Scholar
- Mitchell RM, Medley GF, Collins MT, Schukken YH (2012) A meta-analysis of the effect of dose and age at exposure on shedding of Mycobacterium avium subspecies paratuberculosis (MAP) in experimentally infected calves and cows. Epidemiol Infect 140:231–246View ArticlePubMedGoogle Scholar
- Nielsen SS (2008) Transitions in diagnostic tests used for detection of Mycobacterium avium subsp. paratuberculosis infections in cattle. Vet Microbiol 132:274–282View ArticlePubMedGoogle Scholar
- Matthews HT (1947) On Johne’s disease. Vet Rec 59:397–401PubMedGoogle Scholar
- Ezanno P, van Schaik G, Weber MF, Heesterbeek JA (2005) A modeling study on the sustainability of a certification-and-monitoring program for paratuberculosis in cattle. Vet Res 36:811–826View ArticlePubMedGoogle Scholar
- Bastida F, Juste RA (2011) Paratuberculosis control: a review with a focus on vaccination. J Immune Based Ther Vaccines 9:8View ArticlePubMedPubMed CentralGoogle Scholar
- Beaunée G, Vergu E, Ezanno P (2015) Modelling of paratuberculosis spread between dairy cattle farms at a regional scale. Vet Res 46:111View ArticlePubMedPubMed CentralGoogle Scholar
- Kalis CH, Hesselink JW, Barkema HW, Collins MT (2001) Use of long-term vaccination with a killed vaccine to prevent fecal shedding of Mycobacterium avium subsp. paratuberculosis in dairy herds. Am J Vet Res 62:270–274View ArticlePubMedGoogle Scholar
- Behr MA, Collins DM (2010) Paratuberculosis: organism, disease. Control, vol 1. CABI, Wallingford, UK, Cambridge, USA, p 1–388Google Scholar
- Kirkpatrick BW, Shook GE (2011) Genetic susceptibility to paratuberculosis. Vet Clin N Am Food Anim Pract 27:559–571View ArticleGoogle Scholar
- van Hulzen KJ, Nielen M, Koets AP, de Jong G, van Arendonk JA, Heuven HC (2011) Effect of herd prevalence on heritability estimates of antibody response to Mycobacterium avium subspecies paratuberculosis. J Dairy Sci 94:992–997View ArticlePubMedGoogle Scholar
- Küpper J, Brandt H, Donat K, Erhardt G (2012) Heritability estimates for Mycobacterium avium subspecies paratuberculosis status of German Holstein cows tested by fecal culture. J Dairy Sci 95:2734–2739View ArticlePubMedGoogle Scholar
- Zare Y, Shook GE, Collins MT, Kirkpatrick BW (2014) Short communication: heritability estimates for susceptibility to Mycobacterium avium subspecies paratuberculosis infection defined by ELISA and fecal culture test results in Jersey cattle. J Dairy Sci 97:4562–4567View ArticlePubMedGoogle Scholar
- Purdie AC, Plain KM, Begg DJ, de Silva K, Whittington RJ (2011) Candidate gene and genome-wide association studies of Mycobacterium avium subsp. paratuberculosis infection in cattle and sheep: a review. Comp Immunol Microbiol Infect Dis 34:197–208View ArticlePubMedGoogle Scholar
- Kirkpatrick BW, Shi X, Shook GE, Collins MT (2011) Whole-Genome association analysis of susceptibility to paratuberculosis in Holstein cattle. Anim Genet 42:149–160View ArticlePubMedGoogle Scholar
- van Hulzen KJ, Koets AP, Nielen M, Hoeboer J, van Arendonk JA, Heuven HC (2012) Genetic variation for infection status as determined by a specific antibody response against Mycobacterium avium subspecies paratuberculosis in milk of Dutch dairy goats. J Dairy Sci 95:6145–6151View ArticlePubMedGoogle Scholar
- Alpay F, Zare Y, Kamalludin MH, Huang X, Shi X, Shook GE, Collins MT, Kirkpatrick BW (2014) Genome-wide association study of susceptibility to infection by Mycobacterium avium subspecies paratuberculosis in Holstein cattle. PLoS One 9:e111704View ArticlePubMedPubMed CentralGoogle Scholar
- Zanella R, Settles ML, McKay SD, Schnabel R, Taylor J, Whitlock RH, Schukken Y, Van Kessel JS, Smith JM, Neibergs HL (2011) Identification of loci associated with tolerance to Johne’s disease in Holstein cattle. Anim Genet 42:28–38View ArticlePubMedGoogle Scholar
- Settles M, Zanella R, McKay SD, Schnabel RD, Taylor JF, Whitlock R, Schukken Y, Van Kessel JS, Smith JM, Neibergs H (2009) A whole genome association analysis identifies loci associated with Mycobacterium avium subsp. paratuberculosis infection status in US holstein cattle. Anim Genet 40:655–662View ArticlePubMedGoogle Scholar
- Neibergs HL, Settles ML, Whitlock RH, Taylor JF (2010) GSEA-SNP identifies genes associated with Johne’s disease in cattle. Mamm Genome 21:419–425View ArticlePubMedGoogle Scholar
- van Hulzen KJ, Koets AP, Nielen M, Heuven HC, van Arendonk JA, Klinkenberg D (2014) The effect of genetic selection for Johne’s disease resistance in dairy cattle: results of a genetic-epidemiological model. J Dairy Sci 97:1762–1773View ArticlePubMedGoogle Scholar
- Marcé C, Ezanno P, Weber MF, Seegers H, Pfeiffer DU, Fourichon C (2010) Invited review: modeling within-herd transmission of Mycobacterium avium subspecies paratuberculosis in dairy cattle: a review. J Dairy Sci 93:4455–4470View ArticlePubMedGoogle Scholar
- Lu Z, Schukken YH, Smith RL, Grohn YT (2010) Stochastic simulations of a multi-group compartmental model for Johne’s disease on US dairy herds with test-based culling intervention. J Theor Biol 264:1190–1201View ArticlePubMedGoogle Scholar
- Marcé C, Ezanno P, Seegers H, Pfeiffer DU, Fourichon C (2011) Predicting fadeout versus persistence of paratuberculosis in a dairy cattle herd for management and control purposes: a modelling study. Vet Res 42:36View ArticlePubMedPubMed CentralGoogle Scholar
- Marcé C, Ezanno P, Seegers H, Pfeiffer DU, Fourichon C (2011) Within-herd contact structure and transmission of Mycobacterium avium subspecies paratuberculosis in a persistently infected dairy cattle herd. Prev Vet Med 100:116–125View ArticlePubMedGoogle Scholar
- Robins J, Bogen S, Francis A, Westhoek A, Kanarek A, Lenhart S, Eda S (2015) Agent-based model for Johne’s disease dynamics in a dairy herd. Vet Res 46:68View ArticlePubMedPubMed CentralGoogle Scholar
- Martcheva M, Lenhart S, Eda S, Klinkenberg D, Momotani E, Stabel J (2015) An immuno-epidemiological model for Johne’s disease in cattle. Vet Res 46:69View ArticlePubMedPubMed CentralGoogle Scholar
- Smith RL, Schukken YH, Gröhn YT (2015) A new compartmental model of Mycobacterium avium subsp. paratuberculosis infection dynamics in cattle. Prev Vet Med 122:298–305View ArticlePubMedPubMed CentralGoogle Scholar
- Koets AP, Gröhn YT (2015) Within- and between-host mathematical modeling of Mycobacterium avium subspecies paratuberculosis (MAP) infections as a tool to study the dynamics of host–pathogen interactions in bovine paratuberculosis. Vet Res 46:60View ArticlePubMedPubMed CentralGoogle Scholar
- Al-Mamun MA, Smith RL, Schukken YH, Gröhn YT (2016) Modeling of Mycobacterium avium subsp. paratuberculosis dynamics in a dairy herd: an individual based approach. J Theor Biol 408:105–117View ArticlePubMedGoogle Scholar
- Pradhan AK, Mitchell RM, Kramer AJ, Zurakowski MJ, Fyock TL, Whitlock RH, Smith JM, Hovingh E, Van Kessel JA, Karns JS, Schukken YH (2011) Molecular epidemiology of Mycobacterium avium subsp. paratuberculosis in a longitudinal study of three dairy herds. J Clin Microbiol 49:893–901View ArticlePubMedPubMed CentralGoogle Scholar
- Råberg L, Sim D, Read AF (2007) Disentangling genetic variation for resistance and tolerance to infectious diseases in animals. Science 318:812–814View ArticlePubMedGoogle Scholar
- Schneider DS, Ayres JS (2008) Two ways to survive infection: what resistance and tolerance can teach us about treating infectious diseases. Nat Rev Immunol 8:889–895View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team (2016) R: A language and environment for statistical computing. Vienna, Austria. R Foundation for Statistical. http://www.R-project.org/
- Everitt BS, Hothorn T (2010) Cluster analysis: classifying Romano-British pottery and exoplanets. In: A handbook of statistical analyses using R, 2nd ed. CRC Press, Taylor and Francis Group edition, Boca Raton, p 315–348. http://www.bagualu.net/wordpress/wp-content/uploads/2015/10/A_Handbook_of_Statistical_Analyses_Using_R__Second_Edition.pdf
- Thinsungnoen T, Kaoungku N, Durongdumronchai P, Kerdprasop K, Kerdprasop N (2015) The clustering validity with silhouette and sum of squared errors. In: Proceedings of the 3rd international conference on industrial application engineering, p 44–51. https://www2.ia-engineers.org/iciae/index.php/iciae/iciae2015/paper/viewFile/576/380
- MacQueen J (1967) Some methods for classification and analysis of multivariate observations. In: Proceedings of the fifth berkeley symposium on mathematical statistics and probability, vol 1. University of California Press, Berkeley, California, p 281–297. https://projecteuclid.org/download/pdf_1/euclid.bsmsp/1200512992
- Husson F, Josse J, Le S, Mazet J (2016) FactoMineR: multivariate exploratory data analysis and data mining. R package version 1.32. https://CRAN.R-project.org/package=FactoMineR
- Zare Y, Shook GE, Collins MT, Kirkpatrick BW (2014) Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle. PLoS One 9:e88380View ArticlePubMedPubMed CentralGoogle Scholar
- McSpadden K, Caires K, Zanella R (2013) The effect of Mycobacterium avium subspecies paratuberculosis exposure on animal health. Acta Sci Vet 41:1095Google Scholar
- Pauciullo A, Küpper J, Brandt H, Donat K, Iannuzzi L, Erhardt G (2015) Wingless-type MMTV integration site family member 2 (WNT2) gene is associated with resistance to MAP in faecal culture and antibody response in Holstein cattle. Anim Genet 46:122–132View ArticlePubMedGoogle Scholar
- Carvajal AM, Huircan P, Lepori A (2013) Single nucleotide polymorphisms in immunity-related genes and their association with mastitis in Chilean dairy cattle. Genet Mol Res 12:2702–2711View ArticlePubMedGoogle Scholar
- Rankin JD (1962) The experimental infection of cattle with Mycobacterium johnei. IV. Adult cattle maintained in an infectious environment. J Comp Pathol 72:113–117View ArticlePubMedGoogle Scholar
- Hagan WA (1938) Age as a factor in susceptibility to Johne’s disease. Cornell Vet 28:34–40Google Scholar
- Whitlock RH, Buergelt C (1996) Preclinical and clinical manifestations of paratuberculosis (including pathology). Vet Clin N Am Food Anim Pract 12:345–356View ArticleGoogle Scholar
- Windsor PA, Whittington RJ (2010) Evidence for age susceptibility of cattle to Johne’s disease. Vet J 184:37–44View ArticlePubMedGoogle Scholar
- Begg DJ, Whittington RJ (2008) Experimental animal infection models for Johne’s disease, an infectious enteropathy caused by Mycobacterium avium subsp. paratuberculosis. Vet J 176:129–145View ArticlePubMedGoogle Scholar
- Sweeney RW, Whitlock RH, Rosenberger AE (1992) Mycobacterium paratuberculosis isolated from fetuses of infected cows not manifesting signs of the disease. Am J Vet Res 53:477–480PubMedGoogle Scholar
- Giese SB, Ahrens P (2000) Detection of Mycobacterium avium subsp. paratuberculosis in milk from clinically affected cows by PCR and culture. Vet Microbiol 77:291–297View ArticlePubMedGoogle Scholar
- Rossiter CA, Burhans WS (1996) Farm-specific approach to paratuberculosis (Johne’s disease) control. Vet Clin N Am Food Anim Pract 12:383–415View ArticleGoogle Scholar
- Whittington RJ, Reddacliff LA, Marsh I, McAllister S, Saunders V (2000) Temporal patterns and quantification of excretion of Mycobacterium avium subsp. paratuberculosis in sheep with Johne’s disease. Aust Vet J 78:34–37View ArticlePubMedGoogle Scholar
- Jørgensen JB (1982) An improved medium for culture of Mycobacterium paratuberculosis from bovine faeces. Acta Vet Scand 23:325–335PubMedGoogle Scholar
- Benedictus A, Mitchell RM, Linde-Widmann M, Sweeney R, Fyock T, Schukken YH, Whitlock RH (2008) Transmission parameters of Mycobacterium avium subspecies paratuberculosis infections in a dairy herd going through a control program. Prev Vet Med 83:215–227View ArticlePubMedGoogle Scholar