Efficient genomic prediction based on whole-genome sequence data using split-and-merge Bayesian variable selection

Genetics Selection Evolution

Research Article

Aniek C. Bouwman¹  Mario P. L. Calus¹  Roel F. Veerkamp¹  Chris Schrooten² 

[1] Animal Breeding and Genomics Centre, Wageningen UR Livestock Research, PO Box 338, 6700 AH, Wageningen, The Netherlands;CRV BV, 6800 AL, Arnhem, The Netherlands;

关键词: Quantitative Trait Locus Region;    Effective Sample Size;    Genomic Prediction;    Genomic Estimate Breeding Value;    Somatic Cell Score;   

DOI  :  10.1186/s12711-016-0225-x

 received in 2016-02-17, accepted in 2016-06-16,  发布年份 2016

来源: Springer

【 摘 要 】

BackgroundUse of whole-genome sequence data is expected to increase persistency of genomic prediction across generations and breeds but affects model performance and requires increased computing time. In this study, we investigated whether the split-and-merge Bayesian stochastic search variable selection (BSSVS) model could overcome these issues. BSSVS is performed first on subsets of sequence-based variants and then on a merged dataset containing variants selected in the first step.ResultsWe used a dataset that included 4,154,064 variants after editing and de-regressed proofs for 3415 reference and 2138 validation bulls for somatic cell score, protein yield and interval first to last insemination. In the first step, BSSVS was performed on 106 subsets each containing ~39,189 variants. In the second step, 1060 up to 472,492 variants, selected from the first step, were included to estimate the accuracy of genomic prediction. Accuracies were at best equal to those achieved with the commonly used Bovine 50k-SNP chip, although the number of variants within a few well-known quantitative trait loci regions was considerably enriched. When variant selection and the final genomic prediction were performed on the same data, predictions were biased. Predictions computed as the average of the predictions computed for each subset achieved the highest accuracies, i.e. 0.5 to 1.1 % higher than the accuracies obtained with the 50k-SNP chip, and yielded the least biased predictions. Finally, the accuracy of genomic predictions obtained when all sequence-based variants were included was similar or up to 1.4 % lower compared to that based on the average predictions across the subsets. By applying parallelization, the split-and-merge procedure was completed in 5 days, while the standard analysis including all sequence-based variants took more than three months.ConclusionsThe split-and-merge approach splits one large computational task into many much smaller ones, which allows the use of parallel processing and thus efficient genomic prediction based on whole-genome sequence data. The split-and-merge approach did not improve prediction accuracy, probably because we used data on a single breed for which relationships between individuals were high. Nevertheless, the split-and-merge approach may have potential for applications on data from multiple breeds.

【 授权许可】

CC BY   
© The Author(s) 2016

【 图 表 】

Fig. 3

Fig. 1

40517_2023_273_Article_IEq9.gif

Fig. 2

40517_2023_273_Article_IEq6.gif

40517_2023_273_Article_IEq4.gif

Fig. 1

40517_2023_273_Article_IEq2.gif

Fig. 2

42004_2023_1025_Article_IEq7.gif

Fig. 9

Fig. 9

Fig. 4

Fig. 1

Fig. 4

Fig. 3

Fig. 4

41512_2023_158_Article_IEq26.gif

Fig. 2

Fig. 2

12951_2017_270_Article_IEq3.gif

Fig. 1

41512_2023_158_Article_IEq20.gif

Fig. 1

Fig. 8

Fig. 6

Fig. 2

Fig. 2

Fig. 6

Fig. 1

Fig. 6

12951_2015_155_Article_IEq6.gif

41512_2023_158_Article_IEq9.gif

Fig. 8

12951_2017_323_Article_IEq1.gif

Fig. 5

Fig. 3

41512_2023_158_Article_IEq2.gif

Fig. 7

41512_2023_158_Article_IEq1.gif

Fig. 1

Fig. 5

12936_2017_2051_Article_IEq86.gif

12936_2017_1932_Article_IEq15.gif

12936_2017_2051_Article_IEq85.gif

Fig. 1

Fig. 1

Fig. 7

Fig. 4

Fig. 3

Fig. 3

Fig. 6

Fig. 1

12951_2015_155_Article_IEq4.gif

12936_2017_2014_Article_IEq78.gif

12867_2016_60_Article_IEq1.gif

Fig. 4

Fig. 2

Fig. 9

12867_2016_60_Article_IEq2.gif

Fig. 1

Fig. 1

Fig. 2

Fig. 10

Fig. 1

Fig. 1

Fig. 1

Fig. 7

12951_2015_155_Article_IEq2.gif

Fig. 1

Fig. 9

12936_2017_2051_Article_IEq78.gif

Fig. 5

Fig. 2

Fig. 8