6-Phosphogluconate dehydrogenase (6PGD), a key checkpoint in reprogramming of regulatory T cells metabolism and function
Abstract
Cellular metabolism has key roles in T cells differentiation and function. CD4+ T helper-1 (Th1), Th2, and Th17 subsets are highly glycolytic while regulatory T cells (Tregs) use glucose during expansion but rely on fatty acid oxidation for function. Upon uptake, glucose can enter pentose phosphate pathway (PPP) or be used in glycolysis. Here, we showed that blocking 6-phosphogluconate dehydrogenase (6PGD) in the oxidative PPP resulted in substantial reduction of Tregs suppressive function and shifts toward Th1, Th2, and Th17 phenotypes which led to the development of fetal inflammatory disorder in mice model. These in turn improved anti-tumor responses and worsened the outcomes of colitis model. Metabolically, 6PGD blocked Tregs showed improved glycolysis and enhanced non-oxidative PPP to support nucleotide biosynthesis. These results uncover critical role of 6PGD in modulating Tregs plasticity and function, which qualifies it as a novel metabolic checkpoint for immunotherapy applications.
Introduction
Regulatory T cells (Tregs) are distinct subsets of CD4+ T cells that limit inflammatory responses by suppressing self-reactive T cells to maintain balance between protective and excessive inflammation. Although lipid oxidation and the mevalonate pathway being considered key to their functions (Michalek et al., 2011; Timilshina et al., 2019), cellular metabolism in Tregs is an active area of investigation. Tregs depend on these pathways to maintain their suppressive functions, which are distinct from the high glycolytic signatures of conventional CD4+ effector subsets, Th1, Th2, and Th17 (Michalek et al., 2011). Although glycolysis is not major feature in Tregs, it is critical for Treg expansion (Pacella et al., 2018). Glucose (Glc) enters the cell via glucose transporter 1 (Glut1) and is metabolized to glucose 6-phosphate (G6P), which has alternative fates via glycolysis and as the entry point into the pentose phosphate pathway (PPP) (Figure 1—figure supplement 1).
The PPP, also known as the hexose monophosphate shunt or phosphogluconate pathway, comprises two interdependent branches: the oxidative branch and the non-oxidative branch. Oxidative PPP is composed of three essentially irreversible reactions in which the first is catalyzed by G6PD that converts G6P to 6-phosphogluconolactone (6PGL) while generating reduced nicotinamide adenine dinucleotide phosphate (NADPH). 6PGL then hydrolyzes to 6-phosphogluconate (6PG), which is oxidatively decarboxylated by 6-phosphogluconate dehydrogenase (6PGD), producing ribulose-5-phosphate and a second molecule of NADPH and (Figure 1—figure supplement 1A). Ribulose-5-phosphate is the precursor for ribose-5-phophate (R5P) and thence nucleotides, and NADPH is required both for reductive anabolism and for antioxidative defense.
Tregs are known to have high rates of lipid oxidation (Michalek et al., 2011), which would provoke antioxidant machinery which requires NADPH generated from PPP. We thus investigated the role of PPP in T CD4+ subsets. Our preliminary studies using inhibitors of the PPP enzymes indicated that the third enzyme in the oxidative PPP, 6PGD, has critical roles in modulating T cell functions.
To assess the role of 6PGD in modulating Treg function, we generated mice with flox targeted Pgd allele (Pgdfl/fl), which were crossed with Foxp3YPF-Cre mice to generate Pgdfl/flFoxp3YPF-Cre (thereafter Pgdfl/flFoxp3Cre). Interestingly, Pgdfl/flFoxp3Cre mice showed a distinct phenotype, namely spontaneous fatal autoimmunity and lower lifespan. Isolated Pgd-deficient Tregs had diminished suppressive functions both in vitro and in vivo. They also displayed shift in differentiation toward Th1, Th2, and Th17 CD4+ T cell subtypes. Stable isotope tracing studies on sorted Tregs revealed reprogrammed metabolic pathways consistent with the generation of these subsets. Our results for Tregs directly isolated from pro-inflammatory mice model and derived from naïve CD4+ T cells in vitro uncover a critical role of 6PGD in maintaining Tregs metabolism, phenotype, and function.
Results
PGD deficiency in Tregs induces spontaneous fatal autoimmunity
Our preliminary studies showed that pharmacological inhibitor of 6PGD (6-aminonicotinamide [6-AN]) but not G6PD inhibitor (dehydroepiandrosterone [DHEA]) enhances IFN-γ production by in vitro activated CD4+ T cells (Figure 1—figure supplement 1B), while DHEA but not 6-AN at 10 µM instigated poor T cells viability (Figure 1—figure supplement 1C). Based on these, we generated Pgdfl/fl mice (Figure 1—figure supplement 1D) and crossed with Cd4Cre mice (Figure 1—figure supplement 1F). An enhanced effector phenotype of T CD4+ subsets was observed in Pgdfl/flCd4Cre mice (Figure 1—figure supplement 1I-J) suggesting altered regulation of T cell differentiation and/or activation.
To specifically delete 6PGD in Tregs, Pgdfl/fl mice were crossed with Foxp3Cre mice (Figure 1—figure supplement 1G). In these mice, the Cre recombinase driven by the Foxp3 promoter will delete Pgd only in Tregs upon Foxp3 expression and YFP can be used as expression marker (Rubtsov et al., 2008). Pgd deletion in YFP expressing cells was verified by western blot, real-time PCR (Figure 1A–B, Figure 1—figure supplement 1E), and stable isotope-resolved metabolomics (SIRM) analysis . Notably, disruption of 6PGD in Tregs led to a profound inflammatory disorder while Pgdfl/flCd4Cre mice appeared normal (Figure 1—figure supplement 1K-L). Compared with Foxp3Cre (wild-type [WT]) controls at 20 days of age, Pgdfl/flFoxp3Cre mice exhibited a reduced body size and crusting of ears and eyelids (Figure 1C), enlargement of peripheral lymphoid organs (Figure 1D), and lymphoproliferation (Figure 1E) that was accompanied by a significantly shorter lifespan (Figure 1F). Lymphoproliferation was attributed to a higher frequency of CD4+ and CD8+ T cells in spleen (Figure 1G) with increased frequencies of CD44highCD62Llow effector/memory phenotype (Figure 1H), reduced frequencies of CD44lowCD62Lhigh naïve cells (Figure 1H), higher expressions of CD69 activation marker (Figure 1I), and enhanced production of IFN-γ (Figure 1J) in both CD4+ and CD8+ T cells. In these mice, CD8+ T cells also produced more granzyme-B and higher cytotoxic marker CD107a expression (Figure 1K). The titers of serum antibodies were much higher in Pgdfl/flFoxp3Cre than WT mice (Figure 1L). Moreover, Pgdfl/flFoxp3Cre mice had a significant increase in serum levels of IFN-γ, IL-17A, IL-4, and IL-5 (Figure 1M). Along with splenomegaly and lymphadenopathy in Pgdfl/flFoxp3Cre mice, there was an enhanced immune cell infiltration and tissue damage in main organs (Figure 1N). The development of observed fatal inflammatory disease in Pgdfl/flFoxp3Cre mice suggests a crucial role of 6PGD in modulating the function of Tregs. Such phenotype has also been reported in association with central signaling elements of metabolic reprogramming such as mechanistic target of rapamycin (mTOR) (Zeng et al., 2013) and liver kinase B1 (Yang et al., 2017).
Deletion of 6PGD in Tregs results in loss of their suppressive function
Based on the observed phenotype in Pgdfl/fl Foxp3Cre mice, we further evaluated specific phenotypic markers and suppressive functions of Tregs, which were either isolated from inflamed mice model or generated in vitro from naïve CD4+ T cells. The CD4+Foxp3+ cells in the spleen of 20-day-old Pgdfl/flFoxp3Cre mice were lower than the WT counterparts, both in number (Figure 2A) and in frequency of total CD4+ T cells (Figure 2B). The expression of CD98 (a subunit of the neutral amino acid transporter), CD71 (the transferrin receptor), and CD40L in these 6PGD-deficient Tregs was enhanced (Figure 2C). CD71 and CD98 are key nutrient receptors in Tregs that depend on mTORC1 (Kelly et al., 2007).
Next, YFP+ (indicator of Foxp3 expression) cells (Rubtsov et al., 2008) were sorted from 20-day-old WT (Foxp3YFP-Cre) and Pgdfl/flFoxp3Cre mice and evaluated for phenotypic and functional properties. YFP+ sorted cells are demonstrated in Figure 2D. The level of Foxp3 transcript was lower in 6PGD-deficient cells evaluated by real-time PCR (Figure 2E), which is consistent with lower Foxp3 gMFI in these cells (Figure 2B). When culturing sorted YFP+ Tregs in vitro in the presence of IL-2 (700 IU/ml) and anti-CD3/anti-CD28 coated beads (Treg:beads ratio 1:3), higher number of cells were detected for 6PGD-deficient cells at both 24 and 48 hr time points (Figure 2F). This result showed that 6PGD deficiency was tolerable in Tregs and did not induce cell death in these cells. Evaluation of the culture media demonstrated significantly higher production of cytokines as markers for other CD4+ T helper subsets such as IFN-γ (Th1), IL-13 and IL-5 (Th2), and IL-17A (Th17) due to 6PGD blockade (Figure 2G). RNA sequencing (RNAseq) analysis of isolated Tregs from inflamed mice also demonstrated higher expression levels of Th1 (IFNg, Tbx21), Th2 (Il13, Il5, Il4 and Gata3), and Th17 (Il17a) marker genes in 6PGD-deficient Tregs, shown as heat map in Figure 2—figure supplement 1A and volcano plot in Figure 2—figure supplement 1B. RNAseq data were confirmed by real-time PCR for selected genes (Figure 2—figure supplement 1C) and flow cytometry (Figure 2—figure supplement 1D). It should be noted that although these shifts in gene expression are consequences of 6PGD blockade in Tregs, the resulting severe autoimmunity in mice with pro-inflammatory microenvironment should also be considered. Pro-inflammatory condition of mice could enhance or diminish some of these gene expression changes.
Although Foxp3 expression is a key marker of the Treg formation, additional metabolic and transcriptional regulations can orchestrate Treg plasticity (Shi and Chi, 2019). Pgd-deficient Tregs can be converted to Th1, Th2, and Th17 by modification of metabolic pathways. Notably, the Th2 markers IL13 and IL5 were among the most upregulated genes in Pgd deleted Tregs (Figure 2—figure supplement 1A-B), which points to Treg differentiation into Th2 type cells. These gene expression changes were consistent with significantly reduced suppressive capacity of Pgd-deficient versus WT Treg (Figure 2H–I), which was assayed by co-culturing the corresponding isolated YFP+ cells with CD4+CD45RBhigh effector cells in vitro, as described in supplemental information. Altered Treg function was along with enhancement of CD4+ T cells phenotype shift toward Th1, Th2, and Th17 cells subsets in Pgdfl/flFoxp3Cre mice (Figure 2—figure supplement 1E-G). Pro-inflammatory condition of mice could also contribute in the magnitude of these outcomes.
We also examined Tregs suppressive function in vivo using an inflammatory bowel disease (IBD) model (Figure 2J). YFP+ cells from WT and Pgdfl/flFoxp3Cre mice as well as WT CD4+ T effector (CD4+CD45RBhigh) cells were isolated, mixed, and transferred to Rag1-/- mice (Ostanin et al., 2009; Steinbach et al., 2015). Higher rate of weight loss was observed in Pgdfl/flFoxp3Cre group (Figure 2K) compared to WT. Colons were evaluated on day 45 post cells injection in the IBD model mice, 6PGD deficiency in Tregs led to shorter and thicker colons compared with no injection or injection with WT Treg (Figure 2L). Hematoxylin and eosin staining of colons demonstrated higher immune cell infiltration and destruction of colon structure (Figure 2M). Taken together these results confirm that 6PGD inhibition is associated with loss of Treg suppressive function, both in vitro and in vivo.
6PGD deficiency attenuates Treg suppressive function and induces potent anti-tumor responses
To confirm key observations from YFP+ sorted cells of Pgdfl/flFoxp3Cre mice that displayed inflammatory condition, we characterized changes in the molecular properties of WT Tregs induced by the small molecule inhibitor of 6PGD, 6-AN (Davis and Kauffman, 1987). We also examined the effect of 6PGD blockade in Tregs on tumor formation by tamoxifen treatment of Pgdfl/fl Foxp3EGFP-Cre-ERT2 mice (Figure 1—figure supplement 1H and Figure 3).
First, we used WT Foxp3Cre naïve (CD62Lhigh CD44low) CD4+ T cells to generate Tregs in vitro (as inducible Tregs [iTregs]) in the polarization medium (IL-2+ TGF-β) plus 6-AN or vehicle DMSO. This medium effectively induced Foxp3 expression indicated by the YFP marker in both 6-AN and DMSO treatments (Figure 3A), with lower YFP gMFI in 6-AN treated cell (Figure 3B). Tregs generated under 6-AN treatment in vitro showed a lower level of Foxp3 (Figure 3C), but a higher level of Gata3 and IL5 (Figure 3C) expression. This difference in Foxp3/Gata3 expression in generated Tregs in vitro indicates altered Treg properties in 6PGD blocked cells. Tregs generated in the presence of 6-AN also showed diminished suppressive function in vitro (Figure 3D–E).
To control the timing of Pgd deletion in Tregs and evaluate Tregs suppressive function under controlled inflammation, we generated Pgdfl/flFoxp3EGFP-Cre-ERT2 mice by crossing Pgdfl/fl mice with Foxp3EGFP-Cre-ERT2 mice (Figure 1—figure supplement 1H). In these mice, expression of Cre in Foxp3 promoter is induced by tamoxifen treatment and EGFP can be used as a detection marker (Rubtsov et al., 2010). Suppressive function of targeted Tregs was evaluated in vitro and in vivo in a tumor model. In 6-week-old Pgdfl/flFoxp3EGFP-Cre-ERT2 mice, tamoxifen was injected on days 0, 1, and 3 to induce deletion of 6PGD in Tregs. To induce tumor formation, B16F10 melanoma cells were injected into the flank region of the mice on day 10 post tamoxifen treatment (Figure 3F). Although the tamoxifen induction does not induce deletion in 100% of Tregs in vivo, this model would induce sufficient 6PGD deletion in targeted Tregs to be tracked in vivo (Rubtsov et al., 2010). Foxp3 expression was detected by EGFP expression (CD4+CD25+EGFP+ cells in spleen) (Figure 3G). EGFP+ sorted cells showed reduced suppressive activity in vitro (Figure 3H), as seen in the YFP+ sorted cells from Pgdfl/fl Foxp3Cre mice (Figure 2H–I) or Tregs generated in presence of 6-AN (Figure 3D–E).
In the B16F10 tumor model, induced 6PGD deletion in Tregs resulted in reduced tumor growth (Figure 3I) and smaller final mass (Figure 3J) on day 16 post tumor induction. Assessing cell types in the tumor, we saw more CD8+ and CD4+ T cell infiltration (Figure 3K–L), along with higher IFN-γ production from both CD4+ (Figure 3M) and CD8+ (Figure 3N) T cells in Pgdfl/flFoxp3EGFP-Cre-ERT2 versus WT mice. Thus, in the tumor microenvironment, Tregs were altered by 6PGD blockade to unleash host anti-tumor responses. As such, 6PGD can serve as a metabolic checkpoint for immune activation by modulating Tregs activity, which is considered to be desirable (Wang et al., 2017; Pacella and Piconese, 2019).
6PGD function in Tregs is required for control of allergic (Th2) responses
Tregs are responsible for suppression of T subset effector responses but they also have the plasticity to convert into these subsets (Shi and Chi, 2019). A well-defined organ to follow Th2 responses control by Tregs is lung. In Pgdfl/flFoxp3Cre mice, significant Th2 (allergic) responses were evident in the lung (Figure 3—figure supplement 1). These included higher collagen deposits (Figure 3—figure supplement 1A, top panels), higher mucus accumulation (Figure 3—figure supplement 1A, bottom panels), and higher infiltration of T cells (both CD4+ and CD8+ cells) in the lung (Figure 3—figure supplement 1C) with enhanced production of Th2 markers such as IL-13 (Figure 3—figure supplement 1D), IL-5 (Figure 3—figure supplement 1F), and IL-4 (Figure 3—figure supplement 1H) by CD4+ T cells. The sorted YFP+ Tregs from 6PGD blocked mice also expressed lower levels of Foxp3 and higher levels of Gata3 transcription factors (Figure 3—figure supplement 1B). Moreover, the Th2 markers were higher in the spleen of the Pgdfl/flFoxp3Cre mice (Figure 3—figure supplement 1E,G,I).
High allergic responses in Pgdfl/fl Foxp3Cre mice was also evident with higher infiltration of eosinophils (CD11clow/- Siglect-F+) in the lung (Figure 3—figure supplement 1J; Stevens et al., 2007) and spleen (CD11b+ Siglect-F+) (Figure 3—figure supplement 1K; Hey et al., 2015). These allergic responses were accompanied by infiltration of CD11b+ Ly6G+ cells in the lung and spleen of Pgdfl/flFoxp3Cre mice (Figure 3—figure supplement 1L-M; Xue et al., 2020; Bronte et al., 2016).
Although the shift toward Th2 responses (Figures: 1M; Figure 2—figure supplement 1B-C; Figure 3—figure supplement 1; Figure 2G; and Figure 3C) was dominant in 6PGD-deficient Tregs, enhancement of Th1 and Th17 responses was also noted by higher IFN-γ and IL-17A in the serum (Figure 1M) and higher frequency of CD4+IFN-γ+ and CD4+IL-17A+ cells in the spleen (Figure 2—figure supplement 1E, G). Beside these enhanced Th1 and Th17 responses, YFP+ Tregs isolated from Pgdfl/fl Foxp3Cre mice showed elevated expression of Ifng, Tbx21, and Il17a genes (Figure 2—figure supplement 1A, C) and sorted Tregs cultured in vitro secreted higher amount of IFN-γ and IL-17A into media (Figure 2G).
In other reports, comparable patterns of Treg plasticity have been observed by blocking a key regulator of Tregs metabolism, mTOR. Mtorfl/flFoxp3Cre mice showed 10- to 15-fold increases of Th2 mediators while increases in Th1 and Th17 responses were 2- to 3-fold (Chapman et al., 2018).
6PGD blockade in Tregs induces metabolic reprogramming similar to that of T effector cells
Next, we examined metabolic shifts in the Tregs isolated from WT and Pgdfl/flFoxp3Cre mice. It should be noted that the pro-inflammatory microenvironment of the mice could contribute to Tregs metabolic reprogramming, beside the intrinsic effects of 6PGD. We first assessed glycolysis and mitochondrial respiration by measuring the extracellular acidification rate (ECAR) and oxygen consumption rate (OCR) (Figure 4A–B). Tregs isolated from Pgdfl/flFoxp3Cre mice showed significant increase in both basal and reserve capacity in ECAR (Figure 4A) and OCR (Figure 4B) compared to WT Tregs. We also evaluated Glc uptake and mitochondrial status by fluorescent dye staining (Figure 4C). Compared to WT Tregs, 6PGD-deficient Tregs had increased Glc uptake capacity determined by 2-NBDG uptake (Figure 4C). These cells also showed increase in mitochondrial potential (ΔΨm) as measured by TMRE (Figure 4C) without changes in the mitochondrial mass as examined by MitoTracker Deep Red FM (Figure 4C). Moreover, mitochondrial ROS (assessed by MitoSox Red) was enhanced in 6PGD blocked Tregs (Figure 4C). These results are consistent with the observed skew of the Tregs phenotype and markers toward conventional CD4+ T effector subsets (Michalek et al., 2011; Pacella et al., 2018).
During transition from naïve CD4+ T cells, Tregs switch from Glc utilization to other fuel sources such as fatty acids (Procaccini et al., 2016; Macintyre et al., 2014), which could help sustain their function in low Glc conditions such as in the tumor microenvironment. Glutamine (Gln) is an important fuel source for T cells that is consumed in mitochondria. However, it has been shown that deficiency of the glutamine transporter ASCT2 did not impact Tregs generation (Nakaya et al., 2014), which raises the question of the importance of Gln metabolism in Tregs.
Based on these observations, we tracked simultaneously Glc and Gln metabolism in WT and Pgdfl/flFoxp3Cre Tregs, both as direct isolates from mouse models with inflammatory disorder and also Tregs generated from naïve CD4+ T cells in vitro. In the first set, YFP+ Tregs were sorted from WT and Pgdfl/flFoxp3Cre mice and cultured in the presence of IL-2 (700 IU/ml) and anti-CD3/anti-CD28 coated beads (Treg:beads ratio 1:3) plus D7-Glc and 13C5,15N2-Gln for 48 hr. In the second set, naïve CD4+ from WT and Pgdfl/flFoxp3Cre mice were polarized with IL-2+ TGF-β for 4 days and the medium was changed to that containing D7-Glc + 13C5,15N2-Gln for the last 48 hr. Transformations of D7-Glc and 13C515N2-Gln in sorted Tregs into metabolites of PPP/gluconeogenesis (GNG) (Figure 4D), glycolysis (Figure 4—figure supplement 1A), Krebs cycle (Figure 4—figure supplement 1B), and the nucleotide synthesis pathway (Figure 4—figure supplement 1C) were determined by SIRM analysis. Shift in the same metabolic pathways examined in in vitro generated Tregs is presented in Figure 4—figure supplement 2 and Figure 4—figure supplement 3.
SIRM analysis showed higher buildup of D-labeled glycolytic intermediates (Figure 4—figure supplement 1A), which together with enhanced ECAR (Figure 4A) pointed to increased glycolytic capacity in 6PGD blocked Tregs. These cells also displayed increased levels of D and/or 13C-labeled Krebs cycle metabolites (Figure 4—figure supplement 1B), which together with the observed increased OCR and mitochondrial potential (Figure 4B–C) indicated elevated Krebs cycle activity. Notably, Gln served as the main fuel source for the Krebs cycle (Figure 4—figure supplement 1B), with Gln metabolism being significantly higher in 6PGD-deficient than WT Tregs. Enhanced Gln metabolism via the Krebs cycle led to increased αKG production (Figure 4—figure supplement 1B-F), which could serve as critical ‘metabolic signal’ for CD4+ T subset differentiation (Klysz et al., 2015). αKG level has been shown to control methylation of the Foxp3 gene locus resulting in a shift of balance between Treg/Th1 (Klysz et al., 2015) and Treg/Th17 (Xu et al., 2017). Whether such epigenetic metabolic mechanism is involved in skewing 6PGD blocked Tregs toward Teff subsets awaits further investigation.
6PGD deficiency in Tregs enhances non-oxidative PPP and GNG
SIRM analysis also showed that 6PGD-deficient Tregs extensively accumulated G6P and 6PG (Figure 4D–a and b), which confirmed effective blockade of oxidative PPP at 6PGD (Figure 4—figure supplement 1A). However, instead of depletion, this blockade led to enhanced accumulation of D and/or 13C-labeled R5P (Figure 4D–c), which are downstream products of the 6PGD reaction. Also D and/or 13C-labeled sedoheptulose-7-phosphate (S7P) (Figure 4D–d), erythrose-4-phosphate (Figure 4D–e), and fructose-6-phosphate (F6P) (Figure 4D–f) accumulated in 6PGD-deficient Tregs. These metabolites are all products of the non-oxidative PPP. These results point to enhanced non-oxidative PPP that results in elevated production of R5P from D7-Glc and 13C5,15N2-Gln (Figure 1—figure supplement 1A), despite the blockade at 6PGD. In turn, R5P generated fueled the enhanced biosynthesis of both pyrimidine and purine nucleotides (D, 13C, and 15N labeled species, Figure 4—figure supplement 1C) to fulfill the nucleotide demands of activated T cells.
It should be noted that enhanced incorporation of Gln derived carbon (C* or D*C*/C*NxDx scrambled species) into 6PG (C* and D*C* scrambled species, Figure 4D–b) and ATP/GTP (C*, D*C*, and C*NxDx, Figure 4—figure supplement 1C-m,l) in 6PGD-deficient Tregs are consistent with elevated GNG in these cells. This is because these 13C-labeled species can only be produced via the reaction sequence of glutaminolysis, the Krebs cycle, and GNG initiated by conversion of OAA to PEP via PEPCK, followed by the reversal of the glycolytic reactions. The latter produce 3-phosphoglycerate, glyceraldehyde-3-phosphate, F6P, and G6P, which can fuel both oxidative and non-oxidative PPP as well as serine biosynthesis-one carbon pathway to power purine biosynthesis (Figure 4—figure supplement 1A). Increased accumulation of 13C-G6P, together with that of D-G6P (Figure 4D–a) indicate enhanced activity of G6PD, the first rate-limiting and NADPH-producing enzyme of the oxidative branch of the PPP. This is expected to help compensate for the loss of NADPH production due to blocked 6PGD action. Overall, similar shifts of metabolic patterns in Tregs generated in vitro (iTregs in Figure 4—figure supplement 2 and Figure 4—figure supplement 3) as those sorted from mouse models were evident.
Discussion
In our pilot studies, we found that inhibition of G6PD led to reduced T cell survival with little changes in IFN-γ production (Figure 1—figure supplement 1B-C) but that of 6PGD resulted in higher T cell survival and T subsets activations (Figure 1—figure supplement 1B-C,I-J ). Our G6PD inhibitor study is consistent with a recent study, which showed that G6PD inhibition depleted NADPH and decreased inflammatory cytokine production by T cells (Ghergurovich et al., 2020). To better understand the role of 6PGD in T cells, we generated Pgdfl/flFoxp3Cre mice to examine the effect of 6GPD blockade on Tregs as central immune regulators. Interestingly, Pgdfl/flFoxp3Cre mice showed severe autoimmune disorders highlighting a critical role of 6PGD in Tregs function (Figure 1). Tregs isolated from Pgdfl/flFoxp3Cre mice showed reduced suppressive function in vitro and in vivo (Figure 2), which was accompanied by a shift of gene expression signatures and phenotype markers toward T CD4+ subsets of Th1, Th2, and Th17. Loss of suppressive function and plasticity toward conventional subsets was also reproduced by treatments with 6PGD small molecule inhibitor 6-AN (Figure 3).
Although Glc, but not Gln, contributes to Tregs differentiation and function (Pacella et al., 2018; Procaccini et al., 2016; Nakaya et al., 2014), Glc depletion by blocking the Glc transporter Glut1 did not alter the suppressive function of Tregs (Michalek et al., 2011; Macintyre et al., 2014). However, our study showed that the suppressive function of Tregs was modulated by altered Glc metabolism via the oxidative PPP, that is, 6PGD deficiency in Tregs resulting in significantly improved anti-tumor responses (Figure 3F–N) or changed IBD outcomes (Figure 2). This qualifies 6PGD as an immune checkpoint by modulating the Treg function while implicating the importance of Glc metabolism via the oxidative PPP in Treg activation.
As our knowledge of Glc and Gln metabolism in Tregs is limited, we tracked the metabolic fates of both fuel substrates with stable isotope tracers in 6PGD blocked versus WT Tregs. We found that blocking 6PGD in oxidative PPP shifted Glc flow into glycolysis and non-oxidative PPP to provide sufficient substrates (e.g. one-carbon metabolites and R5P) for fueling nucleotide biosynthesis (Figure 4 and Figure 4—figure supplement 1). This was also accompanied by enhanced Gln flow into GNG to power nucleobase synthesis and NADPH production, which can compensate for NADPH loss due to 6PGD blockade. These reprogrammed metabolic activities co-occurred with the shift of expression in phenotypic and functional markers toward conventional T CD4+ subsets. The importance of glycolytic shift on modulating Tregs fitness has been shown in previous studies (Watson et al., 2021; Zappasodi et al., 2021). Tregs in high Glc environment lose their fitness while upregulation of lactate metabolism is important for maintaining their suppressive phenotype (Watson et al., 2021). Also, the presence of Glc was shown to be a key factor in interfering Treg function under anti-CTLA-4 treatment (Zappasodi et al., 2021). Moreover, although glycolysis is key to supporting the energy and anabolic demands during Treg expansion (Pacella et al., 2018; De Rosa et al., 2015), further increase in glycolysis as a result of mTOR and c-Myc activation suppresses Treg generation (Kabat et al., 2016; Wei et al., 2016). Our metabolic fate analysis not only corroborates with glycolytic control of Treg functions but also reveals the role of 6PGD as a key regulator of Treg stability.
Although disruption of Gln uptake did not impact Treg function (Nakaya et al., 2014), enhanced Gln metabolism in 6PGD-deficient Tregs was evident and co-occurred with the shift from Treg to Th1 (Michalek et al., 2011; Klysz et al., 2015), Th2 (Michalek et al., 2011), or Th17 (Michalek et al., 2011; Xu et al., 2017) phenotypes. Increased Gln oxidation in the Krebs cycle promoted the synthesis of immunomodulatory metabolites such as αKG (Figure 4—figure supplement 1), which has been shown to signal Foxp3 suppression (Klysz et al., 2015; Xu et al., 2017). At the same time, diversion of Gln carbons for enhanced nucleotide synthesis would help meet the metabolic demand of the shift toward conventional T subsets. Thus, Gln metabolism could synergize with 6PGD deficiency in eliciting Treg plasticity toward immune activation. Further studies will be required to substantiate such metabolic mechanism.
In summary, we reported a critical role of 6PGD, the second rate-limiting enzyme of PPP, in skewing Tregs function and metabolism toward those of the Th1, Th2, and Th17 subsets. As such, 6PGD could be a novel metabolic target in clinical applications such as autoimmune disorders or cancer therapies.
Materials and methods
Cell line
Request a detailed protocolB16-F10 cell line was obtained from ATCC (ATCC CRL-6475). The cell line was tested for Mycoplasma contamination.
Mice
Steps to generate of Pgdfl/flCd4Cre, Pgdfl/flFoxp3Cre, and Pgdfl/fl Foxp3EGFP-Cre-ERT2 mice are described in Figure 1—figure supplement 1. Verification of 6PGD deletion was done by PCR, western blot, and alteration in metabolite profile associated with the enzyme. Other mice used during experiments including C57BL/6 J (B6 CD45.2+), B6.SJL-Ptprca Pepcb/BoyJ (B6 CD45.1+), B6.129(Cg)-Foxp3tm4(YFP/icre)Ayr/J (Foxp3YFP-cre), Foxp3tm9(EGFP/cre/ERT2)Ayr/J (Foxp3eGFP-Cre-ERT2), and B6.129S7-Rag1tm1Mom/J (Rag1-/-) mice were purchased from Jackson Laboratories. All mice were kept in specific pathogen-free conditions prior to use. Animal work was done in accordance with the Institutional Animal Care and Use Committee of the University of BIDMC (Protocol number 040–2016).
Detection of cells markers by flow cytometry
Request a detailed protocolFlow cytometry stainings were done on prepared cells. Cells from spleen and lymph nodes were obtained by tissue harvest, mechanically disintegration, passing through 70 μm cell strainer and then red blood cell lysis. Tumor infiltration cells were also prepared using the same steps with additional 30 min collagenase treatment at 37°C. Cultured cells and sorted cells just passed one step of wash with FACS buffer (PBS containing 1% FBS). Cells were suspended in FACS buffer and stained for surface markers for 30 min at 4°C. Aqua LIVE/DEAD Fixable Dead Cell Stain was used in all staining to distinguish dead cells. Intracellular staining for IL-4, IL-13, Il-5, IFN-γ, IL-17A, TGF-β, and granzyme B was performed with the Fixation/Permeabilization Solution Kit with BD GolgiPlug (BD Bioscience) according to the manufacturer’s protocol. To detect cytokines, cells were cultured in RMPI containing 10% FBS and were stimulated with PMA (50 ng/ml), ionomycin (1 μg/ml), and GolgiPlug (1 μl/ml) for 4 hr at 37C before starting the staining. For assessment of Foxp3, cells were stained intracellularly using the Foxp3/Transcription Factor Buffer Set (eBiosciences) following the manufacturer’s protocol. Glc uptake was analyzed by incubation of the cells with 20 μM 2-NBD-Glc (Cayman Chemical) for 30 min at 37°C. Cells were analyzed by a CytoFLEX LX (Beckman Coulter) cytometer and analysis was performed using FlowJo_V10 software. All antibodies are listed in the Key resources table.
Cell isolation by flow cytometry cell sorting
Request a detailed protocolTo obtain purified cells, single cells were first enriched by negative selection using EasySep Mouse CD4+ T cells isolation kit (STEMCELL Technologies) following manufacturer’s protocols. Then cells were stained for desired surface markers and sorted by flow cytometry. Surface markers for gating/sorting of each cell type included: naïve CD4+ T cells (CD4+ CD62Lhigh CD44low), YFP+ Tregs (CD4+ CD25+ YFP+), EFGP+ Tregs (CD4+ CD25+ EGFP+), and T effector cells (CD4+ CD45RBhigh). All the stainings were done at 4°C for 30 min and sorting was done by Becton Dickinson SORP FACSAria II (Beckman Coulter). Isolated cell phenotypes were confirmed by running again with flow cytometer (CytoFLEX LX, Beckman Coulter).
Cell culture and treatment
Naïve CD4+ T cell treatment with inhibitors
Request a detailed protocolNaïve CD4+ T cells (CD4+ CD62Lhigh CD44low) were isolated and cultured in complete RPMI media containing 10% FBS and IL-2 (20 IU/ml) and treated with 6-AN (10 µM), DHEA (10 µM), or vehicle DMSO. Cells stimulation was done by adding soluble anti-CD3 (1 µg/ml) and soluble anti-CD28 (1 µg/ml). Cells were harvested and evaluated 4 days post culture.
Naïve CD4+ T cells polarization toward Tregs and treatment with 6-AN
Request a detailed protocolTo target 6PGD in vitro, naïve CD4+ T cells (CD4+ CD62Lhigh CD44low) were isolated and polarized to Tregs using IL-2 (50 IU/ml) and TGF-β1 (5 ng/ml). Stimulation was done by culturing the cells in plate coated antibodies. For coating anti-CD3 and anti-CD28 (2 µg/ml PBS each) were put on the plate and incubated 4 hr at 37°C and then washed twice with PBS. For inhibitor treatment 6-AN (10 µM) or vehicle DMSO was added.
Culture of YFP+ sorted cells
Request a detailed protocolIsolated Tregs (CD4+ CD25+ YFP+) were cultured in vitro in presence of IL-2 (700 IU/ml) and anti-CD3/anti-CD28 coated beads (Treg:beads ratio 1:3) and cells number was assessed at 24 and 48 hr time points. Media was evaluated for cytokines production.
Suppression assay
Request a detailed protocolSuppression assay was done on prepared Tregs as YFP+ sorted cells from 6PGD+/+Foxp3Cre (WT) and Pgdfl/flFoxp3Cre mice, EGFP+ sorted cells from Pgd+/+ Foxp3EGFP-Cre-ERT2 (WT) and Pgdfl/fl Foxp3EGFP-Cre-ERT2 mice, and also DMSO (control) and 6-AN derived Tregs. These cells express CD45.2 on their surface (B6 CD45.2+). To check suppressive capacity of prepared Tregs, the cells were mixed with WT CD45.1+ effector T cells (CD4+ CD45RBhigh) in different ratio and cultured in U shaped plate. The effecter cells were stained with CFSE (0.5 mM) for 30 min before culture. For optimized stimulation of effector cells, 2.5 × 105 irradiated (30 Gy) splenocytes from CD45.2+ mice (WT B6 CD45.2+) and soluble anti-CD3 (1 µg/ml) was added. After 72 hr of co-culture, amount of CFSE dilution was determined in CD4+ CD45.1+ effector cells using flow cytometer (CytoFLEX LX, Beckman Coulter).
RNA purification, RNAseq, and real-time quantitative PCR
Request a detailed protocolTotal RNA was extracted from YFP+ sorted cells with RNeasy Mini Kit (Qiagen) and RNA concentrations were determined using Nanodrop (Thermo Scientific). Total DNA-free RNA was used for mRNA isolation and library construction. Libraries were sequenced on an Illumina HISEQ 2500 (Illumina). Confirmation of RNAseq was done by examination of key representative genes by real-time quantitative PCR method using ABI 7300 Real-Time PCR system (Applied Biosystems). Used TaqMan probes are listed in the Key resources table. Expression of each gene was normalized to the housekeeping gene expression (18 S rRNA) as ΔCT = CT (gene X) – CT (18 S rRNA). The alteration between groups was calculated as fold change = 2^(ΔΔCT). ΔΔCT is: CT (gene X) – CT (gene X at baseline).
Bioenergetics analysis by Seahorse and mitochondrial activity
Request a detailed protocolPrepared Tregs (YFP+ sorted cells) were seeded on Cell-Tak coated Seahorse XFe24 (Agilent) culture plates (0.8–1 × 106 cells/well) in assay media. The media was DMEM supplemented with 1% BSA and 25 mM Glc, 1 mM pyruvate, and 2 mM Gln. The ECAR and OCR were measured over 75 min. To obtain maximal respiratory and control values, cells were treated with oligomycin (1 μM), FCCP (1.5 μM), and rotenone/antimycin A (0.5 μM).
Mitochondrial activity also was examined by flow cytometry. For this, cells were stained with 200 nM MitoTracker Deep Red FM (ThermoFisher) for mitochondrial mass and 200 nM tetra-methylrhodamine ester (TMRE) (ThermoFisher) for mitochondrial membrane potential. Staining was done in FACS buffer (PBS containing 2% FBS) along with cell surface markers for 30 min at 37°C. To detect mitochondrial reactive oxygen species cells surface markers were added and then 5 μM MitoSOX Red (ThermoFisher) was added in last 10 min incubation at 37°C. Cells were analyzed by a CytoFLEX LX (Beckman Coulter) cytometer and analyzed by FlowJo_V10.
IBD model
Request a detailed protocolTo examine Tregs suppressive activities in vivo, the IBD model was used as described previously (Ostanin et al., 2009; Steinbach et al., 2015). In brief, YFP+ cells from WT and Pgdfl/fl Foxp3Cre mice and also CD4+ effector T cells (CD4+CD45RBhigh) from WT mice were sorted and mixed. Cells were resuspended to 4 × 106 cells/ml effector (CD45high) and 2 × 106 cells/ml Treg (YFP+). Total of 4 × 105 CD45high and 2 × 105 Treg was injected into each Rag1-/- (B6.129S7-Rag1tm1Mom/J) mice intra-peritoneally (i.p.). Control positive group was injection of effector cells alone and negative control group was Rag1-/- mice without injection. Disease progress was monitored by weight loss, change in stools form, and diarrhea. Colons length and thickness along with colon’s histochemistry were evaluated 45 days post cells injection.
SIRM analysis on Tregs
Request a detailed protocolTo trace Glc and Gln in Tregs, YFP+ sorted cells from WT and Pgdfl/flFoxp3Cre mice were cultured in presence of IL-2 (700 IU/ml) and anti-CD3/anti-CD28 coated beads (Treg:beads ratio 1:3) and D7-Glc plus 13C5,15N2-Gln for 48 hr. To prepare labeling media, Glc-free Gln-free RPMI medium was supplemented with 10% dialyzed FBS (Life Technologies), 20 mM HEPES, 0.05 mM 2-mercaptoethanol, and 1% penicillin-streptomycin. Isotope labeling patterns of metabolites of cell extracts were analyzed by Ion Chromatography-UltraHigh Resolution Mass Spectrometry. The procedures were followed as previously described (Fan et al., 2016b; Fan et al., 2016a).
Statistical analysis
Request a detailed protocolAll the statistical analyses were done using GraphPad Prism 6 (GraphPad) and Excel software. Difference between groups was analyzed using the Student’s t test (for two groups) or one-way ANOVA (for multiple groups), followed by the post hoc Tukey test. Differences in tumor growth was calculated by two-way ANOVA. Differences between groups were rated significant at values of p < 0.05. In the figures all data are shown as mean ± SEM and *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.
Data availability
RNAseq data was deposited to NCBI BioSample (BioSample Accession # PRJNA706901).
-
NCBI Sequence Read ArchiveID PRJNA706901. 6-Phosphogluconate dehydrogenase (6PGD) a key checkpoint in reprogramming of Treg metabolism and function.
References
-
A small molecule G6PD inhibitor reveals immune dependence on pentose phosphate pathwayNature Chemical Biology 16:731–739.https://doi.org/10.1038/s41589-020-0533-x
-
Redefining Myeloid Cell Subsets in Murine SpleenFrontiers in Immunology 6:652.https://doi.org/10.3389/fimmu.2015.00652
-
T cell transfer model of chronic colitis: concepts, considerations, and tricks of the tradeAmerican Journal of Physiology. Gastrointestinal and Liver Physiology 296:G135–G146.https://doi.org/10.1152/ajpgi.90462.2008
-
Induction of Murine Intestinal Inflammation by Adoptive Transfer of Effector CD4+ CD45RB high T Cells into Immunodeficient MiceJournal of Visualized Experiments 2015:98.https://doi.org/10.3791/52533
-
Detection and quantitation of eosinophils in the murine respiratory tract by flow cytometryJournal of Immunological Methods 327:63–74.https://doi.org/10.1016/j.jim.2007.07.011
-
Elevated granulocytic myeloid-derived suppressor cells are closely related with elevation of Th17 cells in mice with experimental asthmaTernational Journal of Biological Sciences 16:2072–2083.https://doi.org/10.7150/ijbs.43596
Article and author information
Author details
Funding
Bayer Fund (2014-07-1112)
- Pankaj Seth
NIH Office of the Director (1U24DK097215-01A1)
- Teresa Fan
National Cancer Institute (P30CA177558)
- Teresa Fan
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by BIDMC seed funds to PS, 2014-07-1112 target grant from Bayer to PS, and a pilot grant to PS from 1U24DK097215-01A1 (to RMH and TWMF), and the Markey Cancer Center Redox and Metabolism Shared Resource Facility P30CA177558.
Ethics
Animal work was done in accordance with the Institutional Animal Care and Use Committee of Beth Israel Deaconess Medical Center.(Protocol number 040-2016).
Copyright
© 2021, Daneshmandi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,862
- views
-
- 222
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Physics of Living Systems
The endoplasmic reticulum (ER), the largest cellular compartment, harbours the machinery for the biogenesis of secretory proteins and lipids, calcium storage/mobilisation, and detoxification. It is shaped as layered membranous sheets interconnected with a network of tubules extending throughout the cell. Understanding the influence of the ER morphology dynamics on molecular transport may offer clues to rationalising neuro-pathologies caused by ER morphogen mutations. It remains unclear, however, how the ER facilitates its intra-luminal mobility and homogenises its content. It has been recently proposed that intra-luminal transport may be enabled by active contractions of ER tubules. To surmount the barriers to empirical studies of the minuscule spatial and temporal scales relevant to ER nanofluidics, here we exploit the principles of viscous fluid dynamics to generate a theoretical physical model emulating in silico the content motion in actively contracting nanoscopic tubular networks. The computational model reveals the luminal particle speeds, and their impact in facilitating active transport, of the active contractile behaviour of the different ER components along various time–space parameters. The results of the model indicate that reproducing transport with velocities similar to those reported experimentally in single-particle tracking would require unrealistically high values of tubule contraction site length and rate. Considering further nanofluidic scenarios, we show that width contractions of the ER’s flat domains (perinuclear sheets) generate local flows with only a short-range effect on luminal transport. Only contractions of peripheral sheets can reproduce experimental measurements, provided they are able to contract fast enough.
-
- Cell Biology
- Developmental Biology
How cells regulate the size of their organelles remains a fundamental question in cell biology. Cilia, with their simple structure and surface localization, provide an ideal model for investigating organelle size control. However, most studies on cilia length regulation are primarily performed on several single-celled organisms. In contrast, the mechanism of length regulation in cilia across diverse cell types within multicellular organisms remains a mystery. Similar to humans, zebrafish contain diverse types of cilia with variable lengths. Taking advantage of the transparency of zebrafish embryos, we conducted a comprehensive investigation into intraflagellar transport (IFT), an essential process for ciliogenesis. By generating a transgenic line carrying Ift88-GFP transgene, we observed IFT in multiple types of cilia with varying lengths. Remarkably, cilia exhibited variable IFT speeds in different cell types, with longer cilia exhibiting faster IFT speeds. This increased IFT speed in longer cilia is likely not due to changes in common factors that regulate IFT, such as motor selection, BBSome proteins, or tubulin modification. Interestingly, longer cilia in the ear cristae tend to form larger IFT compared to shorter spinal cord cilia. Reducing the size of IFT particles by knocking down Ift88 slowed IFT speed and resulted in the formation of shorter cilia. Our study proposes an intriguing model of cilia length regulation via controlling IFT speed through the modulation of the size of the IFT complex. This discovery may provide further insights into our understanding of how organelle size is regulated in higher vertebrates.