Toward modeling metabolic state from single-cell transcriptomics

Single-cell omics analysis with genome-scale metabolic modeling

Single-cell technologies have been widely used in biological studies and generated a plethora of single-cell data to be interpreted. Due to the inclusion of the priori metabolic network knowledge as well as gene-protein-reaction associations, genome-scale metabolic models (GEMs) have been a powerful tool to integrate and thereby interpret various omics data mostly from bulk samples. Here, we first review two common ways to leverage bulk omics data with GEMs and then discuss advances on integrative analysis of single-cell omics data with GEMs. We end by presenting our views on current challenges and perspectives in this field.

Resource allocation in mammalian systems

Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Metabolic heterogeneity in cancer

Cancer cells rewire their metabolism to survive during cancer progression. In this context, tumour metabolic heterogeneity arises and develops in response to diverse environmental factors. This metabolic heterogeneity contributes to cancer aggressiveness and impacts therapeutic opportunities. In recent years, technical advances allowed direct characterisation of metabolic heterogeneity in tumours. In addition to the metabolic heterogeneity observed in primary tumours, metabolic heterogeneity temporally evolves along with tumour progression. In this Review, we summarize the mechanisms of environment-induced metabolic heterogeneity. In addition, we discuss how cancer metabolism and the key metabolites and enzymes temporally and functionally evolve during the metastatic cascade and treatment.

Unveiling the future of metabolic medicine: omics technologies driving personalized solutions for precision treatment of metabolic disorders

Metabolic disorders are increasingly prevalent worldwide, leading to high rates of morbidity and mortality. The variety of metabolic illnesses can be addressed through personalized medicine. The goal of personalized medicine is to give doctors the ability to anticipate the best course of treatment for patients with metabolic problems. By analyzing a patient's metabolomic, proteomic, genetic profile, and clinical data, physicians can identify relevant diagnostic, and predictive biomarkers and develop treatment plans and therapy for acute and chronic metabolic diseases. To achieve this goal, real-time modeling of clinical data and multiple omics is essential to pinpoint underlying biological mechanisms, risk factors, and possibly useful data to promote early diagnosis and prevention of complex diseases. Incorporating cutting-edge technologies like artificial intelligence and machine learning is crucial for consolidating diverse forms of data, examining multiple variables, establishing databases of clinical indicators to aid decision-making, and formulating ethical protocols to address concerns. This review article aims to explore the potential of personalized medicine utilizing omics approaches for the treatment of metabolic disorders. It focuses on the recent advancements in genomics, epigenomics, proteomics, metabolomics, and nutrigenomics, emphasizing their role in revolutionizing personalized medicine.

Enablers and challenges of spatial omics, a melting pot of technologies

Spatial omics has emerged as a rapidly growing and fruitful field with hundreds of publications presenting novel methods for obtaining spatially resolved information for any omics data type on spatial scales ranging from subcellular to organismal. From a technology development perspective, spatial omics is a highly interdisciplinary field that integrates imaging and omics, spatial and molecular analyses, sequencing and mass spectrometry, and image analysis and bioinformatics. The emergence of this field has not only opened a window into spatial biology, but also created multiple novel opportunities, questions, and challenges for method developers. Here, we provide the perspective of technology developers on what makes the spatial omics field unique. After providing a brief overview of the state of the art, we discuss technological enablers and challenges and present our vision about the future applications and impact of this melting pot.

FLUXestimator: a webserver for predicting metabolic flux and variations using transcriptomics data

Quantitative assessment of single cell fluxome is critical for understanding the metabolic heterogeneity in diseases. Unfortunately, laboratory-based single cell fluxomics is currently impractical, and the current computational tools for flux estimation are not designed for single cell-level prediction. Given the well-established link between transcriptomic and metabolomic profiles, leveraging single cell transcriptomics data to predict single cell fluxome is not only feasible but also an urgent task. In this study, we present FLUXestimator, an online platform for predicting metabolic fluxome and variations using single cell or general transcriptomics data of large sample-size. The FLUXestimator webserver implements a recently developed unsupervised approach called single cell flux estimation analysis (scFEA), which uses a new neural network architecture to estimate reaction rates from transcriptomics data. To the best of our knowledge, FLUXestimator is the first web-based tool dedicated to predicting cell-/sample-wise metabolic flux and metabolite variations using transcriptomics data of human, mouse and 15 other common experimental organisms. The FLUXestimator webserver is available at http://scFLUX.org/, and stand-alone tools for local use are available at https://github.com/changwn/scFEA. Our tool provides a new avenue for studying metabolic heterogeneity in diseases and has the potential to facilitate the development of new therapeutic strategies.

