Association between a genetic index for digital dermatitis resistance and the presence of digital dermatitis, heel horn erosion, and interdigital hyperplasia in Holstein cows

Digital dermatitis (DD) is a polybacterial disease endemic to most UK dairy farms. It poses a major financial and welfare threat and is characterized by high incidence and recurrence rates. We aimed to investigate the association between the UK EBV for resistance to digital dermatitis, the digital dermatitis index (DDI), and the frequency of DD, heel horn erosion (HHE), and interdigital hyperplasia (IH) in a population of Holstein dairy cows. We enrolled and genotyped 2,352 cows from 4 farms in a prospective cohort study. Foot lesion records were recorded by veterinary surgeons for each animal at 4 time points during a production cycle, starting at approximately 2 mo before calving and ending in late lactation. Importantly, these records were not used in the calculation of the DDI. Lesion records were matched to the animal's own DDI (n = 2,101) and their sire's DDI (n = 1,812). Digital dermatitis index values in our study population ranged from -1.41 to +1.2 and were transformed to represent distance from the mean expressed in SD. The relationship between the DDI and the presence of DD was investigated using a logistic regression model, with farm, parity, and a farm-parity interaction fitted as covariates. A multivariable logistic regression model was fitted to evaluate the relationship between HHE and DDI with farm fitted as a covariate. Finally, a univariable logistic regression model with DDI as explanatory variable was used to investigate the relationship between IH and DDI. The odds ratio of an animal being affected by DD was 0.69 for 1 SD increase in the animal's DDI (95% CI = 0.63-0.76). The odds of HHE and IH were 0.69 (95% CI = 0.62-0.76) and 0.58 (95% CI = 0.49-0.68) respectively for 1 SD increase in DDI. The adjusted probability of DD was 32% (95% CI = 27-36%) for cows with mean DDI value of 0, while it was 24% (95% CI = 20-29%) in cows with a DDI value of +1. Sire DDI breeding values were standardized in the same way and then binned into terciles creating an ordinal variable representing bulls of high, medium, and low genetic merit for DD resistance. The daughters of low genetic merit bulls were at 2.05 (95% CI = 1.60-2.64), 1.96 (95% CI = 1.53-2.50), and 2.85 (95% CI = 1.64-5.16) times greater odds of being affected by DD, HHE, and IH, respectively, compared with the daughters of high genetic merit bulls. The results of this study highlight the potential of digital dermatitis genetic indexes to aid herd management of DD, and suggest that breeding for resistance to DD, alongside environmental and management control practices, could reduce the prevalence of the disease.

A scoring system for grading corkscrew claws in dairy cattle

Corkscrew claw (CC) in dairy cattle is increasingly reported in dairy herds. CC is a progressive deformity of the claw capsule with uncertain aetiology and pathogenesis. Genetics and specific environmental factors are suspected of contributing to the development of this irreversible condition. CC has been found in lame cows; however, the cause and effect has not been established. To perform analysis of risk factors, treatment and pathogenesis, a definition of severity scores is called for. The aim of this study was to measure and analyse CC characteristics from photos of cows' feet to describe and evaluate a scoring system for CC. Width of the visible part of the axial wall, degree of contact between the toe and the floor and angle of the distal part of the abaxial wall as a proxy for the deviation of the abaxial wall was measured from 393 pictures of CC. Based on the measurements on the claws, the parameter "width of the axial wall" was chosen to define the scores. The parameter was divided into three intervals to define either mild CC 0.3-2.0 cm, moderate CC 2.1-3.5 cm or severe CC>3.5 cm and correlation between the parameters; level of contact between the toe and the floor and the angle of the distal abaxial wall was evaluated. There was a significant positive linear correlation between width of the axial wall and angle of the distal part of the abaxial wall (r=0.91), the wider the axial wall, the more the abaxial wall deviated in the distal part. As the width of the axial wall increased the toe increasingly lost contact with the floor, this association was significant for mild CC and moderate CC but not for severe CC. The Interobserver agreement of the CC Scoring system was tested by 30 claw trimmers each scoring 32 cadaver feet and by 2 trained observers on 28 photos of feet using Cohen´s weighted kappa and showed substantial to almost perfect agreement between untrained and trained observers, respectively.

Genetic selection for reduced health treatment costs in Holstein cows: implications from a long-term study

The objective of this study was to estimate genetic parameters of health treatment cost of Holstein cows from producer-recorded health treatments in 8 herds over an 8-yr period of time. Genetic parameters of health treatment cost were estimated in first (n = 2,214), second (n = 1,487) and third (n = 800) parities of US Holstein cows. The health treatments were uniformly defined and consistently recorded by 8 high-performance dairy herds in Minnesota from 2008 to 2015. A fixed treatment cost was assigned to 14 types of health treatments, and the cost included the mean veterinary expense obtained from the veterinary clinics that serviced the 8 herds, pharmaceuticals, and labor cost. The labor cost was $18/h, and the time incurred for each type of health treatment was determined from interviews with the herd owners. The 14 types of health treatment costs were partitioned into 5 categories: mastitis (including mastitis diagnostic test), reproduction (cystic ovary, retained placenta, and metritis), lameness (hoof treatments), metabolic (milk fever, displaced abomasum, ketosis, and digestive), and miscellaneous (respiratory, injury, and other). Health treatment cost for each cow was summed by category within lactation and also across categories within lactation. The estimates of heritability for health treatment cost were 0.13, 0.04, 0.10, 0.12, and 0.04 for the mastitis, reproduction, lameness, metabolic, and miscellaneous categories, respectively, in first parity. Genetic correlations between categories of health treatment cost in first parity were greatest for mastitis and reproduction (r = 0.85); however, phenotypic correlations between all categories were small (r < 0.16). Total health treatment cost had a large genetic correlation with somatic cell score (0.93) and 305-d milk production (0.44) in first parity; however, the genetic correlation (-0.60) between total health treatment cost and udder depth in first parity indicated a genetic relationship exists between shallow udders and less total health treatment cost. Total health treatment cost across categories had a heritability estimate of 0.25 in first parity, 0.16 in second parity, and 0.17 in third parity. Consequently, genetic selection for reduced health treatment cost should be possible by using producer-recorded health treatment records supplemented with treatment costs.

Genetic parameters of sole lesion recovery in Holstein cows

Sole hemorrhage and sole ulcers, referred to as sole lesions, are important causes of lameness in dairy cattle. The objective of this study was to estimate the genetic parameters of a novel trait reflecting how well cows recovered from sole lesions and the genetic correlation of this trait with overall susceptibility to sole lesions. A cohort of Holstein dairy cows was prospectively enrolled on 4 farms and assessed at 4 timepoints: before calving, immediately after calving, in early lactation, and in late lactation. At each timepoint, sole lesions were recorded at the claw level by veterinary surgeons and used to define 2 binary traits: (1) susceptibility to sole lesions-whether animals were affected with sole lesions at least once during the study or were unaffected at every assessment, and (2) sole lesion recovery-whether sole lesions healed between early and late lactation. Animals were genotyped and pedigree details extracted from the national database. Analyses were conducted with BLUPF90 software in a single-step framework; genetic parameters were estimated from animal threshold models using Gibbs sampling. The genetic correlation between both traits was approximated as the correlation between genomic estimated breeding values, adjusting for their reliabilities. A total of 2,025 animals were used to estimate the genetic parameters of sole lesion susceptibility; 44% of animals recorded a sole lesion at least once during the study period. The heritability of sole lesion susceptibility, on the liability scale, was 0.25 (95% highest density interval = 0.16-0.34). A total of 498 animals were used to estimate the genetic parameters of sole lesion recovery; 71% of animals had recovered between the early and late lactation assessments. The heritability of sole lesion recovery, on the liability scale, was 0.27 (95% highest density interval = 0.02-0.52). The approximate genetic correlation between each trait was -0.11 (95% confidence interval = -0.20 to -0.02). Our results indicate that recovery from sole lesions is heritable. If this finding is corroborated in further studies, it may be possible to use selective breeding to reduce the frequency of chronically lame cows. As sole lesion recovery appears to be weakly genetically related to sole lesion susceptibility, successful genetic improvement of sole lesion recovery would benefit from selection on this trait directly.

Evaluating the long-term conformation and hoof growth effects of starting hoof trimming at 5 months of age in New Zealand dairy goats

Hoof overgrowth is associated with poor conformation, an altered weight-bearing surface, and a reduction in the hoof's anatomic and functional integrity. As a result of housing systems that promote hoof overgrowth, hoof trimming is considered a priority in dairy goats. However, there are few data on the effects of the timing of first trimming on hoof conformation and growth rate. The aims of this study were (1) to evaluate the long-term effects of 2 different hoof trimming start times and (2) to investigate the pattern of hoof growth across the first 2 yr of life. Eighty 5-mo-old female Saanen-cross commercially housed dairy goats were allocated randomly to 1 of 2 treatments: (1) early trimmed (trimming beginning at 5 mo old; hooves were trimmed every 4 mo thereafter) and (2) late trimmed (trimming beginning at 13 mo old; hooves were trimmed every 4 mo thereafter). Using a combination of photographs and radiographs, hoof conformation, joint positions, and hoof wall length were assessed before the 13- and 25-mo trimming events. Hoof growth was assessed every 12 wk using caliper measurements. Overall, starting hoof trimming earlier had minor and inconsistent effects. However, detrimental changes in conformation and joint positions occurred between trimming events, particularly in the hind hooves, regardless of trimming treatment. At both assessments, there was a high percentage of overgrown toes and dipped heels, with the hind hooves being more affected compared with the front (overgrown toes at 13 mo, 97.1 vs. 79.1 ± 5.2%; overgrown toes at 25 mo, 91.7 vs. 56.3 ± 6.7%; dipped heels at 13 mo, 98.5 vs. 19.3 ± 5.0%; dipped heels at 25 mo, 88.3 vs. 4.9 ± 4.8%). In addition, at both assessments, the distal interphalangeal joint angle was greater in the hind hooves compared with the front (13 mo, 79.5 vs. 65.2 ± 1.7°; 25 mo, 79.0 vs. 66.7 ± 0.9°), whereas heel angles were less in the hind hooves compared with the front (13 mo, 41.8 vs. 57.1 ± 1.5°; 25 mo, 44.9 vs. 55.9 ± 1.1°). On average, the front hooves grew 4.39 mm/mo and the hind hooves grew 4.20 mm/mo. Early trimming did not have consistent effects on hoof growth rate. Importantly, our results suggest that trimming every 4 mo is not sufficient to prevent hoof overgrowth, the development of poor conformation, and detrimental changes in joint positions, particularly in the hind hooves. Furthermore, the detrimental changes may have masked any long-term treatment effects. Therefore, trimming frequency and age of first trimming should be considered when devising hoof care protocols for dairy goats housed in environments that do not offer opportunities for natural hoof wear.

Benchmarking Based on Regularly Recorded Claw Health Data of Austrian Dairy Cattle for Implementation in the Cattle Data Network (RDV)

While benchmarking is already used for the assessment of performance gaps in cattle herd management and welfare concerns, its application to quantifying claw health performance is relatively new. The goal here was to establish a benchmarking system for claw health in Austrian dairy cattle. We used electronically registered claw health data of cows from 512 dairy herds documented by professional hoof trimmers, culling data from the same herds, and locomotion scores taken at regular milk performance testings in 99 dairy herds during 2020. Mean, median and the 10th, 25th, 75th, and 90th percentiles of the incidences of risk of lameness, 13 common claw lesions, and the annual culling risk directly related to claw and limb disorders were used as key performance indicators. Only validated data sets were used and participating trimmers and locomotion scorers had to pass interobserver reliability tests with weighted Cohen's kappa values ≥ 0.61 indicating substantial interobserver agreement. This claw health benchmarking system is intended to be used henceforth in the transnational cattle data network (RDV) by all participating farmers and is also available for veterinarians and consultants, with the agreement of respective farmers.

Symposium review: Multiple-trait single-step genomic evaluation for hoof health

