E. C. Akanno

Reliability of molecular breeding values for Warner-Bratzler shear force and carcass traits of beef cattle - an independent validation study.

Interest in genetic improvement of carcass and tenderness traits of beef cattle using genome-based selection (GS) and marker-assisted management programs is increasing. The success of such a program depends on the presence of linkage disequilibrium between the observed markers and the underlying QTL as well as on the relationship between the discovery, validation, and target populations. For molecular breeding values (MBV) predicted for a target population using SNP markers, reliabilities of these MBV can be obtained from validation analyses conducted in an independent population distinct from the discovery set. The objective of this study was to test MBV predicted for carcass and tenderness traits of beef cattle in a Canadian-based validation population that is largely independent of a United States-based discovery set. The discovery data set comprised of genotypes and phenotypes from >2,900 multibreed beef cattle while the validation population consisted of 802 crossbred feeder heifers and steers. A bivariate animal model that fitted actual phenotype and MBV was used for validation analyses. The reliability of MBV was defined as square of the genetic correlation (R(2) g) that represents the proportion of the additive genetic variance explained by the SNP markers. Several scenarios involving different starting marker panels (384, 3K, 7K, and 50K) and different sets of SNP selected to compute MBV (50, 100, 200, 375, 400, 600, and 800) were investigated. Validation results showed that the most reliable MBV (R(2) g) were 0.34 for HCW, 0.36 for back fat thickness, 0.28 for rib eye area, 0.30 for marbling score, 0.25 for yield grade, and 0.38 for Warner-Bratzler shear force across the different scenarios explored. The results indicate that smaller SNP panels can be developed for use in genetic improvement of beef carcass and tenderness traits to exploit GS benefits.

Accuracy of genomic predictions for feed efficiency traits of beef cattle using 50K and imputed HD genotypes.

The accuracy of genomic predictions can be used to assess the utility of dense marker genotypes for genetic improvement of beef efficiency traits. This study was designed to test the impact of genomic distance between training and validation populations, training population size, statistical methods, and density of genetic markers on prediction accuracy for feed efficiency traits in multibreed and crossbred beef cattle. A total of 6,794 beef cattle data collated from various projects and research herds across Canada were used. Illumina BovineSNP50 (50K) and imputed Axiom Genome-Wide BOS 1 Array (HD) genotypes were available for all animals. The traits studied were DMI, ADG, and residual feed intake (RFI). Four validation groups of 150 animals each, including Angus (AN), Charolais (CH), Angus-Hereford crosses (ANHH), and a Charolais-based composite (TX) were created by considering the genomic distance between pairs of individuals in the validation groups. Each validation group had 7 corresponding training groups of increasing sizes ( = 1,000, 1,999, 2,999, 3,999, 4,999, 5,998, and 6,644), which also represent increasing average genomic distance between pairs of individuals in the training and validations groups. Prediction of genomic estimated breeding values (GEBV) was performed using genomic best linear unbiased prediction (GBLUP) and Bayesian method C (BayesC). The accuracy of genomic predictions was defined as the Pearsons correlation between adjusted phenotype and GEBV (), unless otherwise stated. Using 50K genotypes, the highest average achieved in purebreds (AN, CH) was 0.41 for DMI, 0.34 for ADG, and 0.35 for RFI, whereas in crossbreds (ANHH, TX) it was 0.38 for DMI, 0.21 for ADG, and 0.25 for RFI. Similarly, when imputed HD genotypes were applied in purebreds (AN, CH), the highest average was 0.14 for DMI, 0.15 for ADG, and 0.14 for RFI, whereas in crossbreds (ANHH, TX) it was 0.38 for DMI, 0.22 for ADG, and 0.24 for RFI. The of GBLUP predictions were greatly reduced with increasing genomic average distance compared to those from BayesC predictions. The results indicate that 50K genotypes, used with BayesC, are more effective for predicting GEBV in purebred cattle. Imputed HD genotypes found utility when dealing with composites and crossbreds. Formulation of a fairly large training set for genomic predictions in beef cattle should consider the genomic distance between the training and target populations.

Genome-wide association for heifer reproduction and calf performance traits in beef cattle.

The aim of this study was to identify SNP markers that associate with variation in beef heifer reproduction and performance of their calves. A genome-wide association study was performed by means of the generalized quasi-likelihood score (GQLS) method using heifer genotypes from the BovineSNP50 BeadChip and estimated breeding values for pre-breeding body weight (PBW), pregnancy rate (PR), calving difficulty (CD), age at first calving (AFC), calf birth weight (BWT), calf weaning weight (WWT), and calf pre-weaning average daily gain (ADG). Data consisted of 785 replacement heifers from three Canadian research herds, namely Brandon Research Centre, Brandon, Manitoba, University of Alberta Roy Berg Kinsella Ranch, Kinsella, Alberta, and Lacombe Research Centre, Lacombe, Alberta. After applying a false discovery rate correction at a 5% significance level, a total of 4, 3, 3, 9, 6, 2, and 1 SNPs were significantly associated with PBW, PR, CD, AFC, BWT, WWT, and ADG, respectively. These SNPs were located on chromosomes 1, 5-7, 9, 13-16, 19-21, 24, 25, and 27-29. Chromosomes 1, 5, and 24 had SNPs with pleiotropic effects. New significant SNPs that impact functional traits were detected, many of which have not been previously reported. The results of this study support quantitative genetic studies related to the inheritance of these traits, and provides new knowledge regarding beef cattle quantitative trait loci effects. The identification of these SNPs provides a starting point to identify genes affecting heifer reproduction traits and performance of their calves (BWT, WWT, and ADG). They also contribute to a better understanding of the biology underlying these traits and will be potentially useful in marker- and genome-assisted selection and management.

Genomic prediction of breed composition and heterosis effects in Angus, Charolais and Hereford crosses using 50K genotypes

Abstract This study examined the feasibility and accuracy of using Illumina BovineSNP50 genotypes to estimate individual cattle breed composition and heterosis relative to estimate from pedigree. First, pedigree was used to compute breed fractions for 1124 crossbred cattle. Given the breed composition of sires and dams, retained heterosis and retained heterozygosity were computed for all individuals. Second, all animals’ genotypes were used to compute individual’s genomic breed fractions by applying a cross-validation method. Average genome-wide heterozygosity and retained heterozygosity based on genomic breed fraction were computed. Lastly, accuracies of breed composition, retained heterozygosity and retained heterosis were assessed as Pearson’s correlation between pedigree- and genome-based predictions. The average breed compositions observed were 0.52 Angus, 0.23 Charolais, and 0.25 Hereford for pedigree-based prediction and 0.46, 0.26, and 0.28 for genome-based prediction, respectively. Correlations of predicted breed composition ranged from 0.94 to 0.96. Genome-based retained heterozygosity and retained heterosis from pedigree were also highly correlated (0.96). A positive association of nonadditive genetic effects was observed for growth traits reflecting the importance of heterosis for these traits. Genomic prediction can aid analyses that depend on knowledge of breed composition and serve as a reliable method to predict heterosis to improve the efficiency of commercial crossbreeding schemes.

Development and validation of a small SNP panel for feed efficiency in beef cattle

The objective of this study was to develop and validate a customized cost-effective single nucleotide polymorphism (SNP) panel for genetic improvement of feed efficiency in beef cattle. The SNPs identified in previous association studies and through extensive analysis of candidate genomic regions and genes, were screened for their functional impact and allele frequency in Angus and Hereford breeds used as validation candidates for the panel. Association analyses were performed on genotypes of 159 SNPs from new samples of Angus (n = 160), Hereford (n = 329), and Angus-Hereford crossbred (n = 382) cattle using allele substitution and genotypic models in ASReml. Genomic heritabilities were estimated for feed efficiency traits using the full set of SNPs, SNPs associated with at least one of the traits (at P ≤ 0.05 and P < 0.10), as well as the Illumina bovine 50K representing a widely used commercial genotyping panel. A total of 63 SNPs within 43 genes showed association (P ≤ 0.05) with at least one trait. The minor alleles of SNPs located in the GHR and CAST genes were associated with decreasing effects on residual feed intake (RFI) and/or RFI adjusted for backfat (RFIf), whereas minor alleles of SNPs within MKI67 gene were associated with increasing effects on RFI and RFIf. Additionally, the minor allele of rs137400016 SNP within CNTFR was associated with increasing average daily gain (ADG). The SNPs genotypes within UMPS, SMARCAL, CCSER1, and LMCD1 genes showed significant over-dominance effects whereas other SNPs located in SMARCAL1, ANXA2, CACNA1G, and PHYHIPL genes showed additive effects on RFI and RFIf. Gene enrichment analysis indicated that gland development, as well as ion and cation transport are important physiological mechanisms contributing to variation in feed efficiency traits. The study revealed the effect of the Jak-STAT signaling pathway on feed efficiency through the CNTFR, OSMR, and GHR genes. Genomic heritability using the 63 significant (P ≤ 0.05) SNPs was 0.09, 0.09, 0.13, 0.05, 0.05, and 0.07 for ADG, dry matter intake, midpoint metabolic weight, RFI, RFIf, and backfat, respectively. These SNPs contributed to genetic variation in the studied traits and thus can potentially be used or tested to generate cost-effective molecular breeding values for feed efficiency in beef cattle.

Modeling heterotic effects in beef cattle using genome-wide SNP-marker genotypes1

An objective of commercial beef cattle crossbreeding programs is to simultaneously optimize use of additive (breed differences) and non-additive (heterosis) effects. A total of 6,794 multibreed and crossbred beef cattle with phenotype and Illumina BovineSNP50 genotype data were used to predict genomic heterosis for growth and carcass traits by applying two methods assumed to be linearly proportional to heterosis. The methods were as follows: 1) retained heterozygosity predicted from genomic breed fractions (HET1) and 2) deviation of adjusted crossbred phenotype from midparent value (HET2). Comparison of methods was based on prediction accuracy from cross-validation. Here, a mutually exclusive random sampling of all crossbred animals (n = 5,327) was performed to form five groups replicated five times with approximately 1,065 animals per group. In each run within a replicate, one group was assigned as a validation set, while the remaining four groups were combined to form the reference set. The phenotype of the animals in the validation set was assumed to be unknown; thus, it resulted in every animal having heterosis values that were predicted without using its own phenotype, allowing their adjusted phenotype to be used for validation. The same approach was used to test the impact of predicted heterosis on accuracy of genomic breeding values (GBV). The results showed positive heterotic effects for growth traits but not for carcass traits that reflect the importance of heterosis for growth traits in beef cattle. Heterosis predicted by HET1 method resulted in less variable estimates that were mostly within the range of estimates generated by HET2. Prediction accuracy was greater for HET2 (0.37-0.98) than HET1 (0.34-0.43). Proper consideration of heterosis in genomic evaluation models has debatable effects on accuracy of EBV predictions. However, opportunity exists for predicting heterosis, improving accuracy of genomic selection, and consequently optimizing crossbreeding programs in beef cattle.

Genomic retained heterosis effects on fertility and lifetime productivity in beef heifers

Abstract: This study evaluated the effects of three genomic indicators of heterosis on female fertility and lifetime productivity, and quantified changes over 11 production cycles in a crossbred cow herd. Pedigree-based breed composition (pBC) was determined and used to calculate retained heterozygosity for 412 replacement heifers born from 2004 to 2014 at the Lacombe Research and Development Centre (AB, Canada). Heifers were followed as cows over 1050 mating opportunities, 11 production cycles, and five parities. Heifers and their sires (51) were genotyped and these genotypes were used to predict each animal’s genomic breed composition (gBC) and three genomic indicators of heterosis: (1) retained heterozygosity (RHETg), (2) heterozygous proportion (H), and (3) retained heterosis (RHg). Correlations between pedigree and genomic breed fractions for Angus, Hereford, and Charolais were high (r p = 0.74–0.94; P < 0.001). Genomic indicators of heterosis were highly related (r p = 0.61 for RHETg vs. H; 0.71 for RHg vs. H; 0.96 for RHETg vs. RHg; P < 0.001). Each 10% change in RHETg resulted in 51 ± 20 d longer survival (P = 0.011) in the herd and 35.7 ± 15.2 kg more (P = 0.019) calf wean weight per cow exposed to breeding when summed over five parities. These differences resulted in an extra

0310 Assessing genetic diversity in Canadian beef cattle populations using Illumina BovineSNP50 chip.

161 per heifer in a year. Optimizing heterosis using genomic tools can be very beneficial for the cow herd if applied correctly.