Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank

Gene-environment interactions within a precision environmental health framework

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

Genetic influences on brain and cognitive health and their interactions with cardiovascular conditions and depression

Approximately 40% of dementia cases could be prevented or delayed by modifiable risk factors related to lifestyle and environment. These risk factors, such as depression and vascular disease, do not affect all individuals in the same way, likely due to inter-individual differences in genetics. However, the precise nature of how genetic risk profiles interact with modifiable risk factors to affect brain health is poorly understood. Here we combine multiple data resources, including genotyping and postmortem gene expression, to map the genetic landscape of brain structure and identify 367 loci associated with cortical thickness and 13 loci associated with white matter hyperintensities (P < 5×10-8), with several loci also showing a significant association with cognitive function. We show that among 220 unique genetic loci associated with cortical thickness in our genome-wide association studies (GWAS), 95 also showed evidence of interaction with depression or cardiovascular conditions. Polygenic risk scores based on our GWAS of inferior frontal thickness also interacted with hypertension in predicting executive function in the Canadian Longitudinal Study on Aging. These findings advance our understanding of the genetic underpinning of brain structure and show that genetic risk for brain and cognitive health is in part moderated by treatable mid-life factors.

Genetic Determinants of Thiazide-Induced Hyperuricemia, Hyperglycemia, and Urinary Electrolyte Disturbances - A Genome-Wide Evaluation of the UK Biobank

Thiazide diuretics, widely used in hypertension, cause a variety of adverse reactions, including hyperglycemia, hyperuricemia, and electrolyte abnormalities. In this study, we aimed to identify genetic variants that interact with thiazide-use to increase the risk of these adverse reactions. Using UK Biobank data, we first performed genomewide variance quantitative trait locus (vQTL) analysis of ~ 6.2 million SNPs on 95,493 unrelated hypertensive White British participants (24,313 on self-reported bendroflumethiazide treatment at recruitment) for 2 blood (glucose and urate) and 2 urine (potassium and sodium) biomarkers. Second, we conducted direct gene-environment interaction (GEI) tests on the significant (P < 2.5 × 10-9) vQTLs, included a second UK Biobank cohort comprising 13,647 unrelated hypertensive White British participants (3,478 on thiazides other than bendroflumethiazide) and set significance at P = 0.05 divided by the number of vQTL SNPs tested for GEIs. The vQTL analysis identified eight statistically significant SNPs for blood glucose (5 SNPs) and serum urate (3 SNPs), with none being identified for the urinary biomarkers. Two of the SNPs (1 glucose SNP: CDKAL1 intron rs35612982, GEI P = 6.24 × 10-3; and 1 serum urate SNP: SLC2A9 intron rs938564, GEI P = 4.51 × 10-4) demonstrated significant GEI effects in the first, but not the second, cohort. Both genes are biologically plausible candidates, with the SLC2A9-mediated interaction having been previously reported. In conclusion, we used a two-stage approach to detect two biologically plausible genetic loci that can interact with thiazides to increase the risk of thiazide-associated biochemical abnormalities. Understanding how environmental exposures (including medications such as thiazides) and genetics interact, is an important step toward precision medicine and improved patient outcomes.

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.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

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.

No gene by stressful life events interaction on individual differences in adults' self-control

Background

Difficulty with self-control, or the ability to alter impulses and behavior in a goal-directed way, predicts interpersonal conflict, lower socioeconomic attainments, and more adverse health outcomes. Etiological understanding, and intervention for low self-control is, therefore, a public health goal. A prominent developmental theory proposes that individuals with high genetic propensity for low self-control that are also exposed to stressful environments may be most at-risk of low levels of self-control. Here we examine if polygenic measures associated with behaviors marked by low self-control interact with stressful life events in predicting self-control.

Methods

Leveraging molecular data from a large population-based Dutch sample (N = 7,090, Mage = 41.2) to test for effects of genetics (i.e., polygenic scores for ADHD and aggression), stressful life events (e.g., traffic accident, violent assault, financial problems), and a gene-by-stress interaction on self-control (measured with the ASEBA Self-Control Scale).

Results

Both genetics (β =.03 -.04, p <.001) and stressful life events (β = .11 -.14, p <.001) were associated with individual differences in self-control. We find no evidence of a gene-by-stressful life events interaction on individual differences in adults' self-control.

Conclusion

Our findings are consistent with the notion that genetic influences and stressful life events exert largely independent effects on adult self-control. However, the small effect sizes of polygenic scores increases the likelihood of null results. Genetically-informed longitudinal research in large samples can further inform the etiology of individual differences in self-control from early childhood into later adulthood and its downstream implications for public health.

Genome-wide analyses of variance in blood cell phenotypes provide new insights into complex trait biology and prediction

Blood cell phenotypes are routinely tested in healthcare to inform clinical decisions. Genetic variants influencing mean blood cell phenotypes have been used to understand disease aetiology and improve prediction; however, additional information may be captured by genetic effects on observed variance. Here, we mapped variance quantitative trait loci (vQTL), i.e. genetic loci associated with trait variance, for 29 blood cell phenotypes from the UK Biobank (N~408,111). We discovered 176 independent blood cell vQTLs, of which 147 were not found by additive QTL mapping. vQTLs displayed on average 1.8-fold stronger negative selection than additive QTL, highlighting that selection acts to reduce extreme blood cell phenotypes. Variance polygenic scores (vPGSs) were constructed to stratify individuals in the INTERVAL cohort (N~40,466), where genetically less variable individuals (low vPGS) had increased conventional PGS accuracy (by ~19%) than genetically more variable individuals. Genetic prediction of blood cell traits improved by ~10% on average combining PGS with vPGS. Using Mendelian randomisation and vPGS association analyses, we found that alcohol consumption significantly increased blood cell trait variances highlighting the utility of blood cell vQTLs and vPGSs to provide novel insight into phenotype aetiology as well as improve prediction.

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Using Genetics to Inform Interventions Related to Sodium and Potassium in Hypertension

BACKGROUND

Hypertension is a key risk factor for major adverse cardiovascular events but remains difficult to treat in many individuals. Dietary interventions are an effective approach to lower blood pressure (BP) but are not equally effective across all individuals. BP is heritable, and genetics may be a useful tool to overcome treatment response heterogeneity. We investigated whether the genetics of BP could be used to identify individuals with hypertension who may receive a particular benefit from lowering sodium intake and boosting potassium levels.

METHODS

In this observational genetic study, we leveraged cross-sectional data from up to 296 475 genotyped individuals drawn from the UK Biobank cohort for whom BP and urinary electrolytes (sodium and potassium), biomarkers of sodium and potassium intake, were measured. Biologically directed genetic scores for BP were constructed specifically among pathways related to sodium and potassium biology (pharmagenic enrichment scores), as well as unannotated genome-wide scores (conventional polygenic scores). We then tested whether there was a gene-by-environment interaction between urinary electrolytes and these genetic scores on BP.

RESULTS

Genetic risk and urinary electrolytes both independently correlated with BP. However, urinary sodium was associated with a larger BP increase among individuals with higher genetic risk in sodium- and potassium-related pathways than in those with comparatively lower genetic risk. For example, each SD in urinary sodium was associated with a 1.47-mm Hg increase in systolic BP for those in the top 10% of the distribution of genetic risk in sodium and potassium transport pathways versus a 0.97-mm Hg systolic BP increase in the lowest 10% (P=1.95×10-3). This interaction with urinary sodium remained when considering estimated glomerular filtration rate and indexing sodium to urinary creatinine. There was no strong evidence of an interaction between urinary sodium and a standard genome-wide polygenic score of BP.

CONCLUSIONS

The data suggest that genetic risk in sodium and potassium pathways could be used in a precision medicine model to direct interventions more specifically in the management of hypertension. Intervention studies are warranted.

Searching for gene-gene interactions through variance quantitative trait loci of 29 continuous Taiwan Biobank phenotypes

Introduction: After the era of genome-wide association studies (GWAS), thousands of genetic variants have been identified to exhibit main effects on human phenotypes. The next critical issue would be to explore the interplay between genes, the so-called "gene-gene interactions" (GxG) or epistasis. An exhaustive search for all single-nucleotide polymorphism (SNP) pairs is not recommended because this will induce a harsh penalty of multiple testing. Limiting the search of epistasis on SNPs reported by previous GWAS may miss essential interactions between SNPs without significant marginal effects. Moreover, most methods are computationally intensive and can be challenging to implement genome-wide. Methods: I here searched for GxG through variance quantitative trait loci (vQTLs) of 29 continuous Taiwan Biobank (TWB) phenotypes. A discovery cohort of 86,536 and a replication cohort of 25,460 TWB individuals were analyzed, respectively. Results: A total of 18 nearly independent vQTLs with linkage disequilibrium measure r 2 < 0.01 were identified and replicated from nine phenotypes. 15 significant GxG were found with p-values <1.1E-5 (in the discovery cohort) and false discovery rates <2% (in the replication cohort). Among these 15 GxG, 11 were detected for blood traits including red blood cells, hemoglobin, and hematocrit; 2 for total bilirubin; 1 for fasting glucose; and 1 for total cholesterol (TCHO). All GxG were observed for gene pairs on the same chromosome, except for the APOA5 (chromosome 11)-TOMM40 (chromosome 19) interaction for TCHO. Discussion: This study provided a computationally feasible way to search for GxG genome-wide and applied this approach to 29 phenotypes.

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification

BACKGROUND

Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored.

METHODS

Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold.

RESULTS

Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10-15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10-46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74).

CONCLUSIONS

Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

Age-dependent genes in adipose stem and precursor cells affect regulation of fat cell differentiation and link aging to obesity via cellular and genetic interactions

BACKGROUND

Age and obesity are dominant risk factors for several common cardiometabolic disorders, and both are known to impair adipose tissue function. However, the underlying cellular and genetic factors linking aging and obesity on adipose tissue function have remained elusive. Adipose stem and precursor cells (ASPCs) are an understudied, yet crucial adipose cell type due to their deterministic adipocyte differentiation potential, which impacts the capacity to store fat in a metabolically healthy manner.

METHODS

We integrated subcutaneous adipose tissue (SAT) bulk (n=435) and large single-nucleus RNA sequencing (n=105) data with the UK Biobank (UKB) (n=391,701) data to study age-obesity interactions originating from ASPCs by performing cell-type decomposition, differential expression testing, cell-cell communication analyses, and construction of polygenic risk scores for body mass index (BMI).

RESULTS

We found that the SAT ASPC proportions significantly decrease with age in an obesity-dependent way consistently in two independent cohorts, both showing that the age dependency of ASPC proportions is abolished by obesity. We further identified 76 genes (72 SAT ASPC marker genes and 4 transcription factors regulating ASPC marker genes) that are differentially expressed by age in SAT and functionally enriched for developmental processes and adipocyte differentiation (i.e., adipogenesis). The 76 age-perturbed ASPC genes include multiple negative regulators of adipogenesis, such as RORA, SMAD3, TWIST2, and ZNF521, form tight clusters of longitudinally co-expressed genes during human adipogenesis, and show age-based differences in cellular interactions between ASPCs and adipose cell types. Finally, our genetic data demonstrate that cis-regional variants of these genes interact with age as predictors of BMI in an obesity-dependent way in the large UKB, while no such gene-age interaction on BMI is observed with non-age-dependent ASPC marker genes, thus independently confirming our cellular ASPC results at the biobank level.

CONCLUSIONS

Overall, we discover that obesity prematurely induces a decrease in ASPC proportions and identify 76 developmentally important ASPC genes that implicate altered negative regulation of fat cell differentiation as a mechanism for aging and directly link aging to obesity via significant cellular and genetic interactions.

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.

X chromosome-wide association study of quantitative biomarkers from the Alzheimer's Disease Neuroimaging Initiative study

Introduction

Alzheimer's disease (AD) is a complex neurodegenerative disease with high heritability. Compared to autosomes, a higher proportion of disorder-associated genes on X chromosome are expressed in the brain. However, only a few studies focused on the identification of the susceptibility loci for AD on X chromosome.

Methods

Using the data from the Alzheimer's Disease Neuroimaging Initiative Study, we conducted an X chromosome-wide association study between 16 AD quantitative biomarkers and 19,692 single nucleotide polymorphisms (SNPs) based on both the cross-sectional and longitudinal studies.

Results

We identified 15 SNPs statistically significantly associated with different quantitative biomarkers of the AD. For the cross-sectional study, six SNPs (rs5927116, rs4596772, rs5929538, rs2213488, rs5920524, and rs5945306) are located in or near to six genes DMD, TBX22, LOC101928437, TENM1, SPANXN1, and ZFP92, which have been reported to be associated with schizophrenia or neuropsychiatric diseases in literature. For the longitudinal study, four SNPs (rs4829868, rs5931111, rs6540385, and rs763320) are included in or near to two genes RAC1P4 and AFF2, which have been demonstrated to be associated with brain development or intellectual disability in literature, while the functional annotations of other five novel SNPs (rs12157031, rs428303, rs5953487, rs10284107, and rs5955016) have not been found.

Discussion

15 SNPs were found statistically significantly associated with the quantitative biomarkers of the AD. Follow-up study in molecular genetics is needed to verify whether they are indeed related to AD. The findings in this article expand our understanding of the role of the X chromosome in exploring disease susceptibility, introduce new insights into the molecular genetics behind the AD, and may provide a mechanistic clue to further AD-related studies.

Cross-ancestry analyses identify new genetic loci associated with 25-hydroxyvitamin D

