Genome-wide gene-diet interaction analysis in the UK Biobank identifies novel effects on hemoglobin A1c

Interaction effect between low birthweight and resistin gene rs1862513 variant on insulin resistance and type 2 diabetes mellitus in adulthood: Toon Genome Study

AIMS/INTRODUCTION

Gene-environment interactions are considered to critically influence type 2 diabetes mellitus development; however, the underlying mechanisms and specific interactions remain unclear. Given the increasing prevalence of low birthweight (LBW) influenced by the intrauterine environment, we sought to investigate genetic factors related to type 2 diabetes development in individuals with LBW.

MATERIALS AND METHODS

The interaction between 20 reported type 2 diabetes susceptibility genes and the development of type 2 diabetes in LBW (<2,500 g) individuals in a population-based Japanese cohort (n = 1,021) was examined by logistic regression and stratified analyses.

RESULTS

Logistic regression analyses showed that only the G/G genotype at the rs1862513 locus of the resistin gene (RETN), an established initiator of insulin resistance, was closely related to the prevalence of type 2 diabetes in individuals with LBW. Age, sex and current body mass index-adjusted stratified analyses showed a significant interaction effect of LBW and the RETN G/G genotype on fasting insulin, homeostatic model assessment 2-insulin resistance, Matsuda index and the prevalence of type 2 diabetes (all P-values for interaction <0.05). The adjusted odds ratio for type 2 diabetes in the LBW + G/G genotype group was 7.33 (95% confidence interval 2.43-22.11; P = 0.002) compared with the non-LBW + non-G/G genotype group. Similar results were obtained after excluding the influence of malnutrition due to World War II.

CONCLUSIONS

Simultaneous assessment of LBW and the RETN G/G genotype can more accurately predict the risk of future type 2 diabetes than assessing each of these factors alone, and provide management strategies, including early lifestyle intervention in LBW population.

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.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

BACKGROUND

The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP.

METHODS

We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses.

RESULTS

We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

CONCLUSIONS

We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

Objective

We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP).

Methods

We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses.

Results

We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

Conclusion

We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

Investigating Gene-Diet Interactions Impacting the Association Between Macronutrient Intake and Glycemic Traits

Few studies have demonstrated reproducible gene-diet interactions (GDIs) impacting metabolic disease risk factors, likely due in part to measurement error in dietary intake estimation and insufficient capture of rare genetic variation. We aimed to identify GDIs across the genetic frequency spectrum impacting the macronutrient-glycemia relationship in genetically and culturally diverse cohorts. We analyzed 33,187 participants free of diabetes from 10 National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program cohorts with whole-genome sequencing, self-reported diet, and glycemic trait data. We fit cohort-specific, multivariable-adjusted linear mixed models for the effect of diet, modeled as an isocaloric substitution of carbohydrate for fat, and its interactions with common and rare variants genome-wide. In main effect meta-analyses, participants consuming more carbohydrate had modestly lower glycemic trait values (e.g., for glycated hemoglobin [HbA1c], -0.013% HbA1c/250 kcal substitution). In GDI meta-analyses, a common African ancestry-enriched variant (rs79762542) reached study-wide significance and replicated in the UK Biobank cohort, indicating a negative carbohydrate-HbA1c association among major allele homozygotes only. Simulations revealed that >150,000 samples may be necessary to identify similar macronutrient GDIs under realistic assumptions about effect size and measurement error. These results generate hypotheses for further exploration of modifiable metabolic disease risk in additional cohorts with African ancestry.

ARTICLE HIGHLIGHTS

We aimed to identify genetic modifiers of the dietary macronutrient-glycemia relationship using whole-genome sequence data from 10 Trans-Omics for Precision Medicine program cohorts. Substitution models indicated a modest reduction in glycemia associated with an increase in dietary carbohydrate at the expense of fat. Genome-wide interaction analysis identified one African ancestry-enriched variant near the FRAS1 gene that may interact with macronutrient intake to influence hemoglobin A1c. Simulation-based power calculations accounting for measurement error suggested that substantially larger sample sizes may be necessary to discover further gene-macronutrient interactions.

The Future of Social Determinants of Health: Looking Upstream to Structural Drivers

Policy Points Policies that redress oppressive social, economic, and political conditions are essential for improving population health and achieving health equity. Efforts to remedy structural oppression and its deleterious effects should account for its multilevel, multifaceted, interconnected, systemic, and intersectional nature. The U.S. Department of Health and Human Services should facilitate the creation and maintenance of a national publicly available, user-friendly data infrastructure on contextual measures of structural oppression. Publicly funded research on social determinants of health should be mandated to (a) analyze health inequities in relation to relevant data on structural conditions and (b) deposit the data in the publicly available data repository.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population

Background

Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins.

Methods

We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8).

Results

The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis.

Conclusion

Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.

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

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

Variance-quantitative trait loci enable systematic discovery of gene-environment interactions for cardiometabolic serum biomarkers

Gene-environment interactions represent the modification of genetic effects by environmental exposures and are critical for understanding disease and informing personalized medicine. These often induce differential phenotypic variance across genotypes; these variance-quantitative trait loci can be prioritized in a two-stage interaction detection strategy to greatly reduce the computational and statistical burden and enable testing of a broader range of exposures. We perform genome-wide variance-quantitative trait locus analysis for 20 serum cardiometabolic biomarkers by multi-ancestry meta-analysis of 350,016 unrelated participants in the UK Biobank, identifying 182 independent locus-biomarker pairs (p < 4.5×10-9). Most are concentrated in a small subset (4%) of loci with genome-wide significant main effects, and 44% replicate (p < 0.05) in the Women's Genome Health Study (N = 23,294). Next, we test each locus-biomarker pair for interaction across 2380 exposures, identifying 847 significant interactions (p < 2.4×10-7), of which 132 are independent (p < 0.05) after accounting for correlation between exposures. Specific examples demonstrate interaction of triglyceride-associated variants with distinct body mass- versus body fat-related exposures as well as genotype-specific associations between alcohol consumption and liver stress at the ADH1B gene. Our catalog of variance-quantitative trait loci and gene-environment interactions is publicly available in an online portal.

Gene–Environment Interaction on Type 2 Diabetes Risk among Chinese Adults Born in Early 1960s

Background: Gene–environment interactions on type 2 diabetes (T2D) risk are studied little among Chinese adults. Aim: This study aimed to explore the interactions among Chinese adults born in early 1960s. Methods: The interaction of single nucleotide polymorphisms (SNPs) and environmental factors on T2D risk were analyzed by multiple linear or logistic regression models, and in total 2216 subjects were included with the age of 49.7 ± 1.5 years. Results: High dietary intake increased the effects of rs340874 on impaired fasting glucose (IFG), rs5015480, rs7612463 on T2D (OR = 2.27, 2.37, 11.37, respectively), and reduced the effects of rs7172432 on IFG, rs459193 on impaired glucose tolerance (IGT) (OR = 0.08, 0.28, respectively). The associations between rs4607517 and T2D, rs10906115 and IGT, rs4607103, rs5015480 and IFG could be modified by drinking/smoking (OR = 2.28, 0.20, 3.27, 2.58, respectively). Physical activity (PA) interacted with rs12970134, rs2191349, rs4607517 on T2D (OR = 0.39, 3.50, 2.35, respectively), rs2796441 and rs4607517 on IGT (OR = 0.42, 0.33, respectively), and rs4430796, rs5215, and rs972283 on IFG (OR = 0.39, 3.05, 7.96, respectively). Significant interactions were identified between socioeconomic status and rs10830963, rs13266634 on T2D (OR = 0.41, 0.44, respectively), rs1470579 and rs2796441 on IGT (OR = 2.13, 2.37, respectively), and rs7202877 and rs7612463 on IFG (OR = 5.64, 9.18, respectively). Conclusion: There indeed existed interactions between environmental factors and genetic variants on T2D risk among Chinese adults.

Interactions between Polygenic Risk Scores, Dietary Pattern, and Menarche Age with the Obesity Risk in a Large Hospital-Based Cohort

Obese Asians are more susceptible to metabolic diseases than obese Caucasians of the same body mass index (BMI). We hypothesized that the genetic variants associated with obesity risk interact with the lifestyles of middle-aged and elderly adults, possibly allowing the development of personalized interventions based on genotype. We aimed to examine this hypothesis in a large city hospital-based cohort in Korea. The participants with cancers, thyroid diseases, chronic kidney disease, or brain-related diseases were excluded. The participants were divided into case and control according to their BMI: ≥25 kg/m2 (case; n = 17,545) and <25 kg/m2 (control; n = 36,283). The genetic variants that affected obesity risk were selected using a genome-wide association study, and the genetic variants that interacted with each other were identified by generalized multifactor dimensionality reduction analysis. The selected genetic variants were confirmed in the Ansan/Ansung cohort, and polygenetic risk scores (PRS)-nutrient interactions for obesity risk were determined. A high BMI was associated with a high-fat mass (odds ratio (OR) = 20.71) and a high skeletal muscle-mass index (OR = 3.38). A high BMI was positively related to metabolic syndrome and its components, including lipid profiles, whereas the initial menstruation age was inversely associated with a high BMI (OR = 0.78). The best model with 5-SNPs included SEC16B_rs543874, DNAJC27_rs713586, BDNF_rs6265, MC4R_rs6567160, and GIPR_rs1444988703. The high PRS with the 5-SNP model was positively associated with an obesity risk of 1.629 (1.475-1.798) after adjusting for the covariates. The 5-SNP model interacted with the initial menstruation age, fried foods, and plant-based diet for BMI risk. The participants with a high PRS also had a higher obesity risk when combined with early menarche, low plant-based diet, and a high fried-food intake than in participants with late menarche, high plant-based diet, and low fried-food intake. In conclusion, people with a high PRS and earlier menarche age are recommended to consume fewer fried foods and a more plant-based diet to decrease obesity risk. This result can be applied to personalized nutrition for preventing obesity.