A cross-platform approach identifies genetic regulators of human metabolism and health

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.

Genetic associations and potential mediators between psychiatric disorders and irritable bowel syndrome: a Mendelian randomization study with mediation analysis

Objective

Potential causal associations between psychiatric disorders and irritable bowel syndrome have been demonstrated in observational studies; however, these studies are susceptible to underlying confounding and reverse causation biases. We aimed to assess the causal effects of psychiatric disorders on irritable bowel syndrome (IBS) and the potential mediators from a genetic perspective by conducting a Mendelian randomization (MR) study with mediation analysis.

Method

Genetic instruments associated with psychiatric disorders, potential mediators, and IBS were obtained from large-scale genome-wide association studies (GWAS). Three MR methods - the inverse-variance weighted (IVW) method, MR-Egger method, and weighted median method, were used to investigate causal association estimates. Heterogeneity among different genetic instrumental variables (IVs) was assessed using Q tests. Additionally, the MR-PRESSO and MR-Pleiotropy methods were used to verify horizontal pleiotropy and detect outliers that might bias the results, which were removed from further analysis. Consequently, we used MR mediation analysis to investigate potential mediators in the causal associations between psychiatric disorders and IBS.

Results

MR provided evidence of the causal effects of genetically predicted broad depression, major depressive disorder (MDD), anxiety disorder, post-traumatic stress disorder (PTSD), and schizophrenia on IBS. The results of MR mediation analysis demonstrated that the reduction in acetate levels mediated 12.6% of the effects of broad depression on IBS; insomnia mediated 16.00%, 16.20%, and 27.14% of the effects of broad depression, MDD, and PTSD on IBS, respectively; and the increase in blood β-hydroxybutyrate levels mediated 50.76% of the effects of schizophrenia on IBS.

Conclusion

Our study confirmed the brain-gut axis involvement and potential modulators in the pathophysiology of psychiatric disorder-induced IBS from a genetic perspective, and suggests potential therapeutic targets for the disrupted brain-gut axis.

Three linked variants have opposing regulatory effects on isovaleryl-CoA dehydrogenase gene expression

While genome-wide association studies (GWAS) and positive selection scans identify genomic loci driving human phenotypic diversity, functional validation is required to discover the variant(s) responsible. We dissected the IVD gene locus-which encodes the isovaleryl-CoA dehydrogenase enzyme-implicated by selection statistics, multiple GWAS, and clinical genetics as important to function and fitness. We combined luciferase assays, CRISPR/Cas9 genome-editing, massively parallel reporter assays (MPRA), and a deletion tiling MPRA strategy across regulatory loci. We identified three regulatory variants, including an indel, that may underpin GWAS signals for pulmonary fibrosis and testosterone, and that are linked on a positively selected haplotype in the Japanese population. These regulatory variants exhibit synergistic and opposing effects on IVD expression experimentally. Alleles at these variants lie on a haplotype tagged by the variant most strongly associated with IVD expression and metabolites, but with no functional evidence itself. This work demonstrates how comprehensive functional investigation and multiple technologies are needed to discover the true genetic drivers of phenotypic diversity.

A Genome-Wide Association Study of Serum Metabolite Profiles in Septic Shock Patients

OBJECTIVES

We sought to assess whether genetic associations with metabolite concentrations in septic shock patients could be used to identify pathways of potential importance for understanding sepsis pathophysiology.

DESIGN

Retrospective multicenter cohort studies of septic shock patients.

SETTING

All participants who were admitted to 27 participating hospital sites in three countries (Australia, New Zealand, and the United Kingdom) were eligible for inclusion.

PATIENTS

Adult, critically ill, mechanically ventilated patients with septic shock (n = 230) who were a subset of the Adjunctive Corticosteroid Treatment in Critically Ill Patients with Septic Shock trial (ClinicalTrials.gov number: NCT01448109).

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

A genome-wide association study was conducted for a range of serum metabolite levels for participants. Genome-wide significant associations (p ≤ 5 × 10-8) were found for the two major ketone bodies (3-hydroxybutyrate [rs2456680] and acetoacetate [rs2213037] and creatinine (rs6851961). One of these single-nucleotide polymorphisms (SNPs) (rs2213037) was located in the alcohol dehydrogenase cluster of genes, which code for enzymes related to the metabolism of acetoacetate and, therefore, presents a plausible association for this metabolite. None of the three SNPs showed strong associations with risk of sepsis, 28- or 90-day mortality, or Acute Physiology and Chronic Health Evaluation score (a measure of sepsis severity).

CONCLUSIONS

We suggest that the genetic associations with metabolites may reflect a starvation response rather than processes involved in sepsis pathophysiology. However, our results require further investigation and replication in both healthy and diseased cohorts including those of different ancestry.

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.

A systematic review of metabolomic findings in adult and pediatric renal disease

Chronic kidney disease (CKD) affects over 0.5 billion people worldwide across their lifetimes. Despite a growingly ageing world population, an increase in all-age prevalence of kidney disease persists. Adult-onset forms of kidney disease often result from lifestyle-modifiable metabolic illnesses such as type 2 diabetes. Pediatric and adolescent forms of renal disease are primarily caused by morphological abnormalities of the kidney, as well as immunological, infectious and inherited metabolic disorders. Alterations in energy metabolism are observed in CKD of varying causes, albeit the molecular mechanisms underlying pathology are unclear. A systematic indexing of metabolites identified in plasma and urine of patients with kidney disease alongside disease enrichment analysis uncovered inborn errors of metabolism as a framework that links features of adult and pediatric kidney disease. The relationship of genetics and metabolism in kidney disease could be classified into three distinct landscapes: (i) Normal genotypes that develop renal damage because of lifestyle and / or comorbidities; (ii) Heterozygous genetic variants and polymorphisms that result in unique metabotypes that may predispose to the development of kidney disease via synergistic heterozygosity, and (iii) Homozygous genetic variants that cause renal impairment by perturbing metabolism, as found in children with monogenic inborn errors of metabolism. Interest in the identification of early biomarkers of onset and progression of CKD has grown steadily in the last years, though it has not translated into clinical routine yet. This systematic review indexes findings of differential concentration of metabolites and energy pathway dysregulation in kidney disease and appraises their potential use as biomarkers.

Multi-omic prediction of incident type 2 diabetes

AIMS/HYPOTHESIS

The identification of people who are at high risk of developing type 2 diabetes is a key part of population-level prevention strategies. Previous studies have evaluated the predictive utility of omics measurements, such as metabolites, proteins or polygenic scores, but have considered these separately. The improvement that combined omics biomarkers can provide over and above current clinical standard models is unclear. The aim of this study was to test the predictive performance of genome, proteome, metabolome and clinical biomarkers when added to established clinical prediction models for type 2 diabetes.

METHODS

We developed sparse interpretable prediction models in a prospective, nested type 2 diabetes case-cohort study (N=1105, incident type 2 diabetes cases=375) with 10,792 person-years of follow-up, selecting from 5759 features across the genome, proteome, metabolome and clinical biomarkers using least absolute shrinkage and selection operator (LASSO) regression. We compared the predictive performance of omics-derived predictors with a clinical model including the variables from the Cambridge Diabetes Risk Score and HbA1c.

RESULTS

Among single omics prediction models that did not include clinical risk factors, the top ten proteins alone achieved the highest performance (concordance index [C index]=0.82 [95% CI 0.75, 0.88]), suggesting the proteome as the most informative single omic layer in the absence of clinical information. However, the largest improvement in prediction of type 2 diabetes incidence over and above the clinical model was achieved by the top ten features across several omic layers (C index=0.87 [95% CI 0.82, 0.92], Δ C index=0.05, p=0.045). This improvement by the top ten omic features was also evident in individuals with HbA1c <42 mmol/mol (6.0%), the threshold for prediabetes (C index=0.84 [95% CI 0.77, 0.90], Δ C index=0.07, p=0.03), the group in whom prediction would be most useful since they are not targeted for preventative interventions by current clinical guidelines. In this subgroup, the type 2 diabetes polygenic risk score was the major contributor to the improvement in prediction, and achieved a comparable improvement in performance when added onto the clinical model alone (C index=0.83 [95% CI 0.75, 0.90], Δ C index=0.06, p=0.002). However, compared with those with prediabetes, individuals at high polygenic risk in this group had only around half the absolute risk for type 2 diabetes over a 20 year period.

CONCLUSIONS/INTERPRETATION

Omic approaches provided marginal improvements in prediction of incident type 2 diabetes. However, while a polygenic risk score does improve prediction in people with an HbA1c in the normoglycaemic range, the group in whom prediction would be most useful, even individuals with a high polygenic burden in that subgroup had a low absolute type 2 diabetes risk. This suggests a limited feasibility of implementing targeted population-based genetic screening for preventative interventions.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Effect of Genetic and Dietary Perturbation of Glycine Metabolism on Atherosclerosis in Humans and Mice

Objective

Epidemiological and genetic studies have reported inverse associations between circulating glycine levels and risk of coronary artery disease (CAD). However, these findings have not been consistently observed in all studies. We sought to evaluate the causal relationship between circulating glycine levels and atherosclerosis using large-scale genetic analyses in humans and dietary supplementation experiments in mice.

Methods

Serum glycine levels were evaluated for association with prevalent and incident CAD in the UK Biobank. A multi-ancestry genome-wide association study (GWAS) meta-analysis was carried out to identify genetic determinants for circulating glycine levels, which were then used to evaluate the causal relationship between glycine and risk of CAD by Mendelian randomization (MR). A glycine feeding study was carried out with atherosclerosis-prone apolipoprotein E deficient (ApoE-/-) mice to determine the effects of increased circulating glycine levels on amino acid metabolism, metabolic traits, and aortic lesion formation.

Results

Among 105,718 subjects from the UK Biobank, elevated serum glycine levels were associated with significantly reduced risk of prevalent CAD (Quintile 5 vs. Quintile 1 OR=0.76, 95% CI 0.67-0.87; P<0.0001) and incident CAD (Quintile 5 vs. Quintile 1 HR=0.70, 95% CI 0.65-0.77; P<0.0001) in models adjusted for age, sex, ethnicity, anti-hypertensive and lipid-lowering medications, blood pressure, kidney function, and diabetes. A meta-analysis of 13 GWAS datasets (total n=230,947) identified 61 loci for circulating glycine levels, of which 26 were novel. MR analyses provided modest evidence that genetically elevated glycine levels were causally associated with reduced systolic blood pressure and risk of type 2 diabetes, but did provide evidence for an association with risk of CAD. Furthermore, glycine-supplementation in ApoE-/- mice did not alter cardiometabolic traits, inflammatory biomarkers, or development of atherosclerotic lesions.

Conclusions

Circulating glycine levels were inversely associated with risk of prevalent and incident CAD in a large population-based cohort. While substantially expanding the genetic architecture of circulating glycine levels, MR analyses and in vivo feeding studies in humans and mice, respectively, did not provide evidence that the clinical association of this amino acid with CAD represents a causal relationship, despite being associated with two correlated risk factors.

GeneMAP: A discovery platform for metabolic gene function

Organisms maintain metabolic homeostasis through the combined functions of small molecule transporters and enzymes. While many of the metabolic components have been well-established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Associations Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions, even pinpointing genes that are distant from the variants implicated by GWAS. In particular, our work identified SLC25A48 as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite, betaine. Rare variant testing and polygenic risk score analyses have elucidated choline-relevant phenomic consequences of SLC25A48 dysfunction. Altogether, our study proposes SLC25A48 as a mitochondrial choline transporter and provides a discovery platform for metabolic gene function.

Genetically predicted plasma levels of amino acids and metabolic dysfunction-associated fatty liver disease risk: a Mendelian randomization study

BACKGROUND

Emerging metabolomics-based studies suggested links between amino acid metabolism and metabolic dysfunction-associated fatty liver disease (MAFLD) risk; however, whether there exists an aetiological role of amino acid metabolism in MAFLD development remains unknown. The aim of the present study was to assess the causal relationship between circulating levels of amino acids and MAFLD risk.

METHODS

We conducted a two-sample Mendelian randomization (MR) analysis using summary-level data from genome-wide association studies (GWAS) to evaluate the causal relationship between genetically predicted circulating levels of amino acids and the risk of MAFLD. In the discovery MR analysis, we used data from the largest MAFLD GWAS (8434 cases and 770,180 controls), while in the replication MR analysis, we used data from a GWAS on MAFLD (1483 cases and 17,781 controls) where MAFLD cases were diagnosed using liver biopsy. We used Wald ratios or inverse variance-weighted (IVW) methods in the MR main analysis and weighted median and MR-Egger regression analyses in sensitivity analyses. Furthermore, we performed a conservative MR analysis by restricting genetic instruments to those directly involved in amino acid metabolism pathways.

RESULTS

We found that genetically predicted higher alanine (OR = 1.43, 95% CI 1.13-1.81) and lower glutamine (OR = 0.83, 95% CI 0.73-0.96) levels were associated with a higher risk of developing MAFLD based on the results from the MR main and conservative analysis. The results from MR sensitivity analyses and complementary analysis using liver proton density fat fraction as a continuous outcome proxying for MAFLD supported the main findings.

CONCLUSIONS

Novel causal metabolites related to MAFLD development were uncovered through MR analysis, suggesting future potential for evaluating these metabolites as targets for MAFLD prevention or treatment.

Assessing the causal effect of genetically predicted metabolites and metabolic pathways on stroke

BACKGROUND

Stroke is a common neurological disorder that disproportionately affects middle-aged and elderly individuals, leading to significant disability and mortality. Recently, human blood metabolites have been discovered to be useful in unraveling the underlying biological mechanisms of neurological disorders. Therefore, we aimed to evaluate the causal relationship between human blood metabolites and susceptibility to stroke.

METHODS

Summary data from genome-wide association studies (GWASs) of serum metabolites and stroke and its subtypes were obtained separately. A total of 486 serum metabolites were used as the exposure. Simultaneously, 11 different stroke phenotypes were set as the outcomes, including any stroke (AS), any ischemic stroke (AIS), large artery stroke (LAS), cardioembolic stroke (CES), small vessel stroke (SVS), lacunar stroke (LS), white matter hyperintensities (WMH), intracerebral hemorrhage (ICH), subarachnoid hemorrhage (SAH), transient ischemic attack (TIA), and brain microbleeds (BMB). A two-sample Mendelian randomization (MR) study was conducted to investigate the causal effects of serum metabolites on stroke and its subtypes. The inverse variance-weighted MR analyses were conducted as causal estimates, accompanied by a series of sensitivity analyses to evaluate the robustness of the results. Furthermore, a reverse MR analysis was conducted to assess the potential for reverse causation. Additionally, metabolic pathway analysis was performed using the web-based MetOrigin.

RESULTS

After correcting for the false discovery rate (FDR), MR analysis results revealed remarkable causative associations with 25 metabolites. Further sensitivity analyses confirmed that only four causative associations involving three specific metabolites passed all sensitivity tests, namely ADpSGEGDFXAEGGGVR* for AS (OR: 1.599, 95% CI 1.283-1.993, p = 2.92 × 10-5) and AIS (OR: 1.776, 95% CI 1.380-2.285, p = 8.05 × 10-6), 1-linoleoylglycerophosph-oethanolamine* for LAS (OR: 0.198, 95% CI 0.091-0.428, p = 3.92 × 10-5), and gamma-glutamylmethionine* for SAH (OR: 3.251, 95% CI 1.876-5.635, p = 2.66 × 10-5), thereby demonstrating a high degree of stability. Moreover, eight causative associations involving seven other metabolites passed both sensitivity tests and were considered robust. The association result of one metabolite (glutamate for LAS) was considered non-robust. As for the remaining metabolites, we speculate that they may potentially possess underlying causal relationships. Notably, no common metabolites emerged from the reverse MR analysis. Moreover, after FDR correction, metabolic pathway analysis identified 40 significant pathways across 11 stroke phenotypes.

CONCLUSIONS

The identified metabolites and their associated metabolic pathways are promising circulating metabolic biomarkers, holding potential for their application in stroke screening and preventive strategies within clinical settings.

Genome-wide association analysis of plasma lipidome identifies 495 genetic associations

The human plasma lipidome captures risk for cardiometabolic diseases. To discover new lipid-associated variants and understand the link between lipid species and cardiometabolic disorders, we perform univariate and multivariate genome-wide analyses of 179 lipid species in 7174 Finnish individuals. We fine-map the associated loci, prioritize genes, and examine their disease links in 377,277 FinnGen participants. We identify 495 genome-trait associations in 56 genetic loci including 8 novel loci, with a considerable boost provided by the multivariate analysis. For 26 loci, fine-mapping identifies variants with a high causal probability, including 14 coding variants indicating likely causal genes. A phenome-wide analysis across 953 disease endpoints reveals disease associations for 40 lipid loci. For 11 coronary artery disease risk variants, we detect strong associations with lipid species. Our study demonstrates the power of multivariate genetic analysis in correlated lipidomics data and reveals genetic links between diseases and lipid species beyond the standard lipids.

The Physiological and Pathological Role of Acyl-CoA Oxidation

Fatty acid metabolism, including β-oxidation (βOX), plays an important role in human physiology and pathology. βOX is an essential process in the energy metabolism of most human cells. Moreover, βOX is also the source of acetyl-CoA, the substrate for (a) ketone bodies synthesis, (b) cholesterol synthesis, (c) phase II detoxication, (d) protein acetylation, and (d) the synthesis of many other compounds, including N-acetylglutamate-an important regulator of urea synthesis. This review describes the current knowledge on the importance of the mitochondrial and peroxisomal βOX in various organs, including the liver, heart, kidney, lung, gastrointestinal tract, peripheral white blood cells, and other cells. In addition, the diseases associated with a disturbance of fatty acid oxidation (FAO) in the liver, heart, kidney, lung, alimentary tract, and other organs or cells are presented. Special attention was paid to abnormalities of FAO in cancer cells and the diseases caused by mutations in gene-encoding enzymes involved in FAO. Finally, issues related to α- and ω- fatty acid oxidation are discussed.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Protein-metabolite association studies identify novel proteomic determinants of metabolite levels in human plasma

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.

Association between Neuroligin-1 polymorphism and plasma glutamine levels in individuals with autism spectrum disorder

BACKGROUND

Unravelling the relationships between candidate genes and autism spectrum disorder (ASD) phenotypes remains an outstanding challenge. Endophenotypes, defined as inheritable, measurable quantitative traits, might provide intermediary links between genetic risk factors and multifaceted ASD phenotypes. In this study, we sought to determine whether plasma metabolite levels could serve as endophenotypes in individuals with ASD and their family members.

METHODS

We employed an untargeted, high-resolution metabolomics platform to analyse 14,342 features across 1099 plasma samples. These samples were collected from probands and their family members participating in the Autism Genetic Resource Exchange (AGRE) (N = 658), compared with neurotypical individuals enrolled in the PrecisionLink Health Discovery (PLHD) program at Boston Children's Hospital (N = 441). We conducted a metabolite quantitative trait loci (mQTL) analysis using whole-genome genotyping data from each cohort in AGRE and PLHD, aiming to prioritize significant mQTL and metabolite pairs that were exclusively observed in AGRE.

FINDINGS

Within the AGRE group, we identified 54 significant associations between genotypes and metabolite levels (P < 5.27 × 10-11), 44 of which were not observed in the PLHD group. Plasma glutamine levels were found to be associated with variants in the NLGN1 gene, a gene that encodes post-synaptic cell-adhesion molecules in excitatory neurons. This association was not detected in the PLHD group. Notably, a significant negative correlation between plasma glutamine and glutamate levels was observed in the AGRE group, but not in the PLHD group. Furthermore, plasma glutamine levels showed a negative correlation with the severity of restrictive and repetitive behaviours (RRB) in ASD, although no direct association was observed between RRB severity and the NLGN1 genotype.

INTERPRETATION

Our findings suggest that plasma glutamine levels could potentially serve as an endophenotype, thus establishing a link between the genetic risk associated with NLGN1 and the severity of RRB in ASD. This identified association could facilitate the development of novel therapeutic targets, assist in selecting specific cohorts for clinical trials, and provide insights into target symptoms for future ASD treatment strategies.

FUNDING

This work was supported by the National Institute of Health (grant numbers: R01MH107205, U01TR002623, R24OD024622, OT2OD032720, and R01NS129188) and the PrecisionLink Biobank for Health Discovery at Boston Children's Hospital.

A genome-wide association study identifies 41 loci associated with eicosanoid levels

Eicosanoids are biologically active derivatives of polyunsaturated fatty acids with broad relevance to health and disease. We report a genome-wide association study in 8406 participants of the Atherosclerosis Risk in Communities Study, identifying 41 loci associated with 92 eicosanoids and related metabolites. These findings highlight loci required for eicosanoid biosynthesis, including FADS1-3, ELOVL2, and numerous CYP450 loci. In addition, significant associations implicate a range of non-oxidative lipid metabolic processes in eicosanoid regulation, including at PKD2L1/SCD and several loci involved in fatty acyl-CoA metabolism. Further, our findings highlight select clearance mechanisms, for example, through the hepatic transporter encoded by SLCO1B1. Finally, we identify eicosanoids associated with aspirin and non-steroidal anti-inflammatory drug use and demonstrate the substantial impact of genetic variants even for medication-associated eicosanoids. These findings shed light on both known and unknown aspects of eicosanoid metabolism and motivate interest in several gene-eicosanoid associations as potential functional participants in human disease.

Large-scale metabolomic profiling and incident non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is a highly prevalent disease with no specific drug therapy. High-throughput metabolomics present an unprecedented opportunity to identify biomarkers and potentially causal risk factors for NAFLD. Here, we determined the impact of 21 circulating metabolites, 17 lipids, and 132 lipoprotein particle characteristics on NAFLD combining prospective observational and two-sample Mendelian randomization (MR) analyses in 121,032 UK Biobank participants. We identified several metabolic factors associated with NAFLD risk in observational and MR analyses including triglyceride-rich and high-density lipoprotein particles composition, as well as the ratio of polyunsaturated fatty acids to total fatty acids. This study, is one of the largest to investigate incident NAFLD, provides concordant observational and genetic evidence that therapies aimed at reducing circulating triglycerides and increasing large HDL particles, as well as interventions aimed at increasing polyunsaturated fatty acid content may warrant further investigation into NAFLD prevention and treatment.

Integrating genetics and metabolomics from multi-ethnic and multi-fluid data reveals putative mechanisms for age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness in older adults. Investigating shared genetic components between metabolites and AMD can enhance our understanding of its pathogenesis. We conduct metabolite genome-wide association studies (mGWASs) using multi-ethnic genetic and metabolomic data from up to 28,000 participants. With bidirectional Mendelian randomization analysis involving 16,144 advanced AMD cases and 17,832 controls, we identify 108 putatively causal relationships between plasma metabolites and advanced AMD. These metabolites are enriched in glycerophospholipid metabolism, lysophospholipid, triradylcglycerol, and long chain polyunsaturated fatty acid pathways. Bayesian genetic colocalization analysis and a customized metabolome-wide association approach prioritize putative causal AMD-associated metabolites. We find limited evidence linking urine metabolites to AMD risk. Our study emphasizes the contribution of plasma metabolites, particularly lipid-related pathways and genes, to AMD risk and uncovers numerous putative causal associations between metabolites and AMD risk.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Biological and functional multimorbidity-from mechanisms to management

Globally, the number of people with multiple co-occurring diseases will increase substantially over the coming decades, with important consequences for patients, carers, healthcare systems and society. Addressing this challenge requires a shift in the prevailing clinical, educational and scientific thinking and organization-with a strong emphasis on the maintenance of generalist skills to balance the specialization trends of medical education and research. Multimorbidity is not a single entity but differs quantitively and qualitatively across life stages, ethnicities, sexes, socioeconomic groups and geographies. Data-driven science that quantifies the impact of disease co-occurrence-beyond the small number of currently well-studied long-term conditions (such as cardiometabolic diseases)-can help illuminate the pathological diversity of multimorbidity and identify common, mechanistically related, and prognostically relevant clusters. Broader access to data opportunities across modalities and disciplines will catalyze vertical and horizontal integration of multimorbidity research, to enable reconfiguring of medical services, clinical trials, guidelines and research in a way that accounts for the complexity of multimorbidity-and provides efficient, joined-up services for patients.

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine

The kidneys operate at the interface of plasma and urine by clearing molecular waste products while retaining valuable solutes. Genetic studies of paired plasma and urine metabolomes may identify underlying processes. We conducted genome-wide studies of 1,916 plasma and urine metabolites and detected 1,299 significant associations. Associations with 40% of implicated metabolites would have been missed by studying plasma alone. We detected urine-specific findings that provide information about metabolite reabsorption in the kidney, such as aquaporin (AQP)-7-mediated glycerol transport, and different metabolomic footprints of kidney-expressed proteins in plasma and urine that are consistent with their localization and function, including the transporters NaDC3 (SLC13A3) and ASBT (SLC10A2). Shared genetic determinants of 7,073 metabolite-disease combinations represent a resource to better understand metabolic diseases and revealed connections of dipeptidase 1 with circulating digestive enzymes and with hypertension. Extending genetic studies of the metabolome beyond plasma yields unique insights into processes at the interface of body compartments.

Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations

Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.

Identification of biomarkers for glycaemic deterioration in type 2 diabetes

We identify biomarkers for disease progression in three type 2 diabetes cohorts encompassing 2,973 individuals across three molecular classes, metabolites, lipids and proteins. Homocitrulline, isoleucine and 2-aminoadipic acid, eight triacylglycerol species, and lowered sphingomyelin 42:2;2 levels are predictive of faster progression towards insulin requirement. Of ~1,300 proteins examined in two cohorts, levels of GDF15/MIC-1, IL-18Ra, CRELD1, NogoR, FAS, and ENPP7 are associated with faster progression, whilst SMAC/DIABLO, SPOCK1 and HEMK2 predict lower progression rates. In an external replication, proteins and lipids are associated with diabetes incidence and prevalence. NogoR/RTN4R injection improved glucose tolerance in high fat-fed male mice but impaired it in male db/db mice. High NogoR levels led to islet cell apoptosis, and IL-18R antagonised inflammatory IL-18 signalling towards nuclear factor kappa-B in vitro. This comprehensive, multi-disciplinary approach thus identifies biomarkers with potential prognostic utility, provides evidence for possible disease mechanisms, and identifies potential therapeutic avenues to slow diabetes progression.

Causal factors underlying diabetes risk informed by Mendelian randomisation analysis: evidence, opportunities and challenges

Diabetes and its complications cause a heavy disease burden globally. Identifying exposures, risk factors and molecular processes causally associated with the development of diabetes can provide important evidence bases for disease prevention and spur novel therapeutic strategies. Mendelian randomisation (MR), an epidemiological approach that uses genetic instruments to infer causal associations between an exposure and an outcome, can be leveraged to complement evidence from observational and clinical studies. This narrative review aims to summarise the evidence on potential causal risk factors for diabetes by integrating published MR studies on type 1 and 2 diabetes, and to reflect on future perspectives of MR studies on diabetes. Despite the genetic influence on type 1 diabetes, few MR studies have been conducted to identify causal exposures or molecular processes leading to increased disease risk. In type 2 diabetes, MR analyses support causal associations of somatic, mental and lifestyle factors with development of the disease. These studies have also identified biomarkers, some of them derived from the gut microbiota, and molecular processes leading to increased disease risk. These studies provide valuable data to better understand disease pathophysiology and explore potential therapeutic targets. Because genetic association studies have mostly been restricted to participants of European descent, multi-ancestry cohorts are needed to examine the role of different types of physical activity, dietary components, metabolites, protein biomarkers and gut microbiome in diabetes development.

Diabetic vascular diseases: molecular mechanisms and therapeutic strategies

Vascular complications of diabetes pose a severe threat to human health. Prevention and treatment protocols based on a single vascular complication are no longer suitable for the long-term management of patients with diabetes. Diabetic panvascular disease (DPD) is a clinical syndrome in which vessels of various sizes, including macrovessels and microvessels in the cardiac, cerebral, renal, ophthalmic, and peripheral systems of patients with diabetes, develop atherosclerosis as a common pathology. Pathological manifestations of DPDs usually manifest macrovascular atherosclerosis, as well as microvascular endothelial function impairment, basement membrane thickening, and microthrombosis. Cardiac, cerebral, and peripheral microangiopathy coexist with microangiopathy, while renal and retinal are predominantly microangiopathic. The following associations exist between DPDs: numerous similar molecular mechanisms, and risk-predictive relationships between diseases. Aggressive glycemic control combined with early comprehensive vascular intervention is the key to prevention and treatment. In addition to the widely recommended metformin, glucagon-like peptide-1 agonist, and sodium-glucose cotransporter-2 inhibitors, for the latest molecular mechanisms, aldose reductase inhibitors, peroxisome proliferator-activated receptor-γ agonizts, glucokinases agonizts, mitochondrial energy modulators, etc. are under active development. DPDs are proposed for patients to obtain more systematic clinical care requires a comprehensive diabetes care center focusing on panvascular diseases. This would leverage the advantages of a cross-disciplinary approach to achieve better integration of the pathogenesis and therapeutic evidence. Such a strategy would confer more clinical benefits to patients and promote the comprehensive development of DPD as a discipline.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

Exploiting the mediating role of the metabolome to unravel transcript-to-phenotype associations

Despite the success of genome-wide association studies (GWASs) in identifying genetic variants associated with complex traits, understanding the mechanisms behind these statistical associations remains challenging. Several methods that integrate methylation, gene expression, and protein quantitative trait loci (QTLs) with GWAS data to determine their causal role in the path from genotype to phenotype have been proposed. Here, we developed and applied a multi-omics Mendelian randomization (MR) framework to study how metabolites mediate the effect of gene expression on complex traits. We identified 216 transcript-metabolite-trait causal triplets involving 26 medically relevant phenotypes. Among these associations, 58% were missed by classical transcriptome-wide MR, which only uses gene expression and GWAS data. This allowed the identification of biologically relevant pathways, such as between ANKH and calcium levels mediated by citrate levels and SLC6A12 and serum creatinine through modulation of the levels of the renal osmolyte betaine. We show that the signals missed by transcriptome-wide MR are found, thanks to the increase in power conferred by integrating multiple omics layer. Simulation analyses show that with larger molecular QTL studies and in case of mediated effects, our multi-omics MR framework outperforms classical MR approaches designed to detect causal relationships between single molecular traits and complex phenotypes.

PPM1K-regulated impaired catabolism of branched-chain amino acids orchestrates polycystic ovary syndrome

BACKGROUND

Polycystic ovary syndrome (PCOS) is one of the most common diseases with the coexistence of reproductive malfunction and metabolic disorders. Previous studies have found increased branched chain amino acid (BCAA) levels in women with PCOS. However, it remains unclear whether BCAA metabolism is causally associated with the risk of PCOS.

METHODS

The changes of BCAA levels in the plasma and follicular fluids of PCOS women were detected. Mendelian randomization (MR) approaches were used to explore the potential causal association between BCAA levels and the risk of PCOS. The function of the gene coding the protein phosphatase Mg2+/Mn2+-dependent 1K (PPM1K) was further explored by using Ppm1k-deficient mouse model and PPM1K down-regulated human ovarian granulosa cells.

FINDINGS

BCAA levels were significantly elevated in both plasma and follicular fluids of PCOS women. Based on MR, a potential direct, causal role for BCAA metabolism was revealed in the pathogenesis of PCOS, and PPM1K was detected as a vital driver. Ppm1k-deficient female mice had increased BCAA levels and exhibited PCOS-like traits, including hyperandrogenemia and abnormal follicle development. A reduction in dietary BCAA intake significantly improved the endocrine and ovarian dysfunction of Ppm1k-/- female mice. Knockdown of PPM1K promoted the conversion of glycolysis to pentose phosphate pathway and inhibited mitochondrial oxidative phosphorylation in human granulosa cells.

INTERPRETATION

Ppm1k deficiency-impaired BCAA catabolism causes the occurrence and development of PCOS. PPM1K suppression disturbed energy metabolism homeostasis in the follicular microenvironment, which provided an underlying mechanism of abnormal follicle development.

FUNDING

This study was supported by the National Key Research and Development Program of China (2021YFC2700402, 2019YFA0802503), the National Natural Science Foundation of China (81871139, 82001503, 92057107), the CAMS Innovation Fund for Medical Sciences (2019-I2M-5-001), Key Clinical Projects of Peking University Third Hospital (BYSY2022043), the China Postdoctoral Science Foundation (2021T140600), and the Collaborative Innovation Program of Shanghai Municipal Health Commission (2020CXJQ01).

Promises and challenges in pharmacoepigenetics

Pharmacogenetics, the study of how interindividual genetic differences affect drug response, does not explain all observed heritable variance in drug response. Epigenetic mechanisms, such as DNA methylation, and histone acetylation may account for some of the unexplained variances. Epigenetic mechanisms modulate gene expression and can be suitable drug targets and can impact the action of nonepigenetic drugs. Pharmacoepigenetics is the field that studies the relationship between epigenetic variability and drug response. Much of this research focuses on compounds targeting epigenetic mechanisms, called epigenetic drugs, which are used to treat cancers, immune disorders, and other diseases. Several studies also suggest an epigenetic role in classical drug response; however, we know little about this area. The amount of information correlating epigenetic biomarkers to molecular datasets has recently expanded due to technological advances, and novel computational approaches have emerged to better identify and predict epigenetic interactions. We propose that the relationship between epigenetics and classical drug response may be examined using data already available by (1) finding regions of epigenetic variance, (2) pinpointing key epigenetic biomarkers within these regions, and (3) mapping these biomarkers to a drug-response phenotype. This approach expands on existing knowledge to generate putative pharmacoepigenetic relationships, which can be tested experimentally. Epigenetic modifications are involved in disease and drug response. Therefore, understanding how epigenetic drivers impact the response to classical drugs is important for improving drug design and administration to better treat disease.

Genetically determined serum serine level has a novel causal effect on multiple sclerosis risk and predicts disability progression

BACKGROUND

There are currently no specific biomarkers for multiple sclerosis (MS). Identifying robust biomarkers for MS is crucial to improve disease diagnosis and management.

METHODS

This study first used six Mendelian randomisation methods to assess causal relationship of 174 metabolites with MS, incorporating data from European-ancestry metabolomics (n=8569-86 507) and MS (n=14 802 MS cases, 26 703 controls) genomewide association studies. Genetic scores for identified causal metabolite(s) were then computed to predict MS disability progression in an independent longitudinal cohort (AusLong study) of 203 MS cases with up to 15-year follow-up.

RESULTS

We found a novel genetic causal effect of serine on MS onset (OR=1.67, 95% CI 1.51 to 1.84, p=1.73×10^-20), such that individuals whose serine level is 1 SD above the population mean will have 1.67 times the risk of developing MS. This is robust across all sensitivity methods (OR ranges from 1.49 to 1.67). In an independent longitudinal MS cohort, we then constructed time-dynamic and time-fixed genetic scores based on serine genetic instrument single-nucleotide polymorphisms, where higher scores for raised serum serine level were associated with increased risk of disability worsening, especially in the time-dynamic model (RR=1.25, 95% CI 1.10 to 1.42, p=7.52×10^-4).

CONCLUSIONS

These findings support investigating serine as an important candidate biomarker for MS onset and disability progression.

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.

Impact of the gut microbiota and associated metabolites on cardiometabolic traits, chronic diseases and human longevity: a Mendelian randomization study

Features of the gut microbiota have been associated with several chronic diseases and longevity in preclinical models as well as in observational studies. Whether these relations underlie causal effects in humans remains to be established. We aimed to determine whether the gut microbiota influences cardiometabolic traits as well as the risk of chronic diseases and human longevity using a comprehensive 2-Sample Mendelian randomization approach. We included as exposures 10 gut-associated metabolites and pathways and 57 microbial taxa abundance. We included as outcomes nine cardiometabolic traits (fasting glucose, fasting insulin, systolic blood pressure, diastolic blood pressure, HDL cholesterol, LDL cholesterol, triglycerides, estimated glomerular filtration rate, body mass index [BMI]), eight chronic diseases previously linked with the gut microbiota in observational studies (Alzheimer's disease, depression, type 2 diabetes, non-alcoholic fatty liver disease, coronary artery disease (CAD), stroke, osteoporosis and chronic kidney disease), as well as parental lifespan and longevity. We found 7 associations with evidence of causality before and after sensitivity analyses, but not after multiple testing correction (1198 tests). Most effect sizes (4/7) were small. The two largest exposure-outcome effects were markedly attenuated towards the null upon inclusion of BMI or alcohol intake frequency in multivariable MR analyses. While finding robust genetic instruments for microbiota features is challenging hence potentially inflating type 2 errors, these results do not support a large causal impact of human gut microbita features on cardiometabolic traits, chronic diseases or longevity. These results also suggest that the previously documented associations between gut microbiota and human health outcomes may not always underly causal relations.

Genetic influences on human blood metabolites in the Japanese population

An increase in ethnic diversity in genetic studies has the potential to provide unprecedented insights into how genetic variations influence human phenotypes. In this study, we conducted a quantitative trait locus (QTL) analysis of 121 metabolites measured using gas chromatography-mass spectrometry with plasma samples from 4,888 Japanese individuals. We found 60 metabolite-gene associations, of which 13 have not been previously reported. Meta-analyses with another Japanese and a European study identified six and two additional unreported loci, respectively. Genetic variants influencing metabolite levels were more enriched in protein-coding regions than in the regulatory regions while being associated with the risk of various diseases. Finally, we identified a signature of strong negative selection for uric acid ( S ˆ  = -1.53, p = 6.2 × 10^-18). Our study expanded the knowledge of genetic influences on human blood metabolites, providing valuable insights into their physiological, pathological, and selective properties.

Comprehensive investigation of pathway enrichment methods for functional interpretation of LC-MS global metabolomics data

BACKGROUND

Global or untargeted metabolomics is widely used to comprehensively investigate metabolic profiles under various pathophysiological conditions such as inflammations, infections, responses to exposures or interactions with microbial communities. However, biological interpretation of global metabolomics data remains a daunting task. Recent years have seen growing applications of pathway enrichment analysis based on putative annotations of liquid chromatography coupled with mass spectrometry (LC-MS) peaks for functional interpretation of LC-MS-based global metabolomics data. However, due to intricate peak-metabolite and metabolite-pathway relationships, considerable variations are observed among results obtained using different approaches. There is an urgent need to benchmark these approaches to inform the best practices.

RESULTS

We have conducted a benchmark study of common peak annotation approaches and pathway enrichment methods in current metabolomics studies. Representative approaches, including three peak annotation methods and four enrichment methods, were selected and benchmarked under different scenarios. Based on the results, we have provided a set of recommendations regarding peak annotation, ranking metrics and feature selection. The overall better performance was obtained for the mummichog approach. We have observed that a ~30% annotation rate is sufficient to achieve high recall (~90% based on mummichog), and using semi-annotated data improves functional interpretation. Based on the current platforms and enrichment methods, we further propose an identifiability index to indicate the possibility of a pathway being reliably identified. Finally, we evaluated all methods using 11 COVID-19 and 8 inflammatory bowel diseases (IBD) global metabolomics datasets.

From 'Omics to Multi-omics Technologies: the Discovery of Novel Causal Mediators

PURPOSE OF REVIEW

'Omics studies provide a comprehensive characterisation of a biological entity, such as the genome, epigenome, transcriptome, proteome, metabolome, or microbiome. This review covers the unique properties of these types of 'omics and their roles as causal mediators in cardiovascular disease. Moreover, applications and challenges of integrating multiple types of 'omics data to increase predictive power, improve causal inference, and elucidate biological mechanisms are discussed.

RECENT FINDINGS

Multi-omics approaches are growing in adoption as they provide orthogonal evidence and overcome the limitations of individual types of 'omics data. Studies with multiple types of 'omics data have improved the diagnosis and prediction of disease states and afforded a deeper understanding of underlying pathophysiological mechanisms, beyond any single type of 'omics data. For instance, disease-associated loci in the genome can be supplemented with other 'omics to prioritise causal genes and understand the function of non-coding variants. Alternatively, techniques, such as Mendelian randomisation, can leverage genetics to provide evidence supporting a causal role for disease-associated molecules, and elucidate their role in disease pathogenesis. As technologies improve, costs for 'omics studies will continue to fall and datasets will become increasingly accessible to researchers. The intrinsically unbiased nature of 'omics data is well-suited to exploratory analyses that discover causal mediators of disease, and multi-omics is an emerging discipline that leverages the strengths of each type of 'omics data to provide insights greater than the sum of its parts.

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.

Comprehensive characterization of putative genetic influences on plasma metabolome in a pediatric cohort

BACKGROUND

The human exposome is composed of diverse metabolites and small chemical compounds originated from endogenous and exogenous sources, respectively. Genetic and environmental factors influence metabolite levels, while the extent of genetic contributions across metabolic pathways is not yet known. Untargeted profiling of human metabolome using high-resolution mass spectrometry (HRMS) combined with genome-wide genotyping allows comprehensive identification of genetically influenced metabolites. As such previous studies of adults discovered and replicated genotype-metabotype associations. However, these associations have not been characterized in children.

RESULTS

We conducted the largest genome by metabolome-wide association study to date of children (N = 441) using 619,688 common genetic variants and 14,342 features measured by HRMS. Narrow-sense heritability (h²) estimates of plasma metabolite concentrations using genomic relatedness matrix restricted maximum likelihood (GREML) method showed a bimodal distribution with high h² (> 0.8) for 15.9% of features and low h² (< 0.2) for most of features (62.0%). The features with high h² were enriched for amino acid and nucleic acid metabolism, while carbohydrate and lipid concentrations showed low h². For each feature, a metabolite quantitative trait loci (mQTL) analysis was performed to identify genetic variants that were potentially associated with plasma levels. Fifty-four associations among 29 features and 43 genetic variants were identified at a genome-wide significance threshold p < 3.5 × 10^-12 (= 5 × 10^-8/14,342 features). Previously reported associations such as UGT1A1 and bilirubin; PYROXD2 and methyl lysine; and ACADS and butyrylcarnitine were successfully replicated in our pediatric cohort. We found potential candidates for novel associations including CSMD1 and a monostearyl alcohol triglyceride (m/z 781.7483, retention time (RT) 89.3 s); CALN1 and Tridecanol (m/z 283.2741, RT 27.6). A gene-level enrichment analysis using MAGMA revealed highly interconnected modules for dADP biosynthesis, sterol synthesis, and long-chain fatty acid transport in the gene-feature network.

CONCLUSION

Comprehensive profiling of plasma metabolome across age groups combined with genome-wide genotyping revealed a wide range of genetic influence on diverse chemical species and metabolic pathways. The developmental trajectory of a biological system is shaped by gene-environment interaction especially in early life. Therefore, continuous efforts on generating metabolomics data in diverse human tissue types across age groups are required to understand gene-environment interaction toward healthy aging trajectories.

Quantifying the role of transcript levels in mediating DNA methylation effects on complex traits and diseases

High-dimensional omics datasets provide valuable resources to determine the causal role of molecular traits in mediating the path from genotype to phenotype. Making use of molecular quantitative trait loci (QTL) and genome-wide association study (GWAS) summary statistics, we propose a multivariable Mendelian randomization (MVMR) framework to quantify the proportion of the impact of the DNA methylome (DNAm) on complex traits that is propagated through the assayed transcriptome. Evaluating 50 complex traits, we find that on average at least 28.3% (95% CI: [26.9%-29.8%]) of DNAm-to-trait effects are mediated through (typically multiple) transcripts in the cis-region. Several regulatory mechanisms are hypothesized, including methylation of the promoter probe cg10385390 (chr1:8'022'505) increasing the risk for inflammatory bowel disease by reducing PARK7 expression. The proposed integrative framework can be extended to other omics layers to identify causal molecular chains, providing a powerful tool to map and interpret GWAS signals.

Genetically personalised organ-specific metabolic models in health and disease

Understanding how genetic variants influence disease risk and complex traits (variant-to-function) is one of the major challenges in human genetics. Here we present a model-driven framework to leverage human genome-scale metabolic networks to define how genetic variants affect biochemical reaction fluxes across major human tissues, including skeletal muscle, adipose, liver, brain and heart. As proof of concept, we build personalised organ-specific metabolic flux models for 524,615 individuals of the INTERVAL and UK Biobank cohorts and perform a fluxome-wide association study (FWAS) to identify 4312 associations between personalised flux values and the concentration of metabolites in blood. Furthermore, we apply FWAS to identify 92 metabolic fluxes associated with the risk of developing coronary artery disease, many of which are linked to processes previously described to play in role in the disease. Our work demonstrates that genetically personalised metabolic models can elucidate the downstream effects of genetic variants on biochemical reactions involved in common human diseases.

Identification of 969 protein quantitative trait loci in an African American population with kidney disease attributed to hypertension

Investigations into the causal underpinnings of disease processes can be aided by the incorporation of genetic information. Genetic studies require populations varied in both ancestry and prevalent disease in order to optimize discovery and ensure generalizability of findings to the global population. Here, we report the genetic determinants of the serum proteome in 466 African Americans with chronic kidney disease attributed to hypertension from the richly phenotyped African American Study of Kidney Disease and Hypertension (AASK) study. Using the largest aptamer-based protein profiling platform to date (6,790 proteins or protein complexes), we identified 969 genetic associations with 900 unique proteins; including 52 novel cis (local) associations and 379 novel trans (distant) associations. The genetic effects of previously published cis-protein quantitative trait loci (pQTLs) were found to be highly reproducible, and we found evidence that our novel genetic signals colocalize with gene expression and disease processes. Many trans- pQTLs were found to reflect associations mediated by the circulating cis protein, and the common trans-pQTLs are enriched for processes involving extracellular vesicles, highlighting a plausible mechanism for distal regulation of the levels of secreted proteins. Thus, our study generates a valuable resource of genetic associations linking variants to protein levels and disease in an understudied patient population to inform future studies of drug targets and physiology.

Metabolic diversity in human populations and correlation with genetic and ancestral geographic distances

DNA polymorphic markers and self-defined ethnicity groupings are used to group individuals with shared ancient geographic ancestry. Here we studied whether ancestral relationships between individuals could be identified from metabolic screening data reported by the California newborn screening (NBS) program. NBS data includes 41 blood metabolites measured by tandem mass spectrometry from singleton babies in 17 parent-reported ethnicity groupings. Ethnicity-associated differences identified for 71% of NBS metabolites (29 of 41, Cohen's d &gt; 0.5) showed larger differences in blood levels of acylcarnitines than of amino acids (P &lt; 1e-4). A metabolic distance measure, developed to compare ethnic groupings based on metabolic differences, showed low positive correlation with genetic and ancient geographic distances between the groups' ancestral world populations. Several outlier group pairs were identified with larger genetic and smaller metabolic distances (Black versus White) or with smaller genetic and larger metabolic distances (Chinese versus Japanese) indicating the influence of genetic and of environmental factors on metabolism. Using machine learning, comparison of metabolic profiles between all pairs of ethnic groupings distinguished individuals with larger genetic distance (Black versus Chinese, AUC = 0.96), while genetically more similar individuals could not be separated metabolically (Hispanic versus Native American, AUC = 0.51). Additionally, we identified metabolites informative for inferring metabolic ancestry in individuals from genetically similar populations, which included biomarkers for inborn metabolic disorders (C10:1, C12:1, C3, C5OH, Leucine-Isoleucine). This work sheds new light on metabolic differences in healthy newborns in diverse populations, which could have implications for improving genetic disease screening.

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.

APOE alleles modulate associations of plasma metabolites with variants from multiple genes on chromosome 19q13.3

The APOE ε2, ε3, and ε4 alleles differentially impact various complex diseases and traits. We examined whether these alleles modulated associations of 94 single-nucleotide polymorphisms (SNPs) harbored by 26 genes in 19q13.3 region with 217 plasma metabolites using Framingham Heart Study data. The analyses were performed in the E2 (ε2ε2 or ε2ε3 genotype), E3 (ε3ε3 genotype), and E4 (ε3ε4 or ε4ε4 genotype) groups separately. We identified 31, 17, and 22 polymorphism-metabolite associations in the E2, E3, and E4 groups, respectively, at a false discovery rate P _FDR < 0.05. These entailed 51 and 19 associations with 20 lipid and 12 polar analytes. Contrasting the effect sizes between the analyzed groups showed 20 associations with group-specific effects at Bonferroni-adjusted P < 7.14E-04. Three associations with glutamic acid or dimethylglycine had significantly larger effects in the E2 than E3 group and 12 associations with triacylglycerol 56:5, lysophosphatidylethanolamines 16:0, 18:0, 20:4, or phosphatidylcholine 38:6 had significantly larger effects in the E2 than E4 group. Two associations with isocitrate or propionate and three associations with phosphatidylcholines 32:0, 32:1, or 34:0 had significantly larger effects in the E4 than E3 group. Nine of 70 SNP-metabolite associations identified in either E2, E3, or E4 groups attained P _FDR < 0.05 in the pooled sample of these groups. However, none of them were among the 20 group-specific associations. Consistent with the evolutionary history of the APOE alleles, plasma metabolites showed higher APOE-cluster-related variations in the E4 than E2 and E3 groups. Pathway enrichment mainly highlighted lipids and amino acids metabolism and citrate cycle, which can be differentially impacted by the APOE alleles. These novel findings expand insights into the genetic heterogeneity of plasma metabolites and highlight the importance of the APOE-allele-stratified genetic analyses of the APOE-related diseases and traits.

The Genetic Basis for the Increased Prevalence of Metabolic Syndrome among Post-Traumatic Stress Disorder Patients

Post-traumatic stress disorder (PTSD) is a highly debilitating psychiatric disorder that can be triggered by exposure to extreme trauma. Even if PTSD is primarily a psychiatric condition, it is also characterized by adverse somatic comorbidities. One illness commonly co-occurring with PTSD is Metabolic syndrome (MetS), which is defined by a set of health risk/resilience factors including obesity, elevated blood pressure, lower high-density lipoprotein cholesterol, higher low-density lipoprotein cholesterol, higher triglycerides, higher fasting blood glucose and insulin resistance. Here, phenotypic association between PTSD and components of MetS are tested on a military veteran cohort comprising chronic PTSD presentation (n = 310, 47% cases, 83% male). Consistent with previous observations, we found significant phenotypic correlation between the various components of MetS and PTSD severity scores. To examine if this observed symptom correlations stem from a shared genetic background, we conducted genetic correlation analysis using summary statistics data from large-scale genetic studies. Our results show robust positive genetic correlation between PTSD and MetS (r_g[SE] = 0.33 [0.056], p = 4.74E-09), and obesity-related components of MetS (r_g = 0.25, SE = 0.05, p = 6.4E-08). Prioritizing genomic regions with larger local genetic correlation implicate three significant loci. Overall, these findings show significant genetic overlap between PTSD and MetS, which may in part account for the markedly increased occurrence of MetS among PTSD patients.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

Background

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

Methods

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

Results

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

Conclusions

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

Funding

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

BACKGROUND

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

METHODS

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

RESULTS

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

CONCLUSIONS

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

FUNDING

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Multi-Omic analyses characterize the ceramide/sphingomyelin pathway as a therapeutic target in Alzheimer's disease

Dysregulation of sphingomyelin and ceramide metabolism have been implicated in Alzheimer's disease. Genome-wide and transcriptome-wide association studies have identified various genes and genetic variants in lipid metabolism that are associated with Alzheimer's disease. However, the molecular mechanisms of sphingomyelin and ceramide disruption remain to be determined. We focus on the sphingolipid pathway and carry out multi-omics analyses to identify central and peripheral metabolic changes in Alzheimer's patients, correlating them to imaging features. Our multi-omics approach is based on (a) 2114 human post-mortem brain transcriptomics to identify differentially expressed genes; (b) in silico metabolic flux analysis on context-specific metabolic networks identified differential reaction fluxes; (c) multimodal neuroimaging analysis on 1576 participants to associate genetic variants in sphingomyelin pathway with Alzheimer's disease pathogenesis; (d) plasma metabolomic and lipidomic analysis to identify associations of lipid species with dysregulation in Alzheimer's; and (e) metabolite genome-wide association studies to define receptors within the pathway as a potential drug target. We validate our hypothesis in amyloidogenic APP/PS1 mice and show prolonged exposure to fingolimod alleviated synaptic plasticity and cognitive impairment in mice. Our integrative multi-omics approach identifies potential targets in the sphingomyelin pathway and suggests modulators of S1P metabolism as possible candidates for Alzheimer's disease treatment.