|
|
||||||||

* School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom, EH9 3JT
Holstein UK, Rickmansworth, Hertfordshire, United Kingdom, WD3 3BB
1 Corresponding author: o.m.onyiro{at}sms.ed.ac.uk
| ABSTRACT |
|---|
|
|
|---|
Key Words: genetic parameter digital dermatitis locomotion trait dairy cow
| INTRODUCTION |
|---|
|
|
|---|
Digital dermatitis has been associated with decreased milk yield, reduced reproductive performance, increased involuntary culling rate, and reduced general well-being of the animals (Nutter and Moffit, 1990; Argaez-Rodriguez et al., 1997; Garbarino et al., 2004). Wells et al. (1999) reported that 9.7% of DD resulted in lameness in affected cows. Hernandez et al. (2001) reported that calving to conception interval was significantly longer in lame cows with multiple lesions including DD, compared with healthy cows. Earlier, Argaez-Rodriguez et al. (1997) reported that cows affected by DD in a Mexican herd had a 20-d increase in calving to conception interval, resulting in reduced fertility and consequently a shorter lifespan. The same study also reported a decrease in the milk production of cows affected with DD compared with healthy cows, although the difference was not significant. Corrie et al. (2000) noted that papillomatous digital dermatitis was common among culled adult cattle with the incidence greater among dairy cattle compared with their beef counterparts.
Estimates of genetic parameters between claw diseases including DD and locomotion traits were reported by van der Waaij et al. (2005). With threshold and linear models, they reported a heritability of 0.10 from both models for DD in Dutch Holstein-Friesian dairy cows. They also reported genetic correlations of 0.16 ± 0.13, –0.22 ± 0.13, –0.67 ± 0.19, and –0.34 ± 0.12 between DD and rear legs, side view (RLS), foot angle (FA), locomotion score (LOCO), and leg and feet composite (L&F), respectively. Approximate genetic correlations reported by Koenig et al. (2005) between DD and linear type traits ranged from 0.03 (hocks) to –0.61 (FA).
Reports of genetic parameters of DD and genetic associations with locomotion, production, fertility, and longevity traits are few in the literature, yet these estimates are required to enable DD to be included in a genetic selection index. The genetic association of DD with longevity has not been reported so far in the literature.
Studies on environmental and management risk factors have associated the incidence of DD with herd size, stage of lactation, parity, damp floors, housing systems, and so on. (Frankena et al., 1991; Rodriguez-Lainz et al., 1996; Somers et al., 2005). Cows housed in cubicles and those with restricted access to pasture were observed to be at greater risk of DD compared with cows that had full access to pasture during summer (Frankena et al., 1991; Somers et al., 2005). Moist and unhygienic floor conditions as well as solid or grooved concrete floors have also been associated with a greater incidence of DD (Rodriguez-Lainz et al., 1996; Wells et al., 1999; Hultgren and Bergsten, 2001).
As part of the UK type classification scheme, field officers now routinely collect information on housing systems and flooring conditions (Onyiro and Brotherstone, 2008). In addition, the field officers record whether cows exhibit signs of DD during type classification. It is, therefore, now possible in the UK to investigate associations between DD and housing conditions and estimate genetic parameters for DD and associated traits using national data.
The objectives of this study were 1) to determine the influence of housing on the occurrence of DD in Holstein-Friesian dairy cows; 2) to estimate the heritability of DD and its genetic correlation with locomotion and other leg and feet traits; and 3) to estimate approximate genetic correlations among locomotion traits, lifespan (LS), production and fertility traits, and DD from sire EBV.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Five type traits related to legs and feet were evaluated. The traits are LOCO, RLS, FA, and bone quality (BONEQ) measured on a scale of 1 to 9, which represents extreme biological values, and L&F, which was measured on a subjective scale of "poor" (65) to "excellent" (95). Table 1
gives a description of the traits scored with their means and standard deviations, and Table 2
shows the locomotion scoring system used for UK national herds. Age of cows at inspection was between 19 and 50 mo.
|
|
Statistical Analysis
Data were first adjusted for field-officer differences in scores by scaling records such that the standard deviation of each field officer was equal to the mean standard deviation of all field officers (Brotherstone, 1994).
The model equation used to describe the data was

where Yi = presence (1) or absence (0) of DD for animal i; µ = overall mean; hyv = random herd-year-visit effect; β1 and β2 = linear and quadratic regression coefficients of DD on age of cow at inspection;
1 and
2 = linear and quadratic regression coefficients of DD on stage of lactation of cow (stage) at inspection;
= linear regression of DD on the proportion of Holstein genes; phols = proportion of Holstein genes; moc = effect of month of calving; hcode.hmns = housing period effect nested in housing type; hfloor = flooring condition at inspection of animal i; ai = random animal effect; and ei = random error term.
The length of time cows spent in a particular housing was classified as 1, 2, 3, 4, 5 mo and then 6 or more due to small sample sizes recorded from 6 mo. The length of time cows spent in slatted or loafing yards was not accounted for because few animals were recorded in these housing types. Housing period was nested within housing type and the combined effect referred to as housing situation. Cow and herd-year-visit were fitted as random effects in the model. Fitting herd-year-visit as a random effect allowed us to estimate the effect of housing systems and flooring conditions as these factors are confounded with herd visit. Herd-year-visit in this model will account for residual differences after the fixed effects have been accounted for. Housing systems and flooring conditions are both indicative of the overall farm management conditions. Age of cow at inspection, stage of lactation at inspection, and mean proportion of Holstein genes (ranging from 91 to 100%) were included as covariates.
Heritability of DD was estimated from a linear model. The heritability obtained from the linear model was then transformed from the 0/1 scale to an assumed underlying continuous normally distributed scale (Robertson and Lerner, 1949). The transformation applied was
![]()
where hn2 = heritability on the continuous normal scale; h0/12 = heritability on the binomial scale; z = ordinate of the standardized normal distribution at the threshold point corresponding to p, and p = incidence of DD; that is, the proportion of cows that had DD.
Genetic and phenotypic correlations between DD and type traits were estimated from a multivariate REML analysis using an animal model. The same model was applied to both DD and the linear traits.
A single multivariate analysis with DD and all 5 type traits was not feasible because insufficient computer memory was available. Therefore, a series of 2-trait analysis (i.e., 5 bivariate analyses) was performed with DD and one of the type traits in turn. No structure was imposed on the variance-covariance matrix. Phenotypic correlations between DD and the type traits (r0/1) were also transformed to an assumed underlying normally distributed scale (rn) by applying the method of Ollausson and Ronningen (1975):
![]()
where p and z are as above. Genetic correlations on the 0/1 and underlying normal scale are expected to be similar, so no transformation is needed (Ollausson and Ronningen, 1975). The ASREML software was used for all analysis (Gilmour et al., 2000).
Estimation of Approximate Genetic Correlations Among Locomotion Traits and Production, Fertility Traits, and Lifespan
Sires EBV were correlated to obtain approximate genetic correlations among type traits (LOCO, RLS, L&F, FA and BONEQ), production traits (milk and fat), fertility traits [calving interval in days (CI) and 56 d nonreturn (NR56), scored 0 if the cow returns to service within 56 d and 1 otherwise], and lifespan. Wall et al. (2003) describe the UK system to predict EBV for fertility and Brotherstone et al. (1997) detail the genetic evaluation of lifespan. The phenotypic lifespan trait is the actual lifespan of the cow, or if she is still alive, her predicted lifespan. Sire EBV for type traits were taken from the output of the analysis described above and their corresponding reliabilities were estimated from the standard errors of the sire solutions as follows:
![]()
where RLi = sire reliability for trait i; SEi = standard error of sire EBV for trait i; and
a(i)2 = additive genetic variance for trait i.
Sire EBV and reliabilities for milk, fat, CI, NR56, and longevity were taken from the UK national evaluation records of May 2007. Sires with missing EBV for any of the traits of interest were eliminated. Minimum reliability for type traits used in estimating the correlations was 0.30, resulting in a minimum reliability for production traits, CI, NR56, and lifespan of 0.50, 0.34, 0.35, and 0.32, respectively. Correlations were based on 973 sires.
Estimation of Approximate Genetic Correlations Between DD and Lifespan, Production, and Fertility Traits
As described above, sire EBV were correlated to obtain approximate genetic correlations between DD and LS, milk, fat, CI, and NR56. Again, sire EBV for DD were taken from sire solutions and their reliabilities calculated from the standard errors of the solutions. The low heritability of DD resulted in many bulls having EBV of low reliability. Therefore, the minimum reliability requirement for DD was reduced from 0.30 to 0.10 which resulted in a lower limit of 0.28 (LS), 0.41 (production traits), and 0.27 (fertility traits), respectively. Correlations were based on 2,461 sires.
All correlations estimated from EBV were adjusted for their reliabilities by the method of Calo et al. (1973):

