Integration of single-cell RNA-seq data into metabolic models to characterize tumour cell populations

Motivation Metabolic reprogramming is a general feature of cancer cells. Regrettably, the comprehensive quantification of metabolites in biological specimens does not promptly translate into knowledge on the utilization of metabolic pathways. Computational models hold the promise to bridge this gap, by estimating fluxes across metabolic pathways. Yet they currently portray the average behavior of intermixed subpopulations, masking their inherent heterogeneity known to hinder cancer diagnosis and treatment. If complemented with the information on single-cell transcriptome, now enabled by RNA sequencing (scRNA-seq), metabolic models of cancer populations are expected to empower the characterization of the mechanisms behind metabolic heterogeneity. To this aim, we propose single-cell Flux Balance Analysis (scFBA) as a computational framework to translate sc-transcriptomes into single-cell fluxomes. Results We show that the integration of scRNA-seq profiles of cells derived from lung ade-nocarcinoma and breast cancer patients, into a multi-scale stoichiometric model of cancer population: 1) significantly reduces the space of feasible single-cell fluxomes; 2) allows to identify clusters of cells with different growth rates within the population; 3) points out the possible metabolic interactions among cells via exchange of metabolites. Availability The scFBA suite of MATLAB functions is available at https://github.com/BIMIB-DISCo/scFBA, as well as the case study datasets. Contact chiara.damiani@unimib.it

[1]  M. V. Heiden,et al.  Targeting cancer metabolism: a therapeutic window opens , 2011, Nature Reviews Drug Discovery.

[2]  T. Sauter,et al.  Constraint Based Modeling Going Multicellular , 2016, Front. Mol. Biosci..

[3]  J. Locasale,et al.  Metabolomics: A Primer. , 2017, Trends in biochemical sciences.

[4]  Anthony T Papenfuss,et al.  How the evolution of multicellularity set the stage for cancer , 2018, British Journal of Cancer.

[5]  Giancarlo Mauri,et al.  popFBA: tackling intratumour heterogeneity with Flux Balance Analysis , 2017, Bioinform..

[6]  S. Croubels,et al.  Metabolic fingerprinting reveals a novel candidate biomarker for prednisolone treatment in cattle , 2015, Metabolomics.

[7]  U. Sauer,et al.  Frontiers of high-throughput metabolomics. , 2017, Current opinion in chemical biology.

[8]  Neil Swainston,et al.  Improving metabolic flux predictions using absolute gene expression data , 2012, BMC Systems Biology.

[9]  D. Sabatini,et al.  Cancer cell metabolism: one hallmark, many faces. , 2012, Cancer discovery.

[10]  R Zenobi,et al.  Single-Cell Metabolomics: Analytical and Biological Perspectives , 2013, Science.

[11]  Daniel Machado,et al.  Systematic Evaluation of Methods for Integration of Transcriptomic Data into Constraint-Based Models of Metabolism , 2014, PLoS Comput. Biol..

[12]  I. Tannock,et al.  Drug resistance and the solid tumor microenvironment. , 2007, Journal of the National Cancer Institute.

[13]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[14]  P. Ward,et al.  Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. , 2012, Cancer cell.

[15]  Jeong Eon Lee,et al.  Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer , 2017, Nature Communications.

[16]  Giancarlo Mauri,et al.  An ensemble evolutionary constraint-based approach to understand the emergence of metabolic phenotypes , 2014, Natural Computing.

[17]  Giancarlo Mauri,et al.  Constraint-Based Modeling and Simulation of Cell Populations , 2016, WIVACE.

[18]  S. Teichmann,et al.  Computational and analytical challenges in single-cell transcriptomics , 2015, Nature Reviews Genetics.

[19]  K. Bachman,et al.  Analysis of glutamine dependency in non-small cell lung cancer , 2012, Cancer biology & therapy.

[20]  Eytan Ruppin,et al.  Modeling cancer metabolism on a genome scale , 2015, Molecular systems biology.

[21]  S. Beloribi-Djefaflia,et al.  Lipid metabolic reprogramming in cancer cells , 2016, Oncogenesis.

[22]  Giancarlo Mauri,et al.  Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models , 2016, Comput. Biol. Chem..

[23]  J. Swinnen,et al.  Cancer Cells Differentially Activate and Thrive on De Novo Lipid Synthesis Pathways in a Low-Lipid Environment , 2014, PloS one.

[24]  Giancarlo Mauri,et al.  A metabolic core model elucidates how enhanced utilization of glucose and glutamine, with enhanced glutamine-dependent lactate production, promotes cancer cell growth: The WarburQ effect , 2017, PLoS Comput. Biol..

[25]  R. Aebersold,et al.  On the Dependency of Cellular Protein Levels on mRNA Abundance , 2016, Cell.

[26]  B. Palsson,et al.  Large-scale in silico modeling of metabolic interactions between cell types in the human brain , 2010, Nature Biotechnology.

[27]  U. Sauer,et al.  Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli , 2007, Molecular systems biology.

[28]  Matteo Mori,et al.  Counting and Correcting Thermodynamically Infeasible Flux Cycles in Genome-Scale Metabolic Networks , 2013, Metabolites.

[29]  P. Fortina,et al.  The reverse Warburg effect: Aerobic glycolysis in cancer associated fibroblasts and the tumor stroma , 2009, Cell cycle.

[30]  M. Junttila,et al.  Influence of tumour micro-environment heterogeneity on therapeutic response , 2013, Nature.

[31]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[32]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

[33]  D. Hanahan,et al.  Hallmarks of Cancer: The Next Generation , 2011, Cell.

[34]  Anne Richelle,et al.  A Systematic Evaluation of Methods for Tailoring Genome-Scale Metabolic Models. , 2017, Cell systems.

[35]  Gurmit Singh,et al.  Inhibitors of glutamate release from breast cancer cells; new targets for cancer-induced bone-pain , 2015, Scientific Reports.

[36]  Travers Ching,et al.  Single-Cell Transcriptomics Bioinformatics and Computational Challenges , 2016, Front. Genet..

[37]  Do-Hyun Nam,et al.  Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells , 2015, Genome Biology.

[38]  Neil Swainston,et al.  Recon 2.2: from reconstruction to model of human metabolism , 2016, Metabolomics.

[39]  F. Bruggeman,et al.  Community Flux Balance Analysis for Microbial Consortia at Balanced Growth , 2013, PloS one.

[40]  Gabor Tigyi,et al.  Lysophosphatidic Acid Inhibits CD8 T-cell Activation and Control of Tumor Progression , 2013, Cancer Immunology Research.

[41]  S. Kauffman,et al.  Cell-cell interaction and diversity of emergent behaviours. , 2011, IET systems biology.

[42]  G. Mauri,et al.  MaREA: Metabolic feature extraction, enrichment and visualization of RNAseq data , 2018, bioRxiv.

[43]  D. Hanahan,et al.  The Hallmarks of Cancer , 2000, Cell.

[44]  Marissa Fessenden,et al.  Metabolomics: Small molecules, single cells , 2016, Nature.

[45]  C. Blanpain,et al.  Unravelling cancer stem cell potential , 2013, Nature Reviews Cancer.