Abstract

The objective of this work was to evaluate the use of regularized quantile regression (RQR) for genomic prediction analyses in traits with or without skewness and with different proportions of epistatic variance. Data were simulated for 2500 individuals. The genome included 5K markers and 50 QTL for a lowly heritable trait (heritability = 0.1) simulated for six scenarios: combinations of trait distributions (symmetric normal or positive skewed) and percentage of epistatic and additive genetic variances (100% epistatic, 50% epistatic and 50% additive, 100% additive). Only additive-by-additive effects were simulated for epistasis. Three replicates were generated per scenario. Estimation of marker effects was performed using 2000 individuals (training). The remaining 500 individuals were used for validation. Accuracy was calculated as the correlation between genomic estimated breeding values and true breeding values. Analyses were performed using Bayesian LASSO (BLASSO) and RQR considering three quantiles (0.15 [RQR0.15], 0.50[RQR0.50], 0.85[RQR0.85]) in R. Results are shown in Table 032. Overall, accuracies of genomic prediction for the normal symmetric trait were variable. There was no relationship between RQR quantiles and accuracies of genomic prediction. In general, accuracies of genomic prediction decreased as the proportion of epistatic variance increased. The highest accuracies of genomic prediction were 0.15 (BLASSO), 0.17 (BLASSO), and 0.20 (RQR0.15) for the proportions of epistatic variance of 100%, 50%, and 0%, respectively. For the positively skewed trait, the greatest accuracies were obtained for RQR0.15, with 0.16, 0.26, and 0.27, followed by RQR0.50, with 0.12, 0.18, and 0.21 for the proportions of epistatic variance of 100%, 50%, and 0%, respectively. Accuracies for RQR0.15 were 77%, 53%, and 50% greater than those obtained using BLASSO. For RQR, accuracies decreased as quantiles increased. Within each method, accuracies were similar for proportions of epistatic variances of 50% and 0%. In conclusion, for the symmetric normal trait, BLASSO had, overall, slightly greater accuracies than for RQR. In contrast, RQR presented greater accuracies of genomic prediction compared to BLASSO when quantiles 0.15 and 0.50 were used for the trait with positive skewness. These results suggest that improved accuracy of genomic prediction for skewed traits can be obtained using RQR.

Table 032.

Accuracies of genomic prediction.

Percentage of epistatic /additive genetic variance
Trait distributionMethodology100% / 0%50% / 50%0% / 100%
Symmetric normalBlasso0.150.170.19
RQR0.150.100.160.13
RQR0.500.140.130.20
RQR0.850.110.160.16
Positive skewnessBLASSO0.090.170.18
RQR0.150.160.260.27
RQR0.500.120.180.21
RQR0.850.050.120.12
Percentage of epistatic /additive genetic variance
Trait distributionMethodology100% / 0%50% / 50%0% / 100%
Symmetric normalBlasso0.150.170.19
RQR0.150.100.160.13
RQR0.500.140.130.20
RQR0.850.110.160.16
Positive skewnessBLASSO0.090.170.18
RQR0.150.160.260.27
RQR0.500.120.180.21
RQR0.850.050.120.12
Table 032.

Accuracies of genomic prediction.

Percentage of epistatic /additive genetic variance
Trait distributionMethodology100% / 0%50% / 50%0% / 100%
Symmetric normalBlasso0.150.170.19
RQR0.150.100.160.13
RQR0.500.140.130.20
RQR0.850.110.160.16
Positive skewnessBLASSO0.090.170.18
RQR0.150.160.260.27
RQR0.500.120.180.21
RQR0.850.050.120.12
Percentage of epistatic /additive genetic variance
Trait distributionMethodology100% / 0%50% / 50%0% / 100%
Symmetric normalBlasso0.150.170.19
RQR0.150.100.160.13
RQR0.500.140.130.20
RQR0.850.110.160.16
Positive skewnessBLASSO0.090.170.18
RQR0.150.160.260.27
RQR0.500.120.180.21
RQR0.850.050.120.12
This content is only available as a PDF.
You do not currently have access to this article.