Systems metabolomics: from metabolomic snapshots to design principles

Altered Cellular Metabolism Is a Consequence of Loss of the Ataxia-Linked Protein Sacsin

Mitochondrial dysfunction is implicated in the pathogenesis of the neurological condition autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS), yet precisely how the mitochondrial metabolism is affected is unknown. Thus, to better understand changes in the mitochondrial metabolism caused by loss of the sacsin protein (encoded by the SACS gene, which is mutated in ARSACS), we performed mass spectrometry-based tracer analysis, with both glucose- and glutamine-traced carbon. Comparing the metabolite profiles between wild-type and sacsin-knockout cell lines revealed increased reliance on aerobic glycolysis in sacsin-deficient cells, as evidenced by the increase in lactate and reduction of glucose. Moreover, sacsin knockout cells differentiated towards a neuronal phenotype had increased levels of tricarboxylic acid cycle metabolites relative to the controls. We also observed disruption in the glutaminolysis pathway in differentiated and undifferentiated cells in the absence of sacsin. In conclusion, this work demonstrates consequences for cellular metabolism associated with a loss of sacsin, which may be relevant to ARSACS.

Emerging Biomarkers in Metabolomics: Advancements in Precision Health and Disease Diagnosis

Metabolomics has come to the fore as an efficient tool in the search for biomarkers that are critical for precision health approaches and improved diagnostics. This review will outline recent advances in biomarker discovery based on metabolomics, focusing on metabolomics biomarkers reported in cancer, neurodegenerative disorders, cardiovascular diseases, and metabolic health. In cancer, metabolomics provides evidence for unique oncometabolites that are important for early disease detection and monitoring of treatment responses. Metabolite profiling for conditions such as neurodegenerative and mental health disorders can offer early diagnosis and mechanisms into the disease especially in Alzheimer's and Parkinson's diseases. In addition to these, lipid biomarkers and other metabolites relating to cardiovascular and metabolic disorders are promising for patient stratification and personalized treatment. The gut microbiome and environmental exposure also feature among the influential factors in biomarker discovery because they sculpt individual metabolic profiles, impacting overall health. Further, we discuss technological advances in metabolomics, current clinical applications, and the challenges faced by metabolomics biomarker validation toward precision medicine. Finally, this review discusses future opportunities regarding the integration of metabolomics into routine healthcare to enable preventive and personalized approaches.

The development of machine learning approaches in two-dimensional NMR data interpretation for metabolomics applications

Metabolomics has been widely applied in human diseases and environmental science to study the systematic changes of metabolites over diverse types of stimuli. NMR-based metabolomics has been widely used, but the peak overlap problems in the one-dimensional (1D) NMR spectrum could limit the accuracy of quantitative analysis for metabolomics applications. Two-dimensional (2D) NMR has been applied to solve the 1D NMR overlap problem, but the data processing is still challenging. In this study, we built an automatic approach to process the 2D NMR data for quantitative applications using machine learning approaches. Partial least square discriminant analysis (PLS-DA), artificial neural network classification (ANN-DA), gradient boosted trees classification (XGBoost-DA), and artificial deep learning neural network classification (ANNDL-DA) were applied in combination with an automatic peak selection approach. Standard mixtures, sea anemone extracts, and mouse fecal samples were tested to demonstrate the approach. Our results showed that ANN-DA and ANNDL-DA have high accuracy in selecting 2D NMR peaks (around 90 %), which have a high potential application in 2D NMR-based metabolomics quantitively study, while PLS-DA and XGBoost-DA showed limitations in either data variation or overfitting. Our study built an automatic approach to applying 2D NMR data to routine quantitative analysis in metabolomics.

The HuMet Repository: Watching human metabolism at work

Metabolism oscillates between catabolic and anabolic states depending on food intake, exercise, or stresses that change a multitude of metabolic pathways simultaneously. We present the HuMet Repository for exploring dynamic metabolic responses to oral glucose/lipid loads, mixed meals, 36-h fasting, exercise, and cold stress in healthy subjects. Metabolomics data from blood, urine, and breath of 15 young, healthy men at up to 56 time points are integrated and embedded within an interactive web application, enabling researchers with and without computational expertise to search, visualize, analyze, and contextualize the dynamic metabolite profiles of 2,656 metabolites acquired on multiple platforms. With examples, we demonstrate the utility of the resource for research into the dynamics of human metabolism, highlighting differences and similarities in systemic metabolic responses across challenges and the complementarity of metabolomics platforms. The repository, providing a reference for healthy metabolite changes to six standardized physiological challenges, is freely accessible through a web portal.

Profiling of serum metabolome of breast cancer: multi-cancer features discriminate between healthy women and patients with breast cancer

Introduction

The progression of solid cancers is manifested at the systemic level as molecular changes in the metabolome of body fluids, an emerging source of cancer biomarkers.

Methods

We analyzed quantitatively the serum metabolite profile using high-resolution mass spectrometry. Metabolic profiles were compared between breast cancer patients (n=112) and two groups of healthy women (from Poland and Norway; n=95 and n=112, respectively) with similar age distributions.

Results

Despite differences between both cohorts of controls, a set of 43 metabolites and lipids uniformly discriminated against breast cancer patients and healthy women. Moreover, smaller groups of female patients with other types of solid cancers (colorectal, head and neck, and lung cancers) were analyzed, which revealed a set of 42 metabolites and lipids that uniformly differentiated all three cancer types from both cohorts of healthy women. A common part of both sets, which could be called a multi-cancer signature, contained 23 compounds, which included reduced levels of a few amino acids (alanine, aspartate, glutamine, histidine, phenylalanine, and leucine/isoleucine), lysophosphatidylcholines (exemplified by LPC(18:0)), and diglycerides. Interestingly, a reduced concentration of the most abundant cholesteryl ester (CE(18:2)) typical for other cancers was the least significant in the serum of breast cancer patients. Components present in a multi-cancer signature enabled the establishment of a well-performing breast cancer classifier, which predicted cancer with a very high precision in independent groups of women (AUC>0.95).

Discussion

In conclusion, metabolites critical for discriminating breast cancer patients from controls included components of hypothetical multi-cancer signature, which indicated wider potential applicability of a general serum metabolome cancer biomarker.

Adjusting for false discoveries in constraint-based differential metabolic flux analysis

One of the critical steps to characterize metabolic alterations in multifactorial diseases, as well as their heterogeneity across different patients, is the identification of reactions that exhibit significantly different usage (or flux) between cohorts. However, since metabolic fluxes cannot be determined directly, researchers typically use constraint-based metabolic network models, customized on post-genomics datasets. The use of random sampling within the feasible region of metabolic networks is becoming more prevalent for comparing these networks. While many algorithms have been proposed and compared for efficiently and uniformly sampling the feasible region of metabolic networks, their impact on the risk of making false discoveries when comparing different samples has not been investigated yet, and no sampling strategy has been so far specifically designed to mitigate the problem. To be able to precisely assess the False Discovery Rate (FDR), in this work we compared different samples obtained from the very same metabolic model. We compared the FDR obtained for different model scales, sample sizes, parameters of the sampling algorithm, and strategies to filter out non-significant variations. To be able to compare the largely used hit-and-run strategy with the much less investigated corner-based strategy, we first assessed the intrinsic capability of current corner-based algorithms and of a newly proposed one to visit all vertices of a constraint-based region. We show that false discoveries can occur at high rates even for large samples of small-scale networks. However, we demonstrate that a statistical test based on the empirical null distribution of Kullback-Leibler divergence can effectively correct for false discoveries. We also show that our proposed corner-based algorithm is more efficient than state-of-the-art alternatives and much less prone to false discoveries than hit-and-run strategies. We report that the differences in the marginal distributions obtained with the two strategies are related to but not fully explained by differences in sample standard deviation, as previously thought. Overall, our study provides insights into the impact of sampling strategies on FDR in metabolic network analysis and offers new guidelines for more robust and reproducible analyses.

Screening of an efficient cholesterol-lowering strain of Lactiplantibacillus plantarum 54-1 and investigation of its degradation molecular mechanism

In this study, an efficient cholesterol-lowering strain of Lactiplantibacillus plantarum 54-1 was screened and its degradation molecular mechanism was investigated. Furthermore, a novel practical MRS medium for screening cholesterol-lowering lactic acid bacteria (LAB) was developed based on ultrasound treatment. L. plantarum 54-1 was found to have the highest ability to eliminate cholesterol (340.69 ± 5.87 µg/mL). According to SEM and the count of viable LAB results, the morphology of LAB in the cholesterol-containing medium developed in this experiment was close to the normal (full and smooth), and it can grow normally. Metabolomics revealed that L. plantarum 54-1 initially converted a portion of cholesterol to 7α-hydroxy-cholesterol and then to the key metabolite taurine, via the phosphotransferase system. These metabolites were further transformed into L-alanine, L-lysine, N6-Acetyl-L-lysine, (R)-b-aminoisobutyric acid, and 2-oxoarginine, through glycine, serine, and threonine metabolism, citrate cycle, D-arginine and D-ornithine metabolism, lysine degradation, and pyruvate metabolism pathways. Prokaryotic reference transcriptomics found that this may be mainly regulated by the bsh, phnE, ptsP, B0667_RS04545, and B0667_RSRS12300 genes, which was further validated by qPCR. Furthermore, molecular docking results demonstrated that 8 differential metabolites might bind to another portion of cholesterol via PI-PI conjugation and hydrophobic interactions and lower cholesterol via co-sedimentation. This study has strategic implications for developing probiotic powder food that lowers cholesterol.

Rapid and automated lipid profiling by nuclear magnetic resonance spectroscopy using neural networks

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for quantitative metabolomics; however, quantification of metabolites from NMR data is often a slow and tedious process requiring user input and expertise. In this study, we propose a neural network approach for rapid, automated lipid identification and quantification from NMR data. Multilayered perceptron (MLP) networks were developed with NMR spectra as the input and lipid concentrations as output. Three large synthetic datasets were generated, each with 55,000 spectra from an original 30 scans of reference standards, by using linear combinations of standards and simulating experimental-like modifications (line broadening, noise, peak shifts, baseline shifts) and common interference signals (water, tetramethylsilane, extraction solvent), and were used to train MLPs for robust prediction of lipid concentrations. The performances of MLPS were first validated on various synthetic datasets to assess the effect of incorporating different modifications on their accuracy. The MLPs were then evaluated on experimentally acquired data from complex lipid mixtures. The MLP-derived lipid concentrations showed high correlations and slopes close to unity for most of the quantified lipid metabolites in experimental mixtures compared with ground-truth concentrations. The most accurate, robust MLP was used to profile lipids in lipophilic hepatic extracts from a rat metabolomics study. The MLP lipid results analyzed by two-way ANOVA for dietary and sex differences were similar to those obtained with a conventional NMR quantification method. In conclusion, this study demonstrates the potential and feasibility of a neural network approach for improving speed and automation in NMR lipid profiling and this approach can be easily tailored to other quantitative, targeted spectroscopic analyses in academia or industry.

Global and Spatial Metabolomics of Individual Cells Using a Tapered Pneumatically Assisted nano-DESI Probe

Single-cell metabolomics has the potential to reveal unique insights into intracellular mechanisms and biological processes. However, the detection of metabolites from individual cells is challenging due to their versatile chemical properties and concentrations. Here, we demonstrate a tapered probe for pneumatically assisted nanospray desorption electrospray ionization (PA nano-DESI) mass spectrometry that enables both chemical imaging of larger cells and global metabolomics of smaller 15 μm cells. Additionally, by depositing cells in predefined arrays, we show successful metabolomics from three individual INS-1 cells per minute, which enabled the acquisition of data from 479 individual cells. Several cells were used to optimize analytical conditions, and 93 or 97 cells were used to monitor metabolome alterations in INS-1 cells after exposure to a low or high glucose concentration, respectively. Our analytical approach offers insights into cellular heterogeneity and provides valuable information about cellular processes and responses in individual cells.

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.

The Expectation and Reality of the HepG2 Core Metabolic Profile

To represent the composition of small molecules circulating in HepG2 cells and the formation of the "core" of characteristic metabolites that often attract researchers' attention, we conducted a meta-analysis of 56 datasets obtained through metabolomic profiling via mass spectrometry and NMR. We highlighted the 288 most commonly studied compounds of diverse chemical nature and analyzed metabolic processes involving these small molecules. Building a complete map of the metabolome of a cell, which encompasses the diversity of possible impacts on it, is a severe challenge for the scientific community, which is faced not only with natural limitations of experimental technologies, but also with the absence of transparent and widely accepted standards for processing and presenting the obtained metabolomic data. Formulating our research design, we aimed to reveal metabolites crucial to the Hepg2 cell line, regardless of all chemical and/or physical impact factors. Unfortunately, the existing paradigm of data policy leads to a streetlight effect. When analyzing and reporting only target metabolites of interest, the community ignores the changes in the metabolomic landscape that hide many molecular secrets.

Tumor heterogeneity: preclinical models, emerging technologies, and future applications

Heterogeneity describes the differences among cancer cells within and between tumors. It refers to cancer cells describing variations in morphology, transcriptional profiles, metabolism, and metastatic potential. More recently, the field has included the characterization of the tumor immune microenvironment and the depiction of the dynamics underlying the cellular interactions promoting the tumor ecosystem evolution. Heterogeneity has been found in most tumors representing one of the most challenging behaviors in cancer ecosystems. As one of the critical factors impairing the long-term efficacy of solid tumor therapy, heterogeneity leads to tumor resistance, more aggressive metastasizing, and recurrence. We review the role of the main models and the emerging single-cell and spatial genomic technologies in our understanding of tumor heterogeneity, its contribution to lethal cancer outcomes, and the physiological challenges to consider in designing cancer therapies. We highlight how tumor cells dynamically evolve because of the interactions within the tumor immune microenvironment and how to leverage this to unleash immune recognition through immunotherapy. A multidisciplinary approach grounded in novel bioinformatic and computational tools will allow reaching the integrated, multilayered knowledge of tumor heterogeneity required to implement personalized, more efficient therapies urgently required for cancer patients.

Closing the gap between in vivo and in vitro omics: using QA/QC to strengthen ex vivo NMR metabolomics

Metabolomics aims to achieve a global quantitation of the pool of metabolites within a biological system. Importantly, metabolite concentrations serve as a sensitive marker of both genomic and phenotypic changes in response to both internal and external stimuli. NMR spectroscopy greatly aids in the understanding of both in vitro and in vivo physiological systems and in the identification of diagnostic and therapeutic biomarkers. Accordingly, NMR is widely utilized in metabolomics and fluxomics studies due to its limited requirements for sample preparation and chromatography, its non-destructive and quantitative nature, its utility in the structural elucidation of unknown compounds, and, importantly, its versatility in the analysis of in vitro, in vivo, and ex vivo samples. This review provides an overview of the strengths and limitations of in vitro and in vivo experiments for translational research and discusses how ex vivo studies may overcome these weaknesses to facilitate the extrapolation of in vitro insights to an in vivo system. The application of NMR-based metabolomics to ex vivo samples, tissues, and biofluids can provide essential information that is close to a living system (in vivo) with sensitivity and resolution comparable to those of in vitro studies. The success of this extrapolation process is critically dependent on high-quality and reproducible data. Thus, the incorporation of robust quality assurance and quality control checks into the experimental design and execution of NMR-based metabolomics experiments will ensure the successful extrapolation of ex vivo studies to benefit translational medicine.

SMART-Miner: A convolutional neural network-based metabolite identification from 1 H-13 C HSQC spectra

The identification of metabolites from complex biofluids and extracts of tissues is an essential process for understanding metabolic profiles. Nuclear magnetic resonance (NMR) spectroscopy is widely used in metabolomics studies for identification and quantification of metabolites. However, the accurate identification of individual metabolites is still a challenging process with higher peak intensity or similar chemical shifts from different metabolites. In this study, we applied a convolutional neural network (CNN) to ¹ H-¹³ C HSQC NMR spectra to achieve accurate peak identification in complex mixtures. The results reveal that the neural network was successfully trained on metabolite identification from these 2D NMR spectra and achieved very good performance compared with other NMR-based metabolomic tools.

Combining denoising of RNA-seq data and flux balance analysis for cluster analysis of single cells

Background

Sophisticated methods to properly pre-process and analyze the increasing collection of single-cell RNA sequencing (scRNA-seq) data are increasingly being developed. On the contrary, the best practices to integrate these data into metabolic networks, aiming at describing metabolic phenotypes within a heterogeneous cell population, have been poorly investigated. In this regard, a critical factor is the presence of false zero values in reactions essential for a fundamental metabolic function, such as biomass or energy production. Here, we investigate the role of denoising strategies in mitigating this problem.

Methods

We applied state-of-the-art denoising strategies - namely MAGIC, ENHANCE, and SAVER - on three public scRNA-seq datasets. We then associated a metabolic flux distribution with every single cell by embedding its noise-free transcriptomics profile in the constraints of the optimization of a core metabolic model. Finally, we used the obtained single-cell optimal metabolic fluxes as features for cluster analysis. We compared the results obtained with different techniques, and with or without the use of denoising. We also investigated the possibility of applying denoising directly on the Reaction Activity Scores, which are metabolic features extracted from the read counts, rather than on the read counts.

Results

We show that denoising of transcriptomics data improves the clustering of single cells. We also illustrate that denoising restores important metabolic properties, such as the correlation between cell cycle phase and biomass accumulation, and between the RAS scores of reactions belonging to the same metabolic pathway. We show that MAGIC performs better than ENHANCE and SAVER, and that, denoising applied directly on the RAS matrix could be an effective alternative in removing false zero values from essential metabolic reactions.

Conclusions

Our results indicate that including denoising as a pre-processing operation represents a milestone to integrate scRNA-seq data into Flux Balance Analysis simulations and to perform single-cell cluster analysis with a focus on metabolic phenotypes.

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.

An Optimized Workflow for the Analysis of Metabolic Fluxes in Cancer Spheroids Using Seahorse Technology

Three-dimensional cancer models, such as spheroids, are increasingly being used to study cancer metabolism because they can better recapitulate the molecular and physiological aspects of the tumor architecture than conventional monolayer cultures. Although Agilent Seahorse XFe96 (Agilent Technologies, Santa Clara, CA, United States) is a valuable technology for studying metabolic alterations occurring in cancer cells, its application to three-dimensional cultures is still poorly optimized. We present a reliable and reproducible workflow for the Seahorse metabolic analysis of three-dimensional cultures. An optimized protocol enables the formation of spheroids highly regular in shape and homogenous in size, reducing variability in metabolic parameters among the experimental replicates, both under basal and drug treatment conditions. High-resolution imaging allows the calculation of the number of viable cells in each spheroid, the normalization of metabolic parameters on a per-cell basis, and grouping of the spheroids as a function of their size. Multivariate statistical tests on metabolic parameters determined by the Mito Stress test on two breast cancer cell lines show that metabolic differences among the studied spheroids are mostly related to the cell line rather than to the size of the spheroid. The optimized workflow allows high-resolution metabolic characterization of three-dimensional cultures, their comparison with monolayer cultures, and may aid in the design and interpretation of (multi)drug protocols.

INTEGRATE: Model-based multi-omics data integration to characterize multi-level metabolic regulation

Metabolism is directly and indirectly fine-tuned by a complex web of interacting regulatory mechanisms that fall into two major classes. On the one hand, the expression level of the catalyzing enzyme sets the maximal theoretical flux level (i.e., the net rate of the reaction) for each enzyme-controlled reaction. On the other hand, metabolic regulation controls the metabolic flux through the interactions of metabolites (substrates, cofactors, allosteric modulators) with the responsible enzyme. High-throughput data, such as metabolomics and transcriptomics data, if analyzed separately, do not accurately characterize the hierarchical regulation of metabolism outlined above. They must be integrated to disassemble the interdependence between different regulatory layers controlling metabolism. To this aim, we propose INTEGRATE, a computational pipeline that integrates metabolomics and transcriptomics data, using constraint-based stoichiometric metabolic models as a scaffold. We compute differential reaction expression from transcriptomics data and use constraint-based modeling to predict if the differential expression of metabolic enzymes directly originates differences in metabolic fluxes. In parallel, we use metabolomics to predict how differences in substrate availability translate into differences in metabolic fluxes. We discriminate fluxes regulated at the metabolic and/or gene expression level by intersecting these two output datasets. We demonstrate the pipeline using a set of immortalized normal and cancer breast cell lines. In a clinical setting, knowing the regulatory level at which a given metabolic reaction is controlled will be valuable to inform targeted, truly personalized therapies in cancer patients.

Human metabolite detection by surface-enhanced Raman spectroscopy

Metabolites are important biomarkers in human body fluids, conveying direct information of cellular activities and physical conditions. Metabolite detection has long been a research hotspot in the field of biology and medicine. Surface-enhanced Raman spectroscopy (SERS), based on the molecular “fingerprint” of Raman spectrum and the enormous signal enhancement (down to a single-molecule level) by plasmonic nanomaterials, has proven to be a novel and powerful tool for metabolite detection. SERS provides favorable properties such as ultra-sensitive, label-free, rapid, specific, and non-destructive detection processes. In this review, we summarized the progress in recent 10 years on SERS-based sensing of endogenous metabolites at the cellular level, in tissues, and in biofluids, as well as drug metabolites in biofluids. We made detailed discussions on the challenges and optimization methods of SERS technique in metabolite detection. The combination of SERS with modern biomedical technology were also anticipated.

Simulating single-cell metabolism using a stochastic flux-balance analysis algorithm

Stochasticity from gene expression in single cells is known to drive metabolic heterogeneity at the level of cellular populations, which is understood to have important consequences for issues such as microbial drug tolerance and treatment of human diseases like cancer. Despite considerable advancements in profiling the genomes, transcriptomes, and proteomes of single cells, it remains difficult to experimentally characterize their metabolism at the genome scale. Computational methods could bridge this gap toward a systems understanding of single-cell biology. To address this challenge, we developed stochastic simulation algorithm with flux-balance analysis embedded (SSA-FBA), a computational framework for simulating the stochastic dynamics of the metabolism of individual cells using genome-scale metabolic models with experimental estimates of gene expression and enzymatic reaction rate parameters. SSA-FBA extends the constraint-based modeling formalism of metabolic network modeling to the single-cell regime, enabling simulation when experimentation is intractable. We also developed an efficient implementation of SSA-FBA that leverages the topology of embedded flux-balance analysis models to significantly reduce the computational cost of simulation. As a preliminary case study, we built a reduced single-cell model of Mycoplasma pneumoniae and used SSA-FBA to illustrate the role of stochasticity on the dynamics of metabolism at the single-cell level.

Multi-Omics Comparison of the Spontaneous Diabetes Mellitus and Diet-Induced Prediabetic Macaque Models

The prevalence of diabetes mellitus has been increasing for decades worldwide. To develop safe and potent therapeutics, animal models contribute a lot to the studies of the mechanisms underlying its pathogenesis. Dietary induction using is a well-accepted protocol in generating insulin resistance and diabetes models. In the present study, we reported the multi-omics profiling of the liver and sera from both peripheral blood and hepatic portal vein blood from Macaca fascicularis that spontaneously developed Type-2 diabetes mellitus with a chow diet (sDM). The other two groups of the monkeys fed with chow diet and high-fat high-sugar (HFHS) diet, respectively, were included for comparison. Analyses of various omics datasets revealed the alterations of high consistency. Between the sDM and HFHS monkeys, both the similar and unique alterations in the lipid metabolism have been demonstrated from metabolomic, transcriptomic, and proteomic data repeatedly. The comparison of the proteome and transcriptome confirmed the involvement of fatty acid binding protein 4 (FABP4) in the diet-induced pathogenesis of diabetes in macaques. Furthermore, the commonly changed genes between spontaneous diabetes and HFHS diet-induced prediabetes suggested that the alterations in the intra- and extracellular structural proteins and cell migration in the liver might mediate the HFHS diet induction of diabetes mellitus.

Transcriptomics and Metabolomics Integration Reveals Redox-Dependent Metabolic Rewiring in Breast Cancer Cells

Rewiring glucose metabolism toward aerobic glycolysis provides cancer cells with a rapid generation of pyruvate, ATP, and NADH, while pyruvate oxidation to lactate guarantees refueling of oxidized NAD⁺ to sustain glycolysis. CtPB2, an NADH-dependent transcriptional co-regulator, has been proposed to work as an NADH sensor, linking metabolism to epigenetic transcriptional reprogramming. By integrating metabolomics and transcriptomics in a triple-negative human breast cancer cell line, we show that genetic and pharmacological down-regulation of CtBP2 strongly reduces cell proliferation by modulating the redox balance, nucleotide synthesis, ROS generation, and scavenging. Our data highlight the critical role of NADH in controlling the oncogene-dependent crosstalk between metabolism and the epigenetically mediated transcriptional program that sustains energetic and anabolic demands in cancer cells.

Metabolomics-Guided Elucidation of Plant Abiotic Stress Responses in the 4IR Era: An Overview

Plants are constantly challenged by changing environmental conditions that include abiotic stresses. These are limiting their development and productivity and are subsequently threatening our food security, especially when considering the pressure of the increasing global population. Thus, there is an urgent need for the next generation of crops with high productivity and resilience to climate change. The dawn of a new era characterized by the emergence of fourth industrial revolution (4IR) technologies has redefined the ideological boundaries of research and applications in plant sciences. Recent technological advances and machine learning (ML)-based computational tools and omics data analysis approaches are allowing scientists to derive comprehensive metabolic descriptions and models for the target plant species under specific conditions. Such accurate metabolic descriptions are imperatively essential for devising a roadmap for the next generation of crops that are resilient to environmental deterioration. By synthesizing the recent literature and collating data on metabolomics studies on plant responses to abiotic stresses, in the context of the 4IR era, we point out the opportunities and challenges offered by omics science, analytical intelligence, computational tools and big data analytics. Specifically, we highlight technological advancements in (plant) metabolomics workflows and the use of machine learning and computational tools to decipher the dynamics in the chemical space that define plant responses to abiotic stress conditions.

The remodel of the "central dogma": a metabolomics interaction perspective

BACKGROUND

In 1957, Francis Crick drew a linear diagram on a blackboard. This diagram is often called the "central dogma." Subsequently, the relationships between different steps of the "central dogma" have been shown to be considerably complex, mostly because of the emerging world of small molecules. It is noteworthy that metabolites can be generated from the diet through gut microbiome metabolism, serve as substrates for epigenetic modifications, destabilize DNA quadruplexes, and follow Lamarckian inheritance. Small molecules were once considered the missing link in the "central dogma"; however, recently they have acquired a central role, and their general perception as downstream products has become reductionist. Metabolomics is a large-scale analysis of metabolites, and this emerging field has been shown to be the closest omics associated with the phenotype and concomitantly, the basis for all omics.

AIM OF REVIEW

Herein, we propose a broad updated perspective for the flux of information diagram centered in metabolomics, including the influence of other factors, such as epigenomics, diet, nutrition, and the gut- microbiome.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Metabolites are the beginning and the end of the flux of information.

Accelerated global sensitivity analysis of genome-wide constraint-based metabolic models

Background

Genome-wide reconstructions of metabolism opened the way to thorough investigations of cell metabolism for health care and industrial purposes. However, the predictions offered by Flux Balance Analysis (FBA) can be strongly affected by the choice of flux boundaries, with particular regard to the flux of reactions that sink nutrients into the system. To mitigate possible errors introduced by a poor selection of such boundaries, a rational approach suggests to focus the modeling efforts on the pivotal ones.

Methods

In this work, we present a methodology for the automatic identification of the key fluxes in genome-wide constraint-based models, by means of variance-based sensitivity analysis. The goal is to identify the parameters for which a small perturbation entails a large variation of the model outcomes, also referred to as sensitive parameters. Due to the high number of FBA simulations that are necessary to assess sensitivity coefficients on genome-wide models, our method exploits a master-slave methodology that distributes the computation on massively multi-core architectures. We performed the following steps: (1) we determined the putative parameterizations of the genome-wide metabolic constraint-based model, using Saltelli’s method; (2) we applied FBA to each parameterized model, distributing the massive amount of calculations over multiple nodes by means of MPI; (3) we then recollected and exploited the results of all FBA runs to assess a global sensitivity analysis.

Results

We show a proof-of-concept of our approach on latest genome-wide reconstructions of human metabolism Recon2.2 and Recon3D. We report that most sensitive parameters are mainly associated with the intake of essential amino acids in Recon2.2, whereas in Recon 3D they are associated largely with phospholipids. We also illustrate that in most cases there is a significant contribution of higher order effects.

Conclusion

Our results indicate that interaction effects between different model parameters exist, which should be taken into account especially at the stage of calibration of genome-wide models, supporting the importance of a global strategy of sensitivity analysis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04002-0.

The Glioblastoma Microenvironment: Morphology, Metabolism, and Molecular Signature of Glial Dynamics to Discover Metabolic Rewiring Sequence

Different functional states determine glioblastoma (GBM) heterogeneity. Brain cancer cells coexist with the glial cells in a functional syncytium based on a continuous metabolic rewiring. However, standard glioma therapies do not account for the effects of the glial cells within the tumor microenvironment. This may be a possible reason for the lack of improvements in patients with high-grade gliomas therapies. Cell metabolism and bioenergetic fitness depend on the availability of nutrients and interactions in the microenvironment. It is strictly related to the cell location in the tumor mass, proximity to blood vessels, biochemical gradients, and tumor evolution, underlying the influence of the context and the timeline in anti-tumor therapeutic approaches. Besides the cancer metabolic strategies, here we review the modifications found in the GBM-associated glia, focusing on morphological, molecular, and metabolic features. We propose to analyze the GBM metabolic rewiring processes from a systems biology perspective. We aim at defining the crosstalk between GBM and the glial cells as modules. The complex networking may be expressed by metabolic modules corresponding to the GBM growth and spreading phases. Variation in the oxidative phosphorylation (OXPHOS) rate and regulation appears to be the most important part of the metabolic and functional heterogeneity, correlating with glycolysis and response to hypoxia. Integrated metabolic modules along with molecular and morphological features could allow the identification of key factors for controlling the GBM-stroma metabolism in multi-targeted, time-dependent therapies.

On the Use of Topological Features of Metabolic Networks for the Classification of Cancer Samples

BACKGROUND

The increasing availability of omics data collected from patients affected by severe pathologies, such as cancer, is fostering the development of data science methods for their analysis.

INTRODUCTION

The combination of data integration and machine learning approaches can provide new powerful instruments to tackle the complexity of cancer development and deliver effective diagnostic and prognostic strategies.

METHODS

We explore the possibility of exploiting the topological properties of sample-specific metabolic networks as features in a supervised classification task. Such networks are obtained by projecting transcriptomic data from RNA-seq experiments on genome-wide metabolic models to define weighted networks modeling the overall metabolic activity of a given sample.

RESULTS

We show the classification results on a labeled breast cancer dataset from the TCGA database, including 210 samples (cancer vs. normal). In particular, we investigate how the performance is affected by a threshold-based pruning of the networks by comparing Artificial Neural Networks, Support Vector Machines and Random Forests. Interestingly, the best classification performance is achieved within a small threshold range for all methods, suggesting that it might represent an effective choice to recover useful information while filtering out noise from data. Overall, the best accuracy is achieved with SVMs, which exhibit performances similar to those obtained when gene expression profiles are used as features.

CONCLUSION

These findings demonstrate that the topological properties of sample-specific metabolic networks are effective in classifying cancer and normal samples, suggesting that useful information can be extracted from a relatively limited number of features.

Dataset on the Effects of Different Pre-Harvest Factors on the Metabolomics Profile of Lettuce (Lactuca sativa L.) Leaves

The study of the relationship between cultivated plants and environmental factors can provide information ranging from a deeper understanding of the plant biological system to the development of more effective management strategies for improving yield, quality, and sustainability of the produce. In this article, we present a comprehensive metabolomics dataset of two phytochemically divergent lettuce (Lactuca sativa L.) butterhead varieties under different growing conditions. Plants were cultivated in hydroponics in a growth chamber with ambient control. The pre-harvest factors that were independently investigated were light intensity (two levels), the ionic strength of the nutrient solutions (three levels), and the molar ratio of three macroelements (K, Mg, and Ca) in the nutrient solution (three levels). We used an untargeted, mass-spectrometry-based approach to characterize the metabolomics profiles of leaves harvested 19 days after transplant. The data revealed the ample impact on both primary and secondary metabolism and its range of variation. Moreover, our dataset is useful for uncovering the complex effects of the genotype, the environmental factor(s), and their interaction, which may deserve further investigation.

Profiling and Targeting of Energy and Redox Metabolism in Grade 2 Bladder Cancer Cells with Different Invasiveness Properties

Bladder cancer is one of the most prevalent deadly diseases worldwide. Grade 2 tumors represent a good window of therapeutic intervention, whose optimization requires high resolution biomarker identification. Here we characterize energy metabolism and cellular properties associated with spreading and tumor progression of RT112 and 5637, two Grade 2 cancer cell lines derived from human bladder, representative of luminal-like and basal-like tumors, respectively. The two cell lines have similar proliferation rates, but only 5637 cells show efficient lateral migration. In contrast, RT112 cells are more prone to form spheroids. RT112 cells produce more ATP by glycolysis and OXPHOS, present overall higher metabolic plasticity and are less sensitive than 5637 to nutritional perturbation of cell proliferation and migration induced by treatment with 2-deoxyglucose and metformin. On the contrary, spheroid formation is less sensitive to metabolic perturbations in 5637 than RT112 cells. The ability of metformin to reduce, although with different efficiency, cell proliferation, sphere formation and migration in both cell lines, suggests that OXPHOS targeting could be an effective strategy to reduce the invasiveness of Grade 2 bladder cancer cells.

Disruption of redox homeostasis for combinatorial drug efficacy in K-Ras tumors as revealed by metabolic connectivity profiling

Abstract

Background

Rewiring of metabolism induced by oncogenic K-Ras in cancer cells involves both glucose and glutamine utilization sustaining enhanced, unrestricted growth. The development of effective anti-cancer treatments targeting metabolism may be facilitated by the identification and rational combinatorial targeting of metabolic pathways.

Methods

We performed mass spectrometric metabolomics analysis in vitro and in vivo experiments to evaluate the efficacy of drugs and identify metabolic connectivity.

Results

We show that K-Ras-mutant lung and colon cancer cells exhibit a distinct metabolic rewiring, the latter being more dependent on respiration. Combined treatment with the glutaminase inhibitor CB-839 and the PI3K/aldolase inhibitor NVP-BKM120 more consistently reduces cell growth of tumor xenografts. Maximal growth inhibition correlates with the disruption of redox homeostasis, involving loss of reduced glutathione regeneration, redox cofactors, and a decreased connectivity among metabolites primarily involved in nucleic acid metabolism.

Conclusions

Our findings open the way to develop metabolic connectivity profiling as a tool for a selective strategy of combined drug repositioning in precision oncology.

Metabolomic Approaches to Study Chemical Exposure-Related Metabolism Alterations in Mammalian Cell Cultures

Biological organisms are constantly exposed to an immense repertoire of molecules that cover environmental or food-derived molecules and drugs, triggering a continuous flow of stimuli-dependent adaptations. The diversity of these chemicals as well as their concentrations contribute to the multiplicity of induced effects, including activation, stimulation, or inhibition of physiological processes and toxicity. Metabolism, as the foremost phenotype and manifestation of life, has proven to be immensely sensitive and highly adaptive to chemical stimuli. Therefore, studying the effect of endo- or xenobiotics over cellular metabolism delivers valuable knowledge to apprehend potential cellular activity of individual molecules and evaluate their acute or chronic benefits and toxicity. The development of modern metabolomics technologies such as mass spectrometry or nuclear magnetic resonance spectroscopy now offers unprecedented solutions for the rapid and efficient determination of metabolic profiles of cells and more complex biological systems. Combined with the availability of well-established cell culture techniques, these analytical methods appear perfectly suited to determine the biological activity and estimate the positive and negative effects of chemicals in a variety of cell types and models, even at hardly detectable concentrations. Metabolic phenotypes can be estimated from studying intracellular metabolites at homeostasis in vivo, while in vitro cell cultures provide additional access to metabolites exchanged with growth media. This article discusses analytical solutions available for metabolic phenotyping of cell culture metabolism as well as the general metabolomics workflow suitable for testing the biological activity of molecular compounds. We emphasize how metabolic profiling of cell supernatants and intracellular extracts can deliver valuable and complementary insights for evaluating the effects of xenobiotics on cellular metabolism. We note that the concepts and methods discussed primarily for xenobiotics exposure are widely applicable to drug testing in general, including endobiotics that cover active metabolites, nutrients, peptides and proteins, cytokines, hormones, vitamins, etc.

An Optimized Dual Extraction Method for the Simultaneous and Accurate Analysis of Polar Metabolites and Lipids Carried out on Single Biological Samples

The functional understanding of metabolic changes requires both a significant investigation into metabolic pathways, as enabled by global metabolomics and lipidomics approaches, and the comprehensive and accurate exploration of specific key pathways. To answer this pivotal challenge, we propose an optimized approach, which combines an efficient sample preparation, aiming to reduce the variability, with a biphasic extraction method, where both the aqueous and organic phases of the same sample are used for mass spectrometry analyses. We demonstrated that this double extraction protocol allows working with one single sample without decreasing the metabolome and lipidome coverage. It enables the targeted analysis of 40 polar metabolites and 82 lipids, together with the absolute quantification of 32 polar metabolites, providing comprehensive coverage and quantitative measurement of the metabolites involved in central carbon energy pathways. With this method, we evidenced modulations of several lipids, amino acids, and energy metabolites in HepaRG cells exposed to fenofibrate, a model hepatic toxicant, and metabolic modulator. This new protocol is particularly relevant for experiments involving limited amounts of biological material and for functional metabolic explorations and is thus of particular interest for studies aiming to decipher the effects and modes of action of metabolic disrupting compounds.

Metabolomic Biomarkers for Detection, Prognosis and Identifying Recurrence in Endometrial Cancer

Metabolic reprogramming is increasingly recognised as one of the defining hallmarks of tumorigenesis. There is compelling evidence to suggest that endometrial cancer develops and progresses in the context of profound metabolic dysfunction. Whilst the incidence of endometrial cancer continues to rise in parallel with the global epidemic of obesity, there are, as yet, no validated biomarkers that can aid risk prediction, early detection, prognostic evaluation or surveillance. Advances in high-throughput technologies have, in recent times, shown promise for biomarker discovery based on genomic, transcriptomic, proteomic and metabolomic platforms. Metabolomics, the large-scale study of metabolites, deals with the downstream products of the other omics technologies and thus best reflects the human phenotype. This review aims to provide a summary and critical synthesis of the existing literature with the ultimate goal of identifying the most promising metabolite biomarkers that can augment current endometrial cancer diagnostic, prognostic and recurrence surveillance strategies. Identified metabolites and their biochemical pathways are discussed in the context of what we know about endometrial carcinogenesis and their potential clinical utility is evaluated. Finally, we underscore the challenges inherent in metabolomic biomarker discovery and validation and provide fresh perspectives and directions for future endometrial cancer biomarker research.