Gene Environment Research Articles

Inclusion of Gene-Gene and Gene-Environment Interactions Unlikely to Dramatically Improve Risk Prediction for Complex Diseases

Genome-wide association studies have identified hundreds of common genetic variants associated with the risk of multifactorial diseases. However, their impact on discrimination and risk prediction is limited. It has been suggested that the identification of gene-gene (G-G) and gene-environment (G-E) interactions would improve disease prediction and facilitate prevention. We conducted a simulation study to explore the potential improvement in discrimination if G-G and G-E interactions exist and are known. We used three diseases (breast cancer, type 2 diabetes, and rheumatoid arthritis) as motivating examples. We show that the inclusion of G-G and G-E interaction effects in risk-prediction models is unlikely to dramatically improve the discrimination ability of these models.

Pseudo semiparametric maximum likelihood estimation exploiting gene environment independence for population-based case-control studies with complex samples

Advances in human genetics have led to epidemiological investigations not only of the effects of genes alone but also of gene-environment (G-E) interaction. A widely accepted design strategy in the study of how G-E relate to disease risks is the population-based case-control study (PBCCS). For simple random samples, semiparametric methods for testing G-E have been developed by Chatterjee and Carroll in 2005. The use of complex sampling in PBCCS that involve differential probabilities of sample selection of cases and controls and possibly cluster sampling is becoming more common. Two complexities, weighting for selection probabilities and intracluster correlation of observations, are induced by the complex sampling. We develop pseudo-semiparametric maximum likelihood estimators (pseudo-SPMLE) that apply to PBCCS with complex sampling. We study the finite sample performance of the pseudo-SPMLE using simulations and illustrate the pseudo-SPMLE with a US case-control study of kidney cancer.

Efficient designs of gene–environment interaction studies: implications of Hardy–Weinberg equilibrium and gene–environment independence

It is important to investigate whether genetic susceptibility variants exercise the same effects in populations that are differentially exposed to environmental risk factors. Here, we assess the power of four two-stage case-control design strategies for assessing multiplicative gene-environment (G-E) interactions or for assessing genetic or environmental effects in the presence of G-E interactions. We considered a di-allelic single nucleotide polymorphism G and a binary environmental variable E under the constraints of G-E independence and Hardy-Weinberg equilibrium and used the Wald statistic for all tests. We concluded that (i) for testing G-E interactions or genetic effects in the presence of G-E interactions when data for E are fully available, it is preferable to ascertain data for G in a subsample of cases with similar numbers of exposed and unexposed and a random subsample of controls; and (ii) for testing G-E interactions or environmental effects in the presence of G-E interactions when data for G are fully available, it is preferable to ascertain data for E in a subsample of cases that has similar numbers for each genotype and a random subsample of controls. In addition, supplementing external control data to an existing case-control sample leads to improved power for assessing effects of G or E in the presence of G-E interactions.

A latent variable approach to study gene-environment interactions in the presence of multiple correlated exposures.

Many existing cohort studies initially designed to investigate disease risk as a function of environmental exposures have collected genomic data in recent years with the objective of testing for gene-environment interaction (G × E) effects. In environmental epidemiology, interest in G × E arises primarily after a significant effect of the environmental exposure has been documented. Cohort studies often collect rich exposure data; as a result, assessing G × E effects in the presence of multiple exposure markers further increases the burden of multiple testing, an issue already present in both genetic and environment health studies. Latent variable (LV) models have been used in environmental epidemiology to reduce dimensionality of the exposure data, gain power by reducing multiplicity issues via condensing exposure data, and avoid collinearity problems due to presence of multiple correlated exposures. We extend the LV framework to characterize gene-environment interaction in presence of multiple correlated exposures and genotype categories. Further, similar to what has been done in case-control G × E studies, we use the assumption of gene-environment (G-E) independence to boost the power of tests for interaction. The consequences of making this assumption, or the issue of how to explicitly model G-E association has not been previously investigated in LV models. We postulate a hierarchy of assumptions about the LV model regarding the different forms of G-E dependence and show that making such assumptions may influence inferential results on the G, E, and G × E parameters. We implement a class of shrinkage estimators to data adaptively trade-off between the most restrictive to most flexible form of G-E dependence assumption and note that such class of compromise estimators can serve as a benchmark of model adequacy in LV models. We demonstrate the methods with an example from the Early Life Exposures in Mexico City to Neuro-Toxicants Study of lead exposure, iron metabolism genes, and birth weight.

Abstract 3742: Gene-environment Interactions in esophageal adenocarcinoma risk: A case-only analysis

