Genetic improvement is a crucial tool to deal with the increasing demand for high quality, sustainably produced dairy. Breeding programs are based on genetic parameters, such as heritability and genetic correlations, for economically important traits in a population. In this study, we estimated population genetic parameters and genetic trends for 67 traits evaluated on heifers and first-lactation Canadian Holstein cows. The data consisted of approximately 500,000 records with pedigree information collected from 1980 to 2019. Genetic parameters were estimated using bivariate linear animal models under a Bayesian approach. Analyses for the 67 traits resulted in 2,211 bivariate combinations, from which the estimated genetic parameters are reported here. The most highly heritable traits were fat percent (0.66) and protein percent (0.69), followed by stature (0.47). Lowest heritabilities (0.01) were observed for disease-related traits, such as lameness and toe ulcer, and calf survival. The genetic correlations between gestation length, calf size, and calving ease measured on both heifer and cows were close to unity. On the other hand, traits such as body condition score and pin width, cystic ovaries and sole ulcer, rear teat placement, and toe ulcer were genetically unrelated. This study reports genetic parameters that have not been previously published for Canadian Holstein cows, and provides updates of those previously estimated. These estimates are useful for building new indexes, updating existing selection indexes, and for predicting correlated responses due to inclusion of novel traits in the breeding programs.

Zebu animals (Bos indicus) are known to take longer to reach puberty compared with taurine animals (Bos taurus), limiting the supply of animals for harvest or breeding and impacting profitability. Genomic information can be a helpful tool to better understand complex traits and improve genetic gains. In this study, we performed a genomewide association study (GWAS) to identify genetic variants associated with reproductive traits in Nelore beef cattle. Heifer pregnancy (HP) was recorded for 1,267 genotyped animals distributed in 12 contemporary groups (CG) with an average pregnancy rate of 0.35 (±0.01). Disregarding one of these CG, the number of antral follicles (NF) was also collected for 937 of these animals, with an average of 11.53 (±4.43). The animals were organized in CG: 12 and 11 for HP and NF, respectively. Genes in linkage disequilibrium (LD) with the associated variants can be considered in a functional enrichment analysis to identify biological mechanisms involved in fertility. Medical Subject Headings (MeSH) were detected using the MESHR package, allowing the extraction of broad meanings from the gene lists provided by the GWAS. The estimated heritability for HP was 0.28 ± 0.07 and for NF was 0.49 ± 0.09, with the genomic correlation being −0.21 ± 0.29. The average LD between adjacent markers was 0.23 ± 0.01, and GWAS identified genomic windows that accounted for >1% of total genetic variance on chromosomes 5, 14, and 18 for HP and on chromosomes 2, 8, 11, 14, 15, 16, and 22 for NF. The MeSH enrichment analyses revealed significant (P < 0.05) terms associated with HP—“Munc18 Proteins,” “Fucose,” and “Hemoglobins”—and with NF—“Cathepsin B,” “Receptors, Neuropeptide,” and “Palmitic Acid.” This is the first study in Nelore cattle introducing the concept of MeSH analysis. The genomic analyses contributed to a better understanding of the genetic control of the reproductive traits HP and NF and provide new selection strategies to improve beef production. Copyright © 2017. American Society of Animal Science .

About 30% of producers use hormone protocols to synchronize ovulation and perform timed artificial insemination (AI) in Canada. Days from calving to first service (CTFS) and first service to conception (FSTC) become masked phenotypes leading to biased genetic evaluations of cows for these fertility traits. The objectives of this study were to (1) demonstrate and quantify the potential amount of bias in genetic evaluations, and (2) find a procedure that could remove the bias. Simulation was used for both objectives. The proposed solution was to identify cows that have been treated by hormone protocols, make their CTFS and FSTC missing, and perform a multiple trait analysis including traits that have high genetic correlations with CTFS and FSTC, and which are not affected by the hormone protocols themselves. A total of 12 scenarios (S1-S12) were tested, changing the percentage of herds and cows that were randomly selected to be under timed AI. Cows that were given hormone protocols had CTFS of 86 d and FSTC of 0, which were used in genetic evaluation. Four criteria were used to indirectly measure the presence of bias: (1) the correlation between true (TBV) and estimated (EBV) breeding values (accuracy); (2) the differences in the mean EBV of top 25, 50, and 75 sires; (3) changes in correlation between TBV and EBV rankings; and (4) the changes in mean EBV over the simulated generations. All criteria changed unfavorably and proportionally to the increased use of timed AI. The accuracy within each class of animals (cows, dams, or sires) decreased proportionally with increased use of timed AI, varying from 0.32 (S12) to 0.52 (S1) for bull EBV for CTFS. The average EBV of the top sires (best 25, 50, 75, or 100 sires) approached population average EBV values when increasing the number of treated animals. The sire rank correlation between EBV and TBV within simulated scenarios was smaller for scenarios with more synchronized animals, going from 0.38 (S12) to 0.67 (S1). The long-term use of hormonal synchronized cows clearly decreased the mean EBV over generations in the population for CTFS and FSTC. The inclusion of genetically correlated traits in a multiple trait model was effective in removing the bias due to the presence of hormonal synchronized cows. However, given the constraints within the simulation, it is important that further investigation with real data is conducted to determine the true effect of including timed AI records within genetic evaluations of fertility traits in dairy cattle.

Understanding how cows respond to heat stress has helped to provide effective herd management practices to tackle this environmental challenge. The possibility of selecting animals that are genetically more heat tolerant may provide additional means to maintain or even improve the productivity of the Canadian dairy industry, which is facing a shifting environment due to climate changes. The objective of this study was to estimate the genetic parameters for heat tolerance of milk, fat, and protein yields in Canadian Holstein cows. A total of 1.3 million test-day records from 195,448 first-parity cows were available. A repeatability test-day model fitting a reaction norm on the temperature-humidity index (THI) was used to estimate the genetic parameters. The estimated genetic correlations between additive genetic effect for production and for heat tolerance ranged from −0.13 to −0.21, indicating an antagonistic relationship between the level of production and heat tolerance. Heritability increased marginally as THI increased above its threshold for milk yield (0.20 to 0.23) and protein yield (0.14 to 0.16) and remained constant for fat yield (0.17). A Spearman rank correlation between the estimated breeding values under thermal comfort and under heat stress showed a potential genotype by environmental interaction. The existence of a genetic variability for heat tolerance allows for the selection of more heat tolerant cows.

Background: Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers. Methods: The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries. Results: A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles’ QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG “Circadian rhythm” and “Neurotrophin signaling pathway” were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones. Conclusions: The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

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.