Vitamin D status-a complex trait influenced by environmental and genetic factors-is tightly associated with skin colour and ancestry. Yet very few studies have investigated the genetic underpinnings of vitamin D levels across diverse ancestries, and the ones that have, relied on small sample sizes, resulting in inconclusive results. Here, we conduct genome-wide association studies (GWAS) of 25 hydroxyvitamin D (25OHD)-the main circulating form of vitamin D-in 442,435 individuals from four broad genetically-determined ancestry groups represented in the UK Biobank: European (N = 421,867), South Asian (N = 9,983), African (N = 8,306) and East Asian (N = 2,279). We identify a new genetic determinant of 25OHD (rs146759773) in individuals of African ancestry, which was not detected in previous analysis of much larger European cohorts due to low minor allele frequency. We show genome-wide significant evidence of dominance effects in 25OHD that protect against vitamin D deficiency. Given that key events in the synthesis of 25OHD occur in the skin and are affected by pigmentation levels, we conduct GWAS of 25OHD stratified by skin colour and identify new associations. Lastly, we test the interaction between skin colour and variants associated with variance in 25OHD levels and identify two loci (rs10832254 and rs1352846) whose association with 25OHD differs in individuals of distinct complexions. Collectively, our results provide new insights into the complex relationship between 25OHD and skin colour and highlight the importance of diversity in genomic studies. Despite the much larger rates of vitamin D deficiency that we and others report for ancestry groups with dark skin (e.g., South Asian), our study highlights the importance of considering ancestral background and/or skin colour when assessing the implications of low vitamin D.

Complex effects of sequence variants on lipid levels and coronary artery disease

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

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.

Applying an evolutionary mismatch framework to understand disease susceptibility

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

Mendelian randomization mixed-scale treatment effect robust identification and estimation for causal inference

Standard Mendelian randomization (MR) analysis can produce biased results if the genetic variant defining an instrumental variable (IV) is confounded and/or has a horizontal pleiotropic effect on the outcome of interest not mediated by the treatment variable. We provide novel identification conditions for the causal effect of a treatment in the presence of unmeasured confounding by leveraging a possibly invalid IV for which both the IV independence and exclusion restriction assumptions may be violated. The proposed Mendelian randomization mixed-scale treatment effect robust identification (MR MiSTERI) approach relies on (i) an assumption that the treatment effect does not vary with the possibly invalid IV on the additive scale; (ii) that the confounding bias does not vary with the possibly invalid IV on the odds ratio scale; and (iii) that the residual variance for the outcome is heteroskedastic with respect to the possibly invalid IV. Although assumptions (i) and (ii) have, respectively, appeared in the IV literature, assumption (iii) has not; we formally establish that their conjunction can identify a causal effect even with an invalid IV. MR MiSTERI is shown to be particularly advantageous in the presence of pervasive heterogeneity of pleiotropic effects on the additive scale. We propose a simple and consistent three-stage estimator that can be used as a preliminary estimator to a carefully constructed efficient one-step-update estimator. In order to incorporate multiple, possibly correlated, and weak invalid IVs, a common challenge in MR studies, we develop a MAny Weak Invalid Instruments (MR MaWII MiSTERI) approach for strengthened identification and improved estimation accuracy. Both simulation studies and UK Biobank data analysis results demonstrate the robustness of the proposed methods.

Quantile Regression for biomarkers in the UK Biobank

Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. An alternative and easy to apply approach is quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest by modeling conditional quantiles within a regression framework. Quantile regression can be applied efficiently at biobank scale using standard statistical packages in much the same way as linear regression, while having some unique advantages such as identifying variants with heterogeneous effects across different quantiles, including non-additive effects and variants involved in gene-environment interactions; accommodating a wide range of phenotype distributions with invariance to trait transformation; and overall providing more detailed information about the underlying genotype-phenotype associations. Here, we demonstrate the value of quantile regression in the context of GWAS by applying it to 39 quantitative traits in the UK Biobank (n > 300 , 000 individuals). Across these 39 traits we identify 7,297 significant loci, including 259 loci only detected by quantile regression. We show that quantile regression can help uncover replicable but unmodelled gene-environment interactions, and can provide additional key insights into poorly understood genotype-phenotype correlations for clinically relevant biomarkers at minimal additional cost.

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

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

A non-coding variant linked to metabolic obesity with normal weight affects actin remodelling in subcutaneous adipocytes

