Application of multiple shrinkage methods to genomic predictions1

Abstract

New challenges have arisen with the development of large marker panels for livestock species. Models easily become overparameterized when all available markers are included. Solutions have led to the development of shrinkage or regularization techniques. The objective of this study was the application and comparison of Bayesian LASSO (B-L), thick-tailed (Student-t), and semiparametric multiple shrinkage methods. The B-L and Student-t methods were also each analyzed within a single shrinkage and a multiple shrinkage framework. Simulated and real data were used to evaluate each method's performance. Real data consisted of SNP genotypes of 4,069 Holstein sires. Traits included in analysis of real data were milk, fat, protein yield, and somatic cell score. The performance of each model was compared based on correlations between true and predicted genomic predicted transmitting abilities. Model performance was also compared with the performance of routinely used methods such as Bayes-A and GBLUP through cross-validation techniques. When using simulated data regardless of shrinkage framework, shrinkage models outperformed genomic BLUP (GBLUP). The average advantage of shrinkage models ranged from 1% to approximately 8% depending on the prior specification. When analyzing real data, shrinkage models slightly outperformed GBLUP for most traits. Shrinkage models were better able to model traits for which 1 or more SNP of large effect have been identified. Overall, results suggested a relatively small advantage in multiple shrinkage models. Multiple shrinkage methods could represent a useful alternative to current methods of prediction; however, their performance in a variety of scenarios needs to be investigated further.