Effectiveness analysis of resistance and tolerance to infection
- Johann C Detilleux^{1}Email author
DOI: 10.1186/1297-9686-43-9
© Detilleux; licensee BioMed Central Ltd. 2011
Received: 13 October 2010
Accepted: 1 March 2011
Published: 1 March 2011
Abstract
Background
Tolerance and resistance provide animals with two distinct strategies to fight infectious pathogens and may exhibit different evolutionary dynamics. However, few studies have investigated these mechanisms in the case of animal diseases under commercial constraints.
Methods
The paper proposes a method to simultaneously describe (1) the dynamics of transmission of a contagious pathogen between animals, (2) the growth and death of the pathogen within infected hosts and (3) the effects on their performances. The effectiveness of increasing individual levels of tolerance and resistance is evaluated by the number of infected animals and the performance at the population level.
Results
The model is applied to a particular set of parameters and different combinations of values. Given these imputed values, it is shown that higher levels of individual tolerance should be more effective than increased levels of resistance in commercial populations. As a practical example, a method is proposed to measure levels of animal tolerance to bovine mastitis.
Conclusions
The model provides a general framework and some tools to maximize health and performances of a population under infection. Limits and assumptions of the model are clearly identified so it can be improved for different epidemiological settings.
Background
The breeding objective in most livestock species is to increase profit by improving performance efficiency. One way to reach this objective is to improve the animals' health, for example, through the implementation of appropriate management methods (e.g. chemotherapy, vaccination, and control of disease vectors). A more sustainable method consists in taking advantage, by selective breeding, of the within-breed variation that exists in the mechanisms of defenses against infectious pathogens [1]. Indeed, hosts have evolved resistance and tolerance defenses [2], thus breeders may choose, as progenitors, animals with the highest levels of resistance, tolerance, or both. One the one hand, resistance is the ability of the host to reduce the success of infection or to increase the rate of clearance of the pathogens. On the other hand, tolerance is the ability to reduce the detrimental effects of the pathogens on the performances of the hosts, either directly or by limiting immunopathological mechanisms [3]. The rate of transmission diminishes naturally among resistant hosts but not necessarily among tolerant ones, as these harbor the pathogen with no or moderate loss in performance [4].
Resistance and tolerance are associated with fitness costs, which arise from the diversion of limiting resources away from biological processes related to performance [5]. If these costs are too high, they may outweigh the effectiveness of the chosen strategy. Direct evidence of such costs can be found in experiments in insects [6], rainbow trout [7], crustaceans [8], wild birds [9] and mice [10].
To decide whether improving resistance, tolerance, neither, or both is the most effective strategy, it is proposed to (1) characterize the dynamics of the pathogens within and between hosts in the population under study, (2) evaluate the impact of the infection on the performances of the population, and (3) choose the most effective strategy. The goal of this study is to illustrate the methodology with a non-lethal micro-parasitic disease in a population where hosts have different levels of resistance to multiplication of the pathogen and different levels of tolerance to damages induced directly by the pathogens.
Methods
Pathogen dynamics
The model chosen here to depict the dynamics of transmission of the infection in a herd is a stochastic version of the SIS (S for susceptible, I for infected) model for the spread of a disease in a closed population of N individuals [11]. This model is appropriate for infections with no permanent immunity after recovery, i.e. individuals are susceptible to the infection, potentially get infected, may recover and become susceptible again. The time-scale of the disease process is assumed to be short compared to the life length of the host and no demographic turnover (natural birth or death) is considered. The area occupied by parasites and hosts is constant, so that numbers and densities coincide. There is only a single non-evolving pathogen species within infected hosts. Once infected, hosts are immediately able to infect other individuals (no latent period). Within the host, the number of pathogens increases following a sigmoidal growth curve and is directly related to the number of immune constituents of the host response to the pathogen, with no distinction between innate and specific immunity. Recovered hosts are as susceptible to infection as naïve hosts and re-exposure does not accelerate development of the disease.
In mathematical terms, the process is described by a continuous time Markov chain, {C^{i}_{t}; t = 0 to T, i = 1 to N}, where C^{i}_{t} denotes the number of pathogens in the i^{th} host at time t. Units of time are chosen arbitrarily. The chain has three transition probabilities (over a small time interval Δt) reflecting the three events, i.e. invasion of a new host by the pathogen, its multiplication and its killing by the immune response of the host.
where o(Δt) tends to 0 when Δt is small, C_{min} is the minimum number of pathogens necessary to have infection, β is the per-capita rate of successful transmission of C_{min} pathogens from an infectious host to a susceptible host upon contact with an infectious individual and during Δt, and I_{t} is the number of infectious hosts with which the i^{th} susceptible has contact.
where γ is the pathogen growth rate. Right after becoming infected, pathogen growth in a host is approximately exponential but it slows down as it reaches a maximum (C^{Max}), at which it stops.
This equation follows from the dynamics between pathogens and immune factors, as observed in experimental studies [12, 13]. Parameter μ represents the maximum number of pathogens killed for each unit of R^{i}_{t}, with R^{i}_{t} being a generic index related to the number, at time t, of the different types of immune factors specific for that pathogen. Because the main interest is on the number of pathogens, the complexity of the immune response is greatly simplified when R^{i}_{t} increases at a rate (h^{i}) that is constant across time. The scaling parameter ρ^{i} varies from 0 to 1 and represents the extra investment in resistance of the i^{th} host with respect to μ. When C^{i}_{t} drops below C_{min}, the infection is assumed to be cleared.
The Markov chain was simulated using the Gillespie algorithm [14], which essentially uses exponential waiting times between events. For all simulations, it was assumed that two individuals in the population were initially infected. Simulation steps were executed until t reaches T units of time (= one replicate) and repeated over 50 replicates. Each cycle took around 4 hours to complete, so the population size was limited at 30 individuals, which is the average size of most dairy herds in the Walloon region of Belgium.
Individual performance
where P^{i}_{t} is the performance of the i^{th} host at time t, when it is infected with C^{i}_{t} pathogens, ω is the maximum amount of performance lost per pathogen (virulence). The parameter λ^{i} is a scaling parameter representing the extra investment in tolerance. If λ^{i} = 1, the host is completely tolerant and produces at a level identical to the one without infection. If λ^{i} = 0, the host is not tolerant to the deleterious effects of the pathogen.
where P^{Max} is the totality of the resources available to the host to insure performance (e.g., production, reproduction, work) and to cope with an infection (resistance and tolerance). If no extra-investments are put in resistance and tolerance, all resources are allocated to insure the highest achievable level of performance in the absence of infection. Parameter c^{i}_{ρ} is the marginal cost of resistance and c^{i}_{λ} is the marginal cost of tolerance (in units of performance). Values for both costs are constrained such that the factor within brackets remains positive (ρ^{i} c^{i}_{ρ} + λ^{i} c^{i}_{λ} ≤ 1). A constraint was also set to insure P^{i}_{t} (equation 4) remains positive or null in totally non-tolerant individuals infected with K pathogens: ω ≤ P^{Max}(1 - ρ^{i} c^{i}_{ρ})/K.
Effectiveness analysis
The most profitable strategy, i.e. the one that will insure the lowest number of infected animals or the highest performance of the population, or both, was identified by weighing the allowed extra investments in resistance, tolerance, or both, against the effectiveness of each of these alternatives.
Effectiveness was computed by comparing populations under the same infection process but in which animals invest ('yes' population) or not ('no' populations) in resistance, tolerance, or both. To do so, the number of infected hosts (I_{t}) and the overall performance (P_{t} = Σ_{i = 1,N} P^{i}_{t}) were followed across time, and the area under the curves of P_{t} (AUC_{P}) and I_{t} (AUC_{I}) were obtained for t = 0 to T with the spline method of the procedure Expand of SAS^{®}[16]. Subsequently, the incremental effects (ΔE_{I} and ΔE_{P}) were computed as the difference between corresponding 'yes' and 'no' populations: ΔE_{I} = AUC_{I}^{no} - AUC_{I}^{yes}, and ΔE_{P} = AUC_{P}^{yes} - AUC_{P}^{no}. Then, the most effective alternative was identified as the one with the highest values for ΔE_{I} and ΔE_{P}.
Model parameters and their values
Symbol | Description | Values |
---|---|---|
P_{t}^{i} | Performance for animal i at time t | |
C_{t}^{i} | Pathogen number in animal i at time t | |
R_{t}^{i} | Immune response in animal i at time t | |
I_{t} | Number of infected animals in the population at time t | |
P_{t} | Population performance at time t | |
ΔE_{P} | Incremental effectiveness for P_{t} over the T period | |
ΔE_{I} | Incremental effectiveness for I_{t} over the T period | |
N | Population size | 30 |
T | Time duration | 300 |
C_{min} | Pathogen number necessary for infection | 10 |
C^{Max} | Maximum number of pathogens | 500 |
P^{Max} | Maximum performance | 100 |
γ | Per-capita pathogen growth rate | 0.5 |
β | Transmission rate | 0.1; 0.5 |
μ | Maximum per-capita pathogen killing rate | 0.25; 1.0 |
ω | Maximum performance loss per pathogen | 0.1; 0.2 |
c^{i}_{ρ} | Marginal costs of resistance | U [0, 0.1]; U [0.1, 0.2]; U [0.2, 0.5] |
c^{i}_{λ} | Marginal costs of tolerance | U [0, 0.1]; U [0.1, 0.2]; U [0.2, 0.5] |
ρ^{i} | Extra investment in resistance | U [0, 1]; U [0,0.5]; U [0.9, 1] |
λ^{i} | Extra investment in tolerance | U [0, 1]; U [0, 0.5]; U [0.9, 1] |
h^{i} | Rate of increase for the immune index | U [0, 0.001]; U [0, 0.01]; U [0, 0.1] |
Finally, effects of low, average and high levels of extra-investments in resistance and tolerance on ΔE_{I} and ΔE_{P} were quantified using fixed linear models (proc GLM on SAS^{®}[16]) that also contained the effects of β, μ, and ω for the characteristics of the pathogen, the averages at the population level of h^{i}, c^{i}_{ρ} and c^{i}_{λ} for the characteristics of the hosts, and all first-order interactions. The resulting least-squares estimates were used to identify epidemiological situations for which investments in tolerance, resistance or both were effective.
Results
Within-host pathogen dynamics
Between-host pathogen dynamics
When individuals invested more in resistance, only a fraction of the population got infected and AUC_{I} was low. For example, AUC_{I} decreased from 22,720 to 19,379 and 13,851 infected hosts in Figures 4b (ρ = 0.22), 4e (ρ = 0.46), and 4h (ρ = 0.94), respectively. When the average level of extra investments in tolerance was high (around 0.95), the impact of I_{t} on P_{t} was almost zero (Figures 4d,g and 4j). Otherwise, P_{t} decreased as I_{t} increased, especially for low levels of tolerance (Figures 4b,e and 4h). This should have translated in an increase in AUC_{P} but, in this particular population, costs associated with tolerance were high (around 0.15) and initial performance was low. For example, P_{0} averaged 79.8 in Figure 4g (λ = 0.95; AUC_{P} = 23,509) and 90.4 in Figure 4e (λ = 0.25; AUC_{P} = 24,779).
Effectiveness analyses
Incremental effectiveness of the performance of the population (ΔE_{P}) associated to different investments in individual tolerance (λ^{i}) and for selected values of c^{i}_{λ}, β and ω, as defined in Table 1
λ^{i}~ U[0,0.5] | λ^{i}~ U[0.9,1] | ||||||
---|---|---|---|---|---|---|---|
c ^{ i } _{ λ } | β | ω | ΔE _{ P } | c ^{ i } _{ λ } | β | ω | ΔE _{ P } |
0 | 0.1 | 0.1 | 1623 | 0 | 0.1 | 0.1 | 2974 |
U[0.2,0.5] | 0.1 | 0.1 | -628 | U[0.2,0.5] | 0.1 | 0.1 | 722 |
0 | 0.1 | 0.2 | 5635 | 0 | 0.1 | 0.2 | 6986 |
U[0.2,0.5] | 0.1 | 0.2 | 3383 | U[0.2,0.5] | 0.1 | 0.2 | 4734 |
0 | 0.5 | 0.1 | 4198 | 0 | 0.5 | 0.1 | 5549 |
U[0.2,0.5] | 0.5 | 0.1 | 1946 | U[0.2,0.5] | 0.5 | 0.1 | 3297 |
0 | 0.5 | 0.2 | 8210 | 0 | 0.5 | 0.2 | 9561 |
U[0.2,0.5] | 0.5 | 0.2 | 5958 | U[0.2,0.5] | 0.5 | 0.2 | 7309 |
Incremental effectiveness of the number of infected (ΔE_{I}) associated to different investments in individual resistance (ρ^{i}) and for selected values of c^{i}_{ρ}, μ and h^{i}, as defined in Table 1
ρ^{i}~ U[0,0.5] | ρ^{i}~ U[0.9,1] | ||||||
---|---|---|---|---|---|---|---|
c^{i}_{ρ} | h^{i} | μ | ΔE_{I} | c^{i}_{ρ} | h^{i} | μ | ΔE_{I} |
U[0.2, 0.5] | U [0, 0.001] | 0.25 | -808 | U[0.2, 0.5] | U [0, 0.001] | 0.25 | 1509 |
0 | U [0, 0.001] | 0.25 | -1215 | 0 | U [0, 0.001] | 0.25 | 1103 |
U[0.2, 0.5] | U [0, 0.001] | 1 | -127 | U[0.2, 0.5] | U [0, 0.001] | 1 | 2191 |
0 | U [0, 0.001] | 1 | -533 | 0 | U [0, 0.001] | 1 | 1784 |
U[0.2, 0.5] | U [0, 0.1] | 0.25 | 5607 | U[0.2, 0.5] | U [0, 0.1] | 0.25 | 7925 |
0 | U [0, 0.1] | 0.25 | 5201 | 0 | U [0, 0.1] | 0.25 | 7518 |
U[0.2, 0.5] | U [0, 0.1] | 1 | 6289 | U[0.2, 0.5] | U [0, 0.1] | 1 | 8607 |
0 | U [0, 0.1] | 1 | 5882 | 0 | U [0, 0.1] | 1 | 8200 |
Discussion
A general framework is proposed to provide insights into the effects of improved resistance and tolerance on the performance and size of an infected population. A clear distinction is made between effects of resistance on multiplication of the pathogen and effects of tolerance on damages induced by the pathogens. Hosts differ in the costs they incur to insure their particular levels of resistance and tolerance, and in the intensity of the response they mount against pathogens. Pathogens differ in their speed of spread between hosts, in virulence, and in the intensity of the response they elicit in the hosts. However, to be useful, the model must be validated and its limits and assumptions must be clarified, as will be discussed in the following, with examples mostly related to bovine mastitis.
Validation of the model
Model validation usually takes the form of a comparison between model outputs and real data but this was not possible here because reliable field data are scarce, difficult to measure or imprecisely defined [17, 18]. For example, estimates of costs associated with resistance and tolerance are limited in animals, in contrast to plants (see review by [19]). Tolerance has often been measured imprecisely as the overall ability to maintain fitness in the face of infection, irrespective of parasite burden. For example, cows infected with E. coli have been classified as moderate and severe responders according to milk production loss in the non-challenged quarters [20]. In this case, it is in reality a measure of the combined effects of resistance and tolerance [4]. It was also a deliberate choice to present a generic model because parameters values are different among disease and host populations, so model outputs for one specific disease may not apply to another disease. For example, transmission rates have been estimated at 0.20 to 1.50 per 1000 quarter-days at risk for S. uberis mastitis [21] but at 7 to 50 for S. aureus mastitis [22]. Similarly, killing rates have been estimated at 0.67 to 1.33 × 10^{-8} mL/cell per min in milk of cows [23, 24] and at 1.64 to 1.76 × 10^{-8} mL/neutrophil per min in dermis of rats inoculated with E. coli[13]. Model outputs will also depend on the virulence of the invading pathogens (ω), as exemplified by the different amount of milk loss at the first occurrence of clinical mastitis depending on bacteria species [25], and on the type of performance (e.g., yield, quality of products, or capacity for work) considered.
As an alternative form of validation, the dynamics of C^{i}_{t} and P^{i}_{t} at the individual, and of I_{t} and P_{t} at the herd levels were evaluated. For instance, as expected, C^{i}_{t} was lowest in resistant and P^{i}_{t} was highest in tolerant hosts (Figure 3), P_{t} remained stable across time when tolerance of the hosts was at its highest level, and I_{t} decreased faster when resistance of hosts was at its highest level (Figure 4). Results from the analysis of variance also validated the model. The null value of ΔE_{I} for h^{i}~U [0, 0.001] was sensible because, at this low rate, pathogens cannot be killed, regardless of how much was invested in resistance. The fact that β did not affect ΔE_{I} may also be explained by the same transmission in 'yes' and 'no' populations, so AUC_{I}^{yes} was close to AUC_{I}^{no} for any value of β. As a final example, ΔE_{P} was higher for β = 0.5 than for β = 0.1 because only few animals got infected with β = 0.1, so improving tolerance of these few hosts was not beneficial at the population level.
Limits and assumptions of the model
The strategy to build this model followed the current trend in epidemiology to begin with simple models and to add complexity only if the model fails to reproduce plausible epidemiological behaviors [26]. Several assumptions were made, some of which have been confirmed previously. One assumption was that available resources are partitioned between performance, resistance and tolerance. Indeed, experiences in poultry [27] and other species [28] have shown that individuals differ in their ability to allocate resources to their needs. This is also one of the factors evoked to explain the increased susceptibility of high yielding dairy cows to mastitis [29]. Lack of resources may lead to vicious cycles because hosts in poor condition are more susceptible to higher pathogen occurrence and infection intensity, which further weaken the condition of the host [30]. Another assumption is that investments in resistance and tolerance are linked through the constraint in equation 4 and this has been confirmed by [2], where a negative relationship was found between resistance and tolerance in rodent malaria.
Some assumptions of the model could also be relaxed with more complex equations that have been used in models examining the effects of mixed infection [21], infectious dose [31] and vaccination/treatment [32] on transmission dynamics. Resistance could vary as a function of exposure to disease [33]. Availability of external resource can vary across time, as in Doesch-Wilson et al. [34]. In the model used here, individual infectious contacts were assumed independent and at random but models with heterogeneous mixing [35] and that consider genetic susceptibility among relatives [36, 37] may be more appropriate. The course of infection within hosts can also be modelled more accurately, in line with the characteristics of the disease under study. For example, models with increasing complexity have been proposed to describe the fate of mastitis-causing E. coli in infected cows [23, 24]. Models for co-evolutionary mechanisms between host and pathogens should be considered [38] if the time scale is longer than the one used in this study.
Other assumptions may be difficult to verify. For such assumptions, a set of arbitrary standard values for the parameters and different forms for equations should be tested in so-called sensitivity analyses. For example, the amount of loss in performance was assumed directly associated with pathogen load, although the most dramatic changes may occur at low or subclinical levels of disease, with diminishing effects of each additional parasite [39].
Effectiveness analyses
Two results from the effectiveness analyses are noteworthy, although they must be further evaluated in empirical studies. One is that the range of possible values of ΔE_{P} and ΔE_{I} for the different input parameters (Figure 5) is wide. This emphasizes the need to accurately model the infection process and its impact on the population before deciding on the most effective strategy. For example, increasing host tolerance is theoretically less effective for improving performance of populations infected with pathogens that cause minor rather than major mastitis. Indeed, pathogens causing minor mastitis are less virulent (ω) and less transmissible (β) than those causing major mastitis [40], so modest advantages of high tolerance would be offset by the associated costs. Likewise, selecting for better resistance to mastitis would be effective to restrict the size of a population epidemic if animals are infected with bacterial strains that are likely to be killed by neutrophils [41], i.e. μ>0 in equation 3.
Another noteworthy observation is that least-squares means for ΔE_{P} were highest in highly tolerant populations, while ΔE_{I} did not change between different tolerance levels. This suggests that selection for increased tolerance would be effective under commercial constraints. This is different from models applied to natural populations that predict an increase in the overall incidence of infection as the frequency of tolerant hosts increases [38]. In natural populations, tolerant hosts survive longer than non-tolerant ones, thus keeping the disease longer in the population and increasing the risk of exposure to disease. Here, the model is for an endemic disease in a population under commercial contraints, in which non-tolerant animals are kept even if they are sick (no natural death, no culling). Consequently, the risk of exposure to disease does not change, even if the pathogen population size (C) increases.
In general, little is known about tolerance mechanisms in animals but their study should provide a good foundation for insuring health over the long term. Indeed, in the long term, advantages of being tolerant should be greater than those associated with resistance. For example, in non-evolving pathogen populations, advantages of being resistant decrease in parallel with the decline in disease frequency, while the advantages of being tolerant are maintained, or even increase if disease frequency rises [42]. In evolving pathogen populations, improved host resistance will pressure pathogens to evolve better mechanisms to evade host defense processes, potentially resulting in cyclical co-evolutionary dynamics. In contrast, tolerance does not interact directly with the pathogen and should not induce selection for counter-adaptations, although elevated levels of tolerance may allow pathogens to be more virulent [43].
where y_{ij}^{t} is the milk yield at time t of the i^{th} cow (yield corrected for fixed and non-genetic random effects estimated from the genetic evaluation model) infected with I_{ij}^{t}, i.e. the bacterial load for bacterial species j; b_{j} is the average tolerance against bacteria of strain j; B_{ij} describes individual random deviations from the average tolerance with B_{ij} ~ IID N(0, σ²_{b}); and e_{ij}^{t} are residuals with e_{ij}^{t} ~ N(0, V_{e}), where V_{e} accounts for the non-independence between repeated e_{ij}^{t}. Such information could be collected from quarters of experimentally infected cows, as was done in the study of [44].
Conclusions
In summary, this paper presents a novel epidemic model to explore the effects of tolerance and resistance on performance and disease spread in a population. Although more research is necessary to validate the model and more empirical studies are needed to obtain values for the input parameters, the analytic approach can be used to find optimal strategies of disease control in commercial populations.
Declarations
Acknowledgements
This study was supported by EADGENE (European Animal Disease Genomics Network of Excellence for Animal Health and Food Safety) and the University of Liege.
Authors’ Affiliations
References
- Jovanović S, Savić M, Živković D: Genetic variation in disease resistance among farm animals. Biotech Anim Husbandry. 2009, 25: 339-347.View Article
- Råberg L, Graham AL, Read AF: Decomposing health: tolerance and resistance to parasites in animal. Phil Trans R Soc B. 2009, 364: 37-49.PubMed CentralView ArticlePubMed
- Medzhitov R: Damage control in host-pathogen interactions. Proc Natl Acad Sci. 2009, 106: 5525-15526. 10.1073/pnas.0908451106.View Article
- Råberg L, Graham AL, Read AF: Decomposing health: tolerance and resistance to parasites in animal. Phil Trans R Soc B. 2009, 364: 37-49.PubMed CentralView ArticlePubMed
- Smith VH, Holt RD: Resource competition and within-host disease dynamics. Tree. 1996, 11: 386-389.PubMed
- Baer B, Armitage SA, Boomsma JJ: Sperm storage induces an immunity cost in ants. Nature. 2006, 441: 872-875. 10.1038/nature04698.View ArticlePubMed
- Kuukka-Anttila H, Peuhkuri N, Kolari I, Paananen T, Kause A: Quantitative genetic architecture of parasite-induced cataract in rainbow trout, Oncorhynchus mykiss. Heredity. 2009, 104: 20-27. 10.1038/hdy.2009.123.View Article
- Little TJ, Killick SC: Evidence for a cost of immunity when the crustacean Daphnia magna is exposed to the bacterial pathogen Pasteuria ramos. J Anim Ecol. 2007, 76: 1202-1207. 10.1111/j.1365-2656.2007.01290.x.View ArticlePubMed
- Read AF, Allen JE: Evolution and immunology: The economics of immunity. Science. 2000, 290: 1104-1105. 10.1126/science.290.5494.1104.View ArticlePubMed
- Coltherd JC, Bünger L, Kyriazakis I, Houdijk JGM: Genetic growth potential interacts with nutrition on the ability of mice to cope with Heligmosomoides bakeri infection. Parasitology. 2009, 136: 1043-1055. 10.1017/S0031182009006428.View ArticlePubMed
- Roberts MG, Heesterbeek JAP: Mathematical models in epidemiology. Mathematical models. Edited by: Filar JA, Krawczyk JB. 2004, UNESCO-EOLSS Oxford, [http://www.eolss.net]
- Detilleux JC: Neutrophils in the war against Staphylococcus aureus: Predator-prey models to the rescue. J Dairy Sci. 2004, 87: 3716-3724. 10.3168/jds.S0022-0302(04)73510-9.View ArticlePubMed
- Li Y, Karlin A, Loike JD, Silverstein SC: A critical concentration of neutrophils is required for effective bacterial killing in suspension. Proc Natl Acad Sci USA. 2002, 99: 8289-8294. 10.1073/pnas.122244799.PubMed CentralView ArticlePubMed
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81: 2340-2361. 10.1021/j100540a008.View Article
- Timms R, Colegrave N, Chan BHK, Read AF: The effect of parasite dose on disease severity in the rodent malaria Plasmodium Chabaudi. Parasitology. 2001, 123: 1-11. 10.1017/S0031182001008083.View ArticlePubMed
- SAS® Institute: SAS^{®} User's Guide. 1998, Cary: SAS® Institute, Version 8.12
- Lanzas C, Ayscue P, Ivanek R, Gröhn YT: Model or meal? Farm animal populations as models for infectious diseases of humans. Nat Rev Microbiol. 2010, 8: 139-148.PubMed
- Miller MR, White A, Boots M: The evolution of parasites in response to tolerance in their hosts: the good, the bad, and apparent commensalism. Evolution. 2006, 60: 945-956.View ArticlePubMed
- Brown JKM: Yield penalties of disease resistance in crops. Curr Opin Plant Biol. 2002, 5: 339-344. 10.1016/S1369-5266(02)00270-4.View ArticlePubMed
- Vandeputte-Van Messom G, Burvenich C, Roets E, Massart-Leën A, Heyneman R, Kremer WDJ, Brand A: Classification of newly calved cows into moderate and severe responders to experimentally induced Escherichia coli mastitis. J Dairy Res. 1993, 60: 19-29. 10.1017/S002202990002731X.View ArticlePubMed
- Zadoks RN, Allore HG, Barkema HW, Sampimom OC, Gröhn YT, Schukken YH: Analysis of an outbreak of Streptococcus uberis mastitis. J Dairy Sci. 2001, 84: 590-599. 10.3168/jds.S0022-0302(01)74512-2.View ArticlePubMed
- Zadoks RN, Allore HG, Hagenaars TJ, Barkema HW, Schukken YH: A mathematical model of Staphylococcus aureus control in dairy herds. Epidemiol Infect. 2002, 129: 397-416. 10.1017/S0950268802007483.PubMed CentralView ArticlePubMed
- Detilleux JC, Vangroenweghe F, Burvenich C: Mathematical model of the acute inflammatory response to Escherichia coli in intramammary challenge. J Dairy Sci. 2006, 89: 3455-3465. 10.3168/jds.S0022-0302(06)72383-9.View ArticlePubMed
- White LJ, Schukken YH, Dogan B, Green L, Döpfer D, Chappell MJ, Medley GF: Modelling the dynamics of intramammary E. coli infections in dairy cows: understanding mechanisms that distinguish transient from persistent infections. Vet Res. 2010, 41: 13-10.1051/vetres/2009061.PubMed CentralView ArticlePubMed
- Gröhn YT, Wilson DJ, González RN, Hertl JA, Schulte H, Bennett G, Schukken YH: Effect of pathogen-specific clinical mastitis on milk yield in dairy cows. J Dairy Sci. 2004, 87: 3358-3374.View ArticlePubMed
- Gilligan CA, van den Bosch F: Epidemiological models for invasion and persistence of pathogens. Ann Rev Phytopathol. 2008, 46: 385-418. 10.1146/annurev.phyto.45.062806.094357.View Article
- Gross WP, Siegel PB, Pierson FW: Effects of genetic selection for high and low antibody response on resistance to a variety of disease challenge and the relationship of resource allocation. Avian Dis. 2002, 46: 1007-1010. 10.1637/0005-2086(2002)046[1007:EOGSFH]2.0.CO;2.View ArticlePubMed
- Zuk M, Stoehr AM: Immune defense and host life history. Am Nat. 2002, 160: S09-S28. 10.1086/342131.View Article
- Ingvartsen KL, Dewhurst RJ, Friggens NC: On the relationship between lactational performance and health: is it yield or metabolic imbalance that cause production diseases in dairy cattle? A position paper. Livest Prod Sci. 2003, 83: 277-308. 10.1016/S0301-6226(03)00110-6.View Article
- Beldomenico PM, Begon M: Disease spread, susceptibility and infection intensity: vicious circles?. Trends Ecol Evol. 2010, 25: 21-27. 10.1016/j.tree.2009.06.015.View ArticlePubMed
- Pujol JM, Eisenberg JE, Haas CN, Koopman JS: The effect of ongoing exposure dynamics in dose response relationships. PLoS Comput Biol Epub. 2009, 5:
- van der Goot JA, Koch G, de Jong M: Quantification of the effect of vaccination on transmission of avian influenza (H7N7) in chickens. Proc Natl Acad Sci USA. 2005, 102: 18141-18146. 10.1073/pnas.0505098102.PubMed CentralView ArticlePubMed
- Dybiec B: SIR model of epidemic spread with accumulated exposure. Eur Phys J. 2009, 67: 377-383.View Article
- Doeschl-Wilson AB, Brindle W, Emmans GC, Kyriazakis I: Unraveling the relationship between animal growth and immune response during micro-parasitic infections. PLoS ONE. 2009, 4: e7508-10.1371/journal.pone.0007508.PubMed CentralView ArticlePubMed
- Brauer F: Epidemic models with heterogeneous mixing and treatment. Bull Math Biol. 2008, 70: 1869-1885. 10.1007/s11538-008-9326-1.View ArticlePubMed
- Detilleux JC: Genetic management of infectious diseases: a heterogeneous epidemio-genetic model illustrated with S. aureus mastitis. Genet Sel Evol. 2005, 37: 437-453. 10.1186/1297-9686-37-5-437.PubMed CentralView ArticlePubMed
- Nieuwhof GJ, Conington J, Bishop SC: A genetic epidemiological model to describe resistance to an endemic bacterial disease in livestock: application to footrot in sheep. Gen Sel Evol. 2009, 41: 19-10.1186/1297-9686-41-19.View Article
- Roy BA, Kirchner JW: Evolution dynamics of pathogen resistance and tolerance. Evolution. 2000, 54: 51-63.View ArticlePubMed
- Hawkins CD, Morris RS: Depression of productivity in sheep infected with Fasciola hepatica. Vet Parasitol. 1978, 4: 341-351. 10.1016/0304-4017(78)90020-1.View Article
- White LJ, Lam TJ, Schukken YH, Green L, Medley G, Chappell MJ: The transmission and control of mastitis in dairy cows: a theoretical approach. Prev Vet Med. 2006, 17: 67-83. 10.1016/j.prevetmed.2006.01.008.View Article
- Rambeaud M, Almeida RA, Pighetti GM, Oliver SP: Dynamics of leukocytes and cytokines during experimentally induced Streptococcus uberis mastitis. Vet Immunol Immunopathol. 2003, 96: 193-205. 10.1016/j.vetimm.2003.08.008.View ArticlePubMed
- Boots M: Fight or learn to live with the consequences. Trends Ecol Evol. 2008, 23: 248-250. 10.1016/j.tree.2008.01.006.View ArticlePubMed
- Svensson E, Råberg L: Resistance and tolerance in animal enemy-victim coevolution. Trends Ecol Evol. 2010, 25: 267-274. 10.1016/j.tree.2009.12.005.View ArticlePubMed
- Vangroenweghe F, Rainart P, Paape M, Duchateau L, Burvenich C: Increase of Escherichia coli inoculum doses induces faster innate immune response in primiparous cows. J Dairy Sci. 2004, 87: 4132-4144. 10.3168/jds.S0022-0302(04)73556-0.View ArticlePubMed
Copyright
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.