|
|
||||||||
Department of Dairy Science, University of Wisconsin, Madison 53706
Corresponding author: Kent Weigel; e-mail: kweigel{at}wisc.edu.
| ABSTRACT |
|---|
|
|
|---|
Key Words: longevity somatic cell count mastitis survival analysis
Abbreviation key: CM = clinical mastitis, PL = length of productive life, RR = relative risk of culling
| INTRODUCTION |
|---|
|
|
|---|
Individual cow SCC values are routinely used as a tool for involuntary culling decisions on commercial dairies, and PTA of dairy sires for SCS are commonly used for genetic selection decisions and routinely incorporated into economic indices, such as Lifetime Net Merit. Estimates of the genetic correlation between SCC and clinical mastitis (CM) are generally large and positive. For example, Rupp and Boichard (1999) reported an estimated correlation of 0.72 in French Holsteins. In Israeli Holsteins, Weller et al. (1992) reported an estimated genetic correlation of 0.994 between SCC and the incidence of bacterial subclinical infection; however, the estimated correlation between SCC and producer-recorded CM was only 0.299. Cranford and Pearson (2001) noted that a strong relationship exists between sires PTA for SCS and the incidence rate of CM in their daughters. Furthermore, SCC is routinely recorded on all cows in an objective manner, using a continuous scale of measurement. On the contrary, CM is scored in a binary manner, and both the diagnosis and (decision to apply) treatment are subjective. For these reasons, heritability estimates for SCC are generally greater in magnitude than corresponding estimates for CM. For example, heritability estimates reported by Rupp and Boichard (1999) were 0.17 for SCC and 0.02 for CM.
The longevity or survival rate of dairy cows can be influenced by SCC through the death or culling of clinically affected animals, as well as the culling of subclinical animals with high SCC (to achieve milk quality premiums). Cranford and Pearson (2001) found significant, unfavorable correlations between sire PTA for SCS and the number of lactations, total (lifetime) DIM, and length of productive life (PL) of their daughters. The magnitude of these relationships increased considerably when sire PTA for SCS were adjusted to a constant level of PTA for milk yield. For example, a one-point increase in sire PTA for SCS corresponded to a decrease of 87 d in PL before standardization to a constant level of milk yield, and 132 d after this standardization.
A potential complication in statistical analyses is that SCC is a time-dependent variable. The SCC of an individual cow changes over time, as does the likelihood that its current SCC level will lead to removal from the herd. Survival analysis methodology effectively uses information from time-dependent covariates and is capable of handling censored observations, so data from animals that are still alive at the time of analysis can be used (Ducrocq and Sölkner, 1998a; Vukasinovic, 1999). Furthermore, the distribution of longevity data is often skewed, and methods based on the assumption of normality have limited applicability to the analysis of longevity data (Egger-Danner, 1993).
Gröhn et al. (1998) used a Cox proportional hazards model, with time-dependent covariates, to assess the impact of CM and other diseases on culling rates in New York State. The number of days from first calving until culling or censoring was reduced significantly by a CM infection, regardless of whether the infection occurred in early, middle, or late lactation. However, the influence of CM on the risk of culling was influenced only minimally by the inclusion of milk yield as a covariate.
Likewise, Neerhof et al. (2000) used survival analysis methodology to investigate the impact of mastitis on functional longevity in Danish Black and White cattle. These authors noted the importance of considering the long-term impact of CM on culling decisions; models in which a CM episode was assumed to influence longevity until the end of the lactation provided the best fit to the data. After adjusting for parity, stage of lactation, previous lactation milk yield, age at first calving, and proportion of Holstein genes, the relative risk of culling for a cow with CM was 1.69 times that of an uninfected cow. Among bulls with at least 100 uncensored daughters, the product-moment correlation between sires genetic evaluations for CM and the relative risk of culling (RR) of their daughters was 0.40.
The objective of this study was to use individual cow SCC data to estimate the impact of clinical and subclinical mastitis on the functional longevity of US Holstein and Jersey cows using survival analysis methodology.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Longevity, also commonly known as herd life or PL, was defined as the number of days from first calving until culling or censoring. Records from cows that were sold for dairy purposes were considered as censored, as were records from cows that resided in herds that discontinued milk recording, cows that were still alive after 5 completed lactations, and cows that were still alive when the data were extracted. Cows that completed a 305-d lactation but did not calve again within 6 mo were considered dead and were treated as uncensored. Functional longevity is the ability to delay involuntary culling, because voluntary culling for low production can be an important reason for disposal. Functional longevity can be approximated by correcting true longevity for production level (Ducrocq and Sölkner, 1998a). Therefore, we created quintiles for within herd-year 305-d mature equivalent milk yield (in Jerseys) or 305-d mature equivalent fat plus protein yield (in Holsteins, for which low component percentages are more likely to be a problem).
The statistical model for analyzing the impact of SCC on functional longevity was as follows:
![]() |
where hijklm(t) = hazard function (instantaneous probability of culling) for a given cow at t days since first calving; h0(t) = Weibull baseline hazard function with scale parameter
and shape parameter
; Pi(t1) = time-dependent fixed parity-stage of lactation effect, assumed to be piecewise constant with change points at t1 = 0, 45, and 270 d postpartum in lactations 1, 2, 3, 4, and 5; Aj = time-independent effect of age at first calving, treated as a continuous covariate with regression coefficient ß; Mk(t2) = time-dependent effect of within-herd-year quintile ranking for mature equivalent 305 d milk yield (Jerseys) or combined fat and protein yield (Holsteins), assumed to be piecewise constant with change points at t2 = calendar dates of calving in lactations 1, 2, 3, 4, and 5; hysl(t2) = time-dependent random effect of herd-year-season, assumed to be independently distributed, following a log-gamma distribution with parameter
, and assumed to be piecewise constant with change points at t2 = January 1, May 1, and September 1 of each calendar year; SCCm(t2) = time-dependent effect of lactation average SCC, grouped into 15 classes (to nearest 50,000 cells/ mL) and assumed to be piecewise constant with changes at t2 = calendar dates of calving in lactations 1, 2, 3, 4, and 5.
This model was applied to data from each category of herds independently, such that separate estimates were obtained for
(the shape parameter of the Weibull distribution) and
(the parameter of the log-gamma distribution of herd-year-season effects). In addition to overall culling, the DHI culling codes corresponding to individual cows were used to identify cows that (according to the herd owner) were culled due to mastitis. We subsequently repeated each analysis, but in this case longevity records were considered completed only if cows were culled due to mastitis; longevity records from cows that were culled for any other reason were considered censored. In this manner, it was possible to determine the impact of SCC on cows overall risk of culling, as well as the impact of SCC on cows risk of culling for mastitis. One must recognize, however, that assignment of reasons for culling by dairy farmers is often imprecise. For example, a nonpregnant cow with poor milk production may be culled after a CM episode, but the farmer could report the corresponding reason for culling as mastitis, low production, or infertility. The present study used lactation average SCC data for each cow, rather than individual test-day SCC records, primarily due to computational limitations, concern about month-to-month variation in test-day SCC measurements (e.g., if the DHI tester happens to arrive on the day a particular cow has CM), and lack of monthly SCC data for herds enrolled in labor-efficient milk recording programs (e.g., herds with bimonthly or quarterly SCC options). However, it is likely that test-day SCC data would have provided more precise information about the udder health status of individual cows at specific times. The Survival Kit Version 3.12, a set of FORTRAN programs by Ducrocq and Sölkner (1998b), was used for the analysis. Cows with lactation average SCC of 200,000 to 250,000 cells/mL were considered "average", and their RR estimates were constrained to unity.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
The relative risk of overall culling or culling due to mastitis for Holstein cows with varying lactation average SCC in low (111,000 to 279,000 cells/mL), medium-low (279,000 to 305,000 cells/mL), medium (305,000 to 329,000 cells/mL), medium-high (329,000 to 358,000 cells/mL), or high (358,000 to 540,000 cells/mL) SCC Holstein herds is shown in Table 2
. The relative risk of culling for individual Holstein cows with high SCC, relative to "average" Holstein cows with lactation SCC values of 200,000 to 250,000 cells/mL and risk ratios of 1.0, was higher in herds with (otherwise) low average SCC. For example, Holstein cows with lactation average SCC > 700,000 cells/mL were roughly 4.4 times more likely to be culled than average cows in herds with low SCC. On the other hand, Holstein cows with lactation average SCC > 700,000 cells/mL were only 2.7 times more likely to be culled than average cows in herds with high SCC. Likewise, the relative risk of culling for Holstein cows with lactation average SCC of 600,000 to 650,000 cells/mL ranged from 2.0 in herds with low average SCC to approximately 1.6 in herds with high average SCC. Interestingly, Holstein cows with extremely low lactation average SCC (<50,000 cells/mL) had a greater relative risk of culling than average cows, particularly in herds with medium-high or high SCC.
|
The relative risk of overall culling and culling due to mastitis for Jersey cows with varying lactation average SCC in low (139,000 to 306,000 cells/mL), medium (306,000 to 349,000 cells/mL), or high (349,000 to 491,000 cells/mL) SCC Jersey herds is shown in Table 3
. Relative to Jersey cows with a lactation average SCC of 200,000 to 250,000 cells/mL, Jersey cows with high lactation average SCC had a much greater risk for culling, particularly in herds with low average SCC. For example, Jersey cows with lactation average SCC > 700,000 cells/mL had RR of 4.7 in herds with low SCC (compared with Jersey cows with SCC of 200,000 to 250,000 cells/mL) and RR of 2.7 in herds with high SCC. As in Holsteins, the risk of culling for Jersey cows with lactation average SCC of 600,000 to 650,000 cells/mL, as compared with average cows, ranged from 1.9 in low SCC herds to 1.6 in high SCC herds. This seems to indicate more stringent culling of cows with clinical or subclinical mastitis in herds that (apart from this cow) tend to have few mastitis problems. However, a similar pattern could be observed if a constant percentage of cows were culled due to CM or high SCC in all categories. For example, the poorest 10% of cows for SCC in a herd with excellent udder health would have lower average SCC than the poorest 10% of cows in herd with a significant mastitis problem, and culling the same proportion from each herd would necessitate more stringent culling for mastitis in the former (based on mean SCC of cows that were removed). As in Holsteins, individual Jersey cows that had extremely low lactation average SCC (e.g., <50,000 cells/mL) appeared to be at greater risk for culling than cows in the reference class (with SCC of 200,000 to 250,000 cells/mL).
|
Previous authors (e.g., Coffey et al., 1986; Shukken et al., 1994) have expressed concern that continuously decreasing SCC through genetic selection could lead to cows with impaired ability to recruit leukocytes and adequately respond to an IMI. Our results in Holsteins and Jerseys appear to indicate that cows with extremely low SCC may suffer from a reduced capacity to resist mastitis, particularly in herds with poor environmental conditions or poor udder health management practices (and, hence, high likelihood of exposure to mastitis pathogens). However, it is important to note (as shown in Figures 1
and 2
) that the number of cows in the lowest SCC category was limited, particularly in herds with high average SCC. The results of Rupp and Boichard (2000), who compared SCC in the first month of first lactation with the subsequent time to CM infection, contradicted the results of the present study. In French Holsteins, the probability of CM increased continuously as the initial SCC level increased, and this pattern occurred in both low SCC and high SCC herds. These authors concluded that low initial SCC values are desirable, and that no "intermediate optimum" exists, with respect to SCC.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication December 8, 2003. Accepted for publication August 29, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. J. Wilson, Y. T. Grohn, G. J. Bennett, R. N. Gonzalez, Y. H. Schukken, and J. Spatz Milk Production Change Following Clinical Mastitis and Reproductive Performance Compared Among J5 Vaccinated and Control Dairy Cattle J Dairy Sci, October 1, 2008; 91(10): 3869 - 3879. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. D. Bannerman, A. C. W. Kauf, M. J. Paape, H. R. Springer, and J. P. Goff Comparison of Holstein and Jersey Innate Immune Responses to Escherichia coli Intramammary Infection J Dairy Sci, June 1, 2008; 91(6): 2225 - 2235. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Ferguson, G. Azzaro, M. Gambina, and G. Licitra Prevalence of Mastitis Pathogens in Ragusa, Sicily, from 2000 to 2006 J Dairy Sci, December 1, 2007; 90(12): 5798 - 5813. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sewalem, F. Miglior, G. J. Kistemaker, and B. J. Van Doormaal Analysis of the relationship between somatic cell score and functional longevity in canadian dairy cattle. J Dairy Sci, September 1, 2006; 89(9): 3609 - 3614. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Reksen, L. Solverod, A. J. Branscum, and O. Osteras Relationships between milk culture results and treatment for clinical mastitis or culling in Norwegian dairy cattle. J Dairy Sci, August 1, 2006; 89(8): 2928 - 2937. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |