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.

Normal range CAG repeat size variations in the HTT gene are associated with an adverse lipoprotein profile partially mediated by body mass index.

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

Circulating metabolites modulated by diet are associated with depression.

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Metabolomics dissection of depression heterogeneity and related cardiometabolic risk.

BACKGROUND: A recent hypothesis postulates the existence of an 'immune-metabolic depression' (IMD) dimension characterized by metabolic dysregulations. Combining data on metabolomics and depressive symptoms, we aimed to identify depressions associated with an increased risk of adverse metabolic alterations. METHOD: Clustering data were from 1094 individuals with major depressive disorder in the last 6 months and measures of 149 metabolites from a 1H-NMR platform and 30 depressive symptoms (IDS-SR30). Canonical correlation analyses (CCA) were used to identify main independent metabolite-symptom axes of variance. Then, for the replication, we examined the association of the identified dimensions with metabolites from the same platform and cardiometabolic diseases in an independent population-based cohort (n = 6572). RESULTS: CCA identified an overall depression dimension and a dimension resembling IMD, in which symptoms such as sleeping too much, increased appetite, and low energy level had higher relative loading. In the independent sample, the overall depression dimension was associated with lower cardiometabolic risk, such as (i.e. per s.d.) HOMA-1B -0.06 (95% CI -0.09 - -0.04), and visceral adipose tissue -0.10 cm2 (95% CI -0.14 - -0.07). In contrast, the IMD dimension was associated with well-known cardiometabolic diseases such as higher visceral adipose tissue 0.08 cm2 (95% CI 0.04-0.12), HOMA-1B 0.06 (95% CI 0.04-0.09), and lower HDL-cholesterol levels -0.03 mmol/L (95% CI -0.05 - -0.01). CONCLUSIONS: Combining metabolomics and clinical symptoms we identified a replicable depression dimension associated with adverse metabolic alterations, in line with the IMD hypothesis. Patients with IMD may be at higher cardiometabolic risk and may benefit from specific treatment targeting underlying metabolic dysregulations.

The association between adiposity and atypical energy-related symptoms of depression: A role for metabolic dysregulations.

BACKGROUND: Adiposity has been shown to be linked with atypical energy-related symptoms (AES) of depression. We used genomics to separate the effect of adiposity from that of metabolic dysregulations to examine whether the link between obesity and AES is dependent on the presence of metabolic dysregulations. METHOD: Data were from NEO (n = 5734 individuals) and NESDA (n = 2238 individuals) cohorts, in which the Inventory of Depressive Symptomatology (IDS-SR30) was assessed. AES profile was based on four symptoms: increased appetite, increased weight, low energy level, and leaden paralysis. We estimated associations between AES and two genetic risk scores (GRS) indexing increasing total body fat with (metabolically unhealthy adiposity, GRS-MUA) and without (metabolically healthy adiposity, GRS-MHA) metabolic dysregulations. RESULTS: We validated that both GRS-MUA and GRS-MHA were associated with higher total body fat in NEO study, but divergently associated with biomarkers of metabolic health (e.g., fasting glucose and HDL-cholesterol) in both cohorts. In the pooled results, per standard deviation, GRS-MUA was specifically associated with a higher AES score (ß = 0.03, 95%CI: 0.01; 0.05), while there was no association between GRS-MHA and AES (ß = -0.01, 95%CI: -0.03; 0.01). CONCLUSION: These results suggest that the established link between adiposity and AES profile emerges in the presence of metabolic dysregulations, which may represent the connecting substrate between the two conditions.

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.

Evaluation of the Value of Waist Circumference and Metabolomics in the Estimation of Visceral Adipose Tissue.

Visceral adipose tissue (VAT) is a strong prognostic factor for cardiovascular disease and a potential target for cardiovascular risk stratification. Because VAT is difficult to measure in clinical practice, we estimated prediction models with predictors routinely measured in general practice and VAT as outcome using ridge regression in 2,501 middle-aged participants from the Netherlands Epidemiology of Obesity study, 2008-2012. Adding waist circumference and other anthropometric measurements on top of the routinely measured variables improved the optimism-adjusted R2 from 0.50 to 0.58 with a decrease in the root-mean-square error (RMSE) from 45.6 to 41.5 cm2 and with overall good calibration. Further addition of predominantly lipoprotein-related metabolites from the Nightingale platform did not improve the optimism-corrected R2 and RMSE. The models were externally validated in 370 participants from the Prospective Investigation of Vasculature in Uppsala Seniors (PIVUS, 2006-2009) and 1,901 participants from the Multi-Ethnic Study of Atherosclerosis (MESA, 2000-2007). Performance was comparable to the development setting in PIVUS (R2 = 0.63, RMSE = 42.4 cm2, calibration slope = 0.94) but lower in MESA (R2 = 0.44, RMSE = 60.7 cm2, calibration slope = 0.75). Our findings indicate that the estimation of VAT with routine clinical measurements can be substantially improved by incorporating waist circumference but not by metabolite measurements.

Plasma Proteomics of Renal Function: A Transethnic Meta-Analysis and Mendelian Randomization Study.

BACKGROUND: Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European populations, and causality of the associations between plasma proteins and kidney function has never been addressed. METHODS: A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR. RESULTS: In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (SPOCK2) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR. CONCLUSIONS: In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

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.

The effect of physical activity level and exercise training on the association between plasma branched-chain amino acids and intrahepatic lipid content in participants with obesity.

AIMS: To evaluate whether the association between plasma branched-chain amino acids (BCAA) and intrahepatic lipid (IHL) was affected by physical activity level. Furthermore, to investigate if a conventional exercise training program, a subcategory of physical activity, could lower plasma BCAA along with alterations in IHL content in patients with type 2 diabetes (T2DM) and people with nonalcoholic fatty liver (NAFL). METHODS: To investigate the effect of physical activity on the association between plasma BCAA and IHL content, linear regression analyses were performed in 1983 individuals from the Netherlands Epidemiology of Obesity (NEO) stratified by physical activity frequency. Furthermore, the effect of a 12-week supervised combined aerobic resistance-exercise program on plasma BCAA, insulin sensitivity (hyperinsulinemic-euglycemic clamp), and IHL (proton-magnetic resonance spectroscopy (1H-MRS)) was investigated in seven patients with T2DM, seven individuals with NAFL and seven BMI-matched control participants (CON). RESULTS: We observed positive associations between plasma valine, isoleucine and leucine level, and IHL content (1.29 (95% CI: 1.21, 1.38), 1.52 (95% CI: 1.43, 1.61), and 1.54 (95% CI: 1.44, 1.64) times IHL, respectively, per standard deviation of plasma amino acid level). Similar associations were observed in less active versus more active individuals. Exercise training did not change plasma BCAA levels among groups, but reduced IHL content in NAFL (from 11.6 ± 3.0% pre-exercise to 8.1 ± 2.0% post exercise, p < 0.05) and CON (from 2.4 ± 0.6% pre-exercise to 1.6 ± 1.4% post exercise, p < 0.05), and improved peripheral insulin sensitivity in NAFL as well by ~23% (p < 0.05). CONCLUSIONS: The association between plasma BCAA levels and IHL is not affected by physical activity level. Exercise training reduced IHL without affecting plasma BCAA levels in individuals with NAFL and CON. We conclude that exercise training-induced reduction in IHL content is not related to changes in plasma BCAA levels. TRIAL REGISTRATION: Trial registry number: NCT01317576.

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.

Deep molecular phenotypes link complex disorders and physiological insult to CpG methylation.

Epigenetic regulation of cellular function provides a mechanism for rapid organismal adaptation to changes in health, lifestyle and environment. Associations of cytosine-guanine di-nucleotide (CpG) methylation with clinical endpoints that overlap with metabolic phenotypes suggest a regulatory role for these CpG sites in the body's response to disease or environmental stress. We previously identified 20 CpG sites in an epigenome-wide association study (EWAS) with metabolomics that were also associated in recent EWASs with diabetes-, obesity-, and smoking-related endpoints. To elucidate the molecular pathways that connect these potentially regulatory CpG sites to the associated disease or lifestyle factors, we conducted a multi-omics association study including 2474 mass-spectrometry-based metabolites in plasma, urine and saliva, 225 NMR-based lipid and metabolite measures in blood, 1124 blood-circulating proteins using aptamer technology, 113 plasma protein N-glycans and 60 IgG-glyans, using 359 samples from the multi-ethnic Qatar Metabolomics Study on Diabetes (QMDiab). We report 138 multi-omics associations at these CpG sites, including diabetes biomarkers at the diabetes-associated TXNIP locus, and smoking-specific metabolites and proteins at multiple smoking-associated loci, including AHRR. Mendelian randomization suggests a causal effect of metabolite levels on methylation of obesity-associated CpG sites, i.e. of glycerophospholipid PC(O-36: 5), glycine and a very low-density lipoprotein (VLDL-A) on the methylation of the obesity-associated CpG loci DHCR24, MYO5C and CPT1A, respectively. Taken together, our study suggests that multi-omics-associated CpG methylation can provide functional read-outs for the underlying regulatory response mechanisms to disease or environmental insults.

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.

Metabolomics analyses in non-diabetic middle-aged individuals reveal metabolites impacting early glucose disturbances and insulin sensitivity.

INTRODUCTION: Several plasma metabolites have been associated with insulin resistance and type 2 diabetes mellitus. OBJECTIVES: We aimed to identify plasma metabolites associated with different indices of early disturbances in glucose metabolism and insulin sensitivity. METHODS: This cross-sectional study was conducted in a subsample of the Leiden Longevity Study comprising individuals without a history of diabetes mellitus (n = 233) with a mean age of 63.3 ± 6.7 years of which 48.1% were men. We tested for associations of fasting glucose, fasting insulin, HOMA-IR, Matsuda Index, Insulinogenic Index and glycated hemoglobin with metabolites (Swedish Metabolomics Platform) using linear regression analysis adjusted for age, sex and BMI. Results were validated internally using an independent metabolomics platform (Biocrates platform) and replicated externally in the independent Netherlands Epidemiology of Obesity (NEO) study (Metabolon platform) (n = 545, mean age of 55.8 ± 6.0 years of which 48.6% were men). Moreover, in the NEO study, we replicated our analyses in individuals with diabetes mellitus (cases: n = 36; controls = 561). RESULTS: Out of the 34 metabolites, a total of 12 plasma metabolites were associated with different indices of disturbances in glucose metabolism and insulin sensitivity in individuals without diabetes mellitus. These findings were validated using a different metabolomics platform as well as in an independent cohort of non-diabetics. Moreover, tyrosine, alanine, valine, tryptophan and alpha-ketoglutaric acid levels were higher in individuals with diabetes mellitus. CONCLUSION: We found several plasma metabolites that are associated with early disturbances in glucose metabolism and insulin sensitivity of which five were also higher in individuals with diabetes mellitus.

Proteome-wide assessment of diabetes mellitus in Qatari identifies IGFBP-2 as a risk factor already with early glycaemic disturbances.

BACKGROUND: Proteomics is expected to provide novel insights in the underlying pathophysiology of type 2 diabetes mellitus. In the present study, we aimed to identify and biochemically characterize proteins associated with diabetes mellitus in a Qatari population. METHODS: In a diabetes case-control study (175 cases, 164 controls; Arab, South Asian and Philippine ethnicities), we conducted a discovery study to screen 1141 blood protein levels for associations with diabetes mellitus. Additional analyses were done in controls in relation to Hb1Ac, and biochemical characterization of the main findings was performed with metabolomics (501 metabolites). We performed two-sample Mendelian Randomization to provide evidence of potential causality using data from European descent of the DIAGRAM consortium (74,124 cases of diabetes mellitus and 824,006 controls) for the identified proteins for T2D and Hb1Ac. RESULTS: After accounting for multiple testing, 30 protein levels were different (p-values<8.6e-5) between cases and controls. Of these, a higher Hb1Ac in controls was associated with a lower IGFBP-2 level (p-value = 4.1e-6). IGFBP-2 protein level was found lower among cases compared with controls across all ethnicities. In controls, IGFBP-2 was associated with 21 metabolite levels, but specifically connected to the metabolite citrulline in network analyses. We observed no evidence, however, that the association between IGFBP-2 and diabetes mellitus was causal. CONCLUSIONS: We specifically identified IGFBP-2 to be associated with diabetes mellitus, although with no evidence for causality, which was specifically connected to citrulline metabolism.

Genome-wide meta-analysis of macronutrient intake of 91,114 European ancestry participants from the cohorts for heart and aging research in genomic epidemiology consortium.

Macronutrient intake, the proportion of calories consumed from carbohydrate, fat, and protein, is an important risk factor for metabolic diseases with significant familial aggregation. Previous studies have identified two genetic loci for macronutrient intake, but incomplete coverage of genetic variation and modest sample sizes have hindered the discovery of additional loci. Here, we expanded the genetic landscape of macronutrient intake, identifying 12 suggestively significant loci (P < 1 × 10-6) associated with intake of any macronutrient in 91,114 European ancestry participants. Four loci replicated and reached genome-wide significance in a combined meta-analysis including 123,659 European descent participants, unraveling two novel loci; a common variant in RARB locus for carbohydrate intake and a rare variant in DRAM1 locus for protein intake, and corroborating earlier FGF21 and FTO findings. In additional analysis of 144,770 participants from the UK Biobank, all identified associations from the two-stage analysis were confirmed except for DRAM1. Identified loci might have implications in brain and adipose tissue biology and have clinical impact in obesity-related phenotypes. Our findings provide new insight into biological functions related to macronutrient intake.

Associations between outdoor temperature and bright sunlight with metabolites in two population-based European cohorts.

BACKGROUND AND AIMS: Outdoor temperature and bright sunlight may directly and/or indirectly modulate systemic metabolism. We assessed the associations between outdoor temperature and bright sunlight duration with metabolomics. METHODS AND RESULTS: Cross-sectional analyses were undertaken in non-diabetic individuals from the Oxford BioBank (OBB; N = 6368; mean age 47.0 years, males 44%) and the Netherlands Epidemiology of Obesity (NEO; N = 5916; mean age 55.6 years, males 43%) study. Data on mean outdoor bright sunlight and temperature were collected from local weather stations in the week prior to blood sampling. Fasting serum levels of 148 metabolites, including 14 lipoprotein subclasses, were measured using NMR spectroscopy. Linear regression analyses were performed to assess the associations between mean outdoor temperature and bright sunlight duration with metabolomics adjusted for age, sex, body mass index, season and either outdoor temperature or bright sunlight. A higher mean outdoor temperature was associated with increased serum concentrations of lipoprotein (sub)particles (ß (SE) = 0.064 (0.018) SD per 5 °C, p = 5.03e-4) and certain amino acids such as phenylalanine (0.066 (0.016) SD, p = 6.44e-05) and leucine (0.111 (0.018) SD, p = 1.25e-09). In contrast, longer duration of bright sunlight was specifically associated with lower concentrations of very low-density lipoprotein (sub)particles (e.g., VLDL cholesterol (-0.024 (0.005) SD per 1-h bright sunlight, p = 8.06e-6)). The direction of effects was generally consistent between the OBB and NEO, although effect sizes were generally larger in the OBB. CONCLUSIONS: Increased bright sunlight duration is associated with an improved metabolic profile whilst higher outdoor temperature may adversely impact cardiometabolic health.

Whole-Exome Sequencing Identifies Loci Associated with Blood Cell Traits and Reveals a Role for Alternative GFI1B Splice Variants in Human Hematopoiesis.

Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis.