Antti Kause1* and Jørgen Ødegård2- 1Biotechnology and Food Research, Biometrical Genetics, MTT Agrifood Research Finland, Jokioinen, Finland
- 2Nofima, Ås, Norway
Tolerance to infections is defined as the ability of a host to limit the impact of a given pathogen burden on host performance. Uncoupling resistance and tolerance is a challenge, and there is a need to be able to separate them using specific trait recording or statistical methods. We present three statistical methods that can be used to investigate genetics of tolerance-related traits. Firstly, using random regressions, tolerance can be analyzed as a reaction norm slope in which host performance (y-axis) is regressed against an increasing pathogen burden (x-axis). Genetic variance in tolerance slopes is the genetic variance for tolerance. Variation in tolerance can induce genotype re-ranking and changes in genetic and phenotypic variation in host performance along the pathogen burden trajectory, contributing to environment-dependent genetic responses to selection. Such genotype-by-environment interactions can be quantified by combining random regressions and covariance functions. To apply random regressions, pathogen burden of individuals needs to be recorded. Secondly, when pathogen burden is not recorded, the cure model for time-until-death data allows separating two traits, susceptibility and endurance. Susceptibility is whether or not an individual was susceptible to an infection, whereas endurance denotes how long time it took until the infection killed a susceptible animal (influenced by tolerance). Thirdly, the normal mixture model can be used to classify continuously distributed host performance, such as growth rate, into different sub-classes (e.g., non-infected and infected), which allows estimation of host performance reduction specific to infected individuals. Moreover, genetics of host performance can be analyzed separately in healthy and affected animals, even in the absence of pathogen burden and survival data. These methods provide novel tools to increase our understanding on the impact of parasites, pathogens, and production diseases on host traits.
Introduction
Tolerance and resistance are two different defense mechanisms to defend against pathogens and parasites. Resistance is the ability of a host to prevent pathogen entry and to control pathogen life cycle in a way to reduce pathogen burden within a host individual. Tolerance to infections, in turn, is defined as the ability of the host to limit the impact of a given pathogen burden on host health, performance, and ultimately on host fitness (Clunies-Ross, 1932; Painter, 1958; Albers et al., 1987; Simms and Triplett, 1994; Simms, 2000) (Figure 1).
Figure 1. Tolerance to infections. Tolerance is the reaction norm slope of host performance regressed against individual's pathogen burden. The lines represent performance of three genotypes with a different degree of tolerance.
Being able to uncouple resistance and tolerance is essential for several reasons. Firstly, they have different impact on the arms-race co-evolution between the host and the pathogen (Mauricio et al., 1997; Rausher, 2001; Bishop and MacKenzie, 2003; Best et al., 2008). Moreover, both in animals and plants, tolerance and resistance are weakly genetically correlated, and thus they are genetically different traits (Leimu and Koricheva, 2006; Ødegård et al., 2011b; Kause et al., 2012). Finally, animal and plant breeders should exploit both increased resistance and tolerance to ensure global food security.
In addition to pathogens, tolerance can be assessed against abiotic factors such as temperature, heavy metals, or against production diseases causing damage to body tissues (Ravagnolo and Misztal, 2000a,b; Schat et al., 2002; Bloemhof et al., 2012; Kause et al., 2012). Naturally, production diseases, such as ascites, are not standard disease traits caused by a pathogen or parasite infection. Thus, there is no co-evolution between a host and a production disease, and the production disease does not evolve in response to the evolution of the host. Nevertheless, improved resistance and tolerance can be both used to reduce the harmful effects of production diseases on farmed animals, motivating their tolerance analysis (Kause et al., 2012). From hereon in this paper, pathogen burden is used as a general term to refer to a pathogen load of an individual, for instance, number or biomass of ecto- and endoparasites, number of pathogens in a blood sample, or severity of a production disease. In plants, pathogen burden may refer to the biomass or number of herbivores, or percentage of leaf area lost to herbivores.
The objective of this paper is to present recent statistical advances in the genetic analysis of tolerance-related traits. Firstly, random regression models have been applied to tolerance analysis. They allow a sophisticated genetic analysis of traits defined as functions as well as the quantification of genotype-by-environment interactions (G × E) induced by infections (Kause, 2011; Kause et al., 2012). Secondly, Ødegård et al. (2011b,c) introduced a cure model to separate “susceptibility” and “endurance” from challenge test data with time-until-death observations (without having any knowledge about infection status of the animals). The first trait is comparable to resistance, while endurance may be influenced by tolerance. “Susceptibility” can be defined as whether or not the animal is liable to die as a result of an infection (i.e., long-term survival, which is likely associated with resistance), while “endurance” is defined as how long time it takes before a potential infection kills the animal (which is likely associated with tolerance). Both endurance and susceptibility may show genetic variation, and may be viewed as different genetic factors affecting survival under an infection. Finally, normal mixture models can be extended to involve responses in host performance traits (e.g., growth curves) specific to healthy and affected individuals (Wang and Bodner, 2007; Madsen et al., 2008).
Random Regression Models
Tolerance is by definition the change in host performance as a function of pathogen burden (Simms, 2000), and hence, it is natural to apply random regression models to estimate genetic parameters and breeding values for tolerance (Kause, 2011). Using random regressions, tolerance can be analyzed as a reaction norm in which host performance (on y-axis) is regressed against pathogen burden of individuals (on x-axis) (Box 1). It is important to note that pathogen burden is measured separately from each individual, and it is not a general environmental characteristic. The slope of such a regression is consistent with the definition of tolerance (Figure 1), and hence genetic variance in regression slopes is the genetic variance for tolerance (Kause, 2011).
Box 1. A random regression model.
The intercept of the tolerance regression is interpreted as the host performance in a pathogen-free environment, and the genetic correlation between the slope and the intercept quantifies the degree to which host performance under no infection is genetically traded off with tolerance. Moreover, genetic correlations of the slope and intercept with third-party traits can be estimated by extending the random regression model to multitrait animal or sire model (Kause et al., 2012).
In animals, pathogen burden is typically a continuously distributed trait, especially when a population is under a natural pathogen infection (Stear et al., 1995; Kuukka-Anttila et al., 2010). Even in a challenge test in which all individuals are exposed to the same initial pathogen load, variation among individuals in resistance creates continuous variation in pathogen burden. Random regression models allow genetic analysis of tolerance along a continuous pathogen burden trajectory. In animal breeding, random regression models have been commonly applied to the reaction norm analysis of G × E (Henderson, 1982; Meyer and Hill, 1997; Calus et al., 2004; Schaeffer, 2004; Lillehammer et al., 2009).
Tolerance-Induced Variation in Host Performance
Genetic variation in tolerance may induce G × E in host performance, leading to changes in genetic variation of host performance along an increasing pathogen burden trajectory. For instance, in Figure 1, genetic variance in host performance is elevated along increased pathogen burden due to diverging tolerance reaction norms. In poultry, pigs, and aquaculture species, breeding nucleuses may be held infection-free due to biosecurity reasons, whereas commercial production and/or collection of sib and progeny information for breeding value estimation occurs at field farms with diverse diseases present. Such a design may induce G × E due to variation in the level of tolerance, which should be accounted for in breeding value evaluations.
In an infection-free environment, individual variation in host performance, e.g., in growth rate, is due to variation in genetic potential for growth and unexplained environmental variation. Under infection, in turn, individual variation in both resistance and tolerance induce additional variation into host performance. Some individuals are fully resistant or are not exposed to an infection, and thus their growth is not influenced by the infection. Some individuals are infected, and the degree to which their growth rate is reduced depends on their pathogen burden and the level of tolerance. Growth of fully tolerant individuals is not affected, whereas growth of very sensitive ones is greatly reduced.
Despite the large number of studies dealing with the changes induced by biotic (e.g., diet) and abiotic factors in general (Hoffmann and Merilä, 1999; Kause and Morin, 2001; Charmantier and Garant, 2005), there has been only a limited focus on infection-induced changes in genetic parameters and the consequent environment-specific genetic responses to selection (van der Waaij et al., 2000). Infections are indeed known to induce changes in heritability of host performance traits (Charmantier et al., 2004; Pakdel et al., 2005; Zerehdaran et al., 2006; Kause et al., 2007, 2012; Vehviläinen et al., 2008; Lewis et al., 2009). Yet, currently we do not know how much of the phenotypic variation in host performance is in fact created by infections and the associated tolerance. A study by Kause et al. (2012) showed that coefficient of phenotypic variation in broiler body weight was elevated from 11.5% when birds were healthy, to 19.1% when birds were severely affected by ascites. Similarly, coefficient of genetic variation was increased from 4.9% to 7.9%, implying the changes in variance can be extensive (Figure 2). It is hypothesized that in populations exposed to infections, a large proportion of phenotypic variance in host traits is induced by infections and the associated individual variation in resistance and tolerance.
Figure 2. Tolerance analysis using random regressions and covariance functions illustrated using data on 7-week body weight and heart ratio of broilers [reproduced from Kause et al. (2012); http://creativecommons.org/licenses/by/3.0/]. Heart ratio, the ratio of right ventricular weight to total heart weight, is an indicator of ascites resistance, the birds with higher than 27–30% heart ratio typically being ascitic (Wideman et al., 1998). Tolerance is the change in body weight as a function of increasing ascites severity, measured as the hearth ratio. (A) Population average tolerance curve with body weight on y-axis and ascites severity on x-axis. (B) Frequency distribution of estimated breeding values for tolerance slopes of individuals, showing sensitive (steep negative slope) and more tolerant (weak negative slope) genotypes. (C) Increased coefficients of phenotypic and genetic variation in body weight as a function of ascites severity, showing ascites molds trait variation. (D) Genetic correlation between healthy birds and birds with different degree of ascites severity, showing ascites creates genotype re-ranking (Kause et al., 2012).
Random regression models combined with covariance functions (Kirkpatrick et al., 1990; Meyer and Hill, 1997) provide means to quantify the changes in phenotypic and genetic variances in host traits along a continuous pathogen burden trajectory (Kause, 2011; Kause et al., 2012). Given the genetic (co)variance estimates of tolerance slope and intercept estimated using random regressions, the changes in genetic variance in host performance can be calculated using formulas (Box 1; Figure 2). The same logic can be applied to the maternal and environmental components of (co)variance.
Tolerance-Induced Genotype Re-Ranking in Host Performance
Crossing tolerance reaction norms create genotype re-ranking in host performance traits across pathogen burden trajectory. This is similar to any genotype re-ranking across environmental gradients (Via and Lande, 1985), with the difference that now the environment is pathogen burden of individuals (Kause et al., 2012). The two forms of G × E, scaling effect and genotype re-ranking, facilitate environment-dependent genetic responses, yet the re-ranking is more severe issue for selective breeding because genotypes in one environment are not necessarily the best ones in the other environments. Re-ranking across environments can be quantified by a genetic correlation between measurements in two environments for a given trait (Falconer, 1952).
The degree of re-ranking between any two pathogen burden levels can be calculated by combining random regression results with covariance functions (Box 1). For instance, ascites induced moderate genotype re-ranking in broiler body weight, the genetic correlation of healthy birds with weakly affected birds being unity but with severely affected birds 0.45 (Kause et al., 2012; Figure 2). In field data sets with multiple environments, infection pressure is typically not the only environmental factor varying across environments, yet the effect of pathogen burden on G × E could be revealed using a combined reaction norm and multi-trait model (Windig et al., 2011), in which pathogen burden is modeled as a continuous reaction norm and the discrete environments capturing other environmental factors are modeled as separate discrete traits. Performing extensive infection-challenge tests is impractical in many farm animal species, but the combined reaction norm and multi-trait model may be an effective additional method for revealing the degree of G × E induced by infections.
Infections do not induce only genotype re-ranking and a change in variance but also changes in the correlation structure of resistance, growth, and reproduction traits (de Greef et al., 2001; Kause et al., 2005, 2012; Zerehdaran et al., 2006; Kuukka-Anttila et al., 2010). The modification of genetic architecture of host traits by pathogens, parasites, and production diseases, mediated by tolerance genetics, may play a more fundamental role in animal breeding and microevolution than has been previously thought.
Data Requirements for Random Regression
Obtaining a solid x-axis is a major challenge for the tolerance analysis in animals because the x-axis should consists of individual-level quantitative data on pathogen burden (e.g., number of parasites, pathogen biomass). Qualitative data on burden (infected vs. non-infected individuals) creates biased estimates of genetic variance for tolerance (Kause, 2011). Moreover, if the x-axis consists of the average burdens of each environment, rather than individual-level burden measurements, then high host performance of a genotype at a given pathogen burden can be a result of high resistance and/or high tolerance, impeding a proper tolerance analysis. The analyzed host performance trait, in turn, can be feed intake, growth, reproduction, survival or a physiological trait, which together can be used to reveal mechanisms contributing to variation among genotypes in tolerance.
A split-family design with both an infection-free control and an experimental challenge test is the most effective design for tolerance analysis. In this way, infected animals are a random sample of their family and thus there will be a real causal relationship between host performance and pathogen burden (Tiffin and Inouye, 2000; Kause, 2011). This requires, however, that all the challenged individuals get the same pathogen burden level. This rarely is the case because individuals have innate individual variation in resistance, creating variation in pathogen burden even in a challenge test. Variation in resistance can be potentially related to the host performance traits used on the y-axis in the tolerance regression, biasing the estimate of genetic variation for tolerance (Tiffin and Inouye, 2000; Kause, 2011). As an alternative to the control-and-challenge test design, all individuals can be first recorded under infection-free conditions (e.g., for mature body weight), and then re-recorded after experimental exposition to equal pathogen burden level. However, such an analysis is unjustified in cases in which host performance shows natural temporal variation (e.g., variation in growth curves), which is thus confounded with tolerance (Albers et al., 1987; Bisset and Morris, 1996; Woolaston and Windon, 2001). Trypanotolerance of African cattle has been analyzed as a change in body weight in response to an experimental infection by Trypanosoma congolense, but although the number of parasites in the blood of individuals was recorded, it was not used to standardize the host performance changes of individuals (Hanotte et al., 2003; van der Waaij et al., 2003).
Under naturally occurring infection, it is possible that either high (or low) performing individuals are infected, leading to biased estimates of genetic variation for tolerance (Tiffin and Inouye, 2000; Kause, 2011). This is a major weakness of field data sets, because it is well established that individuals with initially different growth or life-history trait levels may be differently exposed to infections, parasites and production diseases (Arendt, 1997; Rauw et al., 1998), confounding the cause-and-effect relation between pathogen burden and reduction in host performance.
Random regression models require large sample sizes, e.g., within sire families. Decrease in family size leads to upward-biased genetic variance estimates for tolerance slope (Kause, 2011). This can be illustrated in a sire model set up. When a small number of individuals are sampled for each sire family, the sample is no longer representative of the true distribution and single observations have strong impact on the slope estimate. For some families the slope is underestimated, for others overestimated, and thus genetic variance estimate for slope is artificially increased. With heritability of 0.3 for tolerance slope, more than 50 sibs per family are required to obtain unbiased estimates of slope variance using a sire model analysis (Kause, 2011). Moreover, genetic correlation between tolerance slope and intercept is easily biased downward when family size is low. An upward (downward) bias in the slope of a family pushes the intercept downward (upward), creating an artificial negative genetic trade-off when it does not exist in reality. This can be avoided by using high family sizes and high number of non-infected individuals that force the intercept of a genotype to be placed close to the real value (Mauricio et al., 1997; Kause, 2011).
When each host individual has only a single performance record, it is possible to estimate genetic variance and breeding values for tolerance slope, but not its residual variance. Heritabilities for tolerance slope can be estimated when each individual has several performance observations, e.g., the initial performance under conditions of no infection and thereafter the performance after an infection. By using regression slopes of individuals as raw observations in the genetic analysis, both environmental and genetic components of slope variance and heritability can be estimated (Schaeffer, 2004).
Random regression can be applied to non-linear reaction norms (Kirkpatrick et al., 1990; Meyer and Hill, 1997; Schaeffer, 2004) and plateau-linear regression models (Ravagnolo and Misztal, 2000a,b; Kause et al., 2012), and thus the impact of pathogens on host performance does not need to be analyzed as a linear relationship.
A Cure Model for Time-Until-Death Data
The random regression approach requires individual-level data on pathogen burden which may be challenging to record. The cure model for time-until-death data provides a possibility to analyze genetics of resistance (or susceptibility) and endurance without a need for pathogen burden recording.
Many studies, especially on aquaculture species, have analyzed survival or time-until-death in a challenge test in which individuals are experimentally exposed to a specific pathogen (Ødegård et al., 2011a). Moreover, survival analysis has been applied to time-until-death data when mortality factors remain unknown (e.g., Ducrocq and Casella, 1996; Serenius and Stalder, 2004; Vehviläinen et al., 2010). A typical assumption in such analyses is that individuals with high probability of survival are resistant. However, an individual can survive if it has either high resistance, or low resistance but high tolerance (Figure 3), or was never exposed to a pathogen. The cure survival models are used for modeling of time-until-death data which include a fraction of non-susceptible animals, i.e., animals that are not liable to die as a result of the infection (Farewell, 1982). Ødegård et al. (2011c) developed a cure model aiming to distinguish two traits, “susceptibility” and “endurance,” from time-until-death data. These two concepts may be comparable with resistance and tolerance.
Figure 3. Contribution of resistance and tolerance to mortality due to a specific pathogen. Only individuals without resistance and tolerance will eventually die given a sufficiently long follow-up period. When having a limited follow-up period, individuals with high tolerance may still be alive at the end of an experiment.
In a survival analysis the infection status of each animal is typically unknown. Under pathogen attack, some animals may be fully capable of avoiding death (non-susceptible), either by resisting the infection, or by a successful recovery after the initial infection due to high tolerance (Figure 3). Furthermore, the degree of tolerance may also vary among the susceptible individuals, potentially causing variation in their expected time-until-death. As mortality is usually recorded over a limited follow-up period, a fraction of susceptible animals are also likely to be alive at the time of recording. For susceptible animals, the ability to survive depends on the expected time-until-death of the animal, which may show genetic variation. Hence, analogy of the terms “endurance” and “susceptibility” with tolerance and resistance are not necessarily clear-cut in a survival analysis, due to the fact that one only observes the extreme outcomes of an infection (whether or not an animal dies). Although “endurance” and “susceptibility” are impossible to separate on individual survivors, these two factors may still be distinguished on a family level using longitudinal survival analysis (i.e., short-term mortality rates vs. long-term survival).
A classical survival analysis of time-until-death assumes that all individuals are at risk and that all will eventually die given a sufficiently long follow-up period. When studying lifespan in general this is necessarily true, but may not hold when testing for mortality due to a specific pathogen. For non-susceptible animals time-until-death will necessarily be censored, irrespective of the follow-up time, and survival time may thus be a poor indicator for specific pathogen resistance. The endurance reflects the expected mortality per time-unit among susceptible individuals, but will have no effect on survival of the non-susceptible individuals (Farewell, 1982).
The survivors are likely a mixture of non-susceptible long-term survivors and a fraction of susceptible (but highly endure) animals being still alive, and the true condition of each animal is unknown (unless the animal dies). In the cure model, probabilities of the alternative settings (non-susceptible or susceptible but still alive) can be estimated while simultaneously taking into account variation in endurance among the surviving animals (Ødegård et al., 2011c; Box 2).
Box 2. A cure model.
The cure model has been applied to time-until-death data in farmed shrimp challenge-tested with the Taura syndrome virus (Ødegård et al., 2011b). It was estimated that although 72% of the shrimp survived, only 62% could be considered non-susceptible. The underlying heritability (±SE) for susceptibility was high (0.41 ± 0.07), while the heritability of endurance was low, albeit significant (0.07 ± 0.03). The most striking result was that endurance and susceptibility were seemingly distinct genetic traits (rG = 0.22 ± 0.25). The low genetic variation for endurance and the genetic independency of endurance and susceptibility are in line with the results on other animal species (Kause et al., 2012). These results have substantial impact on how disease challenge-testing should be performed. If the aim is to improve long-term survival under an infection pressure, selective breeding should focus on susceptibility. This implies that the follow-up period should continue until the vast majority of susceptible animals have died, ensuring that the observed end-survival largely resembles the fraction of non-susceptible animals in the population.
Normal Mixture Models
Normal mixture models can be used to analyze genetics of host performance, e.g., growth rate, within a population consisting of individuals affected and unaffected by a pathogen, even in the absence of pathogen burden and time-until-death data.
Finite normal mixture models have earlier been proposed for analysis of infection-affected, continuously distributed phenotypes, assuming that the true infection statuses of individuals are unknown (Detilleux and Leroy, 2000; Ødegård et al., 2003, 2005; Gianola et al., 2004). The mixture model attempts to identify hidden categories (e.g., non-infected and infected) among the observations, assuming that the continuous scale observations originate from two normal distributions differing in mean and (potentially) variance (Box 3; Figure 4). For instance, the broiler ascites example given in Figure 2 can be analyzed using a mixture model analysis assuming that the heart ratio has two underlying distributions, one for non-infected and one for ascitic birds (Zerehdaran et al., 2006). Another example of a mixture trait is somatic cell scores in milk of dairy cattle (Madsen et al., 2008). Somatic cell score is at low level in non-infected cows, but increase to high levels in cases of (unobserved) subclinical mastitis. Hence, the observed somatic cell scores may be viewed as a mixture of two normal distributions (non-infected and mastitic). In Atlantic salmon Salmo salar L., diseases such as infectious pancreas necrosis and pancreas disease can kill a fraction of the animals, but may also reduce subsequent growth of affected survivors. Hence, after an outbreak, observed growth of survivors may be viewed as a mixture trait depending on the individuals' previous health status.
Box 3. A normal mixture model.
Figure 4. An example of a two-component mixture distribution. The dotted lines are the unobserved distributions of non-infected “healthy” individuals (70% of the observation) with ~N(−1.0, 1.0) and infected “diseased” individuals (30%) with ~N(1.0, 1.0). The solid line represents the resulting distribution of the observed phenotypes. Trait values are given on x-axis and the frequencies of observations on y-axis.
Classical selection aims at changing a trait in the desired direction. However, for mixture traits the variation is partly explained by mixing of the two (or more) sub-distributions with different means, and partly by variation within each sub-distribution (Figure 4). Hence, if the aim is to reduce the incidence of the infection rather than altering the observed continuous host trait itself, simple directional selection for the latter (e.g., for somatic cell score) may not be optimal. The mixture model opens new possibilities for selection, and can be used to directly select for reduced infection risk. Additionally, the trait recorded on infected and non-infected animals may be viewed as two distinct sub-traits whose genetic variances and their genetic correlation can be estimated. This resembles the G × E analysis performed with random regression models (Figure 2) with the difference that the mixture model does not take into account that infected individuals may have different pathogen burdens.
Normal mixture models typically assume that an individual is either infected or not, and that infection has a certain effect on the phenotype (Figure 4). However, variation in environmental pathogen load and in individual tolerance for the infection imply that the effect of an infection may vary substantially among individuals and environments. The proposed mixture models may be extended to allow for individual responses to infection (Madsen et al., 2008). Alternatively, the model may be extended to a growth mixture model (Wang and Bodner, 2007). The growth mixture models assume that the observations come from different latent trajectories, i.e., health status does not only affect the expectation of individual observations, but also the slope of a phenotypic trajectory (growth curves). Infected and non-infected animals could show different trajectories, with the non-infected ones being unaffected by the pathogen, while the infected individuals being variably affected by the pathogen burden. Such models may be useful to analyze resistance and infection-affected traits observed on animals with unknown infection status and in environments with variable pathogen loads.
Application of the Methods in Breeding Programmes
Random regression models are routinely applied in farm animal breeding programs, e.g., for milk test-day models in dairy cows and for growth curves (Schaeffer, 2004). Similarly, random regression models can be implemented to select for tolerance, given suitable data are available. The cure model approach for the analysis of time-until-death data (Veerkamp et al., 2001; Ødegård et al., 2011a) have been implemented in the DMU software, allowing the estimation of genetic parameters and breeding values for practical breeding (Madsen and Jensen, 2010). To our knowledge, the cure model has not been implemented in routine genetic evaluations in any breeding program. Árnason (1999) and Urioste et al. (2007) have proposed a bivariate linear-threshold model which can be used to analyze whether an animal survived (a threshold trait) and how long it took until death (a linear trait). Such a model resembles the cure model and is straightforward to apply in multi-trait breeding value evaluations. Also the normal mixture model has been implemented in the DMU software (Madsen and Jensen, 2010), and is therefore available for multi-trait genetic evaluations, but to our knowledge, has not yet been implemented in routine genetic evaluations.
The cure model has been applied to survival data in aquaculture species, leading to altered recommendations for routine disease-challenge testing (Ødegård et al., 2011b). Historically, challenge tests in aquaculture species have been terminated at intermediate cumulative mortalities to ensure maximum variation in binary survival data. However, this approach is only proper given that endurance and susceptibility are equivalent traits, which is not necessarily the case. The current advice is to continue testing until mortality naturally ceases, even at levels above 50% mortality (Ødegård et al., 2011a).
So far, only a limited number of breeding programs have considered selecting for tolerance. Some African cattle breeding programs are specifically selecting for trypanotolerance-related traits, the tolerance being a major breeding objective trait (Hanotte et al., 2003; van der Waaij et al., 2003). In contrast, regardless of the extensive studies conducted in Australia and New Zealand on nematode tolerance in sheep, a decision has been made not to record and select for tolerance because of the need to let animals to suffer and production to be reduced for tolerance to be expressed (Albers et al., 1987; Bisset and Morris, 1996; Woolaston and Windon, 2001). The novel statistical methods and the increasing awareness of the detailed physiological mechanisms of tolerance (Medzhitov et al., 2012) may provide more opportunities for tolerance selection in farm animals.
Conclusions
The recent statistical developments provide tools to increase our understanding of genetics of alternative strategies to defend against parasites, pathogens, and production diseases. Most of the statistical methods can be applied in breeding value evaluations to breed for tolerance. Moreover, the methods presented here provide tools to quantify genotype-by-pathogen burden interactions that may explain a significant proportion of phenotypic variation in traits within populations that are exposed to various infections and production diseases. The traits whose variation is affected are typically production traits that are selected for in breeding programs. To be able to unambiguously select for the genetic potential of a production trait, the effects of resistance and tolerance should be separated from it. The methods presented in this paper provide potential to construct more effective breeding programs to increase both productivity and animal health.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Albers, G. A., Gray, G. D., Piper, L. R., Barker, J. S., Le Jambre, L. F., and Barger, I. A. (1987). The genetics resistance and resilience to Haemonchus contortus infection in young merino sheep. Int. J. Parasitol. 17, 1355–1363.
Arendt, J. D. (1997). Adaptive intrinsic growth rates: an integration across taxa. Q. Rev. Biol. 72, 149–177.
Árnason, T. (1999). Genetic evaluation of Swedish standard-bred trotters for racing performance traits and racing status. J. Anim. Breed. Genet. 116, 387–398.
Best, A., White, A., and Boots, M. (2008). Maintenance of host variation in tolerance to pathogens and parasites. Proc. Natl. Acad. Sci. U.S.A. 105, 20786–20791.
Bishop, S. C., and MacKenzie, K. M. (2003). Genetic management strategies for controlling infectious diseases in livestock populations. Genet. Sel. Evol. 35, S3–S17.
Bisset, S. A., and Morris, C. A. (1996). Feasibility and implications of breeding sheep for resilience to nematode challenge. Int. J. Parasitol. 26, 857–868.
Bloemhof, S., Kause, A., Knol, E. F., van Arendonk, J. A. M., and Misztal, I. (2012). Heat stress effects on farrowing rate in sows: genetic parameter estimation using within-line and crossbred models. J. Anim. Sci. 90, 2109–2119.
Calus, M. P. L., Bijma, P., and Veercamp, R. F. (2004). Effects of data structure on the estimation of covariance functions to describe genotype by environment interactions in a reaction norm model. Genet. Sel. Evol. 36, 489–507.
Charmantier, A., and Garant, D. (2005). Environmental quality and evolutionary potential: lessons from wild populations. Proc. R. Soc. Ser. B 272, 1415–1425.
Charmantier, A., Kruuk, L. E., and Lambrechts, M. M. (2004). Parasitism reduces the potential for evolution in a wild bird population. Evolution 58, 203–206.
Clunies-Ross, I. (1932). Observations on the resistance of sheep to infestations by the stomach worm Haemonchus contortus. J. Coun. Sci. Ind. Res. 5, 73–80.
de Greef, K. H., Janss, L. L., Vereijken, A. L., Pit, R., and Gerritsen, C. L. (2001). Disease-induced variability of genetic correlations: ascites in broilers as a case study. J. Anim. Sci. 79, 1723–1733.
Detilleux, J., and Leroy, P. L. (2000). Application of a mixed normal mixture model for the estimation of mastitis-related parameters. J. Dairy Sci. 83, 2341–2349.
Ducrocq, V., and Casella, G. (1996). A Bayesian analysis of mixed survival models. Genet. Sel. Evol. 28, 505–529.
Farewell, V. T. (1982). The use of mixture-models for the analysis of survival-data with long-term survivors. Biometrics 38, 1041–1046.
Gianola, D., Ødegård, J., Heringstad, B., Klemetsdal, G., Sorensen, D., Madsen, P., et al. (2004). Mixture model for inferring susceptibility to mastitis in dairy cattle: a procedure for likelihood-based inference. Genet. Sel. Evol. 36, 3–27.
Hanotte, O., Ronin, Y., Agaba, M., Nilsson, P., Gelhaus, A., Horstmann, R., et al. (2003). Mapping of quantitative trait loci controlling trypanotolerance in a cross of tolerant West African N'Dama and susceptible East African Boran cattle. Proc. Natl. Acad. Sci. U.S.A. 100, 7443–7448.
Henderson, C. R. (1982). Analysis of covariance in the mixed model: higher-level, non homogeneous, and random regressions. Biometrics 38, 623–640.
Hoffmann, A. A., and Merilä, J. (1999). Heritable variation and evolution under favourable and unfavourable conditions. Trends Ecol. Evol. 14, 96–101.
Kause, A. (2011). Genetic analysis of tolerance to infections using random regressions: a simulation study. Genet. Res. 93, 291–302.
Kause, A., and Morin, J. P. (2001). Seasonality and genetic architecture of development time and body size in the birch feeding sawfly Priophorus pallipes. Genet. Res. 78, 31–40.
Kause, A., Ritola, O., and Paananen, T. (2007). Changes in the expression of genetic characteristics across cohorts in skeletal deformations of farmed salmonids. Genet. Sel. Evol. 39, 529–543.
Kause, A., Ritola, O., Paananen, T., Wahlroos, H., and Mäntysaari, E. A. (2005). Genetic trends in growth, sexual maturity and skeletal deformations, and rate of inbreeding in a breeding programme for rainbow trout. Aquaculture 247, 177–187.
Kause, A., van Dalen, S., and Bovenhuis, H. (2012). Genetics of ascites resistance and tolerance in chicken: a random regression approach. G3 2, 527–535.
Kirkpatrick, M., Lofsvold, D., and Bulmer, M. (1990). Analysis of the inheritance, selection and evolution of growth trajectories. Genetics 124, 979–993.
Kolmodin, R., and Bijma, P. (2004). Response to mass selection when the genotype by environment interaction is modelled as a linear reaction norm. Genet. Sel. Evol. 36, 435–454.
Kuukka-Anttila, H., Peuhkuri, N., Kolari, I., Paananen, T., and Kause, A. (2010). Quantitative genetic architecture of parasite-induced cataract in rainbow trout, Oncorhynchus mykiss. Heredity 104, 20–27.
Leimu, R., and Koricheva, J. (2006). A meta-analysis of trade-offs between plant tolerance and resistance to herbivores: combining the evidence from ecological and agricultural studies. Oikos 112, 1–9.
Lewis, C. R. G., Torremorell, M., Galina-Pantoja, L., and Bishop, S. C. (2009). Genetic parameters for performance traits in commercial sows estimated before and after an outbreak of porcine reproductive and respiratory syndrome. J. Anim. Sci. 87, 876–884.
Lillehammer, M., Ødegård, J., and Meuwissen, T. H. E. (2009). Reducing the bias of estimates of genotype-by-environment interactions in random regression sire models. Genet. Sel. Evol. 41, 30.
Madsen, P., and Jensen, J. (2010). DMU: A user's guide. A package for analysing multivariate mixed models. Version 6, release 5.0. Available Online at: http://dmu.agrsci.dk/dmuv6_guide.5.0.pdf
Madsen, P., Shariati, M. M., and Ødegård, J. (2008). Genetic analysis of somatic cell score in Danish Holsteins using a liability-normal mixture model. J. Dairy Sci. 91, 4355–4364.
Mauricio, R., Rausher, M. D., and Burdick, D. S. (1997). Variation in the defense strategies of plants: are resistance and tolerance mutually exclusive? Ecology 78, 1301–1311.
Medzhitov, R., Schneider, D. S., and Soares, M. P. (2012). Disease tolerance as a defense strategy. Science 6071, 936–941.
Meyer, K., and Hill, W. G. (1997). Estimation of genetic and phenotypic covariance functions for longitudinal or ‘repeated’ records by restricted maximum likelihood. Livest. Prod. Sci. 47, 185–200.
Ødegård, J., Baranski, M., Gjerde, B., and Gjedrem, T. (2011a). Methodology for genetic evaluation of disease resistance in aquaculture species: challenges and future prospects. Aquacult. Res. 42, 103–114.
Ødegård, J., Gitterle, T., Madsen, P., Meuwissen, T. H. E., Yazdi, H. Y., Gjerde, B., et al. (2011b). Quantitative genetics of taura syndrome resistance in pacific white shrimp (Penaeus vannamei): a cure model approach. Genet. Sel. Evol. 43, 14.
Ødegård, J., Madsen, P., Labouriau, R., Gjerde, B., and Meuwissen, T. H. E. (2011c). A sequential threshold cure model for genetic analysis of time-to-event data. J. Anim. Sci. 89, 943–950.
Ødegård, J., Jensen, J., Madsen, P., Gianola, D., Klemetsdal, G., and Herinstad, B. (2003). Detection of mastitis in dairy cattle by use of mixture models for repeated somatic cell scores: a Bayesian approach via Gibbs sampling. J. Dairy Sci. 86, 3694–3703.
Ødegård, J., Madsen, P., Gianola, D., Klemetsdal, G., Jensen, J., Heringstad, B., et al. (2005). A Bayesian threshold-normal mixture model for analysis of a continuous mastitis-related trait. J. Dairy Sci. 88, 2652–2659.
Pakdel, A., van Arendonk, J. A. M., Vereijken, A. L., and Bovenhuis, H. (2005). Genetic parameters of ascites-related traits in broilers: effect of cold and normal temperature conditions. Br. Poult. Sci. 46, 35–42.
Rauw, W. M., Kanis, E., Noordhuizen-Stassen, E. N., and Grommers, F. J. (1998). Undesirable side effects of selection for high production efficiency in farm animals: a review. Livest. Prod. Sci. 56, 15–33.
Ravagnolo, O., and Misztal, I. (2000a). Genetic component of heat stress in dairy cattle: parameter estimation. J. Dairy Sci. 83, 2126–2130.
Ravagnolo, O., and Misztal, I. (2000b). Effect of heat stress on nonreturn rate in Holstein cows: genetic analyses. J. Dairy Sci. 85, 3092–3100.
Schaeffer, L. R. (2004). Application of random regression models in animal breeding. Livest. Prod. Sci. 86, 35–45.
Schat, H., Llugany, M., Vooijs, R., Hartley-Whitaker, J., and Bleeker, P. M. (2002). The role of phytochelatins in constitutive and adaptive heavy metal tolerances in hyperaccumulator and non-hyperaccumulator metallophytes. J. Exp. Bot. 53, 2381–2392.
Serenius, T., and Stalder, K. J. (2004). Genetics of length of productive life and lifetime prolificacy in the Finnish Landrace and Large White pig populations. J. Anim. Sci. 82, 3111–3117.
Simms, E. L., and Triplett, J. (1994). Costs and benefits of plant responses to disease: resistance and tolerance. Evolution 48, 1973–1985.
Stear, M. J., Bairden, K., Duncan, J. L., Gettinby, G., McKellar, Q. A., Murray, M., et al. (1995). The distribution of faecal nematode egg counts in Scottish Blackface lambs following natural, predominantly Ostertagia circumcincta infection. Parasitology 110, 573–581.
Tiffin, P., and Inouye, B. D. (2000). Measuring tolerance to herbivory: accuracy and precision of estimates made using natural versus imposed damage. Evolution 54, 1024–1029.
Urioste, J. I., Misztal, I., and Bertrand, J. K. (2007). Fertility traits in spring-calving Aberdeen Angus cattle. 1. Model development and genetic parameters. J. Anim. Sci. 85, 2854–2860.
van der Waaij, E., Hanotte, O., van Arendonk, J. A. M., Kemp, S. J., Kennedy, D., Gibson, A., et al. (2003). Population parameters for traits defining trypanotolerance in an F2 cross of N'Dama and Boran cattle. Livest. Prod. Sci. 84, 219–230.
van der Waaij, E. H., Bijma, P., Bishop, S. C., and van Arendonk, J. A. M. (2000). Modeling selection for production traits under constant infection pressure. J. Anim. Sci. 78, 2809–2820.
Veerkamp, R. F., Brotherstone, S., Engel, B., and Meuwissen, T. H. E. (2001). Analysis of censored survival data using random regression models. Anim. Sci. 72, 1–10.
Vehviläinen, H., Kause, A., Quinton, C., Koskinen, H., and Paananen, T. (2008). Survival of the currently fittest – Genetics of rainbow trout survival across time and space. Genetics 180, 507–516.
Vehviläinen, H., Kause, A., Quinton, C., Koskinen, H., and Paananen, T. (2010). Genetic architecture of rainbow trout survival from egg to adult. Genet. Res. 92, 1–11.
Via, S., and Lande, R. (1985). Genotype-environment interaction and the evolution of phenotypic plasticity. Evolution 39, 505–522.
Wideman, R. F. Jr., Wing, T., Kirby, Y. K., Forman, M. F., Marson, N., Tackett, C. D., et al. (1998). Evaluation of minimally invasive indices for predicting ascites susceptibility in three successive hatches of broilers exposed to cool temperatures. Poult. Sci. 77, 1565–1573.
Windig, J. J., Mulder, H. A., Bohthe-Wilhelmus, D. I., and Veerkamp, R. F. (2011). Simultaneous estimation of genotype-by-environment interaction accounting for discrete and continuous environmental descriptors in Irish dairy cattle. J. Dairy Sci. 94, 3137–3147.
Woolaston, R. R., and Windon, R. G. (2001). Selection of sheep for response to Trichostrongylus colubriformis larvae: genetic parameters. Anim. Sci. 73, 41–48.
Keywords: cure model, genotype-by-environment interaction, mixture model, quantitative genetics, random regression, resistance, statistical methods, tolerance
Citation: Kause A and Ødegård J (2012) The genetic analysis of tolerance to infections: a review. Front. Gene. 3:262. doi: 10.3389/fgene.2012.00262
Received: 29 March 2012; Accepted: 05 November 2012;
Published online: 14 December 2012.
Edited by:
Andrea B. Doeschl-Wilson, The University of Edinburgh, UKReviewed by:
Mairead L. Bermingham, The University of Edinburgh, UKSunday Peters, Cornell University, USA
Copyright © 2012 Kause and Ødegård. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Antti Kause, Biotechnology and Food Research, Biometrical Genetics, MTT Agrifood Research Finland, FI-31600, Jokioinen, Finland. e-mail: antti.kause@mtt.fi