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.

Detecting latent gene-environment interaction when analyzing binary traits

A Varying Coefficient Model to Jointly Test Genetic and Gene–Environment Interaction Effects

Powerful cocktail methods for detecting genome-wide gene-environment interaction.

Distinct explanations underlie gene-environment interactions in the UK Biobank

Detecting gene–environment interactions from multiple continuous traits

A multivariate approach to joint testing of main genetic and gene-environment interaction effects

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits

Measuring Gene-Gene Interaction Using Kullback-Leibler Divergence

Identification of Gene-Environment Interactions with Marginal Penalization

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

A unified method for rare variant analysis of gene-environment interactions

A general approach to detect gene (G)-environment (E) additive interaction leveraging G-E independence in case-control studies

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

A penalized robust semiparametric approach for gene-environment interactions

A fast algorithm for detecting gene-gene interactions in genome-wide association studies

Structured gene-environment interaction analysis

Discovering non-additive heritability using additive GWAS summary statistics

Detecting Essential and Removable Interactions in Genome-Wide Association Studies

Moderating heritability with genomic data

High-Dimensional Gene–Environment Interaction Analysis