where RL1 and RL2 = reliabilities of traits 1 and 2;
g1,2, = approximate genetic correlation between traits 1 and 2; and r1,2 = correlation between EBV for traits 1 and 2.
The standard errors of correlations were estimated using the formula below:
![]()
where n = number of sires with records (Sokal and Rohlf, 1995).
To investigate whether genetic associations between type and fertility were mediated through production, approximate genetic correlations between these traits were adjusted for milk yield. Similarly, the correlations between DD and fertility traits were adjusted for milk yield. The partial correlation coefficients were estimated using the following formula (Sokal and Rohlf, 1995):

where
1,2,3° = partial correlation coefficient among DD, CI, and NR56 adjusted for milk yield;
g1,3, and
g2,3, = approximate genetic correlations of trait 1 (DD) with milk yield and trait 2 (CI and NR56) with milk yield.
The partial correlations were tested for significance with a t-statistic as shown below (Sokal and Rohlf, 1995):

where m = number of variables kept constant, which is 1 in our case—milk yield.
Standard errors for the partial correlations were estimated in the same way as standard errors of ordinary correlations.
| RESULTS |
|---|
|
|
|---|
|
|
Housing Situation and Flooring Effect
Housing situation and flooring condition removed significant amounts of variation in DD. Predicted effect of housing situation (relative to pasture for all months) and flooring condition (relative to flooring condition 5) on the incidence of DD is given in Table 4
. Cows housed in cubicles had significantly (P < 0.05) increased DD lesions as the housing period increased compared with cows spending similar periods on pasture. This result is consistent with the mean incidence of DD within each housing period (Table 3
). Cows kept in straw yards showed significantly (P < 0.05) greater signs of DD compared with those on pasture when housed for at least 6 mo. The results generally show that cows on pasture were less susceptible to DD than cows in other housing systems.
|
Heritabilities and Correlations Between DD and Type Traits
Estimates of heritabilities of DD and genetic correlations between DD and type traits, together with their standard errors, are given in Table 5
. The genetic variance of DD was significantly greater than zero and the estimate of heritability was 0.011 (±0.003). After transformation to the underlying scale, the heritability of DD increased to 0.029 (±0.007).
|
The phenotypic correlations between DD and the type traits were lower than their corresponding genetic estimates and were all significantly (P < 0.05) different from zero. Only RLS had a positive phenotypic association with DD. The result suggests that, phenotypically, cows with no signs of DD have higher linear type scores and better legs and feet. When transformed from the 0/1 scale to the underlying continuous scale, the phenotypic correlations also increased (Table 5
).
Approximate Genetic Correlations Among Type Traits, LS, Milk, Fat, CI, and NR56
Table 6
gives the approximate genetic correlations of the type traits with LS, production, and fertility, together with standard errors. The correlations between type traits and LS were all significantly different from zero. Lifespan had high approximate genetic correlations with BONEQ (0.50), LOCO (0.66), and L&F (0.69), suggesting that flat and refined bones, good locomotion, and better legs and feet are associated with increased longevity. Foot angle was moderately correlated with LS; evidence that cows with a steeper foot angle had increased longevity. The correlation between RLS and LS (–0.32) was negative, indicating that sickled rear legs as judged from the side are associated with reduced lifespan.
|
Generally, the approximate genetic correlations of type traits with NR56 were low, except for that between BONEQ and NR56 (0.36) but significantly different from zero. The correlations indicate that straighter RLS, steeper FA, better L&F, and flatter, more refined bone quality are associated with better conception rate. Correlations with CI were less intuitive but these correlations are subject to management influences. The partial correlations between type and fertility traits adjusted for milk yield are presented in Table 7
. The partial correlation coefficients were all significant except between L&F and NR56, indicating that only the association between L&F and NR56 is mediated through milk yield.
|
Milk and fat yield had moderate genetic associations with CI, and the correlations were similar (0.27) and significantly different from zero implying that longer CI is associated with increased milk and fat production. Milk, fat, and CI were negatively correlated with NR56, suggesting that high milk or fat yield or longer CI is linked with a lower conception rate.
Approximate Genetic Correlations of DD with LS, Milk, Fat, CI, and NR56
Approximate genetic correlations of DD with LS, milk, and fat were negative and moderate suggesting that the presence of DD is associated with reduced longevity and decreased milk and fat yield. Thus, breeding for resistance to DD should improve both production and longevity.
The association between DD and CI (–0.07) was low and negative, indicating that cows affected by DD had slightly shorter CI than cows with no evidence of DD. From the partial correlation between DD and CI (0.015), it was seen that the shorter CI is mediated through the reduced milk yield of affected cows, because the association was not significantly different from zero. In addition, DD and NR56 were moderately correlated (0.48), indicating an association between DD and greater conception rates. The partial correlation coefficient (0.40) between DD and NR56 (Table 7
) was still significant, suggesting the correlation is not mediated through milk yield.
| DISCUSSION |
|---|
|
|
|---|
Association Between Housing System, Flooring Condition, and DD
Results of this analysis show that housing dairy cows in cubicle houses predisposes them to DD infection, and DD increases as housing periods increase. High incidence of DD has been reported for cows kept in cubicles (Frankena et al., 1991; Somers et al., 2005). Faye and Lescourret (1989) also showed that claw health was worse in indoor cubicle housing than during pasturing. Cubicle houses are characterized by abrasive solid concrete, which causes wear and tear of the cow hooves, thus exposing the inner claws to microorganisms causing infectious diseases such as DD and other injuries. Results from the current study also showed that straw yards predispose cows to DD when housed for extended periods. Straw yards are softer environments than cubicles, and with good management such as daily removal of bedding (which prevents wetness and pathogen build-up), cows should be less vulnerable to DD. There was no significant difference in DD between cows in slatted or loafing yards and those on pasture. However, because the amount of time spent by cows in this system was not accounted for in our analysis, we cannot ascertain whether slatted/loafing yards are compatible with pasture in terms of reduced incidence of DD with longer housing periods. Moreover, a small number of animals were recorded in this system compared with the number of animals at pasture.
Cows with full access to pasture have been found to have a reduced incidence of DD (Wells et al., 1999; Somers et al., 2003, 2005). Our findings corroborate these reports. Figure 1
clearly indicates that cows become less vulnerable to DD as they spend more time on pasture. The high incidence of DD experienced by cows on pasture in the first month may be a result of the housing type where cows were housed in the previous month. The association between flooring condition and incidence of DD was counterintuitive. The lower incidence of DD in dirty, slippery floor conditions compared with the other flooring conditions could be due to difficulty in identifying lesions on dirty feet. As a result, many of the DD cases might go undetected by field officers. Wells et al. (1999) reported a lower incidence of DD on dirt, pasture, or smooth concrete compared with grooved concrete.
Estimates of Genetic Parameters
Heritability
The heritability estimates for DD in this study were lower than previous reports in literature. Using a sire model, van der Waaij et al. (2005) reported an estimate of 0.10 (linear and threshold) for DD in Dutch Holstein-Friesians whereas Koenig et al. (2005) obtained an estimate of 0.073 from a linear logistic mixed model. The disparity in these heritability estimates could be attributed to the differences in the number of records as well as the models used in data analysis. Although our estimates of heritability were low, they were significantly different (P < 0.05) from zero, implying that DD is heritable, and genetic improvement for resistance to DD is possible through selection.
Genetic Correlations
The association between DD and LOCO (–0.67) is similar to the estimate of van der Waaij et al. (2005). The genetic relationships between DD and LOCO and L&F were of similar magnitude in our study and this is not surprising, because both traits have been reported to be highly genetically and phenotypically correlated (van der Waaij et al., 2005; Onyiro and Brotherstone, 2008). The correlations among DD and type traits relating to legs and feet suggest better legs and feet if cows are bred for increased resistance to DD. Conversely, selection for improved type traits would lead to increased resistance to DD. Nonsignificant associations among DD, FA, and RLS were noted by van der Waaij et al. (2005).
Approximate Genetic Correlations
Correlations among Type, LS, Production, and Fertility Traits
The strong association between the type traits and LS could be linked to the fact that improved conformation particularly related to legs and feet will make cows less vulnerable to locomotive disorders. This is consistent with the correlations between DD and type traits (Table 5
). Dekkers et al. (1994) reported moderate genetic association among BONEQ, L&F, set of rear legs on a desirability scale and functional herd life.
The relationships between FA, RLS, L&F, and BONEQ suggest that lower FA, more sickled rear legs, poorer L&F score, and coarser bones are associated with increased milk production. Perez-Cabal et al. (2006) also found that sickled rear legs were associated with increased milk yield. However, they detected a positive association between better legs and feet and greater production, whereas we found that poorer legs and feet score was associated with greater milk yield. In this analysis, we found no significant association between LOCO and milk yield.
In general, the relationships between type and CI are quite surprising but could be due to management effects. Pryce et al. (2000) obtained a negative correlation of FA with CI (–0.20), but positive with RLS (0.19). Because fertility remains an economically important trait in the breeding goals of farmers, and selection for productive life is expected to result in improved health and fertility, the associations between type and fertility traits warrant further research.
Correlations Among LS, Production, and Fertility Traits
The unfavorable approximate genetic correlation between LS and milk yield is an indication that daughters of sires with high EBV for milk yield will have reduced longevity. The correlations among LS, CI and NR56 were favorable (–0.49 and 0.39, respectively). Cows having shorter CI and successful conception will stay longer in the herd. The antagonistic relationship between CI and NR56 (–0.34) means that a longer CI is associated with a reduced pregnancy rate.
The association between milk, fat, and fertility traits suggests that high-yielding cows will have longer CI and impaired conception leading to reduced reproductive fitness. Poor fertility and health along with low production have been identified as reasons for involuntary culling (Esslemont and Kossaibati, 1997). Many studies have reported similar genetic relationships between production and fertility traits. Wall et al. (2003) reported antagonistic correlations between milk yield and CI (0.27) and NR56 (–0.45). Brotherstone et al. (2002) estimated unfavorable genetic correlations between combined 305-d fat and protein and CI (0.40) and NR56 (–0.29). These findings support the inclusion of fitness related traits with production in national selection indices (Miglior et al., 2005).
Correlations Among DD, LS, and Production.
Correlation of the EBV adjusted for their reliabilities among DD, LS, and milk and fat (–0.16, –0.31, –0.43, respectively), suggest that increased incidence of DD is associated with reduced LS and decreased milk and fat production. Contrary to our estimate, Koenig et al. (2005) obtained a positive correlation between DD and milk yield in the first stage of lactation (0.24) but the analysis involved a small data set.
As the presence of DD was identified by field officers rather than veterinarians and was scored only once in first lactation, the incidence of the disease as estimated from our data may be biased downwards. However, recording the presence of DD on the national population (around 60,000 heifers per year are type classified in the UK) should give sufficient information to predict sire breeding values for DD of reasonable reliability. Results from our analyses indicate that breeding for increased resistance to DD is associated with an increase in both longevity and production. Thus, DD should be considered a useful disease trait to be recorded and included in a national selection index for increased productive life.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication March 18, 2008. Accepted for publication June 17, 2008.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |