Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle.

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

Genome-wide association analyses reveals copy number variant regions associated with reproduction and disease traits in Canadian Holstein cattle.

This study aimed to evaluate the impact of copy number variants (CNVs) on 13 reproduction and 12 disease traits in Holstein cattle. Intensity signal files containing Log R ratio and B allele frequency information from 13,730 Holstein animals genotyped with a 95K SNP panel, and 8,467 Holstein animals genotyped with a 50K SNP panel were used to identify the CNVs. Subsequently, the identified CNVs were validated using whole genome sequence data from 126 animals, resulting in 870 high-confidence CNV regions (CNVRs) on 12,131 animals. Out of these, 54 CNVRs had frequencies higher than or equal to 1% in the population and were used in the genome-wide association analysis (one CNVR at a time, including the G matrix). Results revealed that 4 CNVRs were significantly (p-value < 3.7 × 10-5) associated with at least one of the traits analyzed in this study. Specifically, 2 CNVRs were associated with 3 reproduction traits (i.e., calf survival, first service to conception, and non-return rate), and 2 CNVRs were associated with 2 disease traits (i.e., metritis and retained placenta). These CNVRs harbored genes implicated in immune response, cellular signaling, and neuronal development, supporting their potential involvement in these traits. Further investigations to unravel the mechanistic and functional implications of these CNVRs on the mentioned traits are warranted.

The Resilient Dairy Genome Project-A general overview of methods and objectives related to feed efficiency and methane emissions.

The Resilient Dairy Genome Project (RDGP) is an international large-scale applied research project that aims to generate genomic tools to breed more resilient dairy cows. In this context, improving feed efficiency and reducing greenhouse gases from dairy is a high priority. The inclusion of traits related to feed efficiency (e.g., dry matter intake [DMI]) or greenhouse gases (e.g., methane emissions [CH4]) relies on available genotypes as well as high quality phenotypes. Currently, 7 countries (i.e., Australia, Canada, Denmark, Germany, Spain, Switzerland, and United States) contribute with genotypes and phenotypes including DMI and CH4. However, combining data are challenging due to differences in recording protocols, measurement technology, genotyping, and animal management across sources. In this study, we provide an overview of how the RDGP partners address these issues to advance international collaboration to generate genomic tools for resilient dairy. Specifically, we describe the current state of the RDGP database, data collection protocols in each country, and the strategies used for managing the shared data. As of February 2022, the database contains 1,289,593 DMI records from 12,687 cows and 17,403 CH4 records from 3,093 cows and continues to grow as countries upload new data over the coming years. No strong genomic differentiation between the populations was identified in this study, which may be beneficial for eventual across-country genomic predictions. Moreover, our results reinforce the need to account for the heterogeneity in the DMI and CH4 phenotypes in genomic analysis.

Underlying genetic architecture of resistance to mastitis in dairy cattle: A systematic review and gene prioritization analysis of genome-wide association studies.

Mastitis, the most frequent disease in dairy cattle. Resistance to mastitis is a complex, polygenic trait controlled by several genes, each with small effects. Genome-wide association studies have been widely used to identify genomic variants associated with complex traits, including resistance to mastitis, to elucidate the underlying genetic architecture of the trait. However, no systematic review and gene prioritization analysis have been conducted to date on GWAS results for resistance to mastitis in dairy cattle. Hence, the objective was to perform a systematic review and gene prioritization analysis of GWAS studies to identify potential functional candidate genes associated with resistance to mastitis-related traits in dairy cattle. Four electronic databases were searched from inception to December 2020, supplemented with multiple sources of gray literature, to identify eligible articles. Annotation for genes and quantitative trait loci (QTL), and QTL enrichment analysis were conducted using GALLO. Gene prioritization analysis was performed by a guilty-by-association approach using GUILDify and ToppGene. From 52 articles included within this systematic review, 30 articles were used for further functional analyses. Gene and QTL annotation resulted in 9,125 and 43,646 unique genes and QTL, respectively, from 39 studies. In general, overlapping of genes across studies was very low (mean ± SD = 0.02% ± 0.07%). Most annotated genes were associated with somatic cell count-related traits and the Holstein breed. Within all annotated genes, 74 genes were shared among Holstein, Jersey, and Ayrshire breeds. Approximately 7.5% of annotated QTL were related to QTL class "health." Within the health QTL class, 2.6 and 2.2% of QTL were associated with clinical mastitis and somatic cell count-related traits. Enrichment analysis of QTL demonstrated that many enriched QTL were associated with somatic cell score located in Bos taurus autosomes 5, 6, 16, and 20. The prioritization analysis resulted in 427 significant genes after multiple test correction (false discovery rate of 5%) from 26 studies. Most prioritized genes were located in Bos taurus autosomes 19 and 7, and most top-ranked genes were from the cytokine superfamily (e.g., chemokines, interleukins, transforming growth factors, and tumor necrosis factor genes). Although most prioritized genes (397) were associated with somatic cell count-related traits, only 54 genes were associated with clinical mastitis-related traits. Twenty-four genes (ABCC9, ACHE, ADCYAP1, ARC, BCL2L1, CDKN1A, EPO, GABBR2, GDNF, GNRHR, IKBKE, JAG1, KCNJ8, KCNQ1, LIFR, MC3R, MYOZ3, NFKB1, OSMR, PPP3CA, PRLR, SHARPIN, SLC1A3, and TNFRSF25) were reported for both somatic cell count and clinical mastitis-related traits. Prioritized genes were mainly associated with immune response, regulation of secretion, locomotion, cell proliferation, and development. In conclusion, this study provided a fine-mapping of previously identified genomic regions associated with resistance to mastitis and identified key functional candidate genes for resistance to mastitis, which can be used to develop enhanced genomic strategies to combat mastitis by increasing mastitis resistance through genetic selection.

Sharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency.

BACKGROUND: Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows. RESULTS: GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (rg) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase. CONCLUSIONS: The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.

Reproductive tract size and position score: Estimation of genetic parameters for a novel fertility trait in dairy cows.

The dairy industry is moving toward selecting animals with better fertility to decrease the economic losses linked to reproductive issues. The reproductive tract size and position score (SPS) was recently developed in physiological studies as an indicator of pregnancy rate and the number of services to conception. Cows are scored as SPS 1, 2, or 3 based on the size of their reproductive tract and its position in the pelvis, as determined by transrectal palpation. The objective of this study was to estimate genetic parameters for SPS to assess its potential as a novel fertility trait. Phenotypes were collected at the University of British Columbia's research herd from 2017 to 2020, consisting of 3,247 within- and across-lactation SPS records from 490 Holstein cows. A univariate animal model was used to estimate the variance components for SPS. Both threshold and linear models were fit under a Bayesian approach and the results were compared using the Spearman rank correlation (r) between the estimated breeding values. The 2 models ranked the animals very similarly (r = 0.99), and the linear model was selected for further analysis. Genetic correlations with other currently evaluated traits were estimated using a bivariate animal model. The posterior means (± posterior standard deviation) for heritability and repeatability within- and across-lactation were 0.113 (± 0.013), 0.242 (± 0.012), and 0.134 (± 0.014), respectively. The SPS showed null correlations with production traits and favorable correlations with traditional fertility traits, varying from -0.730 (nonreturn rate) to 0.931 (number of services). Although preliminary, these results are encouraging because SPS seems to be more heritable than and strongly genetically correlated with number of services, nonreturn rate, and first service to conception, indicating potential for effective indirect selection response on these traits from SPS genetic selection. Therefore, further studies with larger data sets to validate these findings are warranted.

Predicting dry matter intake in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks.

Dry matter intake (DMI) is a fundamental component of the animal's feed efficiency, but measuring DMI of individual cows is expensive. Mid-infrared reflectance spectroscopy (MIRS) on milk samples could be an inexpensive alternative to predict DMI. The objectives of this study were (1) to assess if milk MIRS data could improve DMI predictions of Canadian Holstein cows using artificial neural networks (ANN); (2) to investigate the ability of different ANN architectures to predict unobserved DMI; and (3) to validate the robustness of developed prediction models. A total of 7,398 milk samples from 509 dairy cows distributed over Canada, Denmark, and the United States were analyzed. Data from Denmark and the United States were used to increase the training data size and variability to improve the generalization of the prediction models over the lactation. For each milk spectra record, the corresponding weekly average DMI (kg/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d), metabolic body weight (MBW), age at calving, year of calving, season of calving, days in milk, lactation number, country, and herd were available. The weekly average DMI was predicted with various ANN architectures using 7 predictor sets, which were created by different combinations MY, FY, PY, MBW, and MIRS data. All predictor sets also included age of calving and days in milk. In addition, the classification effects of season of calving, country, and lactation number were included in all models. The explored ANN architectures consisted of 3 training algorithms (Bayesian regularization, Levenberg-Marquardt, and scaled conjugate gradient), 2 types of activation functions (hyperbolic tangent and linear), and from 1 to 10 neurons in hidden layers). In addition, partial least squares regression was also applied to predict the DMI. Models were compared using cross-validation based on leaving out 10% of records (validation A) and leaving out 10% of cows (validation B). Superior fitting statistics of models comprising MIRS information compared with the models fitting milk, fat and protein yields suggest that other unknown milk components may help explain variation in weekly average DMI. For instance, using MY, FY, PY, and MBW as predictor variables produced a predictive accuracy (r) ranging from 0.510 to 0.652 across ANN models and validation sets. Using MIRS together with MY, FY, PY, and MBW as predictors resulted in improved fitting (r = 0.679-0.777). Including MIRS data improved the weekly average DMI prediction of Canadian Holstein cows, but it seems that MIRS predicts DMI mostly through its association with milk production traits and its utility to estimate a measure of feed efficiency that accounts for the level of production, such as residual feed intake, might be limited and needs further investigation. The better predictive ability of nonlinear ANN compared with linear ANN and partial least squares regression indicated possible nonlinear relationships between weekly average DMI and the predictor variables. In general, ANN using Bayesian regularization and scaled conjugate gradient training algorithms yielded slightly better weekly average DMI predictions compared with ANN using the Levenberg-Marquardt training algorithm.

Predicting methane emission in Canadian Holstein dairy cattle using milk mid-infrared reflectance spectroscopy and other commonly available predictors via artificial neural networks.

Interest in reducing eructed CH4 is growing, but measuring CH4 emissions is expensive and difficult in large populations. In this study, we investigated the effectiveness of milk mid-infrared spectroscopy (MIRS) data to predict CH4 emission in lactating Canadian Holstein cows. A total of 181 weekly average CH4 records from 158 Canadian cows and 217 records from 44 Danish cows were used. For each milk spectra record, the corresponding weekly average CH4 emission (g/d), test-day milk yield (MY, kg/d), fat yield (FY, g/d), and protein yield (PY, g/d) were available. The weekly average CH4 emission was predicted using various artificial neural networks (ANN), partial least squares regression, and different sets of predictors. The ANN architectures consisted of 3 training algorithms, 1 to 10 neurons with hyperbolic tangent activation function in the hidden layer, and 1 neuron with linear (purine) activation function in the hidden layer. Random cross-validation was used to compared the predictor sets: MY (set 1); FY (set 2); PY (set 3); MY and FY (set 4); MY and PY (set 5); MY, FY, and PY (set 6); MIRS (set 7); and MY, FY, PY, and MIRS (set 8). All predictor sets also included age at calving and days in milk, in addition to country, season of calving, and lactation number as categorical effects. Using only MY (set 1), the predictive accuracy (r) ranged from 0.245 to 0.457 and the root mean square error (RMSE) ranged from 87.28 to 99.39 across all prediction models and validation sets. Replacing MY with FY (set 2; r = 0.288-0.491; RMSE = 85.94-98.04) improved the predictive accuracy, but using PY (set 3; r = 0.260-0.468; RMSE = 86.95-98.47) did not. Adding FY to MY (set 4; r = 0.272-0.469; RMSE = 87.21-100.76) led to a negligible improvement compared with sets 1 and 3, but it slightly decreased accuracy compared with set 2. Adding PY to MY (set 5; r = 0.250-0.451; RMSE = 87.66-100.94) did not improve prediction ability. Combining MY, FY, and PY (set 6; r = 0.252-0.455; RMSE = 87.74-101.93) yielded accuracy slightly lower than sets 2 and 3. Using only MIRS data (set 7; r = 0.586-0.717; RMSE = 69.09-96.20) resulted in superior accuracy compared with all previous sets. Finally, the combination of MIRS data with MY, FY, and PY (set 8; r = 0.590-0.727; RMSE = 68.02-87.78) yielded similar accuracy to set 7. Overall, sets including the MIRS data yielded significantly better predictions than the other sets. To assess the predictive ability in a new unseen herd, a limited block cross-validation was performed using 20 cows in the same Canadian herd, which yielded r = 0.229 and RMSE = 154.44, which were clearly much worse than the average r = 0.704 and RMSE = 70.83 when predictions were made by random cross-validation. These results warrant further investigation when more data become available to allow for a more comprehensive block cross-validation before applying the calibrated models for large-scale prediction of CH4 emissions.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages.

BACKGROUND: Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne's disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (-); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip. RESULTS: The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/-) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10- 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids. CONCLUSION: In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins.

BACKGROUND: The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. METHODS: In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. RESULTS: In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a - 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. CONCLUSIONS: ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations.

Genetic analysis of pathogen-specific intramammary infections in dairy cows.

Mastitis is one of the most common diseases in dairy cattle, causing severe economic losses to dairy farmers. Mastitis usually occurs due to intramammary infection (IMI) caused by a variety of pathogenic bacteria. Although good progress has been made in understanding genetics of pathogen-specific clinical mastitis, studies involving genetic analysis of pathogen-specific IMI are scarce. The overall objective of this study was, therefore, to assess genetic variation of overall and pathogen-specific IMI in nonclinical primiparous and multiparous cows using bacterial culture. Data and milk samples were collected over a 2-yr interval as part of the Canadian Bovine Mastitis Research Network. The final data set contained records of 46,900 quarter milk samples from 3,382 clinically healthy primiparous and multiparous Holstein cows from 84 dairy herds. For the genetic analysis, we considered the following 7 traits: overall IMI, non-aureus staphylococci (NAS) IMI, contagious pathogen IMI, environmental pathogen IMI, major pathogen IMI, minor pathogen IMI and somatic cell score (SCS). Data were analyzed at the quarter level using a threshold-probit model via Gibbs sampling in BLUPF90. Prevalence of IMI traits at the quarter level in multiparous cow from 0 to 400 DIM ranged from 6.8 to 45.5%. Posterior mean of quarter heritability estimates (on the underlying scale, posterior SD in brackets) of overall IMI and pathogen-specific IMI traits ranged from 0.017 to 0.073 (±0.009 to 0.030). Weak to strong genetic correlations [ranging from 0.18 to 0.97 (±0.01 to 0.29)] among pathogen-specific IMI traits and with overall IMI indicated that not all of these traits were genetically similar. Weak to moderate Spearman rank correlations between estimated breeding values for overall IMI and pathogen-specific IMI traits (from 0.31 to 0.87) indicated possible substantial reranking of sires. The percentage of daughters with IMI caused by various pathogen groups ranged from 13 to 80% and from 38 to 94% for the best (10% decile) and worst sires (90% decile) according to their IMI trait-specific estimated breeding values, respectively. Pathogen-specific IMI traits and overall IMI had weak to moderate positive genetic correlations [ranging from 0.11 to 0.81 (±0.11 to 0.22)] with SCS. Therefore, selection for lower SCS will improve resistance to IMI. However, based on the observed weak to moderate rank correlations (0.04 to 0.47) between pathogen-specific IMI traits and SCS, selection for lower SCC will not improve resistance to IMI from every pathogen-specific IMI group in the same manner. Therefore, despite low heritability estimates, there was sizeable genetic variation for pathogen-specific IMI traits, indicating that long-term direct genetic selection for pathogen-specific IMI can improve pathogen-specific IMI resistance.

Genome-wide association study between copy number variants and hoof health traits in Holstein dairy cattle.

Genome-wide association studies based on SNP have been completed for multiple traits in dairy cattle; however, copy number variants (CNV) could add genomic information that has yet to be harnessed. The objectives of this study were to identify CNV in genotyped Holstein animals and assess their association with hoof health traits using deregressed estimated breeding values as pseudophenotypes. A total of 23,256 CNV comprising 1,645 genomic regions were identified in 5,845 animals. Fourteen genomic regions harboring structural variations, including 9 deletions and 5 duplications, were associated with at least 1 of the studied hoof health traits. This group of traits included digital dermatitis, interdigital dermatitis, heel horn erosion, sole ulcer, white line lesion, sole hemorrhage, and interdigital hyperplasia; no regions were associated with toe ulcer. Twenty candidate genes overlapped with the regions associated with these traits including SCART1, NRXN2, KIF26A, GPHN, and OR7A17. In this study, an effect on infectious hoof lesions could be attributed to the PRAME (Preferentially Expressed Antigen in Melanoma) gene. Almost all genes detected in association with noninfectious hoof lesions could be linked to known metabolic disorders. The knowledge obtained considering information of associated CNV to the traits of interest in this study could improve the accuracy of estimated breeding values. This may further increase the genetic gain for these traits in the Canadian Holstein population, thus reducing the involuntary animal losses due to lameness.

Effect of recent and ancient inbreeding on production and fertility traits in Canadian Holsteins.

BACKGROUND: Phenotypic performances of livestock animals decline with increasing levels of inbreeding, however, the noticeable decline known as inbreeding depression, may not be due only to the total level of inbreeding, but rather could be distinctly associated with more recent or more ancient inbreeding. Therefore, splitting inbreeding into different age classes could help in assessing detrimental effects of different ages of inbreeding. Hence, this study sought to investigate the effect of recent and ancient inbreeding on production and fertility traits in Canadian Holstein cattle with both pedigree and genomic records. Furthermore, inbreeding coefficients were estimated using traditional pedigree measure (FPED) and genomic measures using segment based (FROH) and marker-by-marker (FGRM) based approaches. RESULTS: Inbreeding depression was found for all production and most fertility traits, for example, every 1% increase in FPED, FROH and FGRM was observed to cause a - 44.71, - 40.48 and - 48.72 kg reduction in 305-day milk yield (MY), respectively. Similarly, an extension in first service to conception (FSTC) of 0.29, 0.24 and 0.31 day in heifers was found for every 1% increase in FPED, FROH and FGRM, respectively. Fertility traits that did not show significant depression were observed to move in an unfavorable direction over time. Splitting both pedigree and genomic inbreeding into age classes resulted in recent age classes showing more detrimental inbreeding effects, while more distant age classes caused more favorable effects. For example, a - 1.56 kg loss in 305-day protein yield (PY) was observed for every 1% increase in the most recent pedigree age class, whereas a 1.33 kg gain was found per 1% increase in the most distant pedigree age class. CONCLUSIONS: Inbreeding depression was observed for production and fertility traits. In general, recent inbreeding had unfavorable effects, while ancestral inbreeding had favorable effects. Given that more negative effects were estimated from recent inbreeding when compared to ancient inbreeding suggests that recent inbreeding should be the primary focus of selection programs. Also, further work to identify specific recent homozygous regions negatively associated with phenotypic traits could be investigated.

Association of genetic polymorphisms related to Johne's disease with estimated breeding values of Holstein sires for milk ELISA test scores.

BACKGROUND: Johne's disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls' estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle. RESULTS: Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls' EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively. CONCLUSIONS: The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs.

Effect of genomic selection on rate of inbreeding and coancestry and effective population size of Holstein and Jersey cattle populations.

Genetic diversity in livestock populations is a significant contributor to the sustainability of animal production. Also, genetic diversity allows animal production to become more responsive to environmental changes and market demands. The loss of genetic diversity can result in a plateau in production and may also result in loss of fitness or viability in animal production. In this study, we investigated the rate of inbreeding (ΔF), rate of coancestry (Δf), and effective population size (Ne) as important quantitative indicators of genetic diversity and evaluated the effect of the recent implementation of genomic selection on the loss of genetic diversity in North American Holstein and Jersey dairy cattle. To estimate the rate of inbreeding and coancestry, inbreeding and coancestry coefficients were calculated using the traditional pedigree method and genomic methods estimated from segment- and marker-based approaches. Furthermore, we estimated Ne from the rate of inbreeding and coancestry and extent of linkage disequilibrium. A total of 205,755 and 89,238 pedigreed and genotyped animals born between 1990 and 2018 inclusively were available for Holsteins and Jerseys, respectively. The estimated average pedigree inbreeding coefficients were 7.74 and 7.20% for Holsteins and Jerseys, respectively. The corresponding values for the segment and marker-by-marker genomic inbreeding coefficients were 13.61, 15.64, and 31.40% for Holsteins and 21.16, 22.54, and 42.62% for Jerseys, respectively. The average coancestry coefficients were 8.33 and 15.84% for Holsteins and 9.23 and 23.46% for Jerseys with pedigree and genomic measures, respectively. Generation interval for the whole 29-yr time period averaged approximately 5 yr for all selection pathways combined. The ΔF per generation based on pedigree, segment, and marker-by-marker genomic measures for the entire 29-yr period was estimated to be 0.75, 1.10, 1.16, and 1.02% for Holstein animals and 0.67, 0.62, 0.63, and 0.59% for Jersey animals, respectively. The Δf was estimated to be 0.98 and 0.98% for Holsteins and 0.73 and 0.78% for Jerseys with pedigree and genomic measures, respectively. These ΔF and Δf translated to an Ne that ranged from 43 to 66 animals for Holsteins and 64 to 85 animals for Jerseys. In addition, the Ne based on linkage disequilibrium was 58 and 120 for Holsteins and Jerseys, respectively. The 10-yr period that involved the application of genomic selection resulted in an increased ΔF per generation with ranges from 1.19 to 2.06% for pedigree and genomic measures in Holsteins. Given the rate at which inbreeding is increasing after the implementation of genomic selection, there is a need to implement measures and means for controlling the rate of inbreeding per year, which will help to manage and maintain farm animal genetic resources.

Estimation of genetic parameters for mid-infrared-predicted lactoferrin and milk fat globule size in Holstein cattle.

Lactoferrin (LF) and milk fat globule (MFG) are 2 biologically active components of milk with great economical and nutritional value in the dairy industry. The objectives of this study were to estimate (1) the heritability of mid-infrared (MIR)-predicted LF and MFG size (MFGS) and (2) the genetic correlations between predicted LF and MFGS with milk, fat, and protein yields, fat and protein percentages, and somatic cell score in first-parity Canadian Holstein cattle. A total of 109,029 test-day records from 22,432 cows and 1,572 farms for MIR-predicted LF and 109,212 test-day records from 22,424 cows and 1,559 farms for MIR-predicted MFGS were used in the analyses. Four separate 5-trait random regression test-day models were used. The models included days in milk, herd test date, and a polynomial regression on DIM nested in age-season of calving classes as fixed effects, random polynomial regressions on DIM nested in herd-year of calving, animal additive genetic and permanent environment classes, and a residual effect. Regression curves were modeled using orthogonal Legendre polynomials of order 4 for the fixed age-season of calving effect and of order 5 for the random effects. Moderate overall heritability estimates of 0.34 and 0.46 were estimated for the MIR-predicted LF and MIR-predicted MFGS, respectively. These heritability estimates were similar to the ones estimated for the direct measure of MFGS in a previous study. The genetic correlations between predicted MFGS and fat percentage (0.53) and between predicted LF and protein percentage (0.41) were both moderate and positive. Predicted LF and somatic cell score showed a weaker correlation (0.06) compared with other studies. The moderate genetic correlation between MIR-predicted MFGS and fat percentage and between MIR-predicted LF and protein percentage suggests that MIR predictions of MFGS and LF are not simply a function of the amount of fat and protein percentage, respectively, in the milk (i.e., the prediction equations are not simply predicting fat or protein percentages). Thus, these MIR-predicted values may provide additional information for selecting for fine milk components in Holstein cattle.

Symposium review: The genomic architecture of inbreeding: How homozygosity affects health and performance.

Inbreeding depression is a growing concern in livestock because it can detrimentally affect animal fitness, health, and production levels. Genomic information can be used to more effectively capture variance in Mendelian sampling, thereby enabling more accurate estimation of inbreeding, but further progress is still required. The calculation of inbreeding for herd management purposes is largely still done using pedigree information only, although inbreeding coefficients calculated in this manner have been shown to be less accurate than genomic inbreeding measures. Continuous stretches of homozygous genotypes, so called runs of homozygosity, have been shown to provide a better estimate of autozygosity at the genomic level than conventional measures based on inbreeding coefficients calculated through conventional pedigree information or even genomic relationship matrices. For improved and targeted management of genomic inbreeding at the population level, the development of methods that incorporate genomic information in mate selection programs may provide a more precise tool for reducing the detrimental effects of inbreeding in dairy herds. Additionally, a better understanding of the genomic architecture of inbreeding and incorporating that knowledge into breeding programs could significantly refine current practices. Opportunities to maintain high levels of genetic progress in traits of interest while managing homozygosity and sustaining acceptable levels of heterozygosity in highly selected dairy populations exist and should be examined more closely for continued sustainability of both the dairy cattle population as well as the dairy industry. The inclusion of precise genomic measures of inbreeding, such as runs of homozygosity, inbreeding, and mating programs, may provide a path forward. In this symposium review article, we describe traditional measures of inbreeding and the recent developments made toward more precise measures of homozygosity using genomic information. The effects of homozygosity resulting from inbreeding on phenotypes, the identification and mapping of detrimental homozygosity haplotypes, management of inbreeding with genomic data, and areas in need of further research are discussed.

Genetic analysis of subclinical mastitis in early lactation of heifers using both linear and threshold models.

Subclinical mastitis (SCM) causes economic losses for dairy producers by reducing milk production and leading to higher incidence of clinical mastitis and premature culling. The prevalence of SCM in first-lactation heifers is highest during early lactation. The objective of this study was to estimate genetic parameters for SCM in early lactation in first-parity Holsteins. Somatic cell count test-day records were collected monthly in 91 Canadian herds participating in the National Cohort of Dairy Farms of the Canadian Bovine Mastitis Research Network. Only the first test-day record available between 5 and 30 d in milk was considered for analysis. The final data set contained 8,518 records from first lactation Holstein heifers. Six alternative traits were defined as indicators of SCM, using various cutoff values of SCC, ranging from 150,000 to 400,000 cells/mL. Both linear and threshold animal models were used. Overall prevalence of SCM using the 6 traits ranged from 13 to 24%. Heritability estimates (standard error) from linear and threshold models ranged from 0.037 to 0.057 (0.015 to 0.018) and from 0.040 to 0.051 (0.017 to 0.020), respectively. We found strong genetic correlations (standard error) among alternative SCC traits, ranging from 0.90 to 0.99 (0.013 to 0.069), indicating that these 6 traits were genetically similar. Despite low heritability, based on estimated breeding values (EBV) predicted from both models, we noted exploitable genetic variation among sires. Higher EBV of SCM resistance corresponded to sires with a higher percentage of daughters without SCM. Based on a linear model (all 6 traits), percentage of daughters with SCM ranged from 5 to 13% and from 19 to 33% for the top 10% and worst 10% of 69 sires with minimum 20 daughters in at least 5 herds, respectively. Spearman's rank correlations among EBV of sires predicted from linear (from 0.75 to 0.95) and threshold (from 0.74 to 0.95) models were moderate to high, respectively. Very high rank correlations (0.98 to 0.99) between EBV predicted for the same trait from linear and threshold model indicated that reranking of sires based on model used was minimal. In conclusion, despite low heritability, we found utilizable genetic variation in early lactation of heifers. Hence, genetic selection to improve genetic resistance to SCM in early lactation of heifers was deemed possible.

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.

A genome-wide association study of immune response traits in Canadian Holstein cattle.

BACKGROUND: Breeding for enhanced immune response (IR) has been suggested as a tool to improve inherent animal health. Dairy cows with superior antibody-mediated (AMIR) and cell-mediated immune responses (CMIR) have been demonstrated to have a lower occurrence of many diseases including mastitis. Adaptive immune response traits are heritable, and it is, therefore, possible to breed for improved IR, decreasing the occurrence of disease. The objective of this study was to perform genome-wide association studies to determine differences in genetic profiles among Holstein cows classified as High or Low for AMIR and CMIR. From a total of 680 cows with immune response phenotypes, 163 cows for AMIR (81 High and 82 Low) and 140 for CMIR (75 High and 65 Low) were selectively genotyped using the Illumina Bovine SNP50 BeadChip. Results were validated using an unrelated population of 164 Holstein bulls IR phenotyped for AMIR and 146 for CMIR. RESULTS: A generalized quasi likelihood score method was used to determine single nucleotide polymorphisms (SNP) and chromosomal regions associated with immune response. After applying a 5% chromosomal false discovery rate, 186 SNPs were significantly associated with AMIR. The majority (93%) of significant markers were on chromosome 23, with a similar peak found in the bull population. For CMIR, 21 SNP markers remained significant. Candidate genes within 250,000 base pairs of significant SNPs were identified to determine biological pathways associated with AMIR and CMIR. Various pathways were identified, including the antigen processing and presentation pathway, important in host defense. Candidate genes included those within the bovine Major Histocompatability Complex such as BoLA-DQ, BoLA-DR and the non-classical BoLA-NC1 for AMIR and BoLA-DQ for CMIR, the complement system including C2 and C4 for AMIR and C1q for CMIR, and cytokines including IL-17A, IL17F for AMIR and IL-17RA for CMIR and tumor necrosis factor for both AMIR and CMIR. Additional genes associated with CMIR included galectins 1, 2 and 3, BCL2 and ß-defensin. CONCLUSIONS: The significant genetic variation associated with AMIR and CMIR in this study may imply feasibility to include immune response in genomic breeding indices as an approach to improve inherent animal health.