Abstract:Genome-wide association studies (GWAS) are commonly employed to study the genetic basis of complex traits/diseases, and a key question is how much heritability could be explained by all single nucleotide polymorphisms (SNPs) in GWAS. One widely used approach that relies on summary statistics only is linkage disequilibrium score regression (LDSC); however, this approach requires certain assumptions about the effects of SNPs (e.g., all SNPs contribute to heritability and each SNP contributes equal variance). More flexible modeling methods may be useful. We previously developed an approach recovering the "true" effect sizes from a set of observed z-statistics with an empirical Bayes approach, using only summary statistics. However, methods for standard error (SE) estimation are not available yet, limiting the interpretation of our results and the applicability of the approach. In this study, we developed several resampling-based approaches to estimate the SE of SNP-based heritability, including two jackknife and three parametric bootstrap methods. The resampling procedures are performed at the SNP level as it is most common to estimate heritability from GWAS summary statistics alone. Simulations showed that the delete-d-jackknife and parametric bootstrap approaches provide good estimates of the SE. In particular, the parametric bootstrap approaches yield the lowest root-mean-squared-error (RMSE) of the true SE. We also explored various methods for constructing confidence intervals (CIs). In addition, we applied our method to estimate the SNP-based heritability of 12 immune-related traits (levels of cytokines and growth factors) to shed light on their genetic architecture. We also implemented the methods to compute the sum of heritability explained and the corresponding SE in an R package SumVg. In conclusion, SumVg may provide a useful alternative tool for calculating SNP heritability and estimating SE/CI, which does not rely on distributional assumptions of SNP effects.

Subsampling Technique to Estimate Variance Component for UK-Biobank Traits

Analytical and computational solution for the estimation of SNP-heritability in biobank-scale and distributed datasets

UK BioCoin: Swift Trait-Specific Summary Statistics Regression for UK Biobank

Dissecting the High-Resolution Genetic Architecture of Complex Phenotypes by Accurately Estimating Gene-Based Conditional Heritability.

Total Heritability Explained by All Variants in Genome-Wide Association Studies Based on Summary Statistics with Standard Error Estimates

Scalable summary statistics-based heritability estimation method with individual genotype level accuracy

Genome-wide discovery for biomarkers using quantile regression at biobank scale

Reliable Heritability Estimation Using Sparse Regularization in Ultrahigh Dimensional Genome-Wide Association Studies

Bias in Estimates of Quantitative-Trait–Locus Effect in Genome Scans: Demonstration of the Phenomenon and a Method-of-Moments Procedure for Reducing Bias

Association Test Between Haplotypes and Longitudinal Traits in Complex Pedigrees.

Estimating trans-ancestry genetic correlation with unbalanced data resources

SumVg: Total heritability explained by all variants in genome-wide association studies based on summary statistics with standard error estimates

Probabilistic inference of the genetic architecture underlying functional enrichment of complex traits

Quantitative Trait Linkage Analysis by Generalized Estimating Equations: Unification of Variance Components and Haseman‐Elston Regression

Accurate and Efficient Estimation of Local Heritability using Summary Statistics and LD Matrix

Optimal Estimation of Co-heritability in High-dimensional Linear Models

Will Big Data Close the Missing Heritability Gap?

Estimating Genetic Variance Contributed by a Quantitative Trait Locus: Removing Nuisance Parameters

Subset scanning for multi-trait analysis using GWAS summary statistics

Leveraging LD eigenvalue regression to improve the estimation of SNP heritability and confounding inflation

A Scalable Adaptive Quadratic Kernel Method for Interpretable Epistasis Analysis in Complex Traits