Abstract:Post-GWAS analysis, in many cases, focuses on fine-mapping targeted genetic regions discovered at GWAS-stage; that is, the aim is to pinpoint potential causal variants and susceptibility genes for complex traits and disease outcomes using next-generation sequencing (NGS) technologies. Large-scale GWAS cohorts are necessary to identify target regions given the typically modest genetic effect sizes. In this context, two-phase sampling design and analysis is a cost-reduction technique that utilizes data collected during phase 1 GWAS to select an informative subsample for phase 2 sequencing. The main goal is to make inference for genetic variants measured via NGS by efficiently combining data from phases 1 and 2. We propose two approaches for selecting a phase 2 design under a budget constraint. The first method identifies sampling fractions that select a phase 2 design yielding an asymptotic variance covariance matrix with certain optimal characteristics, for example, smallest trace, via Lagrange multipliers (LM). The second relies on a genetic algorithm (GA) with a defined fitness function to identify exactly a phase 2 subsample. We perform comprehensive simulation studies to evaluate the empirical properties of the proposed designs for a genetic association study of a quantitative trait. We compare our methods against two ranked designs: residual-dependent sampling and a recently identified optimal design. Our findings demonstrate that the proposed designs, GA in particular, can render competitive power in combined phase 1 and 2 analysis compared with alternative designs while preserving type 1 error control. These results are especially evident under the more practical scenario where design values need to be defined a priori and are subject to misspecification. We illustrate the proposed methods in a study of triglyceride levels in the North Finland Birth Cohort of 1966. R code to reproduce our results is available at github.com/egosv/TwoPhase_postGWAS.

2dGBH: Two-dimensional group Benjamini–Hochberg procedure for false discovery rate control in two-way multiple testing of genomic data

EBT: a Statistic Test Identifying Moderate Size of Significant Features with Balanced Power and Precision for Genome-Wide Rate Comparisons

False discovery rate control: Moving beyond the Benjamini–Hochberg method

Directional false discovery rate control in large-scale multiple comparisons

A grouped, selectively weighted false discovery rate procedure

Asymptotic false discovery control of the Benjamini-Hochberg procedure for pairwise comparisons

Bayesian Hierarchical Hypothesis Testing in Large-Scale Genome-Wide Association Analysis

Optimal False Discovery Rate Control for Large Scale Multiple Testing with Auxiliary Information

2dFDR: a new approach to confounder adjustment substantially increases detection power in omics association studies

Two-phase sample selection strategies for design and analysis in post-genome-wide association fine-mapping studies

Exploratory data analysis for large-scale multiple testing problems and its application in gene expression studies

Optimal control of false discovery criteria in the two-group model

Genome-wide interaction-based association analysis identified multiple new susceptibility Loci for common diseases.

Two-stage Hypothesis Tests for Variable Interactions with FDR Control

Multiple Testing of Submatrices of a Precision Matrix with Applications to Identification of Between Pathway Interactions

Multiple testing in genome-wide association studies via hierarchical hidden Markov models

A more practical approach for the Benjamini-Hochberg FDR controlling procedure for huge-scale testing problems

Accurate cross-platform GWAS analysis via two-stage imputation

A generalized correlated binomial distribution with application in multiple testing problems

Controlling the Rate of GWAS False Discoveries

Simultaneous high-probability bounds on the false discovery proportion in structured, regression, and online settings