|
|
||||||||
1 Danish Institute of Agricultural Sciences, Department of Animal Health and Welfare, PO Box 50, DK-8830 Tjele, Denmark
2 Department of Clinical Studies, Large Animal Medicine, and
3 Department of Animal Science and Animal Health, Division of Epidemiology, The Royal Veterinary and Agricultural University, DK-1870 Frederiksberg C, Denmark
4 Landeskontrollverband Schleswig-Holstein e.V., D-24106 Kiel, Germany
Corresponding author: Ilka C. Klaas; e-mail: ilka.klaas{at}agrsci.dk.
| ABSTRACT |
|---|
|
|
|---|
Key Words: milk leakage teat shape peak milk flow
Abbreviation key: DevMilk = deviation from the population mean milk yield adjusted for stage of lactation and parity, DevPMF = deviation from the population mean peak milk flow adjusted for stage of lactation and parity, MeanResMilk = mean residual of observed expected milk yield
| INTRODUCTION |
|---|
|
|
|---|
-blocking agents (Vandeputte-Van Messom et al., 1984). High-yielding and easy milking cows are reported to be at greater risk for milk leakage (Wendt et al., 1994). Most of our understanding about causes for milk leakage is derived from experimental studies under controlled conditions using small numbers of cows. The aim of our study was to identify risk factors for milk leakage on an individual cow level and to examine the relationships between milk yield, peak milk flow, teat shape, and teat orifice condition under field conditions. Further, we wanted to estimate the variance components resulting from herd and individual cows.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Cows were examined clinically at the morning milking of d 2. Milk leakage was recorded in the milking parlor while the cows entered the milking parlor from the holding area. Milk leakage was recorded for each quarter and was defined as milk dropping or flowing from the teat. Immediately after detachment of the milking cluster, udders and teats of all cows from 9 to 100 DIM and during late lactation (>250 DIM) were examined (visual inspection and palpation). Absence or presence of edema was recorded. Classification of teat and teat-end shape was based on Rosenberger (1979). Teat orifice was classified as normal, teat canal protrusion, white smooth callosity ring, or rough callosity ring. Definitions of teat characteristics are summarized in Table 2
. Teat characteristics were recorded for each quarter. Udder size was classified for each cow by assessing the height of the ventral floor (teats not included) relative to the hocks (Table 2
).
|
Milk leakage was treated as a binary response variable in a random effect logistic regression model using the GLIMMIX macro of SAS (Littell et al., 1996). The initial hierarchical structure was: herds, cows within herd, and quarters within cow. These models resulted in serious underdispersion. Dispersion parameters for the different models were in the interval 0.6 to 0.13. The dispersion parameter allows the residual variances to deviate from their expected values. A value >1 indicates more variation than expected by chance and is referred to as overdispersion. Similarly, a value <1 indicates less than expected variation and is referred to as underdispersion (Brown and Prescott, 1999). A solution to this problem could have been to aggregate quarter observations for each cow, because only 9 of 131 cows having milk leakage had <4 affected quarters. Because the levels in the categorical variables could not be regarded as linear, an aggregation on a cow level as a cow-median or cow-mean was not reasonable. Therefore, one quarter per cow was selected randomly. This approach still permitted evaluation of quarter-level risk factors. The model building process continued with herd as a random effect.
To screen for multicollinearity, all variables were checked using Pearsons correlation coefficients (continuous variables only) or Spearmans rank correlation coefficients. We defined the occurrence of multicollinearity at coefficients >0.5. According to this definition, multicollinearity was present between peak milk flow and milk yield, milk yield and parity, milk yield and stage of lactation, and parity and teat shape.
To control for multicollinearity, milk yield and peak milk flow were corrected for parity and stage of lactation. Instead of milk yield measured during the Lacto-Corder recordings, the energy-corrected total milk yield on the official test-day was used to model lactation curves for all cows. By applying a piecewise-linear model, a lactation curve for 3 parity groups (1, 2, and >2) was fitted. The method was suggested by Enevoldsen et al. (2000) and applied by Nielsen et al. (2002). A cows mean residual derived from the piecewise-linear model characterized milk yield as low, average, or high. Mean residuals (MeanResMilk) were checked for outliers (4 observations with >4 times the standard error were excluded) and normal distribution. As a continuous variable, MeanResMilk entered the random effect logistic regression model with milk leakage as the binary outcome. A second approach to test whether volume of milk affected the risk of milk leakage was to transform the milk yield in deviation from the mean milk yield (DevMilk) within stage of lactation (see Table 3
for grouping) and parity (grouped as 1, 2, and >2) across herds.
|
The model building process started with the baseline random intercept model to evaluate the random herd effect. Risk factors were first analyzed with univariable models, each one tested for its effect on the random herd effect. The final model was generated using stepwise backwards elimination. In each step, model quality was checked by the log-likelihood ratio test and by evaluation of Akaikes information criteria (Brown and Prescott, 1999). Coefficient estimates in different models were checked for changes indicating confounding or poor data quality. Risk factors tested in the initial model were DevPMF, MeanResMilk, DevMilk, udder edema, breed, teat orifice, teat shape, teat-end shape, udder size, season (JanuaryApril, MaySeptember, OctoberDecember in 1997 and 1998), and DIM.
Continuous variables, peak milk flow, DevPMF, MeanResMilk, DevMilk, and DIM were checked for linear relationship with the outcome variable by testing the model with polynomials of different degrees. Furthermore, linearity was tested by transforming the continuous variables into categorical variables. Results indicated a linear relationship. No first-order interactions between continuous variables were found. Interactions between continuous variables and categorical variables were not significant. Models including interaction terms between teat-end shape, teat shape, and teat orifice did not converge.
The final model chosen for primiparous cows included DevPMF, MeanResMilk, teat orifice, and DIM as fixed effects and herd as a random effect. The final model for multiparous cows included DevPMF, MeanResMilk, teat orifice, teat shape, teat-end shape, quarter position of the randomly chosen udder quarter, and DIM as fixed effects, and herd as a random effect. The random logistic regression model was as follows:
![]() | ([1]) |
where P is the probability for milk leakage, ß0 = intercept, ßi = regression coefficient for the continuous risk factor i, RFi = value of risk factor i, FRF is the fixed effect of the categorical risk factor j on level k, and RHm = random herd effect for herd m.
Variances in the models were measured on a logit scale, whereas the variance on residual level was measured on a binomial scale. Under the assumption that no extrabinomial variation is permitted, level 1 (cow) variance was constrained to be equal to 1 on the binomial scale. Estimates of the proportion of variation explained by herd were computed by assuming that level 1 variance on a logit scale was
2/3 (Dohoo et al., 2001).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
When analyzing udder-related risk factors, problems arise concerning the assessment of different quarters within a cow. In our study, the response variable (milk leakage), as we measured and defined it, was more likely to affect the whole cow than individual quarters. Of cows having teats leaking milk, only 9 showed leakage in less than 4 quarters. This may have been the main reason for underdispersion occurring in a hierarchical model in which all quarters were included as a random effect. When there are uniform effect categories, the dispersion parameter will be <1, and removing observations corresponding to the uniform effects is preferable (Brown and Prescott, 1999). Consequently, inclusion of the quarter level in the hierarchical model did not provide much additional information compared with the cow-level model with herd as a random effect. Furthermore, quarter position was insignificant as a fixed effect in the cow-level model for primiparous cows, and of borderline significance in the model for multiparous cows. Different random selections of one quarter per cow were analyzed to check robustness. The overall conclusion, however, regarding significance of the evaluated risk factors was unchanged (data not shown).
The disadvantage of randomly selecting one quarter per cow is the decrease in sample size. A further decrease in sample size was necessary to control the correlation between parity and teat shape (i.e., the risk factor analysis had to be done separately for primiparous and multiparous cows). The decrease in sample size could have affected the results of the model. On the other hand, risk factors with low P values such as DevPMF or teat orifice remained unchanged when analyzed in different random quarter samples. The prevalence of acute lesions and scar tissue at the teat orifice was low, and only one quarter showed milk leakage. Consequently, we could not assess the effect of acute pathological lesions or scar tissue at the teat orifice on the quarter risk and on the cow risk for milk leakage.
As expected and reported for cows milked by an automatic milking system (Persson Waller et al., 2003), larger peak milk flow rates were associated with greater risks for milk leakage in primiparous and multiparous cows. Peak milk flow is determined by size of the teat orifice, intramammary pressure, milk yield, vacuum applied during milking, and other characteristics of the milking machine (Blake and McDaniel, 1978). Greater vacuum pressure stretches the teat orifice. The vacuum and other milk machine characteristics are herd-specific and were not included as fixed effects in the analysis. Only 10% of the variation in the data was at the herd level, indicating that differences between herds must have been small compared with individual cow effects. From a practical point of view, cows with high peak flows are preferable because of their shorter milking duration. Based on our results showing a linear increase in the risk of milk leakage with increasing DevPMF, selection of cows with the greatest peak flow rate cannot be recommended without consideration of the resulting side effects of increased risk for mastitis and hygiene problems. Leaked milk may enhance bacterial growth in the bedding material, and thus increase the risk for environmental mastitis in herd mates.
In both primiparous and the multiparous cows, characteristics of the teat orifice were associated with the risk for milk leakage. In both groups, cows with teat canal protrusions were at greater risk for milk leakage than those having normal teat orifices (Tables 5
and 6
) and all other teat orifice levels (results not shown). Diameter of the teat canal may be larger in cows with teat canal protrusions than in cows having normal teat orifices or in cows having teat-end callosity. Cows having teat canal protrusions also may have had less sphincter muscle tone before milking.
Days in milk were a significant risk factor for milk leakage in all cows with nearly identical estimated coefficients in both parity models. A linear association was detected between DIM and the risk for milk leakage, but it has to be considered that we had no information for cows between 101 and 249 DIM. An increase in DIM by 100 d corresponded to an odds ratio for milk leakage of 0.55 in primiparous cows and 0.60 in multiparous cows. This is in contrast to other observations (Persson Waller et al., 2003), in which no relationship was detected between stage of lactation and milk leakage. In their study, cows were observed at 2-h intervals in the barn, whereas in our study, cows were observed at the time of morning milking. Differences in intramammary pressure and udder filling may have decreased the risk for milk leakage during the course of lactation-intramammary baseline pressure before milking as well as the ejection pressure is greatest during the first 4 mo of lactation (Mayer et al., 1991; Bruckmaier and Hilger, 2001) and declines as lactation advances. As the intramammary pressure decreases, the degree of udder filling before milking also decreases. The smaller the degree of udder filling, the smaller the amount of milk stored in the teat cistern, which reduces the pressure on the teat sphincter before milking is initiated.
The DevMilk was not significant in the models for multiparous and primiparous cows. In multiparous cows, no evidence of a relationship between milk leakage and MeanResMilk was detected. But primiparous cows were at less risk for milk leakage at greater MeanResMilk. Depending on the intensity of milk leakage, the cisternal milk fraction of affected cows may have been reduced before milking. Consequently, the amount of milk harvested during milking may have been less. Higher-yielding primiparous cows might have had more udder storage capacity (and thus lower intramammary pressure) than low-yielding primiparous cows. Two morphological features may be important: size of the teat or gland cistern and size of the udder. In primiparous cows, the size of the cisternal fraction is smaller than that in multiparous cows. The cisternal fraction decreases during lactation, but the decrease is less in primiparous cows than in multiparous cows (Bruckmaier and Blum, 1998). In our study, primiparous cows with tiny udders tended to have greater risk for milk leakage (results not shown). Further investigation is necessary to enlighten this relationship.
Multiparous cows with short teats were at greater risk for milk leakage, whereas teat shape did not affect the risk of milk leakage in primiparous cows. This difference may have been due in part to the smaller sample size of primiparous cows. It is possible that teat shapes are not a relevant risk factor in primiparous cows because, in general, these cows had smaller and narrower teat canals than multiparous cows. In primiparous cows, teat length, and length and diameter of teat canals increased during the first months of lactation (McDonald, 1973), making it difficult to model the effect when only one observation within lactation was analyzed. In multiparous cows, short teats increased the risk for milk leakage compared with longer and thicker teats. Short teats may have had shorter teat canals, resulting in an impaired sphincter muscle function when intramammary pressure before milking is elevated. Cows having inverted or flat teat ends seem to have a greater diameter or a shorter teat canal that may be stretched wider when the intramammary pressure is elevated.
The proportion of the total variation explained in milk leakage at the herd level was similarly small in primiparous (10.8%) and multiparous cows (10.4%), relative to that at the cow level. This suggests that changes in management practices and other herd-level related factors might have small impacts on the risk for milk leakage.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication November 14, 2003. Accepted for publication August 11, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Steeneveld, H. Hogeveen, H. W. Barkema, J. van den Broek, and R. B. M. Huirne The Influence of Cow Factors on the Incidence of Clinical Mastitis in Dairy Cows J Dairy Sci, April 1, 2008; 91(4): 1391 - 1402. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Svennersten-Sjaunja and G. Pettersson Pros and cons of automatic milking in Europe J Anim Sci, March 1, 2008; 86(13_suppl): 37 - 46. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Rovai, M. T. Kollmann, and R. M. Bruckmaier Incontinentia Lactis: Physiology and Anatomy Conducive to Milk Leakage in Dairy Cows J Dairy Sci, February 1, 2007; 90(2): 682 - 690. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |