Nonparametric Reduced-Rank Regression for Multi-SNP, Multi-Trait Association Mapping

Abstract

Genome-wide association studies have proven to be essential for understanding the genetic basis of disease. However, many complex traits---personality traits, facial features, disease subtyping---are inherently high-dimensional, impeding simple approaches to association mapping. We developed a nonparametric Bayesian reduced rank regression model for multi-SNP, multi-trait association mapping that does not require the rank of the linear subspace to be specified. We show in simulations and real data that our model shares strength over SNPs and over correlated traits, improving statistical power to identify genetic associations with an interpretable, SNP-supervised low-dimensional linear projection of the high-dimensional phenotype. On the HapMap phase 3 gene expression QTL study data, we identify pleiotropic expression QTLs that classical univariate tests are underpowered to find and that two step approaches cannot recover. Our Python software, BERRRI, is publicly available at GitHub: https://github.com/ashlee1031/BERRRI.

References in corpus (9)

• Pathwise coordinate optimization
• Inferring Parameters and Structure of Latent Variable Models by Variational Bayes
• A statistical framework for joint eQTL analysis in multiple tissues
• Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs
• A latent factor model with a mixture of sparse and dense factors to model gene expression data with confounding effects
• Nonparametric Reduced Rank Regression
• Bayesian Structured Sparsity from Gaussian Fields
• Differential gene co-expression networks via Bayesian biclustering models
• Bayesian group latent factor analysis with structured sparsity