GxEsum: a novel approach to estimate the phenotypic variance explained by genome-wide GxE interaction based on GWAS summary statistics for biobank-scale data

Genotype × environment interactions in gene regulation and complex traits

Genotype × environment interactions (GxE) have long been recognized as a key mechanism underlying human phenotypic variation. Technological developments over the past 15 years have dramatically expanded our appreciation of the role of GxE in both gene regulation and complex traits. The richness and complexity of these datasets also required parallel efforts to develop robust and sensitive statistical and computational approaches. Although our understanding of the genetic architecture of molecular and complex traits has been maturing, a large proportion of complex trait heritability remains unexplained. Furthermore, there are increasing efforts to characterize the effect of environmental exposure on human health. We therefore review GxE in human gene regulation and complex traits, advocating for a comprehensive approach that jointly considers genetic and environmental factors in human health and disease.

Gene-environment interactions in human health

Gene-environment interactions (G × E), the interplay of genetic variation with environmental factors, have a pivotal impact on human complex traits and diseases. Statistically, G × E can be assessed by determining the deviation from expectation of predictive models based solely on the phenotypic effects of genetics or environmental exposures. Despite the unprecedented, widespread and diverse use of G × E analytical frameworks, heterogeneity in their application and reporting hinders their applicability in public health. In this Review, we discuss study design considerations as well as G × E analytical frameworks to assess polygenic liability dependent on the environment, to identify specific genetic variants exhibiting G × E, and to characterize environmental context for these dynamics. We conclude with recommendations to address the most common challenges and pitfalls in the conceptualization, methodology and reporting of G × E studies, as well as future directions.

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.

Mitigating type 1 error inflation and power loss in GxE PRS: Genotype-environment interaction in polygenic risk score models

The use of polygenic risk score (PRS) models has transformed the field of genetics by enabling the prediction of complex traits and diseases based on an individual's genetic profile. However, the impact of genotype-environment interaction (GxE) on the performance and applicability of PRS models remains a crucial aspect to be explored. Currently, existing genotype-environment interaction polygenic risk score (GxE PRS) models are often inappropriately used, which can result in inflated type 1 error rates and compromised results. In this study, we propose novel GxE PRS models that jointly incorporate additive and interaction genetic effects although also including an additional quadratic term for nongenetic covariates, enhancing their robustness against model misspecification. Through extensive simulations, we demonstrate that our proposed models outperform existing models in terms of controlling type 1 error rates and enhancing statistical power. Furthermore, we apply the proposed models to real data, and report significant GxE effects. Specifically, we highlight the impact of our models on both quantitative and binary traits. For quantitative traits, we uncover the GxE modulation of genetic effects on body mass index by alcohol intake frequency. In the case of binary traits, we identify the GxE modulation of genetic effects on hypertension by waist-to-hip ratio. These findings underscore the importance of employing a robust model that effectively controls type 1 error rates, thus preventing the occurrence of spurious GxE signals. To facilitate the implementation of our approach, we have developed an innovative R software package called GxEprs, specifically designed to detect and estimate GxE effects. Overall, our study highlights the importance of accurate GxE modeling and its implications for genetic risk prediction, although providing a practical tool to support further research in this area.

Statistical methods for gene-environment interaction analysis

Most human complex phenotypes result from multiple genetic and environmental factors and their interactions. Understanding the mechanisms by which genetic and environmental factors interact offers valuable insights into the genetic architecture of complex traits and holds great potential for advancing precision medicine. The emergence of large population biobanks has led to the development of numerous statistical methods aiming at identifying gene-environment interactions (G × E). In this review, we present state-of-the-art statistical methodologies for G × E analysis. We will survey a spectrum of approaches for single-variant G × E mapping, followed by various techniques for polygenic G × E analysis. We conclude this review with a discussion on the future directions and challenges in G × E research.

Re-analysis and meta-analysis of summary statistics from gene-environment interaction studies

MOTIVATION

statistics from genome-wide association studies enable many valuable downstream analyses that are more efficient than individual-level data analysis while also reducing privacy concerns. As growing sample sizes enable better-powered analysis of gene-environment interactions, there is a need for gene-environment interaction-specific methods that manipulate and use summary statistics.

RESULTS

We introduce two tools to facilitate such analysis, with a focus on statistical models containing multiple gene-exposure and/or gene-covariate interaction terms. REGEM (RE-analysis of GEM summary statistics) uses summary statistics from a single, multi-exposure genome-wide interaction study to derive analogous sets of summary statistics with arbitrary sets of exposures and interaction covariate adjustments. METAGEM (META-analysis of GEM summary statistics) extends current fixed-effects meta-analysis models to incorporate multiple exposures from multiple studies. We demonstrate the value and efficiency of these tools by exploring alternative methods of accounting for ancestry-related population stratification in genome-wide interaction study in the UK Biobank as well as by conducting a multi-exposure genome-wide interaction study meta-analysis in cohorts from the diabetes-focused ProDiGY consortium. These programs help to maximize the value of summary statistics from diverse and complex gene-environment interaction studies.

AVAILABILITY AND IMPLEMENTATION

REGEM and METAGEM are open-source projects freely available at https://github.com/large-scale-gxe-methods/REGEM and https://github.com/large-scale-gxe-methods/METAGEM.

Response to anti-IL17 therapy in inflammatory disease is not strongly impacted by genetic background

Response to the anti-IL17 monoclonal antibody secukinumab is heterogeneous, and not all participants respond to treatment. Understanding whether this heterogeneity is driven by genetic variation is a key aim of pharmacogenetics and could influence precision medicine approaches in inflammatory diseases. Using changes in disease activity scores across 5,218 genotyped individuals from 19 clinical trials across four indications (psoriatic arthritis, psoriasis, ankylosing spondylitis, and rheumatoid arthritis), we tested whether genetics predicted response to secukinumab. We did not find any evidence of association between treatment response and common variants, imputed HLA alleles, polygenic risk scores of disease susceptibility, or cross-disease components of shared genetic risk. This suggests that anti-IL17 therapy is equally effective regardless of an individual's genetic background, a finding that has important implications for future genetic studies of biological therapy response in inflammatory diseases.

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.

Genotype-by-environment interactions in chronic back pain

BACKGROUND CONTEXT

Chronic back pain (CBP) is a common debilitating condition with substantial societal impact. While understanding genotype-by-environment (GxE) interactions may be crucial to achieving the goals of personalized medicine, there are few large-scale studies investigating this topic for CBP. None of them systematically explore multiple CBP risk factors.

PURPOSE

To estimate the extent to which genetic effects on CBP are modified by known demographic and clinical risk factors.

RESEARCH DESIGN

Case-control study, genome-wide GxE interaction study.

PATIENT SAMPLE

Data on up to 331,610 unrelated participants (57,881 CBP cases and 273,729 controls) from the UK Biobank cohort were used. UK Biobank is a prospective cohort with collected deep genetic and phenotypic data on approximately 500,000 individuals across the UK.

OUTCOME MEASURES

Self-reported chronic back pain.

METHODS

We applied a whole-genome approach to estimate the proportion of phenotypic variance explained by interactions between genotype and 12 known risk factors. We also analyzed if effects of common single-nucleotide polymorphisms on CBP are changed in presence of known risk factors.

RESULTS

The results indicate a modest, if any, modification of genetic effects by examined risk factors in CBP. Our estimates suggest that detecting such weak effects would require a sample size of millions of individuals.

CONCLUSIONS

The GxE interactions with examined common risk factors for CBP are either weak or absent. Interactions of such magnitude are unlikely to have the potential to inform and influence treatment strategies. Risk estimation models may use common genetic variation and the considered risk factors as independent predictors, without accounting for GxE.

Heuristic Approach Uncovering Biological Significance of Gene-Lifestyle Interactions in Cardiometabolic Traits

BACKGROUND

Gene-lifestyle interaction studies using genome-wide association studies (GWAS) data contribute to a better understanding of individual responses to environmental exposures.

OBJECTIVES

Herein, we aimed at assessing the biological significance of overlapping genes reported in gene-lifestyle interaction studies in cardiometabolic health.

METHOD

A heuristic analysis of genes reporting significant interactions related to cardiometabolic traits was performed to determine the biological pathways common to the different traits.

RESULTS

A total of 873 genes were analyzed. Fine and condensed phenotypic solutions were obtained from overlapping genes common to more than one trait.

CONCLUSIONS

This study revealed significant metabolic pathways associated with the impact of gene-environment interactions on cardiometabolic risk. Graphical Abstract: Publicly available data in cloud-based repositories were used to perform enrichment analyses of genes previously described in GWAS studies that showed interaction with lifestyles. From the enriched pathways, cluster analysis was performed to group enriched metabolic disorders.

Quality control and analytic best practices for testing genetic models of sex differences in large populations

Phenotypic sex-based differences exist for many complex traits. In other cases, phenotypes may be similar, but underlying biology may vary. Thus, sex-aware genetic analyses are becoming increasingly important for understanding the mechanisms driving these differences. To this end, we provide a guide outlining the current best practices for testing various models of sex-dependent genetic effects in complex traits and disease conditions, noting that this is an evolving field. Insights from sex-aware analyses will not only teach us about the biology of complex traits but also aid in achieving the goals of precision medicine and health equity for all.

Gene-environment interaction explains a part of missing heritability in human body mass index

