Pathway analysis in metabolomics: Recommendations for the use of over-representation analysis

Machine learning of cellular metabolic rewiring

Metabolic rewiring allows cells to adapt their metabolism in response to evolving environmental conditions. Traditional metabolomics techniques, whether targeted or untargeted, often struggle to interpret these adaptive shifts. Here, we introduce MetaboLiteLearner, a machine learning framework that harnesses the detailed fragmentation patterns from electron ionization (EI) collected in scan mode during gas chromatography/mass spectrometry (GC/MS) to predict abundance changes in metabolically adapted cells. When tested on breast cancer cells with different preferences to metastasize to specific organs, MetaboLiteLearner predicted the impact of metabolic rewiring on metabolites withheld from the training dataset using only the EI spectra, without metabolite identification or pre-existing knowledge of metabolic networks. The model learned captures shared and unique metabolomic shifts between brain- and lung-homing metastatic lineages, suggesting potential organ-tailored cellular adaptations. Integrating machine learning and metabolomics paves the way for new insights into complex cellular adaptations.

Metabolic Differences in Diabetic Kidney Disease Patients with Normoalbuminuria versus Moderately Increased Albuminuria

Key Points

The pathophysiological mechanisms of diabetic kidney disease (DKD) with normal (nonalbuminuric DKD) versus moderately increased albuminuria (A-DKD) are not well-understood. Fatty acid biosynthesis and oxydation, gluconeogenesis, TCA cycle, and glucose-alanine cycle were more disturbed in patients with A-DKD compared with those with nonalbuminuric DKD with identical eGFR. DKD patients with and without microalbuminuria could represent different clinical phenotypes.

Background

The pathophysiological mechanisms of diabetic kidney disease (DKD) with normal versus moderately increased albuminuria (nonalbuminuric DKD [NA-DKD] and A-DKD) are currently not well-understood and could have implications for diagnosis and treatment.

Methods

Fourteen patients with NA-DKD with urine albumin–creatinine ratio <3 mg/mmol, 26 patients with A-DKD with albumin–creatinine ratio 3–29 mg/mmol, and 60 age- and sex-matched healthy controls were randomly chosen from a population-based cohort study (Nord-Trøndelag Health Study-3, Norway). Seventy-four organic acids, 21 amino acids, 21 biogenic acids, 40 acylcarnitines, 14 sphingomyelins, and 88 phosphatidylcholines were quantified in urine. One hundred forty-six patients with diabetes from the US-based Chronic Renal Insufficiency Cohort study were used to verify main findings.

Results

Patients with NA-DKD and A-DKD had similar age, kidney function, diabetes treatment, and other traditional risk factors. Still, partial least-squares discriminant analysis showed strong metabolite-based separation (R2, 0.82; Q2, 0.52), with patients with NA-DKD having a metabolic profile positioned between the profiles of healthy controls and patients with A-DKD. Seventy-five metabolites contributed significantly to separation between NA-DKD and A-DKD (variable importance in projection scores ≥1.0) with propionylcarnitine (C3), phosphatidylcholine C38:4, medium-chained (C8) fatty acid octenedioic acid, and lactic acid as the top metabolites (variable importance in projection scores, 2.7–2.2). Compared with patients with NA-DKD, those with A-DKD had higher levels of short-chained acylcarnitines, higher long-chained fatty acid levels with more double bounds, higher branched-chain amino acid levels, and lower TCA cycle intermediates. The main findings were similar by random forest analysis and in the Chronic Renal Insufficiency Cohort study. Formal enrichment analysis indicated that fatty acid biosynthesis and oxydation, gluconeogenesis, TCA cycle, and glucose-alanine cycle were more disturbed in patients with A-DKD compared with those with NA-DKD with identical eGFR. We also found indications of a Warburg-like effect in patients with A-DKD (i.e. , metabolism of glucose to lactate despite adequate oxygen).

Conclusion

DKD patients with normoalbuminuria differ substantially in their metabolic disturbances compared with patients with moderately increase albuminuria and could represent different clinical phenotypes.

Differential Correlations Informed Metabolite Set Enrichment Analysis to Decipher Metabolic Heterogeneity of Disease

Metabolic pathways are regarded as functional and basic components of the biological system. In metabolomics, metabolite set enrichment analysis (MSEA) is often used to identify the altered metabolic pathways (metabolite sets) associated with phenotypes of interest (POI), e.g., disease. However, in most studies, MSEA suffers from the limitation of low metabolite coverage. Random walk (RW)-based algorithms can be used to propagate the perturbation of detected metabolites to the undetected metabolites through a metabolite network model prior to MSEA. Nevertheless, most of the existing RW-based algorithms run on a general metabolite network constructed based on public databases, such as KEGG, without taking into consideration the potential influence of POI on the metabolite network, which may reduce the phenotypic specificities of the MSEA results. To solve this problem, a novel pathway analysis strategy, namely, differential correlation-informed MSEA (dci-MSEA), is proposed in this paper. Statistically, differential correlations between metabolites are used to evaluate the influence of POI on the metabolite network, so that a phenotype-specific metabolite network is constructed for RW-based propagation. The experimental results show that dci-MSEA outperforms the conventional RW-based MSEA in identifying the altered metabolic pathways associated with colorectal cancer. In addition, by incorporating the individual-specific metabolite network, the dci-MSEA strategy is easily extended to disease heterogeneity analysis. Here, dci-MSEA was used to decipher the heterogeneity of colorectal cancer. The present results highlight the clustering of colorectal cancer samples with their cluster-specific selection of differential pathways and demonstrate the feasibility of dci-MSEA in heterogeneity analysis. Taken together, the proposed dci-MSEA may provide insights into disease mechanisms and determination of disease heterogeneity.

Epilepsies of presumed genetic etiology show enrichment of rare variants that occur in the general population

Previous studies suggested that severe epilepsies, e.g., developmental and epileptic encephalopathies (DEEs), are mainly caused by ultra-rare de novo genetic variants. For milder disease, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 control individuals. Here, we separately analyzed three different groups of epilepsies: severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in control individuals with an allele count ≥ 1 and a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD ≥ 20), and to have an odds ratio in individuals with epilepsy ≥ 2. We identified genes enriched with QRVs primarily in NAFE (n = 72), followed by GGE (n = 32) and DEE (n = 21). This suggests that rare variants may play a more important role for causality of NAFE than for DEE. Moreover, we found that genes harboring QRVs, e.g., HSGP2, FLNA, or TNC, encode proteins that are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE that occur also in the general population, while in DEE and GGE, the contribution of such variants appears more limited.

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".

Statistical analysis of high-dimensional biomedical data: a gentle introduction to analytical goals, common approaches and challenges

BACKGROUND

In high-dimensional data (HDD) settings, the number of variables associated with each observation is very large. Prominent examples of HDD in biomedical research include omics data with a large number of variables such as many measurements across the genome, proteome, or metabolome, as well as electronic health records data that have large numbers of variables recorded for each patient. The statistical analysis of such data requires knowledge and experience, sometimes of complex methods adapted to the respective research questions.

METHODS

Advances in statistical methodology and machine learning methods offer new opportunities for innovative analyses of HDD, but at the same time require a deeper understanding of some fundamental statistical concepts. Topic group TG9 "High-dimensional data" of the STRATOS (STRengthening Analytical Thinking for Observational Studies) initiative provides guidance for the analysis of observational studies, addressing particular statistical challenges and opportunities for the analysis of studies involving HDD. In this overview, we discuss key aspects of HDD analysis to provide a gentle introduction for non-statisticians and for classically trained statisticians with little experience specific to HDD.

RESULTS

The paper is organized with respect to subtopics that are most relevant for the analysis of HDD, in particular initial data analysis, exploratory data analysis, multiple testing, and prediction. For each subtopic, main analytical goals in HDD settings are outlined. For each of these goals, basic explanations for some commonly used analysis methods are provided. Situations are identified where traditional statistical methods cannot, or should not, be used in the HDD setting, or where adequate analytic tools are still lacking. Many key references are provided.

CONCLUSIONS

This review aims to provide a solid statistical foundation for researchers, including statisticians and non-statisticians, who are new to research with HDD or simply want to better evaluate and understand the results of HDD analyses.

Chewing the fat: How lipidomics is changing our understanding of human health and disease in 2022

Lipids are biological molecules that play vital roles in all living organisms. They perform many cellular functions, such as 1) forming cellular and subcellular membranes, 2) storing and using energy, and 3) serving as chemical messengers during intra- and inter-cellular signal transduction. The large-scale study of the pathways and networks of cellular lipids in biological systems is called "lipidomics" and is one of the fastest-growing omics technologies of the last two decades. With state-of-the-art mass spectrometry instrumentation and sophisticated data handling, clinical studies show how human lipid composition changes in health and disease, thereby making it a valuable medium to collect for clinical applications, such as disease diagnostics, therapeutic decision-making, and drug development. This review gives a comprehensive overview of current workflows used in clinical research, from sample collection and preparation to data and clinical interpretations. This is followed by an appraisal of applications in 2022 and a perspective on the exciting future of clinical lipidomics.

Interpreting omics data with pathway enrichment analysis

Pathway enrichment analysis is indispensable for interpreting omics datasets and generating hypotheses. However, the foundations of enrichment analysis remain elusive to many biologists. Here, we discuss best practices in interpreting different types of omics data using pathway enrichment analysis and highlight the importance of considering intrinsic features of various types of omics data. We further explain major components that influence the outcomes of a pathway enrichment analysis, including defining background sets and choosing reference annotation databases. To improve reproducibility, we describe how to standardize reporting methodological details in publications. This article aims to serve as a primer for biologists to leverage the wealth of omics resources and motivate bioinformatics tool developers to enhance the power of pathway enrichment analysis.

The role of rare genetic variants enrichment in epilepsies of presumed genetic etiology

Previous studies suggested that severe epilepsies e.g., developmental and epileptic encephalopathies (DEE) are mainly caused by ultra-rare de novo genetic variants. For milder phenotypes, rare genetic variants could contribute to the phenotype. To determine the importance of rare variants for different epilepsy types, we analyzed a whole-exome sequencing cohort of 9,170 epilepsy-affected individuals and 8,436 controls. Here, we separately analyzed three different groups of epilepsies : severe DEEs, genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We required qualifying rare variants (QRVs) to occur in controls at a minor allele frequency ≤ 1:1,000, to be predicted as deleterious (CADD≥20), and to have an odds ratio in epilepsy cases ≥2. We identified genes enriched with QRVs in DEE (n=21), NAFE (n=72), and GGE (n=32) - the number of enriched genes are found greatest in NAFE and least in DEE. This suggests that rare variants may play a more important role for causality of NAFE than in DEE. Moreover, we found that QRV-carrying genes e.g., HSGP2, FLNA or TNC are involved in structuring the brain extracellular matrix. The present study confirms an involvement of rare variants for NAFE, while in DEE and GGE, the contribution of such variants appears more limited.

Differences in uterine and serum metabolome associated with metritis in dairy cows

Objectives were to evaluate differences in the uterine and serum metabolomes associated with metritis in dairy cows. Vaginal discharge was evaluated using a Metricheck device (Simcro) at 5, 7, and 11 d in milk (DIM; herd 1) or 4, 6, 8, 10, and 12 DIM (herd 2). Cows with reddish or brownish, watery, and fetid discharge were diagnosed with metritis (n = 24). Cows with metritis were paired with herdmates without metritis (i.e., clear mucous vaginal discharge or clear lochia with ≤50% of pus) based on DIM and parity (n = 24). Day of metritis diagnosis was considered study d 0. All cows diagnosed with metritis received antimicrobial therapy. The metabolome of uterine lavage collected on d 0 and 5, and serum samples collected on d 0 were evaluated using untargeted gas chromatography time-of-flight mass spectrometry. Normalized data were subjected to multivariate canonical analysis of population using the MultBiplotR and MixOmics packages in R Studio. Univariate analyses including t-test, principal component analyses, partial least squares discriminant analyses, and pathway analyses were conducted using Metaboanalyst. The uterine metabolome differed between cows with and without metritis on d 0. Differences in the uterine metabolome associated with metritis on d 0 were related to the metabolism of butanoate, amino acids (i.e., glycine, serine, threonine, alanine, aspartate, and glutamate), glycolysis and gluconeogenesis, and the tricarboxylic acid cycle. No differences in the serum metabolome were observed between cows diagnosed with metritis and counterparts without metritis on d 0. Similarly, no differences in uterine metabolome were observed between cows with metritis and counterparts not diagnosed with metritis on d 5. These results indicate that the establishment of metritis in dairy cows is associated with local disturbances in amino acid, lipid, and carbohydrate metabolism in the uterus. The lack of differences in the uterine metabolome on d 5 indicates that processes implicated with the disease are reestablished by d 5 after diagnosis and treatment.

Metabolomics of asthma, COPD, and asthma-COPD overlap: an overview

The two common progressive lung diseases, asthma and chronic obstructive pulmonary disease (COPD), are the leading causes of morbidity and mortality worldwide. Asthma-COPD overlap, referred to as ACO, is another complex pulmonary disease that manifests itself with features of both asthma and COPD. The disease has no clear diagnostic or therapeutic guidelines, thereby making both diagnosis and treatment challenging. Though a number of studies on ACO have been documented, gaps in knowledge regarding the pathophysiologic mechanism of this disorder exist. Addressing this issue is an urgent need for improved diagnostic and therapeutic management of the disease. Metabolomics, an increasingly popular technique, reveals the pathogenesis of complex diseases and holds promise in biomarker discovery. This comprehensive narrative review, comprising 99 original research articles in the last five years (2017-2022), summarizes the scientific advances in terms of metabolic alterations in patients with asthma, COPD, and ACO. The analytical tools, nuclear magnetic resonance (NMR), gas chromatography-mass spectrometry (GC-MS), and liquid chromatography-mass spectrometry (LC-MS), commonly used to study the expression of the metabolome, are discussed. Challenges frequently encountered during metabolite identification and quality assessment are highlighted. Bridging the gap between phenotype and metabotype is envisioned in the future.

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.

Ten quick tips for computational analysis of medical images

Medical imaging is a great asset for modern medicine, since it allows physicians to spatially interrogate a disease site, resulting in precise intervention for diagnosis and treatment, and to observe particular aspect of patients' conditions that otherwise would not be noticeable. Computational analysis of medical images, moreover, can allow the discovery of disease patterns and correlations among cohorts of patients with the same disease, thus suggesting common causes or providing useful information for better therapies and cures. Machine learning and deep learning applied to medical images, in particular, have produced new, unprecedented results that can pave the way to advanced frontiers of medical discoveries. While computational analysis of medical images has become easier, however, the possibility to make mistakes or generate inflated or misleading results has become easier, too, hindering reproducibility and deployment. In this article, we provide ten quick tips to perform computational analysis of medical images avoiding common mistakes and pitfalls that we noticed in multiple studies in the past. We believe our ten guidelines, if taken into practice, can help the computational-medical imaging community to perform better scientific research that eventually can have a positive impact on the lives of patients worldwide.

Longitudinal associations between metabolites and long-term exposure to ambient air pollution: Results from the KORA cohort study

BACKGROUND

Long-term exposure to air pollution has been associated with cardiopulmonary diseases, while the underlying mechanisms remain unclear.

OBJECTIVES

To investigate changes in serum metabolites associated with long-term exposure to air pollution and explore the susceptibility characteristics.

METHODS

We used data from the German population-based Cooperative Health Research in the Region of Augsburg (KORA) S4 survey (1999-2001) and two follow-up examinations (F4: 2006-08 and FF4: 2013-14). Mass-spectrometry-based targeted metabolomics was used to quantify metabolites among serum samples. Only participants with repeated metabolites measurements were included in the current analysis. Land-use regression (LUR) models were used to estimate annual average concentrations of ultrafine particles, particulate matter (PM) with an aerodynamic diameter less than 10 μm (PM₁₀), coarse particles (PM_coarse), fine particles, PM_2.5 absorbance (a proxy of elemental carbon related to traffic exhaust, PM_2.5abs), nitrogen oxides (NO₂, NO_x), and ozone at individuals' residences. We applied confounder-adjusted mixed-effects regression models to examine the associations between long-term exposure to air pollution and metabolites.

RESULTS

Among 9,620 observations from 4,261 KORA participants, we included 5,772 (60.0%) observations from 2,583 (60.6%) participants in this analysis. Out of 108 metabolites that passed stringent quality control across three study points in time, we identified nine significant negative associations between phosphatidylcholines (PCs) and ambient pollutants at a Benjamini-Hochberg false discovery rate (FDR) corrected p-value < 0.05. The strongest association was seen for an increase of 0.27 μg/m³ (interquartile range) in PM_2.5abs and decreased phosphatidylcholine acyl-alkyl C36:3 (PC ae C36:3) concentrations [percent change in the geometric mean: -2.5% (95% confidence interval: -3.6%, -1.5%)].

CONCLUSIONS

Our study suggested that long-term exposure to air pollution is associated with metabolic alterations, particularly in PCs with unsaturated long-chain fatty acids. These findings might provide new insights into potential mechanisms for air pollution-related adverse outcomes.

Single sample pathway analysis in metabolomics: performance evaluation and application

BACKGROUND

Single sample pathway analysis (ssPA) transforms molecular level omics data to the pathway level, enabling the discovery of patient-specific pathway signatures. Compared to conventional pathway analysis, ssPA overcomes the limitations by enabling multi-group comparisons, alongside facilitating numerous downstream analyses such as pathway-based machine learning. While in transcriptomics ssPA is a widely used technique, there is little literature evaluating its suitability for metabolomics. Here we provide a benchmark of established ssPA methods (ssGSEA, GSVA, SVD (PLAGE), and z-score) alongside the evaluation of two novel methods we propose: ssClustPA and kPCA, using semi-synthetic metabolomics data. We then demonstrate how ssPA can facilitate pathway-based interpretation of metabolomics data by performing a case-study on inflammatory bowel disease mass spectrometry data, using clustering to determine subtype-specific pathway signatures.

RESULTS

While GSEA-based and z-score methods outperformed the others in terms of recall, clustering/dimensionality reduction-based methods provided higher precision at moderate-to-high effect sizes. A case study applying ssPA to inflammatory bowel disease data demonstrates how these methods yield a much richer depth of interpretation than conventional approaches, for example by clustering pathway scores to visualise a pathway-based patient subtype-specific correlation network. We also developed the sspa python package (freely available at https://pypi.org/project/sspa/ ), providing implementations of all the methods benchmarked in this study.

CONCLUSION

This work underscores the value ssPA methods can add to metabolomic studies and provides a useful reference for those wishing to apply ssPA methods to metabolomics data.

Exploring Dynamic Metabolome of the HepG2 Cell Line: Rise and Fall

Both biological and technical variations can discredit the reliability of obtained data in omics studies. In this technical note, we investigated the effect of prolonged cultivation of the HepG2 hepatoma cell line on its metabolomic profile. Using the GC × GC-MS approach, we determined the degree of metabolic variability across HepG2 cells cultured in uniform conditions for 0, 5, 10, 15, and 20 days. Post-processing of obtained data revealed substantial changes in relative abundances of 110 metabolites among HepG2 samples under investigation. Our findings have implications for interpreting metabolomic results obtained from immortal cells, especially in longitudinal studies. There are still plenty of unanswered questions regarding metabolomics variability and many potential areas for future targeted and panoramic research. However, we suggest that the metabolome of cell lines is unstable and may undergo significant transformation over time, even if the culture conditions remain the same. Considering metabolomics variability on a relatively long-term basis, careful experimentation with particular attention to control samples is required to ensure reproducibility and relevance of the research results when testing both fundamentally and practically significant hypotheses.

Human Plasma Proteome During Normal Pregnancy

The human plasma proteome is underexplored despite its potential value for monitoring health and disease. Herein, using a recently developed aptamer-based platform, we profiled 7288 proteins in 528 plasma samples from 91 normal pregnancies (Gene Expression Omnibus identifier GSE206454). The coefficient of variation was <20% for 93% of analytes (median 7%), and a cross-platform correlation for selected key angiogenic and anti-angiogenic proteins was significant. Gestational age was associated with changes in 953 proteins, including highly modulated placenta- and decidua-specific proteins, and they were enriched in biological processes including regulation of growth, angiogenesis, immunity, and inflammation. The abundance of proteins corresponding to RNAs specific to populations of cells previously described by single-cell RNA-Seq analysis of the placenta was highly modulated throughout gestation. Furthermore, machine learning-based prediction of gestational age and of time from sampling to term delivery compared favorably with transcriptomic models (mean absolute error of 2 weeks). These results suggested that the plasma proteome may provide a non-invasive readout of placental cellular dynamics and serve as a blueprint for investigating obstetrical disease.

Survey for Computer-Aided Tools and Databases in Metabolomics

Metabolomics has advanced from innovation and functional genomics tools and is currently a basis in the big data-led precision medicine era. Metabolomics is promising in the pharmaceutical field and clinical research. However, due to the complexity and high throughput data generated from such experiments, data mining and analysis are significant challenges for researchers in the field. Therefore, several efforts were made to develop a complete workflow that helps researchers analyze data. This paper introduces a review of the state-of-the-art computer-aided tools and databases in metabolomics established in recent years. The paper provides computational tools and resources based on functionality and accessibility and provides hyperlinks to web pages to download or use. This review aims to present the latest computer-aided tools, databases, and resources to the metabolomics community in one place.

Evaluation of the metabolomic profile through 1H-NMR spectroscopy in ewes affected by postpartum hyperketonemia

Ketosis is one of the most important health problems in dairy sheep. The aim of this study was to evaluate the metabolic alterations in hyperketonemic (HYK) ewes. Forty-six adult Sardinian ewes were enrolled between 7 ± 3 days post-partum. Blood samples were collected from the jugular vein using Venosafe tubes containing clot activator from jugular vein after clinical examination. The concentration of β-hydroxybutyrate (BHB) was determined in serum and used to divide ewes into assign ewes into: Non-HYK (serum BHB < 0.80 mmol/L) and HYK (serum BHB ≥ 0.80 mmol/L) groups. Animal data and biochemical parameters of groups were examined with one-way ANOVA, and metabolite differences were tested using a t-test. A robust principal component analysis model and a heatmap were used to highlight common trends among metabolites. Over-representation analysis was performed to investigate metabolic pathways potentially altered in connection with BHB alterations. The metabolomic analysis identified 54 metabolites with 14 different between groups. These metabolites indicate altered ruminal microbial populations and fermentations; an interruption of the tricarboxylic acid cycle; initial lack of glucogenic substrates; mobilization of body reserves; the potential alteration of electron transport chain; influence on urea synthesis; alteration of nervous system, inflammatory response, and immune cell function.

A targeted metabolic analysis of football players and its association to player load: Comparison between women and men profiles

Professional athletes undertake a variety of training programs to enhance their physical performance, technical-tactical skills, while protecting their health and well-being. Regular exercise induces widespread changes in the whole body in an extremely complex network of signaling, and evidence indicates that phenotypical sex differences influence the physiological adaptations to player load of professional athletes. Despite that there remains an underrepresentation of women in clinical studies in sports, including football. The objectives of this study were twofold: to study the association between the external load (EPTS) and urinary metabolites as a surrogate of the adaptation to training, and to assess the effect of sex on the physiological adaptations to player load in professional football players. Targeted metabolic analysis of aminoacids, and tryptophan and phenylalanine metabolites detected progressive changes in the urinary metabolome associated with the external training load in men and women’s football teams. Overrepresentation analysis and multivariate analysis of metabolic data showed significant differences of the effect of training on the metabolic profiles in the men and women teams analyzed. Collectively, our results demonstrate that the development of metabolic models of adaptation in professional football players can benefit from the separate analysis of women and men teams, providing more accurate insights into how adaptation to the external load is related to changes in the metabolic phenotypes. Furthermore, results support the use of metabolomics to understand changes in specific metabolic pathways provoked by the training process.

Untargeted approach to investigating the metabolomics profile of type 2 diabetes emiratis

Type 2 Diabetes (T2D) is expected to be the seventh most significant cause of death worldwide by 2030. Although research into its mechanism has received the attention it deserves, our understanding of T2D is still limited. This case-control study employs untargeted metabolomics to explore novel T2D plasma biomarkers in the Emirati population. Ninety-two UAE nationals were included in the cohort, with fifty T2D and forty-two non-T2D profiles. Participants were then stratified into three groups based on metabolic profiles, clinically verified diabetic status, and current HbA1c values: namely controlled diabetics, uncontrolled diabetics and prediabetics, and non-diabetics. The study identified fifteen significant differentially abundant metabolites between the uncontrolled diabetics group and the prediabetics or controlled diabetics group. Interestingly, some metabolites essential for the corticosteroid and thyroid signaling pathways were found to be significantly elevated in poorly controlled T2D, including cortisol, glycocholic acid, bile acids, thyroxine, and the tryptophan metabolite, 5-hydroxyindoleacetic acid. These findings align with those from prior western cohorts and suggest an intriguing linkage between T2D glycemic control and thyroid and adrenal signaling that may provide new diagnostic and prognostic indicators. RESEARCH SIGNIFICANCE: This study investigates the underlooked metabolomic role and correlation with T2D in the UAE population. The report indicates fifteen significant differentially abundant metabolites between on diabetics, uncontrolled diabetics and or controlled diabetics or prediabetics. This panel of metabolites such as thyroxine and corticosteroids should be considered further as potential diagnostic or prognostic biomarkers for T2D in the region.

Nine quick tips for pathway enrichment analysis

Pathway enrichment analysis (PEA) is a computational biology method that identifies biological functions that are overrepresented in a group of genes more than would be expected by chance and ranks these functions by relevance. The relative abundance of genes pertinent to specific pathways is measured through statistical methods, and associated functional pathways are retrieved from online bioinformatics databases. In the last decade, along with the spread of the internet, higher availability of computational resources made PEA software tools easy to access and to use for bioinformatics practitioners worldwide. Although it became easier to use these tools, it also became easier to make mistakes that could generate inflated or misleading results, especially for beginners and inexperienced computational biologists. With this article, we propose nine quick tips to avoid common mistakes and to out a complete, sound, thorough PEA, which can produce relevant and robust results. We describe our nine guidelines in a simple way, so that they can be understood and used by anyone, including students and beginners. Some tips explain what to do before starting a PEA, others are suggestions of how to correctly generate meaningful results, and some final guidelines indicate some useful steps to properly interpret PEA results. Our nine tips can help users perform better pathway enrichment analyses and eventually contribute to a better understanding of current biology.

Metabolomics in pulmonary medicine - extracting the most from your data

The metabolome enables unprecedented insight into biochemistry, providing an integrated signature of the genome, transcriptome, proteome and exposome. Measurement requires rigorous protocols combined with specialised data analysis to achieve its promise.https://bit.ly/3yPiYkQ

On the influence of several factors on pathway enrichment analysis

Pathway enrichment analysis has become a widely used knowledge-based approach for the interpretation of biomedical data. Its popularity has led to an explosion of both enrichment methods and pathway databases. While the elegance of pathway enrichment lies in its simplicity, multiple factors can impact the results of such an analysis, which may not be accounted for. Researchers may fail to give influential aspects their due, resorting instead to popular methods and gene set collections, or default settings. Despite ongoing efforts to establish set guidelines, meaningful results are still hampered by a lack of consensus or gold standards around how enrichment analysis should be conducted. Nonetheless, such concerns have prompted a series of benchmark studies specifically focused on evaluating the influence of various factors on pathway enrichment results. In this review, we organize and summarize the findings of these benchmarks to provide a comprehensive overview on the influence of these factors. Our work covers a broad spectrum of factors, spanning from methodological assumptions to those related to prior biological knowledge, such as pathway definitions and database choice. In doing so, we aim to shed light on how these aspects can lead to insignificant, uninteresting or even contradictory results. Finally, we conclude the review by proposing future benchmarks as well as solutions to overcome some of the challenges, which originate from the outlined factors.

KATP channel dependent heart multiome atlas

Plasmalemmal ATP sensitive potassium (K_ATP) channels are recognized metabolic sensors, yet their cellular reach is less well understood. Here, transgenic Kir6.2 null hearts devoid of the K_ATP channel pore underwent multiomics surveillance and systems interrogation versus wildtype counterparts. Despite maintained organ performance, the knockout proteome deviated beyond a discrete loss of constitutive K_ATP channel subunits. Multidimensional nano-flow liquid chromatography tandem mass spectrometry resolved 111 differentially expressed proteins and their expanded network neighborhood, dominated by metabolic process engagement. Independent multimodal chemometric gas and liquid chromatography mass spectrometry unveiled differential expression of over one quarter of measured metabolites discriminating the Kir6.2 deficient heart metabolome. Supervised class analogy ranking and unsupervised enrichment analysis prioritized nicotinamide adenine dinucleotide (NAD⁺), affirmed by extensive overrepresentation of NAD⁺ associated circuitry. The remodeled metabolome and proteome revealed functional convergence and an integrated signature of disease susceptibility. Deciphered cardiac patterns were traceable in the corresponding plasma metabolome, with tissue concordant plasma changes offering surrogate metabolite markers of myocardial latent vulnerability. Thus, Kir6.2 deficit precipitates multiome reorganization, mapping a comprehensive atlas of the K_ATP channel dependent landscape.

Studying Metabolism by NMR-Based Metabolomics

During the past few decades, the direct analysis of metabolic intermediates in biological samples has greatly improved the understanding of metabolic processes. The most used technologies for these advances have been mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy. NMR is traditionally used to elucidate molecular structures and has now been extended to the analysis of complex mixtures, as biological samples: NMR-based metabolomics. There are however other areas of small molecule biochemistry for which NMR is equally powerful. These include the quantification of metabolites (qNMR); the use of stable isotope tracers to determine the metabolic fate of drugs or nutrients, unravelling of new metabolic pathways, and flux through pathways; and metabolite-protein interactions for understanding metabolic regulation and pharmacological effects. Computational tools and resources for automating analysis of spectra and extracting meaningful biochemical information has developed in tandem and contributes to a more detailed understanding of systems biochemistry. In this review, we highlight the contribution of NMR in small molecule biochemistry, specifically in metabolic studies by reviewing the state-of-the-art methodologies of NMR spectroscopy and future directions.

A Review on Differential Abundance Analysis Methods for Mass Spectrometry-Based Metabolomic Data

This review presents an overview of the statistical methods on differential abundance (DA) analysis for mass spectrometry (MS)-based metabolomic data. MS has been widely used for metabolomic abundance profiling in biological samples. The high-throughput data produced by MS often contain a large fraction of zero values caused by the absence of certain metabolites and the technical detection limits of MS. Various statistical methods have been developed to characterize the zero-inflated metabolomic data and perform DA analysis, ranging from simple tests to more complex models including parametric, semi-parametric, and non-parametric approaches. In this article, we discuss and compare DA analysis methods regarding their assumptions and statistical modeling techniques.

NMR in Metabolomics: From Conventional Statistics to Machine Learning and Neural Network Approaches

NMR measurements combined with chemometrics allow achieving a great amount of information for the identification of potential biomarkers responsible for a precise metabolic pathway. These kinds of data are useful in different fields, ranging from food to biomedical fields, including health science. The investigation of the whole set of metabolites in a sample, representing its fingerprint in the considered condition, is known as metabolomics and may take advantage of different statistical tools. The new frontier is to adopt self-learning techniques to enhance clustering or classification actions that can improve the predictive power over large amounts of data. Although machine learning is already employed in metabolomics, deep learning and artificial neural networks approaches were only recently successfully applied. In this work, we give an overview of the statistical approaches underlying the wide range of opportunities that machine learning and neural networks allow to perform with accurate metabolites assignment and quantification.Various actual challenges are discussed, such as proper metabolomics, deep learning architectures and model accuracy.