Genotype-covariate correlation and interaction disentangled by a whole-genome multivariate reaction norm model

The integration of quantile regression with 3VmrMLM identifies more QTNs and QTN-by-environment interactions using SNP- and haplotype-based markers

Current methods used in genome-wide association studies frequently lack power owing to their inability to detect heterogeneous associations and rare and multiallelic variants. To address these issues, quantile regression is integrated with a three (compressed) variance component multi-locus random-SNP-effect mixed linear model (3VmrMLM) to propose q3VmrMLM for detecting heterogeneous quantitative trait nucleotides (QTNs) and QTN-by-environment interactions (QEIs), and then design haplotype-based q3VmrMLM (q3VmrMLM-Hap) for identifying multiallelic haplotypes and rare variants. In Monte Carlo simulation studies, q3VmrMLM had higher power than 3VmrMLM, sequence kernel association test (SKAT), and integrated quantile rank test (iQRAT). In a re-analysis of 10 traits in 1439 rice hybrids, 261 known genes were identified only by q3VmrMLM and q3VmrMLM-Hap, whereas 175 known genes were detected by both the new and existing methods. Of all the significant QTNs with known genes, q3VmrMLM (179: 140 variance heterogeneity and 157 quantile effect heterogeneity) found more heterogeneous QTNs than 3VmrMLM (123), SKAT (27), and iQRAT (29); q3VmrMLM-Hap (121) mapped more low-frequency (<0.05) QTNs than q3VmrMLM (51), 3VmrMLM (43), SKAT (11), and iQRAT (12); and q3VmrMLM-Hap (12), q3VmrMLM (16), and 3VmrMLM (12) had similar power in identifying gene-by-environment interactions. All significant and suggested QTNs achieved the highest predictive accuracy (r = 0.9045). In conclusion, this study describes a new and complementary approach to mining genes and unraveling the genetic architecture of complex traits in crops.

Evolution, genetic diversity, and health

Human genetic diversity in today's world has been shaped by evolutionary history, demographic shifts and environmental exposures, influencing complex traits, disease susceptibility and drug responses. Capturing this diversity is essential for advancing precision medicine and promoting equitable healthcare. Despite the great progress achieved with initiatives such as the human Pangenome and large biobanks that aim for a better representation of human diversity, important challenges remain. In this Perspective, we discuss the importance of diversity in clinical genomics through an evolutionary lens. We highlight progress and challenges and outline key clinical applications of diverse genetic data. We argue that diversifying both datasets and methodologies-integrating ancestral and environmental factors-is crucial for fully understanding the genetic basis of human health and disease.

Detecting Non-linear Dependence through Genome Wide Analysis

In the current study we introduce statistical methods based on trigonometry, to infer the shape of a (non)-linear bivariate genetic relationship. We do this based on a series of piecemeal GWASs of segments of a target (continuous) trait distribution, and the genetic correlations between those GWASs and a second trait. Simulations confirm that we are able to retrieve the shape of the relationship given certain assumptions about the nature of the relationship between the traits. We applied the method to the genetic relationship between BMI, sleep duration, and height, and psychiatric disorders (ADHD, anorexia nervosa, and depression) using data from approximately 450K individuals from UK Biobank. In the relationship between BMI and psychiatric traits, we found that the expected value of depression is a nonlinear function of BMI i.e. there is a nonlinear genetic relationship between both traits. We observed similar findings for the genetic relationship between BMI and anorexia, sleep duration and depression, and sleep duration and ADHD. We observed no underlying nonlinearity in the genetic relationship between height and psychiatric traits. Using a novel statistical approach, we show that nonlinear genetic relationships between traits are detectable and genetic associations as quantified using global estimators like genetic correlations are not informative about underlying complexities in these relationships. Our findings challenge assumptions of linearity in genetic epidemiology and suggest that bivariate genetic associations are not uniform across the phenotypic spectrum, which may have implications for the development of targeted interventions.

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.

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.

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.

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.

Association between hepatitis A seropositivity and bone mineral density in adolescents and adults: a cross-sectional study using NHANES data

BACKGROUND