Hoof lesions represent an important issue in modern dairy herds, with reported prevalence in different countries ranging from 40 to 70%. This high prevalence of hoof lesions has both economic and social consequences, resulting in increased labor expenses and decreasing animal production, longevity, reproduction, health, and welfare. Therefore, a key goal of dairy herds is to reduce the incidence of hoof lesions, which can be achieved both by improving management practices and through genetic selection. The Canadian dairy industry has recently released a hoof health sub-index. This national genetic evaluation program for hoof health was achieved by creating a centralized data collection system that routinely transfers data recorded by hoof trimmers into a coherent and sustainable national database. The 8 most prevalent lesions (digital dermatitis, interdigital dermatitis, interdigital hyperplasia, heel horn erosion, sole hemorrhage, sole ulcer, toe ulcer, and white line lesion) in Canada are analyzed with a multiple-trait model using a single-step genomic BLUP method. Estimated genomic breeding values for each lesion are combined into a sub-index according to their economic value and prevalence. In addition, data recorded within this system were used to create an interactive management report for dairy producers by Canadian DHI, including the prevalence of lesions on farm, their trends over time, and benchmarks with provincial and national averages.

Two-stage genome-wide association study for the identification of causal variants underlying hoof disorders in cattle

Feet and legs disorders influence dairy cattle breeding by their effect on animal welfare, economic losses due to lower production and fertility, costs of treatment, and problems with herd management. In our study, we estimated heritabilities and performed a 2-step GWAS for 3 traits describing hoof health: hoof health status defined by a veterinarian (HSV), hoof health status defined by a claw trimmer (HSC), and the total number of hoof disorders (NHD), scored in 1,998 Fleckvieh and 979 Braunvieh cows. The individuals were genotyped with a high-density (HD) panel consisting of 76,934 SNP. For significant genomic regions, the SNP information was enhanced by SNP imputed from the whole-genome sequence of Fleckvieh and Braunvieh bulls from the 1000 Bulls Genome project. The heritabilities were estimated to be 0.035 for HSV, 0.249 for HSC, and 0.279 for NHD. Based on the first-stage GWAS with SNP from the HD panel, 7 significant genomic regions on 6 chromosomes were defined: (1) 120 SNP spanning 15,522 bp on BTA1, including the TOPBP1 gene; (2) 4,139 SNP spanning 1,426,046 bp on BTA7, including the RIOK2 and RGMB genes; (3) 167 SNP spanning 167,352 bp on BTA13, including the C13H20orf194 gene; (4) 2 regions on BTA14, one harboring 1,071 SNP spanning 380,024 bp, including RRM2B and NCALD, and the other comprising 632 SNP spanning 385,111 bp, including STK3; (5) 328 SNP on BTA15, spanning 235,567 bp between FAM168A and PLEKHB1; and (6) 1,549 SNP on BTA22, spanning 596,101 bp in the neighborhood of PTPRG. Then, we conducted a second-stage GWAS based on SNP from whole-genome sequences within the significant regions obtained in the first stage of the analysis. For HSV, the highest additive effect was estimated for 23 SNP located within a region on BTA15, close to FAM168A, corresponding to a predicted gene sequence. For HSC, the highest additive effect was attributed to 44 SNP located within a region of BTA22 corresponding to 4 predicted gene sequences, with rs135082893 within a sequence encoding a microRNA. Another potential causal mutation for HSC was rs134142607 on BTA13, within the exon of C13H20orf194. For NHD, 33 SNP with the highest estimated effect were located on BTA7 within a region of a predicted gene positioned between RIOK2 and RGMB. On BTA14, all significant SNP were located in introns of STK3, which is responsible for the "abnormal gait" phenotype in mice.

Distribution analysis of single nucleotide polymorphisms related to feet and legs and their association with lameness in Vrindavani cattle

Lameness is one of the most serious economic problems in dairy animals. The present study was undertaken to identify the SNPs associated with lameness in Vrindavani, a synthetic crossbred cattle developed in India. A total of 187 cows were screened for lameness at Cattle and Buffalo Farm, Indian Veterinary Research Institute, Izatnagar, India and 16.04% prevalence was noticed. The effect of age, body weight, parity, and peak yield was not significant on lameness score. Heritability estimate of lameness score was low (0.07 ± 0.01). A total of 90 Vrindavani cows were genotyped by PCR-RFLP for 10 SNPs that were reported to have association with feet and legs in cattle. Five SNPs displayed significant departure from Hardy Weinberg equilibrium. Varying levels of heterozygosity, polymorphic information content, and allelic diversity were noticed at different SNPs. Odds ratio and least squares analysis revealed a significant association (P < 0.05) of three SNPs with lameness in the herd. Cows with AG, TT, and GG genotypes respectively at rs41632254, rs41603160, and rs41636945 loci were less prone to lameness. These genetic markers would certainly aid in development of an early detection system for lameness.