Recent large-scale genomic association studies found evidence for a genetic link between increased risk of type 2 diabetes and decreased risk for adiposity-related traits, reminiscent of metabolically obese normal weight (MONW) association signatures. However, the target genes and cellular mechanisms driving such MONW associations remain to be identified. Here, we systematically identify the cellular programmes of one of the top-scoring MONW risk loci, the 2q24.3 risk locus, in subcutaneous adipocytes. We identify a causal genetic variant, rs6712203, an intronic single-nucleotide polymorphism in the COBLL1 gene, which changes the conserved transcription factor motif of POU domain, class 2, transcription factor 2, and leads to differential COBLL1 gene expression by altering the enhancer activity at the locus in subcutaneous adipocytes. We then establish the cellular programme under the genetic control of the 2q24.3 MONW risk locus and the effector gene COBLL1, which is characterized by impaired actin cytoskeleton remodelling in differentiating subcutaneous adipocytes and subsequent failure of these cells to accumulate lipids and develop into metabolically active and insulin-sensitive adipocytes. Finally, we show that perturbations of the effector gene Cobll1 in a mouse model result in organismal phenotypes matching the MONW association signature, including decreased subcutaneous body fat mass and body weight along with impaired glucose tolerance. Taken together, our results provide a mechanistic link between the genetic risk for insulin resistance and low adiposity, providing a potential therapeutic hypothesis and a framework for future identification of causal relationships between genome associations and cellular programmes in other disorders.

Molecular mechanisms of environmental exposures and human disease

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

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.

Evolutionary mismatch and the role of GxE interactions in human disease

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

Gene-environment interactions and their impact on human health

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

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.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

Altered macronutrient composition and genetics influence the complex transcriptional network associated with adiposity in the Collaborative Cross

Background

Obesity is a serious disease with a complex etiology characterized by overaccumulation of adiposity resulting in detrimental health outcomes. Given the liver’s critical role in the biological processes that attenuate adiposity accumulation, elucidating the influence of genetics and dietary patterns on hepatic gene expression is fundamental for improving methods of obesity prevention and treatment. To determine how genetics and diet impact obesity development, mice from 22 strains of the genetically diverse recombinant inbred Collaborative Cross (CC) mouse panel were challenged to either a high-protein or high-fat high-sucrose diet, followed by extensive phenotyping and analysis of hepatic gene expression.

Results

Over 1000 genes differentially expressed by perturbed dietary macronutrient composition were enriched for biological processes related to metabolic pathways. Additionally, over 9000 genes were differentially expressed by strain and enriched for biological process involved in cell adhesion and signaling. Weighted gene co-expression network analysis identified multiple gene clusters (modules) associated with body fat % whose average expression levels were influenced by both dietary macronutrient composition and genetics. Each module was enriched for distinct types of biological functions.

Conclusions

Genetic background affected hepatic gene expression in the CC overall, but diet macronutrient differences also altered expression of a specific subset of genes. Changes in macronutrient composition altered gene expression related to metabolic processes, while genetic background heavily influenced a broad range of cellular functions and processes irrespective of adiposity. Understanding the individual role of macronutrient composition, genetics, and their interaction is critical to developing therapeutic strategies and policy recommendations for precision nutrition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12263-022-00714-x.

Genetic determinants of polygenic prediction accuracy within a population

Genomic risk prediction is on the emerging path toward personalized medicine. However, the accuracy of polygenic prediction varies strongly in different individuals. Based on up to 352,277 European ancestry participants in the UK Biobank, we constructed polygenic risk scores for 15 physiological and biochemical quantitative traits. We identified a total of 185 polygenic prediction variability quantitative trait loci for 11 traits by Levene's test among 254,376 unrelated individuals. We validated the effects of prediction variability quantitative trait loci using an independent test set of 58,927 individuals. For instance, a score aggregating 51 prediction variability quantitative trait locus variants for triglycerides had the strongest Spearman correlation of 0.185 (P-value <1.0 × 10-300) with the squared prediction errors. We found a strong enrichment of complex genetic effects conferred by prediction variability quantitative trait loci compared to risk loci identified in genome-wide association studies, including 89 prediction variability quantitative trait loci exhibiting dominance effects. Incorporation of dominance effects into polygenic risk scores significantly improved polygenic prediction for triglycerides, low-density lipoprotein cholesterol, vitamin D, and platelet. In conclusion, we have discovered and profiled genetic determinants of polygenic prediction variability for 11 quantitative biomarkers. These findings may assist interpretation of genomic risk prediction in various contexts and encourage novel approaches for constructing polygenic risk scores with complex genetic effects.

Association of demographic and clinical factors with risk of acute pancreatitis: An exposure-wide Mendelian randomization study

BACKGROUND

The incidence of acute pancreatitis (AP) is increasing over years, which brings enormous economy and health burden. However, the aetiologies of AP and underlying mechanisms are still unclear. Here, we performed a two-sample Mendelian randomization (MR) analysis to investigate the associations between all reported possible risk factors and AP using publicly available genome-wide association study summary statistics.

METHODS

A series of quality control steps were taken in our analysis to select eligible instrumental single nucleotide polymorphisms which were strongly associated with exposures. To make the conclusions more robust and reliable, we utilized several analytical methods (inverse-variance weighting, MR-PRESSO method, weighted median, MR-Egger regression) that are based on different assumptions of two-sample MR analysis. The MR-Egger intercept test, radial regression and leave-one-out sensitivity analysis were performed to evaluate the horizontal pleiotropy, heterogeneities, and stability of these genetic variants on each exposure. A two-step MR method was applied to explore mediators in significant associations.

RESULTS

Genetic predisposition to cholelithiasis (effect estimate: 17.30, 95% CI: 12.25-22.36, p = 1.95 E-11), body mass index (0.32, 95% CI: 0.13-0.51, p < 0.001), body fat percentage (0.57, 95% CI: 0.31-0.83, p = 1.31 E-05), trunk fat percentage (0.36, 95% CI: 0.14-0.59, p < 0.005), ever smoked (1.61, 95% CI: 0.45-2.77, p = 0.007), and limbs fat percentage (0.55, 95% CI: 0.41-0.69, p < 0.001) were associated with an increased risk of AP. In addition, whole-body fat-free mass (-0.32, 95% CI: -0.55 to -0.10, p = 0.004) was associated with a decrease risk of AP.

CONCLUSION

Genetic predisposition to cholelithiasis, obesity and smoking could be causally associated with an increased risk of AP, and whole body fat-free mass could be associated with a decreased risk of AP.

The effect of heteroscedasticity on the prediction efficiency of genome-wide polygenic score for body mass index

Globally, more than 1.9 billion adults are overweight. Thus, obesity is a serious public health issue. Moreover, obesity is a major risk factor for diabetes mellitus, coronary heart disease, and cardiovascular disease. Recently, GWAS examining obesity and body mass index (BMI) have increasingly unveiled many aspects of the genetic architecture of obesity and BMI. Information on genome-wide genetic variants has been used to estimate the genome-wide polygenic score (GPS) for a personalized prediction of obesity. However, the prediction power of GPS is affected by various factors, including the unequal variance in the distribution of a phenotype, known as heteroscedasticity. Here, we calculated a GPS for BMI using LDpred2, which was based on the BMI GWAS summary statistics from a European meta-analysis. Then, we tested the GPS in 354,761 European samples from the UK Biobank and found an effective prediction power of the GPS on BMI. To study a change in the variance of BMI, we investigated the heteroscedasticity of BMI across the GPS via graphical and statistical methods. We also studied the homoscedastic samples for BMI compared to the heteroscedastic sample, randomly selecting samples with various standard deviations of BMI residuals. Further, we examined the effect of the genetic interaction of GPS with environment (GPS×E) on the heteroscedasticity of BMI. We observed the changing variance (i.e., heteroscedasticity) of BMI along the GPS. The heteroscedasticity of BMI was confirmed by both the Breusch-Pagan test and the Score test. Compared to the heteroscedastic sample, the homoscedastic samples from small standard deviation of BMI residuals showed a decreased heteroscedasticity and an improved prediction accuracy, suggesting a quantitatively negative correlation between the phenotypic heteroscedasticity and the prediction accuracy of GPS. To further test the effects of the GPS×E on heteroscedasticity, first we tested the genetic interactions of the GPS with 21 environments and found 8 significant GPS×E interactions on BMI. However, the heteroscedasticity of BMI was not ameliorated after adjusting for the GPS×E interactions. Taken together, our findings suggest that the heteroscedasticity of BMI exists along the GPS and is not affected by the GPS×E interaction.

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

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

A population-wide gene-environment interaction study on how genes, schools, and residential areas shape achievement

A child's environment is thought to be composed of different levels that interact with their individual genetic propensities. However, studies have not tested this theory comprehensively across multiple environmental levels. Here, we quantify the contributions of child, parent, school, neighbourhood, district, and municipality factors to achievement, and investigate interactions between polygenic indices for educational attainment (EA-PGI) and environmental levels. We link population-wide administrative data on children's standardised test results, schools and residential identifiers to the Norwegian Mother, Father, and Child Cohort Study (MoBa), which includes >23,000 genotyped parent-child trios. We test for gene-environment interactions using multilevel models with interactions between EA-PGI and random effects for school and residential environments (thus remaining agnostic to specific features of environments). We use parent EA-PGI to control for gene-environment correlation. We found an interaction between students' EA-PGI and schools suggesting compensation: higher-performing schools can raise overall achievement without leaving children with lower EA-PGI behind. Differences between schools matter more for students with lower EA-PGI, explaining 4 versus 2% of the variance in achievement for students 2 SD below versus 2 SD above the mean EA-PGI. Neighbourhood, district, and municipality variation contribute little to achievement (<2% of the variance collectively), and do not interact with children's individual EA-PGI. Policy to reduce social inequality in achievement in Norway should focus on tackling unequal support across schools for children with difficulties.

Considering strategies for SNP selection in genetic and polygenic risk scores

Genetic risk scores (GRS) and polygenic risk scores (PRS) are weighted sums of, respectively, several or many genetic variant indicator variables. Although they are being increasingly proposed for clinical use, the best ways to construct them are still actively debated. In this commentary, we present several case studies illustrating practical challenges associated with building or attempting to improve score performance when there is expected to be heterogeneity of disease risk between cohorts or between subgroups of individuals. Specifically, we contrast performance associated with several ways of selecting single nucleotide polymorphisms (SNPs) for inclusion in these scores. By considering GRS and PRS as predictors that are measured with error, insights into their strengths and weaknesses may be obtained, and SNP selection approaches play an important role in defining such errors.

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.

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

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

A quantile integral linear model to quantify genetic effects on phenotypic variability

Detecting genetic variants associated with the variance of complex traits, that is, variance quantitative trait loci (vQTLs), can provide crucial insights into the interplay between genes and environments and how they jointly shape human phenotypes in the population. We propose a quantile integral linear model (QUAIL) to estimate genetic effects on trait variability. Through extensive simulations and analyses of real data, we demonstrate that QUAIL provides computationally efficient and statistically powerful vQTL mapping that is robust to non-Gaussian phenotypes and confounding effects on phenotypic variability. Applied to UK Biobank (n = 375,791), QUAIL identified 11 vQTLs for body mass index (BMI) that have not been previously reported. Top vQTL findings showed substantial enrichment for interactions with physical activities and sedentary behavior. Furthermore, variance polygenic scores (vPGSs) based on QUAIL effect estimates showed superior predictive performance on both population-level and within-individual BMI variability compared to existing approaches. Overall, QUAIL is a unified framework to quantify genetic effects on the phenotypic variability at both single-variant and vPGS levels. It addresses critical limitations in existing approaches and may have broad applications in future gene-environment interaction studies.

G × EBLUP: A novel method for exploring genotype by environment interactions and genomic prediction

Genotype by environment (G × E) interaction is fundamental in the biology of complex traits and diseases. However, most of the existing methods for genomic prediction tend to ignore G × E interaction (GEI). In this study, we proposed the genomic prediction method G × EBLUP by considering GEI. Meanwhile, G × EBLUP can also detect the genome-wide single nucleotide polymorphisms (SNPs) subject to GEI. Using comprehensive simulations and analysis of real data from pigs and maize, we showed that G × EBLUP achieved higher efficiency in mapping GEI SNPs and higher prediction accuracy than the existing methods, and its superiority was more obvious when the GEI variance was large. For pig and maize real data, compared with GBLUP, G × EBLUP showed improvement by 3% in the prediction accuracy for backfat thickness, while our findings indicated that the trait of days to 100 kg of pig was not affected by GEI and G × EBLUP did not improve the accuracy of genomic prediction for the trait. A significant advantage was observed for G × EBLUP in maize; the prediction accuracy was improved by ∼5.0 and 7.7% for grain weight and water content, respectively. Furthermore, G × EBLUP was not influenced by the number of environment levels. It could determine a favourable environment using SNP Bayes factors for each environment, implying that it is a robust and useful method for market-specific animal and plant breeding. We proposed G × EBLUP, a novel method for the estimation of genomic breeding value by considering GEI. This method identified the genome-wide SNPs that were susceptible to GEI and yielded higher genomic prediction accuracies and lower mean squared error compared with the GBLUP method.

Exploring Lead loci shared between schizophrenia and Cardiometabolic traits

Individuals with schizophrenia (SCZ) have, on average, a 10- to 20-year shorter expected life span than the rest of the population, primarily due to cardiovascular disease comorbidity. Genome-wide association studies (GWAS) have previously been used to separately identify common variants in SCZ and cardiometabolic traits. However, genetic variants jointly influencing both traits remain to be fully characterised. To assess overlaps (if any) between the genetic architecture of SCZ and cardiometabolic traits, we used conditional false discovery rate (FDR) and local genetic correlation statistical framework analyses. A conjunctional FDR was used to identify shared genetic traits between SCZ and cardiometabolic risk factors. We identified 144 genetic variants which were shared between SCZ and body mass index (BMI), and 15 variants shared between SCZ and triglycerides (TG). Furthermore, we discovered four novel single nucleotide polymorphisms (SNPs) (rs3865350, rs9860913, rs13307 and rs9614186) and four proximate genes (DERL2, SNX4, LY75 and EFCAB6) which were shared by SCZ and BMI. We observed that the novel genetic variant rs13307 and the most proximate gene LY75 exerted potential effects on SCZ and BMI comorbidity. Also, we observed a mixture of concordant and opposite direction associations with shared genetic variants. We demonstrated a moderate to high genetic overlap between SCZ and cardiometabolic traits associated with a pattern of bidirectional associations. Our data suggested a complex interplay between metabolism-related gene pathways in SCZ pathophysiology.

