|
|
||||||||
1 Animal Epidemiology Research Unit, INRA, 63122 Saint Genès Champanelle, France
2 France Contrôle Laitier, 167 rue du Chevaleret, 75013 Paris, France
Corresponding author: Jacques Barnouin; e-mail: barnouin{at}clermont.inra.fr.
| ABSTRACT |
|---|
|
|
|---|
Key Words: clinical mastitis somatic cell score dairy management
Abbreviation key: 36-mo SCS = mean SCS for the 36 mo preceding the beginning of the survey, CM = clinical mastitis, IRCM = incidence rate of clinical mastitis, ZMP = Zero Mastitis Objective program.
| INTRODUCTION |
|---|
|
|
|---|
The very few studies conducted to highlight management factors associated with IRCM in low SCC herds (Schukken et al., 1990; Peeler et al., 2000) indicated a large IRCM variation in this ideal epidemiological situation. A simultaneous evaluation of the risk factors associated with a low SCC and those associated with a lower or higher IRCM is necessary (Barkema et al., 1998), particularly because SCC levels, contrary to IRCM, have fallen over the last 20 yr (Ely et al., 2003; Schukken et al., 2003), whereas the importance of environmental vs. contagious pathogens has increased (Peeler et al., 2003). Nevertheless, if low SCC herds are considered to have higher levels of environmental mastitis (Peeler et al., 2000), clinical mastitis (CM) caused by Staphylococcus aureus has been observed more frequently in herds with a bulk milk cell count lower than 150,000 cells/mL (Elbers et al., 1998).
The main herd risk factors associated with IRCM in low SCC herds (Schukken et al., 1990; Peeler et al., 2000) were: 1) management variables that increased the exposure to environmental pathogens (straw yard housing, poor free stall cleanliness, drinking water from sources other than public water, leaking milk outside the parlor, cleaning the calving area less than once a month, scraping the gathering yard less than twice a day); 2) variables associated with bacterial challenge or host resistance (high percentage of cows leaking milk, postmilking teat disinfection, high frequency of free stall disinfection); 3) practices that could be linked to mastitis problems in the herd (wearing of rubber gloves by the dairyman during milking, changing teat liners more frequently); 4) increased detection and reporting of cases of mastitis (stripping foremilk before attaching the clusters); 5) general management practices (long dry period, high replacement rate); 6) genetic and udder conformation dependent traits (breed, milk yield level); and 7) dietary components (sugar beet pulp in the ration). Nevertheless, it is not always easy to clarify whether a factor is cause or effect: for example, farmers may change liners more frequently because of mastitis, or frequent liner changes may potentially damage the teat (Peeler et al., 2000). Moreover, a herd risk factor associated with IRCM could be due to several reasons: a high culling rate could be an indication of one or more of the following: the relevant herd was not stable; the herd size was increasing; there was increased culling for some reason, including mastitis; or that the introduction of cows into an established herd increased the risk of mastitis (Peeler et al., 2000).
Some studies suggested that low SCC was associated with high IRCM (Elbers et al., 1998; Beaudeau et al., 2002; Green et al., 2004), whereas other works did not reveal any interrelationship (Barkema et al., 1998; Beaudeau et al., 1998). These controversial results could depend on: 1) the situation of the herds in terms of prevailing pathogens (Barkema et al., 1998; Peeler et al., 2003); 2) the accuracy of CM evaluation including the epidemiological unit (Gasqui and Barnouin, 2003; Peeler et al., 2003), SCC measurement level (from most to least accurate: quarter, cow, bulk), and length of the period taken to observe CM cases after SCC measurement (Beaudeau et al., 2002); 3) the stability of management practices in the herds; 4) the fact that the farmers who are better at diagnosing CM might divert more high SCC milk (Elbers et al., 1998); 5) the observed (negative) association between milk yield level and SCC (Barkema et al., 1998), and (positive association) between milk yield level and CM risk (Schukken et al., 1990; Chassagne et al., 1998).
The purpose of the present epidemiological study was to identify, from a large set of potentially relevant variables, dairy herd management practices and characteristics associated with IRCM in French herds with low SCS fulfilling the following criteria: 1) selected on a national basis from DHIA database; 2) experienced a low SCS over a long period; 3) recorded CM cases through a standardized system, which included individual forms collected monthly; and 4) had a stable dairy management system.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Herd Selection
The general herd selection procedure of ZMP has been described previously (Barnouin et al., 2004). The surveyed herds were located in 48 French departments, had at least 20 cows of which 90% were Holstein, Montbéliarde, or Normande breeds, and did not practice vaccination against mastitis. The herd sample was constituted from the National DHIA database according to the herds previous 36-mo history of somatic cell score (36-mo SCS), because SCS (SCS = log2[SCC/100,000] + 3) is, contrary to bulk SCC, a criterion not biased by milk discarding. Herd 36-mo SCS was the arithmetic mean of all monthly cow SCS values determined in a herd during the 36 mo preceding ZMP. All the selected herds of the present study had 36-mo SCS that were among the lowest 5% of SCS for herds within breed, as breed is a key SCS variation factor (Rupp et al., 2000). Consequently, 36-mo SCS was
1.99 (125,000 cells/mL) for Montbéliarde herds,
2.38 for Holstein herds (185,000 cells/mL), and
2.76 (230,000 cells/mL) for Normande herds. Three hundred fifty-seven herds were enrolled, representing 15% of the total number of herds in the National DHIA database that had 36-mo SCS in the lowest 5%. The study did not involve any payment to the selected farmers. Two farmers quit the survey before starting. Five herds were removed because of bovine spongiform encephalopathy. Thirty farmers quit ZMP before it ended for different reasons (DHIA resignation, health problems, lack of time, unknown reasons) or were discarded for incomplete data concerning mastitis recording or questionnaires. Moreover, 23 herds were removed because the farmers did not state the start and end dates of their mastitis recording. In the end, 297 herds were studied, representing 83.2% of the herds initially enrolled.
Clinical Mastitis
Clinical mastitis was identified by the farmers based on clinical signs including udder inflammation or abnormal milk, with or without general clinical signs. Farmers were instructed how to fill out a specific record form and required to describe all the observed CM cases from February 1, 1999 to July 31, 2001. Data included: CM case definition, cow identification, date of CM occurrence, clinical description, affected quarters, and treatment. Data forms were collected and checked monthly by DHIA technicians. To standardize the lag time for a new case, multiple clinical episodes within 8 d were counted as only 1 case (IDF, 1997). Cow-days at risk were calculated as the total number of days that a cow was present in a herd during the observation period, starting 30 d before first calving (IDF, 1997), minus 8 d after each case to correct time at risk for lag time (IDF, 1997). The IRCM was expressed as the reported number of CM episodes during the observation period per 100 cows/yr at risk.
Questionnaire Data
The DHIA technicians collected questionnaire data concerning characteristics and management practices of the farms by interviewing the farmers. Milk samples were also collected through DHIA to determine milk yield and SCC and other milk quality measures. Bacteriological analyses of milk were not available in the study. The questionnaire variables aimed to describe all aspects of the dairy management system. The questionnaires were validated by consensus of an expert working group and pilot interviews of 10 farmers outside ZMP to check its feasibility and comprehensiveness. Questions not easily understood were excluded. In total, the final version of the questionnaire consisted of 1055 herd characteristics and practices, of which 188 were determined by the technician from preexisting information (farm general descriptors, monthly calving numbers, milking machine characteristics, milk production, and quality data). The interviewees did not have to answer every question because some did not apply to their farm (e.g., free-stall characteristics could not be described when the cows were loose-housed in a straw yard). It required 7 to 8 h to collect all the questionnaire data. Two complementary questionnaires were administered. The first of these questionnaires, concerning the characteristics of the dairy cow shed and milking parlor, character traits of the herdsman, and measures of milk quality, was carried out from October to December 1999. The second questionnaire, which concerned housing conditions of the herd, the feeding system, udder and calving hygiene, and herd health, was carried out from October to December 2000. Interviewers were the DHIA technicians of the farms studied and they had a thorough knowledge of the herd management system; at the end of the interview, they checked the farmers answers by requesting confirmation of any information that seemed erroneous. Finally, certain questions concerning 15 herd characteristics and practices, gave rise to answers that were inconsistent with the researchers knowledge of dairy farm management in France. These questions were removed from data analyses. The farmers declared at the end of the study that their dairy management system had not changed during ZMP. This statement was checked by the DHIA technicians and we therefore assumed that herd characteristics and practices were stable between the times of the 2 questionnaires.
Database
All CM and questionnaire data were stored in an Access 2000 (Microsoft France, Courtaboeuf, France) relational database organized through 34 tables. Herd identification number, breed, SCS data, and other milk traits were extracted from the national DHIA database. Checking procedures for data validity were conducted, which consisted of inspecting the CM record forms and questionnaire answers and checking with the interviewer that such answers were normal or incomplete. A set of automatic logical procedures was implemented to detect erroneous entries (e.g., response concerning freestalls if the cows were housed in a straw yard, or number of weeks per year in which the herd was outdoors exceeding 52 wk, mastitis clinical cases without any treatment described). Moreover, to check for typing errors, a double keyboard data entry was performed independently by 2 technicians.
Statistical Model
The association between herd IRCM and potential risk factors was screened using Poisson regression models through the SAS/STAT 8.1 GENMOD procedure (SAS Institute, 1999). Because the variance of IRCM distribution was larger than the mean in the present data, an overdispersed Poisson model with a log link was built (Barkema et al., 1999). With this model, it was assumed that the variance was a constant multiple of the mean instead of equal to the mean (Dohoo et al., 2003). The overdispersion parameter was estimated by the deviance divided by its degrees of freedom. The linearized model used was as follows:
![]() |
where f(NCM) = link function of herd number of cases (NCM) of clinical mastitis during the observation period,
= intercept, LNDAY = natural log of the number of days of observation in a herd during the observation period, calculated as the sum of days of observation for each present cow (used as an offset variable), ß = regression coefficient for RF, RF = herd risk factor, and e = Poisson distributed residual random error. Estimated regression coefficients of the model may be expressed exponentially and interpreted as a relative risk or rate ratio.
Modeling Process
The modeling process included 3 steps. In the first step, all single explanatory variables were screened in a bivariate Poisson model and variables with P
0.25 in the likelihood ratio test were selected, as potential candidates for further multiple Poisson regressions. In the second step, the selected variables were grouped within 8 categories of management practices (Table 1
). The correlations between the selected variables were analyzed for each category of management practices using the CORR procedure. If variables were correlated (correlation coefficient >0.15 or < 0.15), only the variable with the best fit through the bivariate Poisson model of the first step was included in the corresponding multivariate category model. Consequently, it was reasonable to expect that the multicollinearity between variables was notably weak (Barnouin et al., 2004).
|
0.10 selected from the 8 category models using backward elimination were offered to a final Poisson model. The correlations between the offered variables were analyzed as in the second step. The variables that were significant (P
0.05) through the final Poisson model were considered the management variables associated with IRCM. Biologically plausible interactions among significant variables were tested, and none of those interactions was significant. The goodness of fit of the model was assessed by standardized deviance residual analysis (McCullagh and Nelder, 1989). Moreover, the Grubbs test (Grubbs, 1969) was used to detect potential outliers. | RESULTS |
|---|
|
|
|---|
Clinical Mastitis
A total of 5022 CM cases covering 20,084 cows and 37,393 lactations were observed during the program. The average IRCM was 20.1 cases/100 cows per yr (SD = 16.1). This ranged from 0 to 116. Distribution of IRCM per herd is presented in Figure 1
. The IRCM was lower than 20 cases/100 cows per yr in 60% of the surveyed herds and higher than 40 cases/100 cows per yr in 8% of the studied herds.
|
|
|
0.001); b) more than 2 person-years assigned to dairy herd management; c) balanced concentrate in the cow basal diet (b to c: P
0.01). Herds with the following characteristics had higher IRCM: a) milking cows loose-housed in a straw yard; b) no mastitis therapy performed when a single clot was observed in the milk (a to b: P
0.001); c) clusters rinsed using water or soapy water after a high SCC milked cow; d) 305-d milk yield >7435 kg; e) herd located in the South region; f) herd located in the North region (c to f: P
0.01); g) cows with at least 1 nonfunctional quarter; h) premilking holding area with a slippery surface (g to h: P
0.05). Moreover, 23 variables removed from the statistical analysis were correlated with the significant variables of the final model (Table 3
|
| DISCUSSION |
|---|
|
|
|---|
The present study aimed to test many management characteristics. Consequently, an effort was made through study design and statistical procedure to reduce the multicollinearity among explanatory variables, a relatively frequent situation within the dairy management system. The categorization of risk factors was performed to facilitate the detection of interrelated management practices, although categorization is often subjective (Schukken et al., 1990). To lower the collinearity among explanatory variables, a set of correlated variables was removed from the analysis according to statistical reasons. As items for discussion, we summarized (Table 3
) the variables that were removed because they were correlated with the significant variables of the final model. Nevertheless, the fact that a model fits does not ensure the model contains the best biological information. Moreover, in epidemiological studies, the observed association between a factor and a disease did not display an obviously causal relationship (Schukken et al., 1990). But most of the time, it is difficult to conduct rigorous intervention studies on private farms to demonstrate the true relationship between herd management practices and mastitis risk, as has, nevertheless, been implemented by Lam et al. (1997).
From previous ZMP results (Barnouin et al., 2004), 18 major variables allowed discrimination between very low and medium SCS herds, whereas 10 variables were associated with IRCM in the very low SCS herds of the present study. Two variables, which both characterized medium SCS herds vs. low SCS herds and were associated with IRCM in the present study, appeared as mastitis risk factors: 1) milking cows loose-housed in a straw yard; and 2) no mastitis treatment when a single clot was observed in milk at successive milking. Moreover, 2 variables were risk factors for 2 epidemiological situations: 1) cows culled when found with at least 1 damaged teat characterized low SCS herds, whereas cows culled after 3 clinical cases within a lactation characterized low SCS herds that had lower IRCM; and 2) less than 1 person assigned to extra dairy management activities characterized low SCS herds, whereas more than 2 persons assigned to dairy management characterized low SCS herds that had lower IRCM.
Concerning the housing system, which was correlated with numerous variables in the present study (Table 3
), our results agree with those of Peeler et al. (2000), who demonstrated that straw yards resulted in a significant increase in IRCM in low SCC herds. In a study of the incidence of CM in straw yards and freestalls, Fregonesi and Leaver (2001) found less CM cases in freestalls. Straw is a bedding material favoring rapid growth of the main udder pathogens, particularly in humid weather conditions and in bedding kept moist through urine excretion (Ward et al., 2002). As a major part of the udder environment, hygiene and bedding play an important role in the transmission of environmental pathogens implicated in CM (mainly E. coli and Strep. uberis), but also in the development of contagious pathogens (mainly Staph. aureus) with subsequent CM and subclinical cases (Elbers et al., 1998). Moreover, recurrent CM episodes with a transient subclinical stage and caused by E. coli through transmission from one quarter to another, could contribute significantly to increased IRCM in low SCC herds (Döpfer et al., 1999), in which E. coli is considered the predominant pathogen causing CM (Barkema et al., 1998).
Regional differences in SCS have been shown previously in the United States (Miller et al., 1999; Ely et al., 2003) and variation in CM risk in France could depend on meteorological factors (Barnouin et al., 1986). Through the present study, the herds that were located in the North and in the South regions have been shown to have a higher IRCM. In France, these 2 regions are less specialized in bovine production than the West, East, and Center regions. Consequently, fewer well-trained dairy management experts would serve in the North and South regions (Chassagne et al., 2005). Nevertheless, as fresh ryegrass in the diet and farm altitude were correlated with herd location in the present study (Table 3
), some other undetermined reasons could explain the observed regional differences in CM risk.
In ZMP, the presence of a single clot observed in the milk was not considered a mastitis case by herdsmen. Nevertheless, CM therapy was practiced in some farms after observation of a single clot to prevent mastitis. As IRCM were higher in the farms in which such a practice was not performed, the presence of a single clot in the milk could be considered a precursory signal that would precede CM occurrence. Actually, early treatment of mastitis reduced the severity of the disease and more easily brought about a bacteriological cure (Milner et al., 1997). In particular, early intervention using an intramammary antibiotic would achieve a cure for all infections caused by Strep. uberis (Hillerton and Semmens, 1999). Although changes in the milk electrical conductivity and biochemical indicators can help detect early udder infection (Milner et al., 1997), none of the automatic detection methods is, at present, satisfactory, and false alerts are frequent (Pyörälä, 2003). The ZMP results would indicate that a conventional detection method based on identification of clots by the milker could be an effective way to prevent CM and SCM when treatment is initiated after observation of a single clot. However, to be able to control CM, farmers with low SCS herds must be trained not only to detect clots, but also to treat cows efficiently and hygienically (Barnouin et al., 2004), and to separate them from the rest of the herd to reduce spreading contaminated milk from infected to healthy cows.
Culling a cow with 3 CM occurrences within a lactation, that is, within a relatively short time of first becoming infected, would limit the transmission of potentially invasive bacteria able to cause CM (White et al., 2003). Moreover, this practice could have a CM-reducing effect by removing cows that might give rise to cases in the future. The influence on IRCM of such a culling measure, which might depend on the pathogen that caused the CM, has not been previously studied in epidemiological surveys conducted in low SCC herds. Through ZMP, culling practices would be effective to decrease mastitis risk in 2 different ways, as culling cows with damaged teats characterized low SCC herds (Barnouin et al., 2004), whereas culling cows with multiple cases of CM within a lactation was associated with a lower CM risk in the present study.
According to Barnouin et al. (2004), herdsmen must be specialized in dairy farming for long-term control of herd SCS. To decrease IRCM in low SCS farms, at least 2 person-years, in addition to herd specialization, must be assigned to dairy management in French farming conditions, including smaller herd sizes than in some countries. Extra persons are needed because adequate and permanent CM detection, treatment, and management of treated animals require much more time than for herd SCC management. Interaction between the cowman and his animals would define a competent dairy farmer, the regularity of contact being important to develop and maintain satisfactory human-animal interactions (Raussi, 2003). Such satisfactory contacts would decrease stressful breeding conditions, reduce udder injuries, and avoid stress-dependent alteration in immune responsiveness (Cai et al., 1994). Moreover, in low SCC herds, Holtenius et al. (2004) suggested that there were differences in the metabolism and immune status between herds with high or low mastitis incidence indicating increased metabolic stress in high-incidence herds. Finally, frequent contacts between herdsmen, technicians, and cows would promote quick therapeutic decisions following any change in udder or milk appearance, and would favor accurate mastitis control (Barnouin et al., 2004).
A premilking holding area with a slippery surface, which was displayed as a risk factor for CM, could increase udder injuries and stress, with consequences as previously discussed. Avoiding slippery holding areas is a recommended practice in some DHIA and prevention programs, but such a practice was not previously displayed as characterizing very low mastitis risk herds. Moreover, in the high IRCM herds, the clusters were more frequently rinsed using water or soapy water after a cow with a high SCC. The farmers in the present study, as good dairy managers, tried to reduce risk of CM through cluster rinsing when a high SCC was determined exceptionally in a cow of their low SCC herd. Such a practice did not affect the incidence of staphylococcal, streptococcal, and gram-negative bacterial infections, and clinically infected quarters (Smith et al., 1985). Cluster rinsing leads only to a marginal reduction of bacterial load, whether using water or disinfectant (Shearn et al., 1994).
A basal diet containing a commercial concentrate with a balanced energy/protein ratio could contribute, especially in the early lactation, to fewer metabolic disorders that could be linked to risk of mastitis (Holtenius et al., 2004). In low SCC cows, sugar beet pulp as concentrate was associated with a higher IRCM (Schukken et al., 1990), whereas Johnson and Otterby (1981) found higher energy supplies at drying-off in dairy cows that subsequently developed CM. Dietary polyunsaturated fatty acid levels could modulate the leukocyte activity in the milk and inflammation in the udder (Barnouin et al., 1995). Consequently, adequate energy supplies through basal diet and concentrate could be protective factors against CM risk in low SCC herds, although a clear explanation of the etiological mechanism is not available. From the results of ZMP, higher milk production was associated with a high IRCM. Interestingly, this variable also characterized high IRCM herds through univariate (Schukken et al., 1990) or multivariate (Peeler et al., 2000) statistical analyses concerning CM factors in low SCC farms. For Peeler et al. (2000), the risk threshold for an increased CM risk was 7500 kg/cow per yr, similar to current results (7436 kg). The positive association between milk yield level and CM risk was previously demonstrated in French Holstein herds at the cow level (Chassagne et al., 1998), and such a relationship was found in a study conducted in Ontario Holstein herds (Bigras-Poulin et al., 1990). Through the present study, a high correlation has been displayed between milk yield level and importance of corn as a forage (Table 3
), a feedstuff high in calories but rather poor in several amino acids, minerals, trace elements, and vitamins. Moreover, 8 different dietary components, displayed in Tables 1
and 3
as variables selected at the first analysis step or correlated with the significant variables of the final model, might indicate that the relationship between dietary components and CM risk should be more deeply investigated in low SCC herds. Some management practices commonly related to CM risk were not associated with the disease in the present study. In the present study, vaccination against mastitis was not performed in the herds, nor was the use of internal sealants. Moreover, dry cow environment, milking technique, teat dipping, and calving conditions were associated in ZMP to subclinical mastitis risk (Barnouin et al., 2004; Chassagne et al., 2005).
In countries where intensive breeding conditions seem the only way to produce milk profitably, increased CM risk through high milk production could be considered as part of the natural cost of production, as is concentrate feeding or 3x milking in high output herds (Chassagne et al., 1998). Nevertheless, genetic, metabolic, conformation, and behavioral evaluations should be conducted in high-yielding and very low SCC herds to evaluate if production level and CM traits are definitively antagonistic through the udder, an organ so delicate and sensitive, but also efficient. Such studies should focus on clinical cases during the postpartum period and consider dairy management practices characterizing very low mastitis risk from ZMP (Barnouin et al., 2004) and from other mastitis control programs conducted in low SCC herds (Schukken et al., 1990; Peeler et al., 2000). Such studies could profitably involve scientific collaborations including epidemiological studies implemented at the multicountry level and be a further step for a better understanding of mastitis.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication February 24, 2005. Accepted for publication June 21, 2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. G. M. Olde Riekerink, H. W. Barkema, D. F. Kelton, and D. T. Scholl Incidence Rate of Clinical Mastitis on Canadian Dairy Farms J Dairy Sci, April 1, 2008; 91(4): 1366 - 1377. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Rodrigues-Motta, D. Gianola, B. Heringstad, G. J. M. Rosa, and Y. M. Chang A Zero-Inflated Poisson Model for Genetic Analysis of the Number of Mastitis Cases in Norwegian Red Cows J Dairy Sci, November 1, 2007; 90(11): 5306 - 5315. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |