Abstract
Individual tissues perform highly specialized metabolic functions to maintain whole-body homeostasis. Although Drosophila serves as a powerful model for studying human metabolic diseases, a lack of tissue-specific metabolic models makes it challenging to quantitatively assess the metabolic processes of individual tissues and disease models in this organism. To address this issue, we reconstructed 32 tissue-specific genome-scale metabolic models (GEMs) using pseudo-bulk single cell transcriptomics data, revealing distinct metabolic network structures across tissues. Leveraging enzyme kinetics and flux analyses, we predicted tissue-dependent metabolic pathway activities, recapitulating known tissue functions and identifying tissue-specific metabolic signatures, as supported by metabolite profiling. Moreover, to demonstrate the utility of tissue-specific GEMs in a disease context, we examined the effect of a high sugar diet (HSD) on muscle metabolism. Together with (13)C-glucose isotopic tracer studies, we identified glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as a rate-limiting enzyme in response to HSD. Mechanistically, the decreased GAPDH activity was linked to elevated NADH/NAD(+) ratio, caused by disturbed NAD(+) regeneration rates, and oxidation of GAPDH. Furthermore, we introduced a pathway flux index to predict and validate additionally perturbed pathways, including fructose and butanoate metabolism. Altogether, our results represent a significant advance in generating quantitative tissue-specific GEMs and flux analyses in Drosophila, highlighting their use for identifying dysregulated metabolic pathways and their regulation in a human disease model.