Rare and common genetic determinants of metabolic individuality and their effects on human 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 imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Genetic associations with ratios between protein levels detect new pQTLs and reveal protein-protein interactions

Protein quantitative trait loci (pQTLs) are an invaluable source of information for drug target development because they provide genetic evidence to support protein function, suggest relationships between cis- and trans-associated proteins, and link proteins to disease endpoints. Using Olink proteomics data for 1,463 proteins measured in over 54,000 samples of the UK Biobank, we identified 4,248 associations with 2,821 ratios between protein levels (rQTLs). rQTLs were 7.6-fold enriched in known protein-protein interactions, suggesting that their ratios reflect biological links between the implicated proteins. Conducting a GWAS on ratios increased the number of discovered genetic signals by 24.7%. The approach can identify novel loci of clinical relevance, support causal gene identification, and reveal complex networks of interacting proteins. Taken together, our study adds significant value to the genetic insights that can be derived from the UKB proteomics data and motivates the wider use of ratios in large-scale GWAS.

Plasma metabolites and risk of seven cancers: a two-sample Mendelian randomization study among European descendants

BACKGROUND

While circulating metabolites have been increasingly linked to cancer risk, the causality underlying these associations remains largely uninterrogated.

METHODS

We conducted a comprehensive 2-sample Mendelian randomization (MR) study to evaluate the potential causal relationship between 913 plasma metabolites and the risk of seven cancers among European-ancestry individuals. Data on variant-metabolite associations were obtained from a genome-wide association study (GWAS) of plasma metabolites among 14,296 subjects. Data on variant-cancer associations were gathered from large-scale GWAS consortia for breast (N = 266,081), colorectal (N = 185,616), lung (N = 85,716), ovarian (N = 63,347), prostate (N = 140,306), renal cell (N = 31,190), and testicular germ cell (N = 28,135) cancers. MR analyses were performed with the inverse variance-weighted (IVW) method as the primary strategy to identify significant associations at Bonferroni-corrected P < 0.05 for each cancer type separately. Significant associations were subjected to additional scrutiny via weighted median MR, Egger regression, MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO), and reverse MR analyses. Replication analyses were performed using an independent dataset from a plasma metabolite GWAS including 8,129 participants of European ancestry.

RESULTS

We identified 94 significant associations, suggesting putative causal associations between 66 distinct plasma metabolites and the risk of seven cancers. Remarkably, 68.2% (45) of these metabolites were each associated with the risk of a specific cancer. Among the 66 metabolites, O-methylcatechol sulfate and 4-vinylphenol sulfate demonstrated the most pronounced positive and negative associations with cancer risk, respectively. Genetically proxied plasma levels of these two metabolites were significantly associated with the risk of lung cancer and renal cell cancer, with an odds ratio and 95% confidence interval of 2.81 (2.33-3.37) and 0.49 (0.40-0.61), respectively. None of these 94 associations was biased by weak instruments, horizontal pleiotropy, or reverse causation. Further, 64 of these 94 were eligible for replication analyses, and 54 (84.4%) showed P < 0.05 with association patterns consistent with those shown in primary analyses.

CONCLUSIONS

Our study unveils plausible causal relationships between 66 plasma metabolites and cancer risk, expanding our understanding of the role of circulating metabolites in cancer genetics and etiology. These findings hold promise for enhancing cancer risk assessment and prevention strategies, meriting further exploration.

Genetics unravels protein-metabolite relationships

Integrating molecular traits into genetic studies enhances our understanding of how DNA variation influences complex clinical and physiological phenotypes. In a recent article, Benson and colleagues apply this systems genetics approach with proteomics and metabolomics data in plasma from humans to identify and validate several previously unrecognized causal protein-metabolite associations.

Elucidating the role of blood metabolites on pancreatic cancer risk using two-sample Mendelian randomization analysis

Pancreatic ductal adenocarcinoma (PDAC) is an uncommon but highly fatal malignancy. Identifying causal metabolite biomarkers offers an opportunity to facilitate effective risk assessment strategies for PDAC. In this study, we performed a two-sample Mendelian randomization (MR) study to characterize the potential causal effects of metabolites in plasma on PDAC risk. Genetic instruments were determined for a total of 506 metabolites from one set of comprehensive genome-wide association studies (GWAS) involving 913 individuals of European ancestry from the INTERVAL/EPIC-Norfolk cohorts. Another set of genetic instruments was developed for 483 metabolites from an independent GWAS conducted with 8299 individuals of European ancestry from the Canadian Longitudinal Study on Aging (CLSA) cohort. We analyzed GWAS data of the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4), comprising 8275 PDAC cases and 6723 controls of European ancestry. The association of metabolites with PDAC risk was assessed using the inverse-variance weighted (IVW) method, and complemented with sensitivity analyses of MR-Egger and MR-PRESSO tests. Potential side effects of targeting the identified metabolites for PDAC intervention were further evaluated by a phenome-wide MR (Phe-MR) analysis. Forty-four unique metabolites were identified to be significantly associated with PDAC risk, of which four top-ranking metabolites (X: 12798, X: 11787, X: 11308 and X: 19141) showed replication evidence when using instruments developed from both two cohorts. Our results highlight novel blood metabolites related to PDAC risk, which may help prioritize metabolic features for PDAC mechanistic research and further evaluation of their potential role in PDAC risk assessment.

A Metabolome-Wide Mendelian Randomization Study Identifies Dysregulated Arachidonic Acid Synthesis as a Potential Causal Risk Factor for Bipolar Disorder

BACKGROUND

Bipolar disorder (BPD) is a debilitating mood disorder with an unclear etiology. A better understanding of the underlying pathophysiological mechanisms will help to identify novel targets for improved treatment options and prevention strategies. In this metabolome-wide Mendelian randomization study, we screened for metabolites that may have a causal role in BPD.

METHODS

We tested a total of 913 circulating metabolite exposures assessed in 14,296 Europeans using a mass spectrometry-based platform. For the BPD outcome, we used summary data from the largest and most recent genome-wide association study reported to date, including 41,917 BPD cases.

RESULTS

We identified 33 metabolites associated with BPD (padjusted < 5.48 × 10-5). Most of them were lipids, including arachidonic acid (β = -0.154, SE = 0.023, p = 3.30 × 10-11), a polyunsaturated omega-6 fatty acid, along with several complex lipids containing either an arachidonic or a linoleic fatty acid side chain. These associations did not extend to other closely related psychiatric disorders like schizophrenia or depression, although they may be involved in the regulation of lithium response. These lipid associations were driven by genetic variants within the FADS1/2/3 gene cluster, which is a robust BPD risk locus encoding a family of fatty acid desaturase enzymes that are responsible for catalyzing the conversion of linoleic acid into arachidonic acid. Statistical colocalization analyses indicated that 27 of the 33 metabolites shared the same genetic etiology with BPD at the FADS1/2/3 cluster, demonstrating that our findings are not confounded by linkage disequilibrium.

CONCLUSIONS

Overall, our findings support the notion that arachidonic acid and other polyunsaturated fatty acids may represent potential targets for BPD.

Regulation of kynurenine metabolism by blood donor genetics and biology impacts red cell hemolysis in vitro and in vivo

ABSTRACT

In the field of transfusion medicine, the clinical relevance of the metabolic markers of the red blood cell (RBC) storage lesion is incompletely understood. Here, we performed metabolomics of RBC units from 643 donors enrolled in the Recipient Epidemiology and Donor Evaluation Study, REDS RBC Omics. These units were tested on storage days 10, 23, and 42 for a total of 1929 samples and also characterized for end-of-storage hemolytic propensity after oxidative and osmotic insults. Our results indicate that the metabolic markers of the storage lesion poorly correlated with hemolytic propensity. In contrast, kynurenine was not affected by storage duration and was identified as the top predictor of osmotic fragility. RBC kynurenine levels were affected by donor age and body mass index and were reproducible within the same donor across multiple donations from 2 to 12 months apart. To delve into the genetic underpinnings of kynurenine levels in stored RBCs, we thus tested kynurenine levels in stored RBCs on day 42 from 13 091 donors from the REDS RBC Omics study, a population that was also genotyped for 879 000 single nucleotide polymorphisms. Through a metabolite quantitative trait loci analysis, we identified polymorphisms in SLC7A5, ATXN2, and a series of rate-limiting enzymes (eg, kynurenine monooxygenase, indoleamine 2,3-dioxygenase, and tryptophan dioxygenase) in the kynurenine pathway as critical factors affecting RBC kynurenine levels. By interrogating a donor-recipient linkage vein-to-vein database, we then report that SLC7A5 polymorphisms are also associated with changes in hemoglobin and bilirubin levels, suggestive of in vivo hemolysis in 4470 individuals who were critically ill and receiving single-unit transfusions.

The Relationship of Circulating Choline and Choline-Related Metabolite Levels with Health Outcomes: A Scoping Review of Genome-Wide Association Studies and Mendelian Randomization Studies

Choline is essential for proper liver, muscle, brain, lipid metabolism, cellular membrane composition, and repair. Understanding genetic determinants of circulating choline metabolites can help identify new determinants of choline metabolism, requirements, and their link to disease endpoints. We conducted a scoping review to identify studies assessing the association of genetic polymorphisms on circulating choline and choline-related metabolite concentrations and subsequent associations with health outcomes. This study follows the Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement scoping review extension. Literature was searched to September 28, 2022, in 4 databases: Embase, MEDLINE, Web of Science, and the Biological Science Index. Studies of any duration in humans were considered. Any genome-wide association study (GWAS) investigating genetic variant associations with circulating choline and/or choline-related metabolites and any Mendelian randomization (MR) study investigating the association of genetically predicted circulating choline and/or choline-related metabolites with any health outcome were considered. Qualitative evidence is presented in summary tables. From 1248 total reviewed articles, 53 were included (GWAS = 27; MR = 26). Forty-two circulating choline-related metabolites were tested in association with genetic variants in GWAS studies, primarily trimethylamine N-oxide, betaine, sphingomyelins, lysophosphatidylcholines, and phosphatidylcholines. MR studies investigated associations between 52 total unique choline metabolites and 66 unique health outcomes. Of these, 47 significant associations were reported between 16 metabolites (primarily choline, lysophosphatidylcholines, phosphatidylcholines, betaine, and sphingomyelins) and 27 health outcomes including cancer, cardiovascular, metabolic, bone, and brain-related outcomes. Some articles reported significant associations between multiple choline types and the same health outcome. Genetically predicted circulating choline and choline-related metabolite concentrations are associated with a wide variety of health outcomes. Further research is needed to assess how genetic variability influences choline metabolism and whether individuals with lower genetically predicted circulating choline and choline-related metabolite concentrations would benefit from a dietary intervention or supplementation.

IDSL.GOA: gene ontology analysis for interpreting metabolomic datasets

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2393 metabolic GO terms and associated 3144 genes, 1,492 EC annotations, and 2621 metabolites. IDSL.GOA analysis of a case study of older versus young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR < 0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/ .

Mendelian randomization with incomplete measurements on the exposure in the Hispanic Community Health Study/Study of Latinos

Mendelian randomization has been widely used to assess the causal effect of a heritable exposure variable on an outcome of interest, using genetic variants as instrumental variables. In practice, data on the exposure variable can be incomplete due to high cost of measurement and technical limits of detection. In this paper, we propose a valid and efficient method to handle both unmeasured and undetectable values of the exposure variable in one-sample Mendelian randomization analysis with individual-level data. We estimate the causal effect of the exposure variable on the outcome using maximum likelihood estimation and develop an expectation maximization algorithm for the computation of the estimator. Simulation studies show that the proposed method performs well in making inference on the causal effect. We apply our method to the Hispanic Community Health Study/Study of Latinos, a community-based prospective cohort study, and estimate the causal effect of several metabolites on phenotypes of interest.

IDSL.GOA: Gene Ontology Analysis for Interpreting Metabolomic datasets

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2,393 metabolic GO terms and associated 3,144 genes, 1,492 EC annotations, and 2,621 metabolites. IDSL.GOA analysis of a case study of older vs young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR <0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/.

Biological and Genetic Determinants of Glycolysis: Phosphofructokinase Isoforms Boost Energy Status of Stored Red Blood Cells and Transfusion Outcomes

Mature red blood cells (RBCs) lack mitochondria, and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo and during storage in vitro in the blood bank. Here we identify an association between blood donor age, sex, ethnicity and end of storage levels of glycolytic metabolites in 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study. Associations were also observed to ancestry-specific genetic polymorphisms in regions coding for phosphofructokinase 1, platelet - which we detected in mature RBCs - hexokinase 1, and for the ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice, and via multi-omics characterization of 1,929 samples from 643 RBC units during storage. ATP levels, breakdown and deamination into hypoxanthine were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 4,700 heterologous critically ill patients.

Highlights

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry-based genetic polymorphisms in PFKP and CD38/BST1 influence glycolysis in stored human and fresh or stored murine RBCs;PFKP is detected in pure mature RBCs and boosts glycolytic fluxes when ATP is low, such as in stored RBCs;ATP and hypoxanthine associate with RBC hemolysis in vitro and in vivo in mice and thousands of transfusion recipients.

Graphical abstract

Towards personalized genome-scale modeling of inborn errors of metabolism for systems medicine applications

Inborn errors of metabolism (IEMs) are a group of more than 1000 inherited diseases that are individually rare but have a cumulative global prevalence of 50 per 100,000 births. Recently, it has been recognized that like common diseases, patients with rare diseases can greatly vary in the manifestation and severity of symptoms. Here, we review omics-driven approaches that enable an integrated, holistic view of metabolic phenotypes in IEM patients. We focus on applications of Constraint-based Reconstruction and Analysis (COBRA), a widely used mechanistic systems biology approach, to model the effects of inherited diseases. Moreover, we review evidence that the gut microbiome is also altered in rare diseases. Finally, we outline an approach using personalized metabolic models of IEM patients for the prediction of biomarkers and tailored therapeutic or dietary interventions. Such applications could pave the way towards personalized medicine not just for common, but also for rare diseases.

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.

An interactive web application for exploring human plasma and fibroblast metabolomics data from patients with inborn errors of metabolism

Metabolomic profiling is instrumental in understanding the systemic and cellular impact of inborn errors of metabolism (IEMs), monogenic disorders caused by pathogenic genomic variants in genes involved in metabolism. This study encompasses untargeted metabolomics analysis of plasma from 474 individuals and fibroblasts from 67 subjects, incorporating healthy controls, patients with 65 different monogenic diseases, and numerous undiagnosed cases. We introduce a web application designed for the in-depth exploration of this extensive metabolomics database. The application offers a user-friendly interface for data review, download, and detailed analysis of metabolic deviations linked to IEMs at the level of individual patients or groups of patients with the same diagnosis. It also provides interactive tools for investigating metabolic relationships and offers comparative analyses of plasma and fibroblast profiles. This tool emphasizes the metabolic interplay within and across biological matrices, enriching our understanding of metabolic regulation in health and disease. As a resource, the application provides broad utility in research, offering novel insights into metabolic pathways and their alterations in various disorders.

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.

A phenome-wide scan reveals convergence of common and rare variant associations

BACKGROUND

Common and rare variants contribute to the etiology of complex traits. However, the extent to which the phenotypic effects of common and rare variants involve shared molecular mediators remains poorly understood. The question is essential to the basic and translational goals of the science of genomics, with critical basic-science, methodological, and clinical consequences.

METHODS

Leveraging the latest release of whole-exome sequencing (WES, for rare variants) and genome-wide association study (GWAS, for common variants) data from the UK Biobank, we developed a metric, the COmmon variant and RAre variant Convergence (CORAC) signature, to quantify the convergence for a broad range of complex traits. We characterized the relationship between CORAC and effective sample size across phenome-wide association studies.

RESULTS

We found that the signature is positively correlated with effective sample size (Spearman ρ = 0.594, P < 2.2e - 16), indicating increased functional convergence of trait-associated genetic variation, across the allele frequency spectrum, with increased power. Sensitivity analyses, including accounting for heteroskedasticity and varying the number of detected association signals, further strengthened the validity of the finding. In addition, consistent with empirical data, extensive simulations showed that negative selection, in line with enhancing polygenicity, has a dampening effect on the convergence signature. Methodologically, leveraging the convergence leads to enhanced association analysis.

CONCLUSIONS

The presented framework for the convergence signature has important implications for fine-mapping strategies and drug discovery efforts. In addition, our study provides a blueprint for the expectation from future large-scale whole-genome sequencing (WGS)/WES and sheds methodological light on post-GWAS studies.

Metabolic dynamics and prediction of sFGR and adverse fetal outcomes: a prospective longitudinal cohort study

BACKGROUND

Selective fetal growth restriction (sFGR) is an extreme complication that significantly increases the risk of perinatal mortality and long-term adverse neurological outcomes in offspring, affecting approximately 15% of monochorionic diamniotic (MCDA) twin pregnancies. The lack of longitudinal cohort studies hinders the early prediction and intervention of sFGR.

METHODS

We constructed a prospective longitudinal cohort study of sFGR, and quantified 25 key metabolites in 337 samples from maternal plasma in the first, second, and third trimester and from cord plasma. In particular, our study examined fetal growth and brain injury data from ultrasonography and used the Ages and Stages Questionnaire-third edition subscale (ASQ-3) to evaluate the long-term neurocognitive behavioral development of infants aged 2-3 years. Furthermore, we correlated metabolite levels with ultrasound data, including physical development and brain injury indicators, and ASQ-3 data using Spearman's-based correlation tests. In addition, special combinations of differential metabolites were used to construct predictive models for the occurrence of sFGR and fetal brain injury.

RESULTS

Our findings revealed various dynamic patterns for these metabolites during pregnancy and a maximum of differential metabolites between sFGR and MCDA in the second trimester (n = 8). The combination of L-phenylalanine, L-leucine, and L-isoleucine in the second trimester, which were closely related to fetal growth indicators, was highly predictive of sFGR occurrence (area under the curve [AUC]: 0.878). The combination of L-serine, L-histidine, and L-arginine in the first trimester and creatinine in the second trimester was correlated with long-term neurocognitive behavioral development and showed the capacity to identify fetal brain injury with high accuracy (AUC: 0.94).

CONCLUSIONS

The performance of maternal plasma metabolites from the first and second trimester is superior to those from the third trimester and cord plasma in discerning sFGR and fetal brain injury. These metabolites may serve as useful biomarkers for early prediction and promising targets for early intervention in clinical settings.

The Metabolome-Wide Signature of Major Depressive Disorder

Major Depressive Disorder (MDD) is an often-chronic condition with substantial molecular alterations and pathway dysregulations involved. Single metabolite, pathway and targeted metabolomics platforms have indeed revealed several metabolic alterations in depression including energy metabolism, neurotransmission and lipid metabolism. More comprehensive coverage of the metabolome is needed to further specify metabolic dysregulation in depression and reveal previously untargeted mechanisms. Here we measured 820 metabolites using the metabolome-wide Metabolon platform in 2770 subjects from a large Dutch clinical cohort with extensive depression clinical phenotyping (1101 current MDD, 868 remitted MDD, 801 healthy controls) at baseline and 1805 subjects at 6-year follow up (327 current MDD, 1045 remitted MDD, 433 healthy controls). MDD diagnosis was based on DSM-IV psychiatric interviews. Depression severity was measured with the Inventory of Depressive Symptomatology self-report. Associations between metabolites and MDD status and depression severity were assessed at baseline and at the 6-year follow-up. Metabolites consistently associated with MDD status or depression severity on both occasions were examined in Mendelian randomization (MR) analysis using metabolite (N=14,000) and MDD (N=800,000) GWAS results. At baseline, 139 and 126 metabolites were associated with current MDD status and depression severity, respectively, with 79 overlapping metabolites. Six years later, 34 out of the 79 metabolite associations were subsequently replicated. Downregulated metabolites were enriched with long-chain monounsaturated (P=6.7e-07) and saturated (P=3.2e-05) fatty acids and upregulated metabolites with lysophospholipids (P=3.4e-4). Adding BMI to the models changed results only marginally. MR analyses showed that genetically-predicted higher levels of the lysophospholipid 1-linoleoyl-GPE (18:2) were associated with greater risk of depression. The identified metabolome-wide profile of depression (severity) indicated altered lipid metabolism with downregulation of long-chain fatty acids and upregulation of lysophospholipids, for which causal involvement was suggested using genetic tools. This metabolomics signature offers a window on depression pathophysiology and a potential access point for the development of novel therapeutic approaches.

The post-finasteride syndrome: possible etiological mechanisms and symptoms

Finasteride and dutasteride, synthetic 5α-reductase inhibitors (5ARIs) are recommended in many guidelines for the treatment of benign prostatic hyperplasia/lower urinary tract symptoms and alopecia despite a variety of side effects like sexual, neurological, psychiatric, endocrinological, metabolic and ophthalmological dysfunctions and the increased incidence of high grade prostate cancer. The sexual side effects are common during the use of the drug but in a small subgroup of patients, they can persist after stopping the drug. This so-called post-finasteride syndrome has serious implications for the quality of life without a clear etiology or therapy. Three types of 5α-reductases are present in many organs in- and outside the brain where they can be blocked by the two 5ARIs. There is increasing evidence that 5ARIs not only inhibit the conversion of testosterone to 5α-dihydrotestosterone (DHT) in the prostate and the scalp but also in many other tissues. The lipophilic 5ARIs can pass the blood-brain barrier and might block many other neurosteroids in the brain with changes in the neurochemistry and impaired neurogenesis. Further research and therapeutic innovations are urgently needed that might cure or relieve these side effects. More awareness is needed for physicians to outweigh these health risks against the benefits of 5ARIs.

The HuMet Repository: Watching human metabolism at work

The human metabolism constantly responds to stimuli such as food intake, fasting, exercise, and stress, triggering adaptive biochemical processes across multiple metabolic pathways. To understand the role of these processes and disruptions thereof in health and disease, detailed documentation of healthy metabolic responses is needed but still scarce on a time-resolved metabolome-wide level. Here, we present the HuMet Repository, a web-based resource for exploring dynamic metabolic responses to six physiological challenges (exercise, 36 h fasting, oral glucose and lipid loads, mixed meal, cold stress) in healthy subjects. For building this resource, we integrated existing and newly derived metabolomics data measured in blood, urine, and breath samples of 15 young healthy men at up to 56 time points during the six highly standardized challenge tests conducted over four days. The data comprise 1.1 million data points acquired on multiple platforms with temporal profiles of 2,656 metabolites from a broad range of biochemical pathways. By embedding the dataset into an interactive web application, we enable users to easily access, search, filter, analyze, and visualize the time-resolved metabolomic readouts and derived results. Users can put metabolites into their larger context by identifying metabolites with similar trajectories or by visualizing metabolites within holistic metabolic networks to pinpoint pathways of interest. In three showcases, we outline the value of the repository for gaining biological insights and generating hypotheses by analyzing the wash-out of dietary markers, the complementarity of metabolomics platforms in dynamic versus cross-sectional data, and similarities and differences in systemic metabolic responses across challenges. With its comprehensive collection of time-resolved metabolomics data, the HuMet Repository, freely accessible at https://humet.org/, is a reference for normal, healthy responses to metabolic challenges in young males. It will enable researchers with and without computational expertise, to flexibly query the data for their own research into the dynamics of human metabolism.

Beliefs and counseling practices among dermatologists regarding sexual and other adverse effects of finasteride

Finasteride may cause low libido and erectile dysfunction and the product label of finasteride also includes post-marketing reactions of sexual dysfunction that continued after discontinuation of treatment, as well as male infertility and depression. The aim of this study was to evaluate the beliefs and counseling practices among dermatologists regarding adverse effects of finasteride. Anonymous paper surveys were personally distributed to 122 attendees at two annual major dermatology meetings. The participation rate was 82% with 47% women and 77% residents of the United States. 51% of respondents believed that finasteride could cause sexual side effects and 18% believed that it could cause persistent sexual side effects. Fewer than a quarter believed that finasteride could cause depression or lower sperm counts. When initiating finasteride, 69% of respondents counseled at least half of their patients about potential sexual side effects with 52% for persistent sexual side effects and 30% for depression. This study identifies the need for greater awareness of the potential adverse effects of finasteride and identifies opportunities for improvement in counseling practices that reflect finasteride's product labeling.

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.

A Multi-Omic Mosaic Model of Acetaminophen Induced Alanine Aminotransferase Elevation

BACKGROUND

Acetaminophen (APAP) is the most common cause liver injury following alcohol in US patients. Predicting liver injury and subsequent hepatic regeneration in patients taking therapeutic doses of APAP may be possible using new 'omic methods such as metabolomics and genomics. Multi'omic techniques increase our ability to find new mechanisms of injury and regeneration.

METHODS

We used metabolomic and genomic data from a randomized controlled trial of patients administered 4 g of APAP per day for 14 days or longer with blood samples obtained at 0 (baseline), 4, 7, 10, 13 and 16 days. We used the highest ALT as the clinical outcome to be predicted in our integrated analysis. We used penalized regression to model the relationship between genetic variants and day 0 metabolite level, and then performed a metabolite-wide colocalization scan to associate the genetically regulated component of metabolite expression with ALT elevation. Genome-wide association study (GWAS) analyses were conducted for ALT elevation and metabolite level using linear regression, with age, sex, and the first five principal components included as covariates. Colocalization was tested via a weighted sum test.

RESULTS

Out of the 164 metabolites modeled, 120 met the criteria for predictive accuracy and were retained for genetic analyses. After genomic examination, eight metabolites were found to be under genetic control and predictive of ALT elevation due to therapeutic acetaminophen. The metabolites were: 3-oxalomalate, allantoate, diphosphate, L-carnitine, L-proline, maltose, and ornithine. These genes are important in the tricarboxylic acid cycle (TCA), urea breakdown pathway, glutathione production, mitochondrial energy production, and maltose metabolism.

CONCLUSIONS

This multi'omic approach can be used to integrate metabolomic and genomic data allowing identification of genes that control downstream metabolites. These findings confirm prior work that have identified mitochondrial energy production as critical to APAP induced liver injury and have confirmed our prior work that demonstrate the importance of the urea cycle in therapeutic APAP liver injury.

Metabolome-wide Mendelian randomization characterizes heterogeneous and shared causal effects of metabolites on human health

Metabolites are small molecules that are useful for estimating disease risk and elucidating disease biology. Nevertheless, their causal effects on human diseases have not been evaluated comprehensively. We performed two-sample Mendelian randomization to systematically infer the causal effects of 1,099 plasma metabolites measured in 6,136 Finnish men from the METSIM study on risk of 2,099 binary disease endpoints measured in 309,154 Finnish individuals from FinnGen. We identified evidence for 282 causal effects of 70 metabolites on 183 disease endpoints (FDR<1%). We found 25 metabolites with potential causal effects across multiple disease domains, including ascorbic acid 2-sulfate affecting 26 disease endpoints in 12 disease domains. Our study suggests that N-acetyl-2-aminooctanoate and glycocholenate sulfate affect risk of atrial fibrillation through two distinct metabolic pathways and that N-methylpipecolate may mediate the causal effect of N6, N6-dimethyllysine on anxious personality disorder. This study highlights the broad causal impact of plasma metabolites and widespread metabolic connections across diseases.

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.

Red Blood Cell Omics and Machine Learning in Transfusion Medicine: Singularity Is Near

Background

Blood transfusion is a life-saving intervention for millions of recipients worldwide. Over the last 15 years, the advent of high-throughput, affordable omics technologies - including genomics, proteomics, lipidomics, and metabolomics - has allowed transfusion medicine to revisit the biology of blood donors, stored blood products, and transfusion recipients.

Summary

Omics approaches have shed light on the genetic and non-genetic factors (environmental or other exposures) impacting the quality of stored blood products and efficacy of transfusion events, based on the current Food and Drug Administration guidelines (e.g., hemolysis and post-transfusion recovery for stored red blood cells). As a treasure trove of data accumulates, the implementation of machine learning approaches promises to revolutionize the field of transfusion medicine, not only by advancing basic science. Indeed, computational strategies have already been used to perform high-content screenings of red blood cell morphology in microfluidic devices, generate in silico models of erythrocyte membrane to predict deformability and bending rigidity, or design systems biology maps of the red blood cell metabolome to drive the development of novel storage additives.

Key Message

In the near future, high-throughput testing of donor genomes via precision transfusion medicine arrays and metabolomics of all donated products will be able to inform the development and implementation of machine learning strategies that match, from vein to vein, donors, optimal processing strategies (additives, shelf life), and recipients, realizing the promise of personalized transfusion medicine.

The identification of metabolism-related subtypes and potential treatments for idiopathic pulmonary fibrosis

Background: Idiopathic pulmonary fibrosis (IPF) is caused by aberrant repair because of alveolar epithelial injury and can only be effectively treated with several compounds. Several metabolism-related biomolecular processes were found to be involved in IPF. We aimed to identify IPF subtypes based on metabolism-related pathways and explore potential drugs for each subtype. Methods: Gene profiles and clinical information were obtained from the Gene Expression Omnibus (GEO) database (GSE70867 and GSE93606). The enrichment scores for 41 metabolism-related pathways, immune cells, and immune pathways were calculated using the Gene Set Variation Analysis (GSVA) package. The ConsensusClusterPlus package was used to cluster samples. Novel modules and hub genes were identified using weighted correlation network analysis (WGCNA). Receiver operating characteristic (ROC) and calibration curves were plotted, and decision curve analysis (DCA) were performed to evaluate the model in the training and validation cohorts. A connectivity map was used as a drug probe. Results: Two subtypes with significant differences in prognosis were identified based on the metabolism-related pathways. Subtype C1 had a poor prognosis, low metabolic levels, and a unique immune signature. CDS2, LCLAT1, GPD1L, AGPAT1, ALDH3A1, LAP3, ADH5, AHCYL2, and MDH1 were used to distinguish between the two subtypes. Finally, subtype-specific drugs, which can potentially treat IPF, were identified. Conclusion: The aberrant activation of metabolism-related pathways contributes to differential prognoses in patients with IPF. Collectively, our findings provide novel mechanistic insights into subtyping IPF based on the metabolism-related pathway and potential treatments, which would help clinicians provide subtype-specific individualized therapeutic management to patients.

An atlas of genetic scores to predict multi-omic traits

The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables highly cost-effective and powerful analyses for studies that do not have multi-omics1. Here we examine a large cohort (the INTERVAL study2; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank3 to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.

Suicidal risk associated with finasteride versus dutasteride among men treated for benign prostatic hyperplasia: nationwide cohort study

Finasteride, a 5α-reductase inhibitor used in benign prostatic hyperplasia and androgenetic alopecia, has been associated with an increased suicidal risk, whereas it is unclear whether such risk is similar to that for another 5α-reductase inhibitor, dutasteride. We aimed to assess the risk of suicidal behaviours with finasteride relative to dutasteride. A nationwide cohort study was conducted using the French National Health Data System (SNDS). Men aged 50 years or older initiating finasteride 5 mg or dutasteride 0.5 mg in France between 01-01-2012 and 30-06-2016 were included and followed until outcome (suicide death identified from death certificate or self-harm hospitalisation), treatment discontinuation or switch, death, or 31-12-2016. Self-harm by violent means or resulting in admission to an intensive care unit were also examined. Cox proportional hazards models controlled for age and psychiatric and non-psychiatric conditions by inverse probability of treatment weighting (IPTW). Analyses were stratified according to psychiatric history. The study compared 69,786 finasteride new users to 217,577 dutasteride new users (median age: 72.0 years [Q1-Q3 = 64.5-80.2] vs. 71.1 [Q1-Q3 = 65.0-79.2]). During follow-up, 18 suicide deaths (0.57/1000 person-years) and 34 self-harm hospitalisations (1.08/1000) occurred among finasteride users versus 47 deaths (0.43/1000) and 87 hospitalisations (0.79/1000) among dutasteride users. Overall, finasteride was not associated with an increased risk of any suicidal outcome (IPTW-adjusted Hazard Ratio = 1.21 [95% Confidence Interval  .87-1.67]), suicide death or self-harm hospitalisation. However, among individuals with a history of mood disorders, finasteride was associated with an increased risk of any suicidal outcome (25 versus 46 events; HR = 1.64 [95% CI 1.00-2.68]), suicide death (8 versus 10 events; HR = 2.71 [95% CI 1.07-6.91]), self-harm by violent means (6 versus 6 events; HR = 3.11 [95% CI 1.01-9.61]), and self-harm with admission to an intensive care unit (7 versus 5 events; HR = 3.97 [95% CI 1.26-12.5]). None of these risks was significantly increased among individuals without a psychiatric history. These findings do not support an increased risk of suicide with finasteride used in the treatment of benign prostatic hyperplasia. However, an increased risk cannot be excluded among men with a history of mood disorder, but this result based on a limited number of events should be interpreted with caution.

Metabolic individuality: Limitations, challenges, and potential for clinical utility

In Nature Medicine, Surendran and colleagues recently reported the analysis of human plasma metabolomic data for 913 metabolites in ∼20,000 individuals, identifying 2,599 metabolite-genetic variant associations and >400 metabolite signatures comprised of jointly regulated metabolites. This extensive atlas of variant-metabolite relationships reveals novel genomic mechanisms driving metabolic phenotypes.

Identification of PCSK9-like human gene knockouts using metabolomics, proteomics, and whole-genome sequencing in a consanguineous population

Natural human knockouts of genes associated with desirable outcomes, such as PCSK9 with low levels of LDL-cholesterol, can lead to the discovery of new drug targets and treatments. Rare loss-of-function variants are more likely to be found in the homozygous state in consanguineous populations, and deep molecular phenotyping of blood samples from homozygous carriers can help to discriminate between silent and functional variants. Here, we combined whole-genome sequencing with proteomics and metabolomics for 2,935 individuals from the Qatar Biobank (QBB) to evaluate the power of this approach for finding genes of clinical and pharmaceutical interest. As proof-of-concept, we identified a homozygous carrier of a very rare PCSK9 variant with extremely low circulating PCSK9 levels and low LDL. Our study demonstrates that the chances of finding such variants are about 168 times higher in QBB compared with GnomAD and emphasizes the potential of consanguineous populations for drug discovery.