Investigating the genetic background of bovine digital dermatitis using improved definitions of clinical status

Bovine digital dermatitis (DD) is an increasing claw health problem in all cattle production systems worldwide. The objective of this study was to evaluate the use of an improved scoring of the clinical status for DD via M-scores accounting for the dynamics of the disease; that is, the transitions from one stage to another. The newly defined traits were then subjected to a genetic analysis to determine the genetic background for susceptibility to DD. Data consisted of 6,444 clinical observations from 729 Holstein heifers in a commercial dairy herd, collected applying the M-score system. The M-score system is a classification scheme for stages of DD that allows a macroscopic scoring based on clinical inspections of the bovine foot, thus it describes the stages of lesion development. The M-scores were used to define new DD trait definitions with different complexities. Linear mixed models and logistic models were used to identify fixed environmental effects and to estimate variance components. In total, 68% of all observations showed no DD status, whereas 11% were scored as infectious for and affected by DD, and 21% of all observations exhibited an affected but noninfectious status. For all traits, the probability of occurrence and clinical status were associated with age at observation and period of observation. Risk of becoming infected increased with age, and month of observation significantly affected all traits. Identification of the optimal month concerning DD herd status was consistent for all trait definitions; the last month of the trial was identified. In contrast, months exhibiting the highest least squares means of transformed scores differed depending on trait definition. In this respect, traits that can distinguish between healthy, infectious, and noninfectious stages of DD can account for the infectious potential of the herd and can serve as an alert tool. Estimates of heritabilities of traits studied ranged between 0.19 (±0.11) and 0.52 (±0.17), revealing a tendency for higher values for more complex trait definitions. In terms of genetic selection, all trait definitions identified the best (i.e., most resistant) animals, but only the new trait definitions were able to distinguish between animals with average and high predispositions for DD. Considering repeated measurements resulted in heritability estimates ranging between 0.13 (±0.05) and 0.29 (±0.10).

Foot and leg conformation traits have a small effect on genomic predictions of claw disorders in Norwegian Red cows

The aim of this study was to evaluate whether the predictive correlation of genomic breeding values (GEBV) for claw disorders increased by including genetically correlated traits as additional information in the analyses. Predictive correlations of GEBV for claw disorders were calculated based on claw disorders only and by analyzing claw disorders together with genetically correlated foot and leg conformation traits. The claw disorders analyzed were corkscrew claw (CSC); infectious claw disorder, including dermatitis, heel horn erosion, and interdigital phlegmon; and laminitis-related claw disorder, including sole ulcer, white line disorder, and hemorrhage of sole and white line. The foot and leg conformation traits included were hoof quality, foot angle, rear leg rear view new, and rear leg rear view old. The data consisted of 183,728 daughters with claw health records and 421,319 daughters with foot and leg conformation scores. A 25K/54K single nucleotide polymorphism (SNP) data set containing 48,249 SNP was available for the analyses. The number of genotyped sires with daughter information in the analyses was 1,093 including claw disorders and 3,111 including claw disorders and foot and leg conformation traits. Predictive correlations of GEBV for CSC, infectious claw disorder, and laminitis-related claw disorder were calculated from a 10-fold cross-validation and from an additional validation set including the youngest sires. Only sires having daughters with claw health records were in the validation sets, thus increasing the reference population when adding foot and leg conformation traits. The results showed marginal improvement in the predictive correlation of GEBV for CSC when including hoof quality and foot angle, both in 10-fold cross-validation (from 0.35 to 0.37) and in the validation including the youngest sires (from 0.38 to 0.49). For infectious claw disorder and laminitis-related claw disorder, including foot and leg conformation traits had no effect on the predictive correlation of GEBV. Claw disorders are novel traits with a limited amount of historical data and, therefore, a small reference population. Increasing the reference population by including sires with daughter information on foot and leg conformation traits had small effect on the predictive correlation of GEBV. However, the small increase in predictive correlation of GEBV for CSC shows a possible gain when including moderate to high genetically correlated traits.