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.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits

BACKGROUND

Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected.

RESULTS

All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count.

CONCLUSIONS

Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Single-step genome-wide association analyses of claw horn lesions in Holstein cattle using linear and threshold models

BACKGROUND

Lameness in dairy cattle is primarily caused by foot lesions including the claw horn lesions (CHL) of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). This study investigated the genetic architecture of the three CHL based on detailed animal phenotypes of CHL susceptibility and severity. Estimation of genetic parameters and breeding values, single-step genome-wide association analyses, and functional enrichment analyses were performed.

RESULTS

The studied traits were under genetic control with a low to moderate heritability. Heritability estimates of SH and SU susceptibility on the liability scale were 0.29 and 0.35, respectively. Heritability of SH and SU severity were 0.12 and 0.07, respectively. Heritability of WL was relatively lower, indicating stronger environmental influence on the presence and development of WL than the other two CHL. Genetic correlations between SH and SU were high (0.98 for lesion susceptibility and 0.59 for lesion severity), whereas genetic correlations of SH and SU with WL also tended to be positive. Candidate quantitative trait loci (QTL) were identified for all CHL, including some on Bos taurus chromosome (BTA) 3 and 18 with potential pleiotropic effects associated with multiple foot lesion traits. A genomic window of 0.65 Mb on BTA3 explained 0.41, 0.50, 0.38, and 0.49% of the genetic variance for SH susceptibility, SH severity, WL susceptibility, and WL severity, respectively. Another window on BTA18 explained 0.66, 0.41, and 0.70% of the genetic variance for SH susceptibility, SU susceptibility, and SU severity, respectively. The candidate genomic regions associated with CHL harbour annotated genes that are linked to immune system function and inflammation responses, lipid metabolism, calcium ion activities, and neuronal excitability.

CONCLUSIONS

The studied CHL are complex traits with a polygenic mode of inheritance. Most traits exhibited genetic variation suggesting that animal resistance to CHL can be improved with breeding. The CHL traits were positively correlated, which will facilitate genetic improvement for resistance to CHL as a whole. Candidate genomic regions associated with lesion susceptibility and severity of SH, SU, and WL provide insights into a global profile of the genetic background underlying CHL and inform genetic improvement programmes aiming at enhancing foot health in dairy cattle.

Genomic analyses of claw disorders in Holstein cows: Genetic parameters, trait associations, and genome-wide associations considering interactions of SNP and heat stress

The aim of the present study was an in-depth genomic analysis to understand the genomic mechanisms of the 3 claw disorders dermatitis digitalis (DD), interdigital hyperplasia (HYP), and sole ulcer (SU). In this regard, we estimated genetic parameters based on genomic relationship matrices, performed genome-wide association studies, annotated potential candidate genes, and inferred genetic associations with breeding goal traits considering the most important chromosomal segments. As a further novelty of this study, we inferred possible SNP × heat stress interactions for claw disorders. The study consisted of 17,264 first-lactation Holstein Friesian cows kept in 50 large-scale contract herds. The disease prevalence was 15.96, 2.36, and 8.20% for DD, HYP, and SU, respectively. The remaining breeding goal traits consisted of type traits of the feet and leg composite, female fertility, health traits, and 305-d production traits. The final genotype data set included 44,474 SNPs from the 17,264 genotyped cows. Heritabilities for DD, HYP, and SU were estimated in linear and threshold models considering the genomic relationship matrix (G matrix). Genetic correlations with breeding goal traits based on G were estimated in a series of bivariate linear models, which were verified via SNP effect correlations for specific chromosome segments (i.e., segments harboring potential candidate genes for DD, HYP, and SU). Genome-wide association studies were performed for all traits in a case-control design by applying a single SNP linear mixed model. Furthermore, for DD, HYP, and SU, we modeled SNP × heat stress interactions in genome-wide association studies. Single nucleotide polymorphism-based heritabilities were 0.04 and 0.08 for DD, 0.03 and 0.10 for SU, and 0.03 and 0.23 for HYP from linear and threshold models, respectively. The genetic correlations between DD, HYP, and SU with conformation traits from the feet and leg composite were positive throughout, indicating the value of indirect selection on conformation traits to improve claw health. Genetic correlations between DD, SU, and HYP with other breeding goal traits indicated impaired female fertility, impaired udder health status, and productivity decline of diseased cows. Genetic correlations among DD, SU, and HYP were moderate to large, indicating that different claw disorders have similar genetic mechanisms. Nevertheless, we identified disease-specific potential candidate genes, and genetic associations based on the surrounding SNPs partly differed from the genetic correlations. Especially for candidate genes contributing to 2 traits simultaneously, correlations based on SNP effects from the respective chromosome segment were close to 1 or to -1. In this regard, we annotated the candidate genes KRT33A and KRT33B for HYP and DD, KIF27 for HYP and calving to first insemination, and MAN1A1 for SU and the production traits. For SNP × heat stress interactions, we identified significant SNPs on BTA 2, 4, 5, 7, 8, 9, 13, 22, 25, and 28, and we annotated the potential candidate genes FSIP2, CLCN1, ADGRV1, DOP1A, THBD, and RHOBTB1. Results indicate gene-specific mechanisms of the claw disorders only in specific environments.

Association between a genetic index for lameness resistance and the incidence of claw horn lesions in Holstein cows

BACKGROUND

This study aimed to determine the association between the lameness advantage genetic index and four outcomes: sole haemorrhage (SH), sole ulcers (SU), white line lesions (WL), and lameness during mobility scoring.

METHODS

We enrolled 2352 Holstein cows from four predominantly housed dairy herds in the UK. Cows were mobility scored and foot lesions recorded at four time points from before calving to late lactation. Cows were genotyped and genetic indexes were assigned to each cow following national genetic evaluations. Lameness records and genetic indexes were matched for 2107 cows. Four separate multivariable logistic regression models, which included farm and parity as covariables, were used to quantify the association between the lameness advantage index and whether animals were affected by SH, SU, WL, or lameness.

RESULTS

The odds ratios (95% confidence intervals) for one-point increases in the lameness advantage index were 0.79 (0.72-0.86), 0.68 (0.59-0.78), 0.94 (0.84-1.04), and 0.82 (0.74-0.91) for SH, SU, WL, and lameness, respectively. The same trends were present when the sire's lameness advantage index was evaluated in place of the animal's own, although the strength of this association was generally weaker.

CONCLUSION

The lameness advantage index is associated with SH, SU, and lameness, therefore selection on the lameness advantage index could be considered in herds aiming to reduce lameness. Where genomic testing of heifers is not conducted, sire lameness advantage index may still be effective to reduce SH and SU incidence.

Pleiotropic Loci Associated With Foot Disorders and Common Periparturient Diseases in Holstein Cattle

Lameness is an animal welfare issue that incurs substantial financial and environmental costs. This condition is commonly caused by digital dermatitis (DD), sole ulcers (SU), and white line disease (WLD). Susceptibility to these three foot disorders is due in part to genetics, indicating that genomic selection against these foot lesions can be used to reduce lameness prevalence. It is unclear whether selection against foot lesions will lead to increased susceptibility to other common diseases such as mastitis and metritis. Thus, the aim of this study was to determine the genetic correlation between causes of lameness and other common health disorders to identify loci contributing to the correlation. Genetic correlation estimates between SU and DD and between SU and WLD were significantly different from zero (p < 0.05), whereas estimates between DD and mastitis, DD and milk fever, and SU and metritis were suggestive (p < 0.1). All five of these genetic correlation estimates were positive. Two-trait genome-wide association studies (GWAS) for each of these five pairs of traits revealed common regions of association on BTA1 and BTA8 for pairs that included DD or SU as one of the traits, respectively. Other regions of association were unique to the pair of traits and not observed in GWAS for other pairs of traits. The positive genetic correlation estimates between foot disorders and other health disorders imply that selection against foot disorders may also decrease susceptibility to other health disorders. Linkage disequilibrium blocks defined around significant and suggestive SNPs from the two-trait GWAS included genes and QTL that were functionally relevant, supporting that these regions included pleiotropic loci.

Simulating the mechanics behind sub-optimal mobility and the associated economic losses in dairy production

Hoof disorders and sub-optimal mobility (SOM) are economically important health issues in dairy farming. Although the dynamics of hoof disorders have an important effect on cow mobility, they have not been considered in previous simulation models that estimate the economic loss of SOM. Furthermore, these models do not consider the varying severities of SOM. The objective of this study was to develop a novel bio-economic simulation model to simulate the dynamics of 8 hoof disorders: digital dermatitis (DD), interdigital hyperplasia (HYP), interdigital dermatitis/heel-horn erosion (IDHE), interdigital phlegmon (IP), overgrown hoof (OH), sole haemorrhage (SH), sole ulcer (SU) and white-line disease (WLD), their role in SOM, and estimate the economic loss of SOM in a herd of 125 dairy cows. A Reed-Frost model was used for DD and a Greenwood model for the other 7 hoof disorders. Economic analysis was conducted per mobility score according to a 5-point mobility scoring method (1 = perfect mobility; 5 = severely impaired mobility) by comparing a scenario with SOM and one without SOM. Parameters used in the model were based on literature and expert opinion and deemed credible during model validation rounds. Results showed that the mean cumulative incidence for maximum mobility scores 2-5 SOM episodes were respectively 34, 16, 7 and <1 episodes per 100 cows per pasture period and 39, 19, 8, <1 episodes per 100 cows per housing period. The mean total annual economic loss due to SOM resulting from the hoof disorders under study was €15,342: €122 per cow per year. The economic analysis uncovered direct economic losses that could be directly linked to SOM episodes and indirect economic losses that could not be directly linked to SOM episodes but arose due to the presence of SOM. The mean total annual direct economic loss for maximum mobility score 2-5 SOM episodes was €1129, €3098, €4354 and €480, respectively. The mean total annual indirect economic loss varied considerably between the 5th and 95th percentiles: €-6174 and €19,499, and had a mean of €6281. This loss was composed of additional indirect culling due to SOM (∼65%) and changes in the overall herd milk production (∼35%) because of additional younger replacement heifers entering the herd due to increased culling rates. The bio-economic model presented novel results with respect to indirect economic losses arising due to SOM. The results can be used to stimulate farmer awareness and promote better SOM management.

Genome-Wide Association Studies Reveal Susceptibility Loci for Noninfectious Claw Lesions in Holstein Dairy Cattle

Sole ulcers (SUs) and white line disease (WLD) are two common noninfectious claw lesions (NICL) that arise due to a compromised horn production and are frequent causes of lameness in dairy cattle, imposing welfare and profitability concerns. Low to moderate heritability estimates of SU and WLD susceptibility indicate that genetic selection could reduce their prevalence. To identify the susceptibility loci for SU, WLD, SU and/or WLD, and any type of noninfectious claw lesion, genome-wide association studies (GWAS) were performed using generalized linear mixed model (GLMM) regression, chunk-based association testing (CBAT), and a random forest (RF) approach. Cows from five commercial dairies in California were classified as controls having no lameness records and ≥6 years old (n = 102) or cases having SU (n = 152), WLD (n = 117), SU and/or WLD (SU + WLD, n = 198), or any type of noninfectious claw lesion (n = 217). The top single nucleotide polymorphisms (SNPs) were defined as those passing the Bonferroni-corrected suggestive and significance thresholds in the GLMM analysis or those that a validated RF model considered important. Effects of the top SNPs were quantified using Bayesian estimation. Linkage disequilibrium (LD) blocks defined by the top SNPs were explored for candidate genes and previously identified, functionally relevant quantitative trait loci. The GLMM and CBAT approaches revealed the same regions of association on BTA8 for SU and BTA13 common to WLD, SU + WLD, and NICL. These SNPs had effects significantly different from zero, and the LD blocks they defined explained a significant amount of phenotypic variance for each dataset (6.1-8.1%, p < 0.05), indicating the small but notable contribution of these regions to susceptibility. These regions contained candidate genes involved in wound healing, skin lesions, bone growth and mineralization, adipose tissue, and keratinization. The LD block defined by the most significant SNP on BTA8 for SU included a SNP previously associated with SU. The RF models were overfitted, indicating that the SNP effects were very small, thereby preventing meaningful interpretation of SNPs and any downstream analyses. These findings suggested that variants associated with various physiological systems may contribute to susceptibility for NICL, demonstrating the complexity of genetic predisposition.

Laminar Inflammation Responses in the Oligofructose Overload Induced Model of Bovine Laminitis

Bovine laminitis causes substantial economic losses and animal welfare problems in dairy farms worldwide. Previously published studies have reported that the inflammatory response plays a central role in the pathogenesis of the disease. To our knowledge, inflammation associated with bovine laminitis induced by high levels of exposure to oligofructose (OF) has not been reported and characterized. In fact, the disease manifestations in this model closely approximate those of clinical laminitis. The objective of this study was to characterize the inflammatory response in OF-induced bovine laminitis. A total of 12 Chinese Holstein dairy heifers were utilized in this study. The heifers were randomly divided into two groups, treatment (n = 6) and control (n = 6). The treatment group heifers were administered OF solutions via a stomach tube (dose: 17 g/kg of body weight). Upon development of a lameness score of 2 with consecutive positive reactions in the same claw, they would be humanely euthanized. Control heifers were administered deionized water (dose: 2 L/100 kg of body weight) and humanely euthanized at 72 h. Real-time quantitative PCR (qPCR) assays were performed to determine the messenger RNA (mRNA) concentrations of inflammatory mediators in the lamellae. Concentrations of interleukin (IL)-1β, IL-6, IL-8, C-X-C motif chemokine ligand-1 (CXCL-1), macrophage cationic peptide-2 (MCP-2), E-selectin, intercellular adhesion molecule-1 (ICAM-1), cyclooxygenase-2 (COX-2), inducible nitric oxide synthase-1 (iNOS-1), and plasminogen activator inhibitor-1 (PAI-1) were significantly increased (P < 0.05) in the treatment group. No significant difference was found for tumor necrosis factor alpha (TNF-α), IL-10, CXCL-6, and MCP-1. These results demonstrated and characterized the laminar inflammatory response leading to the pathogenesis of bovine laminitis at the early stages.

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.

Genetic parameters and genomic breeding values for digital dermatitis in Holstein Friesian dairy cattle: host susceptibility, infectivity and the basic reproduction ratio

Background

For infectious diseases, the probability that an animal gets infected depends on its own susceptibility, and on the number of infectious herd mates and their infectivity. Together with the duration of the infectious period, susceptibility and infectivity determine the basic reproduction ratio of the disease (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0). \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 is the average number of secondary cases caused by a typical infectious individual in an otherwise uninfected population. An infectious disease dies out when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 is lower than 1. Thus, breeding strategies that aim at reducing disease prevalence should focus on reducing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0, preferably to a value lower than 1. In animal breeding, however, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 has received little attention. Here, we estimate the additive genetic variance in host susceptibility, host infectivity, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 for the endemic claw disease digital dermatitis (DD) in Holstein Friesian dairy cattle, and estimate genomic breeding values (GEBV) for these traits. We recorded DD disease status of both hind claws of 1513 cows from 12 Dutch dairy farms, every 2 weeks, 11 times. The genotype data consisted of 75,904 single nucleotide polymorphisms (SNPs) for 1401 of the cows. We modelled the probability that a cow got infected between recordings, and compared four generalized linear mixed models. All models included a genetic effect for susceptibility; Models 2 and 4 also included a genetic effect for infectivity, while Models 1 and 2 included a farm*period interaction. We corrected for variation in exposure to infectious herd mates via an offset.

Results

GEBV for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 from the model that included genetic effects for susceptibility only had an accuracy of ~ 0.39 based on cross-validation between farms, which is very high given the limited amount of data and the complexity of the trait. Models with a genetic effect for infectivity showed a larger bias, but also a slightly higher accuracy of GEBV. Additive genetic standard deviation for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 was large, i.e. ~ 1.17, while the mean \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 was 2.36.

Conclusions

GEBV for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 showed substantial variation. The mean \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ R_{0} $$\end{document}R0 was only about one genetic standard deviation greater than 1. These results suggest that lowering DD prevalence by selective breeding is promising.

Interdigital Hyperplasia in Holstein Cattle Is Associated With a Missense Mutation in the Signal Peptide Region of the Tyrosine-Protein Kinase Transmembrane Receptor Gene

Bovine interdigital hyperplasia (IH) is a typical disease of the foot with varying prevalence depending on age, breed, and environmental factors resulting in different degrees of lameness. In studies based on assessments of claw health status at time of hoof trimming and applying genetic-statistical models to analyze this data, IH consistently exhibits high estimates of heritability in the range of 0.30–0.40. Although some studies have identified chromosomal regions that could possibly harbor causative genes, a clear identification of molecular causes for IH is lacking. While analyzing the large database of claw health status as documented at time of hoof trimming, we identified one herd with extreme prevalence of IH of > 50% of affected Holstein dairy cows. This herd subsequently was chosen as the object of a detailed study. A total of n = 91 cows was assessed and revealed a prevalence of 59.3% and 38.5% for IH cases, documented as “one-sided” or “two-sided”, respectively. Cows were genotyped using the BovineSNP50 BeadChip. A genome wide association study revealed two significantly associated chromosomal positions (-log₁₀P = 5.57) on bovine chromosome 8 (BTA8) located in intron 5 and downstream of the receptor tyrosine kinase-like orphan receptor 2 (ROR2) gene. As ROR2 plays a key role in ossification of the distal limbs and is associated with brachydactylies in humans, it was a reasonable candidate for IH. A comparative sequencing of the ROR2 gene between cases and controls revealed two missense variants in exon 1 (NC_037335.1:g.85,905,534T > A, ARS-UCD1.2) and exon 9 (NC_037335.1:g.86,140,379A > G, ARS-UCD1.2), respectively. Genotyping of both variants in the cohort of 91 cattle showed that the exon 1 variant (rs377953295) remained significantly associated with IH (p < 0.0001) as a risk factor of the disease. This variant resulted in an amino acid exchange (ENSBTAP00000053765.2:p.Trp9Arg) in the N-terminal region of the ROR2 signal peptide which is necessary for proper topology of the polypeptide during translocation. Quantification of ROR2 mRNA and ROR2 protein showed that the variant resulted in a significant suppression of ROR2 expression in homozygous affected compared to wild type and carrier cows.

Quantitative Trait Loci Mapping for Lameness Associated Phenotypes in Holstein-Friesian Dairy Cattle

Lameness represents a significant challenge for the dairy cattle industry, resulting in economic losses and reduced animal health and welfare. The existence of underlying genomic variation for lameness associated traits has the potential to improve selection strategies by using genomic markers. Therefore, the aim of this study was to identify genomic regions and potential candidate genes associated with lameness traits. Lameness related lesions and digital cushion thickness were studied using records collected by our research team, farm records, and a combination of both. Genome-wide analyses were performed to identify significant genomic effects, and a combination of single SNP association analysis and regional heritability mapping was used to identify associated genomic regions. Significant genomic effects were identified for several lameness related traits: Two genomic regions were identified on chromosome 3 associated with digital dermatitis and interdigital hyperplasia, one genomic region on chromosome 23 associated with interdigital hyperplasia, and one genomic region on chromosome 2 associated with sole haemorrhage. Candidate genes in those regions are mainly related to immune response and fibroblast proliferation. Quantitative trait loci (QTL) identified in this study could enlighten the understanding of lameness pathogenesis, providing an opportunity to improve health and welfare in dairy cattle with the addition of these regions into selection programs.

Udder and Claw-Related Health Traits in Selection of Czech Holstein Cows

The aim of the study was to investigate the impact of the implementation of new health traits into the breeding objective and selection criteria for Czech Holstein cows on the genetic selection response in the breeding objective traits. Incidence of overall claw diseases was included into the current breeding objective for cows (11 traits together). Three traits that indicated claw health (incidence of claw diseases overall and infectious and non-infectious claw diseases) and incidence of clinical mastitis were successively added to the current selection criteria in a cow selection index (a maximum of 19 traits). Selection responses in the breeding objective traits were estimated by applying the general principles of the selection index theory. The required genetic variances for the new traits, the economic weights for all breeding objective traits and the genetic correlations among the selection index traits were estimated within this study. The marginal economic weights, which were calculated for two-year-old cows by applying a bioeconomic model with implemented gene flow, were -193 and -168 € per case for clinical mastitis and overall claw disease incidence, respectively. Using the comprehensive selection index with 19 traits, the reduction in the incidence of both udder and claw diseases was calculated to be 0.004 cases per cow per year. At the same time, a more favourable genetic selection response was obtained for other functional traits, e.g., +0.020% for cow conception rate and +0.010 years for productive life of cow (which represented the profit of 67 € and 367 € per herd and per year, respectively) when compared to the current index. Based on this study, a direct selection of cows for claw and udder health is nowadays recommended to improve the health status of herds and economics in production systems.

A genome-wide association study for susceptibility and infectivity of Holstein Friesian dairy cattle to digital dermatitis

Selection and breeding can be used to fight transmission of infectious diseases in livestock. The prevalence in a population depends on the susceptibility and infectivity of the animals. Knowledge on the genetic background of those traits would facilitate efficient selection for lower disease prevalence. We investigated the genetic background of host susceptibility and infectivity for digital dermatitis (DD), an endemic infectious claw disease in dairy cattle, with a genome-wide association study (GWAS), using either a simple linear mixed model or a generalized linear mixed model based on epidemiological theory. In total, 1,513 Holstein-Friesian cows of 12 Dutch dairy farms were scored for DD infection status and class (M0 to M4.1) every 2 wk for 11 times; 1,401 of these cows were genotyped with a 75k SNP chip. We performed a GWAS with a linear mixed model on 10 host disease status traits, and with a generalized linear mixed model with a complementary log-log link function (GLMM) on the probability that a cow would get infected between 2 scorings. With the GLMM, we fitted SNP effects for host susceptibility and host infectivity, while taking the variation in exposure of the susceptible cow to infectious herd mates into account. With the linear model we detected 4 suggestive SNP (false discovery rate < 0.20), 2 for the fraction of observations a cow had an active lesion on chromosomes 1 and 14, one for the fraction of observations a cow had an M2 lesion on at least one claw on chromosome 1 (the same SNP as for the fraction of observations with an active lesion), and one for the fraction of observations a cow had an M4.1 lesion on at least one claw on chromosome 10. Heritability estimates ranged from 0.09 to 0.37. With the GLMM we did not detect significant nor suggestive SNP. The SNP effects on disease status analyzed with the linear model had a correlation coefficient of only 0.70 with SNP effects on susceptibility of the GLMM, indicating that both models capture partly different effects. Because the GLMM better accounts for the epidemiological mechanisms determining individual disease status and for the distribution of the y-variable, results of the GLMM may be more reliable, despite the absence of suggestive associations. We expect that with an extended GLMM that better accounts for the full genetic variation in infectivity via the environment, the accuracy of SNP effects may increase.

A 100-Year Review: Identification and genetic selection of economically important traits in dairy cattle

Over the past 100 yr, the range of traits considered for genetic selection in dairy cattle populations has progressed to meet the demands of both industry and society. At the turn of the 20th century, dairy farmers were interested in increasing milk production; however, a systematic strategy for selection was not available. Organized milk performance recording took shape, followed quickly by conformation scoring. Methodological advances in both genetic theory and statistics around the middle of the century, together with technological innovations in computing, paved the way for powerful multitrait analyses. As more sophisticated analytical techniques for traits were developed and incorporated into selection programs, production began to increase rapidly, and the wheels of genetic progress began to turn. By the end of the century, the focus of selection had moved away from being purely production oriented toward a more balanced breeding goal. This shift occurred partly due to increasing health and fertility issues and partly due to societal pressure and welfare concerns. Traits encompassing longevity, fertility, calving, health, and workability have now been integrated into selection indices. Current research focuses on fitness, health, welfare, milk quality, and environmental sustainability, underlying the concentrated emphasis on a more comprehensive breeding goal. In the future, on-farm sensors, data loggers, precision measurement techniques, and other technological aids will provide even more data for use in selection, and the difficulty will lie not in measuring phenotypes but rather in choosing which traits to select for.

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).

Genetic relationships between claw health traits of dairy cows in different parities, lactation stages, and herds with different claw disorder frequencies

Claw disorders affect cow welfare and profitability of farms. Not only claw disorders but also the need for trimming, known as trimming status, has been shown to be heritable. Limited knowledge is available on whether claw health traits (claw disorders and trimming status) are genetically the same trait in different parities, lactation stages, or in herds with low or high frequency of claw disorders. The aim of the current study was to estimate frequencies, heritabilities, and genetic correlations of claw health traits measured in different parities (first vs. later parities), in different lactation stages (early vs. late lactation), and in herds with different frequencies of claw disorders (low vs. high frequency). Analyses revealed that heritabilities measured in different parities, lactation stages, or herds with different trait frequencies are similar for most claw health traits. Also, genetic correlations (rg) for most claw health traits were not found to be different from unity for traits in different parities, lactation stages, or herds with different trait frequencies. Sole hemorrhage and infectious lesions were genetically different traits in first or later parities (rg=0.29±0.31 and 0.66±0.15, respectively). White line separation and infectious lesions were genetically different in early and late lactation (rg=0.53±0.20 and 0.69±0.13, respectively), and sole ulcer was genetically different in herds with low or high frequency of sole ulcer (rg=0.75±0.14). In our analysis, we did not find convincing evidence to supports the use of multiple trait models for the analysis of claw health traits treating them as different traits in different parities, lactation stages, and herds with different claw disorder frequency.

Genome-wide association study for claw disorders and trimming status in dairy cattle

Performing a genome-wide association study (GWAS) might add to a better understanding of the development of claw disorders and the need for trimming. Therefore, the aim of the current study was to perform a GWAS on claw disorders and trimming status and to validate the results for claw disorders based on an independent data set. Data consisted of 20,474 cows with phenotypes for claw disorders and 50,238 cows with phenotypes for trimming status. Recorded claw disorders used in the current study were double sole (DS), interdigital hyperplasia (IH), sole hemorrhage (SH), sole ulcer (SU), white line separation (WLS), a combination of infectious claw disorders consisting of (inter-)digital dermatitis and heel erosion, and a combination of laminitis-related claw disorders (DS, SH, SU, and WLS). Of the cows with phenotypes for claw disorders, 1,771 cows were genotyped and these cow data were used for the GWAS on claw disorders. A SNP was considered significant when the false discovery rate≤0.05 and suggestive when the false discovery rate≤0.20. An independent data set of 185 genotyped bulls having at least 5 daughters with phenotypes (6,824 daughters in total) for claw disorders was used to validate significant and suggestive SNP detected based on the cow data. To analyze the trait "trimming status" (i.e., the need for claw trimming), a data set with 327 genotyped bulls having at least 5 daughters with phenotypes (18,525 daughters in total) was used. Based on the cow data, in total 10 significant and 45 suggestive SNP were detected for claw disorders. The 10 significant SNP were associated with SU, and mainly located on BTA8. The suggestive SNP were associated with DS, IH, SU, and laminitis-related claw disorders. Three of the suggestive SNP were validated in the data set of 185 bulls, and were located on BTA13, BTA14, and BTA17. For infectious claw disorders, SH, and WLS, no significant or suggestive SNP associations were detected. For trimming status, 1 significant and 1 suggestive SNP were detected, both located close to each other on BTA15. Some significant and suggestive SNP were located close to SNP detected in studies on feet and leg conformation traits. Genes with major effects could not be detected and SNP associations were spread across the genome, indicating that many SNP, each explaining a small proportion of the genetic variance, influence claw disorders. Therefore, to reduce the incidence of claw disorders by breeding, genomic selection is a promising approach.