Abstract:Abstract Motivation Genetics hold great promise to precision medicine by tailoring treatment to the individual patient based on their genetic profiles. Toward this goal, many large-scale genome-wide association studies (GWAS) have been performed in the last decade to identify genetic variants associated with various traits and diseases. They have successfully identified tens of thousands of disease-related variants. However they have explained only a small proportion of the overall trait heritability for most traits and are of very limited clinical use. This is partly owing to the small effect sizes of most genetic variants, and the common practice of testing association between one trait and one genetic variant at a time in most GWAS, even when multiple related traits are often measured for each individual. Increasing evidence suggests that many genetic variants can influence multiple traits simultaneously, and we can gain more power by testing association of multiple traits simultaneously. It is appealing to develop novel multi-trait association test methods that need only GWAS summary data, since it is generally very hard to access the individual-level GWAS phenotype and genotype data. Results Many existing GWAS summary data-based association test methods have relied on ad hoc approach or crude Monte Carlo approximation. In this article, we develop rigorous statistical methods for efficient and powerful multi-trait association test. We develop robust and efficient methods to accurately estimate the marginal trait correlation matrix using only GWAS summary data. We construct the principal component (PC)-based association test from the summary statistics. PC-based test has optimal power when the underlying multi-trait signal can be captured by the first PC, and otherwise it will have suboptimal performance. We develop an adaptive test by optimally weighting the PC-based test and the omnibus chi-square test to achieve robust performance under various scenarios. We develop efficient numerical algorithms to compute the analytical P-values for all the proposed tests without the need of Monte Carlo sampling. We illustrate the utility of proposed methods through application to the GWAS meta-analysis summary data for multiple lipids and glycemic traits. We identify multiple novel loci that were missed by individual trait-based association test. Availability and implementation All the proposed methods are implemented in an R package available at http://www.github.com/baolinwu/MTAR. The developed R programs are extremely efficient: it takes less than 2 min to compute the list of genome-wide significant single nucleotide polymorphisms (SNPs) for all proposed multi-trait tests for the lipids GWAS summary data with 2.5 million SNPs on a single Linux desktop. Supplementary information Supplementary data are available at Bioinformatics online.

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

Regression-based Approach for Testing the Association Between Multi-Region Haplotype Configuration and Complex Trait

LARGE-SCALE MULTIVARIATE SPARSE REGRESSION WITH APPLICATIONS TO UK BIOBANK

Mapping the genetic architecture of complex traits in experimental populations.

Large-Scale Sparse Regression for Multiple Responses with Applications to UK Biobank

A Fast Multi-Locus Ridge Regression Algorithm for High-Dimensional Genome-Wide Association Studies

A Non-Parametric Method for Building Predictive Genetic Tests on High-Dimensional Data

Identifying Quantitative Trait Loci Via Group-Sparse Multitask Regression and Feature Selection: an Imaging Genetics Study of the Adni Cohort

LLR: a Latent Low-Rank Approach to Colocalizing Genetic Risk Variants in Multiple GWAS.

Group Sparse Reduced Rank Regression for Neuroimaging Genetic Study.

Integrative functional linear model for genome-wide association studies with multiple traits

Accounting for non-genetic factors by low-rank representation and sparse regression for eQTL mapping.

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

Identifying rare variants using a Bayesian regression approach

Nonconvex Penalty Based Low-Rank Representation and Sparse Regression for Eqtl Mapping

Modeling and testing for joint association using a genetic random field model

Integrate multiple traits to detect novel trait–gene association using GWAS summary data with an adaptive test approach

Bayesian Mapping of Quantitative Trait Loci for Multiple Complex Traits with the Use of Variance Components

Family-Based Association Tests for Genomewide Association Scans

Simultaneous Analysis of Common and Rare Variants in Complex Traits: Application to SNPs (Scarvasnp)

Sparse Regression Models for Unraveling Group and Individual Associations in Eqtl Mapping