|
|
||||||||
1 Department of Reproduction, Obstetrics and Herd Health, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
2 Department of Health Management, Atlantic Veterinary College, University of Prince Edward Island, Charlottetown, Canada
3 Department of Physiology, Biochemistry and Biometrics, Faculty of Veterinary Medicine, Ghent University, Merelbeke, Belgium
Corresponding author: Sarne De Vliegher; e-mail: Sarne.Devliegher{at}UGent.be.
| ABSTRACT |
|---|
|
|
|---|
50,000 cells/mL.
Key Words: Cox frailty model culling dairy heifer early lactation somatic cell count
Abbreviation key: HR = hazard ratio, LnSCC = natural log-transformed SCC, LnSCCel = natural log-transformed SCCel, MY = milk yield at test-day (kg; measured after 14 DIM), SCCel = SCC in early lactation (measured between 5 and 14 DIM).
| INTRODUCTION |
|---|
|
|
|---|
Elevated SCC in early lactation (SCCel) in heifers, suggesting presence of IMI around calving, is associated with elevated test-day SCC and higher probabilities of test-day SCC >200,000 cells/mL (De Vliegher et al., 2004a), and with an increased probability of clinical mastitis during the first lactation (Rupp and Boichard, 2000). In addition, elevated SCCel is associated with lower milk production (Coffey et al., 1986; De Vliegher et al., 2005). Furthermore, clinical and subclinical mastitis early postpartum have negative effects on reproductive performance (Barker et al., 1998; Schrick et al., 2001). As a result, dairy herds with a large number of heifers calving with infected udder quarters will suffer substantial economical losses. Moreover, elevated SCCel in heifers could be associated with an increased culling hazard during the first lactation, possibly to some extent because of the aforementioned effects. Disease (including mastitis) has direct and indirect effects on culling. The indirect effects may be reflected by, for instance, milk yield (MY) (Gröhn et al., 1997), as most diseases cause a decline in MY, either temporarily or longer lasting (Gröhn et al., 1998). Comparing models with and without MY helps to estimate the direct and indirect effects of mastitis in general on culling (Rajala-Schultz and Gröhn, 1999b). Survival analysis is often used to assess the effect of covariates that are measured only once (time-independent covariates, e.g., the effect of SCCel in this study). Some covariates, however, are changing over time (e.g., SCC and MY at the different test-days throughout lactation) and can easily be incorporated into the semiparametric Cox models as time-depending or time-varying covariates (Gröhn et al., 1997).
The objectives of the study were 2-fold: 1) to examine the association between SCCel of heifers and culling during the first lactation while accounting for the day of assessment of SCCel and for the variability between herds, and 2) to determine what part of the effect of SCCel on culling acts indirectly through increased test-day SCC and decreased MY, by adjusting the models for test-day SCC and MY, modeled as time-varying covariates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Two data sets were created based on the full data set containing data from 13,835 heifers. The structure of the 2 data sets is given in Table 1
. An event was defined as culling for all reasons combined in the first data set (CULLALLfull; 3204 events), and as culling for udder disorders only in the second data set (CULLUDDfull; 325 events). A similar approach was followed for a subset of data of 7596 heifers with a second test-day SCC (measured after 14 DIM and before 75 DIM) of
50,000 cells/mL. In the first data set, an event was defined as culling for all reasons combined (CULLALLsub; 1587 events), and in the second data set an event was defined as culling for udder disorders only (CULLUDDsub; 136 events).
|
ijk), where tijk represents the time of culling, and
ijk was equal to 0 if the eifer was censored and to 1 in the case of culling. Two x7 different models were fitted using CULLALLfull and CULLUDDfull, respectively, with breed incorporated as a fixed effect in all models. Both natural log-transformed test-day SCC (LnSCC) and test-day MY (measured between 15 and 365 DIM) were considered time-varying covariates when included in the models. To adjust for the fact that LnSCCel was not assessed at the same day after calving for each heifer, the Cox models were stratified according to DIM (5 to 14) on which LnSCCel was measured in the period called "early lactation". This means that for each DIM value (each stratum) another baseline hazard function was assumed.
In the first model, LnSCCel was introduced assuming a constant effect over the different DIM in early lactation (model 1). In the second model, the same relationship between LnSCCel and the culling hazard was studied, but a different relationship between LnSCCel and the culling hazard was allowed according to DIM on which LnSCCel was recorded (model 2). These 2 models were compared with each other based on the likelihoodratio test. Models 3 and 4 contained LnSCC and MY, respectively. Models 5 and 6 contained both LnSCCel (as in model 2) and LnSCC, and LnSCCel (as in model 2) and MY, respectively, to evaluate the changes in LnSCCel when the time-varying covariates were included separately. Model 7 was the full model containing LnSCCel (as in model 2), LnSCC, and MY.
Two additional models (similar to model 2) were fitted using CULLALLsub and CULLUDDsub, respectively, to study the association between the LnSCCel and the culling hazard in the group of heifers with a very low second test-day SCC (
50,000 cells/mL), suggesting no udder health disorders present at that time.
Hazard ratios (HR) with 95% confidence intervals were obtained for all covariates. All models were fitted with S-Plus 6.0 for Windows (Insightful Corp., Seattle, WA).
| RESULTS |
|---|
|
|
|---|
|
50,000 cells/mL were culled, whereas 7% of heifers with a SCCel >1,000,000 cells/mL were culled. At 200 DIM, this was 7 and 13%, respectively. The same trends are present in the heifers that were culled for udder disorders (Figure 2
|
|
2 = 41.2, df = 9, P <0.001) An increase in LnSCCel was associated with an increased culling hazard, except for DIM 5. The association was significant from DIM 10 and onwards (except for DIM 11) (Table 3
2 = 12.8, df = 9, P = 0.17), but all HR were >1. Black Holstein-Friesian heifers were culled less frequently compared with Belgian White-Blue double-purpose heifers and heifers of unknown breed (Tables 3
|
|
Higher MY protected heifers against culling (Tables 3
and 4
, model 4). The magnitude of LnSCCel at the different DIM in early lactation was slightly reduced when MY was taken into account (Tables 3
and 4
, model 6 vs. 2). The changes were smaller compared with the changes due to introducing LnSCC. In addition, incorporating MY into the models took away the breed effect: Black Holstein-Friesian heifers were no longer protected from culling (e.g., Tables 3
and 4
, model 4 vs model 2). The most elaborate model presents the estimates of LnSCCel adjusted for both LnSCC and MY(Tables 3
and 4
, model 7).
Studying the association between LnSCCel and culling in heifers with a second test-day SCC
50,000 cells/mL (CULLALLsub and CULLUDDsub) revealed that, although the association was smaller, an elevated SCCel was, in general, still associated with a higher culling hazard, especially for SCCel measured in the second part of the period called "early lactation" (Table 5
).
|
| DISCUSSION |
|---|
|
|
|---|
Because SCCel decreases substantially in the first 2 wk after calving (Dohoo, 1993, Laevens et al., 1997; De Vliegher et al., 2001, 2004a, De Vliegher et al., b), models that stratified for day of assessment of SCCel were fitted. This implied that no conclusion on a DIM effect on the culling hazard could be drawn, but as no biologically relevant effect could be expected, this was not a problem. Based on our previous findings, an effect was expected related to the infection dynamics during the early lactation period as discussed earlier (De Vliegher et al., 2004a; De Vliegher et al., 2005). The most prevalent group of mastitis pathogens associated with IMI in heifers at parturition are the coagulase-negative staphylococci that are transient in nature (Oliver and Mitchell, 1983), probably because they are colonizing the teat canal rather than the mammary gland and are washed out during milking. In addition, a high rate of spontaneous IMI elimination occurs (Oliver and Jayarao, 1997). Staphylococcus aureus, on the other hand, can also be associated with IMI in heifers at calving (Fox et al., 1995), but tends to persist in lactation (Roberson et al., 1994).
Some of the factors that needed to be evaluated or adjusted for (e.g., test-day SCC and test-day MY) change over the course of lactation. To assess their immediate effect, they were introduced as time-varying covariates, because summary measures (e.g., lactation average SCC or MY) are not able to determine the effect of a covariate throughout lactation (Gröhn et al., 1997). Furthermore, the clustering effects due to the fact that heifers belong to different herds had to be taken into account by introducing a random herd effect in the Cox model, rather than including herd as a fixed effect because the individual farm is not of interest by itself (Duchateau and Janssens, 2004).
Considering only the heifers that were culled for udder health reasons (CULLUDD) vs. all culled heifers (CULLALL) increased the magnitude of association between SCCel and SCC, and the culling hazard. This is comprehensible, as the effect is not diluted by the other reasons why heifers are culled. We wanted to present both approaches because the DHI program only allows farmers to identify one culling reason per heifer, whereas farmers usually consider many factors when deciding to cull an animal (Bascom and Young, 1998). In addition, for 40.5% of the culled heifers, no specific reason was available, although some of them were probably culled because of udder health problems. Still, even when considering all culling reasons, an elevated SCCel predicted a higher culling hazard, confirming the negative economic consequences of heifer mastitis at freshening, even though the HR was only significantly >1 when recorded after 9 DIM. Confidence intervals around the HR were wider for heifers culled for udder health disorders (CULLUDD) compared with all culled heifers (CULLALL). Therefore, fewer HR significantly differed from 1, which also occurred in the analyses in the subsets of data (CULLALLsub and CULLUDDsub). This is comprehensible as the power in survival analysis is a function of the number of events and not of the number of observations (Freedman, 1982).
Beaudeau et al. (1995) included the potential 305-d mature equivalent milk production rather than actual MY in the survival models when studying the effect of disease on culling in French dairy cows. This approach avoided inclusion of part of the impact of disease on culling through their effect on cumulative milk production. However, high-yielding cows, even if they are diseased, are more likely to be kept in the herd (Gröhn et al., 1998). Comparing models with and without MY, can therefore help to estimate the direct and indirect effects of mastitis or disease in general on culling (Rajala-Schultz and Gröhn, 1999b). Hence, MY was included in the models to find out whether this changed the magnitude of the effect of SCCel on culling. Acute mastitis in the first month of lactation had a significant effect on culling throughout lactation (Rajala-Schultz and Gröhn, 1999a), but adding MY to the model, in general, reduced the effect (Rajala-Schultz and Gröhn, 1999b). The same was true in our study, indicating that a small part of the effect of an elevated SCCel on test-day SCC was mediated through MY. Moreover, introducing MY in the models influenced the breed effect. In the models without MY, Black Holstein-Friesian were protected from culling, but this effect was due to their higher MY. According to Rajala-Schultz and Gröhn (1999b), the influence of MY on the culling decision also depends on the lactation stage. We, however, did not distinguish whether the effect differed with differing stages of lactation, as this was not the primary aim of the study.
Test-day SCC was even more associated with culling than was SCCel. Fitting a model with both LnSCCel and LnSCC showed that part of the LnSCCel effect acted through the associated test-day SCC. This was not unexpected as an elevated SCCel increases the odds on elevated test-day SCC (De Vliegher et al., 2004a). Fitting a model containing LnSCCel, LnSCC, and MY results in estimates for LnSCCel that take into account the SCC effect and adjusts for the protective MY effect. Elevated SCCel is associated with elevated test-day SCC and lower MY, thus part of the association between SCCel and culling was mediated through SCC and MY.
Heifers with excellent udder health at the second test-day but with an elevated SCCel were still more at risk for being culled in their first lactation compared with heifers with an equally low second-test-day SCC but a lower SCCel. Most probably this finding is related to the fact that the latter heifers will have fewer test-day SCC >200,000 cells/mL (De Vliegher et al., 2004a) and will out-produce the heifers with the higher SCCel (De Vliegher et al., 2005). This suggests that prevention against elevated SCCel, especially in the second part of the early lactation period (as we defined it), should be preferred over treating an elevated SCCel. Additionally, some of the heifers with (sub)clinical mastitis in early lactation will lose milk production in the affected quarter(s) and will consequently have a low second test-day SCC.
In this data set, no heifers were culled before 32 DIM. This does not reflect the actual situation and is caused by the way the data were collected and handled: only heifers of which the first test-day SCC was measured between 5 and 14 DIM in the year 2000 (n = 14,766) and of which additional test-day SCC were available (n = 14,234; selection procedure outlined in De Vliegher et al., 2004a) were used. Therefore, no heifers that were culled within the first month of lactation were present in the data set. This has resulted in an underestimation of the effect of peripartum udder health problems on the hazard of culling.
The culling decision process is complex and very farmer-, herd-, and time-specific. This is especially true under a quota system (as implemented in Belgium), forcing the farmer to take culling decisions depending on the level of expected fulfillment of the quota. This study ignored some of the points in the complex decision of the farmer as the required herd and time-level information was not available. Economic calculations and implementation of new preventive measures would in addition require cow-specific values. Future studies should collect this information and combine it with data from other studies looking at the association between udder health in early lactation in heifers and udder health (both clinical and subclinical mastitis), production, and fertility in the first and following lactations.
| CONCLUSIONS |
|---|
|
|
|---|
50,000 cells/mL, suggesting that prevention rather than cure of an elevated SCCel is needed.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication August 16, 2004. Accepted for publication November 11, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. Bar, Y. T. Grohn, G. Bennett, R. N. Gonzalez, J. A. Hertl, H. F. Schulte, L. W. Tauer, F. L. Welcome, and Y. H. Schukken Effects of Repeated Episodes of Generic Clinical Mastitis on Mortality and Culling in Dairy Cows J Dairy Sci, June 1, 2008; 91(6): 2196 - 2204. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Green, A. J. Bradley, G. F. Medley, and W. J. Browne Cow, Farm, and Herd Management Factors in the Dry Period Associated with Raised Somatic Cell Counts in Early Lactation J Dairy Sci, April 1, 2008; 91(4): 1403 - 1415. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. I. Parker, C. Compton, F. M. Anniss, A. Weir, C. Heuer, and S. McDougall Subclinical and Clinical Mastitis in Heifers Following the Use of a Teat Sealant Precalving J Dairy Sci, January 1, 2007; 90(1): 207 - 218. [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 |