Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

Personalized metabolic whole-body models for newborns and infants predict growth and biomarkers of inherited metabolic diseases

Comprehensive whole-body models (WBMs) accounting for organ-specific dynamics have been developed to simulate adult metabolism, but such models do not exist for infants. Here, we present a resource of 360 organ-resolved, sex-specific models of newborn and infant metabolism (infant-WBMs) spanning the first 180 days of life. These infant-WBMs were parameterized to represent the distinct metabolic characteristics of newborns and infants, including nutrition, energy requirements, and thermoregulation. We demonstrate that the predicted infant growth was consistent with the recommendation by the World Health Organization. We assessed the infant-WBMs' reliability and capabilities for personalization by simulating 10,000 newborns based on their blood metabolome and birth weight. Furthermore, the infant-WBMs accurately predicted changes in known biomarkers over time and metabolic responses to treatment strategies for inherited metabolic diseases. The infant-WBM resource holds promise for personalized medicine, as the infant-WBMs could be a first step to digital metabolic twins for newborn and infant metabolism.

Cross-attention enables deep learning on limited omics-imaging-clinical data of 130 lung cancer patients

Deep-learning tools that extract prognostic factors derived from multi-omics data have recently contributed to individualized predictions of survival outcomes. However, the limited size of integrated omics-imaging-clinical datasets poses challenges. Here, we propose two biologically interpretable and robust deep-learning architectures for survival prediction of non-small cell lung cancer (NSCLC) patients, learning simultaneously from computed tomography (CT) scan images, gene expression data, and clinical information. The proposed models integrate patient-specific clinical, transcriptomic, and imaging data and incorporate Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathway information, adding biological knowledge within the learning process to extract prognostic gene biomarkers and molecular pathways. While both models accurately stratify patients in high- and low-risk groups when trained on a dataset of only 130 patients, introducing a cross-attention mechanism in a sparse autoencoder significantly improves the performance, highlighting tumor regions and NSCLC-related genes as potential biomarkers and thus offering a significant methodological advancement when learning from small imaging-omics-clinical samples.

Emerging methods for genome-scale metabolic modeling of microbial communities

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels.

Rational strain design with minimal phenotype perturbation

Devising genetic interventions for desired cellular phenotypes remains challenging regarding time and resources. Kinetic models can accelerate this task by simulating metabolic responses to genetic perturbations. However, exhaustive design evaluations with kinetic models are computationally impractical, especially when targeting multiple enzymes. Here, we introduce a framework for efficiently scouting the design space while respecting cellular physiological requirements. The framework employs mixed-integer linear programming and nonlinear simulations with large-scale nonlinear kinetic models to devise genetic interventions while accounting for the network effects of these perturbations. Importantly, it ensures the engineered strain's robustness by maintaining its phenotype close to that of the reference strain. The framework, applied to improve the anthranilate production in E. coli, devises designs for experimental implementation, including eight previously experimentally validated targets. We expect this framework to play a crucial role in future design-build-test-learn cycles, significantly expediting the strain design compared to exhaustive design enumeration.

Generalized metabolic flux analysis framework provides mechanism-based predictions of ophthalmic complications in type 2 diabetes patients

Chronic metabolic diseases arise from changes in metabolic fluxes through biomolecular pathways and gene networks accumulated over the lifetime of an individual. While clinical and biochemical profiles present just real-time snapshots of the patients' health, efficient computation models of the pathological disturbance of biomolecular processes are required to achieve individualized mechanistic insights into disease progression. Here, we describe the Generalized metabolic flux analysis (GMFA) for addressing this gap. Suitably grouping individual metabolites/fluxes into pools simplifies the analysis of the resulting more coarse-grain network. We also map non-metabolic clinical modalities onto the network with additional edges. Instead of using the time coordinate, the system status (metabolite concentrations and fluxes) is quantified as function of a generalized extent variable (a coordinate in the space of generalized metabolites) that represents the system's coordinate along its evolution path and evaluates the degree of change between any two states on that path. We applied GMFA to analyze Type 2 Diabetes Mellitus (T2DM) patients from two cohorts: EVAS (289 patients from Singapore) and NHANES (517) from the USA. Personalized systems biology models (digital twins) were constructed. We deduced disease dynamics from the individually parameterized metabolic network and predicted the evolution path of the metabolic health state. For each patient, we obtained an individual description of disease dynamics and predict an evolution path of the metabolic health state. Our predictive models achieve an ROC-AUC in the range 0.79-0.95 (sensitivity 80-92%, specificity 62-94%) in identifying phenotypes at the baseline and predicting future development of diabetic retinopathy and cataract progression among T2DM patients within 3 years from the baseline. The GMFA method is a step towards realizing the ultimate goal to develop practical predictive computational models for diagnostics based on systems biology. This tool has potential use in chronic disease management in medical practice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13755-023-00218-x.

Cell-specific genome-scale metabolic modeling of SARS-CoV-2-infected lung to identify antiviral enzymes

Computational systems biology plays a key role in the discovery of suitable antiviral targets. We designed a cell-specific, constraint-based modeling technique for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected lungs. We used the gene sequence of the alpha variant of SARS-CoV-2 to build a viral biomass reaction (VBR). We also used the mass proportion of lipids between the viral biomass and its host cell to estimate the stoichiometric coefficients of viral lipids in the reaction. We then integrated the VBR, the gene expression of the alpha variant of SARS-CoV-2, and the generic human metabolic network Recon3D to reconstruct a cell-specific genome-scale metabolic model. An antiviral target discovery (AVTD) platform was introduced using this model to identify therapeutic drug targets for combating COVID-19. The AVTD platform not only identified antiviral genes for eliminating viral replication but also predicted side effects of treatments. Our computational results revealed that knocking out dihydroorotate dehydrogenase (DHODH) might reduce the synthesis rate of cytidine-5'-triphosphate and uridine-5'-triphosphate, which terminate the viral building blocks of DNA and RNA for SARS-CoV-2 replication. Our results also indicated that DHODH is a promising antiviral target that causes minor side effects, which is consistent with the results of recent reports. Moreover, we discovered that the genes that participate in the de novo biosynthesis of glycerophospholipids and ceramides become unidentifiable if the VBR does not involve the stoichiometry of lipids.

Predicting weighted unobserved nodes in a regulatory network using answer set programming

BACKGROUND

The impact of a perturbation, over-expression, or repression of a key node on an organism, can be modelled based on a regulatory and/or metabolic network. Integration of these two networks could improve our global understanding of biological mechanisms triggered by a perturbation. This study focuses on improving the modelling of the regulatory network to facilitate a possible integration with the metabolic network. Previously proposed methods that study this problem fail to deal with a real-size regulatory network, computing predictions sensitive to perturbation and quantifying the predicted species behaviour more finely.

RESULTS

To address previously mentioned limitations, we develop a new method based on Answer Set Programming, MajS. It takes a regulatory network and a discrete partial set of observations as input. MajS tests the consistency between the input data, proposes minimal repairs on the network to establish consistency, and finally computes weighted and signed predictions over the network species. We tested MajS by comparing the HIF-1 signalling pathway with two gene-expression datasets. Our results show that MajS can predict 100% of unobserved species. When comparing MajS with two similar (discrete and quantitative) tools, we observed that compared with the discrete tool, MajS proposes a better coverage of the unobserved species, is more sensitive to system perturbations, and proposes predictions closer to real data. Compared to the quantitative tool, MajS provides more refined discrete predictions that agree with the dynamic proposed by the quantitative tool.

CONCLUSIONS

MajS is a new method to test the consistency between a regulatory network and a dataset that provides computational predictions on unobserved network species. It provides fine-grained discrete predictions by outputting the weight of the predicted sign as a piece of additional information. MajS' output, thanks to its weight, could easily be integrated with metabolic network modelling.

Applications of multi-omics analysis in human diseases

Multi-omics usually refers to the crossover application of multiple high-throughput screening technologies represented by genomics, transcriptomics, single-cell transcriptomics, proteomics and metabolomics, spatial transcriptomics, and so on, which play a great role in promoting the study of human diseases. Most of the current reviews focus on describing the development of multi-omics technologies, data integration, and application to a particular disease; however, few of them provide a comprehensive and systematic introduction of multi-omics. This review outlines the existing technical categories of multi-omics, cautions for experimental design, focuses on the integrated analysis methods of multi-omics, especially the approach of machine learning and deep learning in multi-omics data integration and the corresponding tools, and the application of multi-omics in medical researches (e.g., cancer, neurodegenerative diseases, aging, and drug target discovery) as well as the corresponding open-source analysis tools and databases, and finally, discusses the challenges and future directions of multi-omics integration and application in precision medicine. With the development of high-throughput technologies and data integration algorithms, as important directions of multi-omics for future disease research, single-cell multi-omics and spatial multi-omics also provided a detailed introduction. This review will provide important guidance for researchers, especially who are just entering into multi-omics medical research.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

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.

Integrated regulatory and metabolic networks of the tumor microenvironment for therapeutic target prioritization

Translation of genomic discovery, such as single-cell sequencing data, to clinical decisions remains a longstanding bottleneck in the field. Meanwhile, computational systems biological models, such as cellular metabolism models and cell signaling pathways, have emerged as powerful approaches to provide efficient predictions in metabolites and gene expression levels, respectively. However, there has been limited research on the integration between these two models. This work develops a methodology for integrating computational models of probabilistic gene regulatory networks with a constraint-based metabolism model. By using probabilistic reasoning with Bayesian Networks, we aim to predict cell-specific changes under different interventions, which are embedded into the constraint-based models of metabolism. Applications to single-cell sequencing data of glioblastoma brain tumors generate predictions about the effects of pharmaceutical interventions on the regulatory network and downstream metabolisms in different cell types from the tumor microenvironment. The model presents possible insights into treatments that could potentially suppress anaerobic metabolism in malignant cells with minimal impact on other cell types' metabolism. The proposed integrated model can guide therapeutic target prioritization, the formulation of combination therapies, and future drug discovery. This model integration framework is also generalizable to other applications, such as different cell types, organisms, and diseases.

Machine Learning Methods for Survival Analysis with Clinical and Transcriptomics Data of Breast Cancer

Breast cancer is one of the most common cancers in women worldwide, which causes an enormous number of deaths annually. However, early diagnosis of breast cancer can improve survival outcomes enabling simpler and more cost-effective treatments. The recent increase in data availability provides unprecedented opportunities to apply data-driven and machine learning methods to identify early-detection prognostic factors capable of predicting the expected survival and potential sensitivity to treatment of patients, with the final aim of enhancing clinical outcomes. This tutorial presents a protocol for applying machine learning models in survival analysis for both clinical and transcriptomic data. We show that integrating clinical and mRNA expression data is essential to explain the multiple biological processes driving cancer progression. Our results reveal that machine-learning-based models such as random survival forests, gradient boosted survival model, and survival support vector machine can outperform the traditional statistical methods, i.e., Cox proportional hazard model. The highest C-index among the machine learning models was recorded when using survival support vector machine, with a value 0.688, whereas the C-index recorded using the Cox model was 0.677. Shapley Additive Explanation (SHAP) values were also applied to identify the feature importance of the models and their impact on the prediction outcomes.

Vitamin D and Systems Biology

The biological actions of the vitamin D receptor (VDR) have been investigated intensively for over 100 years and has led to the identification of significant insights into the repertoire of its biological actions. These were initially established to be centered on the regulation of calcium transport in the colon and deposition in bone. Beyond these well-known calcemic roles, other roles have emerged in the regulation of cell differentiation processes and have an impact on metabolism. The purpose of the current review is to consider where applying systems biology (SB) approaches may begin to generate a more precise understanding of where the VDR is, and is not, biologically impactful. Two SB approaches have been developed and begun to reveal insight into VDR biological functions. In a top-down SB approach genome-wide scale data are statistically analyzed, and from which a role for the VDR emerges in terms of being a hub in a biological network. Such approaches have confirmed significant roles, for example, in myeloid differentiation and the control of inflammation and innate immunity. In a bottom-up SB approach, current biological understanding is built into a kinetic model which is then applied to existing biological data to explain the function and identify unknown behavior. To date, this has not been applied to the VDR, but has to the related ERα and identified previously unknown mechanisms of control. One arena where applying top-down and bottom-up SB approaches may be informative is in the setting of prostate cancer health disparities.

Methods for Stratification and Validation Cohorts: A Scoping Review

Personalized medicine requires large cohorts for patient stratification and validation of patient clustering. However, standards and harmonized practices on the methods and tools to be used for the design and management of cohorts in personalized medicine remain to be defined. This study aims to describe the current state-of-the-art in this area. A scoping review was conducted searching in PubMed, EMBASE, Web of Science, Psycinfo and Cochrane Library for reviews about tools and methods related to cohorts used in personalized medicine. The search focused on cancer, stroke and Alzheimer’s disease and was limited to reports in English, French, German, Italian and Spanish published from 2005 to April 2020. The screening process was reported through a PRISMA flowchart. Fifty reviews were included, mostly including information about how data were generated (25/50) and about tools used for data management and analysis (24/50). No direct information was found about the quality of data and the requirements to monitor associated clinical data. A scarcity of information and standards was found in specific areas such as sample size calculation. With this information, comprehensive guidelines could be developed in the future to improve the reproducibility and robustness in the design and management of cohorts in personalized medicine studies.

Toward modeling metabolic state from single-cell transcriptomics

BACKGROUND

Single-cell metabolic studies bring new insights into cellular function, which can often not be captured on other omics layers. Metabolic information has wide applicability, such as for the study of cellular heterogeneity or for the understanding of drug mechanisms and biomarker development. However, metabolic measurements on single-cell level are limited by insufficient scalability and sensitivity, as well as resource intensiveness, and are currently not possible in parallel with measuring transcript state, commonly used to identify cell types. Nevertheless, because omics layers are strongly intertwined, it is possible to make metabolic predictions based on measured data of more easily measurable omics layers together with prior metabolic network knowledge.

SCOPE OF REVIEW

We summarize the current state of single-cell metabolic measurement and modeling approaches, motivating the use of computational techniques. We review three main classes of computational methods used for prediction of single-cell metabolism: pathway-level analysis, constraint-based modeling, and kinetic modeling. We describe the unique challenges arising when transitioning from bulk to single-cell modeling. Finally, we propose potential model extensions and computational methods that could be leveraged to achieve these goals.

MAJOR CONCLUSIONS

Single-cell metabolic modeling is a rising field that provides a new perspective for understanding cellular functions. The presented modeling approaches vary in terms of input requirements and assumptions, scalability, modeled metabolic layers, and newly gained insights. We believe that the use of prior metabolic knowledge will lead to more robust predictions and will pave the way for mechanistic and interpretable machine-learning models.

The role of metabolic ecosystem in cancer progression — metabolic plasticity and mTOR hyperactivity in tumor tissues

Despite advancements in cancer management, tumor relapse and metastasis are associated with poor outcomes in many cancers. Over the past decade, oncogene-driven carcinogenesis, dysregulated cellular signaling networks, dynamic changes in the tissue microenvironment, epithelial-mesenchymal transitions, protein expression within regulatory pathways, and their part in tumor progression are described in several studies. However, the complexity of metabolic enzyme expression is considerably under evaluated. Alterations in cellular metabolism determine the individual phenotype and behavior of cells, which is a well-recognized hallmark of cancer progression, especially in the adaptation mechanisms underlying therapy resistance. In metabolic symbiosis, cells compete, communicate, and even feed each other, supervised by tumor cells. Metabolic reprogramming forms a unique fingerprint for each tumor tissue, depending on the cellular content and genetic, epigenetic, and microenvironmental alterations of the developing cancer. Based on its sensing and effector functions, the mechanistic target of rapamycin (mTOR) kinase is considered the master regulator of metabolic adaptation. Moreover, mTOR kinase hyperactivity is associated with poor prognosis in various tumor types. In situ metabolic phenotyping in recent studies highlights the importance of metabolic plasticity, mTOR hyperactivity, and their role in tumor progression. In this review, we update recent developments in metabolic phenotyping of the cancer ecosystem, metabolic symbiosis, and plasticity which could provide new research directions in tumor biology. In addition, we suggest pathomorphological and analytical studies relating to metabolic alterations, mTOR activity, and their associations which are necessary to improve understanding of tumor heterogeneity and expand the therapeutic management of cancer.

A Practical Guide to Integrating Multimodal Machine Learning and Metabolic Modeling

Complex, distributed, and dynamic sets of clinical biomedical data are collectively referred to as multimodal clinical data. In order to accommodate the volume and heterogeneity of such diverse data types and aid in their interpretation when they are combined with a multi-scale predictive model, machine learning is a useful tool that can be wielded to deconstruct biological complexity and extract relevant outputs. Additionally, genome-scale metabolic models (GSMMs) are one of the main frameworks striving to bridge the gap between genotype and phenotype by incorporating prior biological knowledge into mechanistic models. Consequently, the utilization of GSMMs as a foundation for the integration of multi-omic data originating from different domains is a valuable pursuit towards refining predictions. In this chapter, we show how cancer multi-omic data can be analyzed via multimodal machine learning and metabolic modeling. Firstly, we focus on the merits of adopting an integrative systems biology led approach to biomedical data mining. Following this, we propose how constraint-based metabolic models can provide a stable yet adaptable foundation for the integration of multimodal data with machine learning. Finally, we provide a step-by-step tutorial for the combination of machine learning and GSMMs, which includes: (i) tissue-specific constraint-based modeling; (ii) survival analysis using time-to-event prediction for cancer; and (iii) classification and regression approaches for multimodal machine learning. The code associated with the tutorial can be found at https://github.com/Angione-Lab/Tutorials_Combining_ML_and_GSMM .

Reconstruction of tissue-specific genome-scale metabolic models for human cancer stem cells

Cancer Stem Cells (CSCs) contribute to cancer aggressiveness, metastasis, chemo/radio-therapy resistance, and tumor recurrence. Recent studies emphasized the importance of metabolic reprogramming of CSCs for the maintenance and progression of the cancer phenotype through both the fulfillment of the energetic requirements and the supply of substrates fundamental for fast-cell growth, as well as through metabolite-induced epigenetic regulation. Therefore, it is of paramount importance to develop therapeutic strategies tailored to target the metabolism of CSCs. In this work, we built computational Genome-Scale Metabolic Models (GSMMs) for CSCs of different tissues. Flux simulations were then used to predict metabolic phenotypes, identify potential therapeutic targets, and spot already-known Transcription Factors (TFs), miRNAs and antimetabolites that could be used as part of drug repurposing strategies against cancer. Results were in accordance with experimental evidence, provided insights of new metabolic mechanisms for already known agents, and allowed for the identification of potential new targets and compounds that could be interesting for further in vitro and in vivo validation.

Protocol for hybrid flux balance, statistical, and machine learning analysis of multi-omic data from the cyanobacterium Synechococcus sp. PCC 7002

Combining a computational framework for flux balance analysis with machine learning improves the accuracy of predicting metabolic activity across conditions, while enabling mechanistic interpretation. This protocol presents a guide to condition-specific metabolic modeling that integrates regularized flux balance analysis with machine learning approaches to extract key features from transcriptomic and fluxomic data. We demonstrate the protocol as applied to Synechococcus sp. PCC 7002; we also outline how it can be adapted to any species or community with available multi-omic data. For complete details on the use and execution of this protocol, please refer to Vijayakumar et al. (2020).

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.

Machine learning meets omics: applications and perspectives

The innovation of biotechnologies has allowed the accumulation of omics data at an alarming rate, thus introducing the era of 'big data'. Extracting inherent valuable knowledge from various omics data remains a daunting problem in bioinformatics. Better solutions often need some kind of more innovative methods for efficient handlings and effective results. Recent advancements in integrated analysis and computational modeling of multi-omics data helped address such needs in an increasingly harmonious manner. The development and application of machine learning have largely advanced our insights into biology and biomedicine and greatly promoted the development of therapeutic strategies, especially for precision medicine. Here, we propose a comprehensive survey and discussion on what happened, is happening and will happen when machine learning meets omics. Specifically, we describe how artificial intelligence can be applied to omics studies and review recent advancements at the interface between machine learning and the ever-widest range of omics including genomics, transcriptomics, proteomics, metabolomics, radiomics, as well as those at the single-cell resolution. We also discuss and provide a synthesis of ideas, new insights, current challenges and perspectives of machine learning in omics.

Discovering Essential Multiple Gene Effects Through Large Scale Optimization: An Application to Human Cancer Metabolism

Computational modelling of metabolic processes has proven to be a useful approach to formulate our knowledge and improve our understanding of core biochemical systems that are crucial to maintaining cellular functions. Towards understanding the broader role of metabolism on cellular decision-making in health and disease conditions, it is important to integrate the study of metabolism with other core regulatory systems and omics within the cell, including gene expression patterns. After quantitatively integrating gene expression profiles with a genome-scale reconstruction of human metabolism, we propose a set of combinatorial methods to reverse engineer gene expression profiles and to find pairs and higher-order combinations of genetic modifications that simultaneously optimize multi-objective cellular goals. This enables us to suggest classes of transcriptomic profiles that are most suitable to achieve given metabolic phenotypes. We demonstrate how our techniques are able to compute beneficial, neutral or "toxic" combinations of gene expression levels. We test our methods on nine tissue-specific cancer models, comparing our outcomes with the corresponding normal cells, identifying genes as targets for potential therapies. Our methods open the way to a broad class of applications that require an understanding of the interplay among genotype, metabolism, and cellular behaviour, at scale.

Genome-scale metabolic modelling of SARS-CoV-2 in cancer cells reveals an increased shift to glycolytic energy production

Cancer is considered a high‐risk condition for severe illness resulting from COVID‐19. The interaction between severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) and human metabolism is key to elucidating the risk posed by COVID‐19 for cancer patients and identifying effective treatments, yet it is largely uncharacterised on a mechanistic level. We present a genome‐scale map of short‐term metabolic alterations triggered by SARS‐CoV‐2 infection of cancer cells. Through transcriptomic‐ and proteomic‐informed genome‐scale metabolic modelling, we characterise the role of RNA and fatty acid biosynthesis in conjunction with a rewiring in energy production pathways and enhanced cytokine secretion. These findings link together complementary aspects of viral invasion of cancer cells, while providing mechanistic insights that can inform the development of treatment strategies.

A Convergent Functional Genomics Analysis to Identify Biological Regulators Mediating Effects of Creatine Supplementation

Creatine (Cr) and phosphocreatine (PCr) are physiologically essential molecules for life, given they serve as rapid and localized support of energy- and mechanical-dependent processes. This evolutionary advantage is based on the action of creatine kinase (CK) isozymes that connect places of ATP synthesis with sites of ATP consumption (the CK/PCr system). Supplementation with creatine monohydrate (CrM) can enhance this system, resulting in well-known ergogenic effects and potential health or therapeutic benefits. In spite of our vast knowledge about these molecules, no integrative analysis of molecular mechanisms under a systems biology approach has been performed to date; thus, we aimed to perform for the first time a convergent functional genomics analysis to identify biological regulators mediating the effects of Cr supplementation in health and disease. A total of 35 differentially expressed genes were analyzed. We identified top-ranked pathways and biological processes mediating the effects of Cr supplementation. The impact of CrM on miRNAs merits more research. We also cautiously suggest two dose-response functional pathways (kinase- and ubiquitin-driven) for the regulation of the Cr uptake. Our functional enrichment analysis, the knowledge-based pathway reconstruction, and the identification of hub nodes provide meaningful information for future studies. This work contributes to a better understanding of the well-reported benefits of Cr in sports and its potential in health and disease conditions, although further clinical research is needed to validate the proposed mechanisms.

Undisclosed, unmet and neglected challenges in multi-omics studies

Multi-omics approaches have become a reality in both large genomics projects and small laboratories. However, the multi-omics research community still faces a number of issues that have either not been sufficiently discussed or for which current solutions are still limited. In this Perspective, we elaborate on these limitations and suggest points of attention for future research. We finally discuss new opportunities and challenges brought to the field by the rapid development of single-cell high-throughput molecular technologies.

Abnormal brain cholesterol homeostasis in Alzheimer's disease-a targeted metabolomic and transcriptomic study

The role of brain cholesterol metabolism in Alzheimer's disease (AD) remains unclear. Peripheral and brain cholesterol levels are largely independent due to the impermeability of the blood brain barrier (BBB), highlighting the importance of studying the role of brain cholesterol homeostasis in AD. We first tested whether metabolite markers of brain cholesterol biosynthesis and catabolism were altered in AD and associated with AD pathology using linear mixed-effects models in two brain autopsy samples from the Baltimore Longitudinal Study of Aging (BLSA) and the Religious Orders Study (ROS). We next tested whether genetic regulators of brain cholesterol biosynthesis and catabolism were altered in AD using the ANOVA test in publicly available brain tissue transcriptomic datasets. Finally, using regional brain transcriptomic data, we performed genome-scale metabolic network modeling to assess alterations in cholesterol biosynthesis and catabolism reactions in AD. We show that AD is associated with pervasive abnormalities in cholesterol biosynthesis and catabolism. Using transcriptomic data from Parkinson's disease (PD) brain tissue samples, we found that gene expression alterations identified in AD were not observed in PD, suggesting that these changes may be specific to AD. Our results suggest that reduced de novo cholesterol biosynthesis may occur in response to impaired enzymatic cholesterol catabolism and efflux to maintain brain cholesterol levels in AD. This is accompanied by the accumulation of nonenzymatically generated cytotoxic oxysterols. Our results set the stage for experimental studies to address whether abnormalities in cholesterol metabolism are plausible therapeutic targets in AD.

Exploring the Metabolic Heterogeneity of Cancers: A Benchmark Study of Context-Specific Models

Metabolic heterogeneity is a hallmark of cancer and can distinguish a normal phenotype from a cancer phenotype. In the systems biology domain, context-specific models facilitate extracting physiologically relevant information from high-quality data. Here, to utilize the heterogeneity of metabolic patterns to discover biomarkers of all cancers, we benchmarked thousands of context-specific models using well-established algorithms for the integration of omics data into the generic human metabolic model Recon3D. By analyzing the active reactions capable of carrying flux and their magnitude through flux balance analysis, we proved that the metabolic pattern of each cancer is unique and could act as a cancer metabolic fingerprint. Subsequently, we searched for proper feature selection methods to cluster the flux states characterizing each cancer. We employed PCA-based dimensionality reduction and a random forest learning algorithm to reveal reactions containing the most relevant information in order to effectively identify the most influential fluxes. Conclusively, we discovered different pathways that are probably the main sources for metabolic heterogeneity in cancers. We designed the GEMbench website to interactively present the data, methods, and analysis results.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

Mapping the Metabolic Networks of Tumor Cells and Cancer-Associated Fibroblasts

Metabolism is considered to be the core of all cellular activity. Thus, extensive studies of metabolic processes are ongoing in various fields of biology, including cancer research. Cancer cells are known to adapt their metabolism to sustain high proliferation rates and survive in unfavorable environments with low oxygen and nutrient concentrations. Hence, targeting cancer cell metabolism is a promising therapeutic strategy in cancer research. However, cancers consist not only of genetically altered tumor cells but are interwoven with endothelial cells, immune cells and fibroblasts, which together with the extracellular matrix (ECM) constitute the tumor microenvironment (TME). Cancer-associated fibroblasts (CAFs), which are linked to poor prognosis in different cancer types, are one important component of the TME. CAFs play a significant role in reprogramming the metabolic landscape of tumor cells, but how, and in what manner, this interaction takes place remains rather unclear. This review aims to highlight the metabolic landscape of tumor cells and CAFs, including their recently identified subtypes, in different tumor types. In addition, we discuss various in vitro and in vivo metabolic techniques as well as different in silico computational tools that can be used to identify and characterize CAF-tumor cell interactions. Finally, we provide our view on how mapping the complex metabolic networks of stromal-tumor metabolism will help in finding novel metabolic targets for cancer treatment.

Multiomics analysis of the mdx/mTR mouse model of Duchenne muscular dystrophy

PURPOSE/AIM

Duchenne muscular dystrophy (DMD) is a progressive neuromuscular disease characterized by extensive muscle weakness. Patients with DMD lack a functional dystrophin protein, which transmits force and organizes the cytoskeleton of skeletal muscle. Multiomic studies have been proposed as a way to obtain novel insight about disease processes from preclinical models, and we used this approach to study pathological changes in dystrophic muscles.