Nutrigenetic Interaction of Spontaneously Hypertensive Rat Chromosome 20 Segment and High-Sucrose Diet Sensitizes to Metabolic Syndrome

Several corresponding regions of human and mammalian genomes have been shown to affect sensitivity to the manifestation of metabolic syndrome via nutrigenetic interactions. In this study, we assessed the effect of sucrose administration in a newly established congenic strain BN.SHR20, in which a limited segment of rat chromosome 20 from a metabolic syndrome model, spontaneously hypertensive rat (SHR), was introgressed into Brown Norway (BN) genomic background. We mapped the extent of the differential segment and compared the genomic sequences of BN vs. SHR within the segment in silico. The differential segment of SHR origin in BN.SHR20 spans about 9 Mb of the telomeric portion of the short arm of chromosome 20. We identified non-synonymous mutations e.g., in ApoM, Notch4, Slc39a7, Smim29 genes and other variations in or near genes associated with metabolic syndrome in human genome-wide association studies. Male rats of BN and BN.SHR20 strains were fed a standard diet for 18 weeks (control groups) or 16 weeks of standard diet followed by 14 days of high-sucrose diet (HSD). We assessed the morphometric and metabolic profiles of all groups. Adiposity significantly increased only in BN.SHR20 after HSD. Fasting glycemia and the glucose levels during the oral glucose tolerance test were higher in BN.SHR20 than in BN groups, while insulin levels were comparable. The fasting levels of triacylglycerols were the highest in sucrose-fed BN.SHR20, both compared to the sucrose-fed BN and the control BN.SHR20. The non-esterified fatty acids and total cholesterol concentrations were higher in BN.SHR20 compared to their respective BN groups, and the HSD elicited an increase in non-esterified fatty acids only in BN.SHR20. In a new genetically defined model, we have isolated a limited genomic region involved in nutrigenetic sensitization to sucrose-induced metabolic disturbances.

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.

CNest: A novel copy number association discovery method uncovers 862 new associations from 200,629 whole-exome sequence datasets in the UK Biobank

Copy number variation (CNV) is known to influence human traits, having a rich history of research into common and rare genetic disease, and although CNV is accepted as an important class of genomic variation, progress on copy-number-based genome-wide association studies (GWASs) from next-generation sequencing (NGS) data has been limited. Here we present a novel method for large-scale copy number analysis from NGS data generating robust copy number estimates and allowing copy number GWASs (CN-GWASs) to be performed genome-wide in discovery mode. We provide a detailed analysis in the UK Biobank resource and a specifically designed software package. We use these methods to perform CN-GWAS analysis across 78 human traits, discovering over 800 genetic associations that are likely to contribute strongly to trait distributions. Finally, we compare CNV and SNP association signals across the same traits and samples, defining specific CNV association classes.

Quantitative trait loci, G×E and G×G for glycemic traits: response to metformin and placebo in the Diabetes Prevention Program (DPP)

The complex genetic architecture of type-2-diabetes (T2D) includes gene-by-environment (G×E) and gene-by-gene (G×G) interactions. To identify G×E and G×G, we screened markers for patterns indicative of interactions (relationship loci [rQTL] and variance heterogeneity loci [vQTL]). rQTL exist when the correlation between multiple traits varies by genotype and vQTL occur when the variance of a trait differs by genotype (potentially flagging G×G and G×E). In the metformin and placebo arms of the DPP (n = 1762) we screened 280,965 exomic and intergenic SNPs, for rQTL and vQTL patterns in association with year one changes from baseline in glycemia and related traits (insulinogenic index [IGI], insulin sensitivity index [ISI], fasting glucose and fasting insulin). Significant (p < 1.8 × 10^-7) rQTL and vQTL generated a priori hypotheses of individual G×E tests for a SNP × metformin treatment interaction and secondarily for G×G screens. Several rQTL and vQTL identified led to 6 nominally significant (p < 0.05) metformin treatment × SNP interactions (4 for IGI, one insulin, and one glucose) and 12G×G interactions (all IGI) that exceeded experiment-wide significance (p < 4.1 × 10^-9). Some loci are directly associated with incident diabetes, and others are rQTL and modify a trait's relationship with diabetes (2 diabetes/glucose, 2 diabetes/insulin, 1 diabetes/IGI). rs3197999, an ISI/insulin rQTL, is a possible gene damaging missense mutation in MST1, is associated with ulcerative colitis, sclerosing cholangitis, Crohn's disease, BMI and coronary artery disease. This study demonstrates evidence for context-dependent effects (G×G & G×E) and the complexity of these T2D-related traits.