|
|
||||||||
1 Institute for Animal Science and Health, ID-Lelystad, 8200 AB Lelystad, The Netherlands
2 University of Prince Edward Island, Atlantic Veterinary College, Department of Health Management, Charlottetown, PEI, Canada C1A 4P3
3 Cornell University, College of Veterinary Medicine, Department of Population Medicine and Diagnostic Sciences, Ithaca, NY 14853
Corresponding author: Y. de Haas; e-mail: Y.deHaas{at}id.wag-ur.nl.
| ABSTRACT |
|---|
|
|
|---|
Key Words: somatic cell count pattern lactation curve clinical mastitis
Abbreviation key: CM = clinical mastitis, CSCC = corrected SCC, HF = Holstein-Friesian, OR = odds ratio, SCC200 = lactation average of SCC of more or less than 200,000 cells/ml, SCM = subclinical mastitis
| INTRODUCTION |
|---|
|
|
|---|
Using test-day records of SCC instead of lactation average of SCC increases the possibilities of identifying deviations from the typical curve of SCC during lactation. These deviations might characterize mastitis-causing pathogens. If so, analyzing test-day records of SCC can be more useful in attempting to decrease the prevalence of subclinical mastitis (SCM) and the incidence of CM than the lactation average of SCC. An effective use of SCC test-day records might be achieved by defining patterns of peaks in SCC during the lactation. If these patterns of peaks in SCC provide information on 1) the pathogen distribution on a farm, 2) the success of therapy, and 3) the status of immunity of the cows, these patterns might be better tools to use in mastitis control programs or for genetic selection programs. Management can then be directed specifically on lowering the incidence of pathogen-specific CM, or shortening the duration of infection. Furthermore, decreased genetic susceptibility to the full scope of mastitis-causing pathogens might be accomplished, which cannot completely be achieved by genetic selection on lower lactation average SCC (Nash et al., 2000; De Haas et al., 2002b). Eventually, this all helps to limit the losses due to CM.
Before the patterns of peaks in SCC can be used for health management on farms or for genetics, we first have to find out if patterns of peaks in SCC do differ from the lactation average SCC. Therefore, the overall objectives of this study were: 1) to define several patterns of peaks in SCC and then, 2) to test their sensitivity and specificity, and, 3) to estimate their associations with occurrence of pathogen-specific cases of CM, also in comparison with the lactation average of SCC.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Sampling
Epidemiological studies often depend on farmer-observed disease incidence, introducing a possible bias caused by variation between farmers in diagnostic criteria. Farmer-diagnosed CM was analyzed in this study. Lam et al. (1993) concluded that the diagnostic capability of farmers does not have a negative influence on the validity of farmer-diagnosed CM. Farmer-diagnosed CM, therefore, seems to be a useful tool in epidemiological studies on CM.
Selection of herds and aseptic sampling procedures have been described previously (Barkema et al., 1998). Data collection of milk samples of quarters with CM depended heavily on the willingness of farmers. Therefore, farmers were continuously encouraged, as described by Barkema et al. (1998). During the study period, farmers took milk samples from all quarters that, in their opinion, had clinical signs of mastitis. Samples were stored in a freezer at the farm (at approximately -20°C) and were collected for bacteriological examination at intervals of 6 to 8 wk. Bacteriological culturing of milk samples was performed according to the standards of the National Mastitis Council (Harmon et al., 1990). Briefly, 0.01 mL was cultured and, for each culture, the number of colony-forming units of each of the bacterial species was counted. Collected data contained information on national cow identification number, date of mastitis occurrence, infected quarter, and result of bacteriological culturing of milk sample. Seven groups of pathogens were defined based on their incidence in the data: Staphylococcus aureus, CNS, Escherichia coli, Streptococcus dysgalactiae, Streptococcus uberis, streptococci other than Strep. dysgalactiae and Strep. uberis, and culture-negative samples. Pathogen-specific CM was determined on lactation level. Pathogen-specific cases of CM were scored as categorical traits, so when all cases of CM in a lactation were caused by one of the above-mentioned seven pathogens, a 1 was registered for this particular pathogen, otherwise it was scored as 0. The group of remaining cases of CM consisted of all lactations in which 1) cases of CM were associated with any other pathogen or mixed cultures or 2) consecutive cases of CM were not caused by the same pathogen. In total, 3393 lactations were recorded with cases of CM before 450 DIM (Table 1
).
|
Definition of Patterns
To determine the patterns of peaks in SCC, SCC was expressed relative to SCC in lactations without CM and SCM (De Haas et al., 2002a), by first estimating SCC on each DIM of these lactations with a spline function (Gilmour et al., 2002). Usually, a spline function is used for smoothing data points, and the function allows maximum flexibility and assumes no prescribed curvature (Gilmour et al., 2002). Secondly, these estimates were subtracted from the SCC test-day records for each DIM value. The procedure was carried out separately for heifers and multiparous cows, and this calculated value for SCC is referred to as the corrected SCC (CSCC). Because De Haas et al. (2002a) only used SCC up to 308 DIM, SCC was extrapolated from 309 to 450 DIM by adding the average daily increase in SCC from 100 to 300 DIM to the calculated SCC on the previous day. Here, the assumption was made that milk yield showed a linear association with DIM at the end of the lactation (>309 DIM), so SCC could be linearly extrapolated with DIM.
Five SCC patterns were defined and distinctions were made between lactations with short or longer periods of increased CSCC, and also between lactations with and without recovery within three or five test-day recordings. A recovery is classified by a recording of a low CSCC after a rise in CSCC. Pattern 1 describes a quick rise in CSCC followed by a quick decrease in CSCC and is referred to as "quick recovery pattern." Three consecutive test-day records of CSCC had to be low, high, and low again (Figure 1
). Pattern 2, referred to as "slow recovery pattern," described a slower increase in CSCC, but still with recovery, i.e., low, higher, high, lower, and low again (Figure 1
). Pattern 3 had no restrictions on CSCC on the second and fourth test-day records, but the first one had to be low, the third one had to be high, and the fifth one had to be low again. There is overlap between patterns 2 and 3, because of equal restrictions for the first, third, and fifth test day. Pattern 4 captured a longer increased SCC; i.e., low CSCC followed by four high CSCC, so no recovery took place within four test-day recordings. Pattern 4 is referred to as "no recovery pattern" (Figure 1
). To compare these patterns with the lactation average of SCC, the fifth trait was a binary parameter that indicated whether the lactation average of SCC was higher or lower than 200,000 cells/mL (scored as 1 or 0, respectively) (SCC200).
|
The first criterion of all patterns of peaks in SCC is a low CSCC. However, to be able to identify cows with high CSCC directly after calving, the criteria for the definitions were set differently in the beginning of the lactation. At the end of the lactation, the criteria for the definitions were set differently as well, because cows might have been culled before they were recorded on the required number of test days included in the definition of the patterns of peaks in SCC. Therefore, in the beginning of the lactation, only the second half of the patterns of peaks in SCC was used, and at the end of the lactation the first half of the patterns of peaks in SCC was used. With three consecutive missing recordings on SCC, it was assumed that the cow was either culled or dried off.
Patterns of peaks in SCC were scored as categorical traits, so when a pattern was present in a lactation in 3 or 5 consecutive test-day recordings, it was registered as 1, and otherwise as 0. More than one pattern could be present in one lactation (Table 1
).
Odds Ratio, Sensitivity, Specificity, and Predictive Values
Odds ratios (OR), sensitivities, and specificities were used to describe the associations between presence of patterns and incidence of pathogen-specific CM in lactations (Table 1
). Odds ratio are measures of associations and indicate the change in risk. The sensitivity is the probability that a diseased animal (here defined as a lactation with occurrence of a pathogen-specific case of CM) indeed will be classified as diseased using the test (here defined as a lactation with presence of any of the defined patterns), and the specificity is the probability that a nondiseased animal will be classified as nondiseased with the test (Noordhuizen et al., 1997). Specificity in this study is a relative specificity, since our golden standard is not without error because of variation between farmers in diagnostic criteria (Noordhuizen et al., 1997). The positive predictive value of a test is the probability that in a lactation with a given SCC pattern a case of CM with a given pathogen is diagnosed by the farmer. The negative predictive value is the probability that a case of CM with a given pathogen is not diagnosed in a lactation without a given SCC pattern. Contingency tables are fully filled, indicating that there are also lactations scored with a case of CM but without the presence of any of the patterns and vice versa, i.e., lactations without a case of CM, but with presence of a pattern.
The probabilities for occurrence of pathogen-specific CM in lactations with the presence of a SCC pattern were calculated with single-trait regressions and so were the probabilities for presence of a SCC pattern in lactations with occurrence of pathogen-specific CM. Both single-trait regressions were applied on two datasets, one including all available lactations (n = 24,012) (dataset 1), and one including only the lactations with CM (n = 3393) (dataset 2). These datasets distinguish between probabilities of classifying a diseased cow out of all cows (dataset 1), and to identify the involved pathogen in cows with CM (dataset 2).
The effect of occurrence of pathogen-specific CM on the presence of one of the patterns was also investigated with a multivariate regression. Only one pathogen can occur in a lactation, and, therefore, pathogen-specific cases of CM are fully independent. Odds ratios, therefore, indicate probabilities on occurrence of specific pathogens versus no CM (dataset 1), or versus the group of remaining records (dataset 2).
Multivariate regression was also applied to investigate the effect of presence of any of the patterns in a lactation on occurrence of pathogen-specific CM. More than one pattern can be present in a lactation, so patterns are not independent. Odds ratios indicate the probability of the presence of a certain pattern versus no presence of patterns 1, 2, 3, or 4 or SCC200, independent of the analyzed dataset. Fit of the models was evaluated based on the estimated deviance of the model.
Statistical Analyses
Statistical analyses were carried out to test whether OR differed from 1, using logistic regression in SAS [PROC GENMOD; (SAS/STAT, 2001)]. For investigating the effect of occurrence of pathogen-specific CM on presence of patterns, the analyses were done separately for each pattern (1, 2, 3, and 4), and for SCC200, with this general model:
![]() | ([1]) |
where the outcome variable is a binary indicator variable (0/1) for any of the patterns, with j = 0...5 (pattern 1 to 4, or SCC200),
= intercept, and e = residual random error. The modeled probability was presence of a pattern and included regression coefficients were pathogen-specific CM (ß18), parity (ß9), and the fraction of HF genes (ß10). Eight classes for pathogen-specific CM were defined, i.e., seven classes for defined pathogens, and one group of remaining records. Pathogen-specific CM was recorded as a binary trait (0/1).
For investigating the effect of presence of patterns on occurrence of pathogen-specific CM, the analyses were done separately for each pathogen, with this general model:
![]() | ([2]) |
where the outcome variable is a binary indicator variable (0/1) for pathogen-specific cases of CM, with j = 0...8,
= intercept, and e = residual random error. The modeled probability was occurrence of pathogen-specific CM, and regression coefficients were estimated for the SCC patterns (ß14), SCC200 (ß5), parity (ß6), and the fraction of HF genes (ß7).
Parity had two classes, distinguishing between heifers and multiparous cows, and five classes of fractions of HF genes were defined (0, 1/4, ..., 4/4). Herd was included as a repeated variable (274 classes).
| RESULTS |
|---|
|
|
|---|
Long increased SCC was more frequent in cows with 26 to 75% HF, compared with cows with 76 to 100% HF. Furthermore, more lactations of cows with
75% HF had averages of SCC of > 200,000 cells/mL, than of cows with 76 to 100% HF. Similar results were found in dataset 2; i.e., lactations of cows with 26 to 50% HF showed more often high averages of SCC (>200,000 cells/mL), than lactations of cows with 76 to 100% HF.
Single-Trait Analyses for Associations Between Pattern and Pathogen
All calculated OR in dataset 1 are above 1, indicating that a diseased cow can be classified out of all cows by presence of SCC patterns (Table 2
). Odds ratios for SCC200 were highest for cases of CM associated with Staph. aureus, Strep. dysgalactiae, Strep. uberis, and the group of remaining records. Cases of CM associated with CNS or streptococci other than Strep. dysgalactiae and Strep. uberis showed highest OR for the "no recovery pattern." Chances on the presence of the "quick recovery pattern" were high in lactations with clinical E. coli mastitis compared with the presence of another SCC pattern. No single unambiguous SCC pattern was present for culture-negative samples.
|
The negative predictive values were high, especially when analyzing all lactations, compared with analyzing only lactations with CM (Table 2
). This indicated that in at least 95% of all lactations, absence of pathogen-specific CM was correctly predicted by absence of any of the defined patterns of peaks in SCC. However, probabilities that the presence of a pattern correctly predicted occurrence of pathogen-specific CM were low (
17%) when analyzing dataset 1. When analyzing dataset 2, positive predictive values were higher than with dataset 1. In lactations with CM and with presence of the "no recovery pattern" the predictive value for Staph. aureus is 19%, whereas this is only 9% for E. coli. On the other hand, the "quick recovery pattern" predicted 22% of the cases of E. coli mastitis in lactations with CM, and only 8% of the cases of Staph. aureus mastitis.
Multivariate Analyses for Effect of Pathogens on Pattern
Slightly higher OR were calculated with multivariate analyses (Tables 3
and 4
), compared with single-trait analyses (Table 2
), especially when analyzing dataset 1, and only to some extent when analyzing dataset 2.
|
|
In lactations with CM, the odds of "quick recovery pattern" (versus none of the SCC patterns being present) were roughly 50% when Staph. aureus or CNS were isolated (OR = 0.60 and 0.47, respectively) (Table 4
). Probabilities of observing cases of CM associated with E. coli, Strep. dysgalactiae, or culture-negative samples were low compared with lactations with presence of a "no recovery pattern" and a case of CM associated with the group of remaining records.
High OR were estimated for SCC200 when analyzing dataset 1 (Table 3
), indicating higher chances of finding high lactation averages of SCC in lactations with pathogen-specific CM than in lactations without CM. Among the lactations with CM, Staph. aureus was the only pathogen that was associated with an increased odds of SCC200, compared with the group of remaining records (Table 4
). Most other pathogens were associated with a decreased odds of SCC200, compared with the group of remaining records, although Strep. dysgalactiae and Strep. uberis differed not significantly from 1.
Changing the threshold values had no severe impact on the associations between occurrence of pathogen-specific CM and presence of patterns. Interpretations of OR were similar, independent of the threshold values.
Multivariate Analyses for Effect of Patterns on Pathogen
With multivariate analyses, lower OR were calculated (Table 5
) compared with single-trait analyses (Table 2
), especially when analyzing dataset 1. Different combinations of patterns included in the model resulted in different estimates for OR, and possibilities to elucidate specific mastitis-causing pathogens differed as well. For Staph. aureus, a better fit was found for the model including all patterns versus separate models with each pattern. This indicated that for Staph. aureus a model can be defined including a combination of the patterns, resulting in an optimal fit of the model (results not shown). For other pathogens, including more patterns in the model, did not improve the fit of the model, so a model with only one pattern is preferred (results not shown).
|
| DISCUSSION |
|---|
|
|
|---|
4, respectively. This might explain the higher number of peaks for older cows, compared with cows of first, second, or third parity.
Patterns of peaks in SCC are more often present in cows with relatively low percentages of HF. This might be (partly) explained by the higher incidences of cases of CM in cows with
75% HF, compared with cows with >76% HF. This distribution has been presented before by De Haas et al. (2002a).
Biological Interpretation of Patterns
Intramammary infections with Staph. aureus are often characterized by a long duration and high SCC (Sears et al., 1990; Daley et al., 1991), although a wave-like pattern of increased SCC may be present (Daley et al., 1991). Cases of clinical Staph. aureus mastitis occurred more frequently in later lactation (Vaarst and Enevoldsen, 1997) and may therefore be more easily detected using the patterns. As expected, Staph. aureus was more often associated with the "no recovery pattern," and less often with the "quick recovery pattern."
Phenotypic associations in the current study also showed that clinical E. coli mastitis is associated with short peaks in SCC. Escherichia coli infections are typically acute cases (Vaarst and Enevoldsen, 1997), which agrees with the strong association found between the "quick recovery pattern" SCC and clinical E. coli mastitis in the current study. These results are confirmed by studies that used either induced cases or naturally occurring cases of CM. The effect of naturally occurring cases of clinical E. coli mastitis on the curve of SCC during lactation has been described using the same dataset (De Haas et al., 2002a). Before a case of clinical E. coli mastitis, SCC was close to the SCC level of lactations without CM and SCM, and after the case SCC came down to a level that was only slightly higher than the preinfection level. From experiments inducing E. coli, it was shown that SCC peaks 2 d after inoculation, and preinfection SCC values returned within 3 to 4 wk after challenge (Lohuis et al., 1990; Erskine et al., 1992; Kremer et al., 1993; Pyörälä et al., 1994). In contrast, within 24 h after inoculation with Staph. aureus, SCC increased and remained high for at least 48 d (Shoshani et al., 2000). Culture-negative CM showed strong similarities with clinical coliform mastitis (Vaarst and Enevoldsen, 1997), and associations between patterns and culture-negative CM in the current study were indeed similar to those with E. coli mastitis.
Based on the results in the current study, clinical Strep. dysgalactiae mastitis was not strongly associated with the presence of one pattern in particular, because OR were around 1, which agrees with the specific characteristics of Strep. dysgalactiae reported by Vaarst and Enevoldsen (1997). Leigh (1999) concluded that clinical Strep. uberis mastitis can be associated with longer increased SCC, which was also established in the current study when analyzing dataset 1.
When interpreting the patterns of peaks in SCC, it should be taken into account that the pattern after the rise in SCC is affected by the type of therapy of the cow. However, because all cows with CM were treated, this effect is present for each pattern. In practice, cows are treated as well; hence, the data used in this study correspond with data that are available to farmers. Responses to the therapy might differ per cow, and with the data analyzed in this study it is impossible to specify whether cure occurs spontaneous or as a result of therapy, but that is also unknown in practical situations. To diagnose individual cows by using the SCC patterns might not be the best method, because of the time needed to complete a pattern. For individual cases, bacteriology of milk samples will provide quicker and more accurate results. However, when analyzing these SCC patterns at herd level, a predominant pathogen picture might occur. This can help with introducing pathogen-specific mastitis control programs, consisting of adequate guidelines to control IMI with both contagious and environmental pathogens.
Comparison of Patterns 1 to 4 with SCC200
Although it is not exactly the same trait as the lactation average of SCC, SCC200 was defined to be able to make the comparison with average lactation values for SCC, which is currently used to monitor udder health. The major difference is that lactation average SCC gives as much weight to a reduction of 150,000 to 50,000 cells/mL, as from 900,000 to 800,000 cells/mL, whereas SCC200 only distinguishes between high versus low lactation averages of SCC. Therefore, associations between lactation average SCC and pathogen-specific CM might differ from those found for SCC200. Heuven (1987) compared several SCC variables as predictors of bacteriological status and concluded that cows with high deviations from the typical curve of SCC during lactation were more likely to be treated for CM. Still, lactation average SCC was a more accurate predictor of pathogens than a single SCC measure on bacteriological sample day (Heuven, 1987). Based on comparison of the positive predictive values, occurrence of pathogen-specific CM was equally or even more correctly predicted by presence of one of the SCC patterns than by presence of SCC200. Lactations with presence of the "no recovery pattern" had more chance on occurrence of clinical Staph. aureus mastitis, compared with the base probability of Staph. aureus (2%) (Table 1
). Higher probabilities than the base probability were found for E. coli when the "quick recovery pattern" was present in a lactation. Lactations with CM and with presence of the "no recovery pattern" had high probabilities on clinical Staph. aureus or CNS mastitis, but low probabilities on E. coli or culture-negative CM. The opposite was true for lactations with presence of the "quick recovery pattern". The results support the hypothesis that a different definition of SCC can increase the predictive value of the bacteriological status of a cow, but a combination of several definitions of SCC might be necessary to capture the full scope of mastitis-causing pathogens. Patterns of peaks in SCC may be useful as basic tools for health management advice, as they do provide additional information on, for example, the incidence of (spontaneous or therapeutic) cure, in comparison with the information provided by lactation average SCC alone. Further analyses will be done on the distributions of the SCC patterns on farms and whether this is informative for pathogen-distribution on that farm. Genetic analyses will be carried out as well to establish whether decreased genetic susceptibility to mastitis-causing pathogens can be achieved by genetic selection on patterns of peaks in SCC.
Further Optimization of Trait Definition
Several improvements in the current analyses can be made that are expected to increase the detection of the SCC patterns or to increase the chance of finding an association between SCC patterns and CM. For example, noise in SCC test-day recordings might be removed from the data by correcting for fixed effects, such as herd-test-date, in the statistical model. Chances of finding a peak in SCC were also determined by the frequency of the test-day recordings. Recording an increased SCC on the test day depended on both the day of occurrence of cases of CM in relation to test-day recordings and the duration of increased SCC as a result of pathogen-specific CM. Using this additional information will mainly help to detect acute cases of CM, but not the chronic cases because SCC remains high for some time with a chronic case of CM. Therefore, chances of detecting such a high SCC are higher compared with detecting an increased SCC as a result of an acute case of CM since SCC rises quickly and is low before and after an acute case of CM.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication July 2, 2002. Accepted for publication May 8, 2003.
| REFERENCES |
|---|
|
|
|---|
100.000/ml). Milchwissenschaft 57:183187.
This article has been cited by other articles:
![]() |
O. Reksen, L. Solverod, and O. Osteras Relationships Between Milk Culture Results and Composite Milk Somatic Cell Counts in Norwegian Dairy Cattle J Dairy Sci, August 1, 2008; 91(8): 3102 - 3113. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. P. Sorensen, B. Guldbrandtsen, J. R. Thomasen, and M. S. Lund Pathogen-Specific Effects of Quantitative Trait Loci Affecting Clinical Mastitis and Somatic Cell Count in Danish Holstein Cattle J Dairy Sci, June 1, 2008; 91(6): 2493 - 2500. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. de Haas, W. Ouweltjes, J. t. Napel, J. J. Windig, and G. de Jong Alternative Somatic Cell Count Traits as Mastitis Indicators for Genetic Selection J Dairy Sci, June 1, 2008; 91(6): 2501 - 2511. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
O. Reksen, L. Solverod, and O. Osteras Relationships Between Milk Culture Results and Milk Yield in Norwegian Dairy Cattle J Dairy Sci, October 1, 2007; 90(10): 4670 - 4678. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Lievaart, H. W. Barkema, W. D. J. Kremer, J. van den Broek, J. H. M. Verheijden, and J. A. P. Heesterbeek Effect of Herd Characteristics, Management Practices, and Season on Different Categories of the Herd Somatic Cell Count J Dairy Sci, September 1, 2007; 90(9): 4137 - 4144. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. G. M. O. Riekerink, H. W. Barkema, and H. Stryhn The Effect of Season on Somatic Cell Count and the Incidence of Clinical Mastitis J Dairy Sci, April 1, 2007; 90(4): 1704 - 1715. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ostergaard, M. G. G. Chagunda, N. C. Friggens, T. W. Bennedsgaard, and I. C. Klaas A Stochastic Model Simulating Pathogen-Specific Mastitis Control in a Dairy Herd J Dairy Sci, December 1, 2005; 88(12): 4243 - 4257. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. L. Calus, J. J. Windig, and R. F. Veerkamp Associations Among Descriptors of Herd Management and Phenotypic and Genetic Levels of Health and Fertility J Dairy Sci, June 1, 2005; 88(6): 2178 - 2189. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |