Simultaneously testing for marginal genetic association and gene-environment interaction

A versatile, fast and unbiased method for estimation of gene-by-environment interaction effects on biobank-scale datasets

Identification of gene-by-environment interactions (GxE) is crucial to understand the interplay of environmental effects on complex traits. However, current methods evaluating GxE on biobank-scale datasets have limitations. We introduce MonsterLM, a multiple linear regression method that does not rely on model specification and provides unbiased estimates of variance explained by GxE. We demonstrate robustness of MonsterLM through comprehensive genome-wide simulations using real genetic data from 325,989 individuals. We estimate GxE using waist-to-hip-ratio, smoking, and exercise as the environmental variables on 13 outcomes (N = 297,529-325,989) in the UK Biobank. GxE variance is significant for 8 environment-outcome pairs, ranging from 0.009 - 0.071. The majority of GxE variance involves SNPs without strong marginal or interaction associations. We observe modest improvements in polygenic score prediction when incorporating GxE. Our results imply a significant contribution of GxE to complex trait variance and we show MonsterLM to be well-purposed to handle this with biobank-scale data.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

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

Most human traits are influenced by the interplay between genetic and environmental factors. Many statistical methods have been proposed to screen for gene-environment interaction (GxE) in the post genome-wide association study era. However, most of the existing methods assume a linear interaction between genetic and environmental factors toward phenotypic variations, which diminishes statistical power in the case of nonlinear GxE. In this paper, we present a flexible statistical procedure to detect GxE regardless of whether the underlying relationship is linear or not. By modeling the joint genetic and GxE effects as a varying-coefficient function of the environmental factor, the proposed model is able to capture dynamic trajectories of GxE. We employ a likelihood ratio test with a fast Monte Carlo algorithm for hypothesis testing. Simulations were conducted to evaluate validity and power of the proposed model in various settings. Real data analysis was performed to illustrate its power, in particular, in the case of nonlinear GxE.

Early adversity as the prototype gene × environment interaction in mental disorders?

Childhood adversity (CA) as a significant stressor has consistently been associated with the development of mental disorders. The interaction between CA and genetic variants has been proposed to play a substantial role in disease etiology. In this review, we focus on the gene by environment (GxE) paradigm, its background and interpretation and stress the necessity of its implementation in psychiatric research. Further, we discuss the findings supporting GxCA interactions, ranging from candidate gene studies to polygenic and genome-wide approaches, their strengths and limitations. To illustrate potential underlying epigenetic mechanisms by which GxE effects are translated, we focus on results from FKBP5 × CA studies and discuss how molecular evidence can supplement previous GxE findings. In conclusion, while GxE studies constitute a valuable line of investigation, more harmonized GxE studies in large, deep-phenotyped, longitudinal cohorts, and across different developmental stages are necessary to further substantiate and understand reported GxE findings.

A two-step approach to testing overall effect of gene-environment interaction for multiple phenotypes

MOTIVATION

While gene-environment (GxE) interactions contribute importantly to many different phenotypes, detecting such interactions requires well-powered studies and has proven difficult. To address this, we combine two approaches to improve GxE power: simultaneously evaluating multiple phenotypes and using a two-step analysis approach. Previous work shows that the power to identify a main genetic effect can be improved by simultaneously analyzing multiple related phenotypes. For a univariate phenotype, two-step methods produce higher power for detecting a GxE interaction compared to single step analysis. Therefore, we propose a two-step approach to test for an overall GxE effect for multiple phenotypes.

RESULTS

Using simulations we demonstrate that, when more than one phenotype has GxE effect (i.e., GxE pleiotropy), our approach offers substantial gain in power (18%-43%) to detect an aggregate-level GxE effect for a multivariate phenotype compared to an analogous two-step method to identify GxE effect for a univariate phenotype. We applied the proposed approach to simultaneously analyze three lipids, LDL, HDL and Triglyceride with the frequency of alcohol consumption as environmental factor in the UK Biobank. The method identified two loci with an overall GxE effect on the vector of lipids, one of which was missed by the competing approaches.

AVAILABILITY

We provide an R package MPGE implementing the proposed approach which is available from CRAN: https://cran.r-project.org/web/packages/MPGE/index.html.

Association of Body Mass Index With Colorectal Cancer Risk by Genome-Wide Variants

BACKGROUND

Body mass index (BMI) is a complex phenotype that may interact with genetic variants to influence colorectal cancer risk.

METHODS

We tested multiplicative statistical interactions between BMI (per 5 kg/m2) and approximately 2.7 million single nucleotide polymorphisms with colorectal cancer risk among 14 059 colorectal cancer case (53.2% women) and 14 416 control (53.8% women) participants. All analyses were stratified by sex a priori. Statistical methods included 2-step (ie, Cocktail method) and single-step (ie, case-control logistic regression and a joint 2-degree of freedom test) procedures. All statistical tests were two-sided.

RESULTS

Each 5 kg/m2 increase in BMI was associated with higher risks of colorectal cancer, less so for women (odds ratio [OR] = 1.14, 95% confidence intervals [CI] = 1.11 to 1.18; P = 9.75 × 10-17) than for men (OR = 1.26, 95% CI = 1.20 to 1.32; P = 2.13 × 10-24). The 2-step Cocktail method identified an interaction for women, but not men, between BMI and a SMAD7 intronic variant at 18q21.1 (rs4939827; Pobserved = .0009; Pthreshold = .005). A joint 2-degree of freedom test was consistent with this finding for women (joint P = 2.43 × 10-10). Each 5 kg/m2 increase in BMI was more strongly associated with colorectal cancer risk for women with the rs4939827-CC genotype (OR = 1.24, 95% CI = 1.16 to 1.32; P = 2.60 × 10-10) than for women with the CT (OR = 1.14, 95% CI = 1.09 to 1.19; P = 1.04 × 10-8) or TT (OR = 1.07, 95% CI = 1.01 to 1.14; P = .02) genotypes.

CONCLUSION

These results provide novel insights on a potential mechanism through which a SMAD7 variant, previously identified as a susceptibility locus for colorectal cancer, and BMI may influence colorectal cancer risk for women.

Identification of gene-environment interactions with marginal penalization

Gene-environment (G-E) interaction analysis has been extensively conducted for complex diseases. In marginal analysis, the common practice is to conduct likelihood-based (and other "standard") estimation with each marginal model, and then select significant G-E interactions and main effects based on p values and multiple comparisons adjustment. One limitation of this approach is that the identification results often do not respect the "main effects, interactions" hierarchy, which has been stressed in recent G-E interaction analyses. There is some recent effort tackling this problem, however, with very complex formulations. Another limitation of the common practice is that it may not perform well when regularization is needed, for example, because of "non-normal" distributions. In this article, we propose a marginal penalization approach which adopts a novel penalty to directly tackle the aforementioned problems. The proposed approach has a framework more coherent with that of the recently developed joint analysis methods and an intuitive formulation, and can be effectively realized. In simulation, it outperforms the popular significance-based analysis and simple penalization-based alternatives. Promising findings are made in the analysis of a single-nucleotide polymorphism and a gene expression data.

Subset-Based Analysis Using Gene-Environment Interactions for Discovery of Genetic Associations across Multiple Studies or Phenotypes

OBJECTIVES

Classical methods for combining summary data from genome-wide association studies only use marginal genetic effects, and power can be compromised in the presence of heterogeneity. We aim to enhance the discovery of novel associated loci in the presence of heterogeneity of genetic effects in subgroups defined by an environmental factor.

METHODS

We present a pvalue-assisted subset testing for associations (pASTA) framework that generalizes the previously proposed association analysis based on subsets (ASSET) method by incorporating gene-environment (G-E) interactions into the testing procedure. We conduct simulation studies and provide two data examples.

RESULTS

Simulation studies show that our proposal is more powerful than methods based on marginal associations in the presence of G-E interactions and maintains comparable power even in their absence. Both data examples demonstrate that our method can increase power to detect overall genetic associations and identify novel studies/phenotypes that contribute to the association.

CONCLUSIONS

Our proposed method can be a useful screening tool to identify candidate single nucleotide polymorphisms that are potentially associated with the trait(s) of interest for further validation. It also allows researchers to determine the most probable subset of traits that exhibit genetic associations in addition to the enhancement of power.

Using Bayes model averaging to leverage both gene main effects and G × E interactions to identify genomic regions in genome-wide association studies

Genome-wide association studies typically search for marginal associations between a single-nucleotide polymorphism (SNP) and a disease trait while gene-environment (G × E) interactions remain generally unexplored. More powerful methods beyond the simple case-control (CC) approach leverage either marginal effects or CC ascertainment to increase power. However, these potential gains depend on assumptions whose aptness is often unclear a priori. Here, we review G × E methods and use simulations to highlight performance as a function of main and interaction effects and the association of the two factors in the source population. Substantial variation in performance between methods leads to uncertainty as to which approach is most appropriate for any given analysis. We present a framework that (a) balances the robustness of a CC approach with the power of the case-only (CO) approach; (b) incorporates main SNP effects; (c) allows for incorporation of prior information; and (d) allows the data to determine the most appropriate model. Our framework is based on Bayes model averaging, which provides a principled statistical method for incorporating model uncertainty. We average over inclusion of parameters corresponding to the main and G × E interaction effects and the G-E association in controls. The resulting method exploits the joint evidence for main and interaction effects while gaining power from a CO equivalent analysis. Through simulations, we demonstrate that our approach detects SNPs within a wide range of scenarios with increased power over current methods. We illustrate the approach on a gene-environment scan in the USC Children's Health Study.

Detecting latent exposure in genome-wide association studies using a breakpoint model for logistic regression

Detecting gene-environment (G × E) interactions in the context of genome-wide association studies (GWAS) is a challenging problem since standard methods generally present a lack of power. An additional difficulty arises from the fact that the causal exposure is seldom observed and only a proxy of this exposure is observed. This leads to an additional drop in terms of power and it explains the failure of standard methods in detecting interactions, even very strong ones. In this article, we consider the latent exposure as a source of heterogeneity and we propose a new powerful method, named "Breakpoint Model for Logistic Regression" (BMLR), based on a breakpoint model, in order to detect G × E interactions when causal exposure is unobserved. First, the BMLR method is compared to the ordered-subset analysis for case-control method, which has been developed for the same purpose, through simulations. This highlights the ability of BMLR to detect the heterogeneity, and therefore, to detect interaction with latent exposure. Finally, the BMLR method is compared to standard methods, such as Plink, to perform a GWAS on a published realistic benchmark.

Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases

Recently, many new approaches, study designs, and statistical and analytical methods have emerged for studying gene-environment interactions (G×Es) in large-scale studies of human populations. There are opportunities in this field, particularly with respect to the incorporation of -omics and next-generation sequencing data and continual improvement in measures of environmental exposures implicated in complex disease outcomes. In a workshop called "Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases," held October 17-18, 2014, by the National Institute of Environmental Health Sciences and the National Cancer Institute in conjunction with the annual American Society of Human Genetics meeting, participants explored new approaches and tools that have been developed in recent years for G×E discovery. This paper highlights current and critical issues and themes in G×E research that need additional consideration, including the improved data analytical methods, environmental exposure assessment, and incorporation of functional data and annotations.

Update on the State of the Science for Analytical Methods for Gene-Environment Interactions

The analysis of gene-environment interaction (G×E) may hold the key for further understanding the etiology of many complex traits. The current availability of high-volume genetic data, the wide range in types of environmental data that can be measured, and the formation of consortiums of multiple studies provide new opportunities to identify G×E but also new analytical challenges. In this article, we summarize several statistical approaches that can be used to test for G×E in a genome-wide association study. These include traditional models of G×E in a case-control or quantitative trait study as well as alternative approaches that can provide substantially greater power. The latest methods for analyzing G×E with gene sets and with data in a consortium setting are summarized, as are issues that arise due to the complexity of environmental data. We provide some speculation on why detecting G×E in a genome-wide association study has thus far been difficult. We conclude with a description of software programs that can be used to implement most of the methods described in the paper.

The importance of gene-environment interactions in human obesity

The worldwide obesity epidemic has been mainly attributed to lifestyle changes. However, who becomes obese in an obesity-prone environment is largely determined by genetic factors. In the last 20 years, important progress has been made in the elucidation of the genetic architecture of obesity. In parallel with successful gene identifications, the number of gene-environment interaction (GEI) studies has grown rapidly. This paper reviews the growing body of evidence supporting gene-environment interactions in the field of obesity. Heritability, monogenic and polygenic obesity studies provide converging evidence that obesity-predisposing genes interact with a variety of environmental, lifestyle and treatment exposures. However, some skepticism remains regarding the validity of these studies based on several issues, which include statistical modelling, confounding, low replication rate, underpowered analyses, biological assumptions and measurement precision. What follows in this review includes (1) an introduction to the study of GEI, (2) the evidence of GEI in the field of obesity, (3) an outline of the biological mechanisms that may explain these interaction effects, (4) methodological challenges associated with GEI studies and potential solutions, and (5) future directions of GEI research. Thus far, this growing body of evidence has provided a deeper understanding of GEI influencing obesity and may have tremendous applications in the emerging field of personalized medicine and individualized lifestyle recommendations.

Exposure enriched outcome dependent designs for longitudinal studies of gene-environment interaction

Joint effects of genetic and environmental factors have been increasingly recognized in the development of many complex human diseases. Despite the popularity of case-control and case-only designs, longitudinal cohort studies that can capture time-varying outcome and exposure information have long been recommended for gene-environment (G × E) interactions. To date, literature on sampling designs for longitudinal studies of G × E interaction is quite limited. We therefore consider designs that can prioritize a subsample of the existing cohort for retrospective genotyping on the basis of currently available outcome, exposure, and covariate data. In this work, we propose stratified sampling based on summaries of individual exposures and outcome trajectories and develop a full conditional likelihood approach for estimation that adjusts for the biased sample. We compare the performance of our proposed design and analysis with combinations of different sampling designs and estimation approaches via simulation. We observe that the full conditional likelihood provides improved estimates for the G × E interaction and joint exposure effects over uncorrected complete-case analysis, and the exposure enriched outcome trajectory dependent design outperforms other designs in terms of estimation efficiency and power for detection of the G × E interaction. We also illustrate our design and analysis using data from the Normative Aging Study, an ongoing longitudinal cohort study initiated by the Veterans Administration in 1963. Copyright © 2017 John Wiley & Sons, Ltd.

