Robust Tests for Additive Gene-Environment Interaction in Case-Control Studies Using Gene-Environment Independence

An Evaluation of Gene-Diet Interaction Statistical Methods and Discovery of rs7175421-Whole Grain Interaction in Lung Cancer

Dietary factors show different effects on genetically diverse populations. Scientific research uses gene-environment interaction models to study the effects of dietary factors on genetically diverse populations for lung cancer risk. However, previous study designs have not investigated the degree of type I error inflation and, in some instances, have not corrected for multiple testing. Using a motivating investigation of diet-gene interaction and lung cancer risk, we propose a training and testing strategy and perform real-world simulations to select the appropriate statistical methods to reduce false-positive discoveries. The simulation results show that the unconstrained maximum likelihood (UML) method controls the type I error better than the constrained maximum likelihood (CML). The empirical Bayesian (EB) method can compete with the UML method in achieving statistical power and controlling type I error. We observed a significant interaction between SNP rs7175421 with dietary whole grain in lung cancer prevention, with an effect size (standard error) of -0.312 (0.112) for EB estimate. SNP rs7175421 may interact with dietary whole grains in modulating lung cancer risk. Evaluating statistical methods for gene-diet interaction analysis can help balance the statistical power and type I error.

A Robust Test for Additive Gene-Environment Interaction Under the Trend Effect of Genotype Using an Empirical Bayes-Type Shrinkage Estimator

Evaluating gene by environment (G × E) interaction under an additive risk model (i.e., additive interaction) has gained wider attention. Recently, statistical tests have been proposed for detecting additive interaction, utilizing an assumption on gene-environment (G-E) independence to boost power, that do not rely on restrictive genetic models such as dominant or recessive models. However, a major limitation of these methods is a sharp increase in type I error when this assumption is violated. Our goal was to develop a robust test for additive G × E interaction under the trend effect of genotype, applying an empirical Bayes-type shrinkage estimator of the relative excess risk due to interaction. The proposed method uses a set of constraints to impose the trend effect of genotype and builds an estimator that data-adaptively shrinks an estimator of relative excess risk due to interaction obtained under a general model for G-E dependence using a retrospective likelihood framework. Numerical study under varying levels of departures from G-E independence shows that the proposed method is robust against the violation of the independence assumption while providing an adequate balance between bias and efficiency compared with existing methods. We applied the proposed method to the genetic data of Alzheimer disease and lung cancer.

A Likelihood Ratio Test for Gene-Environment Interaction Based on the Trend Effect of Genotype Under an Additive Risk Model Using the Gene-Environment Independence Assumption

Several statistical methods have been proposed for testing gene-environment (G-E) interactions under additive risk models using data from genome-wide association studies. However, these approaches have strong assumptions from underlying genetic models, such as dominant or recessive effects that are known to be less robust when the true genetic model is unknown. We aimed to develop a robust trend test employing a likelihood ratio test for detecting G-E interaction under an additive risk model, while incorporating the G-E independence assumption to increase power. We used a constrained likelihood to impose 2 sets of constraints for: 1) the linear trend effect of genotype and 2) the additive joint effects of gene and environment. To incorporate the G-E independence assumption, a retrospective likelihood was used versus a standard prospective likelihood. Numerical investigation suggests that the proposed tests are more powerful than tests assuming dominant, recessive, or general models under various parameter settings and under both likelihoods. Incorporation of the independence assumption enhances efficiency by 2.5-fold. We applied the proposed methods to examine the gene-smoking interaction for lung cancer and gene-apolipoprotein E $\varepsilon$4 interaction for Alzheimer disease, which identified 2 interactions between apolipoprotein E $\varepsilon$4 and loci membrane-spanning 4-domains subfamily A (MS4A) and bridging integrator 1 (BIN1) genes at genome-wide significance that were replicated using independent data.

Impact of Gene-Environment Interactions on Cancer Development

Several epidemiological and experimental studies have demonstrated that many human diseases are not only caused by specific genetic and environmental factors but also by gene-environment interactions. Although it has been widely reported that genetic polymorphisms play a critical role in human susceptibility to cancer and other chronic disease conditions, many single nucleotide polymorphisms (SNPs) are caused by somatic mutations resulting from human exposure to environmental stressors. Scientific evidence suggests that the etiology of many chronic illnesses is caused by the joint effect between genetics and the environment. Research has also pointed out that the interactions of environmental factors with specific allelic variants highly modulate the susceptibility to diseases. Hence, many scientific discoveries on gene-environment interactions have elucidated the impact of their combined effect on the incidence and/or prevalence rate of human diseases. In this review, we provide an overview of the nature of gene-environment interactions, and discuss their role in human cancers, with special emphases on lung, colorectal, bladder, breast, ovarian, and prostate cancers.

Estimating Additive Interaction Effect in Stratified Two-Phase Case-Control Design

BACKGROUND AND AIMS

There is considerable interest in epidemiology to estimate an additive interaction effect between two risk factors in case-control studies. An additive interaction is defined as the differential reduction in absolute risk associated with one factor between different levels of the other factor. A stratified two-phase case-control design is commonly used in epidemiology to reduce the cost of assembling covariates. It is crucial to obtain valid estimates of the model parameters by accounting for the underlying stratification scheme to obtain accurate and precise estimates of additive interaction effects. The aim of this paper is to examine the properties of different methods for estimating model parameters and additive interaction effects under a stratified two-phase case-control design.

METHODS

Using simulations, we investigate the properties of three existing methods, namely stratum-specific offset, inverse-probability weighting, and multiple imputation for estimating model parameters and additive interaction effects. We also illustrate these properties using data from two published epidemiology studies.

RESULTS

Simulation studies show that the multiple imputation method performs well when both the true and analysis models are additive (i.e., does not include multiplicative interaction terms) but does not provide a discernible advantage over the offset method when the analysis models are non-additive (i.e., includes multiplicative interaction terms). The offset method exhibits the best overall properties when the analysis model contains multiplicative interaction effects.

CONCLUSION

When estimating additive interaction between risk factors in stratified two-phase case-control studies, we recommend estimating model parameters using multiple imputation when the analysis model is additive, and we recommend the offset method when the analysis model is non-additive.

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

It is increasingly of interest in statistical genetics to test for the presence of an additive interaction between genetic (G) and environmental (E) risk factors. In case-control studies involving a rare disease, a statistical test of no additive G×E interaction typically entails a test of no relative excess risk due to interaction (RERI). It has been shown that a likelihood ratio test of a null RERI incorporating the G-E independence assumption (RERI-LRT) outperforms the standard approach. The RERI-LRT relies on correct specification of a logistic model for the binary outcome, as a function of G, E, and auxiliary covariates. However, when at least one exposure is not categorical or auxiliary covariates are present, nonparametric estimation may not be feasible, while parametric logistic regression will a priori rule out the null hypothesis of no additive interaction in most practical situations, inflating type I error rate. In this paper, we present a general approach to test for G × E additive interaction exploiting G-E independence. Unlike the RERI-LRT, it allows the regression model for the binary outcome to remain unrestricted, and nonetheless still allows for covariate adjustment in order to ensure the G-E independence assumption or to rule out residual confounding. The methods are illustrated through extensive simulation studies and an ovarian cancer study.

A comprehensive gene-environment interaction analysis in Ovarian Cancer using genome-wide significant common variants

As a follow-up to genome-wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well-known ovarian cancer risk factors and their interactions with 28 genome-wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case-control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10-4 ). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46-0.60) compared to 0.71 (95%CI = 0.66-0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1-5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene-environment analysis in ovarian cancer.