Abstract:In genetic association analysis of complex traits, detection of interaction (either GxG or GxE) can help to elucidate the genetic architecture and biological mechanisms underlying the trait. Detection of interaction in a genome-wide interaction study (GWIS) can be methodologically challenging for various reasons,including a high burden of multiple comparisons when testing for epistasis between all possible pairs of a set of genomewide variants, as well as heteroscedasticity effects occurring in the presence of GxG or GxE interaction. In this paper, we address the problem of an even more striking phenomenon that we call the "feast or famine" effect that occurs when testing interaction in a genomewide context. We show that in any given GxE GWIS, the type 1 error of standard interaction tests performed genomewide can vary widely from the nominal level, where the actual type 1 error in any given GWIS varies as a predictable function of the observed trait and environmental values. Using standard methods, some GWISs will have systematically underinflated p-values ("feast"), and others will have systematically overinflated p-values ("famine"), which can lead to false detection of interaction, reduced power, inconsistent results across studies, and failure to replicate true signal. This startling phenomenon is specific to detection of interaction in a GWIS, and it may partly explain why such detection has often proved challenging and difficult to replicate. We show that the feast or famine effect occurs across a wide range of GxE analysis methods, including but not limited to (1) testing interaction in a linear or linear mixed model (LMM) using standard approaches such as t-tests/Wald tests, likelihood ratio tests, or score tests; (2) doing a combined interaction-association test in a linear model or LMM using standard approaches such as F-tests or likelihood ratio tests; (3) testing interaction with multiple environments or multiple SNPs, where these are modeled as random effects in a LMM using standard approaches; (4) performing tests of interaction in a GWIS where significance is assessed using permutation of the trait residuals. We show theoretically that the key cause of this phenomenon is which variables are conditioned on in the analysis, and this suggests an approach to correct the problem by changing the way the conditioning is done. Using this insight, we have developed the TINGA method to adjust the interaction test statistics to make their p-values closer to uniform under the null hypothesis. In simulations we show that TINGA both controls type 1 error and improves power. TINGA allows for covariates and population structure through use of a linear mixed model and accounts for heteroscedasticity. We apply TINGA to detection of epistasis in a study of flowering time in Arabidopsis thaliana.

The Choice of Null Distributions for Detecting Gene-Gene Interactions in Genome-Wide Association Studies

Current Status of Studies on Genome-Wide Gene-Gene Interactions

Detecting Essential and Removable Interactions in Genome-Wide Association Studies

Overcoming the “feast or famine” effect: improved interaction testing in genome-wide association studies

A FAST ALGORITHM FOR DETECTING GENE-GENE INTERACTIONS IN GENOME-WIDE ASSOCIATION STUDIES

Testing Gene-Gene Interactions Based on a Neighborhood Perspective in Genome-wide Association Studies

Testing a Large Number of Composite Null Hypotheses Using Conditionally Symmetric Multidimensional Gaussian Mixtures in Genome-Wide Studies

Measuring Gene-Gene Interaction Using Kullback-Leibler Divergence

Alternative Methods for H1 Simulations in Genome Wide Association Studies

Overcoming the "feast or famine" effect: improved interaction testing in genome-wide association studies

Two-stage Hypothesis Tests for Variable Interactions with FDR Control

A Random Forest Approach to the Detection of Epistatic Interactions in Case-Control Studies

Research on Single Nucleotide Polymorphisms Interaction Detection from Network Perspective

Practical and Theoretical Considerations in Study Design for Detecting Gene-Gene Interactions Using MDR and GMDR Approaches

A Forest-Based Approach to Identifying Gene and Gene-Gene Interactions

Current Status of SNPs Interaction in Genome-Wide Association Study

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

Detecting latent gene-environment interaction when analyzing binary traits

Quick Approximation of Threshold Values for Genome-Wide Association Studies.

Clustermi: Detecting High-Order Snp Interactions Based On Clustering And Mutual Information