Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1,2,3,4,5,6,7. This detailed knowledge of the genetic... Show moreGenome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1,2,3,4,5,6,7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8,9,10,11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases. Show less
Triglyceride (TG)-lowering LPL variants in combination with genetic LDL-C-lowering variants are associated with reduced risk of coronary artery disease (CAD). Genetic variation in the APOA5 gene... Show moreTriglyceride (TG)-lowering LPL variants in combination with genetic LDL-C-lowering variants are associated with reduced risk of coronary artery disease (CAD). Genetic variation in the APOA5 gene encoding apolipoprotein A-V also strongly affects TG levels, but the potential clinical impact and underlying mechanisms are yet to be resolved. Here, we aimed to study the effects of APOA5 genetic variation on CAD risk and plasma lipoproteins through factorial genetic association analyses. Using data from 309,780 European-ancestry participants from the UK Biobank, we evaluated the effects of lower TG levels as a result of genetic variation in APOA5 and/or LPL on CAD risk with or without a background of reduced LDL-C. Next, we compared lower TG levels via APOA5 and LPL variation with over 100 lipoprotein measurements in a combined sample from the Netherlands Epidemiology of Obesity study (N = 4,838) and the Oxford Biobank (N = 6,999). We found that lower TG levels due to combined APOA5 and LPL variation and genetically-influenced lower LDL-C levels afforded the largest reduction in CAD risk (odds ratio: 0.78 (0.73-0.82)). Compared to patients with genetically-influenced lower TG via LPL, genetically-influenced lower TG via APOA5 had similar and independent, but notably larger, effects on the lipoprotein profile. Our results suggest that lower TG levels as a result of APOA5 variation have strong beneficial effects on CAD risk and the lipoprotein profile, which suggest apo A-V may be a potential novel therapeutic target for CAD prevention. Show less
Triglyceride (TG)-lowering LPL variants in combination with genetic LDL-C-lowering variants are associated with reduced risk of coronary artery disease (CAD). Genetic variation in the APOA5 gene... Show moreTriglyceride (TG)-lowering LPL variants in combination with genetic LDL-C-lowering variants are associated with reduced risk of coronary artery disease (CAD). Genetic variation in the APOA5 gene encoding apolipoprotein A-V also strongly affects TG levels, but the potential clinical impact and underlying mechanisms are yet to be resolved. Here, we aimed to study the effects of APOA5 genetic variation on CAD risk and plasma lipoproteins through factorial genetic association analyses. Using data from 309,780 European-ancestry participants from the UK Biobank, we evaluated the effects of lower TG levels as a result of genetic variation in APOA5 and/or LPL on CAD risk with or without a background of reduced LDL-C. Next, we compared lower TG levels via APOA5 and LPL variation with over 100 lipoprotein measurements in a combined sample from the Netherlands Epidemiology of Obesity study (N = 4,838) and the Oxford Biobank (N = 6,999). We found that lower TG levels due to combined APOA5 and LPL variation and genetically-influenced lower LDL-C levels afforded the largest reduction in CAD risk (odds ratio: 0.78 (0.73–0.82)). Compared to patients with genetically-influenced lower TG via LPL, genetically-influenced lower TG via APOA5 had similar and independent, but notably larger, effects on the lipoprotein profile. Our results suggest that lower TG levels as a result of APOA5 variation have strong beneficial effects on CAD risk and the lipoprotein profile, which suggest apo A-V may be a potential novel therapeutic target for CAD prevention. Show less
Background Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed... Show moreBackground Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. Methods Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. Results Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99-1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96-1.04; p value 0.59). Conclusions Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk. Show less
Background and aims: Mendelian randomization studies have shown that triglyceride (TG)- lowering lipoprotein lipase (LPL) alleles and low-density lipoprotein-cholesterol (LDL-C)-lowering alleles... Show moreBackground and aims: Mendelian randomization studies have shown that triglyceride (TG)- lowering lipoprotein lipase (LPL) alleles and low-density lipoprotein-cholesterol (LDL-C)-lowering alleles have independent beneficial associations on cardiovascular disease (CVD) risk. We aimed to provide further insight into this observation by applying Mendelian randomization analyses of genetically-influenced TG and LDL-C levels on plasma metabolomic profiles Methods: We quantified over 100 lipoprotein metabolomic measures in the Netherlands Epidemiology of Obesity (NEO) study (N = 4838) and Oxford Biobank (OBB) (N = 6999) by nuclear magnetic resonance (NMR) spectroscopy. Weighted genetic scores for TG via five LPL alleles and LDL-C via 19 alleles were calculated and dichotomized by the median, resulting in four genotype combinations of high/low TG and high/low LDL-C. We performed linear regression analyses using a two & times; two design with the group with genetically-influenced high TG and LDL-C as a reference. Results: Compared to the individual groups with genetically-influenced lower TG or lower LDL-C only, the group with combined genetically-influenced lower TG and LDL-C showed an overall independent and additive pattern of changes in metabolomic measures. Over 100 measures were different (p < 1.35 & times; 10-3) compared to the reference, with effect sizes and directionality being similar in NEO and OBB. Most notably, levels of all very-low density lipoprotein (VLDL) and LDL sub-particles were lower. Conclusions: Our findings provide evidence that TG-lowering on top of LDL-C-lowering has additive beneficial effects on the lipoprotein profile compared to TG-lowering or LDL-C-lowering only, which is in accordance with reported additive genetic effects on CVD risk reduction. Show less
Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here... Show moreGlycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 x 10(-8)), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.A trans-ancestry meta-analysis of GWAS of glycemic traits in up to 281,416 individuals identifies 99 novel loci, of which one quarter was found due to the multi-ancestry approach, which also improves fine-mapping of credible variant sets. Show less
Bos, M.M.; Goulding, N.J.; Lee, M.A.; Hofman, A.; Bot, M.; Pool, R.; ... ; Lawlor, D.A. 2021
Background Sleep traits are associated with cardiometabolic disease risk, with evidence from Mendelian randomization (MR) suggesting that insomnia symptoms and shorter sleep duration increase... Show moreBackground Sleep traits are associated with cardiometabolic disease risk, with evidence from Mendelian randomization (MR) suggesting that insomnia symptoms and shorter sleep duration increase coronary artery disease risk. We combined adjusted multivariable regression (AMV) and MR analyses of phenotypes of unfavourable sleep on 113 metabolomic traits to investigate possible biochemical mechanisms linking sleep to cardiovascular disease. Methods We used AMV (N = 17,368) combined with two-sample MR (N = 38,618) to examine effects of self-reported insomnia symptoms, total habitual sleep duration, and chronotype on 113 metabolomic traits. The AMV analyses were conducted on data from 10 cohorts of mostly Europeans, adjusted for age, sex, and body mass index. For the MR analyses, we used summary results from published European-ancestry genome-wide association studies of self-reported sleep traits and of nuclear magnetic resonance (NMR) serum metabolites. We used the inverse-variance weighted (IVW) method and complemented this with sensitivity analyses to assess MR assumptions. Results We found consistent evidence from AMV and MR analyses for associations of usual vs. sometimes/rare/never insomnia symptoms with lower citrate (- 0.08 standard deviation (SD)[95% confidence interval (CI) - 0.12, - 0.03] in AMV and - 0.03SD [- 0.07, - 0.003] in MR), higher glycoprotein acetyls (0.08SD [95% CI 0.03, 0.12] in AMV and 0.06SD [0.03, 0.10) in MR]), lower total very large HDL particles (- 0.04SD [- 0.08, 0.00] in AMV and - 0.05SD [- 0.09, - 0.02] in MR), and lower phospholipids in very large HDL particles (- 0.04SD [- 0.08, 0.002] in AMV and - 0.05SD [- 0.08, - 0.02] in MR). Longer total sleep duration associated with higher creatinine concentrations using both methods (0.02SD per 1 h [0.01, 0.03] in AMV and 0.15SD [0.02, 0.29] in MR) and with isoleucine in MR analyses (0.22SD [0.08, 0.35]). No consistent evidence was observed for effects of chronotype on metabolomic measures. Conclusions Whilst our results suggested that unfavourable sleep traits may not cause widespread metabolic disruption, some notable effects were observed. The evidence for possible effects of insomnia symptoms on glycoprotein acetyls and citrate and longer total sleep duration on creatinine and isoleucine might explain some of the effects, found in MR analyses of these sleep traits on coronary heart disease, which warrant further investigation. Show less
Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to similar to 1.3 million participants, we... Show moreGenetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to similar to 1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency <= 0.01) variant BP associations (P < 5 x 10(-8)), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were similar to 8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets. Show less
Smoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been... Show moreSmoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been associated with smoking behaviour-related traits. We tested up to 235,116 single nucleotide variants (SNVs) on the exome-array for association with smoking initiation, cigarettes per day, pack-years, and smoking cessation in a fixed effects meta-analysis of up to 61 studies (up to 346,813 participants). In a subset of 112,811 participants, a further one million SNVs were also genotyped and tested for association with the four smoking behaviour traits. SNV-trait associations withP < 5 x 10(-8)in either analysis were taken forward for replication in up to 275,596 independent participants from UK Biobank. Lastly, a meta-analysis of the discovery and replication studies was performed. Sixteen SNVs were associated with at least one of the smoking behaviour traits (P < 5 x 10(-8)) in the discovery samples. Ten novel SNVs, including rs12616219 nearTMEM182, were followed-up and five of them (rs462779 inREV3L, rs12780116 inCNNM2, rs1190736 inGPR101, rs11539157 inPJA1, and rs12616219 nearTMEM182) replicated at a Bonferroni significance threshold (P < 4.5 x 10(-3)) with consistent direction of effect. A further 35 SNVs were associated with smoking behaviour traits in the discovery plus replication meta-analysis (up to 622,409 participants) including a rare SNV, rs150493199, inCCDC141and two low-frequency SNVs inCEP350andHDGFRP2. Functional follow-up implied that decreased expression ofREV3Lmay lower the probability of smoking initiation. The novel loci will facilitate understanding the genetic aetiology of smoking behaviour and may lead to the identification of potential drug targets for smoking prevention and/or cessation. Show less
Aims/hypothesis Circulating succinate and 12,13-dihydroxy-9Z-octadecenoic acid (12,13-diHOME) were recently shown to promote brown adipocyte thermogenesis and protect against metabolic disorders in... Show moreAims/hypothesis Circulating succinate and 12,13-dihydroxy-9Z-octadecenoic acid (12,13-diHOME) were recently shown to promote brown adipocyte thermogenesis and protect against metabolic disorders in rodents. This study aimed to evaluate the associations between plasma levels of these metabolites and adiposity and metabolic profile in humans. Methods Fasting plasma succinate and 12,13-diHOME levels were quantified using ultra HPLC-tandem MS in 2248 individuals (50% female, mean age 41.3 +/- 5.9 years, mean BMI 26.1 +/- 4.6 kg/m(2)) in addition to fasting plasma biochemistry. Total and regional adiposity were assessed with dual-energy x-ray absorptiometry. An age- and sex-adjusted linear regression model was used to determine the associations between succinate and 12,13-diHOME levels and body composition and metabolic profile. Two-sample Mendelian randomisation was used to assess the associations between genetically determined BMI and metabolic traits with circulating plasma succinate and 12,13-diHOME. Results A one-SD higher plasma succinate and 12,13-diHOME concentration was associated with a 0.15 SD (95% CI 0.28, 0.03) and 0.08 SD (0.15, 0.01) lower total fat mass respectively. Additionally, a one-SD higher plasma 12,13-diHOME level was associated with a 0.09 SD (0.16, 0.02) lower fasting plasma insulin and 0.10 SD (0.17, 0.04) lower plasma triacylglycerol. In Mendelian randomisation analyses, genetically determined higher BMI, fasting hyperinsulinaemia and elevated lipid levels were not associated with changes in either plasma succinate or plasma 12,13-diHOME concentrations. No indications of bias due to directional pleiotropy were detected in the Mendelian randomisation analyses. Conclusions/interpretation Our findings tentatively suggest that plasma succinate and 12,13-diHOME may play a role in the regulation of energy metabolism and brown adipose tissue activation in humans. Further studies encompassing direct assessment of brown adipose tissue activity and dietary supplementation are necessary to investigate the potential beneficial effects of these metabolites on systemic metabolism. Show less
In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear,... Show moreIn many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (F-ROH) for >1.4 million individuals, we show that F-ROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: F-ROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of F-ROH are confirmed within full-sibling pairs, where the variation in F-ROH is independent of all environmental confounding. Show less
Context: Seasonal variation in cold and light exposure may influence metabolic health.Objective: We assessed the associations of bright sunlight and outdoor temperature with measures of glucose and... Show moreContext: Seasonal variation in cold and light exposure may influence metabolic health.Objective: We assessed the associations of bright sunlight and outdoor temperature with measures of glucose and lipid metabolism in two populations of middle-aged European subjects.Design: Cross-sectional study.Setting: Two population-based European cohorts.Participants: Middle-aged nondiabetic subjects from the Oxford Biobank (OBB; N = 4327; mean age, 41.4 years) and the Netherlands Epidemiology of Obesity (NEO) study (N = 5899; mean age, 55.6 years).Interventions: Data on outdoor bright sunlight and temperature collected from local weather stations.Main Outcome Measures: Insulin resistance and fasting lipid levels. Multivariable regression analyses were adjusted for age, sex, percentage body fat, season, and either outdoor temperature or bright sunlight.Results: In the OBB cohort, increased bright sunlight exposure was associated with lower fasting insulin [-1.27% (95% CI, -2.09 to -0.47%) per extra hour of bright sunlight], lower homeostatic model assessment for insulin resistance (-1.36%; 95% CI, -2.23 to -0.50), lower homeostatic model assessment for beta-cell function (-0.80%; 95% CI, -1.31 to -0.30), and lower triglyceride (-1.28%; 95% CI, -2.07 to -0.50) levels. In the NEO cohort generally unidirectional but weaker associations were observed. No associations between outdoor temperature and measures of glucose or lipid metabolism were detected following adjustment for bright sunlight.Conclusions: Bright sunlight, but not outdoor temperature, might be associated with increased insulin sensitivity and lower triglyceride levels. Show less
Liu, D.J.; Peloso, G.M.; Yu, H.; Butterworth, A.S.; Wang, X.; Mahajan, A.; ... ; VA Million Veteran Program 2017
Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this... Show moreWaist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 x 10(-9) to P = 1.8 x 10(-40)) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 x 10(-3) to P = 1.2 x 10(-13)). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions. Show less
Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between... Show moreObesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and similar to 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 x 10(-8)), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less
Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting... Show moreLevels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less
