Detecting latent gene-environment interaction when analyzing binary traits

Abstract

AbstractIn genome-wide association studies (GWAS), it is desirable to test for interactions (GxE) between single-nucleotide polymorphisms (SNPs,G’s) and environmental variables (E’s). However, directly accounting for interaction is often infeasible, becauseEis latent. For quantitative traits (Y) that are approximately normally distributed, it has been shown that indirect testing onGxEcan be done by testing for heteroskedasticity ofYbetween genotypes. However, when traits are binary, the existing methodology based on testing the heteroskedasticity of the trait across genotypes cannot be generalized. In this paper, we propose an approach to indirectly testGxEfor binary traits based on the non-additive effectG, and subsequently propose a joint test that accounts for the main and interaction effects of each SNP during GWAS. We illustrate the statistical features including type-I-error control and power of the proposed method through extensive numerical studies. Applying our method to the UK Biobank dataset, we showcase the practical utility of the proposed method, revealing SNPs and genes with strong potential for latent interaction effects.