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

Systematic discovery of gene-environment interactions underlying the human plasma proteome in UK Biobank

Understanding how gene-environment interactions (GEIs) influence the circulating proteome could aid in biomarker discovery and validation. The presence of GEIs can be inferred from single nucleotide polymorphisms that associate with phenotypic variability - termed variance quantitative trait loci (vQTLs). Here, vQTL association studies are performed on plasma levels of 1463 proteins in 52,363 UK Biobank participants. A set of 677 independent vQTLs are identified across 568 proteins. They include 67 variants that lack conventional additive main effects on protein levels. Over 1100 GEIs are identified between 101 proteins and 153 environmental exposures. GEI analyses uncover possible mechanisms that explain why 13/67 vQTL-only sites lack corresponding main effects. Additional analyses also highlight how age, sex, epistatic interactions and statistical artefacts may underscore associations between genetic variation and variance heterogeneity. This study establishes the most comprehensive database yet of vQTLs and GEIs for the human proteome.

Gene‒Environment Interaction Between CAST Gene and Eye-Rubbing in the Chinese Keratoconus Cohort Study: A Case-Only Study

Purpose

Keratoconus (KC), characterized by progressive corneal protrusion and thinning, is a complex disease influenced by the combination of genetic and environmental factors. The purpose of this study was to explore potential gene‒environment interaction between the calpastatin (CAST) gene and eye-rubbing in KC.

Methods

A case-only study including 930 patients (676 patients with eye-rubbing and 254 patients without eye-rubbing) from the Chinese Keratoconus (CKC) cohort study was performed in the present study. Genotyping of single nucleotide polymorphism (SNP) was conducted using the Illumina Infinium Human Asian Screening Array (ASA) Beadchip. The gene‒environment interactions between CAST gene and eye-rubbing were analyzed using PLINK version 1.90. The interactions between CAST genotypes and eye-rubbing were analyzed by logistic regression models. The SNP-SNP-environment interactions were analyzed using generalized multifactor dimensionality reduction (GMDR).

Results

Three SNPs in CAST gene, namely, rs26515, rs27991, and rs9314177, reached the significance threshold for interactions (defined as P < 2.272 × 10-3). Notably, the minor alleles of these three SNPs exhibited negative interactions with eye-rubbing in KC. The results of logistic regression models revealed that the minor allele homozygotes and heterozygotes of rs26515, rs27991, and rs9314177 also exhibited negative interactions with eye-rubbing. Furthermore, GMDR analysis revealed the significant SNP-SNP-environment interactions among rs26515, rs27991, rs9314177, and eye-rubbing in KC.

Conclusions

This study identified rs26515, rs27991, and rs9314177 in CAST gene existed gene-environment interactions with eye-rubbing in KC, which is highly important for understanding the underlying biological mechanisms of KC and guiding precision prevention and proper management.

Gene-environment interactions within a precision environmental health framework

Understanding the complex interplay of genetic and environmental factors in disease etiology and the role of gene-environment interactions (GEIs) across human development stages is important. We review the state of GEI research, including challenges in measuring environmental factors and advantages of GEI analysis in understanding disease mechanisms. We discuss the evolution of GEI studies from candidate gene-environment studies to genome-wide interaction studies (GWISs) and the role of multi-omics in mediating GEI effects. We review advancements in GEI analysis methods and the importance of large-scale datasets. We also address the translation of GEI findings into precision environmental health (PEH), showcasing real-world applications in healthcare and disease prevention. Additionally, we highlight societal considerations in GEI research, including environmental justice, the return of results to participants, and data privacy. Overall, we underscore the significance of GEI for disease prediction and prevention and advocate for integrating the exposome into PEH omics studies.

Association between biological ageing and periodontitis: Evidence from a cross-sectional survey and multi-omics Mendelian randomization analysis

AIM

To investigate the relationship and potential causality between biological ageing and periodontitis.

MATERIALS AND METHODS

We obtained the National Health and Nutrition Examination Survey (NHANES) and genome-wide association study (GWAS) summary statistics as well as single-cell sequencing data. Multivariate regression analysis based on cross-sectional data, Mendelian randomization (MR) and multi-omics integration analysis were employed to explore the causal association and potential molecular mechanisms between biological ageing and periodontitis. Additionally, two-step MR mediation analysis explored the risk factors in biological ageing-mediated periodontitis.

RESULTS

We analysed data from 3189 participants in the NHANES data and found that higher biological age was associated with increased risk of periodontitis. MR analyses revealed causal associations between biological age measures and periodontitis risk. Frailty (odds ratio [OR] = 2.08, 95% confidence interval [CI]: 1.04-4.18, p = .039) and GrimAge acceleration (OR = 1.16, 95% CI: 1.01-1.32, p = .033) were causally associated with periodontitis risk, and these results were validated in a large-scale meta-periodontitis GWAS dataset. Additionally, the risk effects of body mass index, waist circumference and lifetime smoking on periodontitis were partially mediated by frailty and GrimAge acceleration.

CONCLUSIONS

Evidence from cross-sectional survey and MR analysis suggests that biological ageing increases the risk of periodontitis. Additionally, improving the associated risk factors can help prevent both ageing and periodontitis.

Two genome-wide interaction loci modify the association of nonsteroidal anti-inflammatory drugs with colorectal cancer

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

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.

Data-Driven Characterization of Genetic Variability in Disease Pathways and Pesticide-Induced Nervous System Disease in the United States Population

BACKGROUND

Genetic susceptibility to chemicals is incompletely characterized. However, nervous system disease development following pesticide exposure can vary in a population, implying some individuals may have higher genetic susceptibility to pesticide-induced nervous system disease.

OBJECTIVES

We aimed to build a computational approach to characterize single-nucleotide polymorphisms (SNPs) implicated in chemically induced adverse outcomes and used this framework to assess the link between differential population susceptibility to pesticides and human nervous system disease.

METHODS

We integrated publicly available datasets of Chemical-Gene, Gene-Pathway, and SNP-Disease associations to build Chemical-Pathway-Gene-SNP-Disease linkages for humans. As a case study, we integrated these linkages with spatialized pesticide application data for the US from 1992 to 2018 and spatialized nervous system disease rates for 2018. Through this, we characterized SNPs that may be important in states with high disease occurrence based on the pesticides used there.

RESULTS

We found that the number of SNP hits per pesticide in US states positively correlated with disease incidence and prevalence for Alzheimer's disease, Parkinson disease, and multiple sclerosis. We performed frequent itemset mining to differentiate pesticides used over time in states with high and low disease occurrence and found that only 19% of pesticide sets overlapped between 10 states with high disease occurrence and 10 states with low disease occurrence rates, and more SNPs were implicated in pathways in high disease occurrence states. Through a cross-validation of subsets of five high and low disease occurrence states, we characterized SNPs, genes, pathways, and pesticides more frequently implicated in high disease occurrence states.

DISCUSSION

Our findings support that pesticides contribute to nervous system disease, and we developed priority lists of SNPs, pesticides, and pathways for further study. This data-driven approach can be adapted to other chemicals, diseases, and locations to characterize differential population susceptibility to chemical exposures. https://doi.org/10.1289/EHP14108.

Understanding and applying gene-environment interactions: a guide for nutrition professionals with an emphasis on integration in African research settings

Noncommunicable diseases (NCDs) are influenced by the interplay between genetics and environmental exposures, particularly diet. However, many healthcare professionals, including nutritionists and dietitians, have limited genetic background and, therefore, they may lack understanding of gene-environment interactions (GxEs) studies. Even researchers deeply involved in nutrition studies, but with a focus elsewhere, can struggle to interpret, evaluate, and conduct GxE studies. There is an urgent need to study African populations that bear a heavy burden of NCDs, demonstrate unique genetic variability, and have cultural practices resulting in distinctive environmental exposures compared with Europeans or Americans, who are studied more. Although diverse and rapidly changing environments, as well as the high genetic variability of Africans and difference in linkage disequilibrium (ie, certain gene variants are inherited together more often than expected by chance), provide unparalleled potential to investigate the omics fields, only a small percentage of studies come from Africa. Furthermore, research evidence lags behind the practices of companies offering genetic testing for personalized medicine and nutrition. We need to generate more evidence on GxEs that also considers continental African populations to be able to prevent unethical practices and enable tailored treatments. This review aims to introduce nutrition professionals to genetics terms and valid methods to investigate GxEs and their challenges, and proposes ways to improve quality and reproducibility. The review also provides insight into the potential contributions of nutrigenetics and nutrigenomics to the healthcare sphere, addresses direct-to-consumer genetic testing, and concludes by offering insights into the field's future, including advanced technologies like artificial intelligence and machine learning.

Overcoming the "feast or famine" effect: improved interaction testing in genome-wide association studies

In genetic association analysis of complex traits, detection of interaction (either GxG or GxE) can help to elucidate the genetic architecture and biological mechanisms underlying the trait. Detection of interaction in a genome-wide association study (GWAS) can be methodologically challenging for various reasons, including a high burden of multiple comparisons when testing for epistasis between all possible pairs of a set of genomewide variants, as well as heteroscedasticity effects occurring in the presence of GxG or GxE interaction. In this paper, we address the problem of an even more striking phenomenon that we call the "feast or famine" effect that occurs when testing interaction in a genomewide context. As we verify, even in a simplified setting in which there is no interaction at all (and so no heteroscedasticity), in a GWAS to detect GxG or GxE interaction with a fixed genetic variant or environmental factor, the distribution of the genome-wide p-values under the null hypothesis is not the i.i.d. uniform one that is commonly assumed. Using standard methods, even if all SNPs are independent, some GWASs will have systematically underinflated p-values ("feast"), and others will have systematically overinflated p-values ("famine"), which can lead to false detection of interaction, reduced power, inconsistent results across studies, and failure to replicate true signal. This startling phenomenon is specific to detection of interaction in a GWAS, and it may partly explain why such detection has so far proved challenging and difficult to replicate. We show theoretically that the key cause of this phenomenon is which variables are conditioned on in the analysis, and this suggests an approach to correct the problem by changing the way the conditioning is done. Using this insight, we have developed the TINGA method to adjust the interaction test statistics to make their p-values closer to uniform under the null hypothesis. In simulations we show that TINGA both controls type 1 error and improves power. TINGA allows for covariates and population structure through use of a linear mixed model and accounts for heteroscedasticity. We apply TINGA to detection of epistasis in a study of flowering time in Arabidopsis thaliana.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Stress, Environment and Early Psychosis

Existing literature provides extended evidence of the close relationship between stress dysregulation, environmental insults, and psychosis onset. Early stress can sensitize genetically vulnerable individuals to future stress, modifying their risk for developing psychotic phenomena. Neurobiological substrate of the aberrant stress response to hypothalamic-pituitary-adrenal axis dysregulation, disrupted inflammation processes, oxidative stress increase, gut dysbiosis, and altered brain signaling, provides mechanistic links between environmental risk factors and the development of psychotic symptoms. Early-life and later-life exposures may act directly, accumulatively, and repeatedly during critical neurodevelopmental time windows. Environmental hazards, such as pre- and perinatal complications, traumatic experiences, psychosocial stressors, and cannabis use might negatively intervene with brain developmental trajectories and disturb the balance of important stress systems, which act together with recent life events to push the individual over the threshold for the manifestation of psychosis. The current review presents the dynamic and complex relationship between stress, environment, and psychosis onset, attempting to provide an insight into potentially modifiable factors, enhancing resilience and possibly influencing individual psychosis liability.

Bayesian Approaches in Exploring Gene-environment and Gene-gene Interactions: A Comprehensive Review

Rapid advancements in high-throughput biological techniques have facilitated the generation of high-dimensional omics datasets, which have provided a solid foundation for precision medicine and prognosis prediction. Nonetheless, the problem of missing heritability persists. To solve this problem, it is essential to explain the genetic structure of disease incidence risk and prognosis by incorporating interactions. The development of the Bayesian theory has provided new approaches for developing models for interaction identification and estimation. Several Bayesian models have been developed to improve the accuracy of model and identify the main effect, gene-environment (G×E) and gene-gene (G×G) interactions. Studies based on single-nucleotide polymorphisms (SNPs) are significant for the exploration of rare and common variants. Models based on the effect heredity principle and group-based models are relatively flexible and do not require strict constraints when dealing with the hierarchical structure between the main effect and interactions (M-I). These models have a good interpretability of biological mechanisms. Machine learning-based Bayesian approaches are highly competitive in improving prediction accuracy. These models provide insights into the mechanisms underlying the occurrence and progression of complex diseases, identify more reliable biomarkers, and develop higher predictive accuracy. In this paper, we provide a comprehensive review of these Bayesian approaches.

Pathological imaging-assisted cancer gene-environment interaction analysis

Gene-environment (G-E) interactions have important implications for cancer outcomes and phenotypes beyond the main G and E effects. Compared to main-effect-only analysis, G-E interaction analysis more seriously suffers from a lack of information caused by higher dimensionality, weaker signals, and other factors. It is also uniquely challenged by the "main effects, interactions" variable selection hierarchy. Effort has been made to bring in additional information to assist cancer G-E interaction analysis. In this study, we take a strategy different from the existing literature and borrow information from pathological imaging data. Such data are a "byproduct" of biopsy, enjoys broad availability and low cost, and has been shown as informative for modeling prognosis and other cancer outcomes/phenotypes in recent studies. Building on penalization, we develop an assisted estimation and variable selection approach for G-E interaction analysis. The approach is intuitive, can be effectively realized, and has competitive performance in simulation. We further analyze The Cancer Genome Atlas (TCGA) data on lung adenocarcinoma (LUAD). The outcome of interest is overall survival, and for G variables, we analyze gene expressions. Assisted by pathological imaging data, our G-E interaction analysis leads to different findings with competitive prediction performance and stability.

Diet-Wide Association, Genetic Susceptibility and Colorectal Cancer Risk: A Prospective Cohort Study

BACKGROUND

Both genetic and dietary factors play significant roles in the etiology of colorectal cancer (CRC). To evaluate the relationship between certain food exposures and the risk of CRC, we carried out a large-scale association analysis in the UK Biobank.

METHODS

The associations of 139 foods and nutrients' intake with CRC risk were assessed among 118,210 participants. A polygenic risk score (PRS) of CRC was created to explore any interaction between dietary factors and genetic susceptibility in CRC risk. The hazard ratio (HR) and 95% confidence interval (CI) of CRC risk linked to dietary variables and PRS were estimated using Cox regression models. Multiple comparisons were corrected using the error discovery rate (FDR).

RESULTS

During a mean follow-up of 12.8 years, 1466 incidents of CRC were identified. In the UK Biobank, alcohol and white bread were associated with increased CRC risk, and their HRs were 1.08 (95% CI: 1.03-1.14; FDRP = 0.028) and 1.10 (95% CI: 1.05-1.16; FDRP = 0.003), whereas dietary fiber, calcium, magnesium, phosphorus, and manganese intakes were inversely associated. We found no evidence of any PRS-nutrient interaction relationship in relation to CRC risk.

CONCLUSIONS

Our results show that higher intakes of alcohol and white bread are associated with increased CRC risk, whilst dietary fiber, calcium, magnesium, phosphorus, and manganese are inversely associated.

Joint associations between established genetic susceptibility loci, pesticide exposures, and risk of prostate cancer

More than 200 genetic variants have been independently associated with prostate cancer risk. Studies among farmers have also observed increased prostate cancer risk associated with exposure to specific organophosphate (fonofos, terbufos, malathion, dimethoate) and organochlorine (aldrin, chlordane) insecticides. We examined the joint associations between these pesticides, established prostate cancer loci, and prostate cancer risk among 1,162 cases (588 aggressive) and 2,206 frequency-matched controls nested in the Agricultural Health Study cohort. History of lifetime pesticide use was combined with a polygenic risk score (PRS) generated using 256 established prostate cancer risk variants. Logistic regression models estimated the joint associations of the pesticides, the PRS, and the 256 individual genetic variants with risk of total and aggressive prostate cancer. Likelihood ratio tests assessed multiplicative interaction. We observed interaction between ever use of fonofos and the PRS in relation to total and aggressive prostate cancer risk. Compared to the reference group (never use, PRS < median), men with ever use of fonofos and PRS > median had elevated risks of total (OR 1.35 [1.06-1.73], p-interaction = 0.03) and aggressive (OR 1.49 [1.09-2.04], p-interaction = 0.19) prostate cancer. There was also suggestion of interaction between pesticides and individual genetic variants occurring in regions associated with DNA damage response (CDH3, EMSY genes) and with variants related to altered androgen receptor-driven transcriptional programs critical for prostate cancer. Our study provides evidence that men with greater genetic susceptibility to prostate cancer may be at higher risk if they are also exposed to pesticides and suggests potential mechanisms by which pesticides may increase prostate cancer risk.

Allostatic load in early adolescence: gene / environment contributions and relevance for mental health

Background

Allostatic load is the cumulative "wear and tear" on the body due to chronic adversity. We aimed to test poly-environmental (exposomic) and polygenic contributions to allostatic load and their combined contribution to early adolescent mental health.

Methods

We analyzed data on N = 5,035 diverse youth (mean age 12) from the Adolescent Brain Cognitive Development Study (ABCD). Using dimensionality reduction method, we calculated and overall allostatic load score (AL) using body mass index [BMI], waist circumference, blood pressure, blood glycemia, blood cholesterol, and salivary DHEA. Childhood exposomic risk was quantified using multi-level environmental exposures before age 11. Genetic risk was quantified using polygenic risk scores (PRS) for metabolic system susceptibility (type 2 diabetes [T2D]) and stress-related psychiatric disease (major depressive disorder [MDD]). We used linear mixed effects models to test main, additive, and interactive effects of exposomic and polygenic risk (independent variables) on AL (dependent variable). Mediation models tested the mediating role of AL on the pathway from exposomic and polygenic risk to youth mental health. Models adjusted for demographics and genetic principal components.

Results

We observed disparities in AL with non-Hispanic White youth having significantly lower AL compared to Hispanic and Non-Hispanic Black youth. In the diverse sample, childhood exposomic burden was associated with AL in adolescence (beta=0.25, 95%CI 0.22-0.29, P<.001). In European ancestry participants (n=2,928), polygenic risk of both T2D and depression was associated with AL (T2D-PRS beta=0.11, 95%CI 0.07-0.14, P<.001; MDD-PRS beta=0.05, 95%CI 0.02-0.09, P=.003). Both polygenic scores showed significant interaction with exposomic risk such that, with greater polygenic risk, the association between exposome and AL was stronger. AL partly mediated the pathway to youth mental health from exposomic risk and from MDD-PRS, and fully mediated the pathway from T2D-PRS.

Conclusions

AL can be quantified in youth using anthropometric and biological measures and is mapped to exposomic and polygenic risk. Main and interactive environmental and genetic effects support a diathesis-stress model. Findings suggest that both environmental and genetic risk be considered when modeling stress-related health conditions.

Challenges in studying the interplay of genes and environment. A study of childhood financial distress moderating genetic predisposition for peak smoking

Smoking is one of the leading causes of preventable disease and death in the U.S., and it is strongly influenced both by genetic predisposition and childhood adversity. Using polygenic indices (PGIs) of predisposition to smoking, we evaluate whether childhood financial distress (CFD; a composite measure of financial adversity) moderates genetic risk in explaining peak-cigarette consumption in adulthood. Using the Health and Retirement Study (HRS), we find a substantial reduction in the relationship between genetic risk and peak smoking for those who did not suffer financial adversity in childhood. Among adult smokers who grew up in high-CFD households, a one standard deviation higher PGI is associated with 2.9 more cigarettes smoked per day at peak. By contrast, among smokers who grew up in low-CFD households, this gradient is reduced by 37 percent (or 1.1 fewer). These results are robust to controlling for a host of prime confounders. By contrast, we find no evidence of interactions between the PGI and typical measures of childhood SES such as parental education - a null result that we replicate in the Wisconsin Longitudinal Study (WLS) and the English Longitudinal Study of Aging (ELSA). This suggests the role of childhood financial distress in the relationship with peak smoking is distinct from that of low childhood SES, with high CFD potentially reflecting more acute distress than do measures of low childhood SES. Our evidence also suggests low childhood SES is a weaker proxy for acute distress, providing an alternative explanation for the childhood SES null result.

Cauchy combination methods for the detection of gene-environment interactions for rare variants related to quantitative phenotypes

The characterization of gene-environment interactions (GEIs) can provide detailed insights into the biological mechanisms underlying complex diseases. Despite recent interest in GEIs for rare variants, published GEI tests are underpowered for an extremely small proportion of causal rare variants in a gene or a region. By extending the aggregated Cauchy association test (ACAT), we propose three GEI tests to address this issue: a Cauchy combination GEI test with fixed main effects (CCGEI-F), a Cauchy combination GEI test with random main effects (CCGEI-R), and an omnibus Cauchy combination GEI test (CCGEI-O). ACAT was applied to combine p values of single-variant GEI analyses to obtain CCGEI-F and CCGEI-R and p values of multiple GEI tests were combined in CCGEI-O. Through numerical simulations, for small numbers of causal variants, CCGEI-F, CCGEI-R and CCGEI-O provided approximately 5% higher power than the existing GEI tests INT-FIX and INT-RAN; however, they had slightly higher power than the existing GEI test TOW-GE. For large numbers of causal variants, although CCGEI-F and CCGEI-R exhibited comparable or slightly lower power values than the competing tests, the results were still satisfactory. Among all simulation conditions evaluated, CCGEI-O provided significantly higher power than that of competing GEI tests. We further applied our GEI tests in genome-wide analyses of systolic blood pressure or diastolic blood pressure to detect gene-body mass index (BMI) interactions, using whole-exome sequencing data from UK Biobank. At a suggestive significance level of 1.0 × 10-4, KCNC4, GAR1, FAM120AOS and NT5C3B showed interactions with BMI by our GEI tests.

Plant-based dietary patterns, genetic predisposition and risk of colorectal cancer: a prospective study from the UK Biobank

BACKGROUND

Plant-based dietary patterns may affect colorectal cancer (CRC) related outcomes, while risks differ in the quality of plant foods. We aimed to examine the association of plant-based diet quality with risks of CRC incidence and mortality and whether this association was modified by genetic risk.

METHODS

This prospective cohort study included 186,675 participants free of cancer when the last dietary recall was completed. We calculated three plant-based diet indices (PDIs), i.e., the overall plant-based diet index (PDI), the healthful plant-based diet index (hPDI), and the unhealthful plant-based diet index (uPDI) representing adherence to plant-based diets with diverse quality. Genetic risk was characterized using a weighted polygenic risk score (PRS), capturing overall risk variants associated with CRC. Hazard ratios (HR) and 95% confidential intervals (CI) were estimated by the cause-specific Cox proportional hazards model.

RESULTS

Over a follow-up of 9.5 years, 2163 cases and 466 deaths from CRC were documented. The HR of CRC incidence was 0.88 (95% CI, 0.81-0.96) and 0.91 (95% CI, 0.84-0.99) per 10-score increase in PDI and hPDI, respectively. Compared to the lowest quartile, PDI, hPDI, and uPDI in the highest quartile were associated with a 13% decrease, a 15% decrease, and a 14% increase in risk of incident CRC, respectively. We found a joint association of genetic risk and PDIs with incident CRC, with the highest hazard observed in those carrying higher PRS and adhering to lower-quality PDIs. The inverse association of PDI and hPDI with CRC mortality was pronounced in males.

CONCLUSIONS

Our results suggested that better adherence to overall and healthful plant-based diets was associated with a lower risk of CRC, whereas an unhealthful plant-based diet was associated with a higher CRC risk. Consumption of a higher-quality plant-based diet combined with decreased genetic risk conferred less susceptibility to CRC. Our findings highlighted the importance of food quality when adhering to a plant-based dietary pattern for CRC prevention in the general population.

Genetic risk assessment based on association and prediction studies

The genetic basis of phenotypic emergence provides valuable information for assessing individual risk. While association studies have been pivotal in identifying genetic risk factors within a population, complementing it with insights derived from predictions studies that assess individual-level risk offers a more comprehensive approach to understanding phenotypic expression. In this study, we established personalized risk assessment models using single-nucleotide polymorphism (SNP) data from 200 Korean patients, of which 100 experienced hepatitis B surface antigen (HBsAg) seroclearance and 100 patients demonstrated high levels of HBsAg. The risk assessment models determined the predictive power of the following: (1) genome-wide association study (GWAS)-identified candidate biomarkers considered significant in a reference study and (2) machine learning (ML)-identified candidate biomarkers with the highest feature importance scores obtained by using random forest (RF). While utilizing all features yielded 64% model accuracy, using relevant biomarkers achieved higher model accuracies: 82% for 52 GWAS-identified candidate biomarkers, 71% for three GWAS-identified biomarkers, and 80% for 150 ML-identified candidate biomarkers. Findings highlight that the joint contributions of relevant biomarkers significantly influence phenotypic emergence. On the other hand, combining ML-identified candidate biomarkers into the pool of GWAS-identified candidate biomarkers resulted in the improved predictive accuracy of 90%, demonstrating the capability of ML as an auxiliary analysis to GWAS. Furthermore, some of the ML-identified candidate biomarkers were found to be linked with hepatocellular carcinoma (HCC), reinforcing previous claims that HCC can still occur despite the absence of HBsAg.

Applying an evolutionary mismatch framework to understand disease susceptibility

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of "lifestyle" diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be "mismatched" and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions, with different health effects in "ancestral" versus "modern" environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the "matched" to "mismatched" spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

Misclassification Bias in the Assessment of Gene-by-Environment Interactions

BACKGROUND

Misclassification bias is a common concern in epidemiologic studies. Despite strong bias on main effects, gene-environment interactions have been shown to be biased towards the null under gene-environment independence. In the context of a recent article examining the interaction between nerve agent exposure and paraoxonase-1 gene on Gulf War Illness, we aimed to assess the impact of recall bias-a common misclassfication bias-on the identification of gene-environment interactions when the independence assumption is violated.

METHODS

We derive equations to quantify the bias of the interaction, and numerically illustrate these results by simulating a case-control study of 1000 cases and 1000 controls. Simulation input parameters included exposure prevalence, strength of gene-environment dependence, strength of the main effect, exposure specificity among cases, and strength of the gene-environment interaction.

RESULTS

We show that, even if gene-environment independence is violated, we can bound possible gene-environment interactions by knowing the strength and direction of the gene-environment dependence ( ) and the observed gene-environment interaction ( )-thus often still allowing for the identification of such interactions. Depending on whether is larger or smaller than the inverse of , is a lower (if ) or upper (if ) bound for the true interaction. In addition, the bias magnitude is somewhat predictable by examining other characteristics such as exposure prevalence, the strength of the exposure main effect, and directions of the recall bias and gene-environment dependence.

CONCLUSIONS

Even if gene-environment dependence exists, we may still be able to identify gene-environment interactions even when misclassification bias is present.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Independent and Combined Associations of Sleep Duration, Bedtime, and Polygenic Risk Score with the Risk of Hearing Loss among Middle-Aged and Old Chinese: The Dongfeng-Tongji Cohort Study

Evidence available on the independent and combined associations of sleep duration, bedtime, and genetic predisposition with hearing loss was lacking. The present study included 15,827 participants from the Dongfeng-Tongji cohort study. Genetic risk was characterized by polygenic risk score (PRS) based on 37 genetic loci related to hearing loss. We conducted multivariate logistic regression models to assess the odds ratio (OR) for hearing loss with sleep duration and bedtime, as well as the joint association and interaction with PRS. Results showed that hearing loss was independently associated with sleeping ≥9 h/night compared to the recommended 7 to <8 h/night, and with bedtime ≤9:00 p.m. and >9:00 p.m. to 10:00 p.m. compared to those with bedtime >10:00 p.m. to 11:00 p.m., with estimated ORs of 1.25, 1.27, and 1.16, respectively. Meanwhile, the risk of hearing loss increased by 29% for each 5-risk allele increment of PRS. More importantly, joint analyses showed that the risk of hearing loss was 2-fold in sleep duration ≥9 h/night and high PRS, and 2.18-fold in bedtime ≤9:00 p.m. and high PRS. With significant joint effects of sleep duration and bedtime on hearing loss, we found an interaction of sleep duration with PRS in those with early bedtime and an interaction of bedtime with PRS in those with long sleep duration on hearing loss (P_int <0.05), and such relationships were more evident in high PRS. Similarly, the above relationships were also observed for age-related hearing loss and noise-induced hearing loss, particularly the latter. In addition, age-modified effects of sleep patterns on hearing loss were likewise observed, with stronger estimation among those aged <65 years. Accordingly, longer sleep duration, early bedtime, and high PRS were independently and jointly related to increased risk of hearing loss, suggesting the importance of considering both genetics and sleep pattern for risk assessment of hearing loss.

Gene-environment interaction analysis via deep learning

Gene-environment (G-E) interaction analysis plays an important role in studying complex diseases. Extensive methodological research has been conducted on G-E interaction analysis, and the existing methods are mostly based on regression techniques. In many fields including biomedicine and omics, it has been increasingly recognized that deep learning may outperform regression with its unique flexibility (e.g., in accommodating unspecified nonlinear effects) and superior prediction performance. However, there has been a lack of development in deep learning for G-E interaction analysis. In this article, we fill this important knowledge gap and develop a new analysis approach based on deep neural network in conjunction with penalization. The proposed approach can simultaneously conduct model estimation and selection (of important main G effects and G-E interactions), while uniquely respecting the "main effects, interactions" variable selection hierarchy. Simulation shows that it has superior prediction and feature selection performance. The analysis of data on lung adenocarcinoma and skin cutaneous melanoma overall survival further establishes its practical utility. Overall, this study can advance G-E interaction analysis by delivering a powerful new analysis approach based on modern deep learning.

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 necessity of incorporating non-genetic risk factors into polygenic risk score models

The growing public interest in genetic risk scores for various health conditions can be harnessed to inspire preventive health action. However, current commercially available genetic risk scores can be deceiving as they do not consider other, easily attainable risk factors, such as sex, BMI, age, smoking habits, parental disease status and physical activity. Recent scientific literature shows that adding these factors can improve PGS based predictions significantly. However, implementation of existing PGS based models that also consider these factors requires reference data based on a specific genotyping chip, which is not always available. In this paper, we offer a method naïve to the genotyping chip used. We train these models using the UK Biobank data and test these externally in the Lifelines cohort. We show improved performance at identifying the 10% most at-risk individuals for type 2 diabetes (T2D) and coronary artery disease (CAD) by including common risk factors. Incidence in the highest risk group increases from 3.0- and 4.0-fold to 5.8 for T2D, when comparing the genetics-based model, common risk factor-based model and combined model, respectively. Similarly, we observe an increase from 2.4- and 3.0-fold to 4.7-fold risk for CAD. As such, we conclude that it is paramount that these additional variables are considered when reporting risk, unlike current practice with current available genetic tests.

Evolutionary mismatch and the role of GxE interactions in human disease

Globally, we are witnessing the rise of complex, non-communicable diseases (NCDs) related to changes in our daily environments. Obesity, asthma, cardiovascular disease, and type 2 diabetes are part of a long list of "lifestyle" diseases that were rare throughout human history but are now common. A key idea from anthropology and evolutionary biology-the evolutionary mismatch hypothesis-seeks to explain this phenomenon. It posits that humans evolved in environments that radically differ from the ones experienced by most people today, and thus traits that were advantageous in past environments may now be "mismatched" and disease-causing. This hypothesis is, at its core, a genetic one: it predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions and have differential health effects in ancestral versus modern environments. Here, we discuss how this concept could be leveraged to uncover the genetic architecture of NCDs in a principled way. Specifically, we advocate for partnering with small-scale, subsistence-level groups that are currently transitioning from environments that are arguably more "matched" with their recent evolutionary history to those that are more "mismatched". These populations provide diverse genetic backgrounds as well as the needed levels and types of environmental variation necessary for mapping GxE interactions in an explicit mismatch framework. Such work would make important contributions to our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and sociocultural contexts.

Gene-environment interactions and their impact on human health

The molecular processes underlying human health and disease are highly complex. Often, genetic and environmental factors contribute to a given disease or phenotype in a non-additive manner, yielding a gene-environment (G × E) interaction. In this work, we broadly review current knowledge on the impact of gene-environment interactions on human health. We first explain the independent impact of genetic variation and the environment. We next detail well-established G × E interactions that impact human health involving environmental toxicants, pollution, viruses, and sex chromosome composition. We conclude with possibilities and challenges for studying G × E interactions.

Utilizing a Population-Genetic Framework to Test for Gene-Environment Interactions between Zebrafish Behavior and Chemical Exposure

Individuals within genetically diverse populations display broad susceptibility differences upon chemical exposures. Understanding the role of gene-environment interactions (GxE) in differential susceptibility to an expanding exposome is key to protecting public health. However, a chemical's potential to elicit GxE is often not considered during risk assessment. Previously, we've leveraged high-throughput zebrafish (Danio rerio) morphology screening data to reveal patterns of potential GxE effects. Here, using a population genetics framework, we apportioned variation in larval behavior and gene expression in three different PFHxA environments via mixed-effect modeling to assess significance of GxE term. We estimated the intraclass correlation (ICC) between full siblings from different families using one-way random-effects model. We found a significant GxE effect upon PFHxA exposure in larval behavior, and the ICC of behavioral responses in the PFHxA exposed population at the lower concentration was 43.7%, while that of the control population was 14.6%. Considering global gene expression data, a total of 3746 genes showed statistically significant GxE. By showing evidence that heritable genetics are directly affecting gene expression and behavioral susceptibility of individuals to PFHxA exposure, we demonstrate how standing genetic variation in a heterogeneous population such as ours can be leveraged to test for potential GxE.

Toolkit article: Approaches to measuring social inequities in health in human biology research

Across populations, human morbidity and mortality risks generally follow clear gradients, with socially-disadvantaged individuals and groups tending to have higher morbidity and mortality at all life stages relative to those more socially advantaged. Anthropologists specialize in understanding the proximate and ultimate factors that shape variation in human biological functioning and health and are therefore well-situated to explore the relationships between social position and health in diverse ecological and cultural contexts. While human biologists have developed sophisticated methods for assessing health using minimally-invasive methods, at a disciplinary level, we have room for conceptual and methodological improvement in how we frame, measure, and analyze the social inequities that might shape health inequities. This toolkit paper elaborates on some steps human biologists should take to enhance the quality of our research on health inequities. Specifically, we address: (1) how to frame unequal health outcomes (i.e., inequalities vs. disparities vs. inequities) and the importance of identifying our conceptual models of how these inequities emerge; (2) how to measure various axes of social inequities across diverse cultural contexts, and (3) approaches to community collaboration and dissemination. We end by discussing (4) future directions in human biology research of health inequities, including understanding the ultimate causes of sensitivity to social inequities and transitioning from research to action.

Vitamin C attenuates predisposition to high-fat diet-induced metabolic dysregulation in GLUT10-deficient mouse model

Background

The development of type 2 diabetes mellitus (T2DM) is highly influenced by complex interactions between genetic and environmental (dietary and lifestyle) factors. While vitamin C (ascorbic acid, AA) has been suggested as a complementary nutritional treatment for T2DM, evidence for the significance and beneficial effects of AA in T2DM is thus far inconclusive. We suspect that clinical studies on the topic might need to account for combination of genetic and dietary factors that could influence AA effects on metabolism. In this study, we tested this general idea using a mouse model with genetic predisposition to diet-induced metabolic dysfunction. In particular, we utilized mice carrying a human orthologous GLUT10^G128E variant (GLUT10^G128E mice), which are highly sensitive to high-fat diet (HFD)-induced metabolic dysregulation. The genetic variant has high relevance to human populations, as genetic polymorphisms in glucose transporter 10 (GLUT10) are associated with a T2DM intermediate phenotype in nondiabetic population.

Results

We investigated the impacts of AA supplementation on metabolism in wild-type (WT) mice and GLUT10^G128E mice fed with a normal diet or HFD. Overall, the beneficial effects of AA on metabolism were greater in HFD-fed GLUT10^G128E mice than in HFD-fed WT mice. At early postnatal stages, AA improved the development of compromised epididymal white adipose tissue (eWAT) in GLUT10^G128E mice. In adult animals, AA supplementation attenuated the predisposition of GLUT10^G128E mice to HFD-triggered eWAT inflammation, adipokine dysregulation, ectopic fatty acid accumulation, metabolic dysregulation, and body weight gain, as compared with WT mice.

Conclusions

Taken together, our findings suggest that AA has greater beneficial effects on metabolism in HFD-fed GLUT10^G128E mice than HFD-fed WT mice. As such, AA plays an important role in supporting eWAT development and attenuating HFD-induced metabolic dysregulation in GLUT10^G128E mice. Our results suggest that proper WAT development is essential for metabolic regulation later in life. Furthermore, when considering the usage of AA as a complementary nutrition for prevention and treatment of T2DM, individual differences in genetics and dietary patterns should be taken into account.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12263-022-00713-y.

Multidimensional molecular measurements-environment interaction analysis for disease outcomes

Multiple types of molecular (genetic, genomic, epigenetic, etc.) measurements, environmental risk factors, and their interactions have been found to contribute to the outcomes and phenotypes of complex diseases. In each of the previous studies, only the interactions between one type of molecular measurement and environmental risk factors have been analyzed. In recent biomedical studies, multidimensional profiling, in which data from multiple types of molecular measurements are collected from the same subjects, is becoming popular. A myriad of recent studies have shown that collectively analyzing multiple types of molecular measurements is not only biologically sensible but also leads to improved estimation and prediction. In this study, we conduct an M-E interaction analysis, with M standing for multidimensional molecular measurements and E standing for environmental risk factors. This can accommodate multiple types of molecular measurements and sufficiently account for their overlapping as well as independent information. Extensive simulation shows that it outperforms several closely related alternatives. In the analysis of TCGA (The Cancer Genome Atlas) data on lung adenocarcinoma and cutaneous melanoma, we make some stable biological findings and achieve stable prediction.

Multi-omics analysis: Paving the path toward achieving precision medicine in cancer treatment and immuno-oncology

The acceleration of large-scale sequencing and the progress in high-throughput computational analyses, defined as omics, was a hallmark for the comprehension of the biological processes in human health and diseases. In cancerology, the omics approach, initiated by genomics and transcriptomics studies, has revealed an incredible complexity with unsuspected molecular diversity within a same tumor type as well as spatial and temporal heterogeneity of tumors. The integration of multiple biological layers of omics studies brought oncology to a new paradigm, from tumor site classification to pan-cancer molecular classification, offering new therapeutic opportunities for precision medicine. In this review, we will provide a comprehensive overview of the latest innovations for multi-omics integration in oncology and summarize the largest multi-omics dataset available for adult and pediatric cancers. We will present multi-omics techniques for characterizing cancer biology and show how multi-omics data can be combined with clinical data for the identification of prognostic and treatment-specific biomarkers, opening the way to personalized therapy. To conclude, we will detail the newest strategies for dissecting the tumor immune environment and host–tumor interaction. We will explore the advances in immunomics and microbiomics for biomarker identification to guide therapeutic decision in immuno-oncology.

Diverse environmental perturbations reveal the evolution and context-dependency of genetic effects on gene expression levels

There is increasing appreciation that, in addition to being shaped by an individual's genotype and environment, most complex traits are also determined by poorly understood interactions between these two factors. So-called "genotype × environment" (G×E) interactions remain difficult to map at the organismal level but can be uncovered using molecular phenotypes. To do so at large scale, we used TM3'seq to profile transcriptomes across 12 cellular environments in 544 immortalized B cell lines from the 1000 Genomes Project. We mapped the genetic basis of gene expression levels across environments and revealed a context-dependent genetic architecture: The average heritability of gene expression levels increased in treatment relative to control conditions, and on average, each treatment revealed new expression quantitative trait loci (eQTLs) at 11% of genes. Across our experiments, 22% of all identified eQTLs were context-dependent, and this group was enriched for trait- and disease-associated loci. Further, evolutionary analyses suggested that positive selection has shaped G×E loci involved in responding to immune challenges and hormones but not to man-made chemicals. We hypothesize that this reflects a reduced opportunity for selection to act on responses to molecules recently introduced into human environments. Together, our work highlights the importance of considering an exposure's evolutionary history when studying and interpreting G×E interactions, and provides new insight into the evolutionary mechanisms that maintain G×E loci in human populations.

Schizophrenia Polygenic Risk and Experiences of Childhood Adversity: A Systematic Review and Meta-analysis

BACKGROUND AND HYPOTHESIS

Schizophrenia has been robustly associated with multiple genetic and environmental risk factors. Childhood adversity is one of the most widely replicated environmental risk factors for schizophrenia, but it is unclear if schizophrenia genetic risk alleles contribute to this association.

STUDY DESIGN

In this systematic review and meta-analysis, we assessed the evidence for gene-environment correlation (genes influence likelihood of environmental exposure) between schizophrenia polygenic risk score (PRS) and reported childhood adversity. We also assessed the evidence for a gene-environment interaction (genes influence sensitivity to environmental exposure) in relation to the outcome of schizophrenia and/or psychosis. This study was registered on PROSPERO (CRD42020182812). Following PRISMA guidelines, a search for relevant literature was conducted using Cochrane, MEDLINE, PsycINFO, Web of Science, and Scopus databases until February 2022. All studies that examined the association between schizophrenia PRS and childhood adversity were included.

STUDY RESULTS

Seventeen of 650 identified studies met the inclusion criteria and were assessed against the Newcastle-Ottawa Scale for quality. The meta-analysis found evidence for gene-environment correlation between schizophrenia PRS and childhood adversity (r = .02; 95% CI = 0.01, 0.03; P = .001), but the effect was small and therefore likely to explain only a small proportion of the association between childhood adversity and psychosis. The 4 studies that investigated a gene-environment interaction between schizophrenia PRS and childhood adversity in increasing risk of psychosis reported inconsistent results.

CONCLUSIONS

These findings suggest that a gene-environment correlation could explain a small proportion of the relationship between reported childhood adversity and psychosis.

Gene-environment interactions explain a substantial portion of variability of common neuropsychiatric disorders

In complex diseases, the phenotypic variability can be explained by genetic variation (G), environmental stimuli (E), and interaction of genetic and environmental factors (G-by-E effects), among which the contribution G-by-E remains largely unknown. In this study, we focus on ten major neuropsychiatric disorders using data for 138,383 United States families with 404,475 unique individuals. We show that, while gene-environment interactions account for only a small portion of the total phenotypic variance for a subset of disorders (depression, adjustment disorder, substance abuse), they explain a rather large portion of the phenotypic variation of the remaining disorders: over 20% for migraine and close to or over 30% for anxiety/phobic disorder, attention-deficit/hyperactivity disorder, recurrent headaches, sleep disorders, and post-traumatic stress disorder. In this study, we have incorporated—in the same analysis—clinical data, family pedigrees, the spatial distribution of individuals, their socioeconomic and demographic confounders, and a collection of environmental measurements.

•

G-by-E interactions contribute to the etiology of neuropsychiatric disorders

•

Environmental quality is associated with the risk of neuropsychiatric disorders

•

Geographic locations explain a part of variability of neuropsychiatric disorders

•

Identifying specific G-by-E interactions might be the future of personalized medicine

Substance use and common contributors to morbidity: A genetics perspective

Excessive substance use and substance use disorders (SUDs) are common, serious and relapsing medical conditions. They frequently co-occur with other diseases that are leading contributors to disability worldwide. While heavy substance use may potentiate the course of some of these illnesses, there is accumulating evidence suggesting common genetic architectures. In this narrative review, we focus on four heritable medical conditions - cardiometabolic disease, chronic pain, depression and COVID-19, which are commonly overlapping with, but not necessarily a direct consequence of, SUDs. We find persuasive evidence of underlying genetic liability that predisposes to both SUDs and chronic pain, depression, and COVID-19. For cardiometabolic disease, there is greater support for a potential causal influence of problematic substance use. Our review encourages de-stigmatization of SUDs and the assessment of substance use in clinical settings. We assert that identifying shared pathways of risk has high translational potential, allowing tailoring of treatments for multiple medical conditions.

Funding

SSR acknowledges T29KT0526, T32IR5226 and DP1DA054394; RLK acknowledges AA028292; AA acknowledges DA054869 & K02DA032573. The funders had no role in the conceptualization or writing of the paper.

Gene-diet interactions and cardiovascular diseases: a systematic review of observational and clinical trials

BACKGROUND

Both genetic background and diet are important determinants of cardiovascular diseases (CVD). Understanding gene-diet interactions could help improve CVD prevention and prognosis. We aimed to summarise the evidence on gene-diet interactions and CVD outcomes systematically.

METHODS

We searched MEDLINE^® via Ovid, Embase, PubMed^®, and The Cochrane Library for relevant studies published until June 6th 2022. We considered for inclusion cross-sectional, case-control, prospective cohort, nested case-control, and case-cohort studies as well as randomised controlled trials that evaluated the interaction between genetic variants and/or genetic risk scores and food or diet intake on the risk of related outcomes, including myocardial infarction, coronary heart disease (CHD), stroke and CVD as a composite outcome. The PROSPERO protocol registration code is CRD42019147031.

RESULTS AND DISCUSSION

We included 59 articles based on data from 29 studies; six articles involved multiple studies, and seven did not report details of their source population. The median sample size of the articles was 2562 participants. Of the 59 articles, 21 (35.6%) were qualified as high quality, while the rest were intermediate or poor. Eleven (18.6%) articles adjusted for multiple comparisons, four (7.0%) attempted to replicate the findings, 18 (30.5%) were based on Han-Chinese ethnicity, and 29 (49.2%) did not present Minor Allele Frequency. Fifty different dietary exposures and 52 different genetic factors were investigated, with alcohol intake and ADH1C variants being the most examined. Of 266 investigated diet-gene interaction tests, 50 (18.8%) were statistically significant, including CETP-TaqIB and ADH1C variants, which interacted with alcohol intake on CHD risk. However, interactions effects were significant only in some articles and did not agree on the direction of effects. Moreover, most of the studies that reported significant interactions lacked replication. Overall, the evidence on gene-diet interactions on CVD is limited, and lack correction for multiple testing, replication and sample size consideration.

Accounting for diet and age

The diet and age of mice can modulate how different genetic variants impact body weight, demonstrating the need to take context into account when performing genetic studies.

Cancer as microenvironmental, systemic and environmental diseases: opportunity for transdisciplinary microbiomics science

Cancer is generally regarded as a localised disease, with the well-established role of the tumour microenvironment. However, the realm of cancer goes beyond the tumour microenvironment, and cancer should also be regarded as a systemic and environmental disease. The exposome (ie, the totality of exposures), which encompasses diets, supplements, smoking, alcohol, other lifestyle factors, medications, etc, likely alters the microbiome (inclusive of bacteria, viruses, archaea, fungi, parasites, etc) and immune system in various body sites and influences tumour phenotypes. The systemic metabolic/inflammatory status, which is likely influenced by exposures and intestinal physiological changes, may affect tissue microenvironment of colorectum and any other organs. Germline genomic factors can modify disease phenotypes via gene-by-environment interactions. Although challenges exist, it is crucial to advance not only basic experimental research that can analyse the effects of exposures, microorganisms and microenvironmental components on tumour evolution but also interdisciplinary human population research that can dissect the complex pathogenic roles of the exposome, microbiome and immunome. Metagenomic, metatranscriptomic and metabolomic analyses should be integrated into well-designed population research combined with advanced methodologies of artificial intelligence and molecular pathological epidemiology. Ideally, a prospective cohort study design that enables biospecimen (such as stool) collection before disease detection should be considered to address reverse causation and recall biases. Robust experimental and observational research together can provide insights into dynamic interactions between environmental exposures, microbiota, tumour and immunity during carcinogenesis processes, thereby helping us develop precision prevention and therapeutic strategies to ultimately reduce the cancer burden.

Open problems in human trait genetics

Genetic studies of human traits have revolutionized our understanding of the variation between individuals, and yet, the genetics of most traits is still poorly understood. In this review, we highlight the major open problems that need to be solved, and by discussing these challenges provide a primer to the field. We cover general issues such as population structure, epistasis and gene-environment interactions, data-related issues such as ancestry diversity and rare genetic variants, and specific challenges related to heritability estimates, genetic association studies, and polygenic risk scores. We emphasize the interconnectedness of these problems and suggest promising avenues to address them.

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.

Social Determinants of Health Factors for Gene-Environment COVID-19 Research: Challenges and Opportunities

The characteristics of a person's health status are often guided by how they live, grow, learn, their genetics, as well as their access to health care. Yet, all too often, studies examining the relationship between social determinants of health (behavioral, sociocultural, and physical environmental factors), the role of demographics, and health outcomes poorly represent these relationships, leading to misinterpretations, limited study reproducibility, and datasets with limited representativeness and secondary research use capacity. This is a profound hurdle in what questions can or cannot be rigorously studied about COVID-19. In practice, gene-environment interactions studies have paved the way for including these factors into research. Similarly, our understanding of social determinants of health continues to expand with diverse data collection modalities as health systems, patients, and community health engagement aim to fill the knowledge gaps toward promoting health and wellness. Here, a conceptual framework is proposed, adapted from the population health framework, socioecological model, and causal modeling in gene-environment interaction studies to integrate the core constructs from each domain with practical considerations needed for multidisciplinary science.

Overlapping group screening for detection of gene-environment interactions with application to TCGA high-dimensional survival genomic data

BACKGROUND

In the context of biomedical and epidemiological research, gene-environment (G-E) interaction is of great significance to the etiology and progression of many complex diseases. In high-dimensional genetic data, two general models, marginal and joint models, are proposed to identify important interaction factors. Most existing approaches for identifying G-E interactions are limited owing to the lack of robustness to outliers/contamination in response and predictor data. In particular, right-censored survival outcomes make the associated feature screening even challenging. In this article, we utilize the overlapping group screening (OGS) approach to select important G-E interactions related to clinical survival outcomes by incorporating the gene pathway information under a joint modeling framework.

RESULTS

Simulation studies under various scenarios are carried out to compare the performances of our proposed method with some commonly used methods. In the real data applications, we use our proposed method to identify G-E interactions related to the clinical survival outcomes of patients with head and neck squamous cell carcinoma, and esophageal carcinoma in The Cancer Genome Atlas clinical survival genetic data, and further establish corresponding survival prediction models. Both simulation and real data studies show that our method performs well and outperforms existing methods in the G-E interaction selection, effect estimation, and survival prediction accuracy.

CONCLUSIONS

The OGS approach is useful for selecting important environmental factors, genes and G-E interactions in the ultra-high dimensional feature space. The prediction ability of OGS with the Lasso penalty is better than existing methods. The same idea of the OGS approach can apply to other outcome models, such as the proportional odds survival time model, the logistic regression model for binary outcomes, and the multinomial logistic regression model for multi-class outcomes.

GEInfo: an R package for gene-environment interaction analysis incorporating prior information

SUMMARY

Gene-environment (G-E) interactions have important implications for many complex diseases. With higher dimensionality and weaker signals, G-E interaction analysis is more challenged than the analysis of main G (and E) effects. The accumulation of published literature makes it possible to borrow strength from prior information and improve analysis. In a recent study, a 'quasi-likelihood + penalization' approach was developed to effectively incorporate prior information. Here, we first extend it to linear, logistic and Poisson regressions. Such models are much more popular in practice. More importantly, we develop the R package GEInfo, which realizes this approach in a user-friendly manner. To facilitate direct comparison and routine data analysis, the package also includes functions for alternative methods and visualization.

AVAILABILITY AND IMPLEMENTATION

The package is available at https://CRAN.R-project.org/package=GEInfo.

SUPPLEMENTARY INFORMATION

Supplementary materials are available at Bioinformatics online.

Increased genetic contribution to wellbeing during the COVID-19 pandemic

Physical and mental health are determined by an interplay between nature, for example genetics, and nurture, which encompasses experiences and exposures that can be short or long-lasting. The COVID-19 pandemic represents a unique situation in which whole communities were suddenly and simultaneously exposed to both the virus and the societal changes required to combat the virus. We studied 27,537 population-based biobank participants for whom we have genetic data and extensive longitudinal data collected via 19 questionnaires over 10 months, starting in March 2020. This allowed us to explore the interaction between genetics and the impact of the COVID-19 pandemic on individuals' wellbeing over time. We observe that genetics affected many aspects of wellbeing, but also that its impact on several phenotypes changed over time. Over the course of the pandemic, we observed that the genetic predisposition to life satisfaction had an increasing influence on perceived quality of life. We also estimated heritability and the proportion of variance explained by shared environment using variance components methods based on pedigree information and household composition. The results suggest that people's genetic constitution manifested more prominently over time, potentially due to social isolation driven by strict COVID-19 containment measures. Overall, our findings demonstrate that the relative contribution of genetic variation to complex phenotypes is dynamic rather than static.

Leveraging a High-Throughput Screening Method to Identify Mechanisms of Individual Susceptibility Differences in a Genetically Diverse Zebrafish Model

Understanding the mechanisms behind chemical susceptibility differences is key to protecting sensitive populations. However, elucidating gene-environment interactions (GxE) presents a daunting challenge. While mammalian models have proven useful, problems with scalability to an enormous chemical exposome and clinical translation faced by all models remain; therefore, alternatives are needed. Zebrafish (Danio rerio) have emerged as an excellent model for investigating GxE. This study used a combined bioinformatic and experimental approach to probe the mechanisms underlying chemical susceptibility differences in a genetically diverse zebrafish population. Starting from high-throughput screening (HTS) data, a genome-wide association study (GWAS) using embryonic fish exposed to 0.6 μM Abamectin revealed significantly different effects between individuals. A hypervariable region with two distinct alleles–one with G at the SNP locus (GG) and one with a T and the 16 bp deletion (TT)–associated with differential susceptibility was found. Sensitive fish had significantly lower sox7 expression. Due to their location and the observed expression differences, we hypothesized that these sequences differentially regulate sox7. A luciferase reporter gene assay was used to test if these sequences, alone, could lead to expression differences. The TT allele showed significantly lower expression than the GG allele in MCF-7 cells. To better understand the mechanism behind these expression differences, predicted transcription factor binding differences between individuals were compared in silico, and several putative binding differences were identified. EMSA was used to test for binding differences in whole embryo protein lysate to investigate these TF binding predictions. We confirmed that the GG sequence is bound to protein in zebrafish. Through a competition EMSA using an untagged oligo titration, we confirmed that the GG oligo had a higher binding affinity than the TT oligo, explaining the observed expression differences. This study identified differential susceptibility to chemical exposure in a genetically diverse population, then identified a plausible mechanism behind those differences from a genetic to molecular level. Thus, an HTS-compatible zebrafish model is valuable and adaptable in identifying GxE mechanisms behind susceptibility differences to chemical exposure.

Polygenic scores, diet quality, and type 2 diabetes risk: An observational study among 35,759 adults from 3 US cohorts

BACKGROUND

Both genetic and lifestyle factors contribute to the risk of type 2 diabetes, but the extent to which there is a synergistic effect of the 2 factors is unclear. The aim of this study was to examine the joint associations of genetic risk and diet quality with incident type 2 diabetes.

METHODS AND FINDINGS

We analyzed data from 35,759 men and women in the United States participating in the Nurses' Health Study (NHS) I (1986 to 2016) and II (1991 to 2017) and the Health Professionals Follow-up Study (HPFS; 1986 to 2016) with available genetic data and who did not have diabetes, cardiovascular disease, or cancer at baseline. Genetic risk was characterized using both a global polygenic score capturing overall genetic risk and pathway-specific polygenic scores denoting distinct pathophysiological mechanisms. Diet quality was assessed using the Alternate Healthy Eating Index (AHEI). Cox models were used to calculate hazard ratios (HRs) for type 2 diabetes after adjusting for potential confounders. With over 902,386 person-years of follow-up, 4,433 participants were diagnosed with type 2 diabetes. The relative risk of type 2 diabetes was 1.29 (95% confidence interval [CI] 1.25, 1.32; P < 0.001) per standard deviation (SD) increase in global polygenic score and 1.13 (1.09, 1.17; P < 0.001) per 10-unit decrease in AHEI. Irrespective of genetic risk, low diet quality, as compared to high diet quality, was associated with approximately 30% increased risk of type 2 diabetes (Pinteraction = 0.69). The joint association of low diet quality and increased genetic risk was similar to the sum of the risk associated with each factor alone (Pinteraction = 0.30). Limitations of this study include the self-report of diet information and possible bias resulting from inclusion of highly educated participants with available genetic data.

CONCLUSIONS

These data provide evidence for the independent associations of genetic risk and diet quality with incident type 2 diabetes and suggest that a healthy diet is associated with lower diabetes risk across all levels of genetic risk.