Iterative hard thresholding in genome-wide association studies: Generalized linear models, prior weights, and double sparsity
Open Access
- 1 June 2020
- journal article
- research article
- Published by Oxford University Press (OUP) in GigaScience
- Vol. 9 (6)
Abstract
Background: Consecutive testing of single nucleotide polymorphisms (SNPs) is usually employed to identify genetic variants associated with complex traits. Ideally one should model all covariates in unison, but most existing analysis methods for genome-wide association studies (GWAS) perform only univariate regression. Results: We extend and efficiently implement iterative hard thresholding (IHT) for multiple regression, treating all SNPs simultaneously. Our extensions accommodate generalized linear models, prior information on genetic variants, and grouping of variants. In our simulations, IHT recovers up to 30% more true predictors than SNP-by-SNP association testing and exhibits a 2–3 orders of magnitude decrease in false-positive rates compared with lasso regression. We also test IHT on the UK Biobank hypertension phenotypes and the Northern Finland Birth Cohort of 1966 cardiovascular phenotypes. We find that IHT scales to the large datasets of contemporary human genetics and recovers the plausible genetic variants identified by previous studies. Conclusions: Our real data analysis and simulation studies suggest that IHT can (i) recover highly correlated predictors, (ii) avoid over-fitting, (iii) deliver better true-positive and false-positive rates than either marginal testing or lasso regression, (iv) recover unbiased regression coefficients, (v) exploit prior information and group-sparsity, and (vi) be used with biobank-sized datasets. Although these advances are studied for genome-wide association studies inference, our extensions are pertinent to other regression problems with large numbers of predictors.Keywords
All Related Versions
Funding Information
- National Human Genome Research Institute (HG006139, HG009120)
- National Institute of General Medical Sciences (GM053275, GM053275)
- National Science Foundation (DMS-1264153)
This publication has 44 references indexed in Scilit:
- PUMA: A Unified Framework for Penalized Multiple Regression Analysis of GWAS DataPLOS Computational Biology, 2013
- Mendel: the Swiss army knife of genetic analysis programsBioinformatics, 2013
- Chapter 11: Genome-Wide Association StudiesPLOS Computational Biology, 2012
- Random-Effects Model Aimed at Discovering Associations in Meta-Analysis of Genome-wide Association StudiesThe American Journal of Human Genetics, 2011
- Association screening of common and rare genetic variants by penalized regressionBioinformatics, 2010
- Stability SelectionJournal of the Royal Statistical Society Series B: Statistical Methodology, 2010
- Prioritizing GWAS Results: A Review of Statistical Methods and Recommendations for Their ApplicationThe American Journal of Human Genetics, 2010
- Genome-wide association analysis by lasso penalized logistic regressionBioinformatics, 2009
- Genome-wide association analysis of metabolic traits in a birth cohort from a founder populationNature Genetics, 2008
- Identification of a Novel Risk Locus for Progressive Supranuclear Palsy by a Pooled Genomewide Scan of 500,288 Single-Nucleotide PolymorphismsThe American Journal of Human Genetics, 2007