{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,5,28]],"date-time":"2026-05-28T01:16:16Z","timestamp":1779930976452,"version":"3.53.1"},"reference-count":49,"publisher":"Oxford University Press (OUP)","issue":"5","license":[{"start":{"date-parts":[[2023,4,11]],"date-time":"2023-04-11T00:00:00Z","timestamp":1681171200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"NSF Research Traineeship Program"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2023,5,4]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>The accurate prediction of complex phenotypes such as metabolic fluxes in living systems is a grand challenge for systems biology and central to efficiently identifying biotechnological interventions that can address pressing industrial needs. The application of gene expression data to improve the accuracy of metabolic flux predictions using mechanistic modeling methods such as flux balance analysis (FBA) has not been previously demonstrated in multi-tissue systems, despite their biotechnological importance. We hypothesized that a method for generating metabolic flux predictions informed by relative expression levels between tissues would improve prediction accuracy.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>Relative gene expression levels derived from multiple transcriptomic and proteomic datasets were integrated into FBA predictions of a multi-tissue, diel model of Arabidopsis thaliana\u2019s central metabolism. This integration dramatically improved the agreement of flux predictions with experimentally based flux maps from 13C metabolic flux analysis compared with a standard parsimonious FBA approach. Disagreement between FBA predictions and MFA flux maps was measured using weighted averaged percent error values, and for parsimonious FBA this was169%\u2013180% for high light conditions and 94%\u2013103% for low light conditions, depending on the gene expression dataset used. This fell to 10%-13% and 9%-11% upon incorporating expression data into the modeling process, which also substantially altered the predicted carbon and energy economy of the plant.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>Code and data generated as part of this study are available from https:\/\/github.com\/Gibberella\/ArabidopsisGeneExpressionWeights.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btad186","type":"journal-article","created":{"date-parts":[[2023,4,11]],"date-time":"2023-04-11T11:18:53Z","timestamp":1681211933000},"source":"Crossref","is-referenced-by-count":11,"title":["Accurate flux predictions using tissue-specific gene expression in plant metabolic modeling"],"prefix":"10.1093","volume":"39","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-1942-0315","authenticated-orcid":false,"given":"Joshua A M","family":"Kaste","sequence":"first","affiliation":[{"name":"Department of Biochemistry and Molecular Biology, Michigan State University , 603 Wilson Rd. , East Lansing, Michigan 48824, United States"},{"name":"Department of Plant Biology, Michigan State University , 612 Wilson Rd. , East Lansing, Michigan 48824, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-8793-5084","authenticated-orcid":false,"given":"Yair","family":"Shachar-Hill","sequence":"additional","affiliation":[{"name":"Department of Plant Biology, Michigan State University , 612 Wilson Rd. , East Lansing, Michigan 48824, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2023,4,11]]},"reference":[{"key":"2023050421464706500_btad186-B1","doi-asserted-by":"crossref","first-page":"317","DOI":"10.1007\/s10295-015-1585-x","article-title":"Methods and advances in metabolic flux analysis: a mini-review","volume":"42","author":"Antoniewicz","year":"2015","journal-title":"J Ind Microbiol Biotechnol"},{"key":"2023050421464706500_btad186-B2","doi-asserted-by":"crossref","first-page":"1380","DOI":"10.1104\/pp.114.235358","article-title":"Bottom-up metabolic reconstruction of Arabidopsis and its application to determining the metabolic costs of enzyme production","volume":"165","author":"Arnold","year":"2014","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B3","doi-asserted-by":"crossref","first-page":"e0177292","DOI":"10.1371\/journal.pone.0177292","article-title":"Metabolic flux analysis of heterotrophic growth in Chlamydomonas reinhardtii","volume":"12","author":"Boyle","year":"2017","journal-title":"PLoS One"},{"key":"2023050421464706500_btad186-B4","doi-asserted-by":"crossref","first-page":"647","DOI":"10.1002\/bit.10803","article-title":"OptKnock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization","volume":"84","author":"Burgard","year":"2003","journal-title":"Biotechnol Bioeng"},{"key":"2023050421464706500_btad186-B5","doi-asserted-by":"crossref","first-page":"38","DOI":"10.1016\/j.ymben.2010.11.004","article-title":"Synergy between 13C-metabolic flux analysis and flux balance analysis for understanding metabolic adaption to anaerobiosis in E. coli","volume":"13","author":"Chen","year":"2011","journal-title":"Metab Eng"},{"key":"2023050421464706500_btad186-B6","doi-asserted-by":"crossref","first-page":"1050","DOI":"10.1111\/tpj.12252","article-title":"A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions","volume":"75","author":"Cheung","year":"2013","journal-title":"Plant J"},{"key":"2023050421464706500_btad186-B7","doi-asserted-by":"crossref","first-page":"1671","DOI":"10.1104\/pp.15.00880","article-title":"A method of accounting for enzyme costs in flux balance analysis reveals alternative pathways and metabolite stores in an illuminated Arabidopsis leaf","volume":"169","author":"Cheung","year":"2015","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B8","doi-asserted-by":"crossref","first-page":"231","DOI":"10.1016\/j.ymben.2021.12.003","article-title":"Metabolic flux analysis of the non-transitory starch tradeoff for lipid production in mature tobacco leaves","volume":"69","author":"Chu","year":"2022","journal-title":"Metab Eng"},{"key":"2023050421464706500_btad186-B27","doi-asserted-by":"crossref","first-page":"579","DOI":"10.1104\/pp.109.148817","article-title":"AraGEM, a genome-scale reconstruction of the primary metabolic network in Arabidopsis","volume":"152","author":"de Oliveira Dal\u2019Molin","year":"2010","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B28","doi-asserted-by":"crossref","first-page":"1871","DOI":"10.1104\/pp.110.166488","article-title":"C4GEM, a genome-scale metabolic model to study C4 plant metabolism","volume":"154","author":"de Oliveira Dal\u2019Molin","year":"2010","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B29","first-page":"1","article-title":"A multi-tissue genome-scale metabolic modeling framework for the analysis of whole plant systems","volume":"6","author":"de Oliveira Dal\u2019Molin","year":"2015","journal-title":"Front Plant Sci"},{"key":"2023050421464706500_btad186-B9","doi-asserted-by":"crossref","first-page":"1255","DOI":"10.1016\/j.cell.2019.11.009","article-title":"Conversion of Escherichia coli to generate all biomass carbon from CO2","volume":"179","author":"Gleizer","year":"2019","journal-title":"Cell"},{"key":"2023050421464706500_btad186-B10","author":"Gurobi Optimization, LLC","year":"2019"},{"key":"2023050421464706500_btad186-B11","doi-asserted-by":"crossref","first-page":"526","DOI":"10.1111\/j.1365-313X.2011.04613.x","article-title":"Metabolic network reconstruction and flux variability analysis of storage synthesis in developing oilseed rape (Brassica napus L.) embryos","volume":"67","author":"Hay","year":"2011","journal-title":"Plant J"},{"key":"2023050421464706500_btad186-B12","doi-asserted-by":"crossref","first-page":"639","DOI":"10.1038\/s41596-018-0098-2","article-title":"Creation and analysis of biochemical constraint-based models using the COBRA toolbox v.3.0","volume":"14","author":"Heirendt","year":"2019","journal-title":"Nat Protoc"},{"key":"2023050421464706500_btad186-B13","doi-asserted-by":"crossref","first-page":"593","DOI":"10.1126\/science.1132067","article-title":"Multiple high-throughput analyses monitor the response of E. coli to perturbations","volume":"316","author":"Ishii","year":"2007","journal-title":"Science"},{"key":"2023050421464706500_btad186-B14","doi-asserted-by":"crossref","first-page":"e1007099","DOI":"10.1371\/journal.pcbi.1007099","article-title":"Transcriptome-guided parsimonious flux analysis improves predictions with metabolic networks in complex environments","volume":"16","author":"Jenior","year":"2020","journal-title":"PLoS Comput Biol"},{"key":"2023050421464706500_btad186-B15","doi-asserted-by":"crossref","first-page":"1058","DOI":"10.1111\/tpj.13312","article-title":"A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling","volume":"88","author":"Klepikova","year":"2016","journal-title":"Plant J"},{"key":"2023050421464706500_btad186-B16","doi-asserted-by":"crossref","first-page":"349","DOI":"10.1016\/j.tplants.2004.05.001","article-title":"Dynamic flexibility in the light reactions of photosynthesis governed by both electron and proton transfer reactions","volume":"9","author":"Kramer","year":"2004","journal-title":"Trends Plant Sci"},{"key":"2023050421464706500_btad186-B17","doi-asserted-by":"crossref","first-page":"1034","DOI":"10.1021\/pr401017h","article-title":"Accurate label-free protein quantitation with high- and low-resolution mass spectrometers","volume":"13","author":"Krey","year":"2014","journal-title":"J Proteome Res"},{"key":"2023050421464706500_btad186-B18","doi-asserted-by":"crossref","first-page":"169","DOI":"10.1016\/S0014-5793(01)02613-8","article-title":"Transcriptome meets metabolome: hierarchical and metabolic regulation of the glycolytic pathway","volume":"500","author":"Kuile","year":"2001","journal-title":"FEBS Lett"},{"key":"2023050421464706500_btad186-B19","doi-asserted-by":"crossref","first-page":"149","DOI":"10.1038\/msb4100196","article-title":"Systems metabolic engineering of Escherichia coli for L-threonine production","volume":"3","author":"Lee","year":"2007","journal-title":"Mol Syst Biol"},{"key":"2023050421464706500_btad186-B20","doi-asserted-by":"crossref","first-page":"390","DOI":"10.1038\/msb.2010.47","article-title":"Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models","volume":"6","author":"Lewis","year":"2010","journal-title":"Mol Syst Biol"},{"key":"2023050421464706500_btad186-B21","doi-asserted-by":"crossref","first-page":"535","DOI":"10.1016\/j.cell.2016.03.014","article-title":"On the dependency of cellular protein levels on mRNA abundance","volume":"165","author":"Liu","year":"2016","journal-title":"Cell"},{"key":"2023050421464706500_btad186-B22","doi-asserted-by":"crossref","first-page":"16967","DOI":"10.1073\/pnas.1319485111","article-title":"Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation","volume":"111","author":"Ma","year":"2014","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023050421464706500_btad186-B23","doi-asserted-by":"crossref","first-page":"e1003580","DOI":"10.1371\/journal.pcbi.1003580","article-title":"Systematic evaluation of methods for integration of transcriptomic data into constraint-based models of metabolism","volume":"10","author":"Machado","year":"2014","journal-title":"PLoS Comput Biol"},{"key":"2023050421464706500_btad186-B24","doi-asserted-by":"crossref","first-page":"264","DOI":"10.1016\/j.ymben.2003.09.002","article-title":"The effects of alternate optimal solutions in constraint-based genome-scale metabolic models","volume":"5","author":"Mahadevan","year":"2003","journal-title":"Metab Eng"},{"key":"2023050421464706500_btad186-B25","doi-asserted-by":"crossref","first-page":"3966","DOI":"10.1016\/j.febslet.2009.10.036","article-title":"Correlation of mRNA and protein in complex biological samples","volume":"583","author":"Maier","year":"2009","journal-title":"FEBS Lett"},{"key":"2023050421464706500_btad186-B26","doi-asserted-by":"crossref","first-page":"409","DOI":"10.1038\/s41586-020-2094-2","article-title":"Mass-spectrometry-based draft of the Arabidopsis proteome","volume":"579","author":"Mergner","year":"2020","journal-title":"Nature"},{"key":"2023050421464706500_btad186-B30","doi-asserted-by":"crossref","first-page":"245","DOI":"10.1038\/nbt.1614","article-title":"What is flux balance analysis?","volume":"28","author":"Orth","year":"2010","journal-title":"Nat Biotechnol"},{"key":"2023050421464706500_btad186-B31","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1371\/journal.pcbi.1007036","article-title":"Enhanced flux prediction by integrating relative expression and relative metabolite abundance into thermodynamically consistent metabolic models","volume":"15","author":"Pandey","year":"2019","journal-title":"PLoS Comput Biol"},{"key":"2023050421464706500_btad186-B32","doi-asserted-by":"crossref","first-page":"7797","DOI":"10.1073\/pnas.0702609104","article-title":"Metabolic engineering of Escherichia coli for the production of L-valine based on transcriptome analysis and in silico gene knockout simulation","volume":"104","author":"Park","year":"2007","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023050421464706500_btad186-B33","doi-asserted-by":"crossref","first-page":"1570","DOI":"10.1104\/pp.109.141267","article-title":"A genome-scale metabolic model of Arabidopsis and some of its properties","volume":"151","author":"Poolman","year":"2009","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B34","doi-asserted-by":"crossref","first-page":"490","DOI":"10.1111\/j.1365-313X.2005.02649.x","article-title":"Measuring multiple fluxes through plant metabolic networks","volume":"45","author":"Ratcliffe","year":"2006","journal-title":"Plant J"},{"key":"2023050421464706500_btad186-B35","doi-asserted-by":"crossref","first-page":"e1009589","DOI":"10.1371\/journal.pcbi.1009589","article-title":"\u0394FBA\u2014predicting metabolic flux alterations using genome-scale metabolic models and differential transcriptomic data","volume":"17","author":"Ravi","year":"2021","journal-title":"PLoS Comput Biol"},{"key":"2023050421464706500_btad186-B36","doi-asserted-by":"crossref","first-page":"39","DOI":"10.1016\/j.copbio.2014.02.011","article-title":"Recent advances in the reconstruction of metabolic models and integration of omics data","volume":"29","author":"Saha","year":"2014","journal-title":"Curr Opin Biotechnol"},{"key":"2023050421464706500_btad186-B37","doi-asserted-by":"crossref","first-page":"e21784","DOI":"10.1371\/journal.pone.0021784","article-title":"Zea mays irs1563: a comprehensive genome-scale metabolic reconstruction of maize metabolism","volume":"6","author":"Saha","year":"2011","journal-title":"PLoS One"},{"key":"2023050421464706500_btad186-B38","doi-asserted-by":"crossref","first-page":"1","DOI":"10.3389\/fpls.2014.00668","article-title":"Transcript abundance on its own cannot be used to infer fluxes in central metabolism","volume":"5","author":"Schwender","year":"2014","journal-title":"Front Plant Sci"},{"key":"2023050421464706500_btad186-B39","doi-asserted-by":"crossref","first-page":"21","DOI":"10.1111\/plb.12869","article-title":"Malate valves: old shuttles with new perspectives","volume":"21","author":"Selinski","year":"2019","journal-title":"Plant Biol J"},{"key":"2023050421464706500_btad186-B40","doi-asserted-by":"crossref","first-page":"489","DOI":"10.1007\/s00018-019-03384-y","article-title":"Multi-tissue to whole plant metabolic modelling","volume":"77","author":"Shaw","year":"2020","journal-title":"Cell Mol Life Sci"},{"key":"2023050421464706500_btad186-B41","doi-asserted-by":"crossref","first-page":"694","DOI":"10.1105\/tpc.112.106989","article-title":"Metabolic fluxes in an illuminated Arabidopsis rosette","volume":"25","author":"Szecowka","year":"2013","journal-title":"Plant Cell"},{"key":"2023050421464706500_btad186-B42","doi-asserted-by":"crossref","first-page":"1596","DOI":"10.1104\/pp.105.062141","article-title":"In vivo respiratory metabolism of illuminated leaves","volume":"138","author":"Tcherkez","year":"2005","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B43","doi-asserted-by":"crossref","first-page":"536","DOI":"10.1093\/bioinformatics\/btp704","article-title":"Predicting metabolic engineering knockout strategies for chemical production: accounting for competing pathways","volume":"26","author":"Tepper","year":"2010","journal-title":"Bioinformatics"},{"key":"2023050421464706500_btad186-B44","doi-asserted-by":"crossref","first-page":"3882","DOI":"10.1093\/bioinformatics\/bty445","article-title":"Integrating proteomic or transcriptomic data into metabolic models using linear bound flux balance analysis","volume":"34","author":"Tian","year":"2018","journal-title":"Bioinformatics"},{"key":"2023050421464706500_btad186-B45","first-page":"1218","article-title":"Seeing the wood for the trees: a Forest of methods for optimization and omic-network integration in metabolic modelling","volume":"19","author":"Vijayakumar","year":"2017","journal-title":"Brief Bioinform"},{"key":"2023050421464706500_btad186-B46","doi-asserted-by":"crossref","first-page":"1","DOI":"10.3389\/fpls.2015.00167","article-title":"The relationship between leaf area growth and biomass accumulation in Arabidopsis thaliana","volume":"6","author":"Weraduwage","year":"2015","journal-title":"Front Plant Sci"},{"key":"2023050421464706500_btad186-B47","doi-asserted-by":"crossref","first-page":"311","DOI":"10.1104\/pp.110.158535","article-title":"A genome-scale metabolic model accurately predicts fluxes in central carbon metabolism under stress conditions","volume":"154","author":"Williams","year":"2010","journal-title":"Plant Physiol"},{"key":"2023050421464706500_btad186-B48","doi-asserted-by":"crossref","first-page":"e2121531119","DOI":"10.1073\/pnas.2121531119","article-title":"Reimport of carbon from cytosolic and vacuolar sugar pools into the Calvin-Benson cycle explains photosynthesis labeling anomalies","volume":"119","author":"Xu","year":"2022","journal-title":"Proc Natl Acad Sci USA"},{"key":"2023050421464706500_btad186-B49","doi-asserted-by":"crossref","first-page":"297","DOI":"10.1093\/plphys\/kiab076","article-title":"The metabolic origins of non-photorespiratory CO2 release during photosynthesis: a metabolic flux analysis","volume":"186","author":"Xu","year":"2021","journal-title":"Plant Physiol"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btad186\/49841467\/btad186.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/5\/btad186\/50204855\/btad186.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/39\/5\/btad186\/50204855\/btad186.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,5,4]],"date-time":"2023-05-04T17:47:39Z","timestamp":1683222459000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btad186\/7114031"}},"subtitle":[],"editor":[{"given":"Inanc","family":"Birol","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2023,4,11]]},"references-count":49,"journal-issue":{"issue":"5","published-print":{"date-parts":[[2023,5,4]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btad186","relation":{"has-preprint":[{"id-type":"doi","id":"10.1101\/2022.09.05.506655","asserted-by":"object"}]},"ISSN":["1367-4811"],"issn-type":[{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2023,5,1]]},"published":{"date-parts":[[2023,4,11]]},"article-number":"btad186"}}