Detecting Gene-Environment Interactions for a Quantitative Trait in a Genome-Wide Association Study

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Sex-dependent placental methylation quantitative trait loci provide insight into the prenatal origins of childhood onset traits and conditions

Molecular quantitative trait loci (QTLs) allow us to understand the biology captured in genome-wide association studies (GWASs). The placenta regulates fetal development and shows sex differences in DNA methylation. We therefore hypothesized that placental methylation QTL (mQTL) explain variation in genetic risk for childhood onset traits, and that effects differ by sex. We analyzed 411 term placentas from two studies and found 49,252 methylation (CpG) sites with mQTL and 2,489 CpG sites with sex-dependent mQTL. All mQTL were enriched in regions that typically affect gene expression in prenatal tissues. All mQTL were also enriched in GWAS results for growth- and immune-related traits, but male- and female-specific mQTL were more enriched than cross-sex mQTL. mQTL colocalized with trait loci at 777 CpG sites, with 216 (28%) specific to males or females. Overall, mQTL specific to male and female placenta capture otherwise overlooked variation in childhood traits.

Whole Person Modeling: a transdisciplinary approach to mental health research

The growing global burden of mental illness has prompted calls for innovative research strategies. Theoretical models of mental health include complex contributions of biological, psychosocial, experiential, and other environmental influences. Accordingly, neuropsychiatric research has self-organized into largely isolated disciplines working to decode each individual contribution. However, research directly modeling objective biological measurements in combination with cognitive, psychological, demographic, or other environmental measurements is only now beginning to proliferate. This review aims to (1) to describe the landscape of modern mental health research and current movement towards integrative study, (2) to provide a concrete framework for quantitative integrative research, which we call Whole Person Modeling, (3) to explore existing and emerging techniques and methods used in Whole Person Modeling, and (4) to discuss our observations about the scarcity, potential value, and untested aspects of highly transdisciplinary research in general. Whole Person Modeling studies have the potential to provide a better understanding of multilevel phenomena, deliver more accurate diagnostic and prognostic tests to aid in clinical decision making, and test long standing theoretical models of mental illness. Some current barriers to progress include challenges with interdisciplinary communication and collaboration, systemic cultural barriers to transdisciplinary career paths, technical challenges in model specification, bias, and data harmonization, and gaps in transdisciplinary educational programs. We hope to ease anxiety in the field surrounding the often mysterious and intimidating world of transdisciplinary, data-driven mental health research and provide a useful orientation for students or highly specialized researchers who are new to this area.

A note on p-value multiple-testing adjustment for two-step genome-wide gene-environment interactions scans

Two-step testing is the state-of-the art approach for performing genome-wide interaction scans (GWIS). It is computationally efficient and yields higher power than standard single-step-based GWIS for virtually all biologically plausible scenarios. However, while two-step tests control the genome-wide type I error rate at the desired level, the lack of associated valid p-values can make it difficult for users to compare with single step-results. We show how multiple-testing adjusted p-values can be defined for two-step test based on standard multiple-testing theory, and how they can be in turn scaled to make valid comparisons with single-step tests possible.

Molecular mechanisms of environmental exposures and human disease

A substantial proportion of disease risk for common complex disorders is attributable to environmental exposures and pollutants. An appreciation of how environmental pollutants act on our cells to produce deleterious health effects has led to advances in our understanding of the molecular mechanisms underlying the pathogenesis of chronic diseases, including cancer and cardiovascular, neurodegenerative and respiratory diseases. Here, we discuss emerging research on the interplay of environmental pollutants with the human genome and epigenome. We review evidence showing the environmental impact on gene expression through epigenetic modifications, including DNA methylation, histone modification and non-coding RNAs. We also highlight recent studies that evaluate recently discovered molecular processes through which the environment can exert its effects, including extracellular vesicles, the epitranscriptome and the mitochondrial genome. Finally, we discuss current challenges when studying the exposome - the cumulative measure of environmental influences over the lifespan - and its integration into future environmental health research.

Improved two-step testing of genome-wide gene-environment interactions

Two-step tests for gene-environment (G × E $G\times E$ ) interactions exploit marginal single-nucleotide polymorphism (SNP) effects to improve the power of a genome-wide interaction scan. They combine a screening step based on marginal effects used to "bin" SNPs for weighted hypothesis testing in the second step to deliver greater power over single-step tests while preserving the genome-wide Type I error. However, the presence of many SNPs with detectable marginal effects on the trait of interest can reduce power by "displacing" true interactions with weaker marginal effects and by adding to the number of tests that need to be corrected for multiple testing. We introduce a new significance-based allocation into bins for Step-2G × E $G\times E$ testing that overcomes the displacement issue and propose a computationally efficient approach to account for multiple testing within bins. Simulation results demonstrate that these simple improvements can provide substantially greater power than current methods under several scenarios. An application to a multistudy collaboration for understanding colorectal cancer reveals a G × Sex interaction located near the SMAD7 gene.

Education interacts with genetic variants near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C to confer susceptibility to myopia

Myopia most often develops during school age, with the highest incidence in countries with intensive education systems. Interactions between genetic variants and educational exposure are hypothesized to confer susceptibility to myopia, but few such interactions have been identified. Here, we aimed to identify genetic variants that interact with education level to confer susceptibility to myopia. Two groups of unrelated participants of European ancestry from UK Biobank were studied. A 'Stage-I' sample of 88,334 participants whose refractive error (avMSE) was measured by autorefraction and a 'Stage-II' sample of 252,838 participants who self-reported their age-of-onset of spectacle wear (AOSW) but who did not undergo autorefraction. Genetic variants were prioritized via a 2-step screening process in the Stage-I sample: Step 1 was a genome-wide association study for avMSE; Step 2 was a variance heterogeneity analysis for avMSE. Genotype-by-education interaction tests were performed in the Stage-II sample, with University education coded as a binary exposure. On average, participants were 58 years-old and left full-time education when they were 18 years-old; 35% reported University level education. The 2-step screening strategy in the Stage-I sample prioritized 25 genetic variants (GWAS P < 1e-04; variance heterogeneity P < 5e-05). In the Stage-II sample, 19 of the 25 (76%) genetic variants demonstrated evidence of variance heterogeneity, suggesting the majority were true positives. Five genetic variants located near GJD2, RBFOX1, LAMA2, KCNQ5 and LRRC4C had evidence of a genotype-by-education interaction in the Stage-II sample (P < 0.002) and consistent evidence of a genotype-by-education interaction in the Stage-I sample. For all 5 variants, University-level education was associated with an increased effect of the risk allele. In this cohort, additional years of education were associated with an enhanced effect of genetic variants that have roles including axon guidance and the development of neuronal synapses and neural circuits.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

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.

Variance-quantitative trait loci enable systematic discovery of gene-environment interactions for cardiometabolic serum biomarkers

Gene-environment interactions represent the modification of genetic effects by environmental exposures and are critical for understanding disease and informing personalized medicine. These often induce differential phenotypic variance across genotypes; these variance-quantitative trait loci can be prioritized in a two-stage interaction detection strategy to greatly reduce the computational and statistical burden and enable testing of a broader range of exposures. We perform genome-wide variance-quantitative trait locus analysis for 20 serum cardiometabolic biomarkers by multi-ancestry meta-analysis of 350,016 unrelated participants in the UK Biobank, identifying 182 independent locus-biomarker pairs (p < 4.5×10-9). Most are concentrated in a small subset (4%) of loci with genome-wide significant main effects, and 44% replicate (p < 0.05) in the Women's Genome Health Study (N = 23,294). Next, we test each locus-biomarker pair for interaction across 2380 exposures, identifying 847 significant interactions (p < 2.4×10-7), of which 132 are independent (p < 0.05) after accounting for correlation between exposures. Specific examples demonstrate interaction of triglyceride-associated variants with distinct body mass- versus body fat-related exposures as well as genotype-specific associations between alcohol consumption and liver stress at the ADH1B gene. Our catalog of variance-quantitative trait loci and gene-environment interactions is publicly available in an online portal.

Gene-lifestyle interactions in the genomics of human complex traits

The role and biological significance of gene-environment interactions in human traits and diseases remain poorly understood. To address these questions, the CHARGE Gene-Lifestyle Interactions Working Group conducted series of genome-wide interaction studies (GWIS) involving up to 610,475 individuals across four ancestries for three lipids and four blood pressure traits, while accounting for interaction effects with drinking and smoking exposures. Here we used GWIS summary statistics from these studies to decipher potential differences in genetic associations and G×E interactions across phenotype-exposure-ancestry combinations, and to derive insights on the potential mechanistic underlying G×E through in-silico functional analyses. Our analyses show first that interaction effects likely contribute to the commonly reported ancestry-specific genetic effect in complex traits, and second, that some phenotype-exposures pairs are more likely to benefit from a greater detection power when accounting for interactions. It also highlighted modest correlation between marginal and interaction effects, providing material for future methodological development and biological discussions. We also estimated contributions to phenotypic variance, including in particular the genetic heritability conditional on the exposure, and heritability partitioned across a range of functional annotations and cell types. In these analyses, we found multiple instances of potential heterogeneity of functional partitions between exposed and unexposed individuals, providing new evidence for likely exposure-specific genetic pathways. Finally, along this work, we identified potential biases in methods used to jointly meta-analyze genetic and interaction effects. We performed simulations to characterize these limitations and to provide the community with guidelines for future G×E studies.

Genome-wide interaction study with major depression identifies novel variants associated with cognitive function

Major Depressive Disorder (MDD) often is associated with significant cognitive dysfunction. We conducted a meta-analysis of genome-wide interaction of MDD and cognitive function using data from four large European cohorts in a total of 3510 MDD cases and 6057 controls. In addition, we conducted analyses using polygenic risk scores (PRS) based on data from the Psychiatric Genomics Consortium (PGC) on the traits of MDD, Bipolar disorder (BD), Schizophrenia (SCZ), and mood instability (MIN). Functional exploration contained gene expression analyses and Ingenuity Pathway Analysis (IPA®). We identified a set of significantly interacting single nucleotide polymorphisms (SNPs) between MDD and the genome-wide association study (GWAS) of cognitive domains of executive function, processing speed, and global cognition. Several of these SNPs are located in genes expressed in brain, with important roles such as neuronal development (REST), oligodendrocyte maturation (TNFRSF21), and myelination (ARFGEF1). IPA® identified a set of core genes from our dataset that mapped to a wide range of canonical pathways and biological functions (MPO, FOXO1, PDE3A, TSLP, NLRP9, ADAMTS5, ROBO1, REST). Furthermore, IPA® identified upstream regulator molecules and causal networks impacting on the expression of dataset genes, providing a genetic basis for further clinical exploration (vitamin D receptor, beta-estradiol, tadalafil). PRS of MIN and meta-PRS of MDD, MIN and SCZ were significantly associated with all cognitive domains. Our results suggest several genes involved in physiological processes for the development and maintenance of cognition in MDD, as well as potential novel therapeutic agents that could be explored in patients with MDD associated cognitive dysfunction.

Two-step hypothesis testing to detect gene-environment interactions in a genome-wide scan with a survival endpoint

Defined by their genetic profile, individuals may exhibit differential clinical outcomes due to an environmental exposure. Identifying subgroups based on specific exposure-modifying genes can lead to targeted interventions and focused studies. Genome-wide interaction scans (GWIS) can be performed to identify such genes, but these scans typically suffer from low power due to the large multiple testing burden. We provide a novel framework for powerful two-step hypothesis tests for GWIS with a time-to-event endpoint under the Cox proportional hazards model. In the Cox regression setting, we develop an approach that prioritizes genes for Step-2 G × E testing based on a carefully constructed Step-1 screening procedure. Simulation results demonstrate this two-step approach can lead to substantially higher power for identifying gene-environment ( G × E ) interactions compared to the standard GWIS while preserving the family wise error rate over a range of scenarios. In a taxane-anthracycline chemotherapy study for breast cancer patients, the two-step approach identifies several gene expression by treatment interactions that would not be detected using the standard GWIS.

Uncovering Evidence for Endocrine-Disrupting Chemicals That Elicit Differential Susceptibility through Gene-Environment Interactions

Exposure to endocrine-disrupting chemicals (EDCs) is linked to myriad disorders, characterized by the disruption of the complex endocrine signaling pathways that govern development, physiology, and even behavior across the entire body. The mechanisms of endocrine disruption involve a complex system of pathways that communicate across the body to stimulate specific receptors that bind DNA and regulate the expression of a suite of genes. These mechanisms, including gene regulation, DNA binding, and protein binding, can be tied to differences in individual susceptibility across a genetically diverse population. In this review, we posit that EDCs causing such differential responses may be identified by looking for a signal of population variability after exposure. We begin by summarizing how the biology of EDCs has implications for genetically diverse populations. We then describe how gene-environment interactions (GxE) across the complex pathways of endocrine signaling could lead to differences in susceptibility. We survey examples in the literature of individual susceptibility differences to EDCs, pointing to a need for research in this area, especially regarding the exceedingly complex thyroid pathway. Following a discussion of experimental designs to better identify and study GxE across EDCs, we present a case study of a high-throughput screening signal of putative GxE within known endocrine disruptors. We conclude with a call for further, deeper analysis of the EDCs, particularly the thyroid disruptors, to identify if these chemicals participate in GxE leading to differences in susceptibility.

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.

Robust semiparametric gene-environment interaction analysis using sparse boosting

For the pathogenesis of complex diseases, gene-environment (G-E) interactions have been shown to have important implications. G-E interaction analysis can be challenging with the need to jointly analyze a large number of main effects and interactions and to respect the "main effects, interactions" hierarchical constraint. Extensive methodological developments on G-E interaction analysis have been conducted in recent literature. Despite considerable successes, most of the existing studies are still limited as they cannot accommodate long-tailed distributions/data contamination, make the restricted assumption of linear effects, and cannot effectively accommodate missingness in E variables. To directly tackle these problems, a semiparametric model is assumed to accommodate nonlinear effects, and the Huber loss function and Qn estimator are adopted to accommodate long-tailed distributions/data contamination. A regression-based multiple imputation approach is developed to accommodate missingness in E variables. For model estimation and selection of relevant variables, we adopt an effective sparse boosting approach. The proposed approach is practically well motivated, has intuitive formulations, and can be effectively realized. In extensive simulations, it significantly outperforms multiple direct competitors. The analysis of The Cancer Genome Atlas data on stomach adenocarcinoma and cutaneous melanoma shows that the proposed approach makes sensible discoveries with satisfactory prediction and stability.

Quantile regression analysis reveals widespread evidence for gene-environment or gene-gene interactions in myopia development

A genetic contribution to refractive error has been confirmed by the discovery of more than 150 associated variants in genome-wide association studies (GWAS). Environmental factors such as education and time outdoors also demonstrate strong associations. Currently however, the extent of gene-environment or gene-gene interactions in myopia is unknown. We tested the hypothesis that refractive error-associated variants exhibit effect size heterogeneity, a hallmark feature of genetic interactions. Of 146 variants tested, evidence of non-uniform, non-linear effects were observed for 66 (45%) at Bonferroni-corrected significance (P < 1.1 × 10-4) and 128 (88%) at nominal significance (P < 0.05). LAMA2 variant rs12193446, for example, had an effect size varying from -0.20 diopters (95% CI -0.18 to -0.23) to -0.89 diopters (95% CI -0.71 to -1.07) in different individuals. SNP effects were strongest at the phenotype extremes and weaker in emmetropes. A parsimonious explanation for these findings is that gene-environment or gene-gene interactions in myopia are pervasive.

A Genome-Wide Association Study for Susceptibility to Visual Experience-Induced Myopia

Purpose

The rapid rise in prevalence over recent decades and high heritability of myopia suggest a role for gene-environment (G × E) interactions in myopia susceptibility. Few such G × E interactions have been discovered to date. We aimed to test the hypothesis that genetic analysis of susceptibility to visual experience-induced myopia in an animal model would identify novel G × E interaction loci.

Methods

Chicks aged 7 days (n = 987) were monocularly deprived of form vision for 4 days. A genome-wide association study (GWAS) was carried out in the 20% of chicks most susceptible and least susceptible to form deprivation (n = 380). There were 304,963 genetic markers tested for association with the degree of induced axial elongation in treated versus control eyes (A-scan ultrasonography). A GWAS candidate region was examined in the following three human cohorts: CREAM consortium (n = 44,192), UK Biobank (n = 95,505), and Avon Longitudinal Study of Parents and Children (ALSPAC; n = 4989).

Results

A locus encompassing the genes PIK3CG and PRKAR2B was genome-wide significantly associated with myopia susceptibility in chicks (lead variant rs317386235, P = 9.54e-08). In CREAM and UK Biobank GWAS datasets, PIK3CG and PRKAR2B were enriched for strongly-associated markers (meta-analysis lead variant rs117909394, P = 1.7e-07). In ALSPAC participants, rs117909394 had an age-dependent association with refractive error (-0.22 diopters [D] change over 8 years, P = 5.2e-04) and nearby variant rs17153745 showed evidence of a G × E interaction with time spent reading (effect size -0.23 D, P = 0.022).

Conclusions

This work identified the PIK3CG-PRKAR2B locus as a mediator of susceptibility to visually induced myopia in chicks and suggests a role for this locus in conferring susceptibility to myopia in human cohorts.

Genome-Wide Association Study (GWAS) on Bilirubin Concentrations in Subjects with Metabolic Syndrome: Sex-Specific GWAS Analysis and Gene-Diet Interactions in a Mediterranean Population

Although, for decades, increased serum bilirubin concentrations were considered a threatening sign of underlying liver disease and had been associated with neonatal jaundice, data from recent years show that bilirubin is a powerful antioxidant and suggest that slightly increased serum bilirubin concentrations are protective against oxidative stress-related diseases, such as cardiovascular diseases. Therefore, a better understanding of the gene-diet interactions in determining serum bilirubin concentrations is needed. None of the previous genome-wide association studies (GWAS) on bilirubin concentrations has been stratified by sex. Therefore, considering the increasing interest in incorporating the gender perspective into nutritional genomics, our main aim was to carry out a GWAS on total serum bilirubin concentrations in a Mediterranean population with metabolic syndrome, stratified by sex. Our secondary aim was to explore, as a pilot study, the presence of gene-diet interactions at the GWAS level. We included 430 participants (188 men and 242 women, aged 55–75 years, and with metabolic syndrome) in the PREDIMED Plus-Valencia study. Global and sex-specific GWAS were undertaken to analyze associations and gene-diet interaction on total serum bilirubin. Adherence (low and high) to the Mediterranean diet (MedDiet) was analyzed as the dietary modulator. In the GWAS, we detected more than 55 SNPs associated with serum bilirubin at p < 5 × 10⁻⁸ (GWAS level). The top-ranked were four SNPs (rs4148325 (p = 9.25 × 10⁻²⁴), rs4148324 (p = 9.48 × 10⁻²⁴), rs6742078 (p = 1.29 × 10⁻²³), rs887829 (p = 1.39 × 10⁻²³), and the rs4148324 (p = 9.48 × 10⁻²⁴)) in the UGT1A1 (UDP glucuronosyltransferase family 1 member A1) gene, which replicated previous findings revealing the UGT1A1 as the major locus. In the sex-specific GWAS, the top-ranked SNPs at the GWAS level were similar in men and women (the lead SNP was the rs4148324-UGT1A1 in both men (p = 4.77 × 10⁻¹¹) and women (p = 2.15 × 10⁻¹⁴), which shows homogeneous genetic results for the major locus. There was more sex-specific heterogeneity for other minor genes associated at the suggestive level of GWAS significance (p < 1 × 10⁻⁵). We did not detect any gene-MedDiet interaction at p < 1 × 10⁻⁵ for the major genetic locus, but we detected some gene-MedDiet interactions with other genes at p < 1 × 10⁻⁵, and even at the GWAS level for the IL17B gene (p = 3.14 × 10⁻⁸). These interaction results, however, should be interpreted with caution due to our small sample size. In conclusion, our study provides new data, with a gender perspective, on genes associated with total serum bilirubin concentrations in men and women, and suggests possible additional modulations by adherence to MedDiet.

Robust gene-environment interaction analysis using penalized trimmed regression

In biomedical and epidemiological studies, gene-environment (G-E) interactions have been shown to importantly contribute to the etiology and progression of many complex diseases. Most existing approaches for identifying G-E interactions are limited by the lack of robustness against outliers/contaminations in response and predictor spaces. In this study, we develop a novel robust G-E identification approach using the trimmed regression technique under joint modeling. A robust data-driven criterion and stability selection are adopted to determine the trimmed subset which is free from both vertical outliers and leverage points. An effective penalization approach is developed to identify important G-E interactions, respecting the "main effects, interactions" hierarchical structure. Extensive simulations demonstrate the better performance of the proposed approach compared to multiple alternatives. Interesting findings with superior prediction accuracy and stability are observed in the analysis of TCGA data on cutaneous melanoma and breast invasive carcinoma.

Genome-wide interaction with the insulin secretion locus MTNR1B reveals CMIP as a novel type 2 diabetes susceptibility gene in African Americans

Although type 2 diabetes (T2D) results from metabolic defects in insulin secretion and insulin sensitivity, most of the genetic risk loci identified to date relates to insulin secretion. We reported that T2D loci influencing insulin sensitivity may be identified through interactions with insulin secretion loci, thereby leading to T2D. Here, we hypothesize that joint testing of variant main effects and interaction effects with an insulin secretion locus increases power to identify genetic interactions leading to T2D. We tested this hypothesis with an intronic MTNR1B SNP, rs10830963, which is associated with acute insulin response to glucose, a dynamic measure of insulin secretion. rs10830963 was tested for interaction and joint (main + interaction) effects with genome-wide data in African Americans (2,452 cases and 3,772 controls) from five cohorts. Genome-wide genotype data (Affymetrix Human Genome 6.0 array) was imputed to a 1000 Genomes Project reference panel. T2D risk was modeled using logistic regression with rs10830963 dosage, age, sex, and principal component as predictors. Joint effects were captured using the Kraft two degrees of freedom test. Genome-wide significant (P < 5 × 10-8 ) interaction with MTNR1B and joint effects were detected for CMIP intronic SNP rs17197883 (Pinteraction  = 1.43 × 10-8 ; Pjoint  = 4.70 × 10-8 ). CMIP variants have been nominally associated with T2D, fasting glucose, and adiponectin in individuals of East Asian ancestry, with high-density lipoprotein, and with waist-to-hip ratio adjusted for body mass index in Europeans. These data support the hypothesis that additional genetic factors contributing to T2D risk, including insulin sensitivity loci, can be identified through interactions with insulin secretion loci.

Robust identification of gene-environment interactions for prognosis using a quantile partial correlation approach

Gene-environment (G-E) interactions have important implications for the etiology and progression of many complex diseases. Compared to continuous markers and categorical disease status, prognosis has been less investigated, with the additional challenges brought by the unique characteristics of survival outcomes. Most of the existing G-E interaction approaches for prognosis data share the limitation that they cannot accommodate long-tailed or contaminated outcomes. In this study, for prognosis data, we develop a robust G-E interaction identification approach using the censored quantile partial correlation (CQPCorr) technique. The proposed approach is built on the quantile regression technique (and hence has a solid statistical basis), uses weights to easily accommodate censoring, and adopts partial correlation to identify important interactions while properly controlling for the main genetic and environmental effects. In simulation, it outperforms multiple competitors with more accurate identification. In the analysis of TCGA data on lung cancer and melanoma, biologically sensible findings different from using the alternatives are made.

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.

Impact of a Eukaryotic Translation Initiation Factor 3a Polymorphism on Susceptibility to Gastric Cancer

OBJECTIVE

To investigate single nucleotide polymorphisms in the eukaryotic translation initiation factor 3a (eIF3a) gene and the risk for gastric cancer within the Chinese population.

SUBJECTS AND METHODS

A total of 322 patients with gastric cancer were selected as the patient group and 340 non-gastric cancer patients were selected as the control group using the case-control method. Polymerase chain reaction-sequence-specific primer technology was leveraged to genotype the rs77382849 single nucleotide polymorphism in the eIF3a gene. The demographic characteristics of the study population and other exposures to risk factors were collected. Unconditional logistic regression analysis was performed to determine the association between the risk factors and gastric cancer.

RESULTS

A higher frequency of the eIF3a rs77382849 GG homozygote genotype was observed in the gastric cancer patients compared with the controls (63.98 vs. 54.41%, p < 0.05). After adjustment of exposure risks, such as age, gender, smoking, and drinking, the rs77382849 single nucleotide polymorphism was still associated with susceptibility to gastric cancer. When the eIF3a rs77382849 GG homozygote genotype was used as the reference group, the GA genotype (GA vs. GG: OR = 0.545, 95% CI: 0.386-0.769, p = 0.001) and AA genotype (AA vs. GG: OR = 0.245, 95% CI: 0.072-0.836, p = 0.025) were both correlated with a significantly decreased risk for gastric cancer development.

CONCLUSION

An association between eIF3a rs77382849 polymorphism and susceptibility to gastric cancer was observed in these Chinese patients.