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.

Plasma Metabolites Related to the Consumption of Different Types of Dairy Products and Their Association with New-Onset Type 2 Diabetes: Analyses in the Fenland and EPIC-Norfolk Studies, United Kingdom.

SCOPE: To identify metabolites associated with habitual dairy consumption and investigate their associations with type 2 diabetes (T2D) risk. METHODS AND RESULTS: Metabolomics assays were conducted in the Fenland (n = 10,281) and EPIC-Norfolk (n = 1,440) studies. Using 82 metabolites assessed in both studies, we developed metabolite scores to classify self-reported consumption of milk, yogurt, cheese, butter, and total dairy (Fenland Study-discovery set; n = 6035). Internal and external validity of the scores was evaluated (Fenland-validation set, n = 4246; EPIC-Norfolk, n = 1440). The study assessed associations between each metabolite score and T2D incidence in EPIC-Norfolk (n = 641 cases; 16,350 person-years). The scores classified low and high consumers for all dairy types with internal validity, and milk, butter, and total dairy with external validity. The scores were further associated with lower incident T2D: hazard ratios (95% confidence interval) per standard deviation: milk 0.71 (0.65, 0.77); butter 0.62 (0.57, 0.68); total dairy 0.66 (0.60, 0.72). These associations persisted after adjustment for known dairy-fat biomarkers. CONCLUSION: Metabolite scores identified habitual consumers of milk, butter, and total dairy products, and were associated with lower T2D risk. These findings hold promise for identifying objective indicators of the physiological response to dairy consumption.

Causal effects of maternal circulating amino acids on offspring birthweight: a Mendelian randomisation study.

BACKGROUND: Amino acids are key to protein synthesis, energy metabolism, cell signaling and gene expression; however, the contribution of specific maternal amino acids to fetal growth is unclear. METHODS: We explored the effect of maternal circulating amino acids on fetal growth, proxied by birthweight, using two-sample Mendelian randomisation (MR) and summary data from a genome-wide association study (GWAS) of serum amino acids levels (sample 1, n = 86,507) and a maternal GWAS of offspring birthweight in UK Biobank and Early Growth Genetics Consortium, adjusting for fetal genotype effects (sample 2, n = 406,063 with maternal and/or fetal genotype effect estimates). A total of 106 independent single nucleotide polymorphisms robustly associated with 19 amino acids (p < 4.9 × 10-10) were used as genetic instrumental variables (IV). Wald ratio and inverse variance weighted methods were used in MR main analysis. A series of sensitivity analyses were performed to explore IV assumption violations. FINDINGS: Our results provide evidence that maternal circulating glutamine (59 g offspring birthweight increase per standard deviation increase in maternal amino acid level, 95% CI: 7, 110) and serine (27 g, 95% CI: 9, 46) raise, while leucine (-59 g, 95% CI: -106, -11) and phenylalanine (-25 g, 95% CI: -47, -4) lower offspring birthweight. These findings are supported by sensitivity analyses. INTERPRETATION: Our findings strengthen evidence for key roles of maternal circulating amino acids during pregnancy in healthy fetal growth. FUNDING: A full list of funding bodies that contributed to this study can be found under Acknowledgments.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases.

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

Development and validation of a metabolite score for red meat intake: an observational cohort study and randomized controlled dietary intervention.

BACKGROUND: Self-reported meat consumption is associated with disease risk but objective assessment of different dimensions of this heterogeneous dietary exposure in observational and interventional studies remains challenging. OBJECTIVES: We aimed to derive and validate scores based on plasma metabolites for types of meat consumption. For the most predictive score, we aimed to test whether the included metabolites varied with change in meat consumption, and whether the score was associated with incidence of type 2 diabetes (T2D) and other noncommunicable diseases. METHODS: We derived scores based on 781 plasma metabolites for red meat, processed meat, and poultry consumption assessed with 7-d food records among 11,432 participants in the EPIC-Norfolk (European Prospective Investigation into Cancer and Nutrition-Norfolk) cohort. The scores were then tested for internal validity in an independent subset (n = 853) of the same cohort. In focused analysis on the red meat metabolite score, we examined whether the metabolites constituting the score were also associated with meat intake in a randomized crossover dietary intervention trial of meat (n = 12, Lyon, France). In the EPIC-Norfolk study, we assessed the association of the red meat metabolite score with T2D incidence (n = 1478) and other health endpoints. RESULTS: The best-performing score was for red meat, comprising 139 metabolites which accounted for 17% of the explained variance of red meat consumption in the validation set. In the intervention, 11 top-ranked metabolites in the red meat metabolite score increased significantly after red meat consumption. In the EPIC-Norfolk study, the red meat metabolite score was associated with T2D incidence (adjusted HR per SD: 1.17; 95% CI: 1.10, 1.24). CONCLUSIONS: The red meat metabolite score derived and validated in this study contains metabolites directly derived from meat consumption and is associated with T2D risk. These findings suggest the potential for objective assessment of dietary components and their application for understanding diet-disease associations.The trial in Lyon, France, was registered at clinicaltrials.gov as NCT03354130.

The Relationship of Maternal Gestational Mass Spectrometry-Derived Metabolites with Offspring Congenital Heart Disease: Results from Multivariable and Mendelian Randomization Analyses.

Background: It is plausible that maternal pregnancy metabolism influences the risk of offspring congenital heart disease (CHD). We sought to explore this through a systematic approach using different methods and data. Methods: We undertook multivariable logistic regression of the odds of CHD for 923 mass spectrometry (MS)-derived metabolites in a sub-sample of a UK birth cohort (Born in Bradford (BiB); N = 2605, 46 CHD cases). We considered metabolites reaching a p-value threshold <0.05 to be suggestively associated with CHD. We sought validation of our findings, by repeating the multivariable regression analysis within the BiB cohort for any suggestively associated metabolite that was measured by nuclear magnetic resonance (NMR) or clinical chemistry (N = 7296, 87 CHD cases), and by using genetic risk scores (GRS: weighted genetic risk scores of single nucleotide polymorphisms (SNPs) that were associated with any suggestive metabolite) in Mendelian randomization (MR) analyses. The MR analyses were performed in BiB and two additional European birth cohorts (N = 38,662, 319 CHD cases). Results: In the main multivariable analyses, we identified 44 metabolites suggestively associated with CHD, including those from the following super pathways: amino acids, lipids, co-factors and vitamins, xenobiotics, nucleotides, energy, and several unknown molecules. Of these 44, isoleucine and leucine were available in the larger BiB cohort (NMR), and for these the results were validated. The MR analyses were possible for 27/44 metabolites and for 11 there was consistency with the multivariable regression results. Conclusions: In summary, we have used complimentary data sources and statistical techniques to construct layers of evidence. We found that pregnancy amino acid metabolism, androgenic steroid lipids, and levels of succinylcarnitine could be important contributing factors for CHD.

Metabolomic profiling reveals extensive adrenal suppression due to inhaled corticosteroid therapy in asthma.

The application of large-scale metabolomic profiling provides new opportunities for realizing the potential of omics-based precision medicine for asthma. By leveraging data from over 14,000 individuals in four distinct cohorts, this study identifies and independently replicates 17 steroid metabolites whose levels were significantly reduced in individuals with prevalent asthma. Although steroid levels were reduced among all asthma cases regardless of medication use, the largest reductions were associated with inhaled corticosteroid (ICS) treatment, as confirmed in a 4-year low-dose ICS clinical trial. Effects of ICS treatment on steroid levels were dose dependent; however, significant reductions also occurred with low-dose ICS treatment. Using information from electronic medical records, we found that cortisol levels were substantially reduced throughout the entire 24-hour daily period in patients with asthma who were treated with ICS compared to those who were untreated and to patients without asthma. Moreover, patients with asthma who were treated with ICS showed significant increases in fatigue and anemia as compared to those without ICS treatment. Adrenal suppression in patients with asthma treated with ICS might, therefore, represent a larger public health problem than previously recognized. Regular cortisol monitoring of patients with asthma treated with ICS is needed to provide the optimal balance between minimizing adverse effects of adrenal suppression while capitalizing on the established benefits of ICS treatment.

Identification of Rare Loss-of-Function Genetic Variation Regulating Body Fat Distribution.

CONTEXT: Biological and translational insights from large-scale, array-based genetic studies of fat distribution, a key determinant of metabolic health, have been limited by the difficulty in linking predominantly noncoding variants to specific gene targets. Rare coding variant analyses provide greater confidence that a specific gene is involved, but do not necessarily indicate whether gain or loss of function (LoF) would be of most therapeutic benefit. OBJECTIVE: This work aimed to identify genes/proteins involved in determining fat distribution. METHODS: We combined the power of genome-wide analysis of array-based rare, nonsynonymous variants in 450 562 individuals in the UK Biobank with exome-sequence-based rare LoF gene burden testing in 184 246 individuals. RESULTS: The data indicate that the LoF of 4 genes (PLIN1 [LoF variants, P = 5.86 × 10-7], INSR [LoF variants, P = 6.21 × 10-7], ACVR1C [LoF + moderate impact variants, P = 1.68 × 10-7; moderate impact variants, P = 4.57 × 10-7], and PDE3B [LoF variants, P = 1.41 × 10-6]) is associated with a beneficial effect on body mass index-adjusted waist-to-hip ratio and increased gluteofemoral fat mass, whereas LoF of PLIN4 (LoF variants, P = 5.86 × 10-7 adversely affects these parameters. Phenotypic follow-up suggests that LoF of PLIN1, PDE3B, and ACVR1C favorably affects metabolic phenotypes (eg, triglycerides [TGs] and high-density lipoprotein [HDL] cholesterol concentrations) and reduces the risk of cardiovascular disease, whereas PLIN4 LoF has adverse health consequences. INSR LoF is associated with lower TG and HDL levels but may increase the risk of type 2 diabetes. CONCLUSION: This study robustly implicates these genes in the regulation of fat distribution, providing new and in some cases somewhat counterintuitive insight into the potential consequences of targeting these molecules therapeutically.

Rare and common genetic determinants of metabolic individuality and their effects on human health.

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

Plasma metabolites to profile pathways in noncommunicable disease multimorbidity.

Multimorbidity, the simultaneous presence of multiple chronic conditions, is an increasing global health problem and research into its determinants is of high priority. We used baseline untargeted plasma metabolomics profiling covering >1,000 metabolites as a comprehensive readout of human physiology to characterize pathways associated with and across 27 incident noncommunicable diseases (NCDs) assessed using electronic health record hospitalization and cancer registry data from over 11,000 participants (219,415 person years). We identified 420 metabolites shared between at least 2 NCDs, representing 65.5% of all 640 significant metabolite-disease associations. We integrated baseline data on over 50 diverse clinical risk factors and characteristics to identify actionable shared pathways represented by those metabolites. Our study highlights liver and kidney function, lipid and glucose metabolism, low-grade inflammation, surrogates of gut microbial diversity and specific health-related behaviors as antecedents of common NCD multimorbidity with potential for early prevention. We integrated results into an open-access webserver ( https://omicscience.org/apps/mwasdisease/ ) to facilitate future research and meta-analyses.

Genetically Predicted Glucose-Dependent Insulinotropic Polypeptide (GIP) Levels and Cardiovascular Disease Risk Are Driven by Distinct Causal Variants in the GIPR Region.

There is considerable interest in GIPR agonism to enhance the insulinotropic and extrapancreatic effects of GIP, thereby improving glycemic and weight control in type 2 diabetes (T2D) and obesity. Recent genetic epidemiological evidence has implicated higher GIPR-mediated GIP levels in raising coronary artery disease (CAD) risk, a potential safety concern for GIPR agonism. We therefore aimed to quantitatively assess whether the association between higher GIPR-mediated fasting GIP levels and CAD risk is mediated via GIPR or is instead the result of linkage disequilibrium (LD) confounding between variants at the GIPR locus. Using Bayesian multitrait colocalization, we identified a GIPR missense variant, rs1800437 (G allele; E354), as the putatively causal variant shared among fasting GIP levels, glycemic traits, and adiposity-related traits (posterior probability for colocalization [PPcoloc] > 0.97; PP explained by the candidate variant [PPexplained] = 1) that was independent from a cluster of CAD and lipid traits driven by a known missense variant in APOE (rs7412; distance to E354 ∼770 Kb; R 2 with E354 = 0.004; PPcoloc > 0.99; PPexplained = 1). Further, conditioning the association between E354 and CAD on the residual LD with rs7412, we observed slight attenuation in association, but it remained significant (odds ratio [OR] per copy of E354 after adjustment 1.03; 95% CI 1.02, 1.04; P = 0.003). Instead, E354's association with CAD was completely attenuated when conditioning on an additional established CAD signal, rs1964272 (R 2 with E354 = 0.27), an intronic variant in SNRPD2 (OR for E354 after adjustment for rs1964272: 1.01; 95% CI 0.99, 1.03; P = 0.06). We demonstrate that associations with GIP and anthropometric and glycemic traits are driven by genetic signals distinct from those driving CAD and lipid traits in the GIPR region and that higher E354-mediated fasting GIP levels are not associated with CAD risk. These findings provide evidence that the inclusion of GIPR agonism in dual GIPR/GLP1R agonists could potentiate the protective effect of GLP-1 agonists on diabetes without undue CAD risk, an aspect that has yet to be assessed in clinical trials.