Volume 4 Supplement 1
Proceedings of the 13th European workshop on QTL mapping and marker assisted selection
Sensitivity of genomic selection to using different prior distributions
 Klara L Verbyla^{1, 2, 3, 4}Email author,
 Philip J Bowman^{2},
 Ben J Hayes^{2} and
 Michael E Goddard^{2, 3, 4}
DOI: 10.1186/175365614S1S5
© Verbyla et al; licensee BioMed Central Ltd. 2010
Published: 31 March 2010
Abstract
Genomic selection describes a selection strategy based on genomic estimated breeding values (GEBV) predicted from dense genetic markers such as single nucleotide polymorphism (SNP) data. Different Bayesian models have been suggested to derive the prediction equation, with the main difference centred around the specification of the prior distributions.
Methods
The simulated dataset of the 13^{th} QTLMAS workshop was analysed using four Bayesian approaches to predict GEBV for animals without phenotypic information. Different prior distributions were assumed to assess their affect on the accuracy of the predicted GEBV.
Conclusion
All methods produced GEBV that were highly correlated with the true breeding values. The models appear relatively insensitive to the choice of prior distributions for QTLMAS data set and this is consistent with uniformity of performance of different methods found in real data.
Background
Genomic selection describes a technique for evaluating an animal's breeding value by simultaneously evaluating and summing marker effects across the genome. It uses panels of SNPs covering the whole genome so that ideally all QTL are in linkage disequilibrium with at least one marker, thereby maximizing the proportion of genetic variance explained by the SNPs.
Meuwissen et al (2001) [1] presented three models to produce GEBV. The first invoked the infinitesimal model assumption such that all SNPs had effects derived from the same normal distribution. The other approaches used a Bayesian framework to apply hierarchical models with different prior distributions assuming unequal variances across the SNP, resulting in a t distribution for prior distribution for the QTL effects. The specification of the prior distributions of the QTL effects has been reported to be important to the accurate prediction of breeding values and when mapping multiple QTL across the entire genome [2].
The aim of this study was to assess the effect that different prior distributions and subsequently the models using these priors, had on the accuracy of estimated GEBV using the 13^{th} QTLMAS simulated data set where we had no prior knowledge of the trait's distribution of QTL effects.
Methods
Model
where y is the vector of phenotypes of the trait being analysed for all n individuals, μ is the mean, 1_{ n } is a vector of ones of length n, X_{ j } is a vector of indicator variables representing the genotypes of the j^{ th } marker for all individuals (x_{ ij }=0,1,2), β_{ j } is the size of the QTL effect associated with marker j, u is the vector of random polygenic effects of length n (Z is the associated design matrix) and is assumed to be normally distributed, u ~ N (0, A) where A is the pedigree derived additive genetic relationship matrix and e is the residual error also assumed to be normally distributed, e ~ N(0, I ) where I is the nxn identity matrix. The prior distributions for the variances of the random polygenic effects and the residual were uninformative flat priors of the form Χ^{2}( 2,0). The GEBV at each time point were calculated as
Prior distributions for SNP effects and algorithms
Prior Distribution Specifications
Method  Prior Distribution 

Bayes BLUP 

Bayes A 

Bayes A/B (Hybrid) 
with probability 1 π with probability π 
Bayes C 
γ_{ i } ~bernoulli(π) 1  p(γ_{ i } = 0) = p(γ_{ i } = 1) = π 
The other two models assumed mixture distributions for the SNP effects reflecting the assumption that there is a large number of SNPs with zero or near zero effects and a second smaller set of SNPs with larger significant effects. A Bayes A/B "hybrid" method was used. This approximation to Bayes B [1] was used to keep computational and time demands reasonable. In this algorithm, after every k Bayes A iterations, Bayes B via the reverse jump algorithm is employed. The Reverse Jump algorithm [3] is run multiple times per SNP and then any SNP with a final state of zero in the current Bayes B iterations is set to zero for the subsequent k iterations of the Bayes A. This maintains the correct transitions between models of differing dimensionality. The prior distributions are identical to that of the original Bayes B using a mixture prior distribution for the SNP variance allowing a proportion, 1π, to be set to zero. The other proportion π is sampled from the same mixture distribution as Bayes A. See Meuwissen et al (2001) for more details of priors and conditional distributions used.
A faster alternative to both the Bayes A/B hybrid and Bayes B is to use Stochastic Search Variable Selection (SSVS) [4] (Bayes C [5, 6]). This avoids the problem of the changing dimensionally of the models by providing a technique to maintain constant dimensionality across all models while still allowing the SNP in the predictive set to change. Instead of removing all nonsignificant parameters, their posterior distributions are limited to values close to zero. The major advantage of this method is that it can be implemented using the Gibbs sampler instead of the more computationally demanding algorithms such as the reverse jump algorithm. The indicator variable (γ_{ i }) determines whether the i^{th} SNP effect is sampled from the larger distribution (i.e. significant effect) or from the small distribution with near zero effects (see Table 1). The prior values of π (the proportion sampled from the nonzero distribution or the larger distribution respectively) for both Bayes A\B and Bayes C was set to 0.05, reflecting the fact that with 435 SNP, it appeared reasonable to expect at least 21 SNP would be associated with a QTL.
The algorithms associated with each model were run for 30,000 iterations with the first 10,000 discarded as burnin.
Results and Discussion
Prediction of breeding values at time point 600
The problem of how to model the time series data and estimate GEBV at time point 600 was explored. However, there was little information available to estimate any inflection points or asymptotic values. The GEBV estimated at time points 265, 397 and 530 were found to have a linear relationship (eg. appeared to form the linear part of the growth curve). Consequently, as there was no other information available after time point 530 to predict asymptotes etc., the GEBV at time point 600 were estimated by fitting a linear regression through the breeding values at the three linear time points (265, 397 and 530).
Breeding values
Correlations Between Estimated GEBV for unphenotyped animals at t=600
Bayes C  Bayes A/B  Bayes BLUP  

Bayes A  0.999  0.991  0.860 
Bayes C  1  0.993  0.863 
Bayes A/B  1  0.893 
Comparison of True and Estimated GEBV.
Method  Correlation  MSE  Rank  Regression 

Bayes.BLUP  0.885  5.479  0.691  0.979 
BayesA  0.857  7.092  0.696  1.162 
BayesA/B  0.889  5.435  0.73  1.081 
BayesC  0.861  6.561  0.71  1.024 
The inclusion of the polygenic effect in the model (not simulated in the data) only slightly reduced the accuracy of prediction (.01) but not significantly (results not shown). It was included in the model as its inclusion has been shown to produce slightly better accuracies of prediction while reducing the bias of the variance components[7].
Bayes BLUP produced a significantly different set of GEBV. This is evident by the much lower correlations with the other methods and difference in regression coefficients between BLUP and the other methods. Despite these differences Bayes BLUP produces good accuracy and a low MSE (Table 3). Hayes et al (2009) [8] reports that New Zealand, Australian, the Netherlands and United States studies all found that BLUP gave lower accuracy of GEBV than Bayesian Methods for traits where there is a single QTL that explains a large proportion of the genetic variance e.g. DGAT1 for Fat Percentage. In the current dataset a finite number of QTL were simulated where the largest amount of genetic variance explained by a single QTL was 10.5%. Despite this, Bayes BLUP is still able to produce very accurate GEBV compared to the other methods. One reason this occurs may be that a number of SNPs are required to pick up the effect of a single QTL, resulting in large numbers of SNPs with small effects, which matches the prior distribution of BLUP. However if the percentage of genetic variance explained by a single QTL was to be larger, Bayes BLUP could be expected to produce worse results. Thus this caveat to using Bayes BLUP should be considered when using this method.
Conclusion
All methods produced GEBV that were highly correlated (greater than 0.85) with the true breeding values despite diverse assumptions and prior distributions. This indicates that the hierarchical model is relatively insensitive to the choice of prior distributions for this data set. Thus all models perform well and this is consistent with the general uniformity of performance found across methods in real data. [8]. Despite the general equality in the performance of the different methods, it is still recommended that any information about a trait's QTL effect distribution and phenotypic data should be used to determine the choice of model, prior distributions and setting of the hyper parameters. This will maximise the likelihood of calculating the most accurate GEBV possible.
Declarations
Acknowledgements
KV was funded by the Marie Curie Host Fellowships for Early Stage Research Training, as part of the 6th Framework Programme of the European Commission. This Publication represents the views of the Authors, not the European Commission, and the Commission is not liable for any use that may be made of the information.
This article has been published as part of BMC Proceedings Volume 4 Supplement 1, 2009: Proceedings of 13th European workshop on QTL mapping and marker assisted selection.
The full contents of the supplement are available online at http://www.biomedcentral.com/17536561/4?issue=S1.
Authors’ Affiliations
References
 Meuwissen THE, Hayes BJ, Goddard ME: Prediction of total genetic value using genomewide dense marker maps. Genetics. 2001, 157 (4): 18191829.PubMed CentralPubMed
 Yi NJ: A unified Markov chain Monte Carlo framework for mapping multiple quantitative trait loci. Genetics. 2004, 167 (2): 967975. 10.1534/genetics.104.026286.PubMed CentralView ArticlePubMed
 Green PJ: Reversible jump Markov chain Monte Carlo computation and Bayesian model determination. Biometrika. 1995, 82 (4): 711732. 10.1093/biomet/82.4.711.View Article
 George EI, McCulloch RE: Variable Selection Via Gibbs Sampling. Journal of the American Statistical Association. 1993, 88 (423): 881889. 10.2307/2290777.View Article
 Meuwissen THE, Goddard ME: Mapping multiple QTL using linkage disequilibrium and linkage analysis information and multitrait data. Genetics Selection Evolution. 2004, 36: 261279. 10.1186/12979686363261.View Article
 Verbyla KL, Hayes BJ, Bowman PJ, Goddard ME: Accuracy of genomic selection using stochastic search variable selection in Australian Holstein Friesian dairy cattle. Genetics Research. 2009, 91 (05): 307311. 10.1017/S0016672309990243.View ArticlePubMed
 Calus MPL, Veerkamp RF: Accuracy of breeding values when using and ignoring the polygenic effect in genomic breeding value estimation with a marker density of one SNP per cM. Journal of Animal Breeding and Genetics. 2007, 124 (6): 362368.View ArticlePubMed
 Hayes BJ, Bowman PJ, Chamberlain AJ, Goddard ME: Invited review: Genomic selection in dairy cattle: Progress and challenges. J Dairy Sci. 2009, 92 (2): 433443. 10.3168/jds.20081646.View ArticlePubMed
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.