Abstract:While genome-wide association studies (GWAS) have discovered thousands of risk loci for heritable disorders, so far even very large meta-analyses have recovered only a fraction of the heritability of most complex traits. Recent work utilizing variance components models has demonstrated that a larger fraction of the heritability of complex phenotypes is captured by the additive effects of SNPs than is evident only in loci surpassing genome-wide significance thresholds, typically set at a Bonferroni-inspired p ≤ 5 x 10 -8 . Procedures that control false discovery rate can be more powerful, yet these are still under-powered to detect the majority of non-null effects from GWAS. The current work proposes a novel Bayesian semi-parametric two-group mixture model and develops a Markov Chain Monte Carlo (MCMC) algorithm for a covariate-modulated local false discovery rate ( cmfdr ). The probability of being non-null depends on a set of covariates via a logistic function, and the non-null distribution is approximated as a linear combination of B-spline densities, where the weight of each B-spline density depends on a multinomial function of the covariates. The proposed methods were motivated by work on a large meta-analysis of schizophrenia GWAS performed by the Psychiatric Genetics Consortium (PGC). We show that the new cmfdr model fits the PGC schizophrenia GWAS test statistics well, performing better than our previously proposed parametric gamma model for estimating the non-null density and substantially improving power over usual fdr. Using loci declared significant at cmfdr ≤ 0.20, we perform follow-up pathway analyses using the Kyoto Encyclopedia of Genes and Genomes (KEGG) homo sapiens pathways database. We demonstrate that the increased yield from the cmfdr model results in an improved ability to test for pathways associated with schizophrenia compared to using those SNPs selected according to usual fdr.

Distance-based analysis of variance: approximate inference and an application to genome-wide association studies

Shrinkage-based regularization tests for high-dimensional data with application to gene set analysis

A Kernel approach for extending nonparametric multivariate analysis of variance in high-dimensional settings

D-MANOVA: fast distance-based multivariate analysis of variance for large-scale microbiome association studies

Equivalence of Kernel Machine Regression and Kernel Distance Covariance for Multidimensional Trait Association Studies

Multiple Comparison Procedures for Neuroimaging Genomewide Association Studies

Novel probabilistic models of spatial genetic ancestry with applications to stratification correction in genome-wide association studies

Finding the basic neighborhood in variable range Markov random fields: application in SNP association studies

Fast and Accurate Approximation to Significance Tests in Genome-Wide Association Studies

SEMI-PARAMETRIC COVARIATE-MODULATED LOCAL FALSE DISCOVERY RATE FOR GENOME-WIDE ASSOCIATION STUDIES

Distance-Based Statistical Inference

MANOCCA: a robust and computationally efficient test of covariance in high-dimension multivariate omics data

ADAPTIVE-WEIGHT BURDEN TEST FOR ASSOCIATIONS BETWEEN QUANTITATIVE TRAITS AND GENOTYPE DATA WITH COMPLEX CORRELATIONS.

Statistical Inference for Genetic Relatedness Based on High-Dimensional Logistic Regression

Genetic association tests in family samples for multi-category phenotypes

High-dimensional MANOVA via Bootstrapping and its Application to Functional and Sparse Count Data

A Bayes Factor Approach with Informative Prior for Rare Genetic Variant Analysis from Next Generation Sequencing Data

High-dimensional multivariate analysis of variance via geometric median and bootstrapping

Variance Components Genetic Association Test for Zero-inflated Count Outcomes

Combination rules for homoscedastic and heteroscedastic MANOVA models from multiply imputed datasets

Analysis of distributional variation through multi-scale Beta-Binomial modeling