Novel likelihood ratio tests for screening gene-gene and gene-environment interactions with unbalanced...

There has been extensive literature on modeling gene-gene interaction (GGI) and gene-environment interaction (GEI) in case-control studies with limited literature on statistical methods for GGI and GEI in longitudinal cohort studies. We borrow ideas from the classical two-way analysis of variance literature to address the issue of robust modeling of interactions in repeated-measures studies. While classical interaction models proposed by Tukey and Mandel have interaction structures as a function of main effects, a newer class of models, additive main effects and multiplicative interaction (AMMI) models, do not have similar restrictive assumptions on the interaction structure. AMMI entails a singular value decomposition of the cell residual matrix after fitting the additive main effects and has been shown to perform well across various interaction structures. We consider these models for testing GGI and GEI from two perspectives: likelihood ratio test based on cell means and a regression-based approach using individual observations. Simulation results indicate that both approaches for AMMI models lead to valid tests in terms of maintaining the type I error rate, with the regression approach having better power properties. The performance of these models was evaluated across different interaction structures and 12 common epistasis patterns. In summary, AMMI model is robust with respect to misspecified interaction structure and is a useful screening tool for interaction even in the absence of main effects. We use the proposed methods to examine the interplay between the hemochromatosis gene and cumulative lead exposure on pulse pressure in the Normative Aging Study.