Estimation of genetic parameters and genome-wide association study for enzootic bovine leukosis resistance in Canadian Holstein cattle.

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis (leukosis) frequently observed in North American dairy herds. Infection with BLV can lead to persistent lymphocytosis and tumors, and is associated with decreased production, immunity and fertility. With no available treatment or vaccine, reducing the prevalence of leukosis through management and culling has not yet been successful. Genetic selection could contribute to permanent improvement in dairy cattle resistance to leukosis. This study aimed to examine the prevalence and impact of leukosis in Canada, and to assess the potential for including leukosis resistance in Canadian national genetic evaluations by characterizing the genetic architecture of leukosis resistance using pedigree and genomic information. A total of 117,349 milk enzyme-linked immunosorbent assay test records on 96,779 Holstein cows from 950 Canadian herds taken between 2007 and 2021 were provided by Lactanet Canada (Guelph, ON, Canada). Each cow was classified as test-positive for leukosis or test-negative for leukosis. Leukosis was present in approximately 77% of herds tested; within those herds, an average of 39% of cows tested were test positive for leukosis. Heritabilities of 0.10 (SE = 0.001) and 0.07 (SE <0.001) were estimated for leukosis resistance using a linear animal model and BLUP or single-step GBLUP methodology, respectively. Breeding value correlations were estimated between leukosis resistance and economically important and phenotypically relevant traits. Most correlations between leukosis resistance and traits already included in Canadian genetic evaluations were favorable, with the exception of somatic cell score. The candidate genes for leukosis resistance identified using a genome-wide association study, were on chromosome 23, with some being part of the major histocompatibility complex. This study showed that genetic evaluation for leukosis resistance is possible, and could be considered for inclusion in Canadian national selection indices.

Investigating the genomic background of calving-related traits in Canadian Jersey cattle.

Traits related to calving have a significant impact on animal welfare and farm profitability in dairy production systems. Identifying genomic regions associated with calving traits could contribute to refining dairy cattle breeding programs and management practices in the dairy industry. Therefore, the primary objectives of this study were to estimate genetic parameters and perform genome-wide association studies (GWAS) and functional enrichment analyses for stillbirth, gestation length, calf size, and calving ease traits in North American Jersey cattle. A total of 40,503 animals with phenotypic records and 5,398 animals genotyped for 45,101 single nucleotide polymorphisms (SNPs) were included in the analyses. Genetic parameters were estimated based on animal models and Bayesian methods. The effects of SNPs were estimated using the Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) method. The heritability (standard error) estimates ranged from 0.01 (0.01) for stillbirths (SB) in heifers to 0.11 (0.01) for gestation length (GL) in cows. The genetic correlations ranged from -0.58 (0.11) between calving ease (CE) and SB in heifers to 0.44 (0.14) between calving ease and calf size (CZ) in cows. CE showed the highest genetic correlation between heifers and cows, 0.8 (0.22) respectively. The candidate genes identified, including MTHFR, SERPINA5, IGFBP3, and ZRANB1, are involved in key biological processes and metabolic pathways related to the studied traits. Reducing environmental variation and identifying novel indicators of reproduction traits in the Jersey breed are needed given the low heritability estimates for most traits evaluated in this study. In conclusion, this study provides a characterization of the genetic background of calving-related traits in Jersey cattle. The estimates obtained can be used to improve or build selection indexes in Jersey cattle breeding programs in North America.

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application.