Osteoporosis, characterized by decreased bone density and increased fracture risk, imposes significant physical, psychosocial, and financial burdens. Early detection and prevention are crucial for managing osteoporosis and reducing the risk of fractures.

OBJECTIVES

To investigate the relationship between Hepatitis A seropositivity and bone mineral density (BMD) in adolescents and adults and to explore the potential link between Hepatitis A infection and osteoporosis risk.

DESIGN AND SETTING

This cross-sectional study used data from the National Health and Nutrition Examination Survey (NHANES) from 2011 to 2018 to evaluate the association between hepatitis A seropositivity and BMD in 15,693 participants.

METHODS

Multivariable regression analysis was used to calculate the mean BMD and standard error for adolescents and adults, followed by an independent z-test to determine whether there was a significant difference between the seropositive and seronegative groups.

RESULTS

Hepatitis A seropositive adolescents and adults had lower BMD than their seronegative counterparts, with significant differences in lumber spine (mean difference = -0.03 g/cm2, P < 0.01 for both age groups) and pelvis BMDs (mean difference = -0.02 g/cm2, P < 0.01 for the adult age groups), after adjusting for various covariates.

CONCLUSIONS

This study confirmed that both adolescent and adult individuals seropositive for Hepatitis A antibodies had reduced BMD among both adolescents and adults, especially in the adult group. This finding suggests a possible link between Hepatitis A infection and risk of osteoporosis.

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

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

Unraveling phenotypic variance in metabolic syndrome through multi-omics

Complex multi-omics effects drive the clustering of cardiometabolic risk factors, underscoring the imperative to comprehend how individual and combined omics shape phenotypic variation. Our study partitions phenotypic variance in metabolic syndrome (MetS), blood glucose (GLU), triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and blood pressure through genome, transcriptome, metabolome, and exposome (i.e., lifestyle exposome) analyses. Our analysis included a cohort of 62,822 unrelated individuals with white British ancestry, sourced from the UK biobank. We employed linear mixed models to partition phenotypic variance using the restricted maximum likelihood (REML) method, implemented in MTG2 (v2.22). We initiated the analysis by individually modeling omics, followed by subsequent integration of pairwise omics in a joint model that also accounted for the covariance and interaction between omics layers. Finally, we estimated the correlations of various omics effects between the phenotypes using bivariate REML. Significant proportions of the MetS variance were attributed to distinct data sources: genome (9.47%), transcriptome (4.24%), metabolome (14.34%), and exposome (3.77%). The phenotypic variances explained by the genome, transcriptome, metabolome, and exposome ranged from 3.28% for GLU to 25.35% for HDL-C, 0% for GLU to 19.34% for HDL-C, 4.29% for systolic blood pressure (SBP) to 35.75% for TG, and 0.89% for GLU to 10.17% for HDL-C, respectively. Significant correlations were found between genomic and transcriptomic effects for TG and HDL-C. Furthermore, significant interaction effects between omics data were detected for both MetS and its components. Interestingly, significant correlation of omics effect between the phenotypes was found. This study underscores omics' roles, interaction effects, and random-effects covariance in unveiling phenotypic variation in multi-omics domains.

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.

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.

Partitioning the forms of genotype-by-environment interaction in the reaction norm analysis of stability

KEY MESSAGE

The reaction norm analysis of stability can be enhanced by partitioning the contribution of different types of G × E to the variation in slope. The slope of regression in a reaction norm model, where the performance of a genotype is regressed over an environmental covariable, is often used as a measure of stability of genotype performance. This method could be developed further by partitioning variation in the slope of regression into the two sources of genotype-by-environment interaction (G × E) which cause it: scale-type G × E (heterogeneity of variance) and rank-type G × E (heterogeneity of correlation). Because the two types of G × E have very different properties, separating their effect would enable a clearer understanding of stability. The aim of this paper was to demonstrate two methods which seek to achieve this in reaction norm models. Reaction norm models were fit to yield data from a multi-environment trial in Barley (Hordeum vulgare), with the adjusted mean yield from each environment used as the environmental covariable. Stability estimated from factor-analytic models, which can disentangle the two types of G × E and estimate stability based on rank-type G × E, was used for comparison. Adjusting the reaction norm slope to account for scale-type G × E using a genetic regression more than tripled the correlation with factor-analytic estimates of stability (0.24-0.26 to 0.80-0.85), indicating that it removed variation in the reaction norm slope that originated from scale-type G × E. A standardisation procedure had a more modest increase (055-0.59) but could be useful when curvilinear reaction norms are required. Analyses which use reaction norms to explore the stability of genotypes could gain additional insight into the mechanisms of stability by applying the methods outlined in this study.

A fast and powerful linear mixed model approach for genotype-environment interaction tests in large-scale GWAS

Genotype-by-environment interaction (GEI or GxE) plays an important role in understanding complex human traits. However, it is usually challenging to detect GEI signals efficiently and accurately while adjusting for population stratification and sample relatedness in large-scale genome-wide association studies (GWAS). Here we propose a fast and powerful linear mixed model-based approach, fastGWA-GE, to test for GEI effect and G + GxE joint effect. Our extensive simulations show that fastGWA-GE outperforms other existing GEI test methods by controlling genomic inflation better, providing larger power and running hundreds to thousands of times faster. We performed a fastGWA-GE analysis of ~7.27 million variants on 452 249 individuals of European ancestry for 13 quantitative traits and five environment variables in the UK Biobank GWAS data and identified 96 significant signals (72 variants across 57 loci) with GEI test P-values < 1 × 10-9, including 27 novel GEI associations, which highlights the effectiveness of fastGWA-GE in GEI signal discovery in large-scale GWAS.

Exploring polygenic-environment and residual-environment interactions for depressive symptoms within the UK Biobank

Substantial advances have been made in identifying genetic contributions to depression, but little is known about how the effect of genes can be modulated by the environment, creating a gene-environment interaction. Using multivariate reaction norm models (MRNMs) within the UK Biobank (N = 61294-91644), we investigate whether the polygenic and residual variance components of depressive symptoms are modulated by 17 a priori selected covariate traits-12 environmental variables and 5 biomarkers. MRNMs, a mixed-effects modelling approach, provide unbiased polygenic-covariate interaction estimates for a quantitative trait by controlling for outcome-covariate correlations and residual-covariate interactions. A continuous depressive symptom variable was the outcome in 17 MRNMs-one for each covariate trait. Each MRNM had a fixed-effects model (fixed effects included the covariate trait, demographic variables, and principal components) and a random effects model (where polygenic-covariate and residual-covariate interactions are modelled). Of the 17 selected covariates, 11 significantly modulate deviations in depressive symptoms through the modelled interactions, but no single interaction explains a large proportion of phenotypic variation. Results are dominated by residual-covariate interactions, suggesting that covariate traits (including neuroticism, childhood trauma, and BMI) typically interact with unmodelled variables, rather than a genome-wide polygenic component, to influence depressive symptoms. Only average sleep duration has a polygenic-covariate interaction explaining a demonstrably nonzero proportion of the variability in depressive symptoms. This effect is small, accounting for only 1.22% (95% confidence interval: [0.54, 1.89]) of variation. The presence of an interaction highlights a specific focus for intervention, but the negative results here indicate a limited contribution from polygenic-environment interactions.

From Mendel to quantitative genetics in the genome era: the scientific legacy of W. G. Hill

The quantitative geneticist W. G. ('Bill') Hill, awardee of the 2018 Darwin Medal of the Royal Society and the 2019 Mendel Medal of the Genetics Society (United Kingdom), died on 17 December 2021 at the age of 81 years. Here, we pay tribute to his multiple key scientific contributions, which span population and evolutionary genetics, animal and plant breeding and human genetics. We discuss his theoretical research on the role of linkage disequilibrium (LD) and mutational variance in the response to selection, the origin of the widely used LD metric r² in genomic association studies, the genetic architecture of complex traits, the quantification of the variation in realized relationships given a pedigree relationship and much more. We demonstrate that basic theoretical research in quantitative and statistical genetics has led to profound insights into the genetics and evolution of complex traits and made predictions that were subsequently empirically validated, often decades later.

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

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

Genomic analysis of the slope of the reaction norm for body weight in Australian sheep

Background

Selection of livestock based on their robustness or sensitivity to environmental variation could help improve the efficiency of production systems, particularly in the light of climate change. Genetic variation in robustness arises from genotype-by-environment (G × E) interactions, with genotypes performing differently when animals are raised in contrasted environments. Understanding the nature of this genetic variation is essential to implement strategies to improve robustness. In this study, our aim was to explore the genetics of robustness in Australian sheep to different growth environments using linear reaction norm models (RNM), with post-weaning weight records of 22,513 lambs and 60 k single nucleotide polymorphisms (SNPs). The use of scale-corrected genomic estimated breeding values (GEBV) for the slope to account for scale-type G × E interactions was also investigated.

Results

Additive genetic variance was observed for the slope of the RNM, with genetic correlations between low- and high-growth environments indicating substantial re-ranking of genotypes (0.44–0.49). The genetic variance increased from low- to high-growth environments. The heritability of post-weaning body weight ranged from 0.28 to 0.39. The genetic correlation between intercept and slope of the reaction norm for post-weaning body weight was low to moderate when based on the estimated (co)variance components but was much higher when based on back-solved SNP effects. An initial analysis suggested that a region on chromosome 11 affected both the intercept and the slope, but when the GEBV for the slope were conditioned on the GEBV for the intercept to remove the effect of scale-type G × E interactions on SNP effects for robustness, a single genomic region on chromosome 7 was found to be associated with robustness. This region included genes previously associated with growth traits and disease susceptibility in livestock.

Conclusions

This study shows a significant genetic variation in the slope of RNM that could be used for selecting for increased robustness of sheep. Both scale-type and rank-type G × E interactions contributed to variation in the slope. The correction for scale effects of GEBV for the slope should be considered when analysing robustness using RNM. Overall, robustness appears to be a highly polygenic trait.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00734-6.

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.

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.

Adiposity and cancer: a Mendelian randomization analysis in the UK biobank

BACKGROUND

Observational and Mendelian randomization (MR) studies link obesity and cancer, but it remains unclear whether these depend upon related metabolic abnormalities.

METHODS

We used information from 321,472 participants in the UK biobank, including 30,561 cases of obesity-related cancer. We constructed three genetic instruments reflecting higher adiposity together with either "unfavourable" (82 SNPs), "favourable" (24 SNPs) or "neutral" metabolic profile (25 SNPs). We looked at associations with 14 types of cancer, previously suggested to be associated with obesity.

RESULTS

All genetic instruments had a strong association with BMI (p < 1 × 10^-300 for all). The instrument reflecting unfavourable adiposity was also associated with higher CRP, HbA1c and adverse lipid profile, while instrument reflecting metabolically favourable adiposity was associated with lower HbA1c and a favourable lipid profile. In MR-inverse-variance weighted analysis unfavourable adiposity was associated with an increased risk of non-hormonal cancers (OR = 1.22, 95% confidence interval [CI]:1.08, 1.38), but a lower risk of hormonal cancers (OR = 0.80, 95%CI: 0.72, 0.89). From individual cancers, MR analyses suggested causal increases in the risk of multiple myeloma (OR = 1.36, 95%CI: 1.09, 1.70) and endometrial cancer (OR = 1.77, 95%CI: 1.16, 2.68) by greater genetically instrumented unfavourable adiposity but lower risks of breast and prostate cancer (OR = 0.72, 95%CI: 0.61, 0.83 and OR = 0.81, 95%CI: 0.68, 0.97, respectively). Favourable or neutral adiposity were not associated with the odds of any individual cancer.

CONCLUSIONS

Higher adiposity associated with a higher risk of non-hormonal cancer but a lower risk of some hormone related cancers. Presence of metabolic abnormalities might aggravate the adverse effects of higher adiposity on cancer. Further studies are warranted to investigate whether interventions on adverse metabolic health may help to alleviate obesity-related cancer risk.

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.

The genetic basis of major depression

Major depressive disorder (MDD) is a common, debilitating, phenotypically heterogeneous disorder with heritability ranges from 30% to 50%. Compared to other psychiatric disorders, its high prevalence, moderate heritability, and strong polygenicity have posed major challenges for gene-mapping in MDD. Studies of common genetic variation in MDD, driven by large international collaborations such as the Psychiatric Genomics Consortium, have confirmed the highly polygenic nature of the disorder and implicated over 100 genetic risk loci to date. Rare copy number variants associated with MDD risk were also recently identified. The goal of this review is to present a broad picture of our current understanding of the epidemiology, genetic epidemiology, molecular genetics, and gene-environment interplay in MDD. Insights into the impact of genetic factors on the aetiology of this complex disorder hold great promise for improving clinical care.

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.

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.

Gene-environment dependencies lead to collider bias in models with polygenic scores

The application of polygenic scores has transformed our ability to investigate whether and how genetic and environmental factors jointly contribute to the variation of complex traits. Modelling the complex interplay between genes and environment, however, raises serious methodological challenges. Here we illustrate the largely unrecognised impact of gene-environment dependencies on the identification of the effects of genes and their variation across environments. We show that controlling for heritable covariates in regression models that include polygenic scores as independent variables introduces endogenous selection bias when one or more of these covariates depends on unmeasured factors that also affect the outcome. This results in the problem of conditioning on a collider, which in turn leads to spurious associations and effect sizes. Using graphical and simulation methods we demonstrate that the degree of bias depends on the strength of the gene-covariate correlation and of hidden heterogeneity linking covariates with outcomes, regardless of whether the main analytic focus is mediation, confounding, or gene × covariate (commonly gene × environment) interactions. We offer potential solutions, highlighting the importance of causal inference. We also urge further caution when fitting and interpreting models with polygenic scores and non-exogenous environments or phenotypes and demonstrate how spurious associations are likely to arise, advancing our understanding of such results.

THI Modulation of Genetic and Non-genetic Variance Components for Carcass Traits in Hanwoo Cattle

The phenotype of carcass traits in beef cattle are affected by random genetic and non-genetic effects, which both can be modulated by an environmental variable such as Temperature-Humidity Index (THI), a key environmental factor in cattle production. In this study, a multivariate reaction norm model (MRNM) was used to assess if the random genetic and non-genetic (i.e., residual) effects of carcass weight (CW), back fat thickness (BFT), eye muscle area (EMA), and marbling score (MS) were modulated by THI, using 9,318 Hanwoo steers (N = 8,964) and cows (N = 354) that were genotyped on the Illumina Bovine SNP50 BeadChip (50K). THI was measured based on the period of 15–45 days before slaughter. Both the correlation and the interaction between THI and random genetic and non-genetic effects were accounted for in the model. In the analyses, it was shown that the genetic effects of EMA and the non-genetic effects of CW and MS were significantly modulated by THI. No significant THI modulation of such effects was found for BFT. These results highlight the relevance of THI changes for the genetic and non-genetic variation of CW, EMA, and MS in Hanwoo beef cattle. Importantly, heritability estimates for CW, EMA, and MS from additive models without considering THI interactions were underestimated. Moreover, the significance of interaction can be biased if not properly accounting for the correlation between THI and genetic and non-genetic effects. Thus, we argue that the estimation of genetic parameters should be based on appropriate models to avoid any potential bias of estimates. Our finding should serve as a basis for future studies aiming at revealing genotype by environment interaction in estimation and genomic prediction of breeding values.

Estimation of dynamic SNP-heritability with Bayesian Gaussian process models

Motivation

Improved DNA technology has made it practical to estimate single-nucleotide polymorphism (SNP)-heritability among distantly related individuals with unknown relationships. For growth- and development-related traits, it is meaningful to base SNP-heritability estimation on longitudinal data due to the time-dependency of the process. However, only few statistical methods have been developed so far for estimating dynamic SNP-heritability and quantifying its full uncertainty.

Results

We introduce a completely tuning-free Bayesian Gaussian process (GP)-based approach for estimating dynamic variance components and heritability as their function. For parameter estimation, we use a modern Markov Chain Monte Carlo method which allows full uncertainty quantification. Several datasets are analysed and our results clearly illustrate that the 95% credible intervals of the proposed joint estimation method (which ‘borrows strength’ from adjacent time points) are significantly narrower than of a two-stage baseline method that first estimates the variance components at each time point independently and then performs smoothing. We compare the method with a random regression model using MTG2 and BLUPF90 software and quantitative measures indicate superior performance of our method. Results are presented for simulated and real data with up to 1000 time points. Finally, we demonstrate scalability of the proposed method for simulated data with tens of thousands of individuals.

Availability and implementation

The C++ implementation dynBGP and simulated data are available in GitHub: https://github.com/aarjas/dynBGP. The programmes can be run in R. Real datasets are available in QTL archive: https://phenome.jax.org/centers/QTLA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Adaptability and stability analyses of plants using random regression models

The evaluation of cultivars using multi-environment trials (MET) is an important step in plant breeding programs. One of the objectives of these evaluations is to understand the genotype by environment interaction (GEI). A method of determining the effect of GEI on the performance of cultivars is based on studies of adaptability and stability. Initial studies were based on linear regression; however, these methodologies have limitations, mainly in trials with genetic or statistical unbalanced, heterogeneity of residual variances, and genetic covariance. An alternative would be the use of random regression models (RRM), in which the behavior of the genotypes is characterized as a reaction norm using longitudinal data or repeated measurements and information regarding a covariance function. The objective of this work was the application of RRM in the study of the behavior of common bean cultivars using a MET, based on Legendre polynomials and genotype-ideotype distances. We used a set of 13 trials, which were classified as unfavorable or favorable environments. The results revealed that RRM enables the prediction of the genotypic values of cultivars in environments where they were not evaluated with high accuracy values, thereby circumventing the unbalanced of the experiments. From these values, it was possible to measure the genotypic adaptability according to ideotypes, according to their reaction norms. In addition, the stability of the cultivars can be interpreted as variation in the behavior of the ideotype. The use of ideotypes based on real data allowed a better comparison of the performance of cultivars across environments. The use of RRM in plant breeding is a good alternative to understand the behavior of cultivars in a MET, especially when we want to quantify the adaptability and stability of genotypes.

Detecting Genotype-Population Interaction Effects by Ancestry Principal Components

Heterogeneity in the phenotypic mean and variance across populations is often observed for complex traits. One way to understand heterogeneous phenotypes lies in uncovering heterogeneity in genetic effects. Previous studies on genetic heterogeneity across populations were typically based on discrete groups in populations stratified by different countries or cohorts, which ignored the difference of population characteristics for the individuals within each group and resulted in loss of information. Here, we introduce a novel concept of genotype-by-population (G × P) interaction where population is defined by the first and second ancestry principal components (PCs), which are less likely to be confounded with country/cohort-specific factors. We applied a reaction norm model fitting each of 70 complex traits with significant SNP-heritability and the PCs as covariates to examine G × P interactions across diverse populations including white British and other white Europeans from the UK Biobank (N = 22,229). Our results demonstrated a significant population genetic heterogeneity for behavioral traits such as age at first sexual intercourse and academic qualification. Our approach may shed light on the latent genetic architecture of complex traits that underlies the modulation of genetic effects across different populations.

Whole-Genome Approach Discovers Novel Genetic and Nongenetic Variance Components Modulated by Lifestyle for Cardiovascular Health

Background

Both genetic and nongenetic factors can predispose individuals to cardiovascular risk. Finding ways to alter these predispositions is important for cardiovascular disease prevention.

Methods and Results

We used a novel whole‐genome approach to estimate the genetic and nongenetic effects on—and hence their predispositions to—cardiovascular risk and determined whether they vary with respect to lifestyle factors such as physical activity, smoking, alcohol consumption, and dietary intake. We performed analyses on the ARIC (Atherosclerosis Risk in Communities) Study (N=6896–7180) and validated findings using the UKBB (UK Biobank, N=14 076–34 538). Lifestyle modulation was evident for many cardiovascular traits such as body mass index and resting heart rate. For example, alcohol consumption modulated both genetic and nongenetic effects on body mass index, whereas smoking modulated nongenetic effects on heart rate, pulse pressure, and white blood cell count. We also stratified individuals according to estimated genetic and nongenetic effects that are modulated by lifestyle factors and showed distinct phenotype–lifestyle relationships across the stratified groups. Finally, we showed that neglecting lifestyle modulations of cardiovascular traits would on average reduce single nucleotide polymorphism heritability estimates of these traits by a small yet significant amount, primarily owing to the overestimation of residual variance.

Conclusions

Lifestyle changes are relevant to cardiovascular disease prevention. Individual differences in the genetic and nongenetic effects that are modulated by lifestyle factors, as shown by the stratified group analyses, implies a need for personalized lifestyle interventions. In addition, single nucleotide polymorphism–based heritability of cardiovascular traits without accounting for lifestyle modulations could be underestimated.

A Robust Method Uncovers Significant Context-Specific Heritability in Diverse Complex Traits

Gene-environment interactions (GxE) can be fundamental in applications ranging from functional genomics to precision medicine and is a conjectured source of substantial heritability. However, unbiased methods to profile GxE genome-wide are nascent and, as we show, cannot accommodate general environment variables, modest sample sizes, heterogeneous noise, and binary traits. To address this gap, we propose a simple, unifying mixed model for gene-environment interaction (GxEMM). In simulations and theory, we show that GxEMM can dramatically improve estimates and eliminate false positives when the assumptions of existing methods fail. We apply GxEMM to a range of human and model organism datasets and find broad evidence of context-specific genetic effects, including GxSex, GxAdversity, and GxDisease interactions across thousands of clinical and molecular phenotypes. Overall, GxEMM is broadly applicable for testing and quantifying polygenic interactions, which can be useful for explaining heritability and invaluable for determining biologically relevant environments.

GWEHS: A Genome-Wide Effect Sizes and Heritability Screener

During the last decade, there has been a huge development of Genome-Wide Association Studies (GWAS), and thousands of loci associated to complex traits have been detected. These efforts have led to the creation of public databases of GWAS results, making a huge source of information available on the genetic background of many diverse traits. Here we present GWEHS (Genome-Wide Effect size and Heritability Screener), an open-source online application to screen loci associated to human complex traits and diseases from the NHGRI-EBI GWAS Catalog. This application provides a way to explore the distribution of effect sizes of loci affecting these traits, as well as their contribution to heritability. Furthermore, it allows for making predictions on the change in the expected mean effect size, as well as in the heritability as new loci are found. The application enables inferences on whether the additive contribution of loci expected to be discovered in the future will be able to explain the estimates of familial heritability for the different traits. We illustrate the use of this tool, compare some of the results obtained with those from a previous meta-analysis, and discuss its uses and limitations.

Inferring the Nature of Missing Heritability in Human Traits Using Data from the GWAS Catalog

Thousands of genes responsible for many diseases and other common traits in humans have been detected by Genome Wide Association Studies (GWAS) in the last decade. However, candidate causal variants found so far usually explain only a small fraction of the heritability estimated by family data. The most common explanation for this observation is that the missing heritability corresponds to variants, either rare or common, with very small effect, which pass undetected due to a lack of statistical power. We carried out a meta-analysis using data from the NHGRI-EBI GWAS Catalog in order to explore the observed distribution of locus effects for a set of 42 complex traits and to quantify their contribution to narrow-sense heritability. With the data at hand, we were able to predict the expected distribution of locus effects for 16 traits and diseases, their expected contribution to heritability, and the missing number of loci yet to be discovered to fully explain the familial heritability estimates. Our results indicate that, for 6 out of the 16 traits, the additive contribution of a great number of loci is unable to explain the familial (broad-sense) heritability, suggesting that the gap between GWAS and familial estimates of heritability may not ever be closed for these traits. In contrast, for the other 10 traits, the additive contribution of hundreds or thousands of loci yet to be found could potentially explain the familial heritability estimates, if this were the case. Computer simulations are used to illustrate the possible contribution from nonadditive genetic effects to the gap between GWAS and familial estimates of heritability.