MATERIALS AND METHODS

We evaluated hindlimb muscles of male mdx/mTR mice, which lack a functional dystrophin protein and have deficits in satellite cell abundance and proliferative capacity. Wild type (WT) C57BL/6 J mice served as controls. Muscle fiber contractility was measured, along with changes in the transcriptome using RNA sequencing, and in the proteome, metabolome, and lipidome using mass spectrometry.

RESULTS

While mdx/mTR mice displayed gross pathological changes and continued cycles of degeneration and regeneration, we found no differences in permeabilized fiber contractility between strains. However, there were numerous changes in the transcriptome and proteome related to protein balance, contractile elements, extracellular matrix, and metabolism. There was only a 53% agreement in fold-change data between the proteome and transcriptome. Numerous changes in markers of skeletal muscle metabolism were observed, with dystrophic muscles exhibiting elevated glycolytic metabolites such as 6-phosphoglycerate, fructose-6-phosphate and glucose-6-phosphate, fructose bisphosphate, phosphorylated hexoses, and phosphoenolpyruvate.

CONCLUSIONS

These findings highlight the utility of multiomics in studying muscle disease, and provide additional insight into the pathological changes in dystrophic muscles that might help to indirectly guide evidence-based nutritional or exercise prescription in DMD patients.

Overtraining Syndrome as a Complex Systems Phenomenon

The phenomenon of reduced athletic performance following sustained, intense training (Overtraining Syndrome, and OTS) was first recognized more than 90 years ago. Although hundreds of scientific publications have focused on OTS, a definitive diagnosis, reliable biomarkers, and effective treatments remain unknown. The present review considers existing models of OTS, acknowledges the individualized and sport-specific nature of signs/symptoms, describes potential interacting predisposing factors, and proposes that OTS will be most effectively characterized and evaluated via the underlying complex biological systems. Complex systems in nature are not aptly characterized or successfully analyzed using the classic scientific method (i.e., simplifying complex problems into single variables in a search for cause-and-effect) because they result from myriad (often non-linear) concomitant interactions of multiple determinants. Thus, this review 1) proposes that OTS be viewed from the perspectives of complex systems and network physiology, 2) advocates for and recommends that techniques such as trans-omic analyses and machine learning be widely employed, and 3) proposes evidence-based areas for future OTS investigations, including concomitant multi-domain analyses incorporating brain neural networks, dysfunction of hypothalamic-pituitary-adrenal responses to training stress, the intestinal microbiota, immune factors, and low energy availability. Such an inclusive and modern approach will measurably help in prevention and management of OTS.

A Hybrid Flux Balance Analysis and Machine Learning Pipeline Elucidates Metabolic Adaptation in Cyanobacteria

Machine learning has recently emerged as a promising tool for inferring multi-omic relationships in biological systems. At the same time, genome-scale metabolic models (GSMMs) can be integrated with such multi-omic data to refine phenotypic predictions. In this work, we use a multi-omic machine learning pipeline to analyze a GSMM of Synechococcus sp. PCC 7002, a cyanobacterium with large potential to produce renewable biofuels. We use regularized flux balance analysis to observe flux response between conditions across photosynthesis and energy metabolism. We then incorporate principal-component analysis, k-means clustering, and LASSO regularization to reduce dimensionality and extract key cross-omic features. Our results suggest that combining metabolic modeling with machine learning elucidates mechanisms used by cyanobacteria to cope with fluctuations in light intensity and salinity that cannot be detected using transcriptomics alone. Furthermore, GSMMs introduce critical mechanistic details that improve the performance of omic-based machine learning methods.

A Quantitative Proteome Map of the Human Body

Determining protein levels in each tissue and how they compare with RNA levels is important for understanding human biology and disease as well as regulatory processes that control protein levels. We quantified the relative protein levels from over 12,000 genes across 32 normal human tissues. Tissue-specific or tissue-enriched proteins were identified and compared to transcriptome data. Many ubiquitous transcripts are found to encode tissue-specific proteins. Discordance of RNA and protein enrichment revealed potential sites of synthesis and action of secreted proteins. The tissue-specific distribution of proteins also provides an in-depth view of complex biological events that require the interplay of multiple tissues. Most importantly, our study demonstrated that protein tissue-enrichment information can explain phenotypes of genetic diseases, which cannot be obtained by transcript information alone. Overall, our results demonstrate how understanding protein levels can provide insights into regulation, secretome, metabolism, and human diseases.

Modeling tissue-relevant Caenorhabditis elegans metabolism at network, pathway, reaction, and metabolite levels

Metabolism is a highly compartmentalized process that provides building blocks for biomass generation during development, homeostasis, and wound healing, and energy to support cellular and organismal processes. In metazoans, different cells and tissues specialize in different aspects of metabolism. However, studying the compartmentalization of metabolism in different cell types in a whole animal and for a particular stage of life is difficult. Here, we present MEtabolic models Reconciled with Gene Expression (MERGE), a computational pipeline that we used to predict tissue-relevant metabolic function at the network, pathway, reaction, and metabolite levels based on single-cell RNA-sequencing (scRNA-seq) data from the nematode Caenorhabditis elegans. Our analysis recapitulated known tissue functions in C. elegans, captured metabolic properties that are shared with similar tissues in human, and provided predictions for novel metabolic functions. MERGE is versatile and applicable to other systems. We envision this work as a starting point for the development of metabolic network models for individual cells as scRNA-seq continues to provide higher-resolution gene expression data.

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.

Deep learning meets metabolomics: a methodological perspective

Deep learning (DL), an emerging area of investigation in the fields of machine learning and artificial intelligence, has markedly advanced over the past years. DL techniques are being applied to assist medical professionals and researchers in improving clinical diagnosis, disease prediction and drug discovery. It is expected that DL will help to provide actionable knowledge from a variety of 'big data', including metabolomics data. In this review, we discuss the applicability of DL to metabolomics, while presenting and discussing several examples from recent research. We emphasize the use of DL in tackling bottlenecks in metabolomics data acquisition, processing, metabolite identification, as well as in metabolic phenotyping and biomarker discovery. Finally, we discuss how DL is used in genome-scale metabolic modelling and in interpretation of metabolomics data. The DL-based approaches discussed here may assist computational biologists with the integration, prediction and drawing of statistical inference about biological outcomes, based on metabolomics data.

Metabolomics and Multi-Omics Integration: A Survey of Computational Methods and Resources

As researchers are increasingly able to collect data on a large scale from multiple clinical and omics modalities, multi-omics integration is becoming a critical component of metabolomics research. This introduces a need for increased understanding by the metabolomics researcher of computational and statistical analysis methods relevant to multi-omics studies. In this review, we discuss common types of analyses performed in multi-omics studies and the computational and statistical methods that can be used for each type of analysis. We pinpoint the caveats and considerations for analysis methods, including required parameters, sample size and data distribution requirements, sources of a priori knowledge, and techniques for the evaluation of model accuracy. Finally, for the types of analyses discussed, we provide examples of the applications of corresponding methods to clinical and basic research. We intend that our review may be used as a guide for metabolomics researchers to choose effective techniques for multi-omics analyses relevant to their field of study.

Ameliorative activity of Adansonia digitata fruit on high sugar/high fat diet-simulated Metabolic Syndrome model in male Wistar rats

Metabolic syndrome is a complex of metabolic disorders characterized by oxidative stress which compromises cell functions and entails multiple organs pathologies. We investigated the therapeutic and protective potential of Adansonia digitata fruit -a potent antioxidant- in high sugar/high fat diet-simulated metabolic syndrome in Wistar rats. 42 male rats (140-200 g) were randomly divided into 7 groups. G1 was kept on standard laboratory diet (SLD) for all 9 weeks (negative control). 5 groups were fed high Sugar/high fat diet for 6 weeks then switched to SLD for another 3 weeks + oral treatment as follows: G2+ no treatment (positive control), G3-G5 + 200, 400 and 800 mg/kg/day aqueous A. digitata fruit respectively, G6 + 10 mg/kg/day Simvastatin. G7 + HS/HFD + 400 mg/kg/day A. digitata fruit simultaneously and was terminated at W6. Our results showed that G2-G6 develops dyslipidemia, hyperglycaemia, weight gain, elevated hepatic biomarkers, elevated creatinine and urea plus pathological derangements in the heart, liver and kidney tissues compared to negative control at W6. 200 mg/kg/day A. digitata fruit significantly ameliorated the induced dyslipidemia (P ≤ 0.001), hyperglycaemia (P ≤ 0.001) with a significant reduction in the Atherogenic Index of Plasma (P ≤ 0.000) after 3 weeks treatment. The fruit normalized the elevated hepatic biomarkers as well as creatinine and urea. A dose dependent partial reduction in lesion intensity was observed in the hepatic tissue while the heart and kidney showed mostly reversed to normal histology. The inflammatory infiltration was eliminated. Relevant results were observed for the two higher doses. The simultaneous treatment showed significant lower levels in all biomarkers investigated compared to positive control which could be interpreted as protective activity. A reduction of 4-11% in whole body weight was achieved. CONCLUSION: MetS was successfully simulated with a HS/HFD formula in male Wistar rats. Treatment with aqueous A. digitata fruit showed anti-Metabolic Syndrome potential reflected by weight loss, anti-inflammatory, hypolipidemic, hypoglycaemic, renal, hepatic and cardio-protective activities.

In silico disease model: from simple networks to complex diseases

In this chapter we consider in silico modeling of diseases starting from some simple to some complex (and mathematical) concepts. Examples and applications of in silico modeling for some important categories of diseases (such as for cancers, infectious diseases, and neuronal diseases) are also given.