Abstract The incidence of esophageal adenocarcinoma (EA) has increased approximately 500% in Western countries over the last four decades. The rapidly increasing incidence and sporadically development feature of EA suggest that gene-environment (G-E) interactions dominate the etiology. However, how G-E contributes to EA carcinogenesis is still poorly understood. We used a case-only approach to test the effect of G-E interactions in the occurrence of EA. 1330 functional and/or tagging SNPs selected from 354 cancer-associated genes were genotyped in 335 Caucasian EA cases. G-E (gastroesophageal reflux symptoms, GERD; BMI; and smoking) interactions were assessed by a two-step approach. First, we applied random forest (RF) to screen for important SNPs that had main effects as well as interaction effects. Second, we used case-only logistic regression (LR) models to estimate the G-E interaction effect, adjusting for covariates and false-discovery rate (FDR). RF analyses identified three sets of SNPs (94 SNPs-GERD, 38 SNPs-smoking, and 44 SNPs-BMI, respectively) that had the highest importance scores and lowest classification error rates. Further case-only LR analysis revealed that multiple interaction markers were significantly associated with EA risk: GERD*rs2237051 of EGF(OR = 1.39, P = 3.50E-07), GERD*rs2440 of XRCC5 (OR = 1,80, P = 0.0008), GERD*rs2237051*rs2440 (OR = 1.85, P = 0.004), smoking*rs10842514 of KRAS (OR = 0.85, P = 0.012), BMI*rs2305742 of IL12RB1 (OR = 0.84, P = 0.003), BMI*rs11568820 of VDR (OR = 0.97, P = 0.0005), BMI*rs11244142 of ABL1 (OR = 1.47, P = 2.09E-05), BMI*rs2800975 of XPA*rs743572 of CYP17A (OR = 3.62, P = 0.0009), and BMI*rs11244142*rs20547 of IL13 (OR = 1.89, P = 0.023). These results indicate that genetic variability may contribute to EA development through interactions with environmental factors. Our data also suggest that G-E interactions may be G- and E-specific in modulating the risk of EA. Supported by grants: NIH grants CA92824, CA74386, CA90578, and CA119650); Flight Attendant Medical Research Institute (FAMRI) grant 062459_YCSA; the Kevin Jackson Memorial Fund and Alan Brown Chair of Molecular Genomics. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3742. doi:10.1158/1538-7445.AM2011-3742

Structural Equation Modeling of Gene-Environment Interactions in Coronary Heart Disease

Coronary heart disease (CHD) is a complex disease, which is influenced not only by genetic and environmental factors but also by gene-environment (GE) interactions in interconnected biological pathways or networks. The classical methods are inadequate for identifying GE interactions due to the complex relationships among risk factors, mediating risk factors (e.g., hypertension, blood lipids, and glucose), and CHD. Our aim was to develop a two-level structural equation model (SEM) to identify genes and GE interactions in the progress of CHD to take into account the causal structure among mediating risk factors and CHD (Level 1), and hierarchical family structure (Level 2). The method was applied to the Framingham Heart Study (FHS) Offspring Cohort data. Our approach has several advantages over classical methods: (1) it provides important insight into how genes and contributing factors affect CHD by investigating the direct, indirect, and total effects; and (2) it aids the development of biological models that more realistically reflect the complex biological pathways or networks. Using our method, we are able to detect GE interaction of SERPINE1 and body mass index (BMI) on CHD, which has not been reported. We conclude that SEM modeling of GE interaction can be applied in the analysis of complex epidemiological data sets.

Gene-environment interaction in the etiology of mathematical ability using SNP sets.

Mathematics ability and disability is as heritable as other cognitive abilities and disabilities, however its genetic etiology has received relatively little attention. In our recent genome-wide association study of mathematical ability in 10-year-old children, 10 SNP associations were nominated from scans of pooled DNA and validated in an individually genotyped sample. In this paper, we use a ‘SNP set’ composite of these 10 SNPs to investigate gene-environment (GE) interaction, examining whether the association between the 10-SNP set and mathematical ability differs as a function of ten environmental measures in the home and school in a sample of 1888 children with complete data. We found two significant GE interactions for environmental measures in the home and the school both in the direction of the diathesis-stress type of GE interaction: The 10-SNP set was more strongly associated with mathematical ability in chaotic homes and when parents are negative.

Extended Bayesian Model Averaging in Generalized Linear Mixed Models Applied to Schizophrenia Family Data

The study of disease etiology and the search for susceptible genes of schizophrenia have attracted scientists' attention for decades. Many findings however are inconsistent, possibly due to the higher order interactions involving multi-dimensional genetic and environmental factors or due to the commingling of different ethnic groups. Several studies applied generalized linear mixed models (GLMMs) with family data to identify the genetic contribution to, and environmental influence on, schizophrenia, and to clarify the existence and sources of familial aggregation. Based on an extended Bayesian model averaging (EBMA) procedure, here we estimate the gene-gene (GG) and gene-environment (GE) interactions, and heritability of schizophrenia via variance components of random-effects in GLMMs. Our proposal takes into account the uncertainty in covariates and in genetic model structures, where each competing model includes environmental and genetic covariates, and GE and GG interactions. Simulation studies are conducted to compare the performance of the EBMA approach, permutation test procedure and GEE method. We also illustrate this approach with data from singleton and multiplex schizophrenia families. The results indicate that EBMA is a flexible and stable tool in exploring true candidate genes, and GE and GG interactions, after adjusting for explanatory variables and correlation structures within family members.

Gene–environment interplay in internalizing disorders: consistent findings across six environmental risk factors

Behavior genetic methods can help to elucidate gene-environment (G-E) interplay in the development of internalizing (INT) disorders (i.e., major depression and anxiety disorders). To date, however, no study has conducted a comprehensive analysis examining multiple environmental risk factors with the purpose of delineating general mechanisms of G-E influence in the development of INT disorders. The sample consisted of 1315 male and female twin pairs participating in the age 17 assessment of the Minnesota Twin Family Study. Quantitative G-E interplay models were used to examine how genetic and environmental risk for INT disorders changes as a function of environmental context. Multiple measures and informants were employed to construct composite measures of INT disorders and six environmental risk factors including: stressful life events, mother-child and father-child relationship problems, antisocial and prosocial peer affiliation, and academic achievement and engagement. Significant moderation effects were detected between each environmental risk factor and INT such that in the context of greater environmental adversity, nonshared environmental factors became more important in the etiology of INT symptoms. Our results are consistent with the interpretation that environmental stressors have a causative effect on the emergence of INT disorders. The consistency of our results suggests a general mechanism of environmental influence on INT disorders regardless of the specific form of environmental risk.

Parent–child pair design for detecting gene–environment interactions in complex diseases

It is becoming clear that the etiology of complex diseases involves not only genetic and environmental factors but also gene-environment (GE) interactions. Therefore, it is important to take account of all these factors to improve the power of an epidemiological study design. We propose here a novel parent-child pair (PCP) design for this purpose. In comparison with conventional designs, this approach has the following advantages: (a) PCP is a 4 x 16 design consisting of pairs of parent-child (PC) genotype statuses, PC exposure statuses and PC disease statuses. Therefore, it utilizes more information than the traditional approaches in association studies; (b) It can determine whether findings in studies of association between disease and genetic or environmental factors and their interaction are spurious, arising from Hardy-Weinberg disequilibrium or the other factors; (c) Since the information from both parents and children of the PC pairs are used in this design, it has high power for detecting association of candidate gene, exposure with a complex disease and GE interaction. We also present a set of estimates of relative risks of candidate genes, exposures and GE interactions under the multiplicative model and a method for computing the sample size requirements to test for these relative risks in the context of the PCP design.

Evaluation of two methods for assessing gene‐environment interactions using data from the Danish case‐control study of facial clefts

Epidemiological investigations have begun to consider gene-environment (GE) interactions as potential risk factors for many diseases, including several different birth defects. However, traditional methodological approaches for the analysis of case-control data tend to have low power for detection of interaction effects. A log-linear approach that can impose the assumption that the genotype and exposure of interest occur independently in the population has been proposed as a potentially more powerful method for assessing GE interactions but has not been widely applied in the published literature. The present analyses were undertaken to compare the results obtained when stratified analyses and a log-linear approach were used to assess potential GE interactions. The analyses were conducted using data from a population-based, case-control study conducted in Denmark and considered associations between nonsyndromic cleft lip with or without cleft palate (CL+/-P), infant genotype for variants of RAR-alpha, TGF-alpha, TGF-beta3, and MSX1, and maternal exposure to smoking, alcohol, and multivitamins. Neither the stratified nor the log-linear analyses provided evidence that that risk of CL+/-P is influenced by any of the GE interactions that were evaluated, despite the potential increase in power offered by the latter approach. Further, the analyses highlight concerns regarding the power to reject the assumption of independence of the genetic and environmental factor of interest in the controls and related concerns regarding the validity of results obtained using the log-linear approach when the underlying assumption is violated. The potential increase in power offered by the log-linear approach is offset by concerns regarding the validity of this approach when the independence assumption is violated.

Social Context in Gene–Environment Interactions: Retrospect and Prospect

While many behavioral scientists believe that gene-environment (GE) interactions play an important and perhaps pervasive role in human development and aging, little attention has been devoted to a fundamental conceptual issue: What is it about social context that could alter gene expression? We draw on existing examples of GE interactions to formulate a typology that identifies a set of generic mechanisms by which E moderates G. Empirical studies suggest four ideal types: Social context can trigger a genetic diathesis, compensate for a genetic diathesis, act as a control to prevent behaviors for which there is a genetic predisposition, and enhance adaptation through proximal processes. This typology highlights several problems, however, with prior empirical research, which may explain, in part, why so few GE interactions have actually been observed. These problems include inattention to the dynamic nature of social experience, the manifold, often-intercorrelated dimensions of social context ("EE interactions"), mediators that link social context and the genotype, and analytic models that examine GE interactions as processes that characterize individual development. In turn, these insights call for the integration of life course sociology and behavioral genetics to foster ways of studying genes, context, and aging.