Genomic prediction with whole-genome sequence data in intensely selected pig lines

BackgroundEarly simulations indicated that whole-genome sequence data (WGS) could improve genomic prediction accuracy and its persistence across generations and breeds. However, empirical results have been ambiguous so far. Large data sets that capture most of the genome diversity in a population must be assembled so that allele substitution effects are estimated with high accuracy. The objectives of this study were to use a large pig dataset to assess the benefits of using WGS for genomic prediction compared to using commercial marker arrays, to identify scenarios in which WGS provides the largest advantage, and to identify potential pitfalls for its effective implementation. MethodsWe sequenced 6,931 individuals from seven commercial pig lines with different numerical size. Genotypes of 32.8 million variants were imputed for 396,100 individuals (17,224 to 104,661 per line). We used BayesR to perform genomic prediction for eight complex traits. Genomic predictions were performed using either data from a marker array or variants preselected from WGS based on association tests. ResultsThe prediction accuracy with each set of preselected WGS variants was not robust across traits and lines and the improvements in prediction accuracy that we achieved so far with WGS compared to marker arrays were generally small. The most favourable results for WGS were obtained when the largest training sets were available and used to preselect variants with statistically significant associations to the trait for augmenting the established marker array. With this method and training sets of around 80k individuals, average improvements of genomic prediction accuracy of 0.025 were observed in within-line scenarios. ConclusionsOur results showed that WGS has a small potential to improve genomic prediction accuracy compared to marker arrays in intensely selected pig lines in some settings. Thus, although we expect that more robust improvements could be attained with a combination of larger training sets and optimised pipelines, the use of WGS in the current implementations of genomic prediction should be carefully evaluated on a case-by-case basis against the cost of generating WGS at a large scale.