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.

GWASBrewer: An R Package for Simulating Realistic GWAS Summary Statistics

GWASBrewer: An R Package for Simulating Realistic GWAS Summary Statistics

GAUSS: a summary-statistics-based R package for accurate estimation of linkage disequilibrium for variants, Gaussian imputation, and TWAS analysis of cosmopolitan cohorts

Small-group originating model: Optimized individual-level GWAS simulation featured by SLiM and using open-access data

Simulating genetic risk scores from summary statistics

Multivariate simulation framework reveals performance of multi-trait GWAS methods

CCAFE: Estimating Case and Control Allele Frequencies from GWAS Summary Statistics

WebGWAS: A web server for instant GWAS on arbitrary phenotypes

twas_sim, a Python-based tool for simulation and power analysis of transcriptome-wide association analysis

Simple models of non-random mating and environmental transmission bias standard human genetics statistical methods

A statistical framework for GWAS of high dimensional phenotypes using summary statistics, with application to metabolite GWAS

MungeSumstats: a Bioconductor package for the standardization and quality control of many GWAS summary statistics

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

pgainsim: an R-package to assess the mode of inheritance for quantitative trait loci in GWAS

The goldmine of GWAS summary statistics: a systematic review of methods and tools

A review: simulation tools for genome-wide interaction studies

A statistical simulation model to guide the choices of analytical methods in arrayed CRISPR screen experiments

Alternative Methods for H1 Simulations in Genome Wide Association Studies

Using simulation studies to evaluate statistical methods

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

simmrd: An open‐source tool to perform simulations in Mendelian randomization