New Insights and Potential Therapeutic Interventions in Metabolic Diseases

Endocrine homeostasis and metabolic diseases have been the subject of extensive research in recent years. The development of new techniques and insights has led to a deeper understanding of the mechanisms underlying these conditions and opened up new avenues for diagnosis and treatment. In this review, we discussed the rise of metabolic diseases, especially in Western countries, the genetical, psychological, and behavioral basis of metabolic diseases, the role of nutrition and physical activity in the development of metabolic diseases, the role of single-cell transcriptomics, gut microbiota, epigenetics, advanced imaging techniques, and cell-based therapies in metabolic diseases. Finally, practical applications derived from this information are made.

Systems Immunology Approaches to Metabolism

Over the last decade, immunometabolism has emerged as a novel interdisciplinary field of research and yielded significant fundamental insights into the regulation of immune responses. Multiple classical approaches to interrogate immunometabolism, including bulk metabolic profiling and analysis of metabolic regulator expression, paved the way to appreciating the physiological complexity of immunometabolic regulation in vivo. Studying immunometabolism at the systems level raised the need to transition towards the next-generation technology for metabolic profiling and analysis. Spatially resolved metabolic imaging and computational algorithms for multi-modal data integration are new approaches to connecting metabolism and immunity. In this review, we discuss recent studies that highlight the complex physiological interplay between immune responses and metabolism and give an overview of technological developments that bear the promise of capturing this complexity most directly and comprehensively.

To metabolomics and beyond: a technological portfolio to investigate cancer metabolism

Tumour cells have exquisite flexibility in reprogramming their metabolism in order to support tumour initiation, progression, metastasis and resistance to therapies. These reprogrammed activities include a complete rewiring of the bioenergetic, biosynthetic and redox status to sustain the increased energetic demand of the cells. Over the last decades, the cancer metabolism field has seen an explosion of new biochemical technologies giving more tools than ever before to navigate this complexity. Within a cell or a tissue, the metabolites constitute the direct signature of the molecular phenotype and thus their profiling has concrete clinical applications in oncology. Metabolomics and fluxomics, are key technological approaches that mainly revolutionized the field enabling researchers to have both a qualitative and mechanistic model of the biochemical activities in cancer. Furthermore, the upgrade from bulk to single-cell analysis technologies provided unprecedented opportunity to investigate cancer biology at cellular resolution allowing an in depth quantitative analysis of complex and heterogenous diseases. More recently, the advent of functional genomic screening allowed the identification of molecular pathways, cellular processes, biomarkers and novel therapeutic targets that in concert with other technologies allow patient stratification and identification of new treatment regimens. This review is intended to be a guide for researchers to cancer metabolism, highlighting current and emerging technologies, emphasizing advantages, disadvantages and applications with the potential of leading the development of innovative anti-cancer therapies.

Single-Cell Kinetic Modeling of β-Lapachone Metabolism in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) cells are highly heterogeneous in their metabolism and typically experience elevated reactive oxygen species (ROS) levels such as superoxide and hydrogen peroxide (H2O2) in the tumor microenvironment. Tumor cells survive under these chronic oxidative conditions by upregulating antioxidant systems. To investigate the heterogeneity of cellular responses to chemotherapeutic H2O2 generation in tumor and healthy tissue, we leveraged single-cell RNA-sequencing (scRNA-seq) data to perform redox systems-level simulations of quinone-cycling β-lapachone treatment as a source of NQO1-dependent rapid superoxide and hydrogen peroxide (H2O2) production. Transcriptomic data from 10 HNSCC patient tumors was used to populate over 4000 single-cell antioxidant enzymatic network models of drug metabolism. The simulations reflected significant systems-level differences between the redox states of healthy and cancer cells, demonstrating in some patient samples a targetable cancer cell population or in others statistically indistinguishable effects between non-malignant and malignant cells. Subsequent multivariate analyses between healthy and malignant cellular models pointed to distinct contributors of redox responses between these phenotypes. This model framework provides a mechanistic basis for explaining mixed outcomes of NAD(P)H:quinone oxidoreductase 1 (NQO1)-bioactivatable therapeutics despite the tumor specificity of these drugs as defined by NQO1/catalase expression and highlights the role of alternate antioxidant components in dictating drug-induced oxidative stress.