BACKGROUND: Structural variants (SVs) such as deletions, duplications, and insertions are known to contribute to phenotypic variation but remain challenging to identify and genotype. A more complete, accessible, and assessable collection of SVs will assist efforts to study SV function in cattle and to incorporate SV genotyping into animal evaluation. RESULTS: In this work we produced a large and deeply characterized collection of SVs in Holstein cattle using two popular SV callers (Manta and Smoove) and publicly available Illumina whole-genome sequence (WGS) read sets from 310 samples (290 male, 20 female, mean 20X coverage). Manta and Smoove identified 31 K and 68 K SVs, respectively. In total the SVs cover 5% (Manta) and 6% (Smoove) of the reference genome, in contrast to the 1% impacted by SNPs and indels. SV genotypes from each caller were confirmed to accurately recapitulate animal relationships estimated using WGS SNP genotypes from the same dataset, with Manta genotypes outperforming Smoove, and deletions outperforming duplications. To support efforts to link the SVs to phenotypic variation, overlapping and tag SNPs were identified for each SV, using genotype sets extracted from the WGS results corresponding to two bovine SNP chips (BovineSNP50 and BovineHD). 9% (Manta) and 11% (Smoove) of the SVs were found to have overlapping BovineHD panel SNPs, while 21% (Manta) and 9% (Smoove) have BovineHD panel tag SNPs. A custom interactive database ( https://svdb-dc.pslab.ca ) containing the identified sequence variants with extensive annotations, gene feature information, and BAM file content for all SVs was created to enable the evaluation and prioritization of SVs for further study. Illustrative examples involving the genes POPDC3, ORM1, G2E3, FANCI, TFB1M, FOXC2, N4BP2, GSTA3, and COPA show how this resource can be used to find well-supported genic SVs, determine SV breakpoints, design genotyping approaches, and identify processed pseudogenes masquerading as deletions. CONCLUSIONS: The resources developed through this study can be used to explore sequence variation in Holstein cattle and to develop strategies for studying SVs of interest. The lack of overlapping and tag SNPs from commonly used SNP chips for most of the SVs suggests that other genotyping approaches will be needed (for example direct genotyping) to understand their potential contributions to phenotype. The included SV genotype assessments point to challenges in characterizing SVs, especially duplications, using short-read data and support ongoing efforts to better characterize cattle genomes through long-read sequencing. Lastly, the identification of previously known functional SVs and additional CDS-overlapping SVs supports the phenotypic relevance of this dataset.

Genome-wide association analyses reveal 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 (CNV) 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 copy number variant regions (CNVR) on 12,131 animals. Out of these, 54 CNVR 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 CNVR were significantly associated with at least one of the traits analyzed in this study. Specifically, 2 CNVR were associated with 3 reproduction traits (i.e., calf survival, first service to conception, and nonreturn rate), and 2 CNVR were associated with 2 disease traits (i.e., metritis and retained placenta). These CNVR 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 CNVR 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.

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.

A comprehensive characterization of longevity and culling reasons in Canadian Holstein cattle based on various systematic factors.

The decision of premature culling cows directly impacts the profitability of dairy farms. A comprehensive characterization of the primary causes of culling reasons would greatly improve both management and selection objectives in dairy cattle breeding programs. Therefore, this study aimed to analyze the temporal frequencies of 34 culling reasons in Canadian Holstein cows. After data editing and quality control, records from 3,096,872 cows culled from 9,683 herds spread across Canada were used for the analyses covering the periods from 1996 to 2020. Reproductive issues were the main culling reason accounting for 23.02%, followed by milk production (20.82%), health (20.39%), conformation problems (13.69%), economic factors (13.10%), accidents (5.67%), age-related causes (1.67%), and workability (1.63%). Nearly fifty-eight percent of cows were culled after 47 months of age. The observed frequencies of culling due to economic factors were lower than expected from 1996 to 2014 and higher than expected between 2015 and 2020. Reproduction issues had the highest culling frequencies during fall (24.54%), winter (24.02%), and spring (22.51%), while health issues were the most frequent (22.51%) culling reason in the summer season. Health issues (25.50%) and milk production (27.71%) were the most frequent culling reasons in the provinces of Quebec and Ontario, respectively. Reproductive issues showed the highest frequency across climates based on the Köppen climate classification, except for Csb (Dry-summer subtropical or Mediterranean climate) and Bsk (Middle latitude steppe climate), which correspond to small regions in Canada, where production was the most frequent culling reason (29.42% and 21.56%, respectively). Reproductive and milk performance issues were the two main culling reasons in most ecozones, except in Boreal Shield and Atlantic Marine, where health issues had the highest frequencies (25.12 and 23.75%, respectively). These results will contribute to improving management practices and selective decisions to reduce involuntary culling of Holstein cows.

Genetic Analysis of Methane Emission Traits in Holstein Dairy Cattle.

Genetic selection can be a feasible method to help mitigate enteric methane emissions from dairy cattle, as methane emission-related traits are heritable and genetic gains are persistent and cumulative over time. The objective of this study was to estimate heritability of methane emission phenotypes and the genetic and phenotypic correlations between them in Holstein cattle. We used 1765 individual records of methane emission obtained from 330 Holstein cattle from two Canadian herds. Methane emissions were measured using the GreenFeed system, and three methane traits were analyzed: the amount of daily methane produced (g/d), methane yield (g methane/kg dry matter intake), and methane intensity (g methane/kg milk). Genetic parameters were estimated using univariate and bivariate repeatability animal models. Heritability estimates (±SE) of 0.16 (±0.10), 0.27 (±0.12), and 0.21 (±0.14) were obtained for daily methane production, methane yield, and methane intensity, respectively. A high genetic correlation (rg = 0.94 ± 0.23) between daily methane production and methane intensity indicates that selecting for daily methane production would result in lower methane per unit of milk produced. This study provides preliminary estimates of genetic parameters for methane emission traits, suggesting that there is potential to mitigate methane emission in Holstein cattle through genetic selection.

Parent-offspring genotyped trios unravelling genomic regions with gametic and genotypic epistatic transmission bias on the cattle genome.

Several biological mechanisms affecting the sperm and ova fertility and viability at developmental stages of the reproductive cycle resulted in observable transmission ratio distortion (i.e., deviation from Mendelian expectations). Gene-by-gene interactions (or epistasis) could also potentially cause specific transmission ratio distortion patterns at different loci as unfavorable allelic combinations are under-represented, exhibiting deviation from Mendelian proportions. Here, we aimed to detect pairs of loci with epistatic transmission ratio distortion using 283,817 parent-offspring genotyped trios (sire-dam-offspring) of Holstein cattle. Allelic and genotypic parameterization for epistatic transmission ratio distortion were developed and implemented to scan the whole genome. Different epistatic transmission ratio distortion patterns were observed. Using genotypic models, 7, 19 and 6 pairs of genomic regions were found with decisive evidence with additive-by-additive, additive-by-dominance/dominance-by-additive and dominance-by-dominance effects, respectively. Using the allelic transmission ratio distortion model, more insight was gained in understanding the penetrance of single-locus distortions, revealing 17 pairs of SNPs. Scanning for the depletion of individuals carrying pairs of homozygous genotypes for unlinked loci, revealed 56 pairs of SNPs with recessive epistatic transmission ratio distortion patterns. The maximum number of expected homozygous offspring, with none of them observed, was 23. Finally, in this study, we identified candidate genomic regions harboring epistatic interactions with potential biological implications in economically important traits, such as reproduction.

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.

Unravelling the genetics of non-random fertilization associated with gametic incompatibility.

In the dairy industry, mate allocation is dependent on the producer's breeding goals and the parents' breeding values. The probability of pregnancy differs among sire-dam combinations, and the compatibility of a pair may vary due to the combination of gametic haplotypes. Under the hypothesis that incomplete incompatibility would reduce the odds of fertilization, and complete incompatibility would lead to a non-fertilizing or lethal combination, deviation from Mendelian inheritance expectations would be observed for incompatible pairs. By adding an interaction to a transmission ratio distortion (TRD) model, which detects departure from the Mendelian expectations, genomic regions linked to gametic incompatibility can be identified. This study aimed to determine the genetic background of gametic incompatibility in Holstein cattle. A total of 283,817 genotyped Holstein trios were used in a TRD analysis, resulting in 422 significant regions, which contained 2075 positional genes further investigated for network, overrepresentation, and guilt-by-association analyses. The identified biological pathways were associated with immunology and cellular communication and a total of 16 functional candidate genes were identified. Further investigation of gametic incompatibility will provide opportunities to improve mate allocation for the dairy cattle industry.

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.

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.

The potential for mitigation of methane emissions in ruminants through the application of metagenomics, metabolomics, and other -OMICS technologies.

Ruminant supply chains contribute 5.7 gigatons of CO2-eq per annum, which represents approximately 80% of the livestock sector emissions. One of the largest sources of emission in the ruminant sector is methane (CH4), accounting for approximately 40% of the sectors total emissions. With climate change being a growing concern, emphasis is being put on reducing greenhouse gas emissions, including those from ruminant production. Various genetic and environmental factors influence cattle CH4 production, such as breed, genetic makeup, diet, management practices, and physiological status of the host. The influence of genetic variability on CH4 yield in ruminants indicates that genomic selection for reduced CH4 emissions is possible. Although the microbiology of CH4 production has been studied, further research is needed to identify key differences in the host and microbiome genomes and how they interact with one another. The advancement of "-omics" technologies, such as metabolomics and metagenomics, may provide valuable information in this regard. Improved understanding of genetic mechanisms associated with CH4 production and the interaction between the microbiome profile and host genetics will increase the rate of genetic progress for reduced CH4 emissions. Through a systems biology approach, various "-omics" technologies can be combined to unravel genomic regions and genetic markers associated with CH4 production, which can then be used in selective breeding programs. This comprehensive review discusses current challenges in applying genomic selection for reduced CH4 emissions, and the potential for "-omics" technologies, especially metabolomics and metagenomics, to minimize such challenges. The integration and evaluation of different levels of biological information using a systems biology approach is also discussed, which can assist in understanding the underlying genetic mechanisms and biology of CH4 production traits in ruminants and aid in reducing agriculture's overall environmental footprint.

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.

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.

Effects of Incorporating Dry Matter Intake and Residual Feed Intake into a Selection Index for Dairy Cattle Using Deterministic Modeling.

The inclusion of feed efficiency in the breeding goal for dairy cattle has been discussed for many years. The effects of incorporating feed efficiency into a selection index were assessed by indirect selection (dry matter intake) and direct selection (residual feed intake) using deterministic modeling. Both traits were investigated in three ways: (1) restricting the trait genetic gain to zero, (2) applying negative selection pressure, and (3) applying positive selection pressure. Changes in response to selection from economic and genetic gain perspectives were used to evaluate the impact of including feed efficiency with direct or indirect selection in an index. Improving feed efficiency through direct selection on residual feed intake was the best scenario analyzed, with the highest overall economic response including favorable responses to selection for production and feed efficiency. Over time, the response to selection is cumulative, with the potential for animals to reduce consumption by 0.16 kg to 2.7 kg of dry matter per day while maintaining production. As the selection pressure increased on residual feed intake, the response to selection for production, health, and fertility traits and body condition score became increasingly less favorable. This work provides insight into the potential long-term effects of selecting for feed efficiency as residual feed intake.

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.