Data from: Bayesian methods for estimating GEBVs of threshold traits

Estimation of genomic breeding values is the key step in genomic selection. Many methods have been proposed for continuous traits, but methods for threshold traits are still scarce. Here we introduced threshold model to the framework of genomic selection, and specifically we extended the three Bayesian methods BayesA, BayesB and BayesCπ based on threshold model for estimating genomic breeding values of threshold traits, and the extended methods are correspondingly termed BayesTA, BayesTB and BayesTCπ. Computing procedures of the three BayesT methods using Markov Chain Monte Carlo (MCMC) algorithm were derived. A simulation study was performed to investigate the benefit of the presented methods in accuracy of genomic estimated breeding values (GEBVs) for threshold traits. Factors affecting the performance of the three BayesT methods were addressed. As expected, the three BayesT methods generally performed better than the corresponding normal Bayesian methods, in particular when the number of phenotypic categories was small. In the standard scenario (No. categories = 2, incidence = 30%, No. QTL = 50, h2 = 0.3), the accuracies were improved by 30.4, 2.4, and 5.7 percentage points, respectively. In most scenarios, BayesTB and BayesTCπ generated similar accuracies and both performed better than BayesTA. In conclusion, our work proved that threshold model fits well for predicting GEBVs of threshold traits, and BayesTCπ is supposed to be the method of choice for genomic selection of threshold traits.基因组育种值估计是基因组选择的核心步骤。目前已针对连续性性状提出了多种分析方法，但针对阈性状的相关方法仍较为匮乏。本研究将阈模型引入基因组选择框架，针对阈性状的基因组育种值估计任务，对三种经典贝叶斯方法（BayesA、BayesB及BayesCπ）进行了基于阈模型的拓展，拓展后的方法分别命名为BayesTA、BayesTB及BayesTCπ。本研究推导了这三种BayesT方法基于马尔可夫链蒙特卡洛（Markov Chain Monte Carlo, MCMC）算法的计算流程，并通过模拟试验评估了所提方法对阈性状基因组估计育种值（genomic estimated breeding values, GEBVs）估计精度的提升效果，同时探讨了影响三种BayesT方法性能的相关因素。正如预期，三种BayesT方法的整体表现均优于对应的常规贝叶斯方法，尤其在表型分类数较少时优势更为显著。在标准试验场景（表型分类数=2，发病率=30%，数量性状基因座（Quantitative Trait Locus, QTL）数量=50，遗传力h²=0.3）中，三种方法的估计精度分别提升了30.4、2.4及5.7个百分点。在多数试验场景中，BayesTB与BayesTCπ的估计精度较为相近，且二者的表现均优于BayesTA。综上，本研究证实阈模型可有效用于阈性状的基因组估计育种值预测，且BayesTCπ可作为阈性状基因组选择的首选方法。