Genetic evaluation of mastitis liability and recovery through longitudinal analysis of transition probabilities
- Jessica Franzén^{1, 2}Email author,
- Daniel Thorburn^{2},
- Jorge I Urioste^{1, 3} and
- Erling Strandberg^{1}
DOI: 10.1186/1297-9686-44-10
© Franzén et al; licensee BioMed Central Ltd. 2012
Received: 14 June 2011
Accepted: 4 April 2012
Published: 4 April 2012
Abstract
Background
Many methods for the genetic analysis of mastitis use a cross-sectional approach, which omits information on, e.g., repeated mastitis cases during lactation, somatic cell count fluctuations, and recovery process. Acknowledging the dynamic behavior of mastitis during lactation and taking into account that there is more than one binary response variable to consider, can enhance the genetic evaluation of mastitis.
Methods
Genetic evaluation of mastitis was carried out by modeling the dynamic nature of somatic cell count (SCC) within the lactation. The SCC patterns were captured by modeling transition probabilities between assumed states of mastitis and non-mastitis. A widely dispersed SCC pattern generates high transition probabilities between states and vice versa. This method can model transitions to and from states of infection simultaneously, i.e. both the mastitis liability and the recovery process are considered. A multilevel discrete time survival model was applied to estimate breeding values on simulated data with different dataset sizes, mastitis frequencies, and genetic correlations.
Results
Correlations between estimated and simulated breeding values showed that the estimated accuracies for mastitis liability were similar to those from previously tested methods that used data of confirmed mastitis cases, while our results were based on SCC as an indicator of mastitis. In addition, unlike the other methods, our method also generates breeding values for the recovery process.
Conclusions
The developed method provides an effective tool for the genetic evaluation of mastitis when considering the whole disease course and will contribute to improving the genetic evaluation of udder health.
Background
Mastitis is a common disease in dairy cattle with severe economic consequences [1]. It has been shown that susceptibility to the disease varies between breeds and individuals, with heritabilities ranging from 0.07 to 0.12 [2, 3]. Genetic evaluation of the disease is an important issue and a wide range of methods is available. Methods can be divided into cross-sectional or longitudinal approaches. Cross-sectional methods consider each lactation as a static process, whereas longitudinal methods model changes in disease states during the lactation. The variables mostly used in mastitis analyses are recorded cases of clinical mastitis (CM) or somatic cell counts (SCC). Routine recording of CM is not performed in most countries and SCC can be used as a proxy in the genetic evaluation, due to its high genetic correlation with CM [4–6].
Cross-sectional analyses of recorded cases of CM consider either CM as a binary variable, distinguishing between absence of cases and occurrence of at least one case during lactation, or the lactation average SCC. A shortcoming of these methods is that they ignore the dynamic nature of mastitis, e.g., multiple cases of mastitis or the longitudinal SCC pattern. The dynamic nature of mastitis is taken into account to various extents in longitudinal approaches, for example by counting the number of CM cases during a lactation using a Poisson model [7] or by modeling presence or absence of CM in pre-specified lactation intervals using a longitudinal threshold liability model, which gives as many heritability estimates as intervals [8, 9]. Other longitudinal methods include time in lactation either as a variable in the model, e.g., time to first case in a survival analysis [10], or as an explanatory variable in a random regression model, making heritability a function of time [11–14].
In [15–17], SCC was analyzed longitudinally, using different patterns based on deviations from the typical SCC curve, to investigate how these patterns might be related to pathogen-specific cases of CM. Traits derived from lactation SCC data have also been used to capture longitudinal features such as mastitis intensity through the number of peaks (consecutive SCC test days above a certain level), recovery pace through peak lengths, and lactation fluctuation by the standard deviation of SCC [18–20].
In longitudinal studies, the degree at which the dynamic nature of mastitis is accommodated varies, but research in this field is moving to more sophisticated methods in which the longitudinal approach is one aspect that is stressed. Results from such studies show that considering the longitudinal aspect improves the genetic evaluation of CM.
We present an alternative longitudinal approach in which genetic evaluation of mastitis is performed based on changes in SCC during lactation. These changes are captured by modeling transition probabilities between assumed states of mastitis and non-mastitis. The method simultaneously models the transition probability of developing mastitis and the probability of recovering from an infection. The former has been the focus of many studies but the infection recovery process has been little investigated. In our approach, we model both aspects to capture as much genetic information as possible from the SCC lactation pattern.
Methods
Transition probability model
which gives the probabilities of changing states or remaining in the current state. A transition between two particular states is repeatable, i.e. a cow can be in a specific state more than once during a given lactation. An episode is defined as the duration of each state; after a transition, a new episode begins, leading to multiple episodes within the lactation. The transition probabilities reflect the SCC fluctuations. A widely dispersed pattern, i.e. with many fluctuations between high and low SCC levels, will generate higher values of ${\pi}_{i}^{\left(HD\right)}$ and ${\pi}_{i}^{\left(DH\right)}$ (high transition probabilities between states) and lower values of $1-{{\pi}_{i}}^{\left(HD\right)}$ and $1-{{\pi}_{i}}^{\left(DH\right)}$ (low probabilities of remaining in a current state), compared to values for an individual with a SCC pattern similar to the average lactation curve. A desired structure of the transition matrix would be to have high values of $1-{{\pi}_{i}}^{\left(HD\right)}$ and ${\pi}_{i}^{\left(DH\right)}$ (obviously together with low values of ${\pi}_{i}^{\left(HD\right)}$ and $1-{{\pi}_{i}}^{\left(DH\right)}$), which reflect an individual who rarely develops mastitis but if it does, has a quick recovery.
When modeling survival data with repeated events, such as repeated mastitis cases during lactation, multilevel models are effective and powerful tools [22]. In these models, repeated transitions are viewed as an extra level in a higher level hierarchical structure. Here, we use a three-level structure, in which episodes are nested within cows, and herds and sires are cross-classified on the highest level. Data for mastitis and SCC are most often interval-censored, i.e. the actual time for a transition between states is unknown. When data is collected retrospectively, a state change is only known to have occurred at some point between two data collection times. Although the underlying process is continuous, the structure of the data calls for a model that recognizes its discrete nature. Therefore, a continuous time survival model such as the Cox proportional hazards model [23] was rejected in favor of a multilevel discrete time survival model. Statistical descriptions and evaluations of multilevel survival models with repeated events can be found in [24–26]. By restructuring and expanding the dataset, the multilevel binary response model can be fitted using logistic regression, as in our study, or using other standard methods for discrete response data.
with random effects ${s}_{j}^{\left(DH\right)}~N\left(0,{\left({\sigma}_{S}^{\left(DH\right)}\right)}^{2}\right)$, ${h}_{k}^{\left(DH\right)}~N\left(0,{\left({\sigma}_{h}^{\left(DH\right)}\right)}^{2}\right)$ and ${e}_{ijk}^{\left(DH\right)}~N\left(0,{\left({\sigma}_{e}^{\left(DH\right)}\right)}^{2}\right).$
The two transitions probabilities are expressed conditional on the current state. In a given interval, a cow can only change states in one direction (or stay in the same state).
Simulations and Bayesian inference
As shown by Allison [27] and Browne et al. [22], the likelihood function for a discrete time multilevel model is equivalent to the likelihood for the regression analysis of a dichotomous dependent variable. Therefore, the multilevel discrete time model can be fitted using standard software packages for logistic regression, where the response variable is the binary indicator of the occurrence of an event ${f}_{ijkt}^{\left(HD\right)}$ or ${f}_{ijkt}^{\left(DH\right)}$. The choice of inference of the model is open to both a classical or a Bayesian estimation approach, however, the modular nature of Markov chain Monte Carlo (MCMC) algorithms in Bayesian inference make them an attractive choice to estimate discrete time repeated events [22]. In addition, it has been shown that Bayesian estimates are less biased than maximum likelihood estimates for random-effect logistic regression models [28]. In this paper, Bayesian estimates were obtained with the multilevel software program MLwiN, which is developed by the Centre of Multilevel Modeling, University of Bristol [29]. This software offers the choice of classical or Bayesian estimation and has many features and options to fit multilevel models, including parameter expansion, which is discussed in the next section. The MCMC approach in MLwiN consists of iterative Metropolis-Hastings (MH) and Gibbs sampler steps to successively update the model parameters. Gibbs sampling is used for variances and univariate-update random-walk Metropolis sampling with Gaussian proposal distributions for fixed effects and residuals. For a complete description of the MCMC steps and conditional densities, see [22, 30].
Estimates of parameters in the HD direction were based on 10 000 iterations with a burn in of 500, whereas estimates in the DH direction were based on 100 000 iterations with a burn in of 5000. The strongly reduced dataset for DH estimations called for longer iteration chains to obtain convergence and to estimate variances with the same accuracies as the HD estimates.
Parameter expansion
As previously mentioned, MCMC methods can reduce estimation bias in the discrete time survival model. However, data expansion results in very large datasets for which the MCMC algorithm can be slow and generate chains that exhibit poor mixing. This was apparent for the MCMC algorithms of σ_{ e } , especially for estimation of parameters in the DH direction. One of the main causes of poor mixing in the MCMC algorithm is correlated model parameters. When correlations exist between random effects and their variances, and the variances are close to zero, the MCMC chain can get stuck close to zero, both for the variance and its random effect. Parameter expansion is a reparameterization method, which reformulates the statistical model by replacing certain parameters with others that are not correlated and therefore generate MCMC chains with much better mixing than the original chains [22, 31, 32]. The reparameterization is done in such a way that it is possible to recover the original parameters in the model. The method is a built-in option in MLwiN and is used in the estimations reported in this paper, resulting in much better mixing and therefore faster convergence and less need for long iteration chains.
Simulated data
The process used to generate data on mastitis was an extension of the simulations by Carlén et al. [33] and Schneider et al. [34]. These studies generated data on milk production, interval between calving and first ovulation, conception and mastitis liability. Mastitis history data were simulated as a binary trait until the first case of mastitis or the end of the lactation. For this study, we extended the program to include repeated cases of mastitis and also SCC data connected to every binary response.
Average values for the two mastitis frequency scenarios
Scenario 1 | Scenario 2 | |
---|---|---|
1. Mastitis lactation (%) | 16.5 | 40.0 |
2. Mastitis cases per lactation | 0.28 | 0.95 |
3. Mastitis cases per affected lactation | 1.50 | 2.40 |
Mastitis history data were simulated as a binary trait with underlying normally distributed mastitis and recovery liabilities. The mastitis and recovery liabilities for each cow were modeled by a herd effect, the animal's breeding value, and an environmental value. Breeding values for the cows were simulated by adding half the breeding value of the sire, half the breeding value of the dam and a Mendelian sampling term. The environmental value corresponds to a permanent cow effect, which means that the liability value represents a mean liability over lactation. A weekly random variable was added to allow the mastitis liability to vary from week to week. Before the weekly variations were added, the variance in the liabilities (~N(0, 0.8^{2})) was chosen to give a mastitis heritability of 0.1, which is in line with what is normally found in cross-sectional linear models (on the underlying scale) [5, 37, 38]. Together with the weekly random component (~N(0, 0.6^{2})), variances added up to a phenotypic variance of 1. The additive genetic variance was 0.036, which resulted in a simulated heritability of 0.036 for the weekly data (0.039 if the herd variance of 0.072 is excluded from the phenotypic variance). The recovery liabilities were simulated in the same way and with the same parameters.
For cases in which mastitis did not develop in the previous week, the mastitis liability of the current week generated the next binary outcome. If the resulting mastitis liability was above a defined threshold which corresponded to the targeted mastitis frequency, the cow developed mastitis. If mastitis did develop in the previous week, the current recovery liability decided the next week's binary outcome. If the recovery liability was above the threshold of -0.4 (corresponding to the relative recovery frequency of 0.65), the cow was free from the infection, and if not, the cow remained infected for (at least) another week.
The binary CM data was then used to simulate SCC observations: values for uninfected test days (days when the cow's SCC level is measured) were simulated as random deviations (σ = 0.64 according to previous studies [1]) from a baseline curve and test days with mastitis infections as random deviations from a function with instant SCC increase, followed by a successive decline down to the baseline level. The baseline curve was modeled by a three-phase linear spline function parameterized according to [39]. Another spline function expressed the immediate increase of SCC from the baseline at the time of infection and the successive decline down to the baseline level during the following weeks. This spline function was created according to findings by De Haas et al. [15] by generalizing the effects of Staphylococcus aureus, Escherichia coli, Streptococcus dysgalactiae, and Streptococcus uberis on SCC.
Data expansion
Before applying the method to the data, the sequence of SCC lactation data was converted into a sequence of binary responses, which indicated whether a transition had occurred within each time interval. Cow, sire, and herd indicators were repeated as many times as there were binary responses in the lactation. Furthermore, a time indicator was added, to number the time intervals repeatedly until a transition. After a transition, the counter started over again and ran until the next transition or the end of the lactation. There was no restriction on the cows having the same number of measurements or the same interval between measurements. In theory, a cow only needs two consecutive measurements although, in practice, very short series carry little information. Long intervals between measurements can also lead to information loss due to missed transitions. However, evaluation concerns the performance of the sires, not that of the individual cows, thus longer intervals between measurements can to some extent be compensated by more daughters per sire.
where ${\tau}_{t}$ is the time since calving for the t:th order response of cow i.
for D to H transitions (t = 1,..., T_{ i } -1). For example, a sequence of h_{ ijk } = [0010001110100], containing 13 measurements and three cases of classified mastitis (one which lasted over three measurements) for cow i, generated ${\mathbf{f}}_{ijk}^{\left(HD\right)}=\left[0100110\right]$ and ${\mathbf{f}}_{ijk}^{\left(DH\right)}=\left[10011\right]$, together with two time indicators ${\mathbf{\nu}}_{ijk}^{\left(HD\right)}=\left[1212311\right]$ and ${\mathbf{\nu}}_{ijk}^{\left(DH\right)}=\left[11231\right]$. The time indicators numbered the intervals until a transition occurred (i.e. the elements in f_{ ijk } until a 1 appears) or till the end of the sequence is reached. Since a cow could change states more than once during a lactation, the time indicators could start over several times. For the HD transitions, the counter started at 1 at the first measurement that is classified as healthy (most often the very first one) counting until a transition to a diseased state occurs. The counter started over again when or if the cow returned to a healthy state. For transitions in the other direction, the counter started when the cow developed the first case of classified mastitis and stopped when it returned to a healthy state. If the cow returned to a diseased state, the counter started over again. A cow without a case of classified mastitis did therefore not generate any data for transitions in the D to H direction. At the same time, the duration of a diseased state was usually much shorter than that of a healthy state. This resulted in a strongly reduced data-set for D to H transitions in comparison to the data in the H to D direction. To complement f_{ ijk } and v _{ ijk } , three more sequences k _{ ijk }, s_{ ijk } and h_{ ijk } were created, which contained indicators for cow, sire and herd. Each sequence was just a replicate of the same indicator duplicated as many times as the length of f_{ ijk } and v _{ ijk } .
The expanded data material for the HD direction contained ${\mathbf{f}}_{ijk}^{\left(HD\right)}$, ${\mathbf{\nu}}_{ijk}^{\left(HD\right)}$, ${\mathbf{k}}_{ijk}^{\left(HD\right)}$, ${\mathbf{s}}_{ijk}^{\left(HD\right)}$ and ${\mathbf{h}}_{ijk}^{\left(HD\right)}$ for all cows, while the data in the DH direction contained, ${\mathbf{f}}_{ijk}^{\left(DH\right)}$, ${\mathbf{\nu}}_{ijk}^{\left(DH\right)}$, ${\mathbf{k}}_{ijk}^{\left(DH\right)}$, ${\mathbf{s}}_{ijk}^{\left(DH\right)}$ and ${\mathbf{h}}_{ijk}^{\left(DH\right)}$ for all cows that had at least one case of classified mastitis.
Analysis of simulation results
Breeding values were estimated in separate analyses for the two transitions directions (HD and DH). Correlations between true breeding values (TBV) and estimated breeding values (EBV), i.e. accuracy of selection, for both HD $\left({r}_{BV}^{\left(HD\right)}\right)$ and DH $\left({r}_{BV}^{\left(DH\right)}\right)$ directions, were calculated for all combinations of population sizes, mastitis frequencies, and genetic correlations (r_{ G } ). In addition, correlations between EBV in the HD and DH directions were calculated for the same combinations in order to check if the negative, neutral and positive genetic correlations between mastitis and recovery liabilities in the simulated data are reproduced by the method.
Results
Average correlations between true and estimated breeding values
Daughters/sire | 60 | 150 | |||||||
---|---|---|---|---|---|---|---|---|---|
Mastitis frequency (cases/lactation) | 0.28 Scenario 1 | 0.95 Scenario 2 | 0.28 Scenario 1 | 0.95 Scenario 2 | |||||
Transition direction | HD | DH | HD | DH | HD | DH | HD | DH | |
r_{ G } = 0 | B(τ) = 10 × L(τ) | 0.532 | 0.249 | 0.697 | 0.387 | 0.749 | 0.381 | 0.836 | 0.550 |
(0.012) | (0.018) | (0.012) | (0.021) | (0.009) | (0.016) | (0.005) | (0.013) | ||
Cows | 24000 | 4149 | 24000 | 10151 | 60000 | 10374 | 60000 | 25450 | |
Herds | 1200 | 1095 | 1200 | 1197 | 1200 | 1187 | 1200 | 1200 | |
B(τ) = 15 × L(τ) | 0.549 | 0.271 | 0.702 | 0.436 | 0.751 | 0.410 | 0.838 | 0.619 | |
(0.007) | (0.013) | (0.012) | (0.018) | (0.009) | (0.012) | (0.005) | (0.006) | ||
Cows | 24000 | 4128 | 24000 | 10139 | 60000 | 10328 | 60000 | 25411 | |
Herds | 1200 | 1092 | 1200 | 1197 | 1200 | 1186 | 1200 | 1200 | |
r_{ G } = 0.2 | B(τ) = 10 × L(τ) | 0.589 | 0.268 | 0.684 | 0.407 | 0.760 | 0.386 | 0.834 | 0.533 |
0.0148 | (0.015) | 0.0108 | (0.014) | (0.013) | (0.028) | (0.007) | (0.016) | ||
Cows | 24000 | 4131 | 24000 | 10161 | 60000 | 10366 | 60000 | 25453 | |
Herds | 1200 | 1094 | 1200 | 1198 | 1200 | 1184 | 1200 | 1200 | |
B(τ) = 15 × L(τ) | 0.592 | 0.301 | 0.690 | 0.467 | 0.766 | 0.436 | 0.837 | 0.612 | |
(0.013) | (0.012) | (0.011) | (0.016) | (0.013) | (0.021) | (0.006) | (0.010) | ||
Cows | 24000 | 4115 | 24000 | 10148 | 60000 | 10322 | 60000 | 25425 | |
Herds | 1200 | 1093 | 1200 | 1197 | 1200 | 1183 | 1200 | 1200 | |
r_{ G } = -0.2 | B(τ) = 10 × L(τ) | 0.588 | 0.236 | 0.683 | 0.396 | 0.764 | 0.372 | 0.820 | 0.562 |
(0.017) | (0.019) | (0.013) | (0.021) | (0.011) | (0.021) | (0.015) | (0.012) | ||
Cows | 24000 | 4136 | 24000 | 10189 | 60000 | 10423 | 60000 | 25428 | |
Herds | 1200 | 1087 | 1200 | 1123 | 1200 | 1187 | 1200 | 1200 | |
B(τ) = 15 × L(τ) | 0.591 | 0.257 | 0.688 | 0.431 | 0.766 | 0.468 | 0.824 | 0.607 | |
(0.016) | (0.019) | (0.011) | (0.013) | (0.010) | (0.045) | (0.014) | (0.011) | ||
Cows | 24000 | 4118 | 24000 | 10176 | 60000 | 10408 | 60000 | 25396 | |
Herds | 1200 | 1084 | 1200 | 1195 | 1200 | 1187 | 1200 | 1200 |
Average correlations between estimated breeding values in the HD and DH direction
Daughters/sire | 60 | 150 | |||
---|---|---|---|---|---|
Mastitis frequency (cases/lactation) | 0.28 Scenario 1 | 0.95 Scenario 2 | 0.28 Scenario 1 | 0.95 Scenario 2 | |
r_{ G } = 0 | B(τ) = 10 × L(τ) | 0.048 (0.027) | 0.060 (0.025) | 0.020 (0.018) | 0.003 (0.019) |
B(τ) = 15 × L(τ) | 0.074 (0.010) | 0.172 (0.025) | 0.066 (0.017) | 0.134 (0.016) | |
r_{ G } = 0.2 | B(τ) = 10 × L(τ) | 0.125 (0.092) | 0.115 (0.005) | 0.056 (0.021) | 0.134 (0.032) |
B(τ) = 15 × L(τ) | 0.113 (0.026) | 0.167 (0.024) | 0.114 (0.019) | 0.270 (0.031) | |
r_{ G } = -0.2 | B(τ) = 10 × L(τ) | -0.005 (0.008) | -0.003 (0.026) | -0.074 (0.030) | -0.068 (0.042) |
B(τ) = 15 × L(τ) | 0.046 (0.004) | 0.079 (0.016) | -0.037 (0.022) | 0.041 (0.047) |
Discussion
Model
Mastitis has been the focus of several research projects and is considered in breeding programs in many countries (e.g., [40]). However, in the genetic evaluation of udder health, only the mastitis liabilities are taken into account, leaving aside the recovery process.
Our main objective was to evaluate a new method that simultaneously models transitions to and from states of mastitis considering both the mastitis liability and the recovery process. For example, the method can distinguish between two cows, each with mastitis once during lactation but one showing fast recovery and the other suffering from a protracted infection. Naturally, our aim is to have as few mastitis infections as possible. However, mastitis is unavoidable and a relatively frequent problem among dairy cows. Thus, the capacity for fast recovery is also of interest.
The multilevel discrete time survival model is well suited to analyze the repeated nature of mastitis data. Green et al. [41] used a similar model to investigate how cow, farm and management factors during the dry period influence the incidence of clinical mastitis after calving. However, in their study, only the first case of mastitis was considered, leaving out repeated events in the analysis. The method's ability to take repeated mastitis cases within the same lactation into account was shown in our study through higher correlations between TBV and EBV for higher mastitis frequencies.
Besides the extra genetic information that can be captured by considering repeated cases, it has another beneficial aspect. A higher mastitis frequency reduces the proportion of falsely classified cases, which in turn gives more accurate EBVs. For example, a 10% mastitis misclassification rate among a group of 100 cows with 10 real cases leads to nine falsely classified individuals (10% of the 90 that do not have mastitis). This means 47% (9 of 19) of the mastitis classifications are false. The same misclassification rate among an equally large group, but with 30 real cases, results in seven falsely classified individuals (10% of the 70 that do not have mastitis). In this case, 19% (7 of 37) of the mastitis classifications are false. A higher classification error percentage for lower mastitis frequencies also appears if the misclassification goes in the other direction, i.e. mastitis cases are falsely classified as non-mastitis.
Classification
Misclassifications are unavoidable when SCC is used as an indicator of mastitis. If the boundary between H and D is too low, high random fluctuations around "normal" SCC levels will lead to falsely classified cases of mastitis. There is a trade-off between classifying normal but elevated SCC values as infected if the level is too low and missing possible infections if the level is too high. The results in Table 2 show that the higher boundary level gives higher correlations between EBV and TBV. This is especially apparent for DH transitions, for which the reduced data was more sensitive to misclassifications than the complete data. Boundaries higher than $15\times L(\tau )$ were also tested (not reported here) but did not give an increase in accuracies. On the contrary, if the levels are too high, the results deteriorate.
In this simplified simulation study, the setting of the boundaries (and its consequences) are of course easier to determine than in a real-life dataset. However, the SCC response to a mastitis case was based on real-life observations and the boundaries should therefore be reasonably correct to use also in real-life data.
Performance
The simulated data in this study was used previously to evaluate linear models (LM), threshold models (TM) and survival analysis (SA) [33]. Confirmed cases of mastitis were modeled either as a binary variable (separating absence of mastitis and 1 or more cases during lactation) or as time to first case. Correlations between TBV and EBV in that study were 0.53-0.60 and 0.70-0.76 for 60 and 150 daughters, respectively. The mastitis lactation percentage coincides with scenario 1 in this paper, i.e. the lower mastitis frequency. The corresponding correlations for scenario 1 of the transition method are 0.55-0.59 and 0.75-0.77 (Table 2), showing that the results for the transition method are well in line with those in [33]. Considering that the transition method analyzes SCC values as an indicator of mastitis, while the previous methods analyzed confirmed mastitis cases, the results are even more promising. When the mastitis frequency was higher, the transition method generated even higher correlations, i.e. up to 0.7 for 60 daughters and 0.84 for 150 daughters. A higher mastitis frequency leads to a larger number of lactations with more than one case. More individuals with repeated events makes the transition method even more favorable compared to the other three methods. These results show good prospects for second and higher lactations.
An additional advantage of the transition method is the information generated on the recovery process. These results are not comparable to the LM, TM, and SA methods, because of their one-way approach. However, survival analysis could be used in a similar manner as the transition method. By analyzing the time between infection and recovery, the results could be compared with the recovery results from the transition method. Nevertheless, traditional survival analysis does not take repeated cases into account, neither in the DH nor in the HD direction, which means that a cow can only develop mastitis once and recover once. However, there are a few attempts to analyze repeated cases, where one suggestion is the survival score model [42].
A positive genetic correlation between mastitis and recovery liabilities has two characteristics: cows that easily develop an infection also recover easily (Profile 1), while cows with high disease resistance have a longer recovery time (Profile 2). Cows with Profile 1 have a better chance to generate accurate EBV for their sires than cows with Profile 2 because cows with Profile 1 show their true capacity of moving between states without getting stuck in one of the two states. Cows with a high resistance to the disease, cannot show there aptitude for recovery, simply because they rarely develop mastitis. They are "stuck" in the healthy state.
A negative genetic correlation between mastitis and recovery liabilities is reflected by having cows with low mastitis resistance and long recovery time (Profile 3) and cows with high disease resistance and quick recovery (Profile 4). The inertia factor is even more apparent for data with a negative genetic correlation between mastitis and recovery liabilities, because it is then present for data of both characteristics i.e. Profile 3 and 4. Cows with Profile 3 do not show their inclination for repeated events because they easily get stuck in a diseased state. Among the cows in Profile 4, their inclination for fast recovery will not be evident because these cows hardly develop mastitis. The inertia factor is probably the major reason why the zero and positive genetic correlations are correctly estimated by the method while the negative genetic correlations are not.
Mastitis and recovery liabilities could be estimated jointly in a bivariate model with a prior on the correlation between the two liabilities. In this study, we decided against this for practical reasons and because of the few transitions that occur for the same cow (0.28 (0.95) mastitis cases per lactation), together with the inertia factor, could easily make the estimates of the correlations very imprecise. Even for the high frequency scenario, half of the cows had no mastitis and it is impossible (or at least very model dependent) to determine what the recovery rate is for cows that do not develop mastitis. Genetic correlation between mastitis and recovery liability could be estimated more accurately with multiple lactations per cow and/or from covariances on a sire level. Whether a genetic correlation exists between mastitis and recovery liability is yet to be investigated but the inertia factor should be recognized in such studies.
One should remember that the method used to generate the simulated data is not logistic-normal but uses a probit model. In spite of this, the logit analysis gave estimates close to the values used in data generation. This indicates that our method is not very sensitive to model specification errors. Use of a probit model for analysis is simpler than the logit model, in the sense that is does not require a MH step in the Gibbs sampler. However, the logistic model is less sensitive to outliers for explanatory variables. In practice, the generating mechanism behind the data is unknown, making it difficult to specify the correct model for data analysis. However, the two models usually give very similar results [43].
Further developments
Different possibilities for developing and improving the method could be investigated. One possibility is to develop the method by adding a new state representing subclinical mastitis cases. We would then model transition probabilities between three possible states. The method could also be used to analyze repeated confirmed mastitis cases instead of SCC, excluding the need to classify cows based on SCC. However, complete data of confirmed or veterinary-treated cases are rather rare in field data but sometimes a combination of the two data types is available. The method could accommodate all available data in the genetic evaluation by combining SCC data with incomplete mastitis data. Generalizing the model to allow for time-dependent transition probabilities is another possibility which could be done by including a lactation stage effect in the analysis of the transitions.
Mastitis classification may be improved by going from a strict limit between states to a more flexible and realistic classification. Mixture models have been used successfully to classify mastitis types, e.g. [44, 45], but not in combination with the transition probability model. Multivariate mixture models give the possibility to classify mastitis on the basis of more than one variable and to model overlapping groups, which may improve classification even further.
Associations between pathogen-specific mastitis and SCC pattern have been demonstrated [15, 16]. Different patterns distinguish between long or short increase in SCC and also between lactations with or without recovery. Transition probabilities could be used to describe characteristic patterns of SCC and to identify pathogen-specific mastitis. The transition probability method may also be suitable for a wide range of other diseases for which individuals fall in and out of two or more states and the states are either confirmed or classified by one or more classification variables.
Conclusions
This paper presents and evaluates a novel longitudinal model for genetic evaluation of mastitis. The model captures the dynamic nature of the disease by modeling mastitis liability and by including the recovery process and repeated cases into the analysis. Although a more complete evaluation of the method is necessary, especially on field data, the results point towards a significant gain when broadening the genetic evaluation of udder health to include the whole disease course.
Declarations
Acknowledgements
This work was carried out as part of the RobustMilk project that is financially supported by the European Commission under the Seventh Research Framework Programme, Grant Agreement KBBE-211708. The content of this paper is the sole responsibility of the authors, and it does not necessarily represent the views of the Commission or its services. The authors thank the reviewers for their many suggestions that contributed to a much improved manuscript.
Authors’ Affiliations
References
- Nielsen C: Economic impact of mastitis in dairy cows. Doctoral Thesis No. 2009:29. 2009, Uppsala, Sweden: Swedish University of Agricultural SciencesGoogle Scholar
- Heringstad B, Rekaya R, Gianola G, Klemetsdal G, Weigel KA: Genetic change for clinical mastitis in Norwegian cattle: a threshold model analysis. J Dairy Sci. 2003, 86: 369-375. 10.3168/jds.S0022-0302(03)73615-7.View ArticlePubMedGoogle Scholar
- Lund MS, Jensen J, Petersen PH: Estimation of genetic and phenotypic parameters for clinical mastitis, somatic cell production deviance, and protein yield in dairy cattle using Gibbs sampling. J Dairy Sci. 1999, 82: 1045-1051. 10.3168/jds.S0022-0302(99)75325-7.View ArticlePubMedGoogle Scholar
- Heringstad B, Klemetsdal G, Ruane J: Selection for mastitis resistance in dairy cattle - a review with focus on the situation in the Nordic countries. Livest Prod Sci. 2000, 64: 95-106. 10.1016/S0301-6226(99)00128-1.View ArticleGoogle Scholar
- Carlén E, Strandberg E, Roth A: Genetic parameters for clinical mastitis, somatic cell score, and production in the first three lactations of Swedish Holstein cows. J Dairy Sci. 2004, 87: 3062-3070. 10.3168/jds.S0022-0302(04)73439-6.View ArticlePubMedGoogle Scholar
- Vallimont JE, Dechow CD, Sattler CG, Clay JS: Heritability estimates associated with alternative definitions of mastitis and correlations with somatic cell score and yield. J Dairy Sci. 2009, 92: 3402-3410. 10.3168/jds.2008-1229.View ArticlePubMedGoogle Scholar
- Rodrigues-Motta M, Gianola D, Heringstad B, Rosa GJM, Chang YM: A zero-inflated poisson model for genetic analysis of the number of mastitis cases in Norwegian Red cows. J Dairy Sci. 2007, 90: 5306-5315. 10.3168/jds.2006-898.View ArticlePubMedGoogle Scholar
- Heringstad B, Chang YM, Gianola D, Klemetsda G: Genetic analyis of longitudinal trajectory of clinical mastitis in first-lactation Norwegian cattle. J Dairy Sci. 2003, 86: 2676-2683. 10.3168/jds.S0022-0302(03)73863-6.View ArticlePubMedGoogle Scholar
- Heringstad B, Chang YM, Gianola D, Klemetsdal G: Multivariate threshold model analysis of clinical mastitis in multiparous Norwegian dairy cattle. J Dairy Sci. 2004, 87: 3038-3046. 10.3168/jds.S0022-0302(04)73436-0.View ArticlePubMedGoogle Scholar
- Carlén E, Schneider M, del P, Strandberg E: Comparison between linear models and survival analysis for genetic evaluation of clinical mastitis in dairy cattle. J Dairy Sci. 2005, 88: 797-803. 10.3168/jds.S0022-0302(05)72744-2.View ArticlePubMedGoogle Scholar
- Rekaya R, Gianola D, Weigel K, Shook G: Longitudinal random effects model for genetic analysis of binary data with application to mastitis in dairy cattle. Genet Sel Evol. 2003, 35: 457-468. 10.1186/1297-9686-35-6-457.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang Y-M, Gianola D, Heringstad B, Klemetsdal G: Longitudinal analysis of clinical mastitis at different stages of lactation in Norwegian cattle. Livest Prod Sci. 2004, 88: 251-261. 10.1016/j.livprodsci.2003.11.007.View ArticleGoogle Scholar
- Ødegård J, Jensen J, Klemetsdal P, Madsen P, Heringstad B: Genetic analysis of somatic cell score in Norwegian cattle using random regression test-day models. J Dairy Sci. 2003, 86: 4103-4114. 10.3168/jds.S0022-0302(03)74024-7.View ArticlePubMedGoogle Scholar
- Carlén E, Grandinson K, Emanuelson U, Strandberg E: Random regression models for genetic evaluation of clinical mastitis in dairy cattle. Animal. 2009, 3: 1100-1108. 10.1017/S1751731109004601.View ArticlePubMedGoogle Scholar
- de Haas Y, Barkema HW, Veerkamp RF: The effect of pathogen specific clinical mastitis on the lactation curve for somatic cell count. J Dairy Sci. 2002, 85: 1314-1323. 10.3168/jds.S0022-0302(02)74196-9.View ArticlePubMedGoogle Scholar
- de Haas Y, Veerkamp RF, Barkema HW, Gröhn YT, Schukken YH: Associations between pathogen-specific cases of clinical mastitis and somatic cell count patterns. J Dairy Sci. 2004, 87: 95-105. 10.3168/jds.S0022-0302(04)73146-X.View ArticlePubMedGoogle Scholar
- de Haas Y, Barkema HW, Schukken YH, Veerkamp RF: Associations between somatic cell count patterns and the incidence of clinical mastitis. Prev Vet Med. 2005, 67: 55-68. 10.1016/j.prevetmed.2004.10.006.View ArticlePubMedGoogle Scholar
- de Haas Y, Ouweltjes W, ten Napel J, Windig J, de Jong G: Alternative somatic cell count traits as mastitis indicators for genetic selection. J Dairy Sci. 2008, 91: 2501-2511. 10.3168/jds.2007-0459.View ArticlePubMedGoogle Scholar
- Urioste JI, Franzén J, Strandberg E: Phenotypic and genetic characterization of novel somatic cell count traits from weekly or monthly observations. J Dairy Sci. 2010, 93: 5930-5941. 10.3168/jds.2010-3301.View ArticlePubMedGoogle Scholar
- Windig JJ, Ouwelties W, ten Napel J, de Jong G, Veerkamp RF, de Haas Y: Combining somatic cell count traits for optimal selection against mastitis. J Dairy Sci. 2010, 93: 1690-1701. 10.3168/jds.2009-2052.View ArticlePubMedGoogle Scholar
- Hagnestam C, Emanuelson U, Berglund B: Yield losses associated with clinical mastitis occurring in different weeks of lactation. J Dairy Sci. 2007, 90: 2260-2270. 10.3168/jds.2006-583.View ArticlePubMedGoogle Scholar
- Browne WJ, Steele F, Golalizadeh M, Green MJ: The use of simple reparameterizations to improve the efficiency of Markov chain Monte Carlo estimation for multilevel models with applications to discrete time survival models. J R Stat Soc Ser A Stat Soc. 2009, 172: 579-598. 10.1111/j.1467-985X.2009.00586.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Kleinbaum DG, Klein M: Survival analysis: A self-Learning text. 2005, New York: Springer Science+Business Media IncGoogle Scholar
- Goldstein H, Pan H, Bynner J: A flexible procedure for analyzing longitudinal event histories using a multilevel model. Understand stat. 2004, 3: 85-99. 10.1207/s15328031us0302_2.View ArticleGoogle Scholar
- Steele F, Golstein H, Browne W: A general multilevel multistate competing risk model for event history data, with an application to a study of contraceptive use dynamics. Stat Model. 2004, 4: 145-159. 10.1191/1471082X04st069oa.View ArticleGoogle Scholar
- Steele F: Multilevel models for longitudinal data. J R Stat Soc Ser A Stat Soc. 2008, 171: 5-19.Google Scholar
- Allison PD: Discrete-time methods for the analysis of event histories. Sociological methodology. Edited by: Leinhardt S. 1982, San Fransisco: Jossey-Bass, 61-98.Google Scholar
- Browne WJ, Draper D: A comparison of Bayesian and likelihood methods for fitting multilevel models. Bayes Anal. 2006, 1: 473-550.View ArticleGoogle Scholar
- Centre for Multilevel Modelling (CMM): [http://www.cmm.bristol.ac.uk]
- Browne WJ: Applying MCMC Methods to Multilevel Models. PhD dissertation. 1998, UK: Department of Mathematical Sciences, University of BathGoogle Scholar
- Van Dyk DA, Meng XL: The art of data augmentation. J Comput Graph Stat. 2001, 10: 1-50. 10.1198/10618600152418584.View ArticleGoogle Scholar
- Gelman A, van Dyk DA, Huang ZY, Boscardin WJ: Using redundant parameterizations to fit hierarchical models. J Comput Graph Stat. 2008, 17: 95-122. 10.1198/106186008X287337.View ArticleGoogle Scholar
- Carlén E, Emanuelson U, Strandberg E: Genetic evaluation of mastitis in dairy cattle using linear models, threshold models, and survival analysis: a simulation study. J Dairy Sci. 2006, 89: 4049-4057. 10.3168/jds.S0022-0302(06)72448-1.View ArticlePubMedGoogle Scholar
- Schneider M, del P, Strandberg E, Ducrocq V, Toth A: Survival analysis applied to genetic evaluation for female fertility in dairy cattle. J Dairy Sci. 2005, 88: 2253-2259. 10.3168/jds.S0022-0302(05)72901-5.View ArticleGoogle Scholar
- Svensk Mjölk 2010, Husdjursstatistik, [Cattle statistics]. [http://www.svenskmjolk.se/Statistik/Mjolkforetaget], Accessed Nov 2010
- Svensk Mjölk [Swedish Dairy Association] 2009, Avelsvärden för mjölkrastjurar [Breeding values for dairy bulls]. Svensk Mjölk SE-631:84Google Scholar
- Heringstad B, Karlsen A, Klemetsdal G, Ruane J: Preliminary results from a genetic analysis of clinical mastitis. Proc. Int. Workshop on genetic improvement of functional traits in cattle. Interbull Bull. 1997, 15: 45-49.Google Scholar
- Kadarmideen HN, Thompson R, Simm G: Linear and threshold model genetic parameters for disease, fertility and milk production in dairy cattle. Anim Sci. 2000, 71: 411-419.Google Scholar
- Hagnestam-Nielsen C, Østergaard S: Economic impact of clinical mastitis in a dairy herd assessed by stochastic simulations using different methods to model yield losses. Animal. 2009, 3: 315-328. 10.1017/S1751731108003352.View ArticlePubMedGoogle Scholar
- Miglior F, Muir BL, van Doormaal BJ: Selection indices in Holstein cattle of various countries. J Dairy Sci. 2005, 88: 1255-1263. 10.3168/jds.S0022-0302(05)72792-2.View ArticlePubMedGoogle Scholar
- Green MJ, Bradley AJ, Medley GF, Browne WJ: Cow, farm, and management factors during the dry period that determine the rate of clinical mastitis after calving. J Dairy Sci. 2007, 90: 3764-3776. 10.3168/jds.2007-0107.PubMed CentralView ArticlePubMedGoogle Scholar
- Veerkamp RF, Brotherstone S, Engel B, Meuwissen THE: Analysis of censored survival data using random regression models. Anim Sci. 2001, 72: 1-10.Google Scholar
- Dobson AJ, Barnett AG: Binary variables and logistic regression. An introduction to generalized linear models. 2008, Boca Raton: Chapman and Hall, 123-147. 3Google Scholar
- Boettcher PJ, Caraviello D, Gianola D: Genetic analysis of somatic cell scores in US Holsteins with a Bayesian mixture model. J Dairy Sci. 2007, 90: 435-443. 10.3168/jds.S0022-0302(07)72645-0.View ArticlePubMedGoogle Scholar
- Ødegård J, Madsen P, Gianola D, Klemetsdal G, Jensen J, Heringstad B, Korsgaard IR: A Bayesian threshold-normal mixture model for analysis of a continous mastitis related trait. J Dairy Sci. 2005, 88: 2652-2659. 10.3168/jds.S0022-0302(05)72942-8.View ArticlePubMedGoogle Scholar
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.