Phenotype-specific estimation of metabolic fluxes using gene expression data

A cell's genome influences its metabolism via the expression of enzyme-related genes, but transcriptome and fluxome are not perfectly correlated as post-transcriptional mechanisms also regulate reaction's kinetics. Here, we addressed the question: given a transcriptome, how unobserved mechanisms of reaction kinetics should be systematically accounted for when inferring the fluxome? To infer the most likely and least biased fluxome, we present Pheflux, a constraint-based model maximizing Shannon's entropy of fluxes per mRNA. Benchmarked against ¹³C fluxes of yeast and bacteria, Pheflux accurately estimates the carbon core metabolism. We applied Pheflux to thousands of normal and tumor cell transcriptomes obtained from The Cancer Genome Atlas. Pheflux showed statistically significantly higher glucose yields on lactate in breast, kidney, and bronchus-lung tumoral cells than their normal counterparts. Results are consistent with the Warburg effect, a hallmark of cancer metabolism, suggesting that Pheflux can be efficiently used to study the metabolism of eukaryotic cells.

Systematic diet composition swap in a mouse genome-scale metabolic model reveals determinants of obesogenic diet metabolism in liver cancer

Dietary nutrient availability and gene expression, together, influence tissue metabolic activity. Here, we explore whether altering dietary nutrient composition in the context of mouse liver cancer suffices to overcome chronic gene expression changes that arise from tumorigenesis and western-style diet (WD). We construct a mouse genome-scale metabolic model and estimate metabolic fluxes in liver tumors and non-tumoral tissue after computationally varying the composition of input diet. This approach, called Systematic Diet Composition Swap (SyDiCoS), revealed that, compared to a control diet, WD increases production of glycerol and succinate irrespective of specific tissue gene expression patterns. Conversely, differences in fatty acid utilization pathways between tumor and non-tumor liver are amplified with WD by both dietary carbohydrates and lipids together. Our data suggest that combined dietary component modifications may be required to normalize the distinctive metabolic patterns that underlie selective targeting of tumor metabolism.

Single-cell and spatial analyses reveal the association between gene expression of glutamine synthetase with the immunosuppressive phenotype of APOE+CTSZ+TAM in cancers

An immunosuppressive state is regulated by various factors in the tumor microenvironment (TME), including, but not limited to, metabolic plasticity of immunosuppressive cells and cytokines secreted by these cells. We used single-cell RNA-sequencing (scRNA-seq) data and applied single-cell flux estimation analysis to characterize the link between metabolism and cellular function within the hypoxic TME of colorectal (CRC) and lung cancer. In terms of metabolic heterogeneity, we found myeloid cells potentially inclined to accumulate glutamine but tumor cells inclined to accumulate glutamate. In particular, we uncovered a tumor-associated macrophage (TAM) subpopulation, APOE+CTSZ+TAM, that was present in high proportions in tumor samples and exhibited immunosuppressive characteristics through upregulating the expression of anti-inflammatory genes. The proportion of APOE+CTSZ+TAM and regulatory T cells (Treg) were positively correlated across CRC scRNA-seq samples. APOE+CTSZ+TAM potentially interacted with Treg via CXCL16-CCR6 signals, as seen by ligand-receptor interactions analysis. Notably, glutamate-to-glutamine metabolic flux score and glutamine synthetase (GLUL) expression were uniquely higher in APOE+CTSZ+TAM, compared with other cell types within the TME. GLUL expression in macrophages was positively correlated with anti-inflammatory score and was higher in high-grade and invasive tumor samples. Moreover, spatial transcriptome and multiplex immunofluorescence staining of samples showed that APOE+CTSZ+TAM and Treg potentially colocalized in the tissue sections from CRC clinical samples. These results highlight the specific role and metabolic characteristic of the APOE+CTSZ+TAM subpopulation and provide a new perspective for macrophage subcluster-targeted therapeutic interventions or metabolic checkpoint-based cancer therapies.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

Multienzyme activity profiling for evaluation of cell‐to‐cell variability of metabolic state

In solid organs, cells of the same “type” can vary in their molecular phenotype. The basis of this state variation is being revealed by characterizing cell features including the expression pattern of mRNAs and the internal distribution of proteins. Here, the variability of metabolic state between cells is probed by enzyme activity profiling. We study individual cells of types that can be identified during the post‐mitotic phase of oogenesis in Xenopus laevis. Whole‐cell homogenates of isolated oocytes are used for kinetic analysis of enzymes, with a focus on the initial reaction rate. For each oocyte type studied, the activity signatures of glyceraldehyde 3‐phosphate dehydrogenase (GAPDH) and malate dehydrogenase 1 (MDH1) vary more between the homogenates of single oocytes than between repeat samplings of control homogenates. Unexpectedly, the activity signatures of GAPDH and MDH1 strongly co‐vary between oocytes of each type and change in strength of correlation during oogenesis. Therefore, variability of the kinetic behavior of these housekeeping enzymes between “identical” cells is physiologically programmed. Based on these findings, we propose that single‐cell profiling of enzyme kinetics will improve understanding of how metabolic state heterogeneity is related to heterogeneity revealed by omics methods including proteomics, epigenomics, and metabolomics.

Single-Cell Metabolomics in Hematopoiesis and Hematological Malignancies

Aberrant metabolism contributes to tumor initiation, progression, metastasis, and drug resistance. Metabolic dysregulation has emerged as a hallmark of several hematologic malignancies. Decoding the molecular mechanism underlying metabolic rewiring in hematological malignancies would provide promising avenues for novel therapeutic interventions. Single-cell metabolic analysis can directly offer a meaningful readout of the cellular phenotype, allowing us to comprehensively dissect cellular states and access biological information unobtainable from bulk analysis. In this review, we first highlight the unique metabolic properties of hematologic malignancies and underscore potential metabolic vulnerabilities. We then emphasize the emerging single-cell metabolomics techniques, aiming to provide a guide to interrogating metabolism at single-cell resolution. Furthermore, we summarize recent studies demonstrating the power of single-cell metabolomics to uncover the roles of metabolic rewiring in tumor biology, cellular heterogeneity, immunometabolism, and therapeutic resistance. Meanwhile, we describe a practical view of the potential applications of single-cell metabolomics in hematopoiesis and hematological malignancies. Finally, we present the challenges and perspectives of single-cell metabolomics development.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

Integrating knowledge and omics to decipher mechanisms via large-scale models of signaling networks

Signal transduction governs cellular behavior, and its dysregulation often leads to human disease. To understand this process, we can use network models based on prior knowledge, where nodes represent biomolecules, usually proteins, and edges indicate interactions between them. Several computational methods combine untargeted omics data with prior knowledge to estimate the state of signaling networks in specific biological scenarios. Here, we review, compare, and classify recent network approaches according to their characteristics in terms of input omics data, prior knowledge and underlying methodologies. We highlight existing challenges in the field, such as the general lack of ground truth and the limitations of prior knowledge. We also point out new omics developments that may have a profound impact, such as single‐cell proteomics or large‐scale profiling of protein conformational changes. We provide both an introduction for interested users seeking strategies to study cell signaling on a large scale and an update for seasoned modelers.

Angiogenesis goes computational – The future way forward to discover new angiogenic targets?

Multi-omics technologies are being increasingly utilized in angiogenesis research. Yet, computational methods have not been widely used for angiogenic target discovery and prioritization in this field, partly because (wet-lab) vascular biologists are insufficiently familiar with computational biology tools and the opportunities they may offer. With this review, written for vascular biologists who lack expertise in computational methods, we aspire to break boundaries between both fields and to illustrate the potential of these tools for future angiogenic target discovery. We provide a comprehensive survey of currently available computational approaches that may be useful in prioritizing candidate genes, predicting associated mechanisms, and identifying their specificity to endothelial cell subtypes. We specifically highlight tools that use flexible, machine learning frameworks for large-scale data integration and gene prioritization. For each purpose-oriented category of tools, we describe underlying conceptual principles, highlight interesting applications and discuss limitations. Finally, we will discuss challenges and recommend some guidelines which can help to optimize the process of accurate target discovery.