Human Gain-of-Function MC4R Variants Show Signaling Bias and Protect against Obesity.

The melanocortin 4 receptor (MC4R) is a G protein-coupled receptor whose disruption causes obesity. We functionally characterized 61 MC4R variants identified in 0.5 million people from UK Biobank and examined their associations with body mass index (BMI) and obesity-related cardiometabolic diseases. We found that the maximal efficacy of ß-arrestin recruitment to MC4R, rather than canonical Gαs-mediated cyclic adenosine-monophosphate production, explained 88% of the variance in the association of MC4R variants with BMI. While most MC4R variants caused loss of function, a subset caused gain of function; these variants were associated with significantly lower BMI and lower odds of obesity, type 2 diabetes, and coronary artery disease. Protective associations were driven by MC4R variants exhibiting signaling bias toward ß-arrestin recruitment and increased mitogen-activated protein kinase pathway activation. Harnessing ß-arrestin-biased MC4R signaling may represent an effective strategy for weight loss and the treatment of obesity-related cardiometabolic diseases.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites.

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

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.

Metabolomic differences in lung function metrics: evidence from two cohorts.

RATIONALE: The biochemical mechanisms underlying lung function are incompletely understood. OBJECTIVES: To identify and validate the plasma metabolome of lung function using two independent adult cohorts: discovery-the European Prospective Investigation into Cancer-Norfolk (EPIC-Norfolk, n=10 460) and validation-the VA Normative Aging Study (NAS) metabolomic cohort (n=437). METHODS: We ran linear regression models for 693 metabolites to identify associations with forced expiratory volume in one second (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC), in EPIC-Norfolk then validated significant findings in NAS. Significance in EPIC-Norfolk was denoted using an effective number of tests threshold of 95%; a metabolite was considered validated in NAS if the direction of effect was consistent and p<0.05. MEASUREMENTS AND MAIN RESULTS: Of 156 metabolites that associated with FEV1 in EPIC-Norfolk after adjustment for age, sex, body mass index, height, smoking and asthma status, 34 (21.8%) validated in NAS, including several metabolites involved in oxidative stress. When restricting the discovery sample to men only, a similar percentage, 18 of 79 significant metabolites (22.8%) were validated. A smaller number of metabolites were validated for FEV1/FVC, 6 of 65 (9.2%) when including all EPIC-Norfolk as the discovery population, and 2 of 34 (5.9%) when restricting to men. These metabolites were characterised by involvement in respiratory track secretants. Interestingly, no metabolites were validated for both FEV1 and FEV1/FVC. CONCLUSIONS: The validation of metabolites associated with respiratory function can help to better understand mechanisms of lung health and may assist the development of biomarkers.

The blood metabolome of incident kidney cancer: A case-control study nested within the MetKid consortium.

BACKGROUND: Excess bodyweight and related metabolic perturbations have been implicated in kidney cancer aetiology, but the specific molecular mechanisms underlying these relationships are poorly understood. In this study, we sought to identify circulating metabolites that predispose kidney cancer and to evaluate the extent to which they are influenced by body mass index (BMI). METHODS AND FINDINGS: We assessed the association between circulating levels of 1,416 metabolites and incident kidney cancer using pre-diagnostic blood samples from up to 1,305 kidney cancer case-control pairs from 5 prospective cohort studies. Cases were diagnosed on average 8 years after blood collection. We found 25 metabolites robustly associated with kidney cancer risk. In particular, 14 glycerophospholipids (GPLs) were inversely associated with risk, including 8 phosphatidylcholines (PCs) and 2 plasmalogens. The PC with the strongest association was PC ae C34:3 with an odds ratio (OR) for 1 standard deviation (SD) increment of 0.75 (95% confidence interval [CI]: 0.68 to 0.83, p = 2.6 × 10-8). In contrast, 4 amino acids, including glutamate (OR for 1 SD = 1.39, 95% CI: 1.20 to 1.60, p = 1.6 × 10-5), were positively associated with risk. Adjusting for BMI partly attenuated the risk association for some-but not all-metabolites, whereas other known risk factors of kidney cancer, such as smoking and alcohol consumption, had minimal impact on the observed associations. A mendelian randomisation (MR) analysis of the influence of BMI on the blood metabolome highlighted that some metabolites associated with kidney cancer risk are influenced by BMI. Specifically, elevated BMI appeared to decrease levels of several GPLs that were also found inversely associated with kidney cancer risk (e.g., -0.17 SD change [ßBMI] in 1-(1-enyl-palmitoyl)-2-linoleoyl-GPC (P-16:0/18:2) levels per SD change in BMI, p = 3.4 × 10-5). BMI was also associated with increased levels of glutamate (ßBMI: 0.12, p = 1.5 × 10-3). While our results were robust across the participating studies, they were limited to study participants of European descent, and it will, therefore, be important to evaluate if our findings can be generalised to populations with different genetic backgrounds. CONCLUSIONS: This study suggests a potentially important role of the blood metabolome in kidney cancer aetiology by highlighting a wide range of metabolites associated with the risk of developing kidney cancer and the extent to which changes in levels of these metabolites are driven by BMI-the principal modifiable risk factor of kidney cancer.

The association between circulating 25-hydroxyvitamin D metabolites and type 2 diabetes in European populations: A meta-analysis and Mendelian randomisation analysis.

BACKGROUND: Prior research suggested a differential association of 25-hydroxyvitamin D (25(OH)D) metabolites with type 2 diabetes (T2D), with total 25(OH)D and 25(OH)D3 inversely associated with T2D, but the epimeric form (C3-epi-25(OH)D3) positively associated with T2D. Whether or not these observational associations are causal remains uncertain. We aimed to examine the potential causality of these associations using Mendelian randomisation (MR) analysis. METHODS AND FINDINGS: We performed a meta-analysis of genome-wide association studies for total 25(OH)D (N = 120,618), 25(OH)D3 (N = 40,562), and C3-epi-25(OH)D3 (N = 40,562) in participants of European descent (European Prospective Investigation into Cancer and Nutrition [EPIC]-InterAct study, EPIC-Norfolk study, EPIC-CVD study, Ely study, and the SUNLIGHT consortium). We identified genetic variants for MR analysis to investigate the causal association of the 25(OH)D metabolites with T2D (including 80,983 T2D cases and 842,909 non-cases). We also estimated the observational association of 25(OH)D metabolites with T2D by performing random effects meta-analysis of results from previous studies and results from the EPIC-InterAct study. We identified 10 genetic loci associated with total 25(OH)D, 7 loci associated with 25(OH)D3 and 3 loci associated with C3-epi-25(OH)D3. Based on the meta-analysis of observational studies, each 1-standard deviation (SD) higher level of 25(OH)D was associated with a 20% lower risk of T2D (relative risk [RR]: 0.80; 95% CI 0.77, 0.84; p < 0.001), but a genetically predicted 1-SD increase in 25(OH)D was not significantly associated with T2D (odds ratio [OR]: 0.96; 95% CI 0.89, 1.03; p = 0.23); this result was consistent across sensitivity analyses. In EPIC-InterAct, 25(OH)D3 (per 1-SD) was associated with a lower risk of T2D (RR: 0.81; 95% CI 0.77, 0.86; p < 0.001), while C3-epi-25(OH)D3 (above versus below lower limit of quantification) was positively associated with T2D (RR: 1.12; 95% CI 1.03, 1.22; p = 0.006), but neither 25(OH)D3 (OR: 0.97; 95% CI 0.93, 1.01; p = 0.14) nor C3-epi-25(OH)D3 (OR: 0.98; 95% CI 0.93, 1.04; p = 0.53) was causally associated with T2D risk in the MR analysis. Main limitations include the lack of a non-linear MR analysis and of the generalisability of the current findings from European populations to other populations of different ethnicities. CONCLUSIONS: Our study found discordant associations of biochemically measured and genetically predicted differences in blood 25(OH)D with T2D risk. The findings based on MR analysis in a large sample of European ancestry do not support a causal association of total 25(OH)D or 25(OH)D metabolites with T2D and argue against the use of vitamin D supplementation for the prevention of T2D.

Characterization of missing values in untargeted MS-based metabolomics data and evaluation of missing data handling strategies.

BACKGROUND: Untargeted mass spectrometry (MS)-based metabolomics data often contain missing values that reduce statistical power and can introduce bias in biomedical studies. However, a systematic assessment of the various sources of missing values and strategies to handle these data has received little attention. Missing data can occur systematically, e.g. from run day-dependent effects due to limits of detection (LOD); or it can be random as, for instance, a consequence of sample preparation. METHODS: We investigated patterns of missing data in an MS-based metabolomics experiment of serum samples from the German KORA F4 cohort (n = 1750). We then evaluated 31 imputation methods in a simulation framework and biologically validated the results by applying all imputation approaches to real metabolomics data. We examined the ability of each method to reconstruct biochemical pathways from data-driven correlation networks, and the ability of the method to increase statistical power while preserving the strength of established metabolic quantitative trait loci. RESULTS: Run day-dependent LOD-based missing data accounts for most missing values in the metabolomics dataset. Although multiple imputation by chained equations performed well in many scenarios, it is computationally and statistically challenging. K-nearest neighbors (KNN) imputation on observations with variable pre-selection showed robust performance across all evaluation schemes and is computationally more tractable. CONCLUSION: Missing data in untargeted MS-based metabolomics data occur for various reasons. Based on our results, we recommend that KNN-based imputation is performed on observations with variable pre-selection since it showed robust results in all evaluation schemes.

Association of Genetic Variants Related to Gluteofemoral vs Abdominal Fat Distribution With Type 2 Diabetes, Coronary Disease, and Cardiovascular Risk Factors.

Importance: Body fat distribution, usually measured using waist-to-hip ratio (WHR), is an important contributor to cardiometabolic disease independent of body mass index (BMI). Whether mechanisms that increase WHR via lower gluteofemoral (hip) or via higher abdominal (waist) fat distribution affect cardiometabolic risk is unknown. Objective: To identify genetic variants associated with higher WHR specifically via lower gluteofemoral or higher abdominal fat distribution and estimate their association with cardiometabolic risk. Design, Setting, and Participants: Genome-wide association studies (GWAS) for WHR combined data from the UK Biobank cohort and summary statistics from previous GWAS (data collection: 2006-2018). Specific polygenic scores for higher WHR via lower gluteofemoral or via higher abdominal fat distribution were derived using WHR-associated genetic variants showing specific association with hip or waist circumference. Associations of polygenic scores with outcomes were estimated in 3 population-based cohorts, a case-cohort study, and summary statistics from 6 GWAS (data collection: 1991-2018). Exposures: More than 2.4 million common genetic variants (GWAS); polygenic scores for higher WHR (follow-up analyses). Main Outcomes and Measures: BMI-adjusted WHR and unadjusted WHR (GWAS); compartmental fat mass measured by dual-energy x-ray absorptiometry, systolic and diastolic blood pressure, low-density lipoprotein cholesterol, triglycerides, fasting glucose, fasting insulin, type 2 diabetes, and coronary disease risk (follow-up analyses). Results: Among 452 302 UK Biobank participants of European ancestry, the mean (SD) age was 57 (8) years and the mean (SD) WHR was 0.87 (0.09). In genome-wide analyses, 202 independent genetic variants were associated with higher BMI-adjusted WHR (n = 660 648) and unadjusted WHR (n = 663 598). In dual-energy x-ray absorptiometry analyses (n = 18 330), the hip- and waist-specific polygenic scores for higher WHR were specifically associated with lower gluteofemoral and higher abdominal fat, respectively. In follow-up analyses (n = 636 607), both polygenic scores were associated with higher blood pressure and triglyceride levels and higher risk of diabetes (waist-specific score: odds ratio [OR], 1.57 [95% CI, 1.34-1.83], absolute risk increase per 1000 participant-years [ARI], 4.4 [95% CI, 2.7-6.5], P < .001; hip-specific score: OR, 2.54 [95% CI, 2.17-2.96], ARI, 12.0 [95% CI, 9.1-15.3], P < .001) and coronary disease (waist-specific score: OR, 1.60 [95% CI, 1.39-1.84], ARI, 2.3 [95% CI, 1.5-3.3], P < .001; hip-specific score: OR, 1.76 [95% CI, 1.53-2.02], ARI, 3.0 [95% CI, 2.1-4.0], P < .001), per 1-SD increase in BMI-adjusted WHR. Conclusions and Relevance: Distinct genetic mechanisms may be linked to gluteofemoral and abdominal fat distribution that are the basis for the calculation of the WHR. These findings may improve risk assessment and treatment of diabetes and coronary disease.

Genetic Predisposition to an Impaired Metabolism of the Branched-Chain Amino Acids and Risk of Type 2 Diabetes: A Mendelian Randomisation Analysis.

BACKGROUND: Higher circulating levels of the branched-chain amino acids (BCAAs; i.e., isoleucine, leucine, and valine) are strongly associated with higher type 2 diabetes risk, but it is not known whether this association is causal. We undertook large-scale human genetic analyses to address this question. METHODS AND FINDINGS: Genome-wide studies of BCAA levels in 16,596 individuals revealed five genomic regions associated at genome-wide levels of significance (p < 5 × 10-8). The strongest signal was 21 kb upstream of the PPM1K gene (beta in standard deviations [SDs] of leucine per allele = 0.08, p = 3.9 × 10-25), encoding an activator of the mitochondrial branched-chain alpha-ketoacid dehydrogenase (BCKD) responsible for the rate-limiting step in BCAA catabolism. In another analysis, in up to 47,877 cases of type 2 diabetes and 267,694 controls, a genetically predicted difference of 1 SD in amino acid level was associated with an odds ratio for type 2 diabetes of 1.44 (95% CI 1.26-1.65, p = 9.5 × 10-8) for isoleucine, 1.85 (95% CI 1.41-2.42, p = 7.3 × 10-6) for leucine, and 1.54 (95% CI 1.28-1.84, p = 4.2 × 10-6) for valine. Estimates were highly consistent with those from prospective observational studies of the association between BCAA levels and incident type 2 diabetes in a meta-analysis of 1,992 cases and 4,319 non-cases. Metabolome-wide association analyses of BCAA-raising alleles revealed high specificity to the BCAA pathway and an accumulation of metabolites upstream of branched-chain alpha-ketoacid oxidation, consistent with reduced BCKD activity. Limitations of this study are that, while the association of genetic variants appeared highly specific, the possibility of pleiotropic associations cannot be entirely excluded. Similar to other complex phenotypes, genetic scores used in the study captured a limited proportion of the heritability in BCAA levels. Therefore, it is possible that only some of the mechanisms that increase BCAA levels or affect BCAA metabolism are implicated in type 2 diabetes. CONCLUSIONS: Evidence from this large-scale human genetic and metabolomic study is consistent with a causal role of BCAA metabolism in the aetiology of type 2 diabetes.

Association Between Low-Density Lipoprotein Cholesterol-Lowering Genetic Variants and Risk of Type 2 Diabetes: A Meta-analysis.

Importance: Low-density lipoprotein cholesterol (LDL-C)-lowering alleles in or near NPC1L1 or HMGCR, encoding the respective molecular targets of ezetimibe and statins, have previously been used as proxies to study the efficacy of these lipid-lowering drugs. Alleles near HMGCR are associated with a higher risk of type 2 diabetes, similar to the increased incidence of new-onset diabetes associated with statin treatment in randomized clinical trials. It is unknown whether alleles near NPC1L1 are associated with the risk of type 2 diabetes. Objective: To investigate whether LDL-C-lowering alleles in or near NPC1L1 and other genes encoding current or prospective molecular targets of lipid-lowering therapy (ie, HMGCR, PCSK9, ABCG5/G8, LDLR) are associated with the risk of type 2 diabetes. Design, Setting, and Participants: The associations with type 2 diabetes and coronary artery disease of LDL-C-lowering genetic variants were investigated in meta-analyses of genetic association studies. Meta-analyses included 50 775 individuals with type 2 diabetes and 270 269 controls and 60 801 individuals with coronary artery disease and 123 504 controls. Data collection took place in Europe and the United States between 1991 and 2016. Exposures: Low-density lipoprotein cholesterol-lowering alleles in or near NPC1L1, HMGCR, PCSK9, ABCG5/G8, and LDLR. Main Outcomes and Measures: Odds ratios (ORs) for type 2 diabetes and coronary artery disease. Results: Low-density lipoprotein cholesterol-lowering genetic variants at NPC1L1 were inversely associated with coronary artery disease (OR for a genetically predicted 1-mmol/L [38.7-mg/dL] reduction in LDL-C of 0.61 [95% CI, 0.42-0.88]; P = .008) and directly associated with type 2 diabetes (OR for a genetically predicted 1-mmol/L reduction in LDL-C of 2.42 [95% CI, 1.70-3.43]; P < .001). For PCSK9 genetic variants, the OR for type 2 diabetes per 1-mmol/L genetically predicted reduction in LDL-C was 1.19 (95% CI, 1.02-1.38; P = .03). For a given reduction in LDL-C, genetic variants were associated with a similar reduction in coronary artery disease risk (I2 = 0% for heterogeneity in genetic associations; P = .93). However, associations with type 2 diabetes were heterogeneous (I2 = 77.2%; P = .002), indicating gene-specific associations with metabolic risk of LDL-C-lowering alleles. Conclusions and Relevance: In this meta-analysis, exposure to LDL-C-lowering genetic variants in or near NPC1L1 and other genes was associated with a higher risk of type 2 diabetes. These data provide insights into potential adverse effects of LDL-C-lowering therapy.