Parameters affecting genome simulation for evaluating genomic selection method

Genomic best linear unbiased prediction method including imprinting effects for genomic evaluation

BACKGROUND

Genomic best linear unbiased prediction (GBLUP) is a statistical method used to predict breeding values using single nucleotide polymorphisms for selection in animal and plant breeding. Genetic effects are often modeled as additively acting marker allele effects. However, the actual mode of biological action can differ from this assumption. Many livestock traits exhibit genomic imprinting, which may substantially contribute to the total genetic variation of quantitative traits. Here, we present two statistical models of GBLUP including imprinting effects (GBLUP-I) on the basis of genotypic values (GBLUP-I1) and gametic values (GBLUP-I2). The performance of these models for the estimation of variance components and prediction of genetic values across a range of genetic variations was evaluated in simulations.

RESULTS

Estimates of total genetic variances and residual variances with GBLUP-I1 and GBLUP-I2 were close to the true values and the regression coefficients of total genetic values on their estimates were close to 1. Accuracies of estimated total genetic values in both GBLUP-I methods increased with increasing degree of imprinting and broad-sense heritability. When the imprinting variances were equal to 1.4% to 6.0% of the phenotypic variances, the accuracies of estimated total genetic values with GBLUP-I1 exceeded those with GBLUP by 1.4% to 7.8%. In comparison with GBLUP-I1, the superiority of GBLUP-I2 over GBLUP depended strongly on degree of imprinting and difference in genetic values between paternal and maternal alleles. When paternal and maternal alleles were predicted (phasing accuracy was equal to 0.979), accuracies of the estimated total genetic values in GBLUP-I1 and GBLUP-I2 were 1.7% and 1.2% lower than when paternal and maternal alleles were known.

CONCLUSIONS

This simulation study shows that GBLUP-I1 and GBLUP-I2 can accurately estimate total genetic variance and perform well for the prediction of total genetic values. GBLUP-I1 is preferred for genomic evaluation, while GBLUP-I2 is preferred when the imprinting effects are large, and the genetic effects differ substantially between sexes.

Genomic best linear unbiased prediction method reflecting the degree of linkage disequilibrium

The degree of linkage disequilibrium (LD) between markers differs depending on the location of the genome; this difference biases genetic evaluation by genomic best linear unbiased prediction (GBLUP). To correct this bias, we used three GBLUP methods reflecting the degree of LD (GBLUP-LD). In the three GBLUP-LD methods, genomic relationship matrices were conducted from single nucleotide polymorphism markers weighted according to local LD levels. The predictive abilities of GBLUP-LD were investigated by estimating variance components and assessing the accuracies of estimated breeding values using simulation data. When quantitative trait loci (QTL) were located at weak LD regions, the predictive abilities of the three GBLUP-LD methods were superior to those of GBLUP and Bayesian lasso except when the number of QTL was small. In particular, the superiority of GBLUP-LD increased with decreasing trait heritability. The rates of QTL at weak LD regions would increase when selection by GBLUP continues; this consequently decreases the predictive ability of GBLUP. Thus, the GBLUP-LD could be applicable for populations selected by GBLUP for a long time. However, if QTL were located at strong LD regions, the accuracies of three GBLUP-LD methods were lower than GBLUP and Bayesian lasso.