|
|
||||||||
1 Danish Institute of Agricultural Sciences, Department of Animal Health and Welfare, DK-8830 Tjele, Denmark
2 The Royal Veterinary and Agricultural University, Department of Clinical Studies, Large Animal Medicine, 1870 Frederiksberg C, Denmark
3 The Royal Veterinary and Agricultural University, Department of Animal Science and Animal Health, 1870 Frederiksberg C, Denmark
Corresponding author: I. C. Klaas; E-mail: ilka.klaas{at}agrsci.dk.
| ABSTRACT |
|---|
|
|
|---|
Key Words: udder health principal component analysis clinical investigation udder characteristics
Abbreviation key: ECM = energy-corrected milk, MLOGSCC = mean of log-transformed SCC, MRESMILK = mean residual of observed-expected milk yield, PCA = principal component analysis
| INTRODUCTION |
|---|
|
|
|---|
Consequently, the development of more practical, precise, and cheaper methods to assess udder health on farms is needed. Systematic clinical examination of udder and teat characteristics is relatively inexpensive and may provide valuable information. Collection of such data is quite feasible, especially in herds already applying a herd health management system including other clinical recordings, e.g., body condition scoring or other paraclinical data (Markusfeld et al., 1997, Enevoldsen et al., 2000). In such herds, the marginal costs related to systematic udder examinations are low.
The overall objective of the present study was to evaluate the applicability of systematic clinical udder assessments of individual cows as a supplement to SCC and udder treatment records. The specific objectives of the study were 1) to identify types of udders of cows by using biologically meaningful combinations of systematic clinical assessments of individual cows during cross-sectional examinations and 2) to estimate the relationship between udder traits and milk yield, as a production parameter, and SCC, as a marker of inflammation.
To identify the major udder types, principal component analysis (PCA) was used to reduce the number of observed variables to a limited number of artificial variables, the "principal components."
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Statistical Analysis
Principal component analysis was performed using the factor procedure of SAS (Proc factor, SAS, 1999). The number of components to retain was determined by evaluating the eigenvalue-one criterion, the scree test, the proportion of variance accounted for, and the interpretability criteria (Hatcher, 1994). Because each observed variable contributes one unit of variance to the total variance in the data set, a component displaying eigenvalues >1 accounts for a greater amount of variance than one variable alone and is worthy to be retained. The scree test displayed 4 components with the largest eigenvalues before a visible break. Further, a component was regarded as important and kept for further analysis if at least 3 variables had high loadings on that component and shared the same conceptual meaning. A factor loading greater than ± 0.35 was considered important.
A loading represents the correlation of a given variable with the underlying factor. Loadings range from 1 to +1. Principal component analysis attempts to identify a minimum number of components explaining the majority of the variation in the original data. The first component extracted in a PCA accounts for a maximal amount of total variance in the observed variables (Hatcher, 1994). Typically, the second component will be correlated with some of the observed variables that did not display strong correlations with the first component. Further, the second component is not correlated with the first component. The successive components identified explain progressively smaller portions of the total sample variance. Linear combinations of observed variables are formed. In contrast to common factor analysis, in PCA, no assumptions about the underlying common parts (communalities) and unique parts of the variables are made (Rummel, 1970).
Because ordinary PCA does not account for the effect of the individual cow, only the last observation per cow was included in the analysis (n = 707). Therefore, only clinical examinations conducted by the 2 technicians were analyzed. Principal component analysis was carried out with categorical and continuous variables (Table 2
). Udder shape, teat shape, parity, and stage of lactation were treated as dummy variables with normal udder and teat shape, second parity, and middle part of lactation as reference level. In cases of these dichotomous variables, the loadings on a factor were interpreted as the probability of its presence (positive loading) or absence (negative loading), given that the factor exists (Rummel, 1970).
Three different PCA were performed. First (PCA 1), only clinical variables were included in the analysis to simulate the situation of the practitioner on farm with no further information about the cows. Second (PCA 2), the variables parity and stage of lactation were added as dichotomous variables. Third (PCA 3), milk yield and previous SCC were added. The non-orthogonal promax rotation with a power of 3 was used to enhance the interpretability of the factors.
Milk Yield
Milk yield from each test day of all cows in the herds during 1999 that had at least 9 DIM was transformed to energy-corrected milk (ECM) using the following formula (Anonymous, 2002):
![]() | ([1]) |
When analyzing ECM from test day records to model the lactation curve, the change in slope was expected to occur around 60 d after calving (Enevoldsen et al., 2000). A standard lactation milk curve of the current lactation was estimated for each cow using a piecewise linear regression model (Proc mixed; SAS, 2000) as described by Enevoldsen et al. (2000). The peak of the lactation curve was set to d 60 after calving. It was the intercept of the model and was set to 0 (DIM 60):
![]() | ([2]) |
The second variable (DIMB60) was defined, which has the value 1 if DIM is below 60, otherwise 0. This variable models the slope from calving until 60 DIM:
![]() | ([3]) |
The model accounted for the effect of parity and stage of lactation and was
![]() | ([4]) |
where
| Yijk | = | estimated ECM
| µ | = | expected mean at 60 DIM among third parity cows
| parityi | = | fixed effect of parity (i = 1, 2, and 3+)
| DIM 60(parity)j | = | fixed effect of DIM within parity, and
| DIMB60(parity)k | = | fixed effect of before or after d 60 within parity.
|
The intercept is an estimate of the peak milk production at d 60 after calving and was 32.4 kg for cows in the third parity group (Table 3
). The estimate of DIM 60(parity) is the estimated difference in milk production between exactly 60 DIM and 305 DIM within the different parity groups. Cows in first lactation had the lowest difference of 4.1 kg ECM between DIM 60 and DIM 305, indicating the highest persistency. The estimated difference in milk yield between the day of calving and 60 DIM in second parity cows was not significant, indicating a falling line from calving until d 305 instead of a peak at 60 DIM, which is common in models of milk yield where milk is not adjusted to ECM.
|
SCC
The previous SCC was the log-transformed (log10) SCC of all test days of the examined cows from calving until the day of their clinical examination. The mean gives information about the infection status in a cow during lactation and was included as MLOGSCC in the third step, PCA 3.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
With the introduction of non-clinical variables in PCA 2 and PCA 3, the biological meaning of the components emerges more clearly. In PCA 1, the component "small udder" starts with high positive loadings (probabilities) of small udder shape and short teats together with high negative loadings of deep udder shape and degree of teat end callosity (PCA 1). The addition of parity in PCA 2 and SCC in PCA 3 explains the "small udder" further as having a higher probability among first parity cows in contrast to cows in second (reference level) or higher parity and in contrast to cows having normal (reference level) or deep udder shape. Furthermore, cows with high loadings on that component had a lower MLOGSCC.
From PCA 1 to PCA 2, the order of the components changed. When only clinical variables were analyzed, the component "distressed udder" explained most of the variation in data. When the variables parity and stage of lactation were added, the component "distressed udder" moved to the second position, and "small udder" explained most of the variation of the data and remained the first extracted component in PCA 3. From PCA 2 to PCA 3 the components "mastitis udder" and "soiled udder" changed their order and meaning. The oblique rotation revealed a low inter-correlation between the components in general (Table 8
). The highest inter-component correlation (0.17) was found between component 1 "small udder" and component 4 "soiled udder". That is, in general, small udders were slightly less soiled.
|
| DISCUSSION |
|---|
|
|
|---|
The small udder was related to lower SCC. Deeper udders indicate a decrease in teat end to floor distance, which is correlated with higher SCC (Rogers et al., 1991; Faye et al., 1998) and is even related to a higher risk of subclinical (Ronningen and Reitan, 1990) and clinical (Slettbakk et al., 1995) mastitis. In general, cows in higher parities have higher SCC (Laevens et al., 1997; De Haas et al., 2002). However, in cows that are uninfected by mastitis pathogens, SCC is uncorrelated with parity (Laevens et al., 1997). In conclusion, cows with a high loading on the component "small udder" may be less susceptible to mastitis than cows with normal or deep udder shape.
Component 2, The Distressed Udder
In comparison with normal udder shape, long udders and backward-bulging udders were related to harder tissue and an impaired teat surface with rough teat skin and wounds on warts and teats. Rough teat skin (Burmeister et al., 1998) and wounds on skin (Larsen et al., 2000) are a reservoir for S. aureus. Even superficial lesions of the teat skin increase the risk of subclinical and clinical mastitis (Agger and Willeberg, 1986). Warts caused by viruses can result in local damage of the natural defense mechanisms of the teat and predispose them for bacterial mastitis (Wellenberg et al., 2002).
Hardness of udder tissue characterizes two different morphological features: the udder parenchyma with aggregations of alveoli (secretory cells) partitioned off into lobules by thin connective tissue septa growing into stronger connective tissue, the second morphologic feature, which builds up the udder attachment with the ligaments. The percentage of secretory tissue differs among individual dairy cows and genotypes and is influenced by parity and stage of lactation (Michel, 1994). In freshly calved first parity cows, the percentage of fat tissue is high, and the percentage of secretory tissue is low. During the first lactation, the percentage of secretory tissue increases and replaces the fat tissue, accompanied by increasing udder hardness. In subsequent parities, number and size of alveoli increase further until the fifth or sixth parity. The secretory tissue decreases again during the declining phase of lactation and during involution of the secretory tissue. In contrast to information from literature, stage of lactation and parity did not contribute to component 2. One explanation could be an interaction between connective and secretory tissue: udders appearing fleshier after milking have a heavier stroma or connective tissue than udders with a large amount of secretory tissue (Smith, 1959). Furthermore, the hardness of the udder was measured from behind with both hands palpating the udder close to the ventral abdomen where ligaments and connective tissue were strongest and percentage of secretory tissue was lowest. Cows with a strong ligament and udder capsule may also have had a stronger connective tissue in the parenchyma, leading to a higher degree in hardness than cows with a weaker ligament structure (Akers, 2000). However, previous SCC showed no relationship with high component 2 values on the "distressed udder," but cows with such "distressed udder" may be at a higher risk for subsequent clinical mastitis because of impaired teat surface.
Component 3, The Mastitis Udder
Because the incidence of truly acute clinical mastitis is of short duration (Smith and Hogan, 1995), we expected to find a low frequency of acute cases in our study. In addition to signs of acute clinical mastitis, knotty tissue and asymmetry of quarters significantly contributed to the "mastitis udder," component 3.
This component appeared to be a mixture of acute and chronic symptoms of clinical mastitis. Hardness of mastitis quarter and diagnosed clinical mastitis are acute symptoms, whereas asymmetry and knotty tissue are chronic alterations of the udder. Chronic mastitis can lead to atrophy of the affected quarter, resulting in asymmetry as a post-inflammatory state. Mastitis can also develop progressively with increases in fibrous tissue so that the size of the affected quarter increases (Fuchs, 1994). Further, the farmer may have decided to dry off a mastitis quarter, resulting in involution of the affected quarter. Because these chronic alterations were longer-lasting conditions, the observed frequency in the presented investigation is much higher than the frequency of acute mastitis. The clinical variables with high loadings point toward the same resulta clinical (visible) mastitisand can be used as additional classification of clinical mastitis. Especially in situations where under-reporting is suspected, this component can improve the quality of data because it decreases the number of false-negative non-mastitis cows. A problem when using this component to classify udder health can be that it only tells that there was previous mastitis but tells nothing about the duration or the onset of disease. Conversely, udders once affected by clinical mastitis are more susceptible to recurrent cases of mastitis with acute symptoms (Hogan et al., 1989). When there is no further information about first cases of clinical mastitis available, a distinction between acute clinical mastitis and chronic mastitis with acute episodes is difficult, especially in older cows.
As expected, the variables describing clinical signs of mastitis were related to a higher SCC and a lower milk yield in comparison with other cows in the same parity. Chronic alterations may have reduced the milk production in the affected quarter. Furthermore, the reduction in milk yield after clinical mastitis seemed to persist during the whole lactation as reported in the literature (Morris, 1973; Bartlett et al., 1991). A systematic screening for cows with high loadings on the component may be beneficial in the assessment of udder health status and in the development of a treatment strategy in the herd.
Component 4, The Soiled Udder
The soiled udder was related to early lactation and is consistent with Ward et al. (2002), who found higher counts of Eschericia coli and Staphylococcus uberis in beds of early lactation cows and a higher degree of soiling in early lactation cows. No other clinical or production variable contributed to the "soiled udder" when stage of lactation was added to the analysis. This result is in contrast to the results of a factor analysis performed by Bruun (2003), where days after calving showed no relation with degree of soiling. Contamination via manure reflects hygiene in the barn (Cook, 2002) and because degree of soiling was measured immediately after milking, hygiene practices before the milking process had an influence on the degree of soiling. These management practices affected the cows independently from their udder and teat morphology, milk production, or previous SCC.
General Discussion
Differences in the clinical examinations performed by 2 different persons might have induced a bias. Houe et al. (2002) have shown that the agreement between clinicians on pathological conditions is high, but there was a higher variation when classifying non-pathological conditions such as udder shape or hardness of udder parenchyma. Because our data recording sheets supplied the 2 clinicians with detailed descriptions and figures on specific conditions, we expect the variation between clinicians to be low. Further, the last observation of each cow was chosen to minimize bias caused by a "learning effect" of the examiners during the project.
When developing and interpreting PCA, the researcher should find biologically meaningful descriptions for the extracted components. The order of the components changed with the introduction of non-clinical variables in PCA 2 and PCA 3. Components "soiled udder" and "mastitis udder " reversed in PCA 2 and PCA 3. But, the amount of variation explained by each component was similar; a change in the order of the components was not an interpretational problem. The clinical variables with the highest loadings on a component kept their high loading and their position within a component throughout the model-building process. The stability of these variables is evidence for a meaningful pattern that was not produced by chance. When non-clinical variables were included in the analysis, the focus changed within the component "small udder," and age became the most important variable related to the size of the udder. In contrast, the component "distressed udder" consisted only of clinical variables that were not related to age, stage of lactation, milk production, or previous SCC. Therefore the "distressed udder" remains invariant from PCA 1 to PCA 3. The factor "mastitis udder" is a combination of clinical variables, previous SCC, and milk yield. However, each part of the combination adds up to the whole component, indicating that each part is of the same importance in its contribution to the "mastitis udder." The fourth component, "soiled udder," is determined only by those clinical variables describing soiling of teats and udder and stage of lactation. The presence of a low inter-component correlation with the "small udder" indicates that there is only a weak relationship between udder size and soiling.
Factor analysis was used in 2 other studies on clinical udder examinations in Denmark. Houe et al. (2002) extracted 3 factors with common factor analysis including clinical variables only: the mastitis udder, the long udder, and the high udder. Variables loading high on the mastitis udder were asymmetric front or hind quarters and dry quarters. Similar to our results, but by using the mastitis factor as a variable in a general linear model, Houe et al. (2002) showed that an increase in factor score was correlated with a decrease in lactation milk yield (305-d energy-corrected milk). As in our study, Bruun (2003) extracted 4 factors, but used common factor analysis with maximized variance. The author labelled the first factor "chronic changes," which was similar to our component "mastitis udder." In contrast to our results, no acute signs of mastitis contributed to the factor and no relationship with milk yield were found. The researcher used the milk yield recorded immediately before clinical examination, whereas the present study worked with a milk model to classify cows as high, average, or low yielders within their parity.
Comparable with Bruun et al. (2003), the total amount of variance explained was around 27%, which indicates that several variables have high degrees of uniqueness. That is, they provide important information on their own. Further, the data structure may not have been optimal; Bruun (2003) rescaled ordinal variables with the method of optimal scaling to improve linearity. When rescaled variables are analyzed, the interpretation is more difficult with regard to the original categories. In the present study, dummy variables were used, e.g., for parity and stage of lactation, when non-linearity was suspected. Regardless of this somewhat different approach, extracted factors or components were similar in both studies when the correlation to the factor (component) was strong. In addition to a different method within factor analysis and treatment of variables, differences may also be caused by the choice of variables that were included in the analysis. Introduction of new variables with important contribution to a factor (component) may shift the focus or enhance the interpretability of a factor (component) as seen in the component "small udder." Variables related to factor "udder shape" (Bruun, 2003) and to the component "small udder" of the present study were identical, which is further evidence that the PCA produced a biologically meaningful pattern and contributes to the understanding of udder health.
Welfare and health assessment protocols, including animal-based measurements, become an issue of increasing importance in health management and quality programs in several countries (Noordhuizen and Metz, 2003; Whay et al., 2003). Increasing herd sizes increase the demand for systematic, effective, and economic tools in udder health management. Inclusion of animal-based measurements could be beneficial in advisory services to understand the udder health situation on the farm and to specify udder health risks. The extracted udder types can be implemented in decision support systems on udder health in dairy farms. The presence of cows with "distressed udder" indicates impaired teat and udder condition and could be used to evaluate the effectiveness of hygienic measures or the necessity of measures in the milking parlor. The "soiled udder" could be used as an indicator of hygiene in the barn. Classifying cows according to the "small udder" can support decision making when selecting cows for mastitis prevention on the herd level. The prevalence of cows with a "mastitis udder" provides information about the number of clinically affected cows and the severity of clinical signs, which is of prognostic value. When results from clinical examination are compared with farm records, management practices, such as blinding of mastitis quarters, will be revealed, and the quality of farm records can be evaluated. Further investigations are necessary to test the applicability and the effect of animal-based measurements in udder health monitoring.
| CONCLUSIONS |
|---|
|
|
|---|
Received for publication August 8, 2003. Accepted for publication December 19, 2003.
| 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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |