Genomic selection (GS) is a genetic breeding method that uses genome-wide marker information to improve the accuracy of the prediction of complex traits. The single-step GBLUP (ssGBLUP) model, which integrates pedigree, phenotypic, and genomic data, has improved genomic prediction. However, ssGBLUP assumes that all markers contribute equally to genetic variance, which can limit its predictive accuracy, especially for traits controlled by major genes. To overcome this limitation, we integrate results from genome-wide association studies (GWAS) into an enhanced ssGBLUP framework, termed single-step genome-wide association assisted BLUP (ssGWABLUP). Our approach assigns differential weights to markers on the basis of their GWAS results, thereby increasing the contribution of effective markers while diminishing the influence of ineffective ones during the construction of the genomic relationship matrix. By incorporating pseudo quantitative trait nucleotides (pQTNs) as covariates, we aim to capture the effects of markers closely associated with major causal variants, leading to the development of the ssGWABLUP_pQTNs. Compared with weighted ssGBLUP (WssGBLUP), the ssGWABLUP model demonstrated superior accuracy and dispersion across different genetic architectures. We then compared the performance of our proposed ssGWABLUP_pQTNs model against both ssGBLUP and ssGWABLUP across various genetic scenarios. Our results demonstrate that ssGWABLUP_pQTNs outperforms other models in terms of prediction accuracy, particularly in scenarios with simpler genetic architectures. Additionally, evaluation using pig dataset confirmed the effectiveness of ssGWABLUP_pQTNs, highlighting its potential for practical breeding applications. The incorporation of pQTNs and a weighted genomic relationship matrix presents a promising and potentially scalable approach to further enhance genomic prediction, with potential implications for improving the accuracy of genomic selection in breeding programs.

Natural or artificial selection could shape genetic architecture, e.g., the relationship between minor allele frequency (MAF) and the effect sizes of causal variants (CVs). This study aimed to investigate the impact of the MAF–effect size relationship (as a selection signature, S) on genomic prediction and heritability estimation in livestock, using both simulated data (Holstein) and real datasets (Holstein and pigs). We evaluated the performance of two models: (1) selection-adjusted genomic best linear unbiased prediction (GBLUP-S), and (2) MAF-stratified selection-adjusted genomic best linear unbiased prediction (GBLUP-SMS). Simulation results demonstrated that for traits under strong negative selection (S < −1), both GBLUP-S and GBLUP-SMS outperformed classic GBLUP. The prediction accuracy of GBLUP-S improved by 0.011–0.031, while GBLUP-SMS achieved a gain of 0.005–0.025. Furthermore, GBLUP-SMS exhibited lower sensitivity to variations in S-values, whereas GBLUP-S heavily relied on accurate S specification. When the true S was matched, GBLUP-SMS generated more unbiased (or comparable) heritability estimates and higher prediction accuracy relative to GBLUP-S. Critically, mismatched S in GBLUP-S led to increased bias in heritability estimates and reduced prediction accuracy. Cross-validation with real phenotypic data from Holsteins and pigs demonstrated that implementing selection-adjust methods improved prediction accuracy by 0.015 for FP in Holsteins and 0.01 for T1 in pigs, while enhancing the unbiasedness of heritability estimates across all traits. Negative selection signatures were identified for cattle (S = −0.5) and pig T1, T2, and T3 (S = −1.5, −1, and −2, respectively). These findings advance the theoretical framework of GBLUP-based genomic prediction and heritability estimation.