Fast parallelized sampling of Bayesian regression models for whole-genome prediction