Comparison of predictive ability of single-trait and multitrait genomic selection models for body growth traits in Maiwa yaks

Yaks are grazed extensively on the Qinghai–Tibet Plateau, which has a long history of semi-domestication. The predicted weight of yaks over consecutive years helps make strategic decisions when selecting yak calves for breeding. To achieve more accurate predictions of genomic estimated breeding values, we used a dataset comprising the genotype and weight records of 396 Maiwa yaks collected from 2015 to 2020. We compared the predictive accuracy of the genome best linear unbiased prediction (GBLUP) model with that of six other models. Based on the GBLUP model, we applied two prediction strategies. In the first strategy, the year was treated as a fixed effect in the GBLUP model, and the kinship from all individuals and the markers were treated as random effects. In the second strategy, all individuals were divided into six age groups, with GBLUP performed for each group, and the phenotypes of the closest age groups were treated as fixed effects. Although the GBLUP model provided better prediction accuracy than other single-trait models, most of the predictive capacity was derived from the best linear unbiased estimation. Additionally, incorporating the phenotype of the closest age group as a factor in multitrait prediction enhanced the accuracy of the model. Our findings provide a robust and credible strategy for predicting continuous economic traits in the presence of strong correlations.