BackgroundGenomic selection is a recently developed technology that is beginning to revolutionize animal breeding. The objective of this study was to estimate marker effects to derive prediction equations for direct genomic values for 16 routinely recorded traits of American Angus beef cattle and quantify corresponding accuracies of prediction.MethodsDeregressed estimated breeding values were used as observations in a weighted analysis to derive direct genomic values for 3570 sires genotyped using the Illumina BovineSNP50 BeadChip. These bulls were clustered into five groups using K-means clustering on pedigree estimates of additive genetic relationships between animals, with the aim of increasing within-group and decreasing between-group relationships. All five combinations of four groups were used for model training, with cross-validation performed in the group not used in training. Bivariate animal models were used for each trait to estimate the genetic correlation between deregressed estimated breeding values and direct genomic values.ResultsAccuracies of direct genomic values ranged from 0.22 to 0.69 for the studied traits, with an average of 0.44. Predictions were more accurate when animals within the validation group were more closely related to animals in the training set. When training and validation sets were formed by random allocation, the accuracies of direct genomic values ranged from 0.38 to 0.85, with an average of 0.65, reflecting the greater relationship between animals in training and validation. The accuracies of direct genomic values obtained from training on older animals and validating in younger animals were intermediate to the accuracies obtained from K-means clustering and random clustering for most traits. The genetic correlation between deregressed estimated breeding values and direct genomic values ranged from 0.15 to 0.80 for the traits studied.ConclusionsThese results suggest that genomic estimates of genetic merit can be produced in beef cattle at a young age but the recurrent inclusion of genotyped sires in retraining analyses will be necessary to routinely produce for the industry the direct genomic values with the highest accuracy.

The impact of 9 single nucleotide polymorphisms (SNP) in the leptin (LEP), leptin receptor (LEPR), growth hormone receptor (GHR), and diacylglycerol acyltransferase (DGAT1) gene loci on daily milk production, feed intake, and feed conversion, and weekly measures of live weight, BCS, and body energy traits was evaluated using genetic and phenotypic data on 571 Holstein cows raised at the Langhill Dairy Cattle Research Center in Scotland. Six SNP were typed on the LEP gene and 1 on each of the other 3 loci. Of the 6 LEP SNP, 3 were in very high linkage disequilibrium, meaning there is little gain in typing all of them in the future. Seven LEP haplotypes were identified by parsimony-based analyses. Random-regression allelesubstitution models were used to assess the impact of each SNP allele or haplotype on the traits of interest. Diacylglycerol acyltransferase had a significant effect on milk yield, whereas GHR significantly affected feed intake, feed conversion, and body energy traits. There was also evidence of dominance in allelic effects on milk yield and BCS. The LEP haplotype CCGTTT (corresponding to leptin SNP C207T, C528T, A1457G, C963T, A252T, and C305T, respectively) significantly affected milk yield and feed and dry matter intake. Animals carrying this haplotype produced 3.13 kg more milk daily and consumed 4.64 kg more feed. Furthermore, they tended to preserve more energy than average. Such results may be used to facilitate genetic selection in animal breeding programs.

BackgroundSeveral methods have recently been developed to identify regions of the genome that have been exposed to strong selection. However, recent theoretical and empirical work suggests that polygenic models are required to identify the genomic regions that are more moderately responding to ongoing selection on complex traits. We examine the effects of multi-trait selection on the genome of a population of US registered Angus beef cattle born over a 50-year period representing approximately 10 generations of selection. We present results from the application of a quantitative genetic model, called Birth Date Selection Mapping, to identify signatures of recent ongoing selection.ResultsWe show that US Angus cattle have been systematically selected to alter their mean additive genetic merit for most of the 16 production traits routinely recorded by breeders. Using Birth Date Selection Mapping, we estimate the time-dependency of allele frequency for 44,817 SNP loci using genomic best linear unbiased prediction, generalized least squares, and BayesCπ analyses. Finally, we reconstruct the primary phenotypes that have historically been exposed to selection from a genome-wide analysis of the 16 production traits and gene ontology enrichment analysis.ConclusionsWe demonstrate that Birth Date Selection Mapping utilizing mixed models corrects for time-dependent pedigree sampling effects that lead to spurious SNP associations and reveals genomic signatures of ongoing selection on complex traits. Because multiple traits have historically been selected in concert and most quantitative trait loci have small effects, selection has incrementally altered allele frequencies throughout the genome. Two quantitative trait loci of large effect were not the most strongly selected of the loci due to their antagonistic pleiotropic effects on strongly selected phenotypes. Birth Date Selection Mapping may readily be extended to temporally-stratified human or model organism populations.

The objectives were to estimate genetic parameters needed to elucidate the relationships of a molecular breeding value (MBV) for marbling, intramuscular fat (IMF) of yearling bulls measured with ultrasound, and marbling score (MRB) of slaughtered steers, and to assess the utility of MBV and IMF in predicting the breeding value for MRB. Records for MRB (n = 38,296) and IMF (n = 6,594) were from the American Angus Association database used for national cattle evaluation. A total of 1,006 records of MBV were used in this study. (Co)variance components were estimated with ASREML, fitting an animal model with fixed contemporary groups for MRB and IMF similar to those used in the Angus national genetic evaluation. The overall mean was the only fixed effect included in the model for MBV. Heritability estimates for carcass measures were 0.48 +/- 0.03, 0.31 +/- 0.03, and 0.98 +/- 0.05 for MRB, IMF, and MBV, respectively. Genetic correlations of IMF and MBV with MRB were 0.56 +/- 0.09 and 0.38 +/- 0.10, respectively. The genetic correlation between IMF and MBV was 0.80 +/- 0.22. These results indicate the MBV evaluated may yield a greater genetic advance of approximately 20% when used as an indicator trait for genetic prediction of MRB compared with IMF. However, neither of these indicators alone provides sufficient information to produce highly accurate prediction of breeding value for the economically relevant trait MRB. Given that the goal is a highly accurate prediction of true breeding value for MRB, results of this work point to the need to 1) continue progeny testing, and 2) continue increasing the genetic correlation between the MBV and MRB.

Summary Imputation of moderate-density genotypes from low-density panels is of increasing interest in genomic selection, because it can dramatically reduce genotyping costs. Several imputation software packages have been developed, but they vary in imputation accuracy, and imputed genotypes may be inconsistent among methods. An AdaBoost-like approach is proposed to combine imputation results from several independent software packages, i.e. Beagle(v3.3), IMPUTE(v2.0), fastPHASE(v1.4), AlphaImpute, findhap(v2) and Fimpute(v2), with each package serving as a basic classifier in an ensemble-based system. The ensemble-based method computes weights sequentially for all classifiers, and combines results from component methods via weighted majority 'voting' to determine unknown genotypes. The data included 3078 registered Angus cattle, each genotyped with the Illumina BovineSNP50 BeadChip. SNP genotypes on three chromosomes (BTA1, BTA16 and BTA28) were used to compare imputation accuracy among methods, and the application involved the imputation of 50K genotypes covering 29 chromosomes based on a set of 5K genotypes. Beagle and Fimpute had the greatest accuracy among the six imputation packages, which ranged from 0·8677 to 0·9858. The proposed ensemble method was better than any of these packages, but the sequence of independent classifiers in the voting scheme affected imputation accuracy. The ensemble systems yielding the best imputation accuracies were those that had Beagle as first classifier, followed by one or two methods that utilized pedigree information. A salient feature of the proposed ensemble method is that it can solve imputation inconsistencies among different imputation methods, hence leading to a more reliable system for imputing genotypes relative to independent methods.

BackgroundArtificial neural networks (ANN) mimic the function of the human brain and are capable of performing massively parallel computations for data processing and knowledge representation. ANN can capture nonlinear relationships between predictors and responses and can adaptively learn complex functional forms, in particular, for situations where conventional regression models are ineffective. In a previous study, ANN with Bayesian regularization outperformed a benchmark linear model when predicting milk yield in dairy cattle or grain yield of wheat. Although breeding values rely on the assumption of additive inheritance, the predictive capabilities of ANN are of interest from the perspective of their potential to increase the accuracy of prediction of molecular breeding values used for genomic selection. This motivated the present study, in which the aim was to investigate the accuracy of ANN when predicting the expected progeny difference (EPD) of marbling score in Angus cattle. Various ANN architectures were explored, which involved two training algorithms, two types of activation functions, and from 1 to 4 neurons in hidden layers. For comparison, BayesCπ models were used to select a subset of optimal markers (referred to as feature selection), under the assumption of additive inheritance, and then the marker effects were estimated using BayesCπ with π set equal to zero. This procedure is referred to as BayesCpC and was implemented on a high-throughput computing cluster.ResultsThe ANN with Bayesian regularization method performed equally well for prediction of EPD as BayesCpC, based on prediction accuracy and sum of squared errors. With the 3K-SNP panel, for example, prediction accuracy was 0.776 using BayesCpC, and ranged from 0.776 to 0.807 using BRANN. With the selected 700-SNP panel, prediction accuracy was 0.863 for BayesCpC and ranged from 0.842 to 0.858 for BRANN. However, prediction accuracy for the ANN with scaled conjugate gradient back-propagation was lower, ranging from 0.653 to 0.689 with the 3K-SNP panel, and from 0.743 to 0.793 with the selected 700-SNP panel.ConclusionsANN with Bayesian regularization performed as well as linear Bayesian regression models in predicting additive genetic values, supporting the idea that ANN are useful as universal approximators of functions of interest in breeding contexts.

No abstract

The association between two single nucleotide polymorphisms (SNPs), T945M and UCP1SNP1, with hot carcass weight (HCW, kg, N = 618), longissimus dorsi muscle area (REA, cm(2), N = 633), and backfat thickness (BF, mm, N = 625), measured in Nellore cattle in Brazil, was evaluated. Likelihood ratio tests were used to evaluate reduced (fixed effects of general mean, contemporary group, yearling weight, age at slaughter, and random effect of infinitesimal genetic value) and full model (reduced model effects plus quantitative trait locus effects). Additive and dominance effects were tested for each SNP. Genotypic and gene frequencies were also obtained for the SNPs and a descriptive phenotype analysis was made. Mean values for HCW, REA and BF were equal to 288.13 +/- 0.55 kg, 73.14 +/- 0.27 cm(2), and 4.28 +/- 0.07 mm, respectively; the coefficients of variation were 4.74, 9.24, and 42.43%, respectively. Gene frequencies for T945M and UCP1SNP1 were f(C) = 0.89, f(T) = 0.11, f(C) = 0.81, and f(G) = 0.19. The SNP T945M had a genotypic frequency of only three animals for TT genotype. Additive effects were observed for T945M on REA and BF, while UCP1SNP1 affected HCW and BF. Based on the significant additive effects of the SNPs and the gene frequencies that we found, we can expect genetic gains with marker assisted selection.