Zengdi Zhang1, Oscar C Salgado2, Bing Liu1, Zahra Moazzami1,3, Kristin A Hogquist2,4, Michael A Farrar2,4, Hai-Bin Ruan1,2. 1. Department of Integrative Biology and Physiology, University of Minnesota, Minneapolis, MN, United States. 2. Center for Immunology, University of Minnesota, Minneapolis, MN, United States. 3. Department of Food Science and Nutrition, University of Minnesota, Minneapolis, MN, United States. 4. Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, MN, United States.
Abstract
The immunosuppressive regulatory T (Treg) cells exert emerging effects on adipose tissue homeostasis and systemic metabolism. However, the metabolic regulation and effector mechanisms of Treg cells in coping with obesogenic insults are not fully understood. We have previously established an indispensable role of the O-linked N-Acetylglucosamine (O-GlcNAc) signaling in maintaining Treg cell identity and promoting Treg suppressor function, via STAT5 O-GlcNAcylation and activation. Here, we investigate the O-GlcNAc transferase (OGT)-STAT5 axis in driving the immunomodulatory function of Treg cells for metabolic homeostasis. Treg cell-specific OGT deficiency renders mice more vulnerable to high-fat diet (HFD)-induced adiposity and insulin resistance. Conversely, constitutive STAT5 activation in Treg cells confers protection against adipose tissue expansion and impaired glucose and insulin metabolism upon HFD feeding, in part by suppressing adipose lipid uptake and redistributing systemic iron storage. Treg cell function can be augmented by targeting the OGT-STAT5 axis to combat obesity and related metabolic disorders.
The immunosuppressive regulatory T (Treg) cells exert emerging effects on adipose tissue homeostasis and systemic metabolism. However, the metabolic regulation and effector mechanisms of Treg cells in coping with obesogenic insults are not fully understood. We have previously established an indispensable role of the O-linked N-Acetylglucosamine (O-GlcNAc) signaling in maintaining Treg cell identity and promoting Treg suppressor function, via STAT5 O-GlcNAcylation and activation. Here, we investigate the O-GlcNAc transferase (OGT)-STAT5 axis in driving the immunomodulatory function of Treg cells for metabolic homeostasis. Treg cell-specific OGT deficiency renders mice more vulnerable to high-fat diet (HFD)-induced adiposity and insulin resistance. Conversely, constitutive STAT5 activation in Treg cells confers protection against adipose tissue expansion and impaired glucose and insulin metabolism upon HFD feeding, in part by suppressing adipose lipid uptake and redistributing systemic iron storage. Treg cell function can be augmented by targeting the OGT-STAT5 axis to combat obesity and related metabolic disorders.
Foxp3+ regulatory T (Treg) cells are specialized immunosuppressive lymphocytes that control self-tolerance, inflammatory responses, and tissue homeostasis (1, 2). Treg cells develop in the thymus through two distinct programs involving CD25+Foxp3- and CD25-Foxp3lo progenitors (3). The upregulation of CD25 during Treg cell development is driven by strong T cell receptor (TCR) signals (4). The upregulation of Foxp3 is TCR-independent, but instead involves IL-2 signaling through its receptor CD25 and the downstream activation of STAT5 (5–7). Besides of their indispensability during development, both TCR and IL-2/STAT5 signaling pathways are required for the suppressive function of effector Treg cells (8–10). Treg cells are essential in enforcing peripheral tolerance and loss of them leads to a fatal scurfy phenotype that is characterized by severe inflammation in endocrine and barrier organs (11–13).In addition to preserving tissue integrity, Treg cells have been increasingly shown to support adipose tissue function and systemic metabolism. Treg cells, along with other anti-inflammatory immune cells, are abundant in adipose tissues of lean animals (14, 15). The accumulation and function of adipose Treg cells depend on antigen presentation by parenchymal MHC class-II molecules and local cytokines such as IL-33 to drive a PPARγ-dependent transcriptional program (16–18). Obesity profoundly remodels the immune cell compartment of fat tissue, leading to Treg cell reduction and insulin resistance. Adoptive transfer or genetic induction of Treg cells have been shown to suppress adipose inflammation and improve insulin sensitivity in obese animals (19–21). On the other hand, loss of Treg cell function by cell depletion, PPARγ knockout, or Vα5/Vβ8.2 TCR transgene leads to no or marginal changes in adipose insulin resistance (14, 18, 22). It has been proposed that Treg cells in obesity and aging may adopt pro-inflammatory properties (15, 23); therefore, restoring Treg cell function, rather than increasing Treg cell number, may offer therapeutic opportunities for metabolic diseases. Nonetheless, mechanisms mediating the insulin-sensitizing effect of adipose Treg cells are not fully understood. Also, it has not been established if Treg cells can directly modulated lipid metabolism and energy balance.O-linked N-Acetylglucosamine (O-GlcNAc) modification at serine or threonine residues, termed O-GlcNAcylation (24), plays a central role in signaling pathways relevant to chronic human diseases such as diabetes, cancer, and immune disorders (25–27). O-GlcNAc transferase (OGT), using UDP-GlcNAc derived from the hexosamine biosynthetic pathway as the substrate, controls diverse biological processes such as gene transcription, protein stability, and cell signaling (28–31). O-GlcNAcase (OGA) mediates the removal of O-GlcNAcylation from proteins. Protein O-GlcNAcylation senses glucose availability (30, 32), hormonal cues (29, 33, 34), cellular stress (33, 35), and immune and bacterial signals (27, 36, 37). Recently, we found in Treg cells that TCR signaling promotes OGT-mediated O-GlcNAcylation of Foxp3 and STAT5, two transcriptional factors essential for Treg cell development and function (38). Overexpression of a constitutive active form of STAT5 rescues lineage instability and suppressor dysfunction found in O-GlcNAc deficient Treg cells, thus ameliorating systemic autoimmunity in OGT-deficient mice (38). In this study, we seek to interrogate roles of the OGT-STAT5 axis in instructing Treg cells to cope with obesity-associated remodeling of the adipose tissue and systemic metabolism.
Materials and Methods
Mice
Foxp3 mice (stock number 016959) and Foxp3 mice (stock number 016961) were purchased from the Jackson Laboratory. Rosa26 mice were provided by Dr. Alexander Rudensky at Memorial Sloan Kettering Cancer Center (10). Ogt Mice (Jackson Laboratory stock number 004860) were kindly provided by Dr. Xiaoyong Yang at Yale University. All mice were on C57BL/6 background. Mice were free to access water and fed on a regular chow, tamoxifen food (Teklad, TD.130860), or 60% high fat diet (ResearchDiets, D12492) as indicated. All procedures were approved by the Institutional Animal Care and Use Committee at the University of Minnesota. All relevant ethical regulations for animal testing and research were complied with.
Metabolic Assays
Body weights were monitored every week. Body composition was assessed using an EchoMRI system (39). For food intake measurement, mice were individually housed and food consumption was weighed every morning for 7 consecutive days. For the metabolic cage study, mice were first acclimated in metabolic chambers (Columbus Instruments), and then physical activity and energy expenditure were measured continuously for at least 2 days (40). For oral glucose-, and insulin- tolerance tests. 16 h fasted mice were given glucose (1.5 g/kg body weight) by gavage; 6 h fasted mice were injected with insulin (1 U/kg body weight) intraperitoneally. Tail-vein blood collected at the designated times was used to measure blood glucose level using a Contour Glucometer (Bayer). For lipid tolerance test, mice were fasted for 4 h, gavaged with 200 ul of olive oil, and tail-vein blood was collected every 1.5 h for quantification of triglyceride using a colorimetric assay kit (Cayman #10010303). Fecal energy was determined with a bomb calorimeter at Dr. Pedro Urriol’s lab. Fecal lipid was extracted with chloroform: methanol mix (41), weighed and calculated as percentage to fecal weight.
Flow Cytometry
For surface markers, cells were stained in PBS containing 0.5% BSA with relevant antibodies at 4°C for 30 min. For analysis of intracellular markers, cells were first fixed with Fixation/Permeabilization buffer (ThermoFisher, catalog no. 00-5123) at 4°C for 30 min and then stained in Permeabilization Buffer (ThermoFisher, catalog no. 00-8333) with relevant antibodies at 4°C for 30 min. Flow cytometry data were acquired on BD Fortessa H0081 or X-20 and analyzed with Flowjo.
Histology
Mouse tissues were dissected and fixed in 10% buffered formalin. Sectioning and haematoxylin & eosin (H&E) staining were performed by the Comparative Pathology Shared Resource of UMN.
RNA and Real-time PCR
RNA was isolated with Trizol and reverse transcribed into cDNA with the iScript™ cDNA Synthesis Kit. Real-time RT-PCR was performed using iTaq™ Universal SYBR® Green Supermix and gene-specific primers on a Bio-Rad C1000 Thermal Cycler. All data were normalized to the expression of 18s or 36b4.
Western Blotting
Tissues were lysed in buffer containing 1% Nonidet P-40, 50 mM Tris 3 HCl, 0.1 mM EDTA, 150 mM NaCl, proteinase inhibitors and an OGT inhibitor, TMG. Equal amounts of protein lysates were electrophoresed on 4-20% SDS-PAGE gels and transferred to Nitrocellulose membranes. Primary antibodies were incubated at 4°C overnight. Western blotting was visualized by peroxidase conjugated secondary antibodies and ECL chemiluminescent substrate with a Bio-Rad ChemiDoc Imaging System or imaged by fluorescent IgG secondary antibodies with a LI-COR Odyssey 9120 imager.
Statistical Analyses
Results are shown as mean ± SEM. The comparisons were carried out using two-tailed Student’s t test, one-way ANOVA followed by multiple comparisons with the Tukey adjustment, two-way ANOVA using GraphPad Prism 7.
Results
OGT-insufficiency in Treg Cells Leads to Insulin Resistance
Both the Foxp3 and Ogt genes are located on the mouse X Chromosome; about 40 centimorgans apart (
). When we generated Treg cell-specific OGT knockout (KO) using the Foxp3 mice with YFP-Cre knocked into the endogenous Foxp3 locus (
), hemizygous male KO mice developed scurfy phenotypes and died of autoimmune disease by the age of 4 weeks (38), thus preventing us to study their whole-body metabolism. Due to random inactivation of the X chromosome in females, Foxp3 female mice possess both OGT-sufficient (WT) and -deficient (KO) Treg cells (
). Despite lineage instability and suppressor dysfunction of KO Treg cells, the remaining WT cells were sufficient to avert autoimmunity in mice (38). Therefore, Foxp3 females and their littermate controls were used to investigate metabolic outcomes. When fed with normal chow (NC), no difference in body weight was observed in 12-week-old mice (
). Compared to controls, Foxp3 females had higher blood glucose levels during insulin tolerance test, indicative of impaired insulin sensitivity (
). Glucose tolerance test revealed similar blood glucose levels between the two genotypes (
), suggesting that pancreatic insulin secretion might be increased to compensate insulin resistance in KO mice. We then treated mice with insulin to assess its downstream signaling, and found similar AKT phosphorylation in gonadal white adipose tissue (gWAT) and liver (
). However, insulin-stimulated Ser473 phosphorylation of AKT in brown adipose tissue (BAT) was diminished in Foxp3 females (
). These results demonstrate that OGT in Treg cells controls insulin pathway activation in BAT and whole-body insulin sensitivity.
Figure 1
Mosaic OGT-deficiency in Treg cells decreases insulin sensitivity. (A) Chromosome positions of mouse Foxp3 and Ogt genes. The Xist gene is shown for reference. cM, centimorgan. (B) Schematic view of allele combinations of the Foxp3 and Ogt genes in Foxp3 and Foxp3 female mice. Y-C, YFP-Cre; F, floxed; +, wildtype (wt). Red “X” indicates random inactivation of one of the two X chromosomes. (C) Body weight of 3-month-old, NC-fed Foxp3 (n = 18) and Foxp3 (n = 15) female mice. (D) Insulin tolerance test of 3-month-old, NC-fed Foxp3 (n = 17) and Foxp3 (n = 14) female mice. Area under curve (AUC) is shown to the right. (E) pS473-Akt and total-Akt immunoblotting of BAT protein lysates from Foxp3 and Foxp3 female mice treated with saline or insulin for 30 min. Data are presented as mean ± SEM. n.s., not significant; **, p < 0.01 by unpaired student’s t-test.
Mosaic OGT-deficiency in Treg cells decreases insulin sensitivity. (A) Chromosome positions of mouse Foxp3 and Ogt genes. The Xist gene is shown for reference. cM, centimorgan. (B) Schematic view of allele combinations of the Foxp3 and Ogt genes in Foxp3 and Foxp3 female mice. Y-C, YFP-Cre; F, floxed; +, wildtype (wt). Red “X” indicates random inactivation of one of the two X chromosomes. (C) Body weight of 3-month-old, NC-fed Foxp3 (n = 18) and Foxp3 (n = 15) female mice. (D) Insulin tolerance test of 3-month-old, NC-fed Foxp3 (n = 17) and Foxp3 (n = 14) female mice. Area under curve (AUC) is shown to the right. (E) pS473-Akt and total-Akt immunoblotting of BAT protein lysates from Foxp3 and Foxp3 female mice treated with saline or insulin for 30 min. Data are presented as mean ± SEM. n.s., not significant; **, p < 0.01 by unpaired student’s t-test.
OGT-insufficiency in Treg Cells Promotes Diet-induced Obesity
We then fed animals with 60% high fat diet (HFD) to induce obesity. Female Foxp3 mice gained more weight than controls (
). Body composition analysis showed that weigh gain was due to increased fat mass (
). Similar to animals fed with chow, HFD Foxp3 females displayed impaired insulin tolerance (
), but normal glucose tolerance (
).
Figure 2
Foxp3 mice are prone to diet-induced obesity. (A) Growth curve of Foxp3 (n = 19) and Foxp3 (n = 17) female mice after HFD feeding. (B) Lean mass and (C) fat mass of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice. (D) Insulin tolerance test of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice. Area under curve (AUC) is shown to the right. (E) Glucose tolerance test of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 10) female mice. Area under curve (AUC) is shown to the right. (F) Daily HFD intake of Foxp3 (n = 10) and Foxp3 (n = 9) female mice. ANCOVA was used to correct food intake for body weight. (G–J) Metabolic cage study of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice showing locomotor activity (G), respiratory exchange ratio (RER, H), and heat (I). Body weight-adjusted heat analyzed with ANCOVA is shown (J). (K) Cold-induced change in heat production of HFD-fed Foxp3 (n = 5) and Foxp3 (n = 5) female mice. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01 by unpaired student’s t-test (B–G, J) or two-way ANOVA (A, K).
Foxp3 mice are prone to diet-induced obesity. (A) Growth curve of Foxp3 (n = 19) and Foxp3 (n = 17) female mice after HFD feeding. (B) Lean mass and (C) fat mass of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice. (D) Insulin tolerance test of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice. Area under curve (AUC) is shown to the right. (E) Glucose tolerance test of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 10) female mice. Area under curve (AUC) is shown to the right. (F) Daily HFD intake of Foxp3 (n = 10) and Foxp3 (n = 9) female mice. ANCOVA was used to correct food intake for body weight. (G–J) Metabolic cage study of HFD-fed Foxp3 (n = 10) and Foxp3 (n = 9) female mice showing locomotor activity (G), respiratory exchange ratio (RER, H), and heat (I). Body weight-adjusted heat analyzed with ANCOVA is shown (J). (K) Cold-induced change in heat production of HFD-fed Foxp3 (n = 5) and Foxp3 (n = 5) female mice. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01 by unpaired student’s t-test (B–G, J) or two-way ANOVA (A, K).We went on to determine what caused diet-induced adiposity in Foxp3 mice. Their daily intake of HFD was similar to control mice, when adjusted to body weight using ANCOVA analysis (
). Metabolic cage study found no difference in locomotor activity (
) or respiratory exchange ratio (
). Rather, heat production, when adjusted to body weight (42), was significantly reduced in Foxp3 mice during both light and dark periods (
). Moreover, when challenged with cold, Foxp3 animals showed reduced heat production than the controls (
), indicative thermogenic defects. Together, these data demonstrate that OGT-controlled Treg cell function enables metabolic adaptation to positive energy balance.
STAT5 Activation in Treg Cells Confers Protection Against Diet-induced Obesity
The regulation of Treg cell suppressive program by OGT is dependent on STAT5 O-GlcNAcylation and activation (38). We then sought to determine whether STAT5 overactivation, to the opposite of OGT knockout, would drive Treg cell action in systemic metabolism. Foxp3 male mice were generated to specifically overexpress a constitutively active (CA-) form of STAT5B in Treg cells (10). Flow cytometric analyses confirmed that STAT5B-CA not only increased the number of FOXP3+ Treg cells in the lymph nodes and the activity of the Foxp3 gene promoter (
), but also drastically expanded the population of effector Treg cells that were CD44+KLRG1+ (
) or CD44+CD103+ (
). When fed with a normal chow, even though with comparable body weight (
), Foxp3 mice showed a slight, but statistically significant reduction in fat mass (
). There was a declining trend of blood glucose levels during the insulin tolerance test (
). Similar locomotor activity (
), energy expenditure (
), and energy or lipid excretion into the feces (
) were observed between control and Foxp3 mice.
Figure 3
Characterization of Foxp3 mice. (A–C) Flow cytometric quantification of FOXP3-YFP+CD4+ Treg cells (A), CD44+KLRG1+
(B) and CD44+CD103+
(C) effector Treg cells in lymph nodes from Foxp3 and Foxp3 mice. (D) Growth curve of Foxp3 (n = 10) and Foxp3 (n = 9) mice fed with NC. (E) Lean and fat mass of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice. (F) Insulin tolerance test of NC-fed Foxp3 (n = 8) and Foxp3 (n = 7) mice. Area under curve (AUC) is shown to the right. (G) Locomotor activity and (H) body weight-adjusted heat production of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice, determined by metabolic cages. (I) Fecal energy and (J) fecal lipid percentage of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; ***, p < 0.001 by unpaired student’s t-test or two-way ANOVA (D, F).
Characterization of Foxp3 mice. (A–C) Flow cytometric quantification of FOXP3-YFP+CD4+ Treg cells (A), CD44+KLRG1+
(B) and CD44+CD103+
(C) effector Treg cells in lymph nodes from Foxp3 and Foxp3 mice. (D) Growth curve of Foxp3 (n = 10) and Foxp3 (n = 9) mice fed with NC. (E) Lean and fat mass of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice. (F) Insulin tolerance test of NC-fed Foxp3 (n = 8) and Foxp3 (n = 7) mice. Area under curve (AUC) is shown to the right. (G) Locomotor activity and (H) body weight-adjusted heat production of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice, determined by metabolic cages. (I) Fecal energy and (J) fecal lipid percentage of NC-fed Foxp3 (n = 7) and Foxp3 (n = 6) mice. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; ***, p < 0.001 by unpaired student’s t-test or two-way ANOVA (D, F).We went on to challenge the mice with HFD. The resistant to HFD-induced weight gain could be readily seen in Foxp3 mice as early as 3 weeks after diet-switching (
). The reduction in weight gain was solely because of reduced fat mass (
). Both BAT and gWAT depots were significantly smaller in Foxp3 mice, while liver weight remained the same (
). We then asked what led to the lean phenotype after STAT5B-CA expression in Treg cells. Continuous metabolic monitoring revealed a substantial reduction of food intake in Foxp3 mice (
), while energy expenditure (
) and locomotor activity (
) were both comparable to wildtype controls. Moreover, there was a tendency that fecal energy excretion was increased after Treg cell hyperactivation (
). We then measured lipid content in the feces and found more lipids were excreted in Foxp3 mice (
). As expected, these HFD-fed but relatively lean Foxp3 mice show much better glucose handling during the glucose tolerance test (
). These data demonstrate that Treg cell expansion and activation stimulated by tonic STAT5 signaling restrain dietary fat intake/absorption and render mice resistant to diet-induced metabolic disturbances.
Figure 4
Protection of Foxp3 mice from diet-induced obesity. (A) Body weight of Foxp3 (n = 9) and Foxp3 (n = 6) mice after HFD feeding. (B) Lean and fat mass of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. (C) Weight of BAT, gonadal WAT, and liver of HFD-fed Foxp3 (n = 5) and Foxp3 (n = 5) mice. (D) Cumulative food intake of Foxp3 (n = 7) and Foxp3 (n = 4) mice when switched to HFD feeding. (E) Body weight-adjusted heat production and (F) locomotor activity of HFD-fed Foxp3 (n = 11) and Foxp3 (n = 6) mice, determined by metabolic cages. (G) Fecal energy and (H) fecal lipid percentage of HFD-fed Foxp3 (n = 8) and Foxp3 (n = 6) mice. (I) Glucose tolerance test of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. Area under curve (AUC) is shown to the right. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test or two-way ANOVA (A, I).
Protection of Foxp3 mice from diet-induced obesity. (A) Body weight of Foxp3 (n = 9) and Foxp3 (n = 6) mice after HFD feeding. (B) Lean and fat mass of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. (C) Weight of BAT, gonadal WAT, and liver of HFD-fed Foxp3 (n = 5) and Foxp3 (n = 5) mice. (D) Cumulative food intake of Foxp3 (n = 7) and Foxp3 (n = 4) mice when switched to HFD feeding. (E) Body weight-adjusted heat production and (F) locomotor activity of HFD-fed Foxp3 (n = 11) and Foxp3 (n = 6) mice, determined by metabolic cages. (G) Fecal energy and (H) fecal lipid percentage of HFD-fed Foxp3 (n = 8) and Foxp3 (n = 6) mice. (I) Glucose tolerance test of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. Area under curve (AUC) is shown to the right. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test or two-way ANOVA (A, I).
Inducible STAT5 Activation in Treg Cells Improves Systemic Glucose Metabolism
Due to stochastic activity of the Foxp3 promoter (43), constitutive overexpression of STAT5B-CA led to progressive splenomegaly and hematopoietic neoplasia in Foxp3 mice (data not shown), similar to that after STAT5 activation in other hematopoietic compartments (44–46). We therefore used the inducible Foxp3 line to generate Foxp3 male mice, fed them with tamoxifen-containing food to induce STAT5-CA expression in Treg cells, and challenged them with HFD feeding (
). A trending reduction in weight gain could be observed in Foxp3 mice (
). When compared to their control counterparts, Foxp3 mice possessed significantly less BAT (
) and slightly smaller WAT (
). No change in muscle mass was noticed (
). Moreover, Foxp3 animals also displayed much improved insulin tolerance (
). Concomitantly, insulin signaling activity, reflected by AKT Ser473 phosphorylation, was augmented in BAT of Foxp3 mice (
). The more profound protection observed in the former constitutive model might stem from stronger and broader expression of STAT5B-CA. Nonetheless, using two complementary models, we demonstrate that STAT5 activation in Treg cells is sufficient to protect mice from diet-induced adiposity and insulin resistance.
Figure 5
Inducible STAT5B-CA expression in Treg cells promotes metabolic adaptation to HFD. (A) Schematic view of the generation of tamoxifen-inducible, Treg cell-specific STAT5B-CA overexpression mice and dietary regimens to the mice. (B) Body weight of HFD-fed Foxp3 (n = 6) and Foxp3 (n = 6) mice. (C–F) Weight of BAT (C), gonadal WAT (D), inguinal WAT (E), and quadricep muscle (F) of Foxp3 (n = 9) and Foxp3 (n = 5) mice. (G) Insulin tolerance test of HFD-fed Foxp3 (n = 8) and Foxp3 (n = 5). Area under curve (AUC) is shown to the right. (H) pS473-Akt and total-Akt immunoblotting of BAT protein lysates from Foxp3 and Foxp3 mice treated with saline or insulin for 30 min. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01 by unpaired student’s t-test, one-way or two-way ANOVA.
Inducible STAT5B-CA expression in Treg cells promotes metabolic adaptation to HFD. (A) Schematic view of the generation of tamoxifen-inducible, Treg cell-specific STAT5B-CA overexpression mice and dietary regimens to the mice. (B) Body weight of HFD-fed Foxp3 (n = 6) and Foxp3 (n = 6) mice. (C–F) Weight of BAT (C), gonadal WAT (D), inguinal WAT (E), and quadricep muscle (F) of Foxp3 (n = 9) and Foxp3 (n = 5) mice. (G) Insulin tolerance test of HFD-fed Foxp3 (n = 8) and Foxp3 (n = 5). Area under curve (AUC) is shown to the right. (H) pS473-Akt and total-Akt immunoblotting of BAT protein lysates from Foxp3 and Foxp3 mice treated with saline or insulin for 30 min. Data are presented as mean ± SEM. n.s., not significant; *, p < 0.05; **, p < 0.01 by unpaired student’s t-test, one-way or two-way ANOVA.
Treg Cell Activation Suppresses Lipid Uptake by Adipose Tissue
We then inquired into the causes of metabolic protection after STAT5-dependent Treg cell overactivation. We first interrogated the involvement of BAT thermogenesis. BAT histology clearly revealed smaller adipocytes and less lipid deposition in Foxp3 mice (
). However, the expression of thermogenic genes remained unchanged in BAT (
), or even reduced in inguinal WAT (iWAT) (
). This suggests that reduced weight gain by Treg cell activation was not a result of increased thermogenesis.
Figure 6
Thermogenesis and lipid metabolism in Foxp3 mice. (A) BAT histology of HFD-fed Foxp3 and Foxp3 mice. (B, C) Expression of thermogenic genes in BAT (B) and iWAT (C) of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (D) Gonadal WAT histology of HFD-fed Foxp3 and Foxp3 mice. (E)
Cd36 and (F)
Lpl gene expression in WAT of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (G) Expression of the Cd36 gene in duodenum of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (H) Oral lipid tolerance test of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. Serum triglyceride levels were determined. Area under curve (AUC) is shown to the right. Data are presented as mean ± SEM. *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test.
Thermogenesis and lipid metabolism in Foxp3 mice. (A) BAT histology of HFD-fed Foxp3 and Foxp3 mice. (B, C) Expression of thermogenic genes in BAT (B) and iWAT (C) of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (D) Gonadal WAT histology of HFD-fed Foxp3 and Foxp3 mice. (E)
Cd36 and (F)
Lpl gene expression in WAT of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (G) Expression of the Cd36 gene in duodenum of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (H) Oral lipid tolerance test of HFD-fed Foxp3 (n = 9) and Foxp3 (n = 6) mice. Serum triglyceride levels were determined. Area under curve (AUC) is shown to the right. Data are presented as mean ± SEM. *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test.Not only was gWAT mass reduced in Foxp3 mice (
), but adipocyte size was also much smaller (
). Crown-like structures were more evident in Foxp3 gWAT (
), suggesting active clearing of dying adipocytes by macrophages (47). Interestingly, the expression of both Cd36 (
) and Lpl (
), two important genes required for cellular uptake of free fatty acid (48), were substantially ablated in gWAT and iWAT of Foxp3 mice. There seemed to be no defect in dietary intake of lipid in the intestine, as intestinal Cd36 expression was intact (
) and blood triglyceride levels after an oral gavage of olive oil were even higher in Foxp3 mice (
). Taking in consideration of the elevated fecal lipid content (
), we postulate that Treg cell activation suppresses lipid uptake into adipose tissues, leading to more energy excretion and less weight gain during HFD feeding.
Treg Cell Activation Remodels Systemic Iron Metabolism
Iron is critically important for adipose tissue function and systemic metabolism. Dietary iron supplementation augments thermogenesis, reduces fat mass, and protects mice from diet-induced weight gain (49, 50). On the other hand, nonanemic iron deficiency impairs thermogenesis and promotes visceral obesity in mice (51). Lowering adipocyte iron uptake via genetic deletion of transferrin receptor 1 (TFR1) leads to impaired thermogenesis, increased insulin resistance, and inflammation (52). We then inquired if the metabolic protection after STAT5-dependent Treg cell activation involves iron. Ferritin is a protein complex that stores cellular iron and its protein levels were elevated in BAT, gWAT and iWAT of HFD-fed Foxp3 mice, compared to Cre-only controls (
), indicative of excess iron storage in fat tissues. In Foxp3 adipose tissues, there was a trending increase in Tfr1 gene expression (
) and TFR1 protein levels (
). Conversely, expression of the Fpn1 gene (encoding iron export protein Ferroportin) was significantly downregulated (
). In contrast to adipose tissue, both liver and spleen showed substantial reduction in Ferritin protein levels (
), accompanied with increased TFR1 expression. The presumptive impairment in hepatic iron storage was correlated with reduced expression of iron-responsive Hamp and Hfe genes (
), which encode Hepcidin and Homeostatic iron regulator, respectively. As a result of reduced Hepcidin, intestinal genes involved in iron uptake including Dmt1, Cybrd1, Fn1 displayed varying degrees of upregulation (
). Collectively, STAT5B-CA-mediated Treg cell overactivation redistributes the iron pool from liver and spleen to adipose tissues, which may contribute to the alleviation of metabolic stress associated with diet-induced obesity.
Figure 7
Iron metabolism modulated by Treg cell activation. (A–C) Immunoblotting of Ferritin, TFR1, FPN1 proteins in BAT (A), gWAT (B), and iWAT (C) of HFD-fed Foxp3 and Foxp3 mice. (D, E) Expression of the Tfr1
(D) and Fpn1
(E) genes in adipose tissues of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (F, G) Immunoblotting of Ferritin, TFR1, FPN1 proteins in liver (F) and spleen (G) of HFD-fed Foxp3 and Foxp3 mice. (H) Expression of the Hamp1 and Hfe genes in the liver of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (I) Expression of the Dmt, Cybdr1, and Fpn1 genes in the duodenum of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. Data are presented as mean ± SEM. *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test.
Iron metabolism modulated by Treg cell activation. (A–C) Immunoblotting of Ferritin, TFR1, FPN1 proteins in BAT (A), gWAT (B), and iWAT (C) of HFD-fed Foxp3 and Foxp3 mice. (D, E) Expression of the Tfr1
(D) and Fpn1
(E) genes in adipose tissues of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (F, G) Immunoblotting of Ferritin, TFR1, FPN1 proteins in liver (F) and spleen (G) of HFD-fed Foxp3 and Foxp3 mice. (H) Expression of the Hamp1 and Hfe genes in the liver of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. (I) Expression of the Dmt, Cybdr1, and Fpn1 genes in the duodenum of HFD-fed Foxp3 (n = 7) and Foxp3 (n = 5) mice. Data are presented as mean ± SEM. *, p < 0.05; **, p < 0.01; ***, p < 0.001 by unpaired student’s t-test.
Discussion
In effector Treg cells, protein O-GlcNAcylation induced by the TCR signaling, modifies and activates STAT5 for Treg cell homeostasis and function (38). In the study, we investigate whether and how the OGT-STAT5 regulatory axis dictates Treg cell function during metabolic adaptation to diet-induced obesity. Taking advantage of random X chromosome inactivation and mosaic OGT-deficiency in Foxp3 female mice, we establish that functional deficiency in about half the Treg cells renders animals more prone to diet-induced adiposity and insulin resistance. Our model is more advantageous than other Treg cell depletion models, such as the Foxp3-DTR mice (14), since the absence of systemic inflammation in Foxp3 females allowed us to perform long-term studies and avoided side effects and caveats of tissue inflammation (38). This is possibly also the reason why we observed obesity and thermogenic defects in Foxp3 mice, which have not been reported in other loss-of-function animals (14, 18, 22). We speculate that other heterozygous loss of Treg function models would show similar metabolic phenotypes as the Foxp3 mice, provided that systemic inflammation is avoided. In addition, we do not imply that Treg function in coping with obesity is specifically controlled by OGT, but rather OGT is essential for most if not all Treg functions, including their regulation of insulin sensitivity and systemic metabolism. For instance, we previously demonstrated that Blimp1+ effector Treg cells are ablated in OGT knockout mice (38). Similar to OGT knockout, Blimp1-deficient mice do not have effector Tregs (53). Recently, it was shown that Blimp1+ Tregs are also involved in the homeostasis of adipose tissue (54).It is important to note that BAT harbors Treg cell populations that have shared and distinct gene signatures compared to WAT-resident Treg cells (55). Cold exposure or β3-adrenergic stimuli enhance Treg induction that is linked by a T cell-specific STAT6-Pten axis (56). When OGT is deficient in Treg cells, BAT insulin resistance and impaired thermogenesis could be observed, contributing to diet-induced obesity. Nonetheless, it is unclear if cold or β3-adrenergic activation could directly modulate the hexosamine biosynthetic pathway, OGT activity, and protein O-GlcNAcylation in Treg cells. Other work in our laboratory has revealed a pivotal role of O-GlcNAcylation on STAT6 signaling and type 2 immunity (57). It is warranted in the future to determine if OGT also regulates the STAT6-Pten axis in tissue-specific Treg cells and adipose tissue function.IL-2 and IL-2R-dependent STAT5 activation are required for Treg cell development (5, 58). Foxp3-mediated STAT5a/b knockout in male mice causes scurfy phenotypes and death by the age of 2 months (data not shown). We also generated mosaic Foxp3 female mice avoiding fatal inflammation. However, we did not observe evident metabolic perturbations in these animals (data not shown), differently from Foxp3 females. It indicates that the metabolic effects of Treg cell O-GlcNAcylation might not be mediated just by STAT5, and other mechanisms such as STAT6 could compensate for the loss of STAT5. Nevertheless, STAT5B-CA overexpression is sufficient to booster Treg cell number and immunosuppression (10). Using these Foxp3 mice, we demonstrate that Treg cell activation renders mice resistant to diet-induced weight gain and insulin resistance, in part by suppressing fatty acid uptake in adipose tissue and remodeling systemic iron metabolism. Intriguingly, no substantial changes in BAT thermogenesis were observed in STAT5B-CA overexpression mice, suggesting that STAT5-independent, but maybe STAT6-dependent mechanisms in Treg cells control adipose tissue lipolysis and thermogenesis (56).While it has been widely accepted that Treg cells are important contributors to metabolic homeostasis, the regulatory mechanisms are not fully understood. Previous studies have revealed inflammation, lipolysis, and thermogenesis as tissue-specific targets of Treg cells. Our current study suggests two additional metabolic actions of Treg cells – fatty acid uptake and iron distribution. While intestinal lipid absorption is not perturbed, STAT5B-CA-driven Treg cell activation suppresses the expression of key genes involved in adipose fatty acid uptake – Lpl and Cd36. It would be important in the future to determine if the obesity protection in STAT5B-CA overexpression mice is dependent on LPL/CD36 downregulation. Meanwhile, the adipose iron pool seems to be expanded in STAT5B-CA overexpression, which we believe contributes to thermogenic activation and fat mass reduction. However, we are unclear whether parenchymal adipocytes or resident immune cells such as macrophages play a predominant role in the change in iron content. We also lack a mechanistic understanding of how Treg cells differentially regulate iron metabolism across different tissues (adipose vs. liver and kidney). Finally, it would be interesting to test if the suppression of lipid uptake is linked to iron redistribution in adipose tissues. In summary, our study has revealed new biology of Treg cells in metabolic regulation and the OGT-STAT5 axis as a target to harness Treg cell function for obesity and diabetes therapeutics.
Data Availability Statement
The original contributions presented in the study are included in the article/
. Further inquiries can be directed to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by University of Minnesota.
Author Contributions
ZZ and MF generated and analyzed STAT5B-CA mice. OS and KH did flow cytometry. BL analyzed OGT-KO mice. ZM helped with Western blot. H-BR conceived the project, analyzed data and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by NIH grant R56 AI162791 to H-BR and R01 AI139420 to H-BR and MF.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Authors: Matthias H Tschöp; John R Speakman; Jonathan R S Arch; Johan Auwerx; Jens C Brüning; Lawrence Chan; Robert H Eckel; Robert V Farese; Jose E Galgani; Catherine Hambly; Mark A Herman; Tamas L Horvath; Barbara B Kahn; Sara C Kozma; Eleftheria Maratos-Flier; Timo D Müller; Heike Münzberg; Paul T Pfluger; Leona Plum; Marc L Reitman; Kamal Rahmouni; Gerald I Shulman; George Thomas; C Ronald Kahn; Eric Ravussin Journal: Nat Methods Date: 2011-12-28 Impact factor: 28.547
Authors: Ming Zhao; Kaiqun Ren; Xiwen Xiong; Yue Xin; Yujie Zou; Jason C Maynard; Angela Kim; Alexander P Battist; Navya Koneripalli; Yusu Wang; Qianyue Chen; Ruyue Xin; Chenyan Yang; Rong Huang; Jiahui Yu; Zan Huang; Zengdi Zhang; Haiguang Wang; Daoyuan Wang; Yihui Xiao; Oscar C Salgado; Nicholas N Jarjour; Kristin A Hogquist; Xavier S Revelo; Alma L Burlingame; Xiang Gao; Jakob von Moltke; Zhaoyu Lin; Hai-Bin Ruan Journal: Immunity Date: 2022-04-05 Impact factor: 43.474
Authors: Barbara Maurer; Harini Nivarthi; Bettina Wingelhofer; Ha Thi Thanh Pham; Michaela Schlederer; Tobias Suske; Reinhard Grausenburger; Ana-Iris Schiefer; Michaela Prchal-Murphy; Doris Chen; Susanne Winkler; Olaf Merkel; Christoph Kornauth; Maximilian Hofbauer; Birgit Hochgatterer; Gregor Hoermann; Andrea Hoelbl-Kovacic; Jana Prochazkova; Cosimo Lobello; Abbarna A Cumaraswamy; Johanna Latzka; Melitta Kitzwögerer; Andreas Chott; Andrea Janikova; Šárka Pospíšilova; Joanna I Loizou; Stefan Kubicek; Peter Valent; Thomas Kolbe; Florian Grebien; Lukas Kenner; Patrick T Gunning; Robert Kralovics; Marco Herling; Mathias Müller; Thomas Rülicke; Veronika Sexl; Richard Moriggl Journal: Haematologica Date: 2020-01-31 Impact factor: 9.941
Authors: Bing Liu; Oscar C Salgado; Sangya Singh; Keli L Hippen; Jason C Maynard; Alma L Burlingame; Lauren E Ball; Bruce R Blazar; Michael A Farrar; Kristin A Hogquist; Hai-Bin Ruan Journal: Nat Commun Date: 2019-01-21 Impact factor: 14.919