Context-dependent modification of PFKFB3 in hematopoietic stem cells promotes anaerobic glycolysis and ensures stress hematopoiesis

Metabolic pathways are plastic and rapidly change in response to stress or perturbation. Current metabolic profiling techniques require lysis of many cells, complicating the tracking of metabolic changes over time after stress in rare cells such as hematopoietic stem cells (HSCs). Here, we aimed to identify the key metabolic enzymes that define differences in glycolytic metabolism between steady-state and stress conditions in murine HSCs and elucidate their regulatory mechanisms. Through quantitative 13C metabolic flux analysis of glucose metabolism using high-sensitivity glucose tracing and mathematical modeling, we found that HSCs activate the glycolytic rate-limiting enzyme phosphofructokinase (PFK) during proliferation and oxidative phosphorylation (OXPHOS) inhibition. Real-time measurement of ATP levels in single HSCs demonstrated that proliferative stress or OXPHOS inhibition led to accelerated glycolysis via increased activity of PFKFB3, the enzyme regulating an allosteric PFK activator, within seconds to meet ATP requirements. Furthermore, varying stresses differentially activated PFKFB3 via PRMT1-dependent methylation during proliferative stress and via AMPK-dependent phosphorylation during OXPHOS inhibition. Overexpression of Pfkfb3 induced HSC proliferation and promoted differentiated cell production, whereas inhibition or loss of Pfkfb3 suppressed them. This study reveals the flexible and multilayered regulation of HSC glycolytic metabolism to sustain hematopoiesis under stress and provides techniques to better understand the physiological metabolism of rare hematopoietic cells.


Introduction
Activities governing nutrient requirements and metabolic pathways in individual cells maintain tissue homeostasis and respond to stress through metabolite production.ATP, produced via cytosolic glycolysis and mitochondrial oxidative phosphorylation (OXPHOS), is the universal energy currency of all organisms; it regulates all anabolic or catabolic cellular activities (Schirmer and Evans, 1990;Denton et al., 1975;Harris et al., 1997).Precise control of intracellular ATP concentrations is critical, as ATP is the rate determiner of many ATP-dependent biochemical reactions (Sols, 1981;Gabriel et al., 1985;Frieden, 1965;Hardie et al., 2012;Lin and Hardie, 2018;Hardie and Carling, 1997).
During HSC division, cell metabolism is reprogrammed to activate fatty acid β-oxidation (FAO) and purine metabolism (Ito et al., 2012;Karigane et al., 2016;Umemoto et al., 2022).Furthermore, Liang et al. reported that activated HSCs mainly rely on glycolysis as their energy source (Liang et al., 2020).However, the mechanisms by which each ATP-producing pathway and their connections are differentially regulated between HSCs and differentiated cells at steady state, during cell cycling, or during stress remain unknown.Recently, it has been shown that deeply quiescent HSCs do not activate cell cycle under stress (Bowling et al., 2020;Fanti et al., 2023;Munz et al., 2023).Therefore, it remains unclear whether metabolic changes such as the individual ATP-producing pathways and their interconnections occur uniformly in all HSCs, including these deeply quiescent HSCs.Furthermore, the underlying hub metabolic enzyme responsible for changes in the metabolic system of HSCs under stress has not been identified.HSCs are essential for cell therapy, including HSC transplantation, and in order to comprehensively elucidate the metabolic systems that have attracted attention as their regulatory mechanisms, recent studies have included metabolomic analyses using rare cell types such as HSCs (Qi et al., 2021;Agathocleous et al., 2017;DeVilbiss et al., 2021;Lengefeld et al., 2021;Schönberger et al., 2022), as well as isotope tracer analyses of undifferentiated hematopoietic cells purified after in vivo administration of isotopic glucose (Jun et al., 2021).Although these approaches are useful for obtaining comprehensive information on intracellular metabolites, they are not suited to track real-time changes in cellular metabolism at high resolution.Therefore, new approaches are necessary to analyze metabolites quantitatively and continuously without disturbing the physiological states of single cells while integrating the recently reported metabolome analysis techniques.In this study, we aimed to identify the key metabolic enzymes that define differences in glycolytic metabolism between steady-state and stress conditions in HSCs and elucidate their regulatory mechanisms using a quantitative and mathematical approach.Our findings provide a platform for quantitative metabolic analysis of rare cells such as HSCs, characterize the overall metabolic reprogramming of HSCs during stress loading, and highlight the key enzyme involved in this process.

HSC cell cycling increases anaerobic glycolytic flux
To determine how cell cycle progression alters HSC metabolism in vivo, we intraperitoneally and intravenously treated mice with 5-fluorouracil (5-FU) to induce HSC cell cycling (Figure 1-figure supplement 1A).For analysis after 5-FU administration, the Lineage (Lin) -Sca-1 + c-Kit + (LSK) gate was expanded to include HSCs with decreased c-Kit expression levels early after 5-FU treatment, for example high Sca-1-expressing cells and c-Kit-high to -dim Lin -cells, based on the previous report (Arai et al., 2004;Umemoto et al., 2022;Figure 1-figure supplement 1B).This expanded LSK gate was consistent with the patterns of c-Kit expression observed in endothelial protein C receptor (EPCR) + Lin -CD150 + CD48 -cells (Figure 1-figure supplement 1C) with high stem cell activity after 5-FU administration (Umemoto et al., 2022).We observed a transient decrease in the number of quiescent HSCs (Ki67 -) and an increase in the number of cell-cycling HSCs (Ki67 + ) on day 6 after 5-FU treatment (Figure 1-figure supplement 1D).Along with the loss of cell quiescence, ATP concentration in HSCs decreased transiently on day 6 (Figure 1-figure supplement 1E).Because the route of administration of 5-FU (intraperitoneal or intravenous) made no difference in the Ki67 positivity rate of HSCs (Figure 1-figure supplement 1F), we administered 5-FU intraperitoneally for remaining experiments.Two methods were used to test whether cell cycle progression of HSCs after 5-FU treatment depends on the expression of EPCR.First, phosphorylation of Rb (pRb), a marker of cell cycle progression (Miller et al., 2018), was analyzed in HSCs after 5-FU treatment.Analysis of EPCR + and EPCR -HSCs showed increased pRb in HSCs from 5-FU-treated mice in both fractions compared to HSCs from phosphate-buffered saline (PBS)-treated mice, regardless of EPCR expression (Figure 1figure supplement 1G-H).Second, we used a G 0 marker mouse line (Fukushima et al., 2019).These mice expressed a fusion protein of the p27 inactivation mutant p27K -and the fluorescent protein mVenus (G 0 marker), allowing prospective identification of G 0 cells.We tested whether the expression of G 0 marker in HSCs was altered after 5-FU administration to the G 0 marker mice (Figure 1-figure supplement 1I) and found that 5-FU treatment reduced the frequency of G 0 marker-positive HSCs,

OXPHOS-inhibited HSCs exhibit compensatory glycolytic flux
Previous studies using mouse models of mitochondrial disease or defects in genes involved in electron transport chain and OXPHOS suggest that mitochondrial energy production is essential for maintaining HSC function (Ansó et al., 2017;Inoue et al., 2010;Yu et al., 2013;Bejarano-García et al., 2016), as is the glycolytic system.However, there have been no quantitative reports on how OXPHOS-inhibited HSCs can adapt their metabolism.To understand HSC metabolism under OXPHOS inhibition, we performed in vitro U-13 C 6 -glucose tracer analysis of oligomycin-treated HSCs (Figure 2A; Supplementary file 3).Similar to 5-FU-treated HSCs (Figure 1), oligomycin-treated HSCs exhibited glycolytic system activation (Figure 2B-C

Phosphofructokinase (PFK) metabolism in HSCs is activated during proliferation and OXPHOS inhibition
To investigate whether glycolytic activation in HSCs after 5-FU treatment and OXPHOS inhibition could be demonstrated through unbiased mathematical simulations, we performed quantitative 13 C metabolic flux analysis ( 13 C-MFA).After generating a metabolic model for isotope labeling enrichment and setting appropriate lactate efflux values, a simulation was conducted using the labeled metabolite abundance data obtained from isotope tracer analysis.The appropriate lactate efflux for quiescent HSC (PBS-treated HSC) was determined to 65 after experimenting with values from 0 to 100.The lactate efflux of 5-FU-or oligomycin-treated HSCs was higher than that of quiescent HSCs based on the observation that labeled glycolytic metabolite levels were particularly elevated in in vitro tracer analysis (see 'Quantitative 13 C-MFA with OpenMebius' under 'Materials and methods' for more information).As a result, the variation in the flux values of all enzymatic reactions calculated in HSCs after 5-FU or oligomycin treatment became smaller compared to quiescent HSCs, suggesting that HSCs strictly regulated their metabolism in response to stress (Figure 3-figure supplement 1A-C).Unlike PBS-treated HSCs, those treated with 5-FU or oligomycin exhibited preferential glycolytic activation rather than TCA-or PPP-based metabolic strategies; the first half of the glycolytic system appeared to be the site of metabolic activation (Figure 3A        Figure 1 continued file 4).This increase in metabolic flux upstream of the glycolytic pathway was also supported by our in vitro tracer analysis (Figure 1B and Figure 2B), suggesting that 13 C-MFA was a valid metabolic simulation.Among the reactions in the first half of glycolysis, phosphorylation of fructose 6-phosphate (F6P) by PFK is the irreversible and rate-limiting reaction (Dunaway, 1983).A detailed review of in vitro isotope tracer analysis results showed that the ratio of fructose 1,6-bisphosphate (F1,6BP; the product of PFK) to F6P (the substrate of PFK) was greatly elevated in HSCs during proliferation and OXPHOS inhibition (Figure 3K-L).Together with the results of quantitative 13 C-MFA, these findings suggested that HSCs exhibit elevated glycolytic flux relative to mitochondrial activity by increasing PFK enzyme activity under various stress conditions.

HSCs under stress exhibit activation of glycolysis-initiated TCA cycle and NAS
To investigate the long-term glucose utilization of HSCs, we performed an in vivo tracer analysis with U-13 C 6 glucose based on recent reports (DeVilbiss et al., 2021;Jun et al., 2021;Figure 3-figure supplement 2A; see 'Preparation and storage of in vivo U-13 C 6 -glucose tracer samples' under 'Materials and methods' for more information).In HSCs from 5-FU-treated mice, we observed increased labeling of glycolytic metabolites such as dihydroxyacetone phosphate, glycerol-3-phosphate, and phosphoenolpyruvate, as well as NAS metabolites such as inosine monophosphate and ATP, and those derived from TCA cycle such as aspartic acid and glutamate, compared to HSCs from PBS-treated mice (Figure 3-figure supplement 2B-I, Supplementary file 5).When the amount of U-13 C 6 -glucosederived labeled metabolites in each pathway was calculated, more glucose-derived metabolites entered TCA cycle in the 5-FU-treated group than PBS-treated group (Figure 3-figure supplement 2J).Thus, although short-term (10-30 min) in vitro tracer analysis showed that HSCs exhibited more potent activation of anaerobic glycolysis than of other pathways in response to 5-FU administration, long-term (approximately 3 hr) labeling by in vivo tracer analysis revealed that glycolysis-initiated TCA cycle and NAS flux were activated in addition to enhanced anaerobic glycolysis.Importantly, despite differences in labeling times and supplementation of U-13 C 6 glucose metabolites from non-HSCs to HSCs in vivo, the activation of the glycolytic system was a common finding.

PFKFB3 accelerates glycolytic ATP production during HSC cell cycling
In vitro and in vivo tracer analysis results collectively suggested that the activation of glycolysis catalyzed by PFK may have been the starting point for the activation of the entire HSC metabolism.To analyze the contribution of PFK to ATP metabolism in steady-state or stressed HSCs, we needed to develop an experimental system that could measure the dynamics of ATP concentrations in HSCs in a non-destructive, real-time manner.To this end, we used knock-in GO-ATeam2 mice as a FRET-based biosensor of ATP concentration (see 'Conversion of GO-ATeam2 fluorescence to ATP concentration' under 'Materials and methods' for more information.).The number of bone marrow mononuclear cells (BMMNCs), as well as the frequency of HSCs (CD150 + CD48 -LSK) and other progenitor cells, in the bone marrow (BM) of GO-ATeam2 + mice were almost unchanged compared to C57BL/6J mice, except for a mild decrease in the Lin -fraction (Figure 4-figure supplement 1A-C).Using BMMNCs derived from GO-ATeam2 + mice, we developed a method to detect changes in ATP concentration with high temporal resolution when the activity of PFK was modulated (Figure 4-figure supplement 1D-F).suppressed this decrease in MyPs; however, other metabolites (e.g.pyruvate, lactate, and fatty acids) had minimal effects, suggesting that ATP levels are glycolysis-dependent in MyPs (Figure 4-figure supplement 1G-H), consistent with previous reports that the aerobic glycolytic enzyme M2 pyruvate kinase isoform (PKM2) is required for progenitor cell function (Wang et al., 2014).Further, we analyzed ATP consumption and metabolic dependency of cell-cycling HSCs after 5-FU administration (Figure 4A).After inhibiting glycolysis using 2-deoxy-D-glucose (2-DG) with other mitochondrial substrates, 5-FU-treated HSCs showed more rapid decreases in ATP concentration than PBS-treated HSCs (Figure 4B-C).In contrast, OXPHOS inhibition by oligomycin without glucose or mitochondrial substrates decreased the ATP concentration to a similar extent in both 5-FU-and PBStreated HSCs, although 5-FU-treated HSCs showed earlier ATP exhaustion (Figure 4D-E).These data suggest that 5-FU-treated-HSCs upregulated ATP production via glycolysis, rather than relying on mitochondria.Apoptosis assay revealed a slight increase in early apoptotic cells (annexin V + propidium iodide [PI] -) after 2-DG treatment and a slight decrease in the number of viable cells (Annexin V -PI -) after oligomycin treatment, both to a very limited extent (approximately 5%) compared to the degree of ATP decrease, suggesting that the decrease in ATP after 2-DG or oligomycin treatment did not simply reflect cell death (Figure 4-figure supplement 1I).Importantly, no metabolic changes in glycolysis or OXPHOS were observed in HSCs without cell cycle progression after 5-FU administration (very early phase: day 3; late phase: day 15) (Figure 4-figure supplement 2A-H).
Finally, we investigated the nutrients that drive OXPHOS in PBS-or 5-FU-treated HSCs.Exposure of PBS-or 5-FU-treated HSCs to either etomoxir, a FAO inhibitor, or 6-diazo-5-oxo-L-norleucine (DON), a glutaminolysis inhibitor, alone or in combination, did not decrease ATP concentrations (Figure 4figure supplement 2J-M).Subsequent assessment of FAO activity using FAOBlue, a fluorescent probe for the FAO activity assay (Uchinomiya et al., 2020,) showed no significant differences between PBSand 5-FU-treated HSCs (Figure 4-figure supplement 2N).Thus, neither FAO nor glutaminolysis appeared to be essential for the short-term maintenance of ATP levels in cell-cycling HSCs after 5-FU administration.Notably, the addition of glucose and a PFKFB3 inhibitor to etomoxir rapidly reduced ATP concentrations in HSCs (Figure 4-figure supplement 2O-P).This suggests that etomoxir may partially mimic the effects of oligomycin, indicating that OXPHOS is primarily driven by FAO, but can be compensated by PFKFB3-accelerated glycolysis in HSCs.Conversely, exposure of HSCs to DON in combination with a PFKFB3 inhibitor did not decrease ATP concentrations (Figure 4-figure supplement 2O-P), suggesting that ATP production via glutaminolysis is limited in HSCs.

OXPHOS inhibition accelerates glycolysis to sustain ATP levels in HSCs, but not in progenitors
To assess differences in metabolic dependence between steady-state or stressed HSCs and naturally proliferating HPCs, we altered ATP metabolism in HSCs and progenitors using 2-DG or oligomycin (Figure 5A).Oligomycin treatment rapidly depleted ATP in HSCs and all HPC fractions (green lines in The ATP reduction induced by 2-DG treatment was particularly low (~15%) in HSCs, multipotent progenitor cells (MPPs), and common lymphoid progenitors (CLPs) relative to that in common myeloid progenitors (CMPs), granulocytes-macrophage progenitors (GMPs), and megakaryocyte-erythrocyte progenitors (MEPs; Figure 5D).
Next, we investigated the role of glycolysis in ATP production during OXPHOS inhibition by combining oligomycin administration and glucose supplementation.ATP concentration remained   5E).Similar to oligomycin treatment, rotenone (complex I inhibitor) and carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone (FCCP, mitochondrial uncoupler) treatments, which inhibit OXPHOS-derived ATP production, also decreased ATP concentrations in HSCs, but not when administered simultaneously with glucose (Figure 5-figure supplement 1E-F).Furthermore, with oligomycin, HSCs, but not HPCs, maintained ATP concentrations at low glucose levels (50 mg/dL) (Figure 5-figure supplement 1G).These analyses suggest that ATP was produced by mitochondrial OXPHOS in steady-state HSCs, and that only HSCs, but not HPCs, maintained ATP production by glycolysis when OXPHOS was compromised.

PFKFB3 accelerates glycolytic ATP production during OXPHOS inhibition
Next, to understand whether PFKFB3 contributes to ATP production in HSCs under OXPHOS inhibition, we evaluated PFKFB3 function under OXPHOS inhibition using the GO-ATeam2 + BMMNCs.In oligomycin-treated HSCs, PFKFB3 inhibition led to rapidly decreased ATP concentration that was not observed in HSCs not treated with oligomycin (Figure 5F-I).We examined the effects of HSPC metabolic regulators on ATP levels in oligomycin-treated HSCs.Inhibiting PKM2, which accelerates glycolysis in steady-state progenitors (Wang et al., 2014), significantly reduced ATP levels in oligomycin-treated HSCs (Figure 5-figure supplement 1H, J).Inhibiting LKB1, a kinase upstream of AMPK (Hardie, 2014;Long and Zierath, 2006), did not affect the ATP concentration in oligomycintreated HSCs (Figure 5-figure supplement 1I, K), whereas levels of adenosine monophosphate (AMP), which also activates AMPK, increased in oligomycin-treated but not in 5-FU-treated HSCs (Figure 5-figure supplement 1L).This may explain differences in AMPK-dependent ATP production between proliferative HSCs and HSCs under OXPHOS inhibition.
Next, we tested the effects of PFKFB3 on ATP concentration in HPCs.Unlike HSCs, HPCs exhibited PFKFB3-dependent ATP production, even without oligomycin (Figure 5-figure supplement 1M-Q).Therefore, ATP production in steady-state HSCs was PFKFB3-independent, and proliferative stimulation or OXPHOS inhibition plastically activated glycolytic ATP production in a PFKFB3-dependent manner to meet ATP demand.

PFKFB3 activity renders HSCs dependent on glycolysis
Next, we investigated whether PFKFB3 activity itself confers glycolytic dependence on HSCs.We retrovirally overexpressed Pfkfb3 in HSCs and performed cell cycle analysis (Figure 5J).Pfkfb3overexpressed HSCs increased the proportion of cells in the S/G2/M phase and decreased the number of G 0 cells compared to mock-overexpressed HSCs (Figure 5K).Next, we retrovirally overexpressed Pfkfb3 in GO-ATeam2 + HSCs and performed real-time ATP measurement (Figure 5J).Pfkfb3overexpressing GO-ATeam2 + HSCs did not show changes in ATP concentrations relative to those in mock-transduced cells (Figure 5L;

PFKFB3 methylation by PRMT1 supports ATP production by cell-cycling HSCs
Next, we investigated how 5-FU-treated-HSCs regulate PFKFB3 independently of AMPK (Figure 4G-J).PFKFB3 activity is regulated at multiple levels (Shi et al., 2017), and PFKFB3 transcript and protein levels in HSCs remained unchanged during 5-FU-induced cell cycling (Figure 6A-B).Phosphorylation can also regulate PFKFB3 activity (Marsin et al., 2002;Novellasdemunt et al., 2013;Okamura and Sakakibara, 1998); however, we observed no change in PFKFB3 phosphorylation in 5-FU-treated-HSCs (Figure 6C).Upon oligomycin exposure, PFKFB3 was phosphorylated by AMPK in the HSCs (Figure 6D).PFKFB3 is also methylated, and its activity is upregulated by protein arginine methyltransferase 1 (PRMT1; Yamamoto et al., 2014).We observed that Prmt1 expression increased in 5-FU-treated-HSCs relative to that in PBS-treated-HSCs (Figure 6E).Furthermore, PFKFB3 methylation was significantly induced in 5-FU-treated-HSCs than in PBS-treated-HSCs (Figure 6F).Treatment of HSCs with a PRMT1 inhibitor decreased PFKFB3 methylation (Figure 6G), suggesting that PRMT1 catalyzed PFKFB3 methylation.In contrast, the number of transcripts regulated by PRMT1 decreased or was unchanged (Figure 6-figure supplement 1), suggesting that the transcriptional regulatory function of PRMT1 is limited.To investigate whether glycolytic activity in HSCs was regulated by methylated-PFKFB3 (m-PFKFB3), mice treated with PBS or 5-FU were injected with 2-NBDG, and m-PFKFB3 levels in HSCs with high and low 2-NBDG uptake were quantified.Regardless of PBS or 5-FU treatment, HSCs with high 2-NBDG uptake exhibited higher m-PFKFB3 levels than those with low uptake (Figure 6H), suggesting that m-PFKFB3 regulated the activity of the glycolytic system in HSCs.
Further, we analyzed the potential effects of PRMT1 inhibition on ATP concentration in GO-ATeam2 + HSCs.Treatment with the PRMT1 inhibitor significantly decreased ATP levels in 5-FU-treated-HSCs than in PBS-treated-HSCs (Figure 6I).In contrast, the retroviral overexpression of Pfkfb3 in GO-ATeam2 + HSCs abolished the effect of the PRMT1 inhibitor on ATP reduction (Figure 6J).These findings indicated that ATP levels in 5-FU-treated-HSCs were supported by PRMT1 methylation-mediated PFKFB3 activation.

PFKFB3 contributes to HSPC pool expansion and stress hematopoiesis maintenance
Finally, we analyzed PFKFB3 function in HSCs during hematopoiesis.We cultured HSCs with a PFKFB3 inhibitor in vitro under quiescence-maintaining or proliferative conditions (Figure 7-figure supplement 1A; Kobayashi et al., 2019).Cell count in HSC-derived colonies decreased following treatment with a PFKFB3 inhibitor under proliferative, but not quiescence-maintaining, conditions (Figure 7figure supplement 1B).We also knocked out Pfkfb3 in HSCs using the less toxic, vector-free CRISPR-Cas9 system and cultured the cells under quiescence-maintaining or proliferative conditions (Figure 7-figure supplement 1A) based on recent reports by Shiroshita et al., 2022.Again, cell numbers in Pfkfb3-knockout (KO) HSC-derived colonies decreased only in proliferative cultures when compared to control cultures (Rosa26-KO HSCs) (Figure 7-figure supplement 1C, E, F).We retrovirally overexpressed Pfkfb3 in HSCs and cultured them under quiescence maintenance or proliferative conditions (Figure 7-figure supplement 1A).Pfkfb3-overexpressing HSC colonies showed increased cell count compared to that of mock-transduced cells, but only under proliferative conditions (Figure 7-figure supplement 1D).
To assess PFKFB3 function in HSCs in vivo, we transplanted Pfkfb3-KO HSCs (Ly5.2 + ) or wild type (WT) control HSCs into lethally irradiated recipients (Ly5.1 + ) as well as Ly5.1 + competitor cells (Figure 7A), and the behavior of Pfkfb3-KO cells was evaluated by Sanger sequencing of peripheral blood (PB) cells (Shiroshita et al., 2022).In the KO group, donor-derived chimerism in PB cells decreased relative to that in the WT control group during the early phase (1 month post-transplant) but recovered thereafter (Figure 7B).Next, we retrovirally transduced Ly5.2 + HSCs with Pfkfb3 S461E (Pfkfb3CA), a constitutively active PFKFB3 mutant, and transplanted them into lethally irradiated recipients (Ly5.2 + ), along with Ly5.1 + competitor cells (Figure 7A, Figure 7-figure supplement 1G).Donor chimerism during the early post-transplant period in the Pfkfb3CA-overexpressing group was significantly higher than that in the mock-transduced group (Figure 7C).These findings suggest that PFKFB3 may play a role in the differentiation and proliferation of HSCs.Therefore, we compared the contribution of PFKFB3 to HSPC function at steady state and after myeloproliferative stimulation.Pfkfb3-or Rosa26-KO HSPCs were transplanted into recipients (Ly5.1 + ).After 2 months, recipients received 5-FU intraperitoneally, and the dynamics of Pfkfb3-or Rosa26-KO cell abundance in PB was assessed (Figure 7D).In PB cells prior to 5-FU administration, Pfkfb3-or Rosa26-KO HSPC-derived blood cells were almost equally present, suggesting a limited involvement of PFKFB3 in steady-state blood cell production (Figure 7E).However, after 5-FU administration, Pfkfb3-KO HSPC-derived blood cell abundance was reduced compared to that in the Rosa26-KO group (Figure 7E).This change occurred on day 6 after 5-FU administration (day 1), when the cell cycle of HSCs was activated (Figure 1-figure supplement 1D), supporting the idea that PFKFB3 contributes to HSC proliferation and differentiation into HSPCs.
To investigate the mechanisms underlying the short-term effects of PFKFB3 on hematopoiesis after bone marrow transplantation (BMT), we evaluated cell cycle and apoptosis of Pfkfb3-KO or    -overexpressing HSPCs on day 2 after BMT (Figure 7F).Cell cycle was analyzed by Ki67/Hoechst33432 staining and in vivo BrdU labeling (Jun et al., 2021), which showed that cell cycle progression was suppressed in Pfkfb3-KO HSPCs (Figure 7G-J).In contrast, Pfkfb3-KO cells did not show increased apoptotic rates or decreased homing efficiency after BMT (Figure 7K; Figure 7-figure supplement 1H).Furthermore, we examined the cell cycle of HSPCs overexpressing Pfkfb3CA on day 2 after BMT (Figure 7L) and found that Pfkfb3CA-overexpressing HSPCs showed accelerated cell cycle compared to mock-overexpressing HSPCs (Figure 7M-N).These data suggest that PFKFB3 contributes to HSC proliferation and differentiates cell production in in vitro and in vivo proliferative environments (cytokine stimulation and transplantation).

Discussion
In this study, by combining metabolomic tracing of U-13 C 6 -labeled glucose and 13 C-MFA, we quantitatively identified the metabolic programs used by HSCs during steady-state, cell-cycling, and OXPHOS inhibition.Under proliferative stress, HSCs uniformly shift from mitochondrial respiration to glycolytic ATP production and PPP activation, which represent hallmarks of cell-cycling mammalian cells (Intlekofer and Finley, 2019).Previous reports have emphasized the importance of glycolysis in maintaining HSC quiescence, but have primarily analyzed HSCs in transplant assays, wherein HSCs must enter the cell cycle (Takubo et al., 2013;Takubo et al., 2010).Prior analysis of repopulation capacity, which is positively correlated with enhanced glycolysis, may have overestimated glycolytic ATP production and overlooked mitochondrial ATP production during native hematopoiesis.In fact, some studies have suggested that OXPHOS activity is important for HSC maintenance and function (Ansó et al., 2017).
Our method was based on recently reported quantitative metabolic analysis techniques for very small numbers of cells (Qi et al., 2021;Agathocleous et al., 2017;DeVilbiss et al., 2021;Lengefeld et al., 2021;Schönberger et al., 2022;Jun et al., 2021), such as HSCs, and expands our knowledge of HSC metabolism during stress hematopoiesis.In our study, 5-FU administration in mice transiently decreased ATP concentration in HSCs in parallel with cell cycle progression, suggesting that HSC differentiation and cell cycle progression are closely related to intracellular metabolism and can be monitored by measuring ATP concentration.We mainly analyzed a mixture of EPCR + and EPCR - HSCs, and we believe that the observed cell cycle progression and promotion of glycolysis in both EPCR + and EPCR -HSCs support the validity of our claims (Figure 1L, Figure 1-figure supplement 1G-K).According to 13 C-MFA enzymatic reaction flux of PFK in 5-FU-treated HSCs indicated a relative increase of approximately 10%.However, the flux value obtained by 13 C-MFA was calculated with glucose uptake as 100.Thus, when combined with the overall increase in the glycolytic pool demonstrated by in vitro isotopic glucose tracer analysis and in vivo NBDG analysis, rapid acceleration of glycolysis becomes evident throughout the HSCs, including subpopulations that were less responsive to stress (Bowling et al., 2020;Fanti et al., 2023;Munz et al., 2023).These findings are consistent with reports suggesting that HSCs have relatively low biosynthetic activity (Signer et al., 2014;Essers et al., 2009) that is rapidly activated in response to cell proliferation stimuli (Karigane et al., 2016;Umemoto et al., 2018).Notably, we found that HSCs could accelerate glycolytic ATP production to fully compensate for mitochondrial ATP production under OXPHOS inhibition, a phenomenon that is difficult to identify without real-time ATP analysis.Thus, HSCs exposed to acute stresses choose to schema (F).Representative plots of Ki67/Hoechst33432 staining of Rosa26-KO (G) or Pfkfb3-KO (H) HSPCs and summary of analysis (I); summary of in vivo BrdU labeling assay (J).Apoptosis assay results (K).n=4-5 biological replicates for each group.(L-N) Cell cycle analysis of Pfkfb3CA or Mockoverexpressing HSPCs on day 2 after BMT.Experimental Schema (L).Representative plot of Ki67/Hoechst33432 staining for both groups (M) and summary of analysis (N).n=5 biological replicates for each group.(O) Models showing ATP production and regulation in quiescent, OXPHOS-inhibited, and cell-cycling HSCs.Note that the GO-ATeam2 system identified plastic acceleration of glycolysis by PFKFB3 in response to different types of stress maintains ATP levels.Data are presented as mean ± SD. * p≤0.05, ** p≤0.01, *** p≤0.001 as determined by Student's t-test (B, C, E, I-K, and N).See also    change the efficiency of glucose utilization (accelerated glycolytic ATP production) rather than other energy sources.In vivo, a completely glucose-deficient environment is improbable.Therefore, even under conditions such as hypoxia, where OXPHOS is inhibited, it is conceivable that glycolysis is accelerated to maintain ATP concentrations.Glucose tracer analysis showed NAS suppression under OXPHOS inhibition, leading to glycolysis without cell proliferation (Figure 2C-F; Figure 1-figure supplement 1N).This suppression can be attributed to several factors: phosphates derived from ATP are added to nucleotide mono-/di-phosphates during NAS; the primary source of ATP production, OXPHOS, is impaired; and the presence of enzymes, such as dihydroorotate dehydrogenase, which are conjugated with OXPHOS (Liu et al., 2000).Such multifactorial effects raise new questions about the relationship between OXPHOS and nucleotide synthesis.On the other hand, we observed that ATP production in steady-state or cell-cycling HSCs and in naturally proliferating HPCs depended more on mitochondrial OXPHOS than on glycolysis; inhibiting glycolysis in steady-state HSCs resulted in only mild ATP decreases, suggesting that OXPHOS is still the major source of ATP production even in a medium saturated with hypoxia mimicking the BM environment.The p50 value of mitochondria (the partial pressure of oxygen at which respiration is half maximal) is less than 0.1 kPa, corresponding to an oxygen concentration of less than 0.1% under atmospheric pressure (Gnaiger et al., 2000), suggesting that even under hypoxic conditions, OXPHOS can maintain some level of activity.Because FAO and the mitochondrial respiratory chain are necessary for HSC self-renewal and quiescence (Ansó et al., 2017;Bejarano-García et al., 2016;Ito et al., 2012;Kobayashi et al., 2019), fatty acids may support mitochondrial ATP production independently of fluxes from glycolysis.FAO and glutaminolysis were not immediately essential for ATP production in HSCs.Given reports on the longterm necessity of FAO and glutaminolysis for HSC maintenance (Ito et al., 2012;Oburoglu et al., 2014), ATP concentrations could be maintained in the short term by compensatory pathways.Furthermore, although glycolysis and TCA cycle are decoupled in steady-state HSCs, in response to cell cycle progression, anaerobic glycolytic metabolism in HSCs is enhanced (Figure 1) and fluxes to TCA cycle and PPP from the glycolytic system are also promoted (Figure 3-figure supplement 2).The degree of glycolysis and TCA cycle coupling observed by in vitro and in vivo tracer analysis differed, likely due to differences in labeling time (10-30 min in vitro and 3 hr in vivo).In particular, in vivo tracer analysis allows all cells to be capable of metabolizing U-13 C 6 -glucose and providing its metabolites to HSCs, and there is a significant amount of time, approximately 120-180 min, after glucose labeling to purify HSCs.Metabolic reactions will continue during this time and subsequent processing on ice, which may increase the influx of labeled carbon into the TCA cycle.This complex dynamic in the in vivo tracer analysis makes it difficult to determine whether the labeled carbon influx is the result of direct influx from glycolysis or the re-uptake of metabolites by HSCs that have been processed by other cells.This is in contrast to in vitro analysis where such extended metabolic processing does not occur.Furthermore, despite an increased carbon influx into the TCA cycle in vivo, ATP production from mitochondria does not show a corresponding increase after 5-FU treatment, as shown by the GO-ATeam2 analysis shown in Figure 4C.Despite these technical differences, an essential common finding from both in vivo and in vitro analyses is the activation of glycolysis and nucleotide synthesis (NAS) in 5-FU-treated HSCs, highlighting critical metabolic changes in response to treatment.Moreover, these data provide direct evidence that glycolysis and TCA cycle become functionally uncoupled in quiescent HSCs (Takubo et al., 2013;Halvarsson et al., 2017).Our findings are also consistent with previous reports of OXPHOS activation associated with HSC proliferation (Takubo et al., 2013;Yu et al., 2013;Maryanovich et al., 2015;Ito et al., 2012).In other words, HSCs exhibit an increased proportion of anaerobic glycolysis-derived ATP by PFKFB3 upon proliferation and OXPHOS inhibition; furthermore, the glycolytic system is the starting point of metabolic activation and is indispensable for the overall enhancement of HSC metabolism (Figure 7H).
HPCs and leukemic cells accelerate glycolytic ATP production using PKM2 for differentiation and transformation, respectively Wang et al., 2014; however, we demonstrated that glycolytic acceleration does not fully compensate for mitochondrial ATP production in HPCs.Mechanistically, PFKFB3 increased glycolytic activity in HSCs to maintain ATP concentrations during proliferation and OXPHOS inhibition.Furthermore, inhibition of PFKFB3 in addition to OXPHOS does not result in a complete loss of ATP in HSCs, suggesting the robustness of HSC metabolism (Figure 5G).Under steady-state conditions, naturally proliferating HPCs rely on PFKFB3 for ATP production, whereas HSCs do not.This may explain the reduction of ECAR after oligomycin treatment in MyPs as shown by the Mito stress test (Figure 2G).In other words, while PFKFB3-dependent active glycolysis and mitochondria must always be coupled in MyPs, this is not necessarily the case in HSCs, even after 5-FU treatment (Figure 1G).Therefore, we can infer that quiescent HSCs at steady state can produce ATP via PFKFB3 activation in response to stress, enabling additional ATP generation.Furthermore, overexpression of Pfkfb3 in HSCs increased glycolytic dependency, suggesting that PFKFB3 itself can modulate metabolic dependency in HSCs.Changes in glycolytic dependency in HSCs overexpressing Pfkfb3 may seem small (0.06-0.13 mM; Figure 5L and M).However, it is noteworthy that the rate of the reaction catalyzed by PFK varies greatly within a very narrow range of ATP concentrations, less than 1 mM.Webb et al. analyzed the factors controlling PFK activity and reported that the reaction rate of PFK varies by approximately 40% in the 0.3-1 mM ATP concentration range (Webb et al., 2015).The reason that differences in glycolytic dependence could be detected in cells overexpressing Pfkfb3 may be that the ATP concentration at the time of analysis was approximately 0.5-0.6 mM, which is within the range where a small change in ATP concentration can dynamically alter PFK activity.
PFKFB3 supports hematopoiesis in contexts that require robust HSPC proliferation in vitro and in vivo.We showed that the positive or negative effect of Pfkfb3 overexpression or KO on differentiated blood cell production is gradually lost after BMT.This is because HSPCs require PFKFB3 for cell cycle progression during stress hematopoiesis in the early phase after BMT (Figure 7F-J and  L-N).However, even during stress hematopoiesis, PFKFB3 is not involved in cell death or homing efficiency (Figure 7K; Figure 7-figure supplement 1H) and appears to contribute primarily to the regulation of transient HSPC proliferation in the BM cavity.HSCs no longer require PFKFB3 for a certain period of time after BMT, probably because they regain a quiescent state.This is consistent with the fact that inhibition of PFKFB3 in quiescent HSCs does not reduce the ATP concentration (Figure 5F and H), suggesting that the activity of PFKFB3 is plastically modified.HSC metabolic plasticity is also illustrated by the mode of PFKFB3 activation, differing depending on stress type.During proliferative stress, PRMT1 methylates PFKFB3 in the HSCs to promote glycolytic ATP production, a modification that increases its activity (Yamamoto et al., 2014).PRMT1 is required for stress hematopoiesis (Zhu et al., 2019), but its downstream targets in HSCs remain unclear.Our results strongly suggest that PRMT1 targets PFKFB3 to stimulate glycolysis in HSCs.In contrast, under OXPHOS inhibition, PFKFB3 phosphorylation by AMPK is induced-another modification that also upregulates its activity.These two PFKFB3 protein modifications allow for flexible regulation of ATP production by glycolysis, even under simultaneous and different stresses.In fact, the constitutively active S461E PFKFB3 mutant, designed to mimic phosphorylation in response to OXPHOS inhibition, enhanced HSC reconstitution capacity after transplantation, suggesting that even if PFKFB3 is activated by one stress (in this case, proliferative), it has the activation capacity to respond to a different stress (i.e.mitochondrial).Therefore, the functions of phosphorylated and methylated forms of PFKFB3 are to some extent interchangeable, and either modification can be used to handle diverse stresses.
In summary, we found that HSCs exhibit a highly dynamic range of glycolytic flux.Our study highlights glycolysis as a pivotal source of energy production in stressed HSCs, and indicates that OXPHOS, although an important source of ATP, can be uncoupled from glycolysis in steady-state HSCs without compromising ATP levels.Because multiple PFKFB3 modifications safeguard HSCs against different stresses by accelerating glycolysis, interventions targeting these might effectively induce or manage stress hematopoiesis.This study provides a platform for comprehensive and quantitative real-time analysis of ATP concentration and its dynamics in HSPCs.Our approach allows for analysis of metabolic programs in rare cells and detection of various metabolic activities within a diverse cell population, making it applicable to the analysis of various tissue systems in normal and diseased states.

Limitations of the study
In this study, 5-FU-treated HSCs were analyzed as cell-cycling HSCs, but if more sensitive and timesaving glucose tracer analysis methods (especially after in vivo labeling with isotopic glucose) are developed, it may be possible to prospectively differentiate and quantitatively analyze HSC metabolism based on the cell surface antigens and cell cycle status.Although our assay uses media that mimic the BM environment, in the near future, in vivo GO-ATeam2 analysis will allow us to measure ATP concentrations in physiologically hypoxic BM.

Mice and genotyping
C57BL/6 mice (7-16 weeks old, Ly5.2 + ) were purchased from Japan SLC (Shizuoka, Japan).C57BL/6 mice (Ly5.1 + ) were purchased from CLEA Japan (Shizuoka, Japan).Knock-in mice harboring GO-ATeam2 (Imamura et al., 2009;Nakano et al., 2011;Yamamoto et al., 2019) in the Rosa26 locus were generated in the Yamamoto laboratory.The GO-ATeam2 mice (8-16 weeks old) were used to analyze HSPCs.Ubc-GFP reporter mice (Ubc-GFP mice) were from the Jackson Laboratory and genotyped using PCR-based assays.GO-ATeam2 mice were genotyped by PCR of tail DNA or by transdermal GFP fluorescence.The PCR protocol was as follows: 94 °C for 5 min; 34 cycles of 94 °C for 30 s, 56 °C for 30 s, 72 °C for 30 s; 72 °C for 5 min; and 4 °C hold.Primers for GO-ATeam2 or Ubc-GFP mice are listed in Supplementary file 6. mVenus-p27K -mice (17-20 weeks old) were provided by Kitamura Laboratory and used for cell cycle analysis (Fukushima et al., 2019;Oki et al., 2014).Mice were genotyped using PCR-based assays of tail DNA or transdermal EGFP fluorescence.All mice were maintained in the animal facility at the National Center for Global Health and Medicine Research Institute under specific pathogen-free conditions and fed ad libitum.Mice were euthanized by cervical dislocation.All animal experiments were approved by the Institutional Animal Care and Use Committee (IACUC) at the National Center for Global Health and Medicine Research Institute.Both male and female mice were used.

Cell preparation
For C57BL/6 mice, bone marrow (BM) cells were isolated from bilateral femurs and tibiae by flushing with PBS + 2% fetal calf serum (FCS) (Gibco) using a 21-gauge needle (Terumo Corporation, Tokyo, Japan) and a 10 mL syringe (Terumo).As an exception, for U-13 C 6 -labeled glucose tracer experiments using C57BL/6 mice, BM was flushed with PBS +0.1% bovine serum albumin (BSA, Cat# A4503).The BM plug was dispersed by refluxing through the needle, and the suspension was centrifuged 680 × g for 5 min at 4 °C.Cells were lysed with lysis buffer (0.17 M NH 4 Cl, 1 mM EDTA, 10 mM NaHCO 3 ) at room temperature (RT) for 5 min, washed with two volumes PBS + 2% FCS (or PBS +0.1% BSA for tracer experiments), and centrifuged at 680 × g for 5 min at 4 °C.Cells were resuspended in PBS +2% FCS (or PBS +0.1% BSA for tracer experiments) and filtered through 40 µm nylon mesh (BD Biosciences).Cells were again centrifuged 680 × g for 5 min at 4 °C and treated with anti-CD16/32 antibody for Fc-receptor block (2 µL/mouse; BD Biosciences, Cat# 553152) for 10 min at 4 °C.Antic-Kit magnetic beads (Miltenyi Biotec,Bergisch Gladbach,Germany, were added at a 1:5 v/v ratio for 15 min at 4 °C.After removing the antibody with two PBS +2% FCS (or PBS +0.1% BSA for tracer experiments) washes, c-Kit-positive cells were isolated using Auto-MACS Pro (Miltenyi Biotec) with the Possel-s or Possel-d2 program.Isolated cells were centrifuged once at 340 × g for 5 min and stained with an antibody cocktail for flow cytometry.
For analysis of the GO-ATeam2 hematopoietic cells, BM from GO-ATeam2 mice was flushed with PBS + 0.1% BSA to minimize exposure to nutrients in FCS.Hemolysis, centrifugation, filtering, and Fc receptor blocking were performed in the same manner as for cell preparation using C57BL/6 mice.Cells were stained for 30 min with an antibody cocktail at 4 °C and then washed and suspended in 1000 µL PBS +0.1% BSA and centrifuged at 340 × g at 4 °C for 5 min.Supernatants were discarded in preparation for flow cytometry.

Intracellular staining for phosphorylated Rb (pRb)
EPCR + or EPCR -LSK-SLAM cells from PBS-or 5-FU-treated C57BL/6 mice were purified separately (see "Flow cytometry and cell sorting" for details).Anti-phospho-Rb (Ser807/811) antibody (CST, Cat# 8516T) was used as the primary antibody and Anti-rabbit IgG (H+L), F(ab') Fragment (Alexa Fluor488 Conjugate) (CST, Cat# 4412) was used as the secondary antibody.Fixation and permeabilization were performed according to the manufacturer protocol.pRb and DNA content (stained with PI) were analyzed by flow cytometry.

Analysis of mVenus-p27K-mouse-derived BM cells
Surface-marker-stained BM mononuclear cells (MNCs) (see 'Flow cytometry and cell sorting' for details) were analyzed by flow cytometry to determine the frequency of G0 marker positivity for EPCR + or EPCR -CD150 + CD48 -LSK or progenitor cells.

Seahorse flux analyzer
The extracellular acidification rate (ECAR) and oxygen consumption rate (OCR) were measured using a Seahorse XFe96 extracellular flux analyzer according to the manufacturer's instructions (Agilent Technologies).Briefly, sorted cells were dispensed to culture plates pre-coated with Cell-Tak (Corning) and then the media was replaced with pre-warmed XF-DMEM medium (Agilent) supplemented with 10 mM glucose, 1 mM pyruvate, and 2 mM glutamine, followed by centrifugation at 200 × g for 5 min.OCR and ECAR were measured at baseline and again after sequential addition of respiratory inhibitors at final concentrations of 1 μM oligomycin (an inhibitor of ATP synthase), 2 µM FCCP (an uncoupling agent of mitochondrial respiration), 0.5 µM rotenone/antimycin (an inhibitor of mitochondrial complex I/III) and 50 mM 2-deoxy-D-glucose (an inhibitor of glycolysis).The experiment was performed by dispensing 75,000 HSCs (PBS or 5-FU treated) or MyPs per well.
For the in vivo 2-NBDG assay, C57BL/6 mice treated with PBS or 5-FU were subjected to an in vivo 2-NBDG assay as reported by Jun et al., 2021.Mice received a bolus dose of 375 µg 2-NBDG intravenously and were euthanized by cervical dislocation after 1 hr.Mice were immediately placed on ice, and all subsequent cell preparation processes were performed while the cells were chilled on ice.The 2-NBDG positive cell fraction was detected by flow cytometry.

Conversion of GO-ATeam2 fluorescence to ATP concentration
The GO-ATeam2 knock-in mice were reported by Yamamoto et al., 2019.Briefly, we used a CAG promoter-based knock-in strategy targeting the Rosa26 locus to generate GO-ATeam2 knock-in mice.A study presenting the significance of measuring the absolute concentration of ATP at the singlecell level is currently in preparation for submission, but briefly, the FRET efficiency was converted to the absolute concentration of ATP using the following method (Watanuki et al., in preparation).To permeabilize BM cells, α-hemolysin stock solution (Sigma-Aldrich, St. Louis, MO, USA) was diluted in permeabilization buffer (140 mM KCl (Wako, Cat# 163-03545), 6 mM NaCl (Wako, Cat# 191-01665), 0.1 mM EGTA (Wako, Cat# QB-6401), and 10 mM HEPES (Wako, Cat# 342-01375) [pH 7.4]) to a final concentration of 50 µg/mL α-hemolysin.GO-ATeam2-knock-in BMMNCs were added to the buffer and permeabilized for 30 min at 37 °C under 5% CO 2 .To calibrate ATP concentration, calibration buffer (140 mM KCl, 6 mM NaCl, 0.5 mM MgCl2 (Wako, Cat# 136-03995), and 10 mM HEPES [pH 7.4]) and Mg-ATP stock solution (Sigma-Aldrich, Cat# A9187) were prepared.After washing GO-ATeam2knock-in BMMNCs with calibration buffer, fresh calibration buffer without ATP was added.Mg-ATP was gradually added to increase ATP concentration in the cell suspension, and FRET values of the GO-ATeam2 biosensor at defined ATP concentrations were analyzed by flow cytometry.The FRET value (relative ratio of FRET to EGFP fluorescence intensities) was calculated by the following equation.
The excitation wavelength of FRET and EGFP was set at 488 nm.
The FRET value was then fitted to Hill's formula (Hill, 1910) as a function of ATP concentration: where θ is the original percentage of receptor proteins occupied by the ligand, [L] is the free (unbound) ligand concentration, K A is the concentration of ligand at half saturation, and n is Hill's coefficient.Equation 2 was transformed as such that θ could be expressed by the FRET value as follows: We estimated parameters n and K A by fitting observed FRET values to the linear regression model represented in Equation 3. In our experiment, n=3.1234, and K A = 0.84699.Using these parameters, cellular ATP concentration, [L], was estimated.

Time-course analysis of FRET values
GO-ATeam2 is a ratiometric biosensor that monitors ATP concentration through Förster resonance energy transfer (FRET) from EGFP to the monomeric version of Kusabira Orange (mKO), regardless of the sensor expression levels (Nakano et al., 2011).Surface-marker-stained BMMNCs from GO-ATeam2 mice were dispensed into a basal medium (Ba-M) containing minimal salts, vitamins, and buffers (HEPES and sodium bicarbonate), but no glucose or mitochondrial substrates (Supplementary file 1), or into a medium containing mitochondrial substrates, pyruvate, lactate, fatty acids, and amino acids, but no glucose (PLFA medium).Depending on the experiment, fresh and surface marker-stained BMMNCs obtained from mice 2, 5, or 14 days after intraperitoneal administration of PBS or 5-FU were dispensed into a Ba-M or PLFA medium.The FRET/EGFP ratio data was imported continuously during analysis in a real-time manner using the BD FACSAria IIIu under ambient pressure.Depending on their purpose, experiments were conducted in the presence or absence of various nutrients or metabolic modulators (Figure 4-figure supplement 1D).For this platform, 2 mL of Ba-M or PLFA medium per tube was pre-saturated with 1% O 2 /5% CO 2 /94% N 2 to stabilize ATP levels of BMMNCs (Figure 4-figure supplement 1E-F) and mimic the hypoxic BM environment; when medium was not pre-saturated, ATP concentrations rapidly decreased, even in the presence of glucose, pyruvate, or lactate (Figure 4-figure supplement 1E-F).
To reduce the effect of autofluorescence as much as possible, the top 40-50% of EGFP and FRET fractions of MFI were used in the analysis (MFI >1000 for EGFP and FRET).Then, data reporting EGFP and FRET fluorescence values in individual cells from each gating (e.g.HSCs, MPPs) were extracted along with time course data.Relevant nutrients and inhibitors were added to medium with samples for analysis.Data acquired by the FACSAria IIIu device and retrieved as FCS files were analyzed by the flowCore package in R software.The FRET/EGFP ratio of each set of single cells was fitted to a generalized additive model using the 'gam' function in the 'mgcv' package with 's', a spline-based smoothing function, in default settings as a function of time, then smoothened using the 'predict' function.Pseudocolor plots of the FRET/EGFP ratio were created using the 'kde2d function'.If needed, fitted data were converted to ATP concentration using the model described above.
To compare changes in ATP concentrations in PBS-and 5-FU-treated groups, we corrected differences in baseline ATP concentrations by multiplying all data from the PBS-treated group by the following value: ATP concentration at 0 s in the 5-FU group/ATP concentration at 0 s in the PBS group.

FAOBlue assay
Surface marker-stained BMMNCs from PBS-or 5-FU-treated mice were dispensed at 3×10 5 cells in 500 µL Ba-M, which had been pre-saturated for 48 hr under 1% O 2 and 5% CO 2 conditions and contained 200 mg/dL glucose and 50 µM verapamil.These cells were then exposed to 5 µM FAOBlue (Funakoshi) for 15 min.As a negative control, BMMNCs were exposed to 100 µM etomoxir simultaneously with FAOBlue.The FAOBlue-stained BMMNCs were then centrifuged at 340 × g for 5 min at 4 °C, and the supernatant was discarded.The fluorescence of FAOBlue was excited at a wavelength of 405 nm and detected in the V 450/50 channel.After analysis, the HSC fraction data were extracted.
Comparison of metabolite levels before and after sorting c-Kit-positive cells were isolated using Auto-MACS Pro (Miltenyi Biotec) with the Possel-s or Possel-d2 program as described above (see 'Cell preparation' for details).Isolated cells were counted, and 1×10 5 viable cells were dispensed into methanol containing an internal standard as a pre-sorting cell sample and stored at −80 °C until IC-MS analysis.To the isolated cell suspension, 0.1% PI was added and samples were sorted using the FACS Aria IIIu.A total of 1×10 5 viable cells (PI -cells) were sorted directly into methanol containing an internal standard as a post-sorting cell sample and stored at −80 °C until IC-MS analysis.The detected metabolites were quantified based on calibration curve data (see 'Ion chromatography mass spectrometry (IC-MS) analysis' for details).
Preparation and storage of in vitro U-13 C 6 -glucose tracer samples For tracer analysis, C57BL/6 mice were euthanized to obtain 25,000-50,000 cells per sample of each fraction (HSC, MPP, GMP, CLP) from BM using the FACSAria IIIu instrument.Numbers of mice used to obtain each fraction were as follows: 30-35 each for steady state HSCs and MPPs, 60-65 each for 5-FU treated HSCs, 10 each for GMPs and CLPs.In addition, bone and BM cells were chilled by placing dishes and tubes on ice during the cell preparation process; samples were washed with icecold buffer throughout the entire process before cell sorting.Experiments and experimental manipulations regarding the sampling of mouse femurs and tibias were also performed in the shortest amount of time possible by skilled personnel.Cells were sorted in 0.1% BSA +PBS and sorted cells were centrifuged at 340 × g and 4 °C for 5 min.After discarding the supernatant, cells were added to 1 mL pre-saturated (under 1% O 2 and 5% CO 2 ) GO-ATeam2 Ba-M +0.1% BSA+200 mg/dL U-13 C 6 -(Sigma-Aldrich, Cat# 389374) or U-12 C 6 -glucose and incubated 10 or 30 min.If the process of presaturation was omitted, ATP levels dropped rapidly within a short time (Figure 4-figure supplement 1F).When using oligomycin (1 µM Cell Signaling Technology, Cat# 9996), exposure time was set to 10 min.Samples were then immediately centrifuged at 1000 × g and 4 °C for 3 min.After discarding supernatants, cells were frozen and stored at −80 °C.
Preparation and storage of in vivo U-13 C 6 -glucose tracer samples U-13 C 6 -glucose administration to C57BL/6 mice was performed based on the methods of Jun et al., 2021, with some modifications.Mice were intraperitoneally administered medetomidine hydrochloride, midazolam, and butorphanol tartrate at 0.75 mg/kg, 4 mg/kg, and 5 mg/kg, respectively.After anesthesia, mice were kept warm on a hot plate set at 37 °C while a 27-gauge needle was placed in the external tail vein and U-13 C 6 -glucose was continuously administered.The dose and duration of U-13 C 6 -glucose administration followed (Jun et al., 2021), and 0.4125 mg/g body mass was administered in 1 min, followed by 0.008 mg/g body mass per minute for 3 hr.After U-13 C 6 -glucose administration, mice were euthanized by cervical dislocation and immediately placed on ice.For in vivo tracer analysis, BMMNCs from the bilateral femur, tibia, pelvis, and sternum of each mouse were used to prepare sufficient numbers of HSCs, and pre-chilled 0.1% BSA +PBS was used for BM flushing and washing.HSCs were directly sorted in methanol and stored at −80 °C until IC-MS analysis.A total of 1-3×10 4 HSCs were purified from one or two mice in the PBS group and from two or three mice in the 5-FU group.
When generating the heat map of labeling rates in each metabolite, 1 was added as a pseudo number to the labeling rate of all metabolites.When calculating the total amount of 13 C labeled metabolites for each pathway, metabolites other than M+0 were summed in each metabolite.

Ion chromatography mass spectrometry (IC-MS) analysis
For metabolome analysis focused on glycolytic metabolites and nucleotides, anionic metabolites were measured using an orbitrap-type MS (Q-Exactive Focus; Thermo Fisher Scientific, Waltham, MA, USA) connected to a high-performance IC system (ICS-5000+, Thermo Fisher Scientific), enabling highly selective and sensitive metabolite quantification owing to the IC-separation and Fourier Transfer MS principle (Miyajima et al., 2017).The IC instrument was equipped with an anion electrolytic suppressor (Dionex AERS 500; Thermo Fisher Scientific) to convert the potassium hydroxide gradient into pure water before the sample entered the mass spectrometer.Separation was performed using a Dionex IonPac AS11-HC-4 μm IC column (Thermo Fisher Scientific).The IC flow rate was 0.25 mL/min supplemented post-column with a 0.18 mL/min makeup flow of MeOH.The potassium hydroxide gradient conditions for IC separation were as follows: 1-100 mM (0-40 min), 100 mM (40-50 min), and 1 mM (50.1-60 min), with a column temperature of 30 °C.The Q-Exactive Focus mass spectrometer was operated under the ESI negative mode for all detections.A full mass scan (m/z 70-900) was performed at a resolution of 70,000.The automatic gain control target was set at 3×10 6 ions, and the maximum ion injection time was 100ms.Source ionization parameters were optimized with a spray voltage of 3 kV, and other parameters were as follows: transfer temperature, 320 °C; S-lens level, 50; heater temperature, 300 °C; sheath gas, 36; and aux gas, 10.Metabolite amounts were quantified from calibration curve data generated based on peak areas and respective metabolite amounts.
Quantitative 13 C-MFA with OpenMebius OpenMebius (Open source software for 13 C-MFA) provides the platform to simulate isotope labeling enrichment from a user-defined metabolic model setup worksheet developed in MATLAB (Math-Works, Natick, MA, USA; Kajihata et al., 2014).Quantitative 13 C-MFA was performed according to a manual prepared by the software developer (http://www-shimizu.ist.osaka-u.ac.jp/hp/en/software/OpenMebius.html), but some metabolic model modifications were made to more faithfully reflect our measured data.Specifically, the model was modified to include (a) the conversion of pyruvate to lactate catalyzed by lactate dehydrogenase, (b) the formation of citrate from acetyl CoA and oxaloacetate catalyzed by citrate synthase, (c) the synthesis of alpha-ketoglutarate from citrate catalyzed by aconitase and isocitrate dehydrogenase, and (d) the synthesis of fumarate from succinate by succinate dehydrogenase.Reactions with pyruvate formate lyase performed by Escherichia coli, Streptococcus spp., and ethanol fermentation of acetyl CoA were excluded from the default metabolic network sheet.
The lactate efflux values in 13 C-MFA were determined using the following trial and error method.First, various values (0-100) were entered as candidate lactate efflux values and simulations were run to determine the optimal lactate efflux.When the lactate efflux value was set low (below 50), either the simulation could not be run and an error occurred, or the simulation resulted in the glycolytic system progressing in the opposite direction.These results suggested that the appropriate solution was not obtained because the lactate efflux was unnatural compared to the level of glycolytic metabolites.This was validated by experimental data showing that isotopic labeling rates for most glycolytic metabolites were close to 100% at short labeling times (Figure 1-figure supplement 2C).Therefore, we ran the simulation with a higher lactate efflux value.Finally, we set the lactate efflux to 65, which yielded reasonably satisfactory results for nearly 100% labeling of glycolytic and PPP metabolites in PBS-or DMSO-treated HSCs.
The rate of lactate efflux 5-FU-treated HSCs with the rate of glucose uptake set to 100 was defined using the following equation, with the flux in stationary phase HSC set to 65: 65×(Percentage of glycolytic metabolites labeled with 13 C in the total 13 C-labelled metabolites [5-FU-treated HSCs])/(Percentage of glycolytic metabolites labeled with 13 C in the total 13 C-labelled metabolites [PBS-treated HSC]) In the metabolic flux measurements of HSCs under mitochondrial stress, the lactate efflux determined by the above method exceeded the maximum value that could be modeled (85>), so we decreased the lactate efflux flux by 5 and adopted the maximum value, 80, at which modeling became possible.For values of efflux other than those of lactate efflux flux, the values specified by the Open-Mebius manual were used to eliminate arbitrary factors as much as possible.
The metabolic substrate used for labeling was set to 100% U-13 C 6 glucose.Metabolites used in the analysis included the first intermediate metabolite produced when U-13 C 6 glucose is metabolized (e.g.G6P or F6P with all carbons labeled, the labeled metabolite of the first cycle of the TCA cycle) and the unlabeled metabolite that was measured.Some of the labeled metabolites in the TCA cycle (e.g.citrate [M2]) and erythrose 4-phosphate (M4) in PPP were detected with non-negligible amounts of natural isotopes (>5% even when labeled with U-12 C 6 glucose compared to U-13 C 6 glucose).The presence of such natural isotopes may result in overestimation of the amount of increased labeling with U-13 C 6 glucose.In such cases, the amount of natural isotope detected when labeled with U-12 C 6 glucose was subtracted from the amount of labeled metabolite detected with U-13 C 6 glucose.If the resulting true labeled isotope abundance was negative, the labeled amount was modeled as zero.When analyzing in MATLAB, the number of modeling cycles was set to 100, and the iteration time was set to a maximum of 2000 cycles.

Luminometric ATP measurement
HSCs were sorted from C57BL/6 mice treated with PBS or 5-FU and dispensed into pre-saturated GO-ATeam2 medium with 0.1% BSA in a 1%O 2 /5%CO 2 incubator.HSCs were then exposed to 15 µM of PFKFB3 inhibitor (AZ PFKFB3 26) or DMSO and placed in a 1%O 2 /5%CO 2 incubator for 10 min.Cells were centrifuged at 4 °C and 340 × g and the supernatant was removed.ATP measurements were performed according to manufacturer instructions using Cell ATP Assay Reagent Ver. 2 (Toyo B-Net Corporation).The amount of ATP per cell was calculated by dividing the amount of ATP detected by the number of cells used for analysis.

CRISPR/Cas9 knockout (KO) of Pfkfb3
Target sequences of single guide RNA (sgRNA) were provided in a previous report (Chu et al., 2016) and identified using the web tool GenScript (https://www.genscript.com) for Pfkfb3.sgRNAs were synthesized using a CUGA7 gRNA Synthesis Kit (Nippon Gene, Tokyo, Japan, Cat# 314-08691) following manufacturer instructions, diluted to 1.5 µg/µL, and cryopreserved at −80 °C until use.CD150 + CD48 -Flt3 -LSK cells sorted by FACSAria IIIu were cultured in SF-O3 medium supplemented with stem cell factor (SCF) (50 ng/mL) (Peprotech, Cat# 250-03) and thrombopoietin (TPO) (Peprotech, Cat# 300-18) (50 ng/mL) (S50T50 medium) and incubated under 20% O 2 /5% CO 2 conditions for 16-24 hr, enabling subsequent HSC-specific gene editing with the CRISPR-Cas9 system.Ribonucleoprotein complex preparation and electroporation were conducted as previously reported (Gundry et al., 2016).Briefly, 3 µg Cas9 protein (TrueCut Cas9 Protein v2, Thermo Fisher Scientific, Cat# A36496) plus 3 µg of sgRNA were incubated in Buffer T (Invitrogen, Cat# MPK10096) for 20 min at RT in a volume 6 µL.Cultured cells were resuspended in 30 µL Buffer T and added to ribonucleoprotein at a total volume of 36 µL.Cells were electroporated using the Neon Transfection System (Thermo Fisher Scientific) at 1700 V for 20ms with one pulse.The cell suspension was transferred to S50T50 medium and cultured under 20% O 2 /5% CO 2 conditions.To evaluate gene editing efficiency, genomic DNA from LSK cells was extracted using the NucleoSpin system (Macherey-Nagel, Dürin, Germany) 2-3 d after electroporation.PCR was performed using the following settings: 95 °C for 2 min; 35 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s; followed by final extension at 72 °C for 5 min.PCR products were purified using Wizard SV Gel and the PCR Clean-Up System (Promega Corporation, Madison, WI, USA, Cat# A9281) following manufacturer instructions.A tracking of indels by decomposition (TIDE) assay (Brinkman et al., 2014) or inference of CRISPR edits analysis (Conant et al., 2022) was performed to analyze the sequence data of each PCR product obtained by Sanger sequencing.Among five sgRNAs, Pfkfb3-sg1 displayed the best editing efficiency and was used for subsequent transplant and culture experiments.

BM transplant of Pfkfb3-KO HSCs
Either Rosa26 (control) or Pfkfb3 sequences in HSCs were targeted using CRISPR/Cas9.After electroporation, HSCs were incubated for 2-3 hr in S50T50 medium under 5%CO 2 /20%O 2 conditions, and then counted using a TC10 Automated Cell Counter (Bio-Rad Laboratories, Inc, Hercules, CA, USA).Subsequently, 500 gene-edited HSCs together with 2×10 6 BM cells from Ly5.1 congenic mice were transplanted retro-orbitally into lethally (9.5 Gy using MBR-1520R with a 125 kV 10 mA, 0.5 mm Al, 0.2 mm Cu filter)-irradiated Ly5.1 mice.During Pfkfb3 KO using the vector-free CRISPR-Cas9 system, the KO efficiency was not 100%, so the transplanted cells were a mixture of Pfkfb3-KO cells and wildtype cells.Therefore, after 2, 8, and 16 weeks, peripheral blood was collected and donor-derived chimerism was assessed by a TIDE assay based on a recent study by Shiroshita et al., 2022.The following oligonucleotides for sgRNA synthesis and primers for post-knockout genomic PCR were used.
The sequences of triple-gRNA and the primer set used to confirm KO efficiency were as follows.

Cell cycle analysis and apoptosis assay of Pfkfb3-KO/overexpressing
HSPCs after non-competitive BMT BMMNCs were collected from the bilateral femur, tibia, pelvic bone, and sternum of each individual recipient mouse on day 2 after noncompetitive BMT.Recipient BMMNCs were then stained with Lineage-marker-PerCP-Cy5.5, Ly5.1-PerCP-Cy5.5, and Ly5.2-PE (cell cycle analysis) or Lineage-marker-FITC, Ly5.1-FITC, and Ly5.2-Alexa Fluor700 (apoptosis assay).For the analysis, all BMMNCs from each recipient were used in one analysis, and all lineage-marker negative Ly5.2 + cells were analyzed.Cell cycle analysis (Ki67/Hoechst33432 staining) was performed as described in the 'Ki67/Hoechst33432 staining' section.In vivo BrdU labeling assays were performed as reported by Jun et al., 2021 using the FITC BrdU Flow Kit (BD Biosciences, Cat# 559619).Apoptosis assays were performed using the PE Annexin V Apoptosis Detection Kit I according to manufacturer instructions.
Cell cycle analysis (Ki67/Hoechst33432 staining) of Pfkfb3-overexpressing HSPCs after transplantation was also performed using all BMMNCs from each recipient mouse, and the analysis was performed on all Pfkfb3-overexpressing cells (GFP + ).

5-FU administration after BM recovery in Pfkfb3-KO HSPCs
PFKFB3 was gene-edited in HSPCs using triple-gRNA as described above, and the equivalent of 3×10 5 LSK cells were transplanted retro-orbitally into lethally-irradiated (8.5 Gy using MBR-1520R-3) recipient Ly5.1 mice noncompetitively.After 2 months, recipient mice were treated with 150 mg/kg of 5-FU intraperitoneally.Peripheral blood was collected on the day of 5-FU administration (day 1), and on days 4, 6, 9, and 16.The dynamics of Pfkfb3-or Rosa26-KO cell abundance (as control group) were analyzed by Sanger sequencing as described above.
Homing assay of Pfkfb3-KO HSPCs PFKFB3 was gene-edited in GFP + HSPCs using triple-gRNA as described above.After editing, 2×10 5 cells were retro-orbitally transplanted into lethally-irradiated (8.5 Gy) C57BL/6 mice.After 16 hours, BMMNCs from recipients were stained for surface antigens and analyzed for the percentage of GFP + cells within the PI-negative cells.

Immunocytochemistry
HSCs from PBS-or 5-FU-treated C57BL/6 mice were subjected to immunocytochemistry using antibodies for PFKFB3 (Abcam, Cat# ab181861), phosphorylated-PFKFB3 (Bioss, Cat# bs-3331R), and methylated-PFKFB3 (developed by Takehiro Yamamoto) (Yamamoto et al., 2014).Purified HSCs were resuspended in 50% FCS-PBS and cytospun using the Thermo Scientific Cytospin 4 system (Thermo Fisher Scientific).When using 2-NBDG-positive or -negative HSCs, C57BL/6 mice were given 2-NBDG intravenously (see 'In vivo 2-NBDG assay' for details) and subjected to cytospinning.Cytospun cells were fixed using 4% paraformaldehyde in PBS pH 7.4 for 10 min at RT.Fixed cells were washed twice with ice-cold PBS.For permeabilization, cells were incubated for 5 min with PBS containing 0.1% Triton X-100.Permeabilized cells were washed once with ice-cold PBS.After blocking with 3% BSA-PBS for 30 min, cells were incubated in the diluted antibody with 0.3% BSA-PBS in a humidified chamber overnight at 4 °C.A dilution factor of 1:100 was used for all antibodies.The next day, cells were incubated with Goat anti-Mouse IgG2a Secondary Antibody, Alexa Fluor 555 (Thermo Fisher Scientific, Cat# A-21137) and DAPI in 0.3% BSA-PBS for 1 hr at RT.After two washes with ice-cold PBS, samples were coverslipped with a drop of mounting medium and imaged with a Zeiss LSM 880 microscope (ZEISS, Jena, Germany).Images were acquired at room temperature under darkened conditions using a 100 x oil immersion lens.The obtained image data was analyzed using Imaris software (Bitplane) to calculate the MFI of the target for each cell.

RNA sequencing
Library preparation for RNA-seq was performed on 3000-3500 HSCs derived from mice after 5-FU or PBS administration.Total RNA was prepared using Rneasy Micro kit (QIAGEN, Hilden, Germany).cDNA was synthesized and amplified using SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Takara Bio, Inc, Shiga, Japan).RNA-seq libraries were prepared using the Nextera XT Kit (Illumina, San Diego, CA, USA).Single-end 75 bp sequencing was performed on a NextSeq 500 platform (Illumina).RNA-seq data were obtained from three independent experiments (biological replicates) for each cell type.TopHat (version 2.0.13; with default parameters) was used for mapping to the reference genome (UCSC/mm10) with annotation data from iGenomes (Illumina).Then, gene expression levels were quantified using Cuffdiff (Cufflinks version 2.2.1; with default parameters).
Megakaryocytes were identified as cells with high forward scatter and side scatter, as well as high CD150 and CD41 expression.

Statistical analysis
Data are presented as means ± SD unless otherwise stated.For multiple comparisons, statistical significance was determined by Tukey's multiple comparison test using the Tukey HSD function in the R×64 4.0.3software (R Core Team, Vienna, Austria).A paired or unpaired two-tailed Student's t-test and two-way ANOVA with Sidak's test were used for experiments with two groups.A p-value < 0.05 was considered statistically significant.
The following dataset was generated:

Figure supplement 1 .
Figure supplement 1. Dependence on glycolysis increases with cell cycle progression of HSCs.

Figure supplement 3 .
Figure supplement 3. Mito Stress test results using the Seahorse flux analyzer.

Figure 3 .
Figure 3. Quantitative 13 C-MFA of quiescent, proliferative, and stressed HSCs.(A-C) Overview of quantitative 13 C-MFA of PBS-treated HSCs (A), 5-FUtreated HSCs (B), and OXPHOS-inhibited HSCs (C).The representative net flux for each reaction with glucose uptake as 100 is shown in the squares below the catalytic enzymes for each reaction listed in green letters.Red arrows indicate reactions with particularly elevated fluxes and blue arrows indicate reactions with particularly decreased fluxes.(D) Heatmap of the relative flux of each enzyme in the 5-FU or oligomycin groups compared to Figure 3 continued on next page

Figure 3 -
figure supplements 1-2.The online version of this article includes the following source data and figure supplement(s) for figure 3: Source data 1.Raw data for Figure 3D-L.
Figure3 continued Figure 5B-C; Figure 5-figure supplement 1A-D).Treatment with 2-DG decreased ATP concentrations for a short amount of time (~12 min) in HSCs and HPCs, but ATP reduction was less evident than that induced by oligomycin (blue lines in Figure 5B-C; Figure 5-figure supplement 1A-D).

Figure 4 .
Figure 4. PFKFB3 activates the glycolytic system in proliferating HSCs.(A) Experimental design used to conduct real-time ATP analysis of HSCs treated with 5-FU or PBS.PLFA medium containing mitochondrial substrates (pyruvate, lactate, fatty acids, and amino acids) but no glucose, was used for experiments with 2-DG; Ba-M containing neither mitochondrial substrates nor glucose was used for experiments with oligomycin, PFKFB3 inhibitor, or AMPK inhibitor.(B-E) Results of real-time ATP analysis of PBS-(red) or 5-FU-treated (blue) HSCs after treatment with 2-DG (B, D), oligomycin (C, E). (F) Normalized mRNA counts of PFKFB isozymes based on the RNA sequencing of HSCs.(G-J) Results of real-time ATP analysis of PBS-(red) or 5-FU-treated (blue) HSCs after treatment with PFKFB3 inhibitor (G, I), or AMPK inhibitor (H, J).Bar graphs show corrected ATP concentrations for the last 2 min (D) of (B), 6-7 min (E) of (C), or the last 1 min (I, J) of (G, H) for PFKFB3 and AMPK inhibitors, respectively.Each group represents at least 60 cells.Data are representative results of pooled samples from three biological replicates.(see 'Time-course analysis of FRET values' in 'Materials and methods' for details of the correction method used to calculate ATP concentration.)Data are presented as mean ± SD. * p≤0.05, ** p≤0.01, *** p≤0.001 as determined by Student's t-test (D, E, I, and J) or a one-way ANOVA followed by Tukey's test (F).See also Figure 4-figure supplements 1-2.The online version of this article includes the following source data and figure supplement(s) for figure 4: Source data 1.Raw data for Figure 4D-F, I and J.

Figure supplement 1 .
Figure supplement 1. Establishment of a real-time ATP concentration analysis system using GO-ATeam2.

Figure supplement 1
Figure supplement 1-source data 1.Raw data for Figure 4-figure supplement 1B, C and I.

Figure supplement 2 .
Figure supplement 2. FAO is not active in proliferating HSCs.

Figure supplement 2
Figure supplement 2-source data 1.Raw data for Figure 4-figure supplement 2C, D, G-I, K, M, N and P.

Figure 5 .
Figure 4 continued stem/progenitor fractions in PLFA medium (red bars) alone or PLFA medium plus 2-DG (blue bars).ATP concentration for the last 2 min of the analysis time is shown.Data is summarized from (B, C) and Figure 5-figure supplement 1.Each group represents at least 110 cells.Data are representative results of pooled samples from three biological replicates.(E) ATP concentration in indicated stem/progenitor fractions in Ba-M plus glucose (dark blue bars) or Ba-M plus glucose and oligomycin (orange bars).ATP concentration for the last 1 min of the analysis period is shown.Data is summarized from (B, C) and Figure 5-figure supplement 1.Each group represents at least 43 cells.Data are representative results of pooled samples from three biological replicates.(F-I) Effects of PFKFB3 or AMPK inhibitors (PFKFB3i or AMPKi, respectively) on ATP concentration in HSCs from GO-ATeam2 mice in Ba-M plus glucose only (F) or Ba-M plus glucose and oligomycin (G).ATP concentrations for the last 1 min of the analysis period are shown in (H) and (I) for glucose only and glucose with oligomycin groups, respectively.Each group represents at least 90 cells.Data are representative results of pooled samples from three biological replicates.(J) Experimental schema for cell cycle assay and real-time ATP concentration analysis after overexpression of Pfkfb3.(K) Cell cycle status of Pfkfb3-overexpressing (Pfkfb3OE) and mock-transduced HSCs.(L-M) Effects of inhibitors on ATP concentration in Pfkfb3overexpressing GO-ATeam2 + HSCs.Cells were exposed to vehicle or 2-DG (L), oligomycin in the presence or absence of glucose 12.5 mg/dL (M), and ATP concentrations for the last 2 min (L) or 1 min (M) of the analysis period were calculated.Data are representative results of pooled samples from three biological replicates.Data are presented as mean ± SD. * p≤0.05, ** p≤0.01, *** p≤0.001 as determined by Student's t-test (D, E, and K) or oneway ANOVA followed by Tukey's test (H, I, L, and M).See also Figure 5-figure supplement 1.The online version of this article includes the following source data and figure supplement(s) for figure 5: Source data 1.Raw data for Figure 5D, E, H, I, K, L and M.

Figure supplement 1 .
Figure supplement 1. Steady-state PFKFB3 activity defines HSC and HPC metabolic kinetics and cell cycle.

Figure
Figure supplement 1-source data 1.Raw data for Figure 5-figure supplement 1F,J-L, Q.

Figure 6 .
Figure 5 continued in each group.n=26-27 single HSCs for each group.The data are representative results from two independent experiments.(C) Quantification of MFI of phosphorylated-PFKFB3 (p-PFKFB3) protein in PBS-or 5-FU-treated HSCs.The lower part of the graph shows representative images of immunocytochemistry of p-PFKFB3 in each group.n=27 single HSCs for each group.The data are representative results from two independent experiments.(D) Quantification of MFI of p-PFKFB3 in HSCs treated with glucose (200 mg/dL); glucose plus oligomycin (1 µM); and glucose, oligomycin, and dorsomorphin (100 µM) for 5 min.The lower part of the graph shows representative images of immunocytochemistry of p-PFKFB3 in each group.n=32-36 for each group.The data are representative results from two independent experiments.(E) Normalized Prmt1 mRNA counts based on RNA sequencing of PBS-treated (red) or 5-FU-treated (blue) HSCs.Data are representative results of pooled samples from three biological replicates.(F) MFI quantification of methylated-PFKFB3 (m-PFKFB3) in PBS-or 5-FU-treated HSCs.The lower part of the graph shows representative images of immunocytochemistry of m-PFKFB3 in each group.n=23-41 for each group.The data are representative results from three independent experiments.(G) Quantification of MFI of m-PFKFB3 in PBS-or 5-FU-treated HSCs or 5-FU-treated HSCs after 15 min treatment with a PRMT1 inhibitor (90 μg/ mL GSK3368715); n=25-35 single HSCs for each group.The lower part of the graph shows representative images showing immunocytochemistry of m-PFKFB3.Data represent a single experiment.(H) Quantitation of m-PFKFB3 in NBDG-positive or -negative HSCs in mice treated with PBS or 5-FU.The lower part of the graph shows representative images of immunocytochemistry of m-PFKFB3 in each group.n=28-41 for each group.The data are representative results from two independent experiments.(I) Corrected ATP levels in PBS-(red) or 5-FU-treated (blue) HSCs 15 min after treatment with vehicle or a PRMT1 inhibitor (90 µg/mL GSK3368715).Each group represents at least 101 cells.Data are representative results of pooled samples of two biological replicates.(see 'Time-course analysis of FRET values' in 'Materials and methods' for details of the correction method used to calculate ATP concentration.)(J) ATP concentration in mock-transduced (Ctl) or Pfkfb3-overexpressed (OE) HSCs after treatment with the PRMT1 inhibitor (90 µg/mL GSK3368715).ATP concentration for the last 1 min of the analysis period is shown.Data are presented as mean ± SD. * p≤0.05, ** p≤0.01, *** p≤0.001 as determined by Student's t-test (A-C, E-F, and I-J) or one-way ANOVA followed by Tukey's test (D, G, and H).See also Figure 6-figure supplement 1.The online version of this article includes the following source data and figure supplement(s) for figure 6: Source data 1.Raw data for Figure 6A-J.

Figure supplement 1 .
Figure supplement 1. Gene expression related to Prmt1 in proliferating HSCs.

Figure
Figure supplement 1-source data 1.Raw data for Figure 6-figure supplement 1.

Figure 7 -
figure supplement 1.The online version of this article includes the following source data and figure supplement(s) for figure 7: Source data 1.Raw data for Figure 7B, C, E, I-K and N.

Figure supplement 1 .
Figure supplement 1. PFKFB3 contributes to HSC proliferation and differentiation in vitro.
Figure 7 continued