Gene-environment (G×E) interaction could partially explain missing heritability in traits; however, the magnitudes of G×E interaction effects remain unclear. Here, we estimate the heritability of G×E interaction for body mass index (BMI) by subjecting genome-wide interaction study data of 331,282 participants in the UK Biobank to linkage disequilibrium score regression (LDSC) and linkage disequilibrium adjusted kinships-software for estimating SNP heritability from summary statistics (LDAK-SumHer) analyses. Among 14 obesity-related lifestyle factors, MET score, pack years of smoking, and alcohol intake frequency significantly interact with genetic factors in both analyses, accounting for the partial variance of BMI. The G×E interaction heritability (%) and standard error of these factors by LDSC and LDAK-SumHer are as follows: MET score, 0.45% (0.12) and 0.65% (0.24); pack years of smoking, 0.52% (0.13) and 0.93% (0.26); and alcohol intake frequency, 0.32% (0.10) and 0.80% (0.17), respectively. Moreover, these three factors are partially validated for their interactions with genetic factors in other obesity-related traits, including waist circumference, hip circumference, waist-to-hip ratio adjusted with BMI, and body fat percentage. Our results suggest that G×E interaction may partly explain the missing heritability in BMI, and two G×E interaction loci identified could help in understanding the genetic architecture of obesity.

The Polygenic Risk Score Knowledge Base offers a centralized online repository for calculating and contextualizing polygenic risk scores

The process of identifying suitable genome-wide association (GWA) studies and formatting the data to calculate multiple polygenic risk scores on a single genome can be laborious. Here, we present a centralized polygenic risk score calculator currently containing over 250,000 genetic variant associations from the NHGRI-EBI GWAS Catalog for users to easily calculate sample-specific polygenic risk scores with comparable results to other available tools. Polygenic risk scores are calculated either online through the Polygenic Risk Score Knowledge Base (PRSKB; https://prs.byu.edu ) or via a command-line interface. We report study-specific polygenic risk scores across the UK Biobank, 1000 Genomes, and the Alzheimer's Disease Neuroimaging Initiative (ADNI), contextualize computed scores, and identify potentially confounding genetic risk factors in ADNI. We introduce a streamlined analysis tool and web interface to calculate and contextualize polygenic risk scores across various studies, which we anticipate will facilitate a wider adaptation of polygenic risk scores in future disease research.

Considering hormone-sensitive cancers as a single disease in the UK biobank reveals shared aetiology

Hormone-related cancers, including cancers of the breast, prostate, ovaries, uterine, and thyroid, globally contribute to the majority of cancer incidence. We hypothesize that hormone-sensitive cancers share common genetic risk factors that have rarely been investigated by previous genomic studies of site-specific cancers. Here, we show that considering hormone-sensitive cancers as a single disease in the UK Biobank reveals shared genetic aetiology. We observe that a significant proportion of variance in disease liability is explained by the genome-wide single nucleotide polymorphisms (SNPs), i.e., SNP-based heritability on the liability scale is estimated as 10.06% (SE 0.70%). Moreover, we find 55 genome-wide significant SNPs for the disease, using a genome-wide association study. Pair-wise analysis also estimates positive genetic correlations between some pairs of hormone-sensitive cancers although they are not statistically significant. Our finding suggests that heritable genetic factors may be a key driver in the mechanism of carcinogenesis shared by hormone-sensitive cancers.

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.

An integrative analysis of genomic and exposomic data for complex traits and phenotypic prediction

Complementary to the genome, the concept of exposome has been proposed to capture the totality of human environmental exposures. While there has been some recent progress on the construction of the exposome, few tools exist that can integrate the genome and exposome for complex trait analyses. Here we propose a linear mixed model approach to bridge this gap, which jointly models the random effects of the two omics layers on phenotypes of complex traits. We illustrate our approach using traits from the UK Biobank (e.g., BMI and height for N ~ 35,000) with a small fraction of the exposome that comprises 28 lifestyle factors. The joint model of the genome and exposome explains substantially more phenotypic variance and significantly improves phenotypic prediction accuracy, compared to the model based on the genome alone. The additional phenotypic variance captured by the exposome includes its additive effects as well as non-additive effects such as genome-exposome (gxe) and exposome-exposome (exe) interactions. For example, 19% of variation in BMI is explained by additive effects of the genome, while additional 7.2% by additive effects of the exposome, 1.9% by exe interactions and 4.5% by gxe interactions. Correspondingly, the prediction accuracy for BMI, computed using Pearson's correlation between the observed and predicted phenotypes, improves from 0.15 (based on the genome alone) to 0.35 (based on the genome and exposome). We also show, using established theories, that integrating genomic and exposomic data can be an effective way of attaining a clinically meaningful level of prediction accuracy for disease traits. In conclusion, the genomic and exposomic effects can contribute to phenotypic variation via their latent relationships, i.e. genome-exposome correlation, and gxe and exe interactions, and modelling these effects has a potential to improve phenotypic prediction accuracy and thus holds a great promise for future clinical practice.