Genes Interacting with Occupational Exposures to Low Molecular Weight Agents and Irritants on Adult-Onset Asthma in Three European Studies

BACKGROUND

The biological mechanisms by which cleaning products and disinfectants-an emerging risk factor-affect respiratory health remain incompletely evaluated. Studying genes by environment interactions (G × E) may help identify new genes related to adult-onset asthma.

OBJECTIVES

We identified interactions between genetic polymorphisms of a large set of genes involved in the response to oxidative stress and occupational exposures to low molecular weight (LMW) agents or irritants on adult-onset asthma.

METHODS

Our data came from three large European cohorts: Epidemiological Family-based Study of the Genetics and Environment of Asthma (EGEA), Swiss Cohort Study on Air Pollution and Lung and Heart Disease in Adults (SAPALDIA), and European Community Respiratory Health Survey in Adults (ECRHS). A candidate pathway-based strategy identified 163 genes involved in the response to oxidative stress and potentially related to exposures to LMW agents/irritants. Occupational exposures were evaluated using an asthma job-exposure matrix and job-specific questionnaires for cleaners and healthcare workers. Logistic regression models were used to detect G × E interactions, adjusted for age, sex, and population ancestry, in 2,599 adults (mean age, 47 years; 60% women, 36% exposed, 18% asthmatics). p-Values were corrected for multiple comparisons.

RESULTS

Ever exposure to LMW agents/irritants was associated with current adult-onset asthma [OR = 1.28 (95% CI: 1.04, 1.58)]. Eight single nucleotide polymorphism (SNP) by exposure interactions at five loci were found at p < 0.005: PLA2G4A (rs932476, chromosome 1), near PLA2R1 (rs2667026, chromosome 2), near RELA (rs931127, rs7949980, chromosome 11), PRKD1 (rs1958980, rs11847351, rs1958987, chromosome 14), and PRKCA (rs6504453, chromosome 17). Results were consistent across the three studies and after accounting for smoking.

CONCLUSIONS

Using a pathway-based selection process, we identified novel genes potentially involved in adult asthma by interaction with occupational exposure. These genes play a role in the NF-κB pathway, which is involved in inflammation. Citation: Rava M, Ahmed I, Kogevinas M, Le Moual N, Bouzigon E, Curjuric I, Dizier MH, Dumas O, Gonzalez JR, Imboden M, Mehta AJ, Tubert-Bitter P, Zock JP, Jarvis D, Probst-Hensch NM, Demenais F, Nadif R. 2017. Genes interacting with occupational exposures to low molecular weight agents and irritants on adult-onset asthma in three European studies. Environ Health Perspect 125:207-214; http://dx.doi.org/10.1289/EHP376.

Tests for Gene-Environment Interactions and Joint Effects With Exposure Misclassification

The number of methods for genome-wide testing of gene-environment (G-E) interactions continues to increase, with the aim of discovering new genetic risk factors and obtaining insight into the disease-gene-environment relationship. The relative performance of these methods, assessed on the basis of family-wise type I error rate and power, depends on underlying disease-gene-environment associations, estimates of which may be biased in the presence of exposure misclassification. This simulation study expands on a previously published simulation study of methods for detecting G-E interactions by evaluating the impact of exposure misclassification. We consider 7 single-step and modular screening methods for identifying G-E interaction at a genome-wide level and 7 joint tests for genetic association and G-E interaction, for which the goal is to discover new genetic susceptibility loci by leveraging G-E interaction when present. In terms of statistical power, modular methods that screen on the basis of the marginal disease-gene relationship are more robust to exposure misclassification. Joint tests that include main/marginal effects of a gene display a similar robustness, which confirms results from earlier studies. Our results offer an increased understanding of the strengths and limitations of methods for genome-wide searches for G-E interaction and joint tests in the presence of exposure misclassification.

Augmented case-only designs for randomized clinical trials with failure time endpoints

Under suitable assumptions and by exploiting the independence between inherited genetic susceptibility and treatment assignment, the case-only design yields efficient estimates for subgroup treatment effects and gene-treatment interaction in a Cox model. However it cannot provide estimates of the genetic main effect and baseline hazards, that are necessary to compute the absolute disease risk. For two-arm, placebo-controlled trials with rare failure time endpoints, we consider augmenting the case-only design with random samples of controls from both arms, as in the classical case-cohort sampling scheme, or with a random sample of controls from the active treatment arm only. The latter design is motivated by vaccine trials for cost-effective use of resources and specimens so that host genetics and vaccine-induced immune responses can be studied simultaneously in a bigger set of participants. We show that these designs can identify all parameters in a Cox model and that the efficient case-only estimator can be incorporated in a two-step plug-in procedure. Results in simulations and a data example suggest that incorporating case-only estimators in the classical case-cohort design improves the precision of all estimated parameters; sampling controls only in the active treatment arm attains a similar level of efficiency.

Mixed Linear Model Approaches of Association Mapping for Complex Traits Based on Omics Variants

Precise prediction for genetic architecture of complex traits is impeded by the limited understanding on genetic effects of complex traits, especially on gene-by-gene (GxG) and gene-by-environment (GxE) interaction. In the past decades, an explosion of high throughput technologies enables omics studies at multiple levels (such as genomics, transcriptomics, proteomics, and metabolomics). The analyses of large omics data, especially two-loci interaction analysis, are very time intensive. Integrating the diverse omics data and environmental effects in the analyses also remain challenges. We proposed mixed linear model approaches using GPU (Graphic Processing Unit) computation to simultaneously dissect various genetic effects. Analyses can be performed for estimating genetic main effects, GxG epistasis effects, and GxE environment interaction effects on large-scale omics data for complex traits, and for estimating heritability of specific genetic effects. Both mouse data analyses and Monte Carlo simulations demonstrated that genetic effects and environment interaction effects could be unbiasedly estimated with high statistical power by using the proposed approaches.

Estimating genetic effect sizes under joint disease-endophenotype models in presence of gene-environment interactions

Effects of genetic variants on the risk of complex diseases estimated from association studies are typically small. Nonetheless, variants may have important effects in presence of specific levels of environmental exposures, and when a trait related to the disease (endophenotype) is either normal or impaired. We propose polytomous and transition models to represent the relationship between disease, endophenotype, genotype and environmental exposure in family studies. Model coefficients were estimated using generalized estimating equations and were used to derive gene-environment interaction effects and genotype effects at specific levels of exposure. In a simulation study, estimates of the effect of a genetic variant were substantially higher when both an endophenotype and an environmental exposure modifying the variant effect were taken into account, particularly under transition models, compared to the alternative of ignoring the endophenotype. Illustration of the proposed modeling with the metabolic syndrome, abdominal obesity, physical activity and polymorphisms in the NOX3 gene in the Quebec Family Study revealed that the positive association of the A allele of rs1375713 with the metabolic syndrome at high levels of physical activity was only detectable in subjects without abdominal obesity, illustrating the importance of taking into account the abdominal obesity endophenotype in this analysis.

The impact of exposure-biased sampling designs on detection of gene-environment interactions in case-control studies with potential exposure misclassification

With limited funding and biological specimen availability, choosing an optimal sampling design to maximize power for detecting gene-by-environment (G-E) interactions is critical. Exposure-enriched sampling is often used to select subjects with rare exposures for genotyping to enhance power for tests of G-E effects. However, exposure misclassification (MC) combined with biased sampling can affect characteristics of tests for G-E interaction and joint tests for marginal association and G-E interaction. Here, we characterize the impact of exposure-biased sampling under conditions of perfect exposure information and exposure MC on properties of several methods for conducting inference. We assess the Type I error, power, bias, and mean squared error properties of case-only, case-control, and empirical Bayes methods for testing/estimating G-E interaction and a joint test for marginal G (or E) effect and G-E interaction across three biased sampling schemes. Properties are evaluated via empirical simulation studies. With perfect exposure information, exposure-enriched sampling schemes enhance power as compared to random selection of subjects irrespective of exposure prevalence but yield bias in estimation of the G-E interaction and marginal E parameters. Exposure MC modifies the relative performance of sampling designs when compared to the case of perfect exposure information. Those conducting G-E interaction studies should be aware of exposure MC properties and the prevalence of exposure when choosing an ideal sampling scheme and method for characterizing G-E interactions and joint effects.

Assessing gene-environment interactions for common and rare variants with binary traits using gene-trait similarity regression

Accounting for gene-environment (G×E) interactions in complex trait association studies can facilitate our understanding of genetic heterogeneity under different environmental exposures, improve the ability to discover susceptible genes that exhibit little marginal effect, provide insight into the biological mechanisms of complex diseases, help to identify high-risk subgroups in the population, and uncover hidden heritability. However, significant G×E interactions can be difficult to find. The sample sizes required for sufficient power to detect association are much larger than those needed for genetic main effects, and interactions are sensitive to misspecification of the main-effects model. These issues are exacerbated when working with binary phenotypes and rare variants, which bear less information on association. In this work, we present a similarity-based regression method for evaluating G×E interactions for rare variants with binary traits. The proposed model aggregates the genetic and G×E information across markers, using genetic similarity, thus increasing the ability to detect G×E signals. The model has a random effects interpretation, which leads to robustness against main-effect misspecifications when evaluating G×E interactions. We construct score tests to examine G×E interactions and a computationally efficient EM algorithm to estimate the nuisance variance components. Using simulations and data applications, we show that the proposed method is a flexible and powerful tool to study the G×E effect in common or rare variant studies with binary traits.

Attributing effects to interactions

A framework is presented that allows an investigator to estimate the portion of the effect of one exposure that is attributable to an interaction with a second exposure. We show that when the 2 exposures are statistically independent in distribution, the total effect of one exposure can be decomposed into a conditional effect of that exposure when the second is absent and also a component due to interaction. The decomposition applies on difference or ratio scales. We discuss how the components can be estimated using standard regression models, and how these components can be used to evaluate the proportion of the total effect of the primary exposure attributable to the interaction with the second exposure. In the setting in which one of the exposures affects the other, so that the 2 are no longer statistically independent in distribution, alternative decompositions are discussed. The various decompositions are illustrated with an example in genetic epidemiology. If it is not possible to intervene on the primary exposure of interest, the methods described in this article can help investigators to identify other variables that, if intervened upon, would eliminate the largest proportion of the effect of the primary exposure.

Atrazine and pregnancy outcomes: a systematic review of epidemiologic evidence

Atrazine (ATR) is a commonly used agricultural herbicide that has been the subject of epidemiologic studies assessing its relation to reproductive health problems. This review evaluates both the consistency and the quality of epidemiologic evidence testing the hypothesis that ATR exposure, at usually encountered levels, is a risk factor for birth defects, small for gestational age birth weight, prematurity, miscarriages, and problems of fetal growth and development. We followed the current methodological guidelines for systematic reviews by using two independent researchers to identify, retrieve, and evaluate the relevant epidemiologic literature on the relation of ATR to various adverse outcomes of birth and pregnancy. Each eligible paper was summarized with respect to its methods and results with particular attention to study design and exposure assessment, which have been cited as the main areas of weakness in ATR research. As a quantitative meta-analysis was not feasible, the study results were categorized qualitatively as positive, null, or mixed. The literature on ATR and pregnancy-related health outcomes is growing rapidly, but the quality of the data is poor with most papers using aggregate rather than individual-level information. Without good quality data, the results are difficult to assess; however, it is worth noting that none of the outcome categories demonstrated consistent positive associations across studies. Considering the poor quality of the data and the lack of robust findings across studies, conclusions about a causal link between ATR and adverse pregnancy outcomes are not warranted.

Integrative modeling of multiple genomic data from different types of genetic association studies

Genome-wide association studies (GWASs) and expression-/methylation-quantitative trait loci (eQTL/mQTL) studies constitute popular approaches for investigating the association of single nucleotide polymorphisms (SNPs) with disease and expression/methylation, respectively. Here, we propose to integrate QTL studies to more powerfully test the SNP effect on disease in GWASs when they are conducted among different subjects. We propose a model for the joint effect of SNPs, methylation, and gene expression on disease risk and obtain the marginal model for SNPs by integrating out methylation and expression. We characterize all possible causal relations among SNPs, methylation, and expression and study the corresponding null hypotheses of no SNP effect in terms of the regression coefficients in the joint model. We develop a score test for variance components of regression coefficients to evaluate the genetic effect. We further propose an omnibus test to accommodate different models. We illustrate the utility of the proposed method in an asthma GWAS study, a brain tumor study, and numerical simulations.

Identification of new genetic susceptibility loci for breast cancer through consideration of gene-environment interactions

Genes that alter disease risk only in combination with certain environmental exposures may not be detected in genetic association analysis. By using methods accounting for gene-environment (G × E) interaction, we aimed to identify novel genetic loci associated with breast cancer risk. Up to 34,475 cases and 34,786 controls of European ancestry from up to 23 studies in the Breast Cancer Association Consortium were included. Overall, 71,527 single nucleotide polymorphisms (SNPs), enriched for association with breast cancer, were tested for interaction with 10 environmental risk factors using three recently proposed hybrid methods and a joint test of association and interaction. Analyses were adjusted for age, study, population stratification, and confounding factors as applicable. Three SNPs in two independent loci showed statistically significant association: SNPs rs10483028 and rs2242714 in perfect linkage disequilibrium on chromosome 21 and rs12197388 in ARID1B on chromosome 6. While rs12197388 was identified using the joint test with parity and with age at menarche (P-values = 3 × 10(-07)), the variants on chromosome 21 q22.12, which showed interaction with adult body mass index (BMI) in 8,891 postmenopausal women, were identified by all methods applied. SNP rs10483028 was associated with breast cancer in women with a BMI below 25 kg/m(2) (OR = 1.26, 95% CI 1.15-1.38) but not in women with a BMI of 30 kg/m(2) or higher (OR = 0.89, 95% CI 0.72-1.11, P for interaction = 3.2 × 10(-05)). Our findings confirm comparable power of the recent methods for detecting G × E interaction and the utility of using G × E interaction analyses to identify new susceptibility loci.

A Tutorial on Interaction

In this tutorial, we provide a broad introduction to the topic of interaction between the effects of exposures. We discuss interaction on both additive and multiplicative scales using risks, and we discuss their relation to statistical models (e.g. linear, log-linear, and logistic models). We discuss and evaluate arguments that have been made for using additive or multiplicative scales to assess interaction. We further discuss approaches to presenting interaction analyses, different mechanistic forms of interaction, when interaction is robust to unmeasured confounding, interaction for continuous outcomes, qualitative or “crossover” interactions, methods for attributing effects to interactions, case-only estimators of interaction, and power and sample size calculations for additive and multiplicative interaction.

Gene-environment interactions in cancer epidemiology: a National Cancer Institute Think Tank report

Cancer risk is determined by a complex interplay of genetic and environmental factors. Genome-wide association studies (GWAS) have identified hundreds of common (minor allele frequency [MAF] > 0.05) and less common (0.01 < MAF < 0.05) genetic variants associated with cancer. The marginal effects of most of these variants have been small (odds ratios: 1.1-1.4). There remain unanswered questions on how best to incorporate the joint effects of genes and environment, including gene-environment (G × E) interactions, into epidemiologic studies of cancer. To help address these questions, and to better inform research priorities and allocation of resources, the National Cancer Institute sponsored a "Gene-Environment Think Tank" on January 10-11, 2012. The objective of the Think Tank was to facilitate discussions on (1) the state of the science, (2) the goals of G × E interaction studies in cancer epidemiology, and (3) opportunities for developing novel study designs and analysis tools. This report summarizes the Think Tank discussion, with a focus on contemporary approaches to the analysis of G × E interactions. Selecting the appropriate methods requires first identifying the relevant scientific question and rationale, with an important distinction made between analyses aiming to characterize the joint effects of putative or established genetic and environmental factors and analyses aiming to discover novel risk factors or novel interaction effects. Other discussion items include measurement error, statistical power, significance, and replication. Additional designs, exposure assessments, and analytical approaches need to be considered as we move from the current small number of success stories to a fuller understanding of the interplay of genetic and environmental factors.

Gene-environment interactions in genome-wide association studies: current approaches and new directions

BACKGROUND

Complex psychiatric traits have long been thought to be the result of a combination of genetic and environmental factors, and gene-environment interactions are thought to play a crucial role in behavioral phenotypes and the susceptibility and progression of psychiatric disorders. Candidate gene studies to investigate hypothesized gene-environment interactions are now fairly common in human genetic research, and with the shift toward genome-wide association studies, genome-wide gene-environment interaction studies are beginning to emerge.

METHODS

We summarize the basic ideas behind gene-environment interaction, and provide an overview of possible study designs and traditional analysis methods in the context of genome-wide analysis. We then discuss novel approaches beyond the traditional strategy of analyzing the interaction between the environmental factor and each polymorphism individually.

RESULTS

Two-step filtering approaches that reduce the number of polymorphisms tested for interactions can substantially increase the power of genome-wide gene-environment studies. New analytical methods including data-mining approaches, and gene-level and pathway-level analyses, also have the capacity to improve our understanding of how complex genetic and environmental factors interact to influence psychologic and psychiatric traits. Such methods, however, have not yet been utilized much in behavioral and mental health research.

CONCLUSIONS

Although methods to investigate gene-environment interactions are available, there is a need for further development and extension of these methods to identify gene-environment interactions in the context of genome-wide association studies. These novel approaches need to be applied in studies of psychology and psychiatry.

Environmental confounding in gene-environment interaction studies

We show that, in the presence of uncontrolled environmental confounding, joint tests for the presence of a main genetic effect and gene-environment interaction will be biased if the genetic and environmental factors are correlated, even if there is no effect of either the genetic factor or the environmental factor on the disease. When environmental confounding is ignored, such tests will in fact reject the joint null of no genetic effect with a probability that tends to 1 as the sample size increases. This problem with the joint test vanishes under gene-environment independence, but it still persists if estimating the gene-environment interaction parameter itself is of interest. Uncontrolled environmental confounding will bias estimates of gene-environment interaction parameters even under gene-environment independence, but it will not do so if the unmeasured confounding variable itself does not interact with the genetic factor. Under gene-environment independence, if the interaction parameter without controlling for the environmental confounder is nonzero, then there is gene-environment interaction either between the genetic factor and the environmental factor of interest or between the genetic factor and the unmeasured environmental confounder. We evaluate several recently proposed joint tests in a simulation study and discuss the implications of these results for the conduct of gene-environment interaction studies.

Case-only method for cause-specific hazards models with application to assessing differential vaccine efficacy by viral and host genetics

Cause-specific proportional hazards models are commonly used for analyzing competing risks data in clinical studies. Motivated by the objective to assess differential vaccine protection against distinct pathogen types in randomized preventive vaccine efficacy trials, we present an alternative case-only method to standard maximum partial likelihood estimation that applies to a rare failure event, e.g. acquisition of HIV infection. A logistic regression model is fit to the counts of cause-specific events (infecting pathogen type) within study arms, with an offset adjusting for the randomization ratio. This formulation of cause-specific hazard ratio estimation permits immediate incorporation of host-genetic factors to be assessed as effect modifiers, an important area of vaccine research for identifying immune correlates of protection, thus inheriting the estimation efficiency, and cost benefits of the case-only estimator commonly used for assessing gene-treatment interactions. The method is used to reassess HIV genotype-specific vaccine efficacy in the RV144 trial, providing nearly identical results to standard Cox methods, and to assess if and how this vaccine efficacy depends on Fc-γ receptor genes.

Gene-environment interactions in asthma and allergic diseases: challenges and perspectives

The concept of gene-environment (GxE) interactions has dramatically evolved in the last century and has now become a central theme in studies that assess the causes of human disease. Despite the numerous efforts to discover genes associated in asthma and allergy through various approaches, including the recent genome-wide association studies, investigation of GxE interactions has been mainly limited to candidate genes, candidate environmental exposures, or both. This review discusses the various strategies from hypothesis-driven strategies to the full agnostic search of GxE interactions with an illustration from recently published articles. Challenges raised by each piece of the puzzle (ie, phenotype, environment, gene, and analysis of GxE interaction) are put forward, and tentative solutions are proposed. New perspectives to integrate various types of data generated by new sequencing technologies and to progress toward a systems biology approach of disease are outlined. The future of a molecular network-based approach of disease to which GxE interactions are related requires space for innovative and multidisciplinary research. Assembling the various parts of a puzzle in a complex system could well occur in a way that might not necessarily follow the rules of logic.

Two-stage testing procedures with independent filtering for genome-wide gene-environment interaction

Several two-stage multiple testing procedures have been proposed to detect gene-environment interaction in genome-wide association studies. In this article, we elucidate general conditions that are required for validity and power of these procedures, and we propose extensions of two-stage procedures using the case-only estimator of gene-treatment interaction in randomized clinical trials. We develop a unified estimating equation approach to proving asymptotic independence between a filtering statistic and an interaction test statistic in a range of situations, including marginal association and interaction in a generalized linear model with a canonical link. We assess the performance of various two-stage procedures in simulations and in genetic studies from Women's Health Initiative clinical trials.

Exploring the interaction between SNP genotype and postmenopausal hormone therapy effects on stroke risk

BACKGROUND

Genome-wide association studies have identified several genomic regions that are associated with stroke risk, but these provide an explanation for only a small fraction of familial stroke aggregation. Genotype by environment interactions may contribute further to such an explanation. The Women's Health Initiative (WHI) clinical trial found increased stroke risk with postmenopausal hormone therapy (HT) and provides an efficient setting for evaluating genotype-HT interaction on stroke risk.

METHODS

We examined HT by genotype interactions for 392 SNPs selected from candidate gene studies, and 2,371 SNPs associated with changes in blood protein concentrations after hormone therapy, in analyses that included 2,045 postmenopausal women who developed stroke during WHI clinical trial and observational study follow-up and one-to-one matched controls. A two-stage procedure was implemented where SNPs passing the first stage screening based on marginal association with stroke risk were tested in the second stage for interaction with HT using case-only analysis.

RESULTS

The two-stage procedure identified two SNPs, rs2154299 and rs12194855, in the coagulation factor XIII subunit A (F13A1) region and two SNPs, rs630431 and rs560892, in the proprotein convertase subtilisin kexin 9 (PCSK9) region, with an estimated false discovery rate <0.05 based on interaction tests. Further analyses showed significant stroke risk interaction between these F13A1 SNPs and estrogen plus progestin (E+P) treatment for ischemic stroke and for ischemic and hemorrhagic stroke combined, and suggested interactions between PCSK9 SNPs with either E+P or estrogen-alone treatment.

CONCLUSIONS

Genotype by environment interaction information may help to define genomic regions relevant to stroke risk. Two-stage analysis among postmenopausal women generates novel hypotheses concerning the F13A1 and PCSK9 genomic regions and the effects of hormonal exposures on postmenopausal stroke risk for subsequent independent validation.