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

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 Pearson's 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.

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.

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.