Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-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-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.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

No evidence linking sleep traits with white blood cell counts: Multivariable-adjusted and Mendelian randomization analyses.

BACKGROUND: Disturbances in habitual sleep have been associated with multiple age-associated diseases. However, the biological mechanisms underpinning these associations remain largely unclear. We assessed the possible involvement of the circulating immune system by determining the associations between sleep traits and white blood cell counts using multivariable-adjusted linear regression and Mendelian randomization. METHODS: Cross-sectional multivariable-adjusted linear regression analyses were done using participants within the normal range of total white blood cell counts (>4.5 × 109 and <11.0 × 109/µL) from UK Biobank. For the sleep traits, we examined (short and long) sleep duration, chronotype, insomnia symptoms and daytime dozing. Two-sample Mendelian randomization analyses were done using instruments for sleep traits derived from European-ancestry participants from UK Biobank (over 410,000 participants) and using SNP-outcome data derived from European-ancestry participants from the Blood Cell Consortium (N = 563,946) to which no data from UK Biobank contributed. RESULTS: Using data from 357,656 participants (mean [standard deviation] age: 56.5 [8.1] years, and 44.4% men), we did not find evidence that disturbances in any of the studied sleep traits were associated with differences in blood cell counts (total, lymphocytes, neutrophiles, eosinophiles and basophiles). Also, we did not find associations between disturbances in any of the studied sleep traits and white blood cell counts using Mendelian Randomization. CONCLUSION: Based on the results from two different methodologies, disturbances in habitual sleep are unlikely to cause changes in blood cell counts and thereby differences in blood cell counts are unlikely to be underlying the observed sleep-disease associations.

PLIS: A metabolomic response monitor to a lifestyle intervention study in older adults.

The response to lifestyle intervention studies is often heterogeneous, especially in older adults. Subtle responses that may represent a health gain for individuals are not always detected by classical health variables, stressing the need for novel biomarkers that detect intermediate changes in metabolic, inflammatory, and immunity-related health. Here, our aim was to develop and validate a molecular multivariate biomarker maximally sensitive to the individual effect of a lifestyle intervention; the Personalized Lifestyle Intervention Status (PLIS). We used 1 H-NMR fasting blood metabolite measurements from before and after the 13-week combined physical and nutritional Growing Old TOgether (GOTO) lifestyle intervention study in combination with a fivefold cross-validation and a bootstrapping method to train a separate PLIS score for men and women. The PLIS scores consisted of 14 and four metabolites for females and males, respectively. Performance of the PLIS score in tracking health gain was illustrated by association of the sex-specific PLIS scores with several classical metabolic health markers, such as BMI, trunk fat%, fasting HDL cholesterol, and fasting insulin, the primary outcome of the GOTO study. We also showed that the baseline PLIS score indicated which participants respond positively to the intervention. Finally, we explored PLIS in an independent physical activity lifestyle intervention study, showing similar, albeit remarkably weaker, associations of PLIS with classical metabolic health markers. To conclude, we found that the sex-specific PLIS score was able to track the individual short-term metabolic health gain of the GOTO lifestyle intervention study. The methodology used to train the PLIS score potentially provides a useful instrument to track personal responses and predict the participant's health benefit in lifestyle interventions similar to the GOTO study.

Hepatic triglyceride content is intricately associated with numerous metabolites and biochemical pathways.

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, 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 × 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: https://tofaquih.github.io/AtlasLiver/. 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.

A novel approach for pharmacological substantiation of safety signals using plasma concentrations of medication and administrative/healthcare databases: A case study using Danish registries for an FDA warning on lamotrigine.

PHARMACOM-EPI is a novel framework to predict plasma concentrations of drugs at the time of occurrence of clinical outcomes. In early 2021, the U.S. Food and Drug Administration (FDA) issued a warning on the antiseizure drug lamotrigine claiming that it has the potential to increase the risk of arrhythmias and related sudden cardiac death due to a pharmacological sodium channel-blocking effect. We hypothesized that the risk of arrhythmias and related death is due to toxicity. We used the PHARMACOM-EPI framework to investigate the relationship between lamotrigine's plasma concentrations and the risk of death in older patients using real-world data. Danish nationwide administrative and healthcare registers were used as data sources and individuals aged 65 years or older during the period 1996 - 2018 were included in the study. According to the PHARMACOM-EPI framework, plasma concentrations of lamotrigine were predicted at the time of death and patients were categorized into non-toxic and toxic groups based on the therapeutic range of lamotrigine (3-15 mg/L). Over 1 year of treatment, the incidence rate ratio (IRR) of all-cause mortality was calculated between the propensities score matched toxic and non-toxic groups. In total, 7286 individuals were diagnosed with epilepsy and were exposed to lamotrigine, 432 of which had at least one plasma concentration measurement The pharmacometric model by Chavez et al. was used to predict lamotrigine's plasma concentrations considering the lowest absolute percentage error among identified models (14.25 %, 95 % CI: 11.68-16.23). The majority of lamotrigine associated deaths were cardiovascular-related and occurred among individuals with plasma concentrations in the toxic range. The IRR of mortality between the toxic group and non-toxic group was 3.37 [95 % CI: 1.44-8.32] and the cumulative incidence of all-cause mortality exponentially increased in the toxic range. Application of our novel framework PHARMACOM-EPI provided strong evidence to support our hypothesis that the increased risk of all-cause and cardiovascular death was associated with a toxic plasma concentration level of lamotrigine among older lamotrigine users.

Delirium in older patients with COVID-19: Prevalence, risk factors and clinical outcomes across the first three waves of the pandemic.

OBJECTIVES: Delirium is a serious condition, which poses treatment challenges during hospitalisation for COVID-19. Improvements in testing, vaccination and treatment might have changed patient characteristics and outcomes through the pandemic. We evaluated whether the prevalence and risk factors for delirium, and the association of delirium with in-hospital mortality changed through the pandemic. METHODS: This study was part of the COVID-OLD study in 19 Dutch hospitals including patients ≥70 years in the first (spring 2020), second (autumn 2020) and third wave (autumn 2021). Multivariable logistic regression models were used to study risk factors for delirium, and in-hospital mortality. Differences in effect sizes between waves were studied by including interaction terms between wave and risk factor in logistic regression models. RESULTS: 1540, 884 and 370 patients were included in the first, second and third wave, respectively. Prevalence of delirium in the third wave (12.7%) was significantly lower compared to the first (22.5%) and second wave (23.5%). In multivariable-adjusted analyses, pre-existing memory problems was a consistent risk factor for delirium across waves. Previous delirium was a risk factor for delirium in the first wave (OR 4.02), but not in the second (OR 1.61) and third wave (OR 2.59, p-value interaction-term 0.028). In multivariable-adjusted analyses, delirium was not associated with in-hospital mortality in all waves. CONCLUSION: Delirium prevalence declined in the third wave, which might be the result of vaccination and improved treatment strategies. Risk factors for delirium remained consistent across waves, although some attenuation was seen in the second wave.

Weight loss, visit-to-visit body weight variability and cognitive function in older individuals.

OBJECTIVE: to investigate the association between variability and loss of body weight with subsequent cognitive performance and activities of daily living in older individuals. DESIGN: cross-sectional cohort study. SETTING: PROspective Study of Pravastatin in the Elderly at Risk, multicentre trial with participants from Scotland, Ireland and the Netherlands. SUBJECTS: 4,309 participants without severe cognitive dysfunction (mean age 75.1 years, standard deviation (SD) = 3.3), at higher risk for cardiovascular disease (CVD). METHODS: body weight was measured every 3 months for 2.5 years. Weight loss was defined as an average slope across all weight measurements and as ≥5% decrease in baseline body weight during follow-up. Visit-to-visit variability was defined as the SD of weight measurements (kg) between visits. Four tests of cognitive function were examined: Stroop test, letter-digit coding test (LDCT), immediate and delayed picture-word learning tests. Two measures of daily living activities: Barthel Index (BI) and instrumental activities of daily living (IADL). All tests were examined at month 30. RESULTS: both larger body weight variability and loss of ≥5% of baseline weight were independently associated with worse scores on all cognitive tests, but minimally with BI and IADL. Compared with participants with stable weight, participants with significant weight loss performed 5.83 seconds (95% CI 3.74; 7.92) slower on the Stroop test, coded 1.72 digits less (95% CI -2.21; -1.13) on the LDCT and remembered 0.71 pictures less (95% CI -0.93; -0.48) on the delayed picture-word learning test. CONCLUSION: in older people at higher risk for CVD, weight loss and variability are independent risk-factors for worse cognitive function.

Apolipoprotein A-V is a potential target for treating coronary artery disease: evidence from genetic and metabolomic analyses.

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 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.

Assessment of the bi-directional relationship between blood mitochondrial DNA copy number and type 2 diabetes mellitus: a multivariable-adjusted regression and Mendelian randomisation study.

AIMS/HYPOTHESIS: Mitochondrial dysfunction, which can be approximated by blood mitochondrial DNA copy number (mtDNA-CN), has been implicated in the pathogenesis of type 2 diabetes mellitus. Thus far, however, insights from prospective cohort studies and Mendelian randomisation (MR) analyses on this relationship are limited. We assessed the association between blood mtDNA-CN and incident type 2 diabetes using multivariable-adjusted regression analyses, and the associations between blood mtDNA-CN and type 2 diabetes and BMI using bi-directional MR. METHODS: Multivariable-adjusted Cox proportional hazard models were used to estimate the association between blood mtDNA-CN and incident type 2 diabetes in 285,967 unrelated European individuals from UK Biobank free of type 2 diabetes at baseline. Additionally, a cross-sectional analysis was performed to investigate the association between blood mtDNA-CN and BMI. We also assessed the potentially causal relationship between blood mtDNA-CN and type 2 diabetes (N=898,130 from DIAGRAM, N=215,654 from FinnGen) and BMI (N=681,275 from GIANT) using bi-directional two-sample MR. RESULTS: During a median follow-up of 11.87 years, 15,111 participants developed type 2 diabetes. Participants with a higher level of blood mtDNA-CN are at lower risk of developing type 2 diabetes (HR 0.90 [95% CI 0.89, 0.92]). After additional adjustment for BMI and other confounders, these results attenuated moderately and remained present. The multivariable-adjusted cross-sectional analyses showed that higher blood mtDNA-CN was associated with lower BMI (-0.12 [95% CI -0.14, -0.10]) kg/m2. In the bi-directional MR analyses, we found no evidence for causal associations between blood mtDNA-CN and type 2 diabetes, and blood mtDNA-CN and BMI in either direction. CONCLUSIONS/INTERPRETATION: The results from the present study indicate that the observed association between low blood mtDNA-CN and higher risk of type 2 diabetes is likely not causal.

Timing of objectively-collected physical activity in relation to body weight and metabolic health in sedentary older people: a cross-sectional and prospective analysis.

BACKGROUND: Little is known about the impact of timing as opposed to frequency and intensity of daily physical activity on metabolic health. Therefore, we assessed the association between accelerometery-based daily timing of physical activity and measures of metabolic health in sedentary older people. METHODS: Hourly mean physical activity derived from wrist-worn accelerometers over a 6-day period was collected at baseline and after 3 months in sedentary participants from the Active and Healthy Ageing study. A principal component analysis (PCA) was performed to reduce the number of dimensions (e.g. define periods instead of separate hours) of hourly physical activity at baseline and change during follow-up. Cross-sectionally, a multivariable-adjusted linear regression analysis was used to associate the principal components, particularly correlated with increased physical activity in data-driven periods during the day, with body mass index (BMI), fasting glucose and insulin, HbA1c and the homeostatic model assessment for insulin resistance (HOMA-IR). For the longitudinal analyses, we calculated the hourly changes in physical activity and change in metabolic health after follow-up. RESULTS: We included 207 individuals (61.4% male, mean age: 64.8 [SD 2.9], mean BMI: 28.9 [4.7]). Higher physical activity in the early morning was associated with lower fasting glucose (-2.22%, 95% CI: -4.19, -0.40), fasting insulin (-13.54%, 95%CI: -23.49, -4.39), and HOMA-IR (-16.07%, 95%CI: -27.63, -5.65). Higher physical activity in the late afternoon to evening was associated with lower BMI (-2.84%, 95% CI: -4.92, -0.70). Higher physical activity at night was associated with higher BMI (2.86%, 95% CI: 0.90, 4.78), fasting glucose (2.57%, 95% CI: 0.70, 4.30), and HbA1c (2.37%, 95% CI: 1.00, 3.82). Similar results were present in the prospective analysis. CONCLUSION: Specific physical activity timing patterns were associated with more beneficial metabolic health, suggesting particular time-dependent physical activity interventions might maximise health benefits.

Genetic Susceptibility Loci in Genomewide Association Study of Cluster Headache.

OBJECTIVE: Identifying common genetic variants that confer genetic risk for cluster headache. METHODS: We conducted a case-control study in the Dutch Leiden University Cluster headache neuro-Analysis program (LUCA) study population (n = 840) and unselected controls from the Netherlands Epidemiology of Obesity Study (NEO; n = 1,457). Replication was performed in a Norwegian sample of 144 cases from the Trondheim Cluster headache sample and 1,800 controls from the Nord-Trøndelag Health Survey (HUNT). Gene set and tissue enrichment analyses, blood cell-derived RNA-sequencing of genes around the risk loci and linkage disequilibrium score regression were part of the downstream analyses. RESULTS: An association was found with cluster headache for 4 independent loci (r2 < 0.1) with genomewide significance (p < 5 × 10-8 ), rs11579212 (odds ratio [OR] = 1.51, 95% confidence interval [CI] = 1.33-1.72 near RP11-815 M8.1), rs6541998 (OR = 1.53, 95% CI = 1.37-1.74 near MERTK), rs10184573 (OR = 1.43, 95% CI = 1.26-1.61 near AC093590.1), and rs2499799 (OR = 0.62, 95% CI = 0.54-0.73 near UFL1/FHL5), collectively explaining 7.2% of the variance of cluster headache. SNPs rs11579212, rs10184573, and rs976357, as proxy SNP for rs2499799 (r2  = 1.0), replicated in the Norwegian sample (p < 0.05). Gene-based mapping yielded ASZ1 as possible fifth locus. RNA-sequencing indicated differential expression of POLR1B and TMEM87B in cluster headache patients. INTERPRETATION: This genomewide association study (GWAS) identified and replicated genetic risk loci for cluster headache with effect sizes larger than those typically seen in complex genetic disorders. ANN NEUROL 2021;90:203-216.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure.

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

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 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 × 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.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure.

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling.

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 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.

Investigating the relationships between unfavourable habitual sleep and metabolomic traits: evidence from multi-cohort multivariable regression and Mendelian randomization analyses.

BACKGROUND: 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.

Differential insulin sensitivity of NMR-based metabolomic measures in a two-step hyperinsulinemic euglycemic clamp study.

BACKGROUND: Insulin is the key regulator of glucose metabolism, but it is difficult to dissect direct insulin from glucose-induced effects. We aimed to investigate the effects of hyperinsulemia on metabolomic measures under euglycemic conditions in nondiabetic participants. METHODS: We assessed concentrations of 151 metabolomic measures throughout a two-step hyperinsulinemic euglycemic clamp procedure. We included 24 participants (50% women, mean age = 62 [s.d. = 4.2] years) and metabolomic measures were assessed under baseline, low-dose (10 mU/m2/min) and high-dose (40 mU/m2/min) insulin conditions. The effects of low- and high-dose insulin infusion on metabolomic measures were analyzed using linear mixed-effect models for repeated measures. RESULTS: After low-dose insulin infusion, 90 metabolomic measures changed in concentration (p < 1.34e-4), among which glycerol (beta [Confidence Interval] = - 1.41 [- 1.54, - 1.27] s.d., p = 1.28e-95) and three-hydroxybutyrate (- 1.22 [- 1.36, - 1.07] s.d., p = 1.44e-61) showed largest effect sizes. After high-dose insulin infusion, 121 metabolomic measures changed in concentration, among which branched-chain amino acids showed the largest additional decrease compared with low-dose insulin infusion (e.g., Leucine, - 1.78 [- 1.88, - 1.69] s.d., P = 2.7e-295). More specifically, after low- and high-dose insulin infusion, the distribution of the lipoproteins shifted towards more LDL-sized particles with decreased mean diameters. CONCLUSION: Metabolomic measures are differentially insulin sensitive and may thus be differentially affected by the development of insulin resistance. Moreover, our data suggests insulin directly affects metabolomic measures previously associated with increased cardiovascular disease risk.

Association of measures of body fat with serum alpha-tocopherol and its metabolites in middle-aged individuals.

BACKGROUND AND AIMS: The accumulation of fat increases the formation of lipid peroxides, which are partly scavenged by alpha-tocopherol (α-TOH). Here, we aimed to investigate the associations between different measures of (abdominal) fat and levels of urinary α-TOH metabolites in middle-aged individuals. METHODS AND RESULTS: In this cross-sectional analysis in the Netherlands Epidemiology of Obesity study (N = 511, 53% women; mean [SD] age of 55 [6.1] years), serum α-TOH and α-TOH metabolites from 24-h urine were measured as alpha-tocopheronolactone hydroquinone (α-TLHQ, oxidized) and alpha-carboxymethyl-hydroxychroman (α-CEHC, enzymatically converted) using liquid-chromatography-tandem mass spectrometry. Body mass index and total body fat were measured, and abdominal subcutaneous and visceral adipose tissue (aSAT and VAT) were assessed using magnetic resonance imaging. Using multivariable-adjusted linear regression analyses, we analysed the associations of BMI, TBF, aSAT and VAT with levels of urinary α-TOH metabolites, adjusted for confounders. We observed no evidence for associations between body fat measures and serum α-TOH. Higher BMI and TBF were associated with lower urinary levels of TLHQ (0.95 [95%CI: 0.90, 1.00] and 0.94 [0.88, 1.01] times per SD, respectively) and with lower TLHQ relative to CEHC (0.93 [0.90, 0.98] and 0.93 [0.87, 0.98] times per SD, respectively). We observed similar associations for VAT (TLHQ: 0.94 [0.89, 0.99] times per SD), but not for aSAT. CONCLUSIONS: Opposite to our research hypothesis, higher abdominal adiposity was moderately associated with lower levels of oxidized α-TOH metabolites, which might reflect lower vitamin E antioxidative activity in individuals with higher abdominal fat instead.

Low-Frequency Synonymous Coding Variation in CYP2R1 Has Large Effects on Vitamin D Levels and Risk of Multiple Sclerosis.

Vitamin D insufficiency is common, correctable, and influenced by genetic factors, and it has been associated with risk of several diseases. We sought to identify low-frequency genetic variants that strongly increase the risk of vitamin D insufficiency and tested their effect on risk of multiple sclerosis, a disease influenced by low vitamin D concentrations. We used whole-genome sequencing data from 2,619 individuals through the UK10K program and deep-imputation data from 39,655 individuals genotyped genome-wide. Meta-analysis of the summary statistics from 19 cohorts identified in CYP2R1 the low-frequency (minor allele frequency = 2.5%) synonymous coding variant g.14900931G>A (p.Asp120Asp) (rs117913124[A]), which conferred a large effect on 25-hydroxyvitamin D (25OHD) levels (-0.43 SD of standardized natural log-transformed 25OHD per A allele; p value = 1.5 × 10-88). The effect on 25OHD was four times larger and independent of the effect of a previously described common variant near CYP2R1. By analyzing 8,711 individuals, we showed that heterozygote carriers of this low-frequency variant have an increased risk of vitamin D insufficiency (odds ratio [OR] = 2.2, 95% confidence interval [CI] = 1.78-2.78, p = 1.26 × 10-12). Individuals carrying one copy of this variant also had increased odds of multiple sclerosis (OR = 1.4, 95% CI = 1.19-1.64, p = 2.63 × 10-5) in a sample of 5,927 case and 5,599 control subjects. In conclusion, we describe a low-frequency CYP2R1 coding variant that exerts the largest effect upon 25OHD levels identified to date in the general European population and implicates vitamin D in the etiology of multiple sclerosis.