RaMP-DB 2.0: a renovated knowledgebase for deriving biological and chemical insight from metabolites, proteins, and genes

SmartGraph API: Programmatic Knowledge Mining in Network-Pharmacology Setting

The recent SmartGraph platform facilitates the execution of complex drug-discovery workflows with ease in the network-pharmacology paradigm. However, at the time of its publication we identified the need for the development of an Application Programming Interface (API) that could promote biomedical data integration and hypothesis generation in an automated manner. This need was magnified at the time of the COVID-19 pandemic. This study addresses the absence of such an API. Accordingly, most functionalities of the original platform were implemented within the SmartGraph API. We demonstrate that by using the API it is possible to transform the original semiautomated workflow behind the Neo4COVID19 database to a fully automated one. The availability of the SmartGraph API lends a significant improvement to the programmatic integration of network-pharmacology-oriented knowledge graphs and analytics, as well as predictive functionalities and workflows.

Introducing "Identification Probability" for Automated and Transferable Assessment of Metabolite Identification Confidence in Metabolomics and Related Studies

Methods for assessing compound identification confidence in metabolomics and related studies have been debated and actively researched for the past two decades. The earliest effort in 2007 focused primarily on mass spectrometry and nuclear magnetic resonance spectroscopy and resulted in four recommended levels of metabolite identification confidence─the Metabolite Standards Initiative (MSI) Levels. In 2014, the original MSI Levels were expanded to five levels (including two sublevels) to facilitate communication of compound identification confidence in high resolution mass spectrometry studies. Further refinement in identification levels have occurred, for example to accommodate use of ion mobility spectrometry in metabolomics workflows, and alternate approaches to communicate compound identification confidence also have been developed based on identification points schema. However, neither qualitative levels of identification confidence nor quantitative scoring systems address the degree of ambiguity in compound identifications in the context of the chemical space being considered. Neither are they easily automated nor transferable between analytical platforms. In this perspective, we propose that the metabolomics and related communities consider identification probability as an approach for automated and transferable assessment of compound identification and ambiguity in metabolomics and related studies. Identification probability is defined simply as 1/N, where N is the number of compounds in a database that matches an experimentally measured molecule within user-defined measurement precision(s), for example mass measurement or retention time accuracy, etc. We demonstrate the utility of identification probability in an in silico analysis of multiproperty reference libraries constructed from a subset of the Human Metabolome Database and computational property predictions, provide guidance to the community in transparent implementation of the concept, and invite the community to further evaluate this concept in parallel with their current preferred methods for assessing metabolite identification confidence.

Pairing omics to decode the diversity of plant specialized metabolism

Plants have evolved complex bouquets of specialized natural products that are utilized in medicine, agriculture, and industry. Untargeted natural product discovery has benefitted from growing plant omics data resources. Yet, plant genome complexity limits the identification and curation of biosynthetic pathways via single omics. Pairing multi-omics types within experiments provides multiple layers of evidence for biosynthetic pathway mining. The extraction of paired biological information facilitates connecting genes to transcripts and metabolites, especially when captured across time points, conditions and chemotypes. Experimental design requires specific adaptations to enable effective paired-omics analysis. Ultimately, metadata standards are required to support the integration of paired and unpaired public datasets and to accelerate collaborative efforts for natural product discovery in the plant research community.

Pharmacometabolomics Approach to Explore Pharmacokinetic Variation and Clinical Characteristics of a Single Dose of Desvenlafaxine in Healthy Volunteers

This study investigated the effects of a single dose of desvenlafaxine via oral administration on the pharmacokinetic parameters and clinical and laboratory characteristics in healthy volunteers using a pharmacometabolomics approach. In order to optimize desvenlafaxine's therapeutic use and minimize potential adverse effects, this knowledge is essential. Methods: Thirty-five healthy volunteers were enrolled after a health trial and received a single dose of desvenlafaxine (Pristiq®, 100 mg). First, liquid chromatography coupled to tandem mass spectrometry was used to determine the main pharmacokinetic parameters. Next, ultra-performance liquid chromatography-quadrupole time-of-flight mass spectrometry was used to identify plasma metabolites with different relative abundances in the metabolome at pre-dose and when the desvenlafaxine peak plasma concentration was reached (pre-dose vs. post-dose). Results: Correlations were observed between metabolomic profiles, such as tyrosine, sphingosine 1-phosphate, and pharmacokinetic parameters, as well as acetoacetic acid and uridine diphosphate glucose associated with clinical characteristics. Our findings suggest that desvenlafaxine may have a broader effect than previously thought by acting on the proteins responsible for the transport of various molecules at the cellular level, such as the solute carrier SLC and adenosine triphosphate synthase binding cassette ABC transporters. Both of these molecules have been associated with PK parameters and adverse events in our study. Conclusions: This altered transporter activity may be related to the reported side effects of desvenlafaxine, such as changes in blood pressure and liver function. This finding may be part of the explanation as to why people respond differently to the drug.

Modulation of anxiety-like behavior in galactooligosaccharide-fed mice: A potential role for bacterial tryptophan metabolites and reduced microglial reactivity

Prebiotic galactooligosaccharides (GOS) reduce anxiety-like behaviors in mice and humans. However, the biological pathways behind these behavioral changes are not well understood. To begin to study these pathways, we utilized C57BL/6 mice that were fed a standard diet with or without GOS supplementation for 3 weeks prior to testing on the open field. After behavioral testing, colonic contents and serum were collected for bacteriome (16S rRNA gene sequencing, colonic contents only) and metabolome (UPLC-MS, colonic contents and serum data) analyses. As expected, GOS significantly reduced anxiety-like behavior (i.e., increased time in the center) and decreased cytokine gene expression (Tnfa and Ccl2) in the prefrontal cortex. Notably, time in the center of the open field was significantly correlated with serum methyl-indole-3-acetic acid (methyl-IAA). This metabolite is a methylated form of indole-3-acetic acid (IAA) that is derived from bacterial metabolism of tryptophan. Sequencing analyses showed that GOS significantly increased Lachnospiraceae UCG006 and Akkermansia; these taxa are known to metabolize both GOS and tryptophan. To determine the extent to which methyl-IAA can affect anxiety-like behavior, mice were intraperitoneally injected with methyl-IAA. Mice given methyl-IAA had a reduction in anxiety-like behavior in the open field, along with lower Tnfa in the prefrontal cortex. Methyl-IAA was also found to reduce TNF-α (as well as CCL2) production by LPS-stimulated BV2 microglia. Together, these data support a novel pathway through which GOS reduces anxiety-like behaviors in mice and suggests that the bacterial metabolite methyl-IAA reduces microglial cytokine and chemokine production, which in turn reduces anxiety-like behavior.

Introducing 'identification probability' for automated and transferable assessment of metabolite identification confidence in metabolomics and related studies

Methods for assessing compound identification confidence in metabolomics and related studies have been debated and actively researched for the past two decades. The earliest effort in 2007 focused primarily on mass spectrometry and nuclear magnetic resonance spectroscopy and resulted in four recommended levels of metabolite identification confidence - the Metabolite Standards Initiative (MSI) Levels. In 2014, the original MSI Levels were expanded to five levels (including two sublevels) to facilitate communication of compound identification confidence in high resolution mass spectrometry studies. Further refinement in identification levels have occurred, for example to accommodate use of ion mobility spectrometry in metabolomics workflows, and alternate approaches to communicate compound identification confidence also have been developed based on identification points schema. However, neither qualitative levels of identification confidence nor quantitative scoring systems address the degree of ambiguity in compound identifications in context of the chemical space being considered, are easily automated, or are transferable between analytical platforms. In this perspective, we propose that the metabolomics and related communities consider identification probability as an approach for automated and transferable assessment of compound identification and ambiguity in metabolomics and related studies. Identification probability is defined simply as 1/N, where N is the number of compounds in a reference library or chemical space that match to an experimentally measured molecule within user-defined measurement precision(s), for example mass measurement or retention time accuracy, etc. We demonstrate the utility of identification probability in an in silico analysis of multi-property reference libraries constructed from the Human Metabolome Database and computational property predictions, provide guidance to the community in transparent implementation of the concept, and invite the community to further evaluate this concept in parallel with their current preferred methods for assessing metabolite identification confidence.

Metabolomic Profiling of Adipose Tissue in Type 2 Diabetes: Associations with Obesity and Insulin Resistance

The global prevalence of Type 2 Diabetes (T2D) poses significant public health challenges due to its associated severe complications. Insulin resistance is central to T2D pathophysiology, particularly affecting adipose tissue function. This cross-sectional observational study investigates metabolic alterations in subcutaneous adipose tissue (SAT) associated with T2D to identify potential therapeutic targets. We conducted a comprehensive metabolomic analysis of SAT from 40 participants (20 T2D, 20 ND-T2D), matched for sex, age, and BMI (Body Mass Index). Metabolite quantification was performed using GC/MS and LC/MS/MS platforms. Correlation analyses were conducted to explore associations between metabolites and clinical parameters. We identified 378 metabolites, including significant elevations in TCA cycle (tricarboxylic acid cycle) intermediates, branched-chain amino acids (BCAAs), and carbohydrates, and a significant reduction in the nucleotide-related metabolites in T2D subjects compared to those without T2D. Obesity exacerbated these alterations, particularly in amino acid metabolism. Adipocyte size negatively correlated with BCAAs, while adipocyte glucose uptake positively correlated with unsaturated fatty acids and glycerophospholipids. Our findings reveal distinct metabolic dysregulation in adipose tissue in T2D, particularly in energy metabolism, suggesting potential therapeutic targets for improving insulin sensitivity and metabolic health. Future studies should validate these findings in larger cohorts and explore underlying mechanisms to develop targeted interventions.

WebGestalt 2024: faster gene set analysis and new support for metabolomics and multi-omics

Enrichment analysis, crucial for interpreting genomic, transcriptomic, and proteomic data, is expanding into metabolomics. Furthermore, there is a rising demand for integrated enrichment analysis that combines data from different studies and omics platforms, as seen in meta-analysis and multi-omics research. To address these growing needs, we have updated WebGestalt to include enrichment analysis capabilities for both metabolites and multiple input lists of analytes. We have also significantly increased analysis speed, revamped the user interface, and introduced new pathway visualizations to accommodate these updates. Notably, the adoption of a Rust backend reduced gene set enrichment analysis time by 95% from 270.64 to 12.41 s and network topology-based analysis by 89% from 159.59 to 17.31 s in our evaluation. This performance improvement is also accessible in both the R package and a newly introduced Python package. Additionally, we have updated the data in the WebGestalt database to reflect the current status of each source and have expanded our collection of pathways, networks, and gene signatures. The 2024 WebGestalt update represents a significant leap forward, offering new support for metabolomics, streamlined multi-omics analysis capabilities, and remarkable performance enhancements. Discover these updates and more at https://www.webgestalt.org.

MetaboAnalyst 6.0: towards a unified platform for metabolomics data processing, analysis and interpretation

We introduce MetaboAnalyst version 6.0 as a unified platform for processing, analyzing, and interpreting data from targeted as well as untargeted metabolomics studies using liquid chromatography - mass spectrometry (LC-MS). The two main objectives in developing version 6.0 are to support tandem MS (MS2) data processing and annotation, as well as to support the analysis of data from exposomics studies and related experiments. Key features of MetaboAnalyst 6.0 include: (i) a significantly enhanced Spectra Processing module with support for MS2 data and the asari algorithm; (ii) a MS2 Peak Annotation module based on comprehensive MS2 reference databases with fragment-level annotation; (iii) a new Statistical Analysis module dedicated for handling complex study design with multiple factors or phenotypic descriptors; (iv) a Causal Analysis module for estimating metabolite - phenotype causal relations based on two-sample Mendelian randomization, and (v) a Dose-Response Analysis module for benchmark dose calculations. In addition, we have also improved MetaboAnalyst's visualization functions, updated its compound database and metabolite sets, and significantly expanded its pathway analysis support to around 130 species. MetaboAnalyst 6.0 is freely available at https://www.metaboanalyst.ca.

Consistent Multi-Omic Relationships Uncover Molecular Basis of Pediatric Asthma IgE Regulation

Serum total immunoglobulin E levels (total IgE) capture the state of the immune system in relation to allergic sensitization. High levels are associated with airway obstruction and poor clinical outcomes in pediatric asthma. Inconsistent patient response to anti-IgE therapies motivates discovery of molecular mechanisms underlying serum IgE level differences in children with asthma. To uncover these mechanisms using complementary metabolomic and transcriptomic data, abundance levels of 529 named metabolites and expression levels of 22,772 genes were measured among children with asthma in the Childhood Asthma Management Program (CAMP, N=564) and the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS, N=309) via the TOPMed initiative. Gene-metabolite associations dependent on IgE were identified within each cohort using multivariate linear models and were interpreted in a biochemical context using network topology, pathway and chemical enrichment, and representation within reactions. A total of 1,617 total IgE-dependent gene-metabolite associations from GACRS and 29,885 from CAMP met significance cutoffs. Of these, glycine and guanidinoacetic acid (GAA) were associated with the most genes in both cohorts, and the associations represented reactions central to glycine, serine, and threonine metabolism and arginine and proline metabolism. Pathway and chemical enrichment analysis further highlighted additional related pathways of interest. The results of this study suggest that GAA may modulate total IgE levels in two independent pediatric asthma cohorts with different characteristics, supporting the use of L-Arginine as a potential therapeutic for asthma exacerbation. Other potentially new targetable pathways are also uncovered.

The Plight of the Metabolite: Oxidative Stress and Tear Film Destabilisation Evident in Ocular Allergy Sufferers across Seasons in Victoria, Australia

Ocular allergy (OA) is characterised by ocular surface itchiness, redness, and inflammation in response to allergen exposure. The primary aim of this study was to assess differences in the human tear metabolome and lipidome between OA and healthy controls (HCs) across peak allergy (spring-summer) and off-peak (autumn-winter) seasons in Victoria, Australia. A total of 19 participants (14 OA, 5 HCs) aged 18-45 were recruited and grouped by allergy questionnaire score. Metabolites and lipids from tear samples were analysed using mass spectrometry. Data were analysed using TraceFinder and Metaboanalyst. Metabolomics analysis showed 12 differentially expressed (DE) metabolites between those with OA and the HCs during the peak allergy season, and 24 DE metabolites were found in the off-peak season. The expression of niacinamide was upregulated in OA sufferers vs. HCs across both seasons (p ≤ 0.05). A total of 6 DE lipids were DE between those with OA and the HCs during the peak season, and 24 were DE in the off-peak season. Dysregulated metabolites affected oxidative stress, inflammation, and homeostasis across seasons, suggesting a link between OA-associated itch and ocular surface damage via eye rubbing. Tear lipidome changes were minimal between but suggested tear film destabilisation and thinning. Such metabolipodome findings may pave new and exciting ways for effective diagnostics and therapeutics for OA sufferers in the future.

WikiPathways 2024: next generation pathway database

WikiPathways (wikipathways.org) is an open-source biological pathway database. Collaboration and open science are pivotal to the success of WikiPathways. Here we highlight the continuing efforts supporting WikiPathways, content growth and collaboration among pathway researchers. As an evolving database, there is a growing need for WikiPathways to address and overcome technical challenges. In this direction, WikiPathways has undergone major restructuring, enabling a renewed approach for sharing and curating pathway knowledge, thus providing stability for the future of community pathway curation. The website has been redesigned to improve and enhance user experience. This next generation of WikiPathways continues to support existing features while improving maintainability of the database and facilitating community input by providing new functionality and leveraging automation.

Nucleotide, Phospholipid, and Kynurenine Metabolites Are Robustly Associated with COVID-19 Severity and Time of Plasma Sample Collection in a Prospective Cohort Study

Understanding the molecular underpinnings of disease severity and progression in human studies is necessary to develop metabolism-related preventative strategies for severe COVID-19. Metabolites and metabolic pathways that predispose individuals to severe disease are not well understood. In this study, we generated comprehensive plasma metabolomic profiles in >550 patients from the Longitudinal EMR and Omics COVID-19 Cohort. Samples were collected before (n = 441), during (n = 86), and after (n = 82) COVID-19 diagnosis, representing 555 distinct patients, most of which had single timepoints. Regression models adjusted for demographics, risk factors, and comorbidities, were used to determine metabolites associated with predisposition to and/or persistent effects of COVID-19 severity, and metabolite changes that were transient/lingering over the disease course. Sphingolipids/phospholipids were negatively associated with severity and exhibited lingering elevations after disease, while modified nucleotides were positively associated with severity and had lingering decreases after disease. Cytidine and uridine metabolites, which were positively and negatively associated with COVID-19 severity, respectively, were acutely elevated, reflecting the particular importance of pyrimidine metabolism in active COVID-19. This is the first large metabolomics study using COVID-19 plasma samples before, during, and/or after disease. Our results lay the groundwork for identifying putative biomarkers and preventive strategies for severe COVID-19.

Metabolomic epidemiology offers insights into disease aetiology

Metabolomic epidemiology is the high-throughput study of the relationship between metabolites and health-related traits. This emerging and rapidly growing field has improved our understanding of disease aetiology and contributed to advances in precision medicine. As the field continues to develop, metabolomic epidemiology could lead to the discovery of diagnostic biomarkers predictive of disease risk, aiding in earlier disease detection and better prognosis. In this Review, we discuss key advances facilitated by the field of metabolomic epidemiology for a range of conditions, including cardiometabolic diseases, cancer, Alzheimer's disease and COVID-19, with a focus on potential clinical utility. Core principles in metabolomic epidemiology, including study design, causal inference methods and multi-omic integration, are briefly discussed. Future directions required for clinical translation of metabolomic epidemiology findings are summarized, emphasizing public health implications. Further work is needed to establish which metabolites reproducibly improve clinical risk prediction in diverse populations and are causally related to disease progression.

Recent advances in mass spectrometry-based computational metabolomics

The computational metabolomics field brings together computer scientists, bioinformaticians, chemists, clinicians, and biologists to maximize the impact of metabolomics across a wide array of scientific and medical disciplines. The field continues to expand as modern instrumentation produces datasets with increasing complexity, resolution, and sensitivity. These datasets must be processed, annotated, modeled, and interpreted to enable biological insight. Techniques for visualization, integration (within or between omics), and interpretation of metabolomics data have evolved along with innovation in the databases and knowledge resources required to aid understanding. In this review, we highlight recent advances in the field and reflect on opportunities and innovations in response to the most pressing challenges. This review was compiled from discussions from the 2022 Dagstuhl seminar entitled "Computational Metabolomics: From Spectra to Knowledge".

Integrative transcriptomic and metabolic analyses of the mammalian hibernating brain identifies a key role for succinate dehydrogenase in ischemic tolerance

Ischemic stroke results in a loss of tissue homeostasis and integrity, the underlying pathobiology of which stems primarily from the depletion of cellular energy stores and perturbation of available metabolites ¹ . Hibernation in thirteen-lined ground squirrels (TLGS), Ictidomys tridecemlineatus , provides a natural model of ischemic tolerance as these mammals undergo prolonged periods of critically low cerebral blood flow without evidence of central nervous system (CNS) damage ² . Studying the complex interplay of genes and metabolites that unfolds during hibernation may provide novel insights into key regulators of cellular homeostasis during brain ischemia. Herein, we interrogated the molecular profiles of TLGS brains at different time points within the hibernation cycle via RNA sequencing coupled with untargeted metabolomics. We demonstrate that hibernation in TLGS leads to major changes in the expression of genes involved in oxidative phosphorylation and this is correlated with an accumulation of the tricarboxylic acid (TCA) cycle intermediates citrate, cis-aconitate, and α-ketoglutarate-αKG. Integration of the gene expression and metabolomics datasets led to the identification of succinate dehydrogenase (SDH) as the critical enzyme during hibernation, uncovering a break in the TCA cycle at that level. Accordingly, the SDH inhibitor dimethyl malonate (DMM) was able to rescue the effects of hypoxia on human neuronal cells in vitro and in mice subjected to permanent ischemic stroke in vivo . Our findings indicate that studying the regulation of the controlled metabolic depression that occurs in hibernating mammals may lead to novel therapeutic approaches capable of increasing ischemic tolerance in the CNS.

Exploiting Multi-Omics Profiling and Systems Biology to Investigate Functions of TOMM34

Although modern biology is now in the post-genomic era with vastly increased access to high-quality data, the set of human genes with a known function remains far from complete. This is especially true for hundreds of mitochondria-associated genes, which are under-characterized and lack clear functional annotation. However, with the advent of multi-omics profiling methods coupled with systems biology algorithms, the cellular role of many such genes can be elucidated. Here, we report genes and pathways associated with TOMM34, Translocase of Outer Mitochondrial Membrane, which plays role in the mitochondrial protein import as a part of cytosolic complex together with Hsp70/Hsp90 and is upregulated in various cancers. We identified genes, proteins, and metabolites altered in TOMM34-/- HepG2 cells. To our knowledge, this is the first attempt to study the functional capacity of TOMM34 using a multi-omics strategy. We demonstrate that TOMM34 affects various processes including oxidative phosphorylation, citric acid cycle, metabolism of purine, and several amino acids. Besides the analysis of already known pathways, we utilized de novo network enrichment algorithm to extract novel perturbed subnetworks, thus obtaining evidence that TOMM34 potentially plays role in several other cellular processes, including NOTCH-, MAPK-, and STAT3-signaling. Collectively, our findings provide new insights into TOMM34's cellular functions.

IntLIM 2.0: identifying multi-omic relationships dependent on discrete or continuous phenotypic measurements

Motivation

IntLIM uncovers phenotype-dependent linear associations between two types of analytes (e.g. genes and metabolites) in a multi-omic dataset, which may reflect chemically or biologically relevant relationships.

Results

The new IntLIM R package includes newly added support for generalized data types, covariate correction, continuous phenotypic measurements, model validation and unit testing. IntLIM analysis uncovered biologically relevant gene-metabolite associations in two separate datasets, and the run time is improved over baseline R functions by multiple orders of magnitude.

Availability and implementation

IntLIM is available as an R package with a detailed vignette (https://github.com/ncats/IntLIM) and as an R Shiny app (see Supplementary Figs S1-S6) (https://intlim.ncats.io/).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.