Zhang et al. tackle the O(nΒ³) computational burden of mixed linear models (MLMs) in GWAS with two complementary methods: compression, which clusters individuals into s kinship-based groups (UPGMA clustering; SAS proc cluster) so the random effect dimension shrinks, giving a cΒ³ speedup where c = n/s; and P3D, which estimates variance components once from a reduced model and fixes them as empirical-Bayes priors for all marker tests, eliminating per-marker iteration .
Evidence spans three species: human GOLDN (n=1,315; 637 markers), dog (n=292; 1,000 of 23,500 SNPs), and maize (n=277; 553 SNPs). P3D P-values correlated rΒ² > 0.96 with full optimization across observed phenotypes and rΒ² > 0.97 across simulated architectures (hΒ² = 0.25β0.5, 5β20 QTNs, dominance d β {0,0.25,0.5,1}, epistasis 0β20% variance; 1,000 replicates) .
Reported values from the human dataset (821 s vs 0.34 s per marker with compression level 3.8; projected 26 years β 2.7 days for 1M markers).
Counterintuitively, moderate compression (levels 1.5β10) improved both model fit and power versus the standard MLM β up to 34%, 42%, and 20% for human, dog, and maize β because group-level kinship buffers noise in individual kinship estimates. QQ plots showed better type-I error control than both GLM and standard MLM at these levels. The reduced-model fit (marker-free) tracked full-model power, so optimal compression can be chosen cheaply once per study .
Downstream literature (e.g., multi-locus MLM extensions) built directly on this framework, and the method remains a foundational, widely implemented GWAS acceleration approach .
Confidence: high for reported within-study benchmarks (exact quantitative data, 1,000-replicate simulations); moderate for generalization to modern biobank-scale data, which the paper did not test. What would change the verdict: independent failure of compression/P3D to preserve power on large diverse cohorts or with imputed-sequence kinship.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.