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
Per- and polyfluoroalkyl substances (PFAS) are widely used and persistent chemicals, leading to ubiquitous exposure. Although high PFAS levels have been associated with an adverse cardiovascular... Show morePer- and polyfluoroalkyl substances (PFAS) are widely used and persistent chemicals, leading to ubiquitous exposure. Although high PFAS levels have been associated with an adverse cardiovascular risk profile, the distribution of levels and relations with cardio-metabolic risk markers in the general population have not been fully characterized. We assessed the association between blood levels of perfluorooctaneic acid (PFOA), perfluorooctane sulfonic acid (PFOS), and perfluorohexanesulfonic acid (PFHxS) and a range of lipoproteins and metabolites as well as clinical lipid measurements. We used data from participants of the Netherlands Epidemiology of Obesity study (NEO) (n = 584) and the Rhineland Study (n = 1962), jointly spanning an age range of 30 to 89 years. PFAS were measured with the Metabolon HD4 platform, and lipoprotein and metabolite profiles were measured using Nightingale's nuclear magnetic resonance-spectroscopy platform, and mainly comprised lipoprotein markers. Using linear regression analyses, we quantified age-, sex-, and education-adjusted associations of PFOA, PFOS, and PFHxS with clinical lipid measurements and 224 lipoproteins and metabolites. Higher levels of PFAS, particularly PFOS and PFHxS, were associated with higher concentrations of total lipid, cholesterol and phospholipid content in most HDL, IDL, LDL, and VLDL subclasses. The effect sizes were age-dependent for the majority of the associations, with the deleterious effects of PFAS being generally stronger in people below compared to those above median age. Our observation that in the general population even low PFAS concentrations are associated with an unfavorable lipid profile, calls for further critical regulation of PFAS substances. Show less
Yuan, L.S.; Li-Gao, R.; Verhoeven, A.; Eyk, H.J. van; Bizino, M.B.; Rensen, P.C.N.; ... ; Berg, B.M. van den 2023
Background: Composition of high-density lipoproteins (HDL) is emerging as an important determinant in the development of microvascular complications in type 2 diabetes mellitus (T2DM). Dutch South... Show moreBackground: Composition of high-density lipoproteins (HDL) is emerging as an important determinant in the development of microvascular complications in type 2 diabetes mellitus (T2DM). Dutch South Asian (DSA) individuals with T2DM display an increased risk of microvascular complications compared with Dutch white Caucasian (DwC) individuals with T2DM. In this study, we aimed to investigate whether changes in HDL composition associate with increased microvascular risk in this ethnic group and lead to new lipoprotein biomarkers. Materials and Methods: Using H-1 nuclear magnetic resonance spectroscopy and Bruker IVDr Lipoprotein Subclass Analysis (B.I.LISA) software, plasma lipoprotein changes were determined in 51 healthy individuals (30 DwC, 21 DSA) and 92 individuals with T2DM (45 DwC, 47 DSA) in a cross-sectional, case-control study. Differential HDL subfractions were investigated using multinomial logistic regression analyses, adjusting for possible confounders including BMI and diabetes duration. Results: We identified HDL compositional differences between healthy and diabetic individuals in both ethnic groups. Specifically, levels of apolipoprotein A2 and HDL-4 subfractions were lower in DSA compared with DwC with T2DM. Apolipoprotein A2 and HDL-4 subfractions also negatively correlated with waist circumference, waistto-hip ratio, haemoglobin A1c, glucose levels and disease duration in DSA with T2DM, and associated with increased incidence of microvascular complications .Conclusion: While HDL composition differed between controls and T2DM in both ethnic groups, the lower levels of lipid content in the smallest HDL subclass (HDL-4) in DSA with T2DM appeared to be more clinically relevant, with higher odds of having diabetes-related pan-microvascular complications such as retinopathy and neuropathy. These typical differences in HDL could be used as ethnicity-specific T2DM biomarkers. Show less
Background and Aims Non-alcoholic fatty liver disease (NAFLD) is characterized by the pathological accumulation of triglycerides in hepatocytes and is associated with insulin resistance,... Show moreBackground and Aims Non-alcoholic fatty liver disease (NAFLD) is characterized by the pathological accumulation of triglycerides in hepatocytes and is associated with insulin resistance, atherogenic dyslipidaemia and cardiometabolic diseases. Thus far, the extent of metabolic dysregulation associated with hepatic triglyceride accumulation has not been fully addressed. In this study, we aimed to identify metabolites associated with hepatic triglyceride content (HTGC) and map these associations using network analysis. Methods: To gain insight in the spectrum of metabolites associated with hepatic triglyceride accumulation, we performed a comprehensive plasma metabolomics screening of 1363 metabolites in apparently healthy middle aged (age 45-65) individuals (N = 496) in whom HTGC was measured by proton magnetic resonance spectroscopy. An atlas of metabolite-HTGC associations, based on univariate results, was created using correlation-based Gaussian graphical model (GGM) and genome scale metabolic model network analyses. Pathways associated with the clinical prognosis marker fibrosis 4 (FIB-4) index were tested using a closed global test. Results: Our analyses revealed that 118 metabolites were univariately associated with HTGC (p-value <6.59 x 10(-5)), including 106 endogenous, 1 xenobiotic and 11 partially characterized/uncharacterized metabolites. These associations were mapped to several biological pathways including branched amino acids (BCAA), diglycerols, sphingomyelin, glucosyl-ceramide and lactosyl-ceramide. We also identified a novel possible HTGC-related pathway connecting glutamate, metabolonic lactone sulphate and X-15245 using the GGM network. These pathways were confirmed to be associated with the FIB-4 index as well. The full interactive metabolite-HTGC atlas is provided online: . Conclusions: The combined network and pathway analyses indicated extensive associations between BCAA and the lipids pathways with HTGC and the FIB-4 index. Moreover, we report a novel pathway glutamate-metabolonic lactone sulphate-X-15245 with a potential strong association with HTGC. These findings can aid elucidating HTGC metabolomic profiles and provide insight into novel drug targets for fibrosis-related outcomes. Show less
Tandem cytosine-adenine-guanine (CAG) repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington's disease (HD). Apart from neuropsychiatric complications, the disease is also... Show moreTandem cytosine-adenine-guanine (CAG) repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington's disease (HD). Apart from neuropsychiatric complications, the disease is also accompanied by metabolic dysregulation and weight loss, which contribute to a progressive functional decline. Recent studies also reported an association between repeats below the pathogenic threshold (<36) for HD and body mass index (BMI), suggesting that HTT repeat sizes in the non-pathogenic range are associated with metabolic dysregulation. In this study, we hypothesized that HTT repeat sizes < 36 are associated with metabolite levels, possibly mediated through reduced BMI. We pooled data from three European cohorts (n = 10 228) with genotyped HTT CAG repeat size and metabolomic measurements. All 145 metabolites were measured on the same targeted platform in all studies. Multilevel mixed-effects analysis using the CAG repeat size in HTT identified 67 repeat size metabolite associations. Overall, the metabolomic profile associated with larger CAG repeat sizes in HTT were unfavorable-similar to those of higher risk of coronary artery disease and type 2 diabetes-and included elevated levels of amino acids, fatty acids, low-density lipoprotein (LDL)-, very low-density lipoprotein- and intermediate density lipoprotein (IDL)-related metabolites while with decreased levels of very large high-density lipoprotein (HDL)-related metabolites. Furthermore, the associations of 50 metabolites, in particular, specific very large HDL-related metabolites, were mediated by lower BMI. However, no mediation effect was found for 17 metabolites related to LDL and IDL. In conclusion, our findings indicate that large non-pathogenic CAG repeat sizes in HTT are associated with an unfavorable metabolomic profile despite their association with a lower BMI. Show less
Li-Gao, R.; Hughes, D.A.; Klinken, J.B. van; Mutsert, R. de; Rosendaal, F.R.; Mook-Kanamori, D.O.; ... ; Dijk, K.W. van 2021
Humans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation... Show moreHumans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation in concentrations of postprandial metabolites (t = 150 min) in response to a liquid mixed meal through genome-wide association studies (GWAS) performed in the Netherlands Epidemiology of Obesity (NEO) study (n = 5,705). The metabolite response GWAS identified an association between glucose change and rs10830963:G in the melatonin receptor 1B (beta [SE] -0.23 [0.03], P = 2.15 x 10(-19)). In addition, the ANKRD55 locus led by rs458741:C showed strong associations with extremely large VLDL (XXLVLDL) particle response (XXLVLDL total cholesterol: beta [SE] 0.17 [0.03], P = 5.76 x 10(-10); XXLVLDL cholesterol ester: beta [SE] 0.17 [0.03], P = 9.74 x 10(-10)), which also revealed strong associations with body composition and diabetes in the UK Biobank (P < 5 x 10(-8)). Furthermore, the associations between XXLVLDL response and insulinogenic index, HOMA-beta, Matsuda insulin sensitivity index, and HbA(1c) in the NEO study implied the role of chylomicron synthesis in diabetes (with false discovery rate-corrected q <0.05). To conclude, genetic studies of metabolomics change after a liquid meal illuminate novel pathways for glucose and lipid metabolism. Further studies are warranted to corroborate biological pathways of the ANKRD55 locus underlying diabetes. 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
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
Fuentes, L. de las; Sung, Y.J.; Noordam, R.; Winkler, T.; Feitosa, M.F.; Schwander, K.; ... ; Lifelines Cohort Study 2020
Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed... Show moreEducational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 x 10(-8)). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP. 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