Genotype-covariate interaction effects and the heritability of adult body mass index

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

Development and validation of an intrinsic capacity score in the UK Biobank study

BACKGROUND

In 2015, the World Health Organization introduced the concept of intrinsic capacity (IC) to define the individual-level characteristics that enable an older person to be and do the things they value. This study developed an intrinsic capacity score for UK Biobank study participants and validated its use as a tool for health outcome prediction, understanding healthy aging trajectories, and genetic research.

METHODS

Our analysis included data from 45,208 UK biobank participants who had a complete record of the ten variables included in the analysis. Factor adequacy was tested using Kaiser-Meyer-Olkin, Barthelt's, and the determinant of matrix tests, and the number of factors was determined by the parallel analysis method. Exploratory and confirmatory factor analyses were employed to determine the structure and dimensionality of indicators. Finally, the intrinsic capacity score was generated, and its construct and predictive validities as well as reliability were assessed.

RESULTS

The factor analysis identified a multidimensional construct comprising one general factor (intrinsic capacity) and five specific factors (locomotor, vitality, cognitive, psychological, and sensory). The bifactor structure showed a better fit (comparative fit index = 0.995, Tucker Lewis index = 0.976, root mean square error of approximation = 0.025, root mean square residual = 0.009) than the conventional five-factor structure. The intrinsic capacity score generated using the bifactor confirmatory factor analysis has good construct validity, as demonstrated by an inverse association with age (lower intrinsic capacity in older age; (β) =-0.035 (95%CI: -0.036, -0.034)), frailty (lower intrinsic capacity score in prefrail participants, β = -0.104 (95%CI: (-0.114, -0.094)) and frail participants, β = -0.227 (95%CI: -0.267, -0.186) than robust participants), and comorbidity (a lower intrinsic capacity score associated with increased Charlson's comorbidity index, β =-0.019 (95%CI: -0.022, -0.015)). The intrinsic capacity score also predicted comorbidity (a one-unit increase in baseline intrinsic capacity score led to a lower Charlson's comorbidity index, β = 0.147 (95%CI: -0.173, -0.121)) and mortality (a one-unit increase in baseline intrinsic capacity score led to 25 % lower risk of death, odds ratio = 0.75(95%CI: 0.663, 0.848)).

CONCLUSION

The bifactor structure showed a better fit in all goodness of fit tests. The intrinsic capacity construct has strong structural, construct, and predictive validities and is a promising tool for monitoring aging trajectories.

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.

Mutational landscape of risk variants in comorbid depression and obesity: a next-generation sequencing approach

Major depression (MD) and obesity are complex genetic disorders that are frequently comorbid. However, the study of both diseases concurrently remains poorly addressed and therefore the underlying genetic mechanisms involved in this comorbidity remain largely unknown. Here we examine the contribution of common and rare variants to this comorbidity through a next-generation sequencing (NGS) approach. Specific genomic regions of interest in MD and obesity were sequenced in a group of 654 individuals from the PISMA-ep epidemiological study. We obtained variants across the entire frequency spectrum and assessed their association with comorbid MD and obesity, both at variant and gene levels. We identified 55 independent common variants and a burden of rare variants in 4 genes (PARK2, FGF21, HIST1H3D and RSRC1) associated with the comorbid phenotype. Follow-up analyses revealed significantly enriched gene-sets associated with biological processes and pathways involved in metabolic dysregulation, hormone signaling and cell cycle regulation. Our results suggest that, while risk variants specific to the comorbid phenotype have been identified, the genes functionally impacted by the risk variants share cell biological processes and signaling pathways with MD and obesity phenotypes separately. To the best of our knowledge, this is the first study involving a targeted sequencing approach toward the study of the comorbid MD and obesity. The framework presented here allowed a deep characterization of the genetics of the co-occurring MD and obesity, revealing insights into the mutational and functional profile that underlies this comorbidity and contributing to a better understanding of the relationship between these two disabling disorders.

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.

Genetic influence on within-person longitudinal change in anthropometric traits in the UK Biobank

The causes of temporal fluctuations in adult traits are poorly understood. Here, we investigate the genetic determinants of within-person trait variability of 8 repeatedly measured anthropometric traits in 50,117 individuals from the UK Biobank. We found that within-person (non-directional) variability had a SNP-based heritability of 2-5% for height, sitting height, body mass index (BMI) and weight (P ≤ 2.4 × 10-3). We also analysed longitudinal trait change and show a loss of both average height and weight beyond about 70 years of age. A variant tracking the Alzheimer's risk APOE- E 4 allele (rs429358) was significantly associated with weight loss ( β  = -0.047 kg per yr, s.e. 0.007, P = 2.2 × 10-11), and using 2-sample Mendelian Randomisation we detected a relationship consistent with causality between decreased lumbar spine bone mineral density and height loss (bxy = 0.011, s.e. 0.003, P = 3.5 × 10-4). Finally, population-level variance quantitative trait loci (vQTL) were consistent with within-person variability for several traits, indicating an overlap between trait variability assessed at the population or individual level. Our findings help elucidate the genetic influence on trait-change within an individual and highlight disease risks associated with these changes.

Identifying latent genetic interactions in genome-wide association studies using multiple traits

The "missing" heritability of complex traits may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. We propose a new kernel-based method called Latent Interaction Testing (LIT) to screen for genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Using simulated data, we demonstrate that LIT increases power to detect latent genetic interactions compared to univariate methods. We then apply LIT to obesity-related traits in the UK Biobank and detect variants with interactive effects near known obesity-related genes (URL: https://CRAN.R-project.org/package=lit ).

Distinct explanations underlie gene-environment interactions in the UK Biobank

The role of gene-environment (GxE) interaction in disease and complex trait architectures is widely hypothesized, but currently unknown. Here, we apply three statistical approaches to quantify and distinguish three different types of GxE interaction for a given trait and E variable. First, we detect locus-specific GxE interaction by testing for genetic correlation r g < 1 across E bins. Second, we detect genome-wide effects of the E variable on genetic variance by leveraging polygenic risk scores (PRS) to test for significant PRSxE in a regression of phenotypes on PRS, E, and PRSxE, together with differences in SNP-heritability across E bins. Third, we detect genome-wide proportional amplification of genetic and environmental effects as a function of the E variable by testing for significant PRSxE with no differences in SNP-heritability across E bins. Simulations show that these approaches achieve high sensitivity and specificity in distinguishing these three GxE scenarios. We applied our framework to 33 UK Biobank traits (25 quantitative traits and 8 diseases; average N = 325 K ) and 10 E variables spanning lifestyle, diet, and other environmental exposures. First, we identified 19 trait-E pairs with r g significantly < 1 (FDR<5%) (average r g = 0.95 ); for example, white blood cell count had r g = 0.95 (s.e. 0.01) between smokers and non-smokers. Second, we identified 28 trait-E pairs with significant PRSxE and significant SNP-heritability differences across E bins; for example, BMI had a significant PRSxE for physical activity (P=4.6e-5) with 5% larger SNP-heritability in the largest versus smallest quintiles of physical activity (P=7e-4). Third, we identified 15 trait-E pairs with significant PRSxE with no SNP-heritability differences across E bins; for example, waist-hip ratio adjusted for BMI had a significant PRSxE effect for time spent watching television (P=5e-3) with no SNP-heritability differences. Across the three scenarios, 8 of the trait-E pairs involved disease traits, whose interpretation is complicated by scale effects. Analyses using biological sex as the E variable produced additional significant findings in each of the three scenarios. Overall, we infer a significant contribution of GxE and GxSex effects to complex trait and disease variance.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Risk factors affecting polygenic score performance across diverse cohorts

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed effects of covariate stratification and interaction on body mass index (BMI) PGS (PGSBMI) across four cohorts of European (N=491,111) and African (N=21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R2 being nearly double between best and worst performing quintiles for certain covariates. 28 covariates had significant PGSBMI-covariate interaction effects, modifying PGSBMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R2 differences among strata and interaction effects - across all covariates, their main effects on BMI were correlated with their maximum R2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGSBMI individuals have highest R2 and increase in PGS effect. Using quantile regression, we show the effect of PGSBMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGSBMI performance and effects, we investigated ways to increase model performance taking into account non-linear effects. Machine learning models (neural networks) increased relative model R2 (mean 23%) across datasets. Finally, creating PGSBMI directly from GxAge GWAS effects increased relative R2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGSBMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Genome-wide analysis reveals extensive genetic overlap between childhood phenotypes and later-life type 2 diabetes

Observational studies have indicated a potential influence of childhood phenotypes on the later development of type 2 diabetes (T2D). However, the underlying biological mechanisms remain unclear. In this study, we conducted a comprehensive genome-wide analysis to investigate the shared genetic architecture and genetic loci between nine childhood phenotypes (N = 4202-620,26) and later-life T2D (N = 80,154) using genetic correlation, mendelian randomization (MR), and conjunctional false discovery rate (conjFDR) statistical frameworks. Our findings demonstrated substantial genetic correlations and pleiotropic enrichment between childhood obesity, body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), and later-life T2D. Childhood obesity exhibited a significant association with increased later-life T2D risk through 10 mediators, 6 of which were adulthood obesity-related phenotypes. Additionally, we identified 69, 83, 3, 5, 10, 5, 3, and 7 loci shared between childhood obesity, BMI, SBP, high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), triglycerides (TG), apolipoprotein A-I (ApoA-I), apolipoprotein B (ApoB), and T2D at conjFDR <0.05, with the majority of these loci being novel discoveries. Overall, our study reveals extensive genetic overlap between childhood obesity-related phenotypes and T2D with concordant effect directions, shedding new light on variants and phenotypes with lifelong effects.

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.

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.

Identifying latent genetic interactions in genome-wide association studies using multiple traits

Genome-wide association studies of complex traits frequently find that SNP-based estimates of heritability are considerably smaller than estimates from classic family-based studies. This 'missing' heritability may be partly explained by genetic variants interacting with other genes or environments that are difficult to specify, observe, and detect. To circumvent these challenges, we propose a new method to detect genetic interactions that leverages pleiotropy from multiple related traits without requiring the interacting variable to be specified or observed. Our approach, Latent Interaction Testing (LIT), uses the observation that correlated traits with shared latent genetic interactions have trait variance and covariance patterns that differ by genotype. LIT examines the relationship between trait variance/covariance patterns and genotype using a flexible kernel-based framework that is computationally scalable for biobank-sized datasets with a large number of traits. We first use simulated data to demonstrate that LIT substantially increases power to detect latent genetic interactions compared to a trait-by-trait univariate method. We then apply LIT to four obesity-related traits in the UK Biobank and detect genetic variants with interactive effects near known obesity-related genes. Overall, we show that LIT, implemented in the R package lit, uses shared information across traits to improve detection of latent genetic interactions compared to standard approaches.

Genetic insights into the age-specific biological mechanisms governing human ovarian aging

There is currently little evidence that the genetic basis of human phenotype varies significantly across the lifespan. However, time-to-event phenotypes are understudied and can be thought of as reflecting an underlying hazard, which is unlikely to be constant through life when values take a broad range. Here, we find that 74% of 245 genome-wide significant genetic associations with age at natural menopause (ANM) in the UK Biobank show a form of age-specific effect. Nineteen of these replicated discoveries are identified only by our modeling framework, which determines the time dependency of DNA-variant age-at-onset associations without a significant multiple-testing burden. Across the range of early to late menopause, we find evidence for significantly different underlying biological pathways, changes in the signs of genetic correlations of ANM to health indicators and outcomes, and differences in inferred causal relationships. We find that DNA damage response processes only act to shape ovarian reserve and depletion for women of early ANM. Genetically mediated delays in ANM were associated with increased relative risk of breast cancer and leiomyoma at all ages and with high cholesterol and heart failure for late-ANM women. These findings suggest that a better understanding of the age dependency of genetic risk factor relationships among health indicators and outcomes is achievable through appropriate statistical modeling of large-scale biobank data.

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.

The Impact of Bariatric Surgery on Quality of Life in Patients with Obesity

Obesity has become a widespread health problem influencing people's health, general well-being, and healthcare costs. It also represents an important risk factor for multiple comorbidities and malignancies.

OBJECTIVES

the primary objective of this study was to provide notable insights to healthcare professionals regarding the management of patients with obesity, to highlight the effectiveness of bariatric surgical methods in losing excess weight, and to establish the relationship between weight loss and changes in quality of life (QoL).

MATERIAL AND METHODS

our study evaluated the QoL of 76 patients following bariatric surgery at different postoperative stages using the 36-Item Short Form Survey (SF-36) and The World Health Organization Quality of Life-BREF (WHOQOL-BREF) questionnaires.

RESULTS

regarding the type of bariatric procedure performed, out of the 76 respondents, 39.47% underwent gastric bypass surgery (RYGB), 56.57% underwent sleeve gastrectomy (LSG), and only 3.94% underwent single anastomosis duodeno-ileal switch (SADI-S). Pertaining to the SF-36 questionnaire, the lowest average scores were found in the energy/fatigue subscales and in the limitations due to mental health, which remained consistent across surgery types with a significant decrease in the SADI-S group. Concerning the WHOQOL-BREF questionnaire, the lowest mean scores were found in the environment (15.03 ± 2.37) and social relations (16.08 ± 2.22) subscales, whilst the highest average scores were in physical health (16.30 ± 2.03) and mental health (16.57 ± 2.16).

CONCLUSIONS

the findings revealed that whilst bariatric surgery significantly improved physical health, it resulted in a decrease in mental health scores. Consequently, the study emphasizes the importance of adopting a holistic approach to managing obesity that considers improving both physical and mental health outcomes.

Prevalence of overweight and obesity in patients with major depressive disorder with anxiety: Mediating role of thyroid hormones and metabolic parameters

BACKGROUNDS

Anxiety is a common comorbidity in major depressive disorder (MDD); however, its role in overweight and obesity in MDD patients remains unclear. We examined the relationship between severe anxiety and overweight and obesity, as well as the mediating role of thyroid hormones and metabolic parameters in MDD patients.

METHODS

This cross-sectional study recruited 1718 first-episode drug-naïve MDD outpatients. All participants were rated on the Hamilton Depression Rating Scale for depression and Hamilton Anxiety Rating Scale for anxiety and measured in thyroid hormones and metabolic parameters.

RESULTS

A total of 218 (12.7 %) individuals had severe anxiety. The prevalence of overweight and obesity in patients with severe anxiety was 62.8 % and 5.5 %, respectively. Severe anxiety symptoms were related to overweight (Odds Ratio [OR]: 1.47, 95 % CI: 1.08, 2.00) and obesity (OR: 2.10, 95 % CI: 1.07, 4.15). The association between severe anxiety and overweight was mainly attenuated by thyroid hormones (40.4 %), blood pressure (31.9 %), and plasma glucose (19.1 %). For obesity, the association with severe anxiety was mainly attenuated by thyroid hormones (48.2 %), blood pressure (39.1 %), and total cholesterol (28.2 %).

LIMITATIONS

Due to the cross-sectional design, no causal relationship could be derived.

CONCLUSIONS

Thyroid hormones and metabolic parameters can explain the risk of overweight and obesity associated with severe anxiety in MDD patients. These findings add to the knowledge of the pathological pathway of overweight and obesity in severe anxiety comorbid MDD.

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.

Phenotypic subtyping via contrastive learning

Defining and accounting for subphenotypic structure has the potential to increase statistical power and provide a deeper understanding of the heterogeneity in the molecular basis of complex disease. Existing phenotype subtyping methods primarily rely on clinically observed heterogeneity or metadata clustering. However, they generally tend to capture the dominant sources of variation in the data, which often originate from variation that is not descriptive of the mechanistic heterogeneity of the phenotype of interest; in fact, such dominant sources of variation, such as population structure or technical variation, are, in general, expected to be independent of subphenotypic structure. We instead aim to find a subspace with signal that is unique to a group of samples for which we believe that subphenotypic variation exists (e.g., cases of a disease). To that end, we introduce Phenotype Aware Components Analysis (PACA), a contrastive learning approach leveraging canonical correlation analysis to robustly capture weak sources of subphenotypic variation. In the context of disease, PACA learns a gradient of variation unique to cases in a given dataset, while leveraging control samples for accounting for variation and imbalances of biological and technical confounders between cases and controls. We evaluated PACA using an extensive simulation study, as well as on various subtyping tasks using genotypes, transcriptomics, and DNA methylation data. Our results provide multiple strong evidence that PACA allows us to robustly capture weak unknown variation of interest while being calibrated and well-powered, far superseding the performance of alternative methods. This renders PACA as a state-of-the-art tool for defining de novo subtypes that are more likely to reflect molecular heterogeneity, especially in challenging cases where the phenotypic heterogeneity may be masked by a myriad of strong unrelated effects in the data.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

Do gene-environment interactions have implications for the precision prevention of type 2 diabetes?

The past decades have seen a rapid global rise in the incidence of type 2 diabetes. This surge has been driven by diabetogenic environmental changes that may act together with a genetic predisposition to type 2 diabetes. It is possible that there is a synergistic gene-environment interaction, where the effects of the diabetogenic environment depend on the genetic predisposition to type 2 diabetes. Randomised trials have shown that it is possible to delay, or even prevent the development of type 2 diabetes in individuals at elevated risk through behavioural modification, focusing on weight loss, physical activity and diet. There is wide heterogeneity between individuals regarding the effectiveness of these interventions, which could, in part, be due to genetic differences. However, the studies of gene-environment interactions performed thus far suggest that behavioural modifications appear equally effective in reducing the incidence of type 2 diabetes from the stage of impaired glucose tolerance, regardless of the known underlying genetic predisposition. Recent studies suggest that there may be several subtypes of type 2 diabetes, which give new opportunities for gaining insight into gene-environment interactions. At present, the role of gene-environment interactions in the development of type 2 diabetes remains unclear. With many puzzle pieces missing in the general picture of type 2 diabetes development, the available evidence of gene-environment interactions is not ready for translation to individualised type 2 diabetes prevention based on genetic profiling.

Genetic and Environmental Contribution to the Co-Occurrence of Endocrine-Metabolic Disorders and Depression: A Nationwide Swedish Study of Siblings

OBJECTIVE

Depression is common in individuals with endocrine-metabolic disorders and vice versa, and a better understanding of the underlying factors contributing to the comorbidity of these disorders is needed. This study investigated the familial coaggregation of depression and endocrine-metabolic disorders and estimated the contribution of genetic and environmental factors to their co-occurrence.

METHODS

This population-based cohort study included 2.2 million individuals born in Sweden between 1973 and 1996, with follow-up through 2013. Participants were linked to their biological parents, allowing identification of full siblings, maternal half siblings, and paternal half siblings. Diagnoses of depression and endocrine-metabolic conditions were investigated, with the latter grouped into autoimmune disorders (autoimmune hypothyroidism, Graves' disease, and type 1 diabetes) and non-autoimmune disorders (type 2 diabetes, obesity, and polycystic ovary syndrome). Logistic regression and Cox regression were used to estimate the associations between endocrine-metabolic disorders and depression within the same individual and across siblings. Quantitative genetic modeling was performed to investigate the relative contribution of genetic and environmental influences.

RESULTS

Individuals with endocrine-metabolic disorders had a significantly higher risk of depression, with odds ratios ranging from 1.43 (95% CI=1.30, 1.57) for Graves' disease to 3.48 (95% CI=3.25, 3.72) for type 2 diabetes. Increased risks extended to full and half siblings. These correlations were mainly explained by shared genetic influences for non-autoimmune conditions, and by nonshared environmental factors for autoimmune disorders, especially for type 1 diabetes.

CONCLUSIONS

These findings provide phenotypic and etiological insights into the co-occurrence of depression and various endocrine-metabolic conditions, which could guide future research aiming at identifying pathophysiological mechanisms and intervention targets.

Mapping the genetic architecture of cortical morphology through neuroimaging: progress and perspectives

Cortical morphology is a key determinant of cognitive ability and mental health. Its development is a highly intricate process spanning decades, involving the coordinated, localized expression of thousands of genes. We are now beginning to unravel the genetic architecture of cortical morphology, thanks to the recent availability of large-scale neuroimaging and genomic data and the development of powerful biostatistical tools. Here, we review the progress made in this field, providing an overview of the lessons learned from genetic studies of cortical volume, thickness, surface area, and folding as captured by neuroimaging. It is now clear that morphology is shaped by thousands of genetic variants, with effects that are region- and time-dependent, thereby challenging conventional study approaches. The most recent genome-wide association studies have started discovering common genetic variants influencing cortical thickness and surface area, yet together these explain only a fraction of the high heritability of these measures. Further, the impact of rare variants and non-additive effects remains elusive. There are indications that the quickly increasing availability of data from whole-genome sequencing and large, deeply phenotyped population cohorts across the lifespan will enable us to uncover much of the missing heritability in the upcoming years. Novel approaches leveraging shared information across measures will accelerate this process by providing substantial increases in statistical power, together with more accurate mapping of genetic relationships. Important challenges remain, including better representation of understudied demographic groups, integration of other 'omics data, and mapping of effects from gene to brain to behavior across the lifespan.

Identification of five genetic variants with differential effects on obesity-related traits based on age

Obesity is a major public health concern, and its prevalence generally increases with age. As the number of elderly people is increasing in the aging population, the age-dependent increase in obesity has raised interest in the underlying mechanism. To understand the genetic basis of age-related increase in obesity, we identified genetic variants showing age-dependent differential effects on obesity. We conducted stratified analyses between young and old groups using genome-wide association studies of 355,335 United Kingom Biobank participants for five obesity-related phenotypes, including body mass index, body fat percentage, waist-hip ratio, waist circumference, and hip circumference. Using t-statistic, we identified five significant lead single nucleotide polymorphisms: rs2258461 with body mass index, rs9861311 and rs429358 with body fat percentage, rs2870099 with waist-hip ratio, and rs145500243 with waist circumference. Among these single nucleotide polymorphisms, rs429358, located in APOE gene was associated with diverse age-related diseases, such as Alzheimer’s disease, coronary artery disease, age-related degenerative macular diseases, and cognitive decline. The C allele of rs429358 gradually decreases body fat percentage as one grows older in the range of 40–69 years. In conclusion, we identified five genetic variants with differential effects on obesity-related phenotypes based on age using a stratified analysis between young and old groups, which may help to elucidate the mechanisms by which age influences the development of obesity.

Importance of Including Non-European Populations in Large Human Genetic Studies to Enhance Precision Medicine

One goal of genomic medicine is to uncover an individual's genetic risk for disease, which generally requires data connecting genotype to phenotype, as done in genome-wide association studies (GWAS). While there may be clinical promise to employing prediction tools such as polygenic risk scores (PRS), it currently stands that individuals of non-European ancestry may not reap the benefits of genomic medicine because of underrepresentation in large-scale genetics studies. Here, we discuss why this inequity poses a problem for genomic medicine and the reasons for the low transferability of PRS across populations. We also survey the ancestry representation of published GWAS and investigate how estimates of ancestry diversity in GWASparticipants might be biased. We highlight the importance of expanding genetic research in Africa, one of the most underrepresented regions in human genomics research, and discuss issues of ethics, resources, and technology for equitable advancement of genomic medicine.

Immuno-Modulatory Effects of Intervertebral Disc Cells

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

Age and diet shape the genetic architecture of body weight in diversity outbred mice

Understanding how genetic variation shapes a complex trait relies on accurately quantifying both the additive genetic and genotype–environment interaction effects in an age-dependent manner. We used a linear mixed model to quantify diet-dependent genetic contributions to body weight measured through adulthood in diversity outbred female mice under five diets. We observed that heritability of body weight declined with age under all diets, except the 40% calorie restriction diet. We identified 14 loci with age-dependent associations and 19 loci with age- and diet-dependent associations, with many diet-dependent loci previously linked to neurological function and behavior in mice or humans. We found their allelic effects to be dynamic with respect to genomic background, age, and diet, identifying several loci where distinct alleles affect body weight at different ages. These results enable us to more fully understand and predict the effectiveness of dietary intervention on overall health throughout age in distinct genetic backgrounds.

Body weight is one trait influenced by genes, age and environmental factors. Both internal and external environmental pressures are known to affect genetic variation over time. However, it is largely unknown how all factors – including age – interact to shape metabolism and bodyweight.

Wright et al. set out to quantify the interactions between genes and diet in ageing mice and found that the effect of genetics on mouse body weight changes with age. In the experiments, Wright et al. weighed 960 female mice with diverse genetic backgrounds, starting at two months of age into adulthood. The animals were randomized to different diets at six months of age. Some mice had unlimited food access, others received 20% or 40% less calories than a typical mouse diet, and some fasted one or two days per week.

Variations in their genetic background explained about 80% of differences in mice’s weight, but the influence of genetics relative to non-genetic factors decreased as they aged. Mice on the 40% calorie restriction diet were an exception to this rule and genetics accounted for 80% of their weight throughout adulthood, likely due to reduced influence from diet and reduced interactions between diet and genes. Several genes involved in metabolism, neurological function, or behavior, were associated with mouse weight.

The experiments highlight the importance of considering interactions between genetics, environment, and age in determining complex traits like body weight. The results and the approaches used by Wright et al. may help other scientists learn more about how the genetic predisposition to disease changes with environmental stimuli and age.

Heritability Analyses Uncover Shared Genetic Effects of Lung Function and Change over Time

Genetic influence on lung functions has been identified in previous studies; however, the relative longitudinal effects of genetic factors and their interactions with smoking on lung function remain unclear. Here, we identified the longitudinal effects of genetic variants on lung function by determining single nucleotide polymorphism (SNP) heritability and genetic correlations, and by analyzing interactions with smoking. Subject-specific means and annual change rates were calculated for eight spirometric measures obtained from 6622 Korean adults aged 40−69 years every two years for 14 years, and their heritabilities were estimated separately. Statistically significant (p < 0.05) heritability for the subject-specific means of all spirometric measures (8~32%) and change rates of forced expiratory volume in 1 s to forced vital capacity ratio (FEV1/FVC; 16%) and post-bronchodilator FEV1/FVC (17%) were detected. Significant genetic correlations of the change rate with the subject-specific mean were observed for FEV1/FVC (ρg = 0.64) and post-bronchodilator FEV1/FVC (ρg = 0.47). Furthermore, post-bronchodilator FEV1/FVC showed significant heritability of SNP-by-smoking interaction (hGXS2 = 0.4) for the annual change rate. The GWAS also detected genome-wide significant SNPs for FEV1 (rs4793538), FEV1/FVC (rs2704589, rs62201158, and rs9391733), and post-bronchodilator FEV1/FVC (rs2445936). We found statistically significant evidence of heritability role on the change in lung function, and this was shared with the effects on cross-sectional measurements. We also found some evidence of interaction with smoking for the change of lung function.

An opportunity for primary prevention research in psychotic disorders

An opportunity has opened for research into primary prevention of psychotic disorders, based on progress in endophenotypes, genetics, and genomics. Primary prevention requires reliable prediction of susceptibility before any symptoms are present. We studied a battery of measures where published data supports abnormalities of these measurements prior to appearance of initial psychosis symptoms. These neurobiological and behavioral measurements included cognition, eye movement tracking, Event Related Potentials, and polygenic risk scores. They generated an acceptably precise separation of healthy controls from outpatients with a psychotic disorder. METHODS: The Bipolar and Schizophrenia Network on Intermediate Phenotypes (B-SNIP) measured this battery in an ancestry-diverse series of consecutively recruited adult outpatients with a psychotic disorder and healthy controls. Participants include all genders, 16 to 50 years of age, 261 with psychotic disorders (Schizophrenia (SZ) 109, Bipolar with psychosis (BPP) 92, Schizoaffective disorder (SAD) 60), 110 healthy controls. Logistic Regression, and an extension of the Linear Mixed Model to include analysis of pairwise interactions between measures (Environmental kernel Relationship Matrices (ERM)) with multiple iterations, were performed to predict case-control status. Each regression analysis was validated with four-fold cross-validation. RESULTS AND CONCLUSIONS: Sensitivity, specificity, and Area Under the Curve of Receiver Operating Characteristic of 85%, 62%, and 86%, respectively, were obtained for both analytic methods. These prediction metrics demonstrate a promising diagnostic distinction based on premorbid risk variables. There were also statistically significant pairwise interactions between measures in the ERM model. The strong prediction metrics of both types of analytic model provide proof-of-principle for biologically-based laboratory tests as a first step toward primary prevention studies. Prospective studies of adolescents at elevated risk, vs. healthy adolescent controls, would be a next step toward development of primary prevention strategies.

A compressed variance component mixed model for detecting QTNs and QTN-by-environment and QTN-by-QTN interactions in genome-wide association studies

Although genome-wide association studies are widely used to mine genes for quantitative traits, the effects to be estimated are confounded, and the methodologies for detecting interactions are imperfect. To address these issues, the mixed model proposed here first estimates the genotypic effects for AA, Aa, and aa, and the genotypic polygenic background replaces additive and dominance polygenic backgrounds. Then, the estimated genotypic effects are partitioned into additive and dominance effects using a one-way analysis of variance model. This strategy was further expanded to cover QTN-by-environment interactions (QEIs) and QTN-by-QTN interactions (QQIs) using the same mixed-model framework. Thus, a three-variance-component mixed model was integrated with our multi-locus random-SNP-effect mixed linear model (mrMLM) method to establish a new methodological framework, 3VmrMLM, that detects all types of loci and estimates their effects. In Monte Carlo studies, 3VmrMLM correctly detected all types of loci and almost unbiasedly estimated their effects, with high powers and accuracies and a low false positive rate. In re-analyses of 10 traits in 1439 rice hybrids, detection of 269 known genes, 45 known gene-by-environment interactions, and 20 known gene-by-gene interactions strongly validated 3VmrMLM. Further analyses of known genes showed more small (67.49%), minor-allele-frequency (35.52%), and pleiotropic (30.54%) genes, with higher repeatability across datasets (54.36%) and more dominance loci. In addition, a heteroscedasticity mixed model in multiple environments and dimension reduction methods in quite a number of environments were developed to detect QEIs, and variable selection under a polygenic background was proposed for QQI detection. This study provides a new approach for revealing the genetic architecture of quantitative traits.

Discrimination of Geographical Origins of Wolfberry (Lycium barbarum L.) Fruits Using Stable Isotopes, Earth Elements, Free Amino Acids, and Saccharides

To develop sophisticated approaches for distinguishing goji origins, 325 wolfberry fruit samples of a certain cultivar, plant age, drying method, and collection season were gathered from 26 producing areas across Northwest China in 2017 and 2018. We employed 49 indices, including stable isotopes, earth elements, soluble amino acids, and saccharides, to identify the regions of origin of these goji fruits. Analysis of variance (ANOVA) and heritability analysis were used to assess the effects of the environment (producing areas), cultivar, plant age, drying process, and collection season. Samples from the same place can be classified and partially discriminated using principal component analysis (PCA). We were able to distinguish fruits produced in Zhongning County from those produced in the other five producing provinces using orthogonal projection to latent structure-discriminant analysis (OPLS-DA). Calcium (Ca), manganese (Mn), ornithine (Orn), cystine (Cys-Cys), glutamate (Glu), phenylalanine (Phe), phosphoserine (Ps), serine (Ser), lysine (Lys), taurine (Tau), proline (Pro), and tyrosine (Tyr) indices were chosen using S-plots and heritability analysis, and their repeatability was established with samples collected in 2018. The indices selected in this study can distinguish goji berries produced in Zhongning County from fruits originating from five other Provinces with high repeatability, which was validated with various cultivars, drying methods, harvest seasons, and plant ages and with heritability analysis.

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.

Body mass index interacts with a genetic-risk score for depression increasing the risk of the disease in high-susceptibility individuals

Depression is strongly associated with obesity among other chronic physical diseases. The latest mega- and meta-analysis of genome-wide association studies have identified multiple risk loci robustly associated with depression. In this study, we aimed to investigate whether a genetic-risk score (GRS) combining multiple depression risk single nucleotide polymorphisms (SNPs) might have utility in the prediction of this disorder in individuals with obesity. A total of 30 depression-associated SNPs were included in a GRS to predict the risk of depression in a large case-control sample from the Spanish PredictD-CCRT study, a national multicentre, randomized controlled trial, which included 104 cases of depression and 1546 controls. An unweighted GRS was calculated as a summation of the number of risk alleles for depression and incorporated into several logistic regression models with depression status as the main outcome. Constructed models were trained and evaluated in the whole recruited sample. Non-genetic-risk factors were combined with the GRS in several ways across the five predictive models in order to improve predictive ability. An enrichment functional analysis was finally conducted with the aim of providing a general understanding of the biological pathways mapped by analyzed SNPs. We found that an unweighted GRS based on 30 risk loci was significantly associated with a higher risk of depression. Although the GRS itself explained a small amount of variance of depression, we found a significant improvement in the prediction of depression after including some non-genetic-risk factors into the models. The highest predictive ability for depression was achieved when the model included an interaction term between the GRS and the body mass index (BMI), apart from the inclusion of classical demographic information as marginal terms (AUC = 0.71, 95% CI = [0.65, 0.76]). Functional analyses on the 30 SNPs composing the GRS revealed an over-representation of the mapped genes in signaling pathways involved in processes such as extracellular remodeling, proinflammatory regulatory mechanisms, and circadian rhythm alterations. Although the GRS on its own explained a small amount of variance of depression, a significant novel feature of this study is that including non-genetic-risk factors such as BMI together with a GRS came close to the conventional threshold for clinical utility used in ROC analysis and improves the prediction of depression. In this study, the highest predictive ability was achieved by the model combining the GRS and the BMI under an interaction term. Particularly, BMI was identified as a trigger-like risk factor for depression acting in a concerted way with the GRS component. This is an interesting finding since it suggests the existence of a risk overlap between both diseases, and the need for individual depression genetics-risk evaluation in subjects with obesity. This research has therefore potential clinical implications and set the basis for future research directions in exploring the link between depression and obesity-associated disorders. While it is likely that future genome-wide studies with large samples will detect novel genetic variants associated with depression, it seems clear that a combination of genetics and non-genetic information (such is the case of obesity status and other depression comorbidities) will still be needed for the optimization prediction of depression in high-susceptibility individuals.

Lifestyle Modifies the Diabetes-Related Metabolic Risk, Conditional on Individual Genetic Differences

Metabolic syndrome is a group of heritable metabolic traits that are highly associated with type 2 diabetes (T2DM). Classical interventions to T2DM include individual self-management of environmental risk factors, such as improving diet quality, increasing physical activity, and reducing smoking and alcohol consumption, which decreases the risk of developing metabolic syndrome. However, it is poorly understood how the phenotypes of diabetes-related metabolic traits change with respect to lifestyle modifications at the individual level. In the analysis, we used 12 diabetes-related metabolic traits and eight lifestyle covariates from the UK Biobank comprising 288,837 white British participants genotyped for 1,133,273 genome-wide single nucleotide polymorphisms. We found 16 GxE interactions. Modulation of genetic effects by physical activity was seen for four traits (glucose, HbA1c, C-reactive protein, systolic blood pressure) and by alcohol and smoking for three (BMI, glucose, waist-hip ratio and BMI and diastolic and systolic blood pressure, respectively). We also found a number of significant phenotypic modulations by the lifestyle covariates, which were not attributed to the genetic effects in the model. Overall, modulation in the metabolic risk in response to the level of lifestyle covariates was clearly observed, and its direction and magnitude were varied depending on individual differences. We also showed that the metabolic risk inferred by our model was notably higher in T2DM prospective cases than controls. Our findings highlight the importance of individual genetic differences in the prevention and management of diabetes and suggest that the one-size-fits-all approach may not benefit all.

A genome-wide association study of childhood adiposity and blood lipids

Background: The rising prevalence of childhood obesity and dyslipidaemia is a major public health concern due to its association with morbidity and mortality in later life. Methods: In this study, we have conducted genome-wide association studies (GWAS) for eight measures of adiposity and lipids in a cohort of young individuals (mean age 9.9) from the Avon Longitudinal Study of Parents and Children (ALSPAC). These measures were body mass index (BMI), systolic and diastolic blood pressure, high- density and low-density lipoprotein cholesterol, triglycerides, apolipoprotein A-I and apolipoprotein B. We next undertook functional enrichment, pathway analyses and linkage disequilibrium (LD) score regression to evaluate genetic correlations with later-life cardiometabolic diseases. Results: Using GWAS we identified 14 unique loci associated with at least one risk factor in this cohort of age 10 individuals (P<5x10-8), with lipoprotein lipid-associated loci being enriched for liver tissue-derived gene expression and lipid synthesis pathways. LD score regression provided evidence of various genetic correlations, such as childhood systolic blood pressure being genetically correlated with later-life coronary artery disease (rG=0.26, 95% CI=0.07 to 0.46, P=0.009) and hypertension (rG=0.37, 95% CI=0.19 to 0.55, P=6.57x10-5), as well as childhood BMI with type 2 diabetes (rG=0.35, 95% CI=0.18 to 0.51, P=3.28x10-5). Conclusions: Our findings suggest that there are genetic variants inherited at birth which begin to exert their effects on cardiometabolic risk factors as early as age 10 in the life course. However, further research is required to assess whether the genetic correlations we have identified are due to direct or indirect effects of childhood adiposity and lipid traits.

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.

Non-linear interaction between physical activity and polygenic risk score of body mass index in Danish and Russian populations

Body mass index (BMI) is a highly heritable polygenic trait. It is also affected by various environmental and behavioral risk factors. We used a BMI polygenic risk score (PRS) to study the interplay between the genetic and environmental factors defining BMI. First, we generated a BMI PRS that explained more variance than a BMI genetic risk score (GRS), which was using only genome-wide significant BMI-associated variants (R² = 13.1% compared to 6.1%). Second, we analyzed interactions between BMI PRS and seven environmental factors. We found a significant interaction between physical activity and BMI PRS, even when the well-known effect of the FTO region was excluded from the PRS, using a small dataset of 6,179 samples. Third, we stratified the study population into two risk groups using BMI PRS. The top 22% of the studied populations were included in a high PRS risk group. Engagement in self-reported physical activity was associated with a 1.66 kg/m² decrease in BMI in this group, compared to a 0.84 kg/m² decrease in BMI in the rest of the population. Our results (i) confirm that genetic background strongly affects adult BMI in the general population, (ii) show a non-linear interaction between BMI genetics and physical activity, and (iii) provide a standardized framework for future gene-environment interaction analyses.

Little Evidence of Modified Genetic Effect of rs16969968 on Heavy Smoking Based on Age of Onset of Smoking

INTRODUCTION

Tobacco smoking is the leading cause of preventable death globally. Smoking quantity, measured in cigarettes per day, is influenced both by the age of onset of regular smoking (AOS) and by genetic factors, including a strong effect of the nonsynonymous single-nucleotide polymorphism rs16969968. A previous study by Hartz et al. reported an interaction between these two factors, whereby rs16969968 risk allele carriers who started smoking earlier showed increased risk for heavy smoking compared with those who started later. This finding has yet to be replicated in a large, independent sample.

METHODS

We performed a preregistered, direct replication attempt of the rs16969968 × AOS interaction on smoking quantity in 128 383 unrelated individuals from the UK Biobank, meta-analyzed across ancestry groups. We fit statistical association models mirroring the original publication as well as formal interaction tests on multiple phenotypic and analytical scales.

RESULTS

We replicated the main effects of rs16969968 and AOS on cigarettes per day but failed to replicate the interaction using previous methods. Nominal significance of the rs16969968 × AOS interaction term depended strongly on the scale of analysis and the particular phenotype, as did associations stratified by early/late AOS. No interaction tests passed genome-wide correction (α = 5e-8), and all estimated interaction effect sizes were much smaller in magnitude than previous estimates.

CONCLUSIONS

We failed to replicate the strong rs16969968 × AOS interaction effect previously reported. If such gene-moderator interactions influence complex traits, they likely depend on scale of measurement, and current biobanks lack the power to detect significant genome-wide associations given the minute effect sizes expected.

IMPLICATIONS

We failed to replicate the strong rs16969968 × AOS interaction effect on smoking quantity previously reported. If such gene-moderator interactions influence complex traits, current biobanks lack the power to detect significant genome-wide associations given the minute effect sizes expected. Furthermore, many potential interaction effects are likely to depend on the scale of measurement employed.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

FTO Obesity Variant-Exercise Interaction on Changes in Body Weight and BMI: The Taiwan Biobank Study

CONTEXT

Gene-exercise interaction on cross-sectional body mass index (BMI) has been extensively studied and is well established. However, gene-exercise interaction on changes in body weight/BMI remains controversial.

OBJECTIVE

To examine the interaction between the FTO obesity variant and regular exercise on changes in body weight/BMI.

PARTICIPANTS

Taiwan Biobank participants aged 30-70 years (N = 20 906) were examined at both baseline and follow-up visit (mean follow-up duration: 3.7 years).

MAIN OUTCOME MEASURES

The interaction between the FTO obesity variant rs1421085 and regular exercise habit (no exercise, ≤20 metabolic equivalent of tasks (METs)/week exercise, >20 METs/week exercise) on changes in body weight/BMI.

RESULTS

Individuals with the risk allele of rs1421085 gained more weight and increased BMI than those without the risk allele if they did not exercise. In contrast, individuals with the risk allele gained less weight and BMI if they exercised regularly, indicating an interaction between rs1421085 and regular exercise habit (P = .030 for Δbody weight and P = .034 for ΔBMI). The effect of exercise on maintaining body weight was larger in those with the risk allele of rs1421085. When we focused on individuals without regular exercise at baseline, individuals with the risk allele again tended to lose more weight than those with a nonrisk allele if they had acquired an exercise habit by the follow-up visit.

CONCLUSION

The beneficial effect of exercise is greater in individuals genetically prone to obesity due to the interaction between the FTO obesity variant rs1421085 and regular exercise on changes in body weight and BMI.

Genetic Influences on Disease Subtypes

Disease classification, or nosology, was historically driven by careful examination of clinical features of patients. As technologies to measure and understand human phenotypes advanced, so too did classifications of disease, and the advent of genetic data has led to a surge in genetic subtyping in the past decades. Although the fundamental process of refining disease definitions and subtypes is shared across diverse fields, each field is driven by its own goals and technological expertise, leading to inconsistent and conflicting definitions of disease subtypes. Here, we review several classical and recent subtypes and subtyping approaches and provide concrete definitions to delineate subtypes. In particular, we focus on subtypes with distinct causal disease biology, which are of primary interest to scientists, and subtypes with pragmatic medical benefits, which are of primary interest to physicians. We propose genetic heterogeneity as a gold standard for establishing biologically distinct subtypes of complex polygenic disease. We focus especially on methods to find and validate genetic subtypes, emphasizing common pitfalls and how to avoid them.

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

Genetic variation in response to the environment, that is, genotype-by-environment interaction (GxE), is fundamental in the biology of complex traits and diseases. However, existing methods are computationally demanding and infeasible to handle biobank-scale data. Here, we introduce GxEsum, a method for estimating the phenotypic variance explained by genome-wide GxE based on GWAS summary statistics. Through comprehensive simulations and analysis of UK Biobank with 288,837 individuals, we show that GxEsum can handle a large-scale biobank dataset with controlled type I error rates and unbiased GxE estimates, and its computational efficiency can be hundreds of times higher than existing GxE methods.

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes

BACKGROUND/OBJECTIVES

Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil.

SUBJECTS/METHODS

Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European, and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of body mass index (BMI) in three Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambuí), and South (Pelotas).

RESULTS

We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p = 2.76e-06). This variant is rare in Europeans but with frequencies of ~3% in West Africa and has a strong female-specific effect (95% CI: 2.32-5.65 kg/m² per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among women with morbid obesity from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects.

CONCLUSIONS

We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.

Strategies to Understand the Weight-Reduced State: Genetics and Brain Imaging

Most individuals with obesity or overweight have difficulty maintaining weight loss. The weight-reduced state induces changes in many physiological processes that appear to drive weight regain. Here, we review the use of cell biology, genetics, and imaging techniques that are being used to begin understanding why weight regain is the normal response to dieting. As with obesity itself, weight regain has both genetic and environmental drivers. Genetic drivers for "thinness" and "obesity" largely overlap, but there is evidence for specific genetic loci that are different for each of these weight states. There is only limited information regarding the genetics of weight regain. Currently, most genetic loci related to weight point to the central nervous system as the organ responsible for determining the weight set point. Neuroimaging tools have proved useful in studying the contribution of the central nervous system to the weight-reduced state in humans. Neuroimaging technologies fall into three broad categories: functional, connectivity, and structural neuroimaging. Connectivity and structural imaging techniques offer unique opportunities for testing mechanistic hypotheses about changes in brain function or tissue structure in the weight-reduced state.

Evidence for shared genetics between physical activity, sedentary behaviour and adiposity-related traits

Observational, cross-sectional and longitudinal studies showed that physical activity and sedentary behaviour are associated with adiposity-related traits, apparently in a bidirectional manner. Physical activity is also suggested to suppress the genetic risk of adiposity. Since phenotypic associations with genetic variants are not subject to reverse causation or confounding, they may be used as tools to shed light on cause and effect in this complex interdependency. We review the evidence for shared genetics of physical activity and adiposity-related traits and for gene-by-physical activity interactions on adiposity-related traits in human studies. We outline limitations, challenges and opportunities in studying and understanding of these relationships. In summary, physical activity and sedentary behaviour are genetically correlated with body mass index and fat percentage but may not be correlated with lean body mass. Mendelian randomisation analyses show that physical activity and sedentary behaviour have bidirectional relationships with adiposity. Several studies suggest that physical activity suppresses genetic risk of adiposity. No studies have yet tested whether adiposity enhances genetic predisposition to sedentariness. The complexity of the comprehensive causal model makes the assessment of the single or combined components challenging. Substantial progress in this field may need long-term intervention studies.

Novelty seeking is associated with increased body weight and orbitofrontal grey matter volume reduction

Novelty seeking (NS) has previously been identified as a personality trait that is associated with elevated body mass index (BMI) and obesity. Of note, both obesity and reduced impulse control - a core feature of NS - have previously been associated with grey matter volume (GMV) reductions in the orbitofrontal cortex (OFC). Yet, it remains unknown, if body weight-related grey matter decline in the OFC might be explained by higher levels of NS. To address this question, we studied associations between NS, BMI and brain structure in 355 healthy subjects. Brain images were pre-processed using voxel-based morphometry (VBM). BMI was calculated from self-reported height and weight. The Tridimensional Personality Questionnaire (TPQ) was used to assess NS. NS and BMI were associated positively (r = .137, p = .01) with NS being a significant predictor of BMI (B = 0.172; SE B = 0.05; ß = 0.184; p = 0.001). Significant associations between BMI and GMV specifically in the OFC (x = -44, y = 56, z = -2, t(350) = 4.34, k = 5, p_FWE = 0.011) did not uphold when correcting for NS in the model. In turn, a significant negative association between NS and OFC GMV was found independent of BMI (x = -2, y = 48, z = -10, t(349) = 4.42, k = 88, p_FWE = 0.008). Body mass-related grey matter decrease outside the OFC could not be attributed to NS. Our results suggest that body-weight-related orbitofrontal grey matter reduction can at least partly be linked to higher levels of NS. Given the pivotal role of the OFC in overweight as well as cognitive domains such as impulse inhibition, executive control and reward processing, its association with NS seems to provide a tenable neurobiological correlate for future research.

Quantifying genetic heterogeneity between continental populations for human height and body mass index

Genome-wide association studies (GWAS) in samples of European ancestry have identified thousands of genetic variants associated with complex traits in humans. However, it remains largely unclear whether these associations can be used in non-European populations. Here, we seek to quantify the proportion of genetic variation for a complex trait shared between continental populations. We estimated the between-population correlation of genetic effects at all SNPs (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{g}$$\end{document}rg) or genome-wide significant SNPs (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{{g\left( {GWS} \right)}}$$\end{document}rgGWS) for height and body mass index (BMI) in samples of European (EUR; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 49,839$$\end{document}n=49,839) and African (AFR; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 17,426$$\end{document}n=17,426) ancestry. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{g}$$\end{document}r^g between EUR and AFR was 0.75 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.035$$\end{document}s.e.=0.035) for height and 0.68 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.062$$\end{document}s.e.=0.062) for BMI, and the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{{g\left( {GWS} \right)}}$$\end{document}r^gGWS was 0.82 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.030$$\end{document}s.e.=0.030) for height and 0.87 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.064$$\end{document}s.e.=0.064) for BMI, suggesting that a large proportion of GWAS findings discovered in Europeans are likely applicable to non-Europeans for height and BMI. There was no evidence that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{g}$$\end{document}r^g differs in SNP groups with different levels of between-population difference in allele frequency or linkage disequilibrium, which, however, can be due to the lack of power.