Abstract Macrophages may acquire a reparative phenotype that supports tissue repair and remodeling in response to tissue injury. However, the metabolic requirements underpinning this process are incompletely understood. Here, we show that posttranslational modification (PTM) of PPARγ regulates lipid synthesis in response to wound microenvironmental cues and that metabolic rewiring orchestrates function of reparative macrophages. In injured tissues, repair signaling leads to decreased macrophage PPARγ threonine 166 (T166) phosphorylation, which results in a partially active PPARγ transcriptional program comprised of increased binding activity to the regulator regions of lipid synthesis-associated genes, thereby increased lipogenesis. The accumulated lipids serve as signaling molecules, triggering STAT3-mediated growth factor expression, and supporting the synthesis of phospholipids for the expansion of the endoplasmic reticulum (ER), which is required for protein secretion. Genetic or pharmacological inhibition of PPARγ T166 phosphorylation promotes the reparative function of macrophages and facilitates tissue regeneration. In summary, our work identifies PPARγ T166-regulated lipid biosynthesis as an essential pathway for meeting the anabolic demands of the activation and function of macrophages and provides a rationale for potential therapeutic targeting of tissue repair. Subject terms: Monocytes and macrophages, Phosphorylation, Lipid signalling, Skin diseases __________________________________________________________________ Macrophages with a reparative phenotype are important for tissue repair and have distinctive metabolic features. Here authors show that metabolic rewiring of macrophages during wound healing involves dephosphorylation of the transcription factor PPARγ, which results in activation of target genes that regulate lipid biosynthesis. Introduction The restoration of tissue integrity and homeostasis after injury is a complex biological process that proceeds in the following sequential phases: inflammation, tissue formation and maturation^[48]1. Among all cell types, macrophages play critical roles in supporting each phase of wound repair by rapidly adopting diverse activation states in response to local tissue microenvironment cues^[49]2,[50]3. When damage occurs, macrophages polarize to a classically activated phenotype and exhibit proinflammatory functions such as phagocytosis, antigen presentation, and proinflammatory cytokine production^[51]4. When the wound-healing response is well organized and controlled, macrophages shift toward a resolution phenotype to inhibit inflammation by releasing anti-inflammatory cytokines^[52]2,[53]5. Although early studies investigating the contribution of macrophages to wound repair focused on their role as immune cells, it is now clear that macrophages exhibit much more complex roles not only in inflammation control but also in the mechanisms of tissue regeneration^[54]6–[55]8. Macrophages are an important source of growth factors, including platelet-derived growth factor (PDGF), transforming growth factor-β1 (TGF-β1), and vascular endothelial growth factor (VEGF), which induce re-epithelialization and neovascularization by promoting cellular proliferation and blood vessel development^[56]4,[57]9. Macrophage depletion leads to less efficient repair and regeneration following injury^[58]10. However, the precise molecular mechanisms that drive and sustain the macrophage reparative phenotype remain poorly understood. There is a growing understanding that cellular metabolic profiles are associated with the balance of macrophage activation and function^[59]11. Macrophages with different activation states exhibit distinct metabolic programs in response to activation to meet their bioenergetic needs^[60]11,[61]12. In addition, metabolic intermediates can function as signaling molecules to facilitate the polarization of macrophages^[62]13,[63]14. For instance, lipopolysaccharide (LPS) strongly induces glycolysis in classically activated macrophages to meet energy demands^[64]12,[65]15. Moreover, this process increases succinate levels, which block hypoxia-inducible factor-1α (HIF-1α) degradation to sustain the proinflammatory phenotype^[66]16,[67]17. However, the metabolic requirements of macrophages in the repair phase remain unclear. Studies of metabolic programs fueling macrophage activation have largely focused on catabolic pathways, such as glycolysis and oxidative phosphorylation (OXPHOS)^[68]18, while anabolic processes, including the synthesis of lipids, amino acids, and nucleotides, are also essential for proliferation and cytokine production in immune cells^[69]19–[70]21. The production and secretion of repair factors are crucial for the function of reparative macrophages. It is imperative to further investigate the anabolic processes of macrophages. To determine the functional interplay between cellular metabolism and specific repair effector functions in injury-associated macrophages, we conduct a comprehensive analysis of these macrophages at the transcriptional and functional levels. The single-cell transcriptome atlas reveals that fatty acid synthesis (FAS) activities are consistently enhanced following stimulation with healing signals in 9 tissue injury models. This metabolic reprogramming is precisely regulated by PPARγ T166 (corresponding to the T136 site in PPARγ1) phosphorylation. Intracellular fatty acids induce STAT3-mediated growth factor expression and support the expansion of the ER, which is required for the secretion of reparative proteins. Genetic and pharmacological inhibition of PPARγ T166 phosphorylation markedly improve tissue repair in vivo and the migration of primary dermal fibroblasts (DF) in vitro. Collectively, our results identify a function and mechanism of PPARγ in reparative macrophages and open avenues for therapeutically modulating macrophage polarization to promote wound repair. Results The increase in lipid anabolic metabolism is a feature of reparative macrophages during wound healing To explore the metabolic pathways that are used by reparative macrophages, we obtained single-cell RNA sequencing (scRNA-seq) data on monocyte-macrophage lineage cells from mouse wound healing models in 9 tissue types, including airway, joint, heart, kidney, liver, lung, skeletal muscle, skin, and spinal cord (Fig. [71]1a). After strict quality control and filtration, we collected a total of 59332 macrophages derived from the inflammatory phase (phase I) and reparative phase (phase II) (Supplementary Fig. [72]1a, b). Considering the heterogeneity of macrophages and batch effects among different datasets, we analyzed each dataset independently. Using skin injury as an example, we found that inflammatory phase and reparative phase macrophages clustered predominantly in disparate populations, as shown by uniform manifold approximation and projection (UMAP) (Fig. [73]1b). Gene ontology (GO) enrichment analysis showed that inflammatory phase macrophages had significantly upregulated genes associated with glycolytic processes, while reparative phase macrophages highly expressed genes associated with lipid metabolic processes (Fig. [74]1c). To verify this metabolic difference, we sorted wound macrophages in the inflammatory and reparative phases in a full-thickness excisional skin injury model. As expected, we found significantly elevated expression of genes associated with lipid metabolism in skin macrophages in the reparative phase (Fig. [75]1d). Intriguingly, we observed consistent upregulation of FAS activity in reparative phase macrophages in all injury models, while the other types of metabolism exhibited various regulation patterns in phase II macrophages (Fig. [76]1e and Supplementary Fig. [77]1c). The differences in metabolic features indicate the flexibility of macrophages in their energy usage to meet the bioenergetic demands of cellular activation, but FAS is a more conserved metabolic pathway for the activation of reparative macrophages. Fig. 1. Fatty acid synthesis is a feature of reparative macrophages. [78]Fig. 1 [79]Open in a new tab a Tissue injury models included in the analysis. b The scRNA-seq analysis of skin wound macrophages of phase I (1 dpi) and phase II (6 dpi). UMAP highlights wound healing phases. c GO enrichment analysis showing the highly enriched metabolic pathway associated with phase I and phase II. Pie charts represent enriched biological processes and the ratio of gene number that enriched in each group, small dots represent genes associated with the processes. The size of each pie indicates the number of genes enriched in a specific term. d The mRNA expression of lipid metabolism-associated genes in sorted wound macrophages isolated at 1 and 6 dpi (n = 3 mice per group). e Heatmaps showing ssGSEA scores of phase I and phase II macrophages using signature genes associated with FAO (mmu00071) and FAS (mmu00061) pathways. f Developmental trajectory of skin wound macrophages. Sample are colored by wound healing phases or macrophage subtypes. g Heatmap showing expression patterns of glycolysis-, tricarboxylic acid cycle (TCA)-, OXPHOS-, pentose phosphate pathway (PPP)-, FAO-, and FAS-associated signature genes among macrophage subtypes from [80]GSE131364 dataset. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. Considering the strong influences of the tissue microenvironment, we further confirmed lipid synthesis activity in reparative macrophages in vitro. Polarized macrophages in vitro are instrumental in analyzing macrophage function, and we tested polarizing conditions to generate reparative macrophages. We used a skin injury model and found a significant increase in the secretion of IL-10 during the tissue repair phase (Supplementary Fig. [81]2a). SingleR revealed that 58.9% of macrophages in the reparative phase exhibited an IL-10-induced phenotype (Fig. [82]1f). IL-10-treated bone marrow-derived macrophages (BMDM) showed upregulated expression of tissue repair-related genes (Supplementary Fig. [83]2b), and significantly promoted the migration of skin DFs (Supplementary Fig. [84]2c). These results suggest that IL-10 polarizes macrophages toward a typical reparative phenotype. We retrieved public transcriptomes of IL-10-induced macrophages M(IL10) and showed that M(IL10) exhibited a distinct metabolic pattern compared to M0 (undifferentiated), M(LPS) (classically activated), or M(IL4) (alternatively activated) macrophages (Fig. [85]1g). Interestingly, M(IL10) and M(IL4) exhibited increased lipid metabolic activity, while M(IL4) showed higher expression levels of fatty acid oxidation (FAO)-related genes, and M(IL10) macrophages had strongly upregulated FAS genes (Fig. [86]1g). We also observed that IL-10 enhanced FAS-related gene expression in human macrophages (Supplementary Fig. [87]1d, e). In addition, previously reported reparative macrophage induction strategies, that used dexamethasone or TGF-β1, induced highly expressions of lipid metabolism-associated genes, especially FAS-related genes in macrophages (Supplementary Fig. [88]1f). Collectively, our data revealed that reparative macrophages possess a distinct metabolic program characterized by a significant increase in lipid synthesis-related genes. Lipid synthesis is required for the activation and function of reparative macrophages To investigate the lipid metabolic impact of reparative signaling in macrophages, we performed lipidomic analysis (Supplementary Fig. [89]3a). We observed that glycerolipids were highly accumulated in M(IL10) (Fig. [90]2a). Pathway enrichment analysis validated the alteration of glycerolipid metabolism in IL-10-induced macrophages (Fig. [91]2b). Consistently, M(IL10) exhibited markedly increased staining with neutral lipid dye (BODIPY 493/503), which is indicative of increased lipid accumulation (Fig. [92]2c, d). Similarly, reparative phase (6 days post injury, 6 dpi) wound macrophages accumulated more lipids than macrophages isolated from the early phase (1 dpi) of skin injury (Supplementary Fig. [93]3b). Stimulating macrophages with IL-10 increased intracellular triglycerides (TG) and nonesterified fatty acids (NEFA) (Supplementary Fig. [94]3c). These results indicate that reparative signaling induces lipid accumulation in macrophages. Fig. 2. Limiting de novo lipid synthesis reduces reparative function of macrophages. [95]Fig. 2 [96]Open in a new tab a Volcano plot showing changes of lipid metabolites in BMDMs with or without IL-10 treatment (n = 4 per group). b GO enrichment analysis showing the pathway alterations in BMDMs with or without IL-10 treatment. c Representative images of BODIPY staining in control or IL-10-treated BMDMs (n = 5 biological replicates). Scale bar, 10 μm. d BODIPY staining analysis by flow cytometry in control or IL-10-treated macrophages. e Protein expression of ACC1 and FASN in BMDMs in response to IL-10. β-Actin was used as loading control. f Lipid synthesis rate in control and IL-10 BMDMs. The data represent the incorporation of ^14C-glucose in the lipid fraction of n = 4 biological replicates per group. g Diagram of de novo synthesis of fatty acids from [U-^13C] glucose. h The percentages of isotopomers of FA 16:0 (palmitate) after trace to [U-^13C] glucose in BMDMs in the presence or absence of IL-10. n = 3 (control) or n = 4 (IL-10) biological replicates. i Flow cytometry assay of CD163 and CD206 on BMDMs with IL-10 or C75 treatment (n = 3 biological replicates). j TGF-β1, PDGFβ and VEGFβ in cell culture medium is determined by ELISA upon stimulation with IL-10 or C75. n = 6 (TGF-β1 and PDGFβ) or n = 4 (VEGFβ) independent experiments. k The scratch wound closure rate of DFs co-cultured with M(IL10) or treatment with C75 or M(IL10) knockout (KO) Fsan for 24 h (n = 3 biological replicates). Data were analyzed by one-way analysis of variance (ANOVA) followed by a Dunnett’s multiple-comparisons test. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. In line with an accumulation of lipidome, the protein expression of metabolic enzymes involved in FAS (ACC1 and FASN) was also strongly induced in M(IL10) and reparative phase macrophages of skin injury (Fig. [97]2e and Supplementary Fig. [98]3d). Considering that glucose is the major substrate for FAS, we traced the fate of glucose in reparative macrophage using isotope labeling. Using C^14-labeled glucose, we observed that IL-10-stimulated macrophages substantially increased C^14-glucose incorporation into lipid components (Fig. [99]2f). We further observed that M(IL10) cells had significantly increased palmitic acid and stearic acid labeling from [U-^13C] glucose tracers (Fig. [100]2g, h and Supplementary Fig. [101]3e). These data confirm that the accumulated lipids in reparative macrophages are derived from de novo lipid synthesis pathway. To further investigate the lipid synthesis activity of reparative macrophages in wound healing processes, we generated cardiotoxin (CTX)-induced skeletal muscle injury and CCl[4]-induced acute liver injury models (Supplementary Fig. [102]4a–d). Consistent with the findings observed in skin injury models, reparative macrophages from both muscle and liver injury models exhibited enhanced lipid accumulation and upregulated expression of lipid synthesis-related genes (Supplementary Fig. [103]4e–h). Together, these results demonstrate that increased de novo lipid synthesis activity is a feature of reparative macrophages. We then examined whether perturbations in the de novo lipid synthesis in macrophages impacted reparative functions. C75, a competitive fatty acid synthase (FASN) inhibitor, significantly inhibited IL-10-induced CD163 and CD206 expression (Fig. [104]2i), and impaired the secretion of reparative cytokines by macrophages (Fig. [105]2j and Supplementary Fig. [106]5a). Knockout of FASN encoding gene Fasn also significantly inhibited the reparative phenotype of M(IL10) (Supplementary Fig. [107]5b, c). Importantly, both FASN inhibition and Fasn knockout in M(IL10) decreased dermal fibroblasts migration in scratch tests (Fig. [108]2k and Supplementary Fig. [109]5d, e). Collectively, these results confirm the strong upregulation of lipid synthesis and demonstrate that de novo lipid synthesis is required for the reparative function of macrophages. PPARγ T166 phosphorylation regulates lipid synthesis in reparative macrophages To investigate the regulator of lipid metabolism in reparative macrophages, we performed transcription factor activity analysis using SCENIC. We found that the transcriptional activity of 3 lipid metabolism-related transcription factors, including CEBPδ, PPARγ and SREBF2, increased in reparative phase macrophages (Fig. [110]3a and Supplementary Fig. [111]6a, b). However, only PPARγ transcriptional activity positively correlated with FAS in macrophages (Fig. [112]3b). The activation of PPARγ in reparative macrophages was validated by a fluorescence reporter assay (Fig. [113]3c). Using PPARγ-knockout macrophages, we confirmed that PPARγ regulated lipid synthesis activity in macrophages (Supplementary Fig. [114]6c, d). In addition, PPARγ deficiency downregulated the expression of reparative markers CD163 and CD206, while upregulated the expression of proinflammatory makers CD80 and MHCII in macrophages (Supplementary Fig. [115]6e), further indicating the strong association between lipid synthesis activity and reparative phenotype in macrophages. These results suggest that lipid synthesis in reparative macrophages is regulated by PPARγ activity. Then, we investigated the molecular mechanism through which PPARγ could regulate lipid synthesis in reparative macrophages. First, we observed no significant changes in the expression levels or cellular localization of PPARγ in macrophages after IL-10 treatment (Supplementary Fig. [116]6f). Considering that PPARγ activity is heavily regulated by posttranslational modification (PTM)^[117]22, we next measured three important phosphorylation sites (S112, T166 and S273) of PPARγ. We found that the levels of T166 phosphorylation gradually increased with increasing inflammatory stimulation time, and decreased in response to IL-10 stimulation (Fig. [118]3d and Supplementary Fig. [119]6g). The levels of the other two phosphorylation sites (S112 and S273) showed no significant changes in response to IL-10 (Supplementary Fig. [120]6h), suggesting that reparative signaling induces PPARγ activity through dephosphorylating T166 site in macrophages. Fig. 3. PPARγ T166 dephosphorylation regulates lipid synthesis activity of skin wound macrophages. [121]Fig. 3 [122]Open in a new tab a UMAP and violin plots showing the PPARγ activity in macrophages from skin wound healing phase I and phase II. b Correlations between FAS-related gene expression pattern and transcriptional activity of CEBPδ, PPARγ and SREBF2 in skin wound macrophages. c Schematic illustration indicating the reporter plasmid structure (left), representative flow cytometry histograms showing the reporter fluorescence levels (right) (n = 4 biological replicates). d Representative immunoblot of phospho-T166-PPARγ (pPPARγ^T166) and total PPARγ in BMDMs treated with IL-10. β-actin was used as a loading control (n = 3 biological replicates). Data were analyzed by a one-way ANOVA followed by a Dunnett’s multiple-comparisons test. e Heatmap showing expression patterns of FAS-associated signature genes between PPARγ^+/+ and PPARγ^A/A BMDMs or skin wound macrophages (6 dpi). f Protein expression of ACC1 and FASN in BMDMs from PPARγ^+/+ or PPARγ^A/A mice in response to IL-10. β-actin was used as loading control. n = 3 biological replicates. g Representative pictures of BODIPY staining in PPARγ^+/+ and PPARγ^A/A BMDMs in response to IL-10. The lipid droplets were quantified (n = 5 biological replicates). Scale bar, 10 μm. Data were analyzed by a one-way ANOVA followed by a Tukey’s multiple comparisons test. h Lipid droplets in PPARγ^+/+ and PPARγ^A/A BMDMs were quantified by flow cytometry. Dashed line indicates peak of the fluorescent signal. i Lipid synthesis rate in PPARγ^+/+ and PPARγ^A/A BMDMs. The data represent the incorporation of ^14C-glucose in the lipid fraction of n = 4 biological replicates per group. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. To further determine the role of PPARγ T166 phosphorylation in macrophages, we generated homozygous knock-in mice (PPARγ^T166A/T166A, hereafter referred to as PPARγ^A/A) to abrogate the phosphorylation of PPARγ at T166. We performed RNA sequencing on wild-type (PPARγ^+/+) and PPARγ^A/A BMDMs. The results revealed distinct transcriptional features in PPARγ^+/+ and PPARγ^A/A macrophages (Supplementary Fig. [123]7a), and PPARγ^A/A macrophages negatively enriched with pro-inflammation-related GO items (Supplementary Fig. [124]7b), indicating T166A induces a typical reparative phenotype in macrophages. A similar result was observed in macrophages isolated from reparative phase (6 dpi) wound (Supplementary Fig. [125]7c, d). Then, we found that both PPARγ^A/A BMDMs and skin wound macrophages exhibited significantly upregulated levels of and FAS activity (Fig. [126]3e and Supplementary Fig. [127]7e). Consistent with these findings, PPARγ T166 dephosphorylation increased the expression of ACC1 and FASN in macrophages (Fig. [128]3f). There were also noticeable increases in intracellular TGs and lipid droplets in PPARγ^A/A macrophages compared to wild-type cells (Fig. [129]3g, h and Supplementary Fig. [130]7f). Moreover, lipogenesis assays using C^14-labeled glucose showed increased incorporation of labeled glucose-derived carbon into lipids in PPARγ^A/A macrophages compared to control cells (Fig. [131]3i). Collectively, we reveal that PPARγ T166 dephosphorylation is critical in promoting lipid generation in reparative macrophages. Blocking PPARγ T166 phosphorylation accelerates wound healing We further examined the role of PPARγ T166 dephosphorylation in the tissue repair function of macrophages. A full-thickness excisional skin injury model was created on the dorsal skin of PPARγ^+/+ and PPARγ^A/A mice. The wound healing rates were monitored over time (Fig. [132]4a). Notably, no histological changes in the skin of PPARγ^A/A mice were observed under homeostatic conditions (Supplementary Fig. [133]8a), while PPARγ^A/A mice had smaller skin wounds than PPARγ^+/+ mice from day 3 to 12 post injury (Fig. [134]4b, c). Histologic assessment of the wounds revealed an increase in re-epithelialization in PPARγ^A/A mice (Supplementary Fig. [135]8b), indicating an accelerated repair process in PPARγ^A/A mice. Importantly, skin inflammation significantly increased PPARγ T166 phosphorylation, and as the wound progressed into the repair phase, T166 phosphorylation gradually decreased (Fig. [136]4d). Furthermore, we tested the healing rates of PPARγ^A/A mice using the other injury models. In muscle injury, the cross-sectional area (CSA) of the regenerating muscle fibers was significantly larger in the PPARγ^A/A mice compared to PPARγ^+/+ mice on 8 dpi (Supplementary Fig. [137]9a, b). In acute liver injury models, CCl[4]-induced serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities were significantly higher in PPARγ^+/+ mice on 4 dpi (Supplementary Fig. [138]9c, d). Besides, wound tissues from PPARγ^A/A mice showed a higher abundance of reparative macrophages in both muscle and liver injury models (Supplementary Fig. [139]9e, f). Overall, these data demonstrate the pathophysiological relevance of PPARγ T166 dephosphorylation and wound healing. Fig. 4. Blockage of PPARγ T166 phosphorylation improves skin wound healing. [140]Fig. 4 [141]Open in a new tab a Schematic illustration indicating the phases of skin inflammation, tissue repair after injury. b Representative time-lapse photographs of skin wounds in PPARγ^+/+ and PPARγ^A/A mice over 12 dpi. c Wound area and closure rate in PPARγ^+/+ and PPARγ^A/A mice over 12 dpi (n = 5 per genotype). d Protein levels of p-PPARγ^T166 and PPARγ were determined by western blots using the wounds from 1, 6 and 12 dpi mice. Control, 1 dpi and 12 dpi, n = 3 biologically independent mice per condition. 6 dpi, n = 4 biologically independent mice. Data were analyzed by a one-way ANOVA followed by a Dunnett’s multiple-comparisons test. e Percentage of CD80^+ and CD206^+ macrophages in wound tissue of 1-, 6- and 12-dpi in PPARγ^+/+ and PPARγ^A/A mice within the CD45^+ immune cell population. n = 3 biological replicates. f The quantification of wound area in PPARγ^+/+ and PPARγ^A/A mice after treatment with CL-lipo and Lipo control. n = 4 (PPARγ^+/+ and PPARγ^+/+ + CL), n = 7 (PPARγ^A/A and PPARγ^A/A + CL) biologically independent samples. Data were analyzed by a one-way ANOVA followed by a Tukey’s multiple comparisons test. g The scratch wound closure rate of DFs co-cultured with PPARγ^+/+ or PPARγ^A/A BMDMs for 24 h (n = 3 per group). Data were analyzed by a one-way ANOVA followed by a Dunnett’s multiple-comparisons test. h Representative flow cytometry histograms and quantifications of surface levels of CD163 and CD206 in PPARγ^+/+ and PPARγ^A/A BMDMs treated with IL-10 for 48 h. n = 3 biological replicates per group. Data were analyzed by a one-way ANOVA followed by a Tukey’s multiple comparisons test. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. Next, we investigated the linkage between tissue repair and the cutaneous immune response in PPARγ^A/A mice. As expected, injury led to the recruitment of neutrophils, monocytes, and macrophages at the wound site, but the absolute numbers of monocytes/macrophages, T lymphocytes and granulocytes were unchanged between the genotypes after injury (Supplementary Fig. [142]10a). Strikingly, PPARγ^A/A mice showed a decrease in MHCII^+ macrophage numbers but a significant increase in the number of CD206^+ macrophages around the wound compared to those in PPARγ^+/+ mice 6 dpi (Supplementary Fig. [143]10b). In addition, macrophages from wild type mice showed higher levels of the inflammatory M1 marker CD80 at 1 dpi, while macrophages from PPARγ^A/A mice showed significantly higher surface expression of the reparative marker CD206 at 1 and 6 dpi (Fig. [144]4e). Genes encoding anti-inflammatory and reparative molecules, including Arg1, Tgfb1, Col2a1 and Vegfb, were upregulated in PPARγ^A/A mice skin, whereas inflammatory cytokines and chemokines, such as Il6 and Ccl2, were downregulated in PPARγ^A/A mice skin in the reparative phase (Supplementary Fig. [145]10c). These data reveal that PPARγ T166 dephosphorylation promotes wounding healing by controlling macrophage reparative phenotype. To determine the contribution of macrophages in these transgenic mice, we used a clodronate liposome-based approach to deplete macrophages from the dorsal skin of mice. As expected, intradermal injection of clodronate-encapsulated liposomes (CL-lipos) slowed the wound healing rate in wild-type mice, suggesting a crucial role for macrophages in regulating skin regeneration (Fig. [146]4f and Supplementary Fig. [147]11a). Importantly, the comparable repair rates between PPARγ^A/A mice and wild-type mice following macrophage depletion suggests a critical role of PPARγ T166 dephosphorylation in macrophages during tissue repair (Fig. [148]4f and Supplementary Fig. [149]11a). Moreover, compared to wild-type macrophages, PPARγ^A/A macrophages more rapidly promote the migration of skin fibroblasts in vitro (Fig. [150]4g and Supplementary Fig. [151]11b). We further examined the reparative function of PPARγ^A/A macrophages in vitro. Intriguingly, we found that IL-10 stimulation did not further increase the levels of CD206 and CD163 in PPARγ^A/A macrophages, suggesting that the reparative phenotype of PPARγ^A/A macrophages are fully polarized (Fig. [152]4h). Then, we mimicked the tissue injury phase transition in vitro by shifting the stimulation from LPS + IFNγ to IL-10 (Supplementary Fig. [153]11c). We found that PPARγ^A/A BMDMs expressed significantly increased levels of reparative and anti-inflammatory genes, while the expressions of inflammatory genes were lower in PPARγ^A/A macrophages (Supplementary Fig. [154]11d), suggesting that PPARγ T166 dephosphorylation polarize macrophages toward a robust reparative phenotype. T166 dephosphorylation promotes PPARγ binding to lipid synthesis-related gene loci To further understand the role of PPARγ T166A-mediated macrophage reparative activity, we performed chromatin-immunoprecipitation sequencing (ChIP-seq) of PPARγ^+/+ and PPARγ^A/A macrophages. We observed no binding of endogenous PPARγ to reparative phenotype-related gene, such as Cd163, Mrc1 and Arg1 (Supplementary Fig. [155]12a). Instead, the occupancy of PPARγ binding sites was highly correlated with lipid metabolism (Fig. [156]5a). Both wild type and T166A PPARγ binding was highly enriched at promoters of FAO-related genes, while the regulation ability of lipid synthesis and localization was significantly increased in the mutated PPARγ (Fig. [157]5a). Representative gene tracks of PPARγ bindings are shown in Fig. [158]5b. Subsequent ChIP-qPCR and CUT&RUN analysis confirmed the increases in PPARγ binding, and further revealed that the occupancy of PPARγ binding sites at Acly, Acaca, Fasn, Dgat1 and Dgat2 regulatory regions increased by T166A mutation (Fig. [159]5c, d). The mRNA levels of these PPARγ target genes were significantly upregulated in T166A macrophages (Fig. [160]5e). These data demonstrate that PPARγ T166A upregulate lipid synthesis by directly promoting lipid synthesis-related gene expression. Moreover, the expression of genes associated with phenotypic markers and reparative factors was upregulated in M(IL10) or T166A macrophages but was attenuated in macrophages exposed to the FASN inhibitor C75 (Fig. [161]5f). Taken together, these results provide evidence that posttranslational modification of T166 plays a key role in PPARγ-mediated transcriptional activity of lipid synthesis-related genes, leading to the acquisition of a reparative phenotype in macrophages. Fig. 5. Blockage of T166 phosphorylation promotes PPARγ binding to lipid synthesis and accumulation-related gene loci. [162]Fig. 5 [163]Open in a new tab a GO enrichment analysis showing the highly enriched metabolic pathway associated with WT or T166A PPARγ binding sites. b Binding peaks of WT or T166A PPARγ around indicated gene region. c ChIP-qPCR analysis of PPARγ binding to the Dbi, Plin2, Dgat1, Dgat2, Acly, Acaca and Fasn promoters in WT and T166A BMDMs. qPCR was performed with primers specific to the PPARγ-binding motifs. Data were normalized to the input. n = 3 biological replicates per group. d CUT&RUN qPCR of PPARγ binding to the indicated promoters in WT and T166A BMDMs (n = 3 biological replicates per group). Data were normalized to the spike in DNA. e The relative mRNA levels of FAS-related genes in WT or T166A BMDMs (n = 3 biological replicates per group) were assessed by RT-qPCR. f The relative mRNA levels of repair-related genes in WT or T166A BMDMs (n = 3 biological replicates per group) were assessed by RT-qPCR after treatment with IL-10 and C75 for 48 h. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. Intracellular lipids promote growth factor expression by activating STAT3 transcriptional activity We further explored the mechanism by which intracellular lipids could affect macrophage reparative functions. Previous studies have revealed that intracellular lipids may serve as signaling molecules to induce the activation of STAT3^[164]23–[165]25, which directly regulates the expression of growth factors, impacting cellular proliferation, survival, differentiation, and other physiological processes^[166]26,[167]27. Thus, we reasoned that the de novo-synthesized lipids in macrophages would be essential for the increased expression of growth factors mediated by STAT3 activation. We showed that STAT3 was substantially enriched in the regulator regions of Vegfb, Pdgfb, and Tgfb1 in macrophages (Fig. [168]6a). A STAT3 inhibitor, Stattic, substantially inhibited the increase in the expression of growth factors and concomitant reparative macrophage activation induced by PPARγ T166A mutation (Fig. [169]6b, c). However, JAK inhibition did not impact the reparative phenotype of T166A macrophages (Supplementary Fig. [170]12b). These data indicate that PPARγ^A/A macrophage exert reparative function through activating STAT3 in a JAK independent manner. Next, we evaluated the role of lipids in STAT3 activation in reparative macrophages. Free fatty acids (FFA) were employed to elevate intracellular lipid levels of macrophages. FFAs addition in macrophage culture significantly induced STAT3 phosphorylation, which was essential for STAT3 activation (Fig. [171]6d). Luciferase reporter assays showed that FFAs strongly induced STAT3-mediated Vegfb, Pdgfb, and Tgfb1 expression (Fig. [172]6e). CUT&RUN assays revealed that IL-10 strongly enhanced the binding of STAT3 to the regulatory regions of Vegfb, Pdgfb, and Tgfb1, and FAS inhibition using C75 significantly decreased this effect (Fig. [173]6f). Furthermore, STAT3 knockout strongly reduced IL-10-induced expression and secretion of reparative factors (Supplementary Fig. [174]12c–e), and the effects of C75 or FFAs were markedly abolished in macrophages lacking STAT3 (Supplementary Fig. [175]12e, f). These results suggest that intracellular lipids induce reparative factor expressions by promoting STAT3 activation. Besides, we also found that unsaturated fatty acid, such as oleic acid, induced higher levels of STAT3 activation (Fig. [176]6g). Collectively, our data indicate that the increased flux of de novo FAS is essential for STAT3 activation and growth factor expression. Fig. 6. Intracellular lipids are required to support the production of growth factor by macrophages. [177]Fig. 6 [178]Open in a new tab a Binding peaks of STAT3 around indicated gene region. b mRNA expression of the reparative markers Cd163, Mrc1, Tgfb1, Pdgfb and Vegfb in PPARγ^+/+ and PPARγ^A/A BMDMs treated with STAT3 inhibitor Stattic (5 μM). n = 3 biological replicates per group. c Representative flow cytometry histograms and quantifications of surface levels of CD163 and CD206 in PPARγ^+/+ and PPARγ^A/A BMDMs treated with stattic for 48 h (n = 3 biological replicates per group). d The protein expression of phosphorylated (Tyr705) STAT3 and total STAT3 in PPARγ^+/+ and PPARγ^A/A BMDMs in response to free fatty acid (FFA). β-Actin was used as a control. e Transcriptional assays of fatty acid on STAT3-driven gene transcription in dual-luciferase reporter assay. Vegfb, Pgdfb and Tgfb1 promoter was transfected in HEK293T cells with or without STAT3 plasmid in response to FFA. n = 4 biological replicates per group. f CUT&RUN qPCR of STAT3 binding to the indicated promoters in BMDMs treated with C75 in response to IL-10 (n = 3 biological replicates per group). Data were normalized to the spike in DNA. g CUT&RUN qPCR of STAT3 binding to the indicated promoters in BMDMs treated with oleate or palmitate (n = 3 biological replicates per group). Data were normalized to the spike in DNA. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. Fatty acids support the expansion of the endoplasmic reticulum, which is required for the secretion of reparative proteins FAS plays a pivotal role in supplying building blocks for the synthesis of phospholipids (Supplementary Fig. [179]13a). Phospholipids are integral components of biological membranes that support cell proliferation, protein transport and secretion^[180]28. We observed significant alterations in phospholipid metabolic pathways in reparative macrophages using lipidomic analysis (Fig. [181]2b). We also found a notable increase in the mRNA levels of genes related to phospholipid metabolism, such as Chka, Chkb, Cds1, Pcyt1a, Pcyt2 and Plpp2, in IL-10-treated or PPARγ T166A-mutated macrophages (Supplementary Fig. [182]13b). In addition, the expression of phospholipid metabolism genes was significantly attenuated in C75-treated macrophages compared with treated controls (Supplementary Fig. [183]13c). We further addressed the role of phospholipids in membrane system remodeling in reparative macrophage. Although macrophages are typically not proliferative in vitro, there is evidence that macrophages possess the capability of proliferating in vivo in response to tissue microenvironment cues^[184]29. However, we observe no proliferation of macrophages decreased during reparative phase (Supplementary Fig. [185]13d). Thus, we hypothesized that the main role of phospholipids was as essential factors for the endoplasmic reticulum (ER) and Golgi to support the synthesis, transport, and secretion of reparative proteins. Consistent with this hypothesis, we found that ER compartments were considerably expanded in macrophages after stimulation with IL-10 or in T166A macrophages (Supplementary Fig. [186]13e, f). Using [U-^13C] glucose tracer, we observed that IL-10 treatment and T166A mutation strongly upregulated the isotope ratio of ^13C/^12C in ER (Supplementary Fig. [187]13g). Collectively, our results demonstrate that FAS supports the expansion of the ER, which is required for the secretion of reparative proteins through the generation of phospholipids. Topical inhibition of PPARγ T166 phosphorylation accelerates wound healing To apply these findings to clinically translatable outcomes, we tested the effects of local PPARγ T166 phosphorylation inhibition on skin wound healing. We previously found that bardoxolone (CDDO) selectively inhibited PPARγ T166 phosphorylation^[188]30. Thus, we further investigated whether CDDO could affect macrophage activation and metabolic regulation by treating BMDMs with CDDO. The results showed that CDDO-stimulated macrophages exhibited increased lipid accumulation and lipogenesis-related gene expression (Fig. [189]7a–c). More importantly, CDDO induced a reparative phenotype in macrophages (Fig. [190]7c, d). We further investigated the specificity of CDDO in macrophages. The lipid synthesis activity and reparative phenotype induced by CDDO were significantly abolished in PPARγ-knockout macrophages (Supplementary Fig. [191]14a–c). However, CDDO treatment did not further enhance the reparative phenotype and lipid anabolism in PPARγ T166A mutated cells (Supplementary Fig. [192]14d–f). Additionally, although CDDO is known to induce Nrf2 activation, the levels of Nrf2 did not affect CDDO-induced lipid synthesis and reparative phenotype (Supplementary Fig. [193]14g–i). These results reveal that CDDO induces macrophages reparative phenotype through PPARγ T166 dephosphorylation. Fig. 7. CDDO accelerates healing of murine full-thickness excisional wounds. [194]Fig. 7 [195]Open in a new tab a BODIPY staining analysis by flow cytometry in control (CTR) or CDDO-treated macrophages. Dashed line indicates peak of the fluorescent signal. b TG content of BMDMs in response to CDDO (n = 4 biological replicates). c The relative mRNA levels of FAS- and repair-related genes in control or CDDO-treated BMDMs were assessed by RT-qPCR. n = 3 biological replicates per group. d Representative flow cytometry histograms and quantifications of surface levels of CD163 and CD206 in control or CDDO-treated BMDMs (n = 3 biological replicates per group). e Schematic illustration of experimental design. f Representative time-lapse photographs of skin wounds in mice over 12 dpi, treated with vehicle or CDDO ointment. g Wound area after vehicle and CDDO treatment of C57BL/6 mice. (n = 10 per group). h Quantitative analysis of wound re-epithelialization of 6 dpi wound sections from vehicle (n = 3) and CDDO-treated (n = 3) mice. i BODIPY staining analysis by flow cytometry in MACS-sorted wound macrophages at 6 dpi from vehicle and CDDO-treated mice. Dashed line indicates peak of the fluorescent signal. j Representative flow cytometry histograms and quantifications of surface levels of CD163 and CD206 in MACS-sorted wound macrophages at 6 dpi from vehicle and CDDO-treated mice, n = 3 biological replicates per group. Unless specified otherwise, the data are presented as means ± s.e.m. (error bar) and compared using the two-tailed Student’s t test. Source data are provided as a Source Data file. Given that CDDO shaped the reparative property of macrophages in vitro, we created mice with full-thickness excisional wounds and topically applied CDDO ointment or vehicle ointment (Fig. [196]7e). CDDO significantly accelerated wound closure, particularly at 3-6 dpi, during which time macrophages transitioned from a proinflammatory to a reparative phenotype (Fig. [197]7f, g). Histological examination showed that CDDO enhanced the re-epithelialization processes during the repair phase (Fig. [198]7h and Supplementary Fig. [199]15a). Furthermore, we isolated F4/80^+ macrophages from wounds and found increased lipid accumulation and reparative phenotype polarization in CDDO-treated mice compared to vehicle-treated mice (Fig. [200]7i, j). Thus, the effects of local PPARγ T166 phosphorylation inhibition are similar to those of the PPARγ T166A mutation on improving skin wound healing. To further examine the therapeutic effect of CDDO to human wounds, we performed scratch tests using human macrophages. The results showed that CDDO treated human macrophages significantly promoted human skin fibroblasts (HSF) and human epidermal keratinocytes (HaCaT) migration (Supplementary Fig. [201]15b, c). Collectively, these results show that topical PPARγ T166 phosphorylation inhibition could accelerate skin wound healing. Discussion Macrophages are present in virtually all tissues in the body, where they act as pivotal players in signal responses and exhibit multiple biological functions, such as immune response, development, tissue remodeling and wound healing^[202]9,[203]31–[204]33. In the context of tissue repair, macrophages not only secrete an array of proinflammatory or anti-inflammatory cytokines to mediate immune response but also release many matrix metalloproteinases and growth factors^[205]9. These bioactive molecules remodel the extracellular matrix, foster angiogenesis, and ultimately facilitate wound healing. However, the molecular mechanisms that drive and sustain the reparative phenotype of macrophages remain unclear. In this study, we identified PPARγ T166 dephosphorylation-mediated lipid anabolic network as a major regulator of reparative function in macrophages. We detected the activation of FAS in reparative macrophages across 9 tissue injury models, as well as IL-10, TGF-β1 and dexamethasone induced in vitro models. We also revealed that this lipid metabolic reprogramming is conserved in mouse and human reparative macrophages, suggesting targeting this pathway may have therapeutic value. Moreover, intracellular lipids were necessary for sustaining the reparative function of macrophages, partly because fatty acids trigger STAT3-mediated growth factor expression and support the expansion of ER, which is required for the secretion of reparative proteins. Our results provide an additional rationale for the efficacy of FAS activation in the treatment of tissue injury, that offers significant potential in the development of innovative therapeutic strategies. Recent advances in immunometabolism have revealed that macrophages with different activation states have distinct metabolic pathways to meet their demands for energy and the production of specific cell-associated factors^[206]18,[207]34. For instance, M1 macrophages primarily rely on glycolysis for energy production, which enables them to rapidly generate ATP in response to inflammation^[208]35. M2 macrophages exhibit an OXPHOS metabolic phenotype. These cells rely on fatty acid oxidation and mitochondrial respiration for energy production^[209]36. Significant progress has been made in metabolic research on classical activation and alternative activation of macrophages, while there is limited understanding regarding the specific metabolic adaptations that are crucial for wound healing and tissue repair. We demonstrate that injury-associated macrophages undergo distinct metabolic adaptations to facilitate phase-specific repair functions after tissue injury. The most striking observations of this study is that reparative macrophages adjust their metabolism by activating a fatty acid synthetic program which is critical for supporting the growth factor production and secretion. Using scRNA-seq atlas, we emphasized the indispensable role of lipid anabolism in macrophage function, while the catabolic processes, such as glycolysis, OXPHOS and FAO, were strongly impacted by tissue microenvironments. This finding is consistent with previous study showing that both FAO and OXPHOS are dispensable for macrophage activation^[210]37,[211]38. By modulating macrophage metabolism, it is possible to manipulate macrophage phenotypes and functions, thereby impacting disease progression and treatment outcomes^[212]39,[213]40. In cancer immunotherapy, reprogramming the metabolism of tumor-associated macrophages (TAM) can shift them toward an antitumor phenotype, thereby enhancing tumor elimination. Strategies such as inhibiting glycolysis and promoting OXPHOS in TAMs have shown promising results in preclinical models of cancer^[214]41,[215]42. In chronic inflammatory diseases such as atherosclerosis, targeting macrophage lipid metabolism regulates inflammatory responses and cholesterol handling, reducing plaque formation and promoting stability^[216]36,[217]43. These findings highlight the potential of manipulating macrophage metabolism to improve disease treatment outcomes. Therefore, targeting the metabolic regulation of macrophages during the repair phase will allow for more precise modulation of tissue repair and regeneration. In this study, we created an ointment to regulate FAS activity of injury-associated macrophages. This strategy was proven to be useful in promoting skin wound healing. Thus, further exploration of the intricate relationship between macrophage metabolism and tissue repair programs will pave the way for the development of therapies that exploit macrophage plasticity and metabolic rewiring. PPARγ is a peroxisome proliferator-activated receptor (PPAR) family protein that has received particular attention as a critical regulator of a multitude of metabolic processes that maintain cellular energy homeostasis^[218]44. Notably, an extensive amount of literature has addressed the association between PTMs and PPARγ function in the regulation of cell differentiation, tissue development, and energy metabolism^[219]45–[220]47. However, most of the previously published work was completed in the context of adipocytes^[221]48,[222]49. Recent studies have revealed that PPARγ also plays a major role in macrophage differentiation and functions^[223]50–[224]52. Conditional ablation of Pparg in macrophages leads to a pronounced delay in skeletal muscle regeneration^[225]53. Despite this, there remains a paucity of data on PPARγ PTM mechanisms that regulate macrophage phenotypes. It is worth noting that the biological function and regulation mechanism of PPARγ PTMs might be different between adipocytes and macrophages. Phosphorylation of PPARγ at S273 in macrophages regulates the expression of GDF3^[226]54, thereby affecting macrophage reparative functions, while S273 phosphorylation primarily influences insulin sensitivity in adipocytes through the regulation of genes related to glucose and lipid metabolism^[227]55,[228]56. These findings indicate that PPARγ PTMs have cell-specific roles. We previously identified T166 as a phosphorylation site of PPARγ in the DNA-binding domain^[229]30. Here, PPARγ T166 phosphorylation was proven to be a molecular regulator of macrophage activation and metabolism. Genetic blockade of T166 phosphorylation by an alanine mutation enhanced the binding of PPARγ to the regulator regions of lipid synthesis-related genes. The accumulation of lipids may allow reparative phase macrophages to increase anti-inflammatory cytokine and reparative growth factor production by activating STAT3 signaling and directly serving as membrane structures to facilitate protein secretion. In conclusion, this work establishes PPARγ as a required metabolic sensor and transcriptional regulator of skin reparative macrophages, but much more work needs to be done to explore the role of PPARγ in macrophages and other types of immune cells to better elucidate the molecular connection between PPARγ, tissue repair and lipid metabolism. This new insight could lead to the identification of PPARγ T166 as a new pharmaceutical target to facilitate the healthy repair of injured tissues. Methods Generation of PPARγ^T166A mice All mouse procedures and experiments for this study were approved by the Institutional Animal Care and Use Committee at Nanjing University (IACUC-2207009). Animals were housed in a temperature-controlled specific pathogen-free facilities with a 12 h light/dark cycle with 45–65% relative humidity. Food and water were available ad libitum. All mice were more than 8 weeks old at the beginning of the experiments and included both males and females. Animal sample sizes were determined based on previous studies and literatures of the field using similar experimental paradigms. Mice with the PPARγ mutation were generated using CRISPR-Cas9 strategy, as described previously^[230]57. Site-directed mutagenesis was used to alter the codon for T166 to alanine. Single-cell zygotes from superovulated C57BL/6J female mice mated to C57BL/6J males were injected with Cas9 mRNA, sgRNA and ssDNA, and cultured overnight. Transfer the two-cell embryos into oviducts of pseudopregnant Swiss Webster surrogate females. The resulting progeny gave rise to highly chimeric offspring. These chimeras were mated to C57BL/6J mice to obtain germline transmission of the mutant allele. These heterozygous mice (PPARγ^T166A/+) were mated to generate homozygous knock-in mice PPARγ^T166A/T166A. Mice were backcrossed onto C57BL/6J for 6 generations, after which mice heterozygous for the PPARγ^T166A allele (PPARγ^T166A/+) were bred to generate experimental cohorts. The PPARγ^T166A allele was detected by PCR genotyping with primers: 5’-CAAGACATGGTCTCCCTCTG-3’, 5’-AGTTGTCCTAGTGCCTACTG-3’ and 5’-CTGGCCTACACGAAGAAGCTC-3’. Wild-type littermates were used for experimental controls. Animal model of skin injury C57BL/6J mice (Strain No. N000013) were purchased from GemPharmatech Co., Ltd. Before surgery, all mice were anesthetized and dorsal area was shaved. The dorsal skins of mice were secured by inelastic medical adhesive tapes to prevent wound contraction. Subsequently, full-thickness skin excisional wounds (10 mm in diameter) were made using surgical scissors under sterile surgical conditions. To deplete the macrophages in vivo, 100 μl of 5 mg/ml of clodronate liposomes or PBS liposomes (YEASEN) were intradermally injected into wound of each mouse every other day. For local treatment of wounds with CDDO (HY-14909, MCE), either vehicle or CDDO ointment (3 mg/kg) was applied directly to the wound at 0, 3, 6, 9 dpi. The vehicle ointment was formulated by mixing sunflower oil (S304663, aladdin) and Vaseline (V105023, aladdin) in a 1:1 ratio (w/w). For the preparation of the ointment containing the PPARγ T166 inhibitor, CDDO was incorporated into the vehicle ointment with a weight ratio of 1:9. The epithelial space was measured on every 3 days. The area of the wound was calculated using ImageJ software. Then the percentage of wound area closed relative to the original wound was calculated using the following formula: relative closure rate = [1 − (wound size/original wound size)] × 100%. 12 days later, the mice were euthanized in a CO[2] chamber adhering IACUC guidelines. For histologic analysis, wounds were excised at different times after injury. Muscle injury Mice were anaesthetized with isoflurane and 50 μl of 12 μM cardiotoxin (CTX) (RP17303, GenScript) was injected in the tibialis anterior muscle. Muscles were collected for flow cytometry analysis or for histology analysis at different time points after injury. Muscles were isolated and digested with RPMI 1640 medium containing 2 mg/ml Collagenase II (HY-E70005B, MCE), 0.5 mg/ml Collagenase IV (HY-E70005D, MCE) and 4 mg/ml Dispase ([231]S25046, Yuanye Bio-Technology) for 1 h at 37 °C. Digested material was subsequently passed through a 40 μm cell strainer. The resultant cell suspension was washed with PBS and centrifuged at 500 × g for 5 min. Macrophages were isolated by F4/80-positive selection using anti-F4/80 magnetic beads (130-110-443, Miltenyi Biotech), as recommended by the manufacturer. Acute liver injury WT or PPARγ^A/A mice were intraperitoneally injected with a single dose of olive oil (vehicle) or 1% CCl[4] (10 ml/kg, dissolved in olive oil) and sacrificed at 1 day or 4 days post treatment. Livers were perfused with HBSS solution and digested with 0.1 mg/ml DNase I (Roche) and 0.5 mg/ml collagenase IV (HY-E70005D, MCE). The cell suspension was filtered through a 70 μm cell strainer and centrifuged at 30 × g for 3 min to eliminate hepatocytes. The non-parenchymal cell fraction was washed once with PBS at 500 × g for 5 min at 4 °C. Kupffer cells were isolated using anti-F4/80 magnetic beads (130-110-443, Miltenyi Biotech), as recommended by the manufacturer. Blood from WT and PPARγ^A/A mice treated with oil (vehicle) or CCl[4] was collected, and sera were extracted from the blood. The levels of serum ALT and AST were measured with the Alanine aminotransferase Assay Kit (C009-2-1) and Aspartate aminotransferase Assay Kit (C010-2-1) from Jiancheng Bioengineering Institute. Cell culture experiments HEK293T (#GNHu17), RAW 264.7 (#TCM13) were purchased from the cell bank of the Committee on Type Culture Collection of the Chinese Academy of Sciences. HSF and HaCaT cells were generous gifts from Prof. JinHui Wu’s Lab at Nanjing University. DFs were isolated from mouse skin. BMDMs were isolated from mouse bone marrow. Human macrophages were isolated from peripheral blood mononuclear cells. HEK293T, RAW 264.7, DF, HSF and HaCaT cells were cultured in DMEM. BMDMs were maintained in RPMI 1640 medium. Human macrophages were maintained in X-VIVO 15 medium (Lonza). Cells were maintained in culture supplemented with 10% FBS (GIBCO) and 1% penicillin/streptomycin at 37 °C with 5% CO[2] in a humidified incubator. BMDMs generation and in vitro treatment BMDMs were isolated from the femurs and tibiase of C57BL/6J or PPARγ^A/A mice. Femurs and tibiase were placed in PBS. Using surgical scissor cuts off the ends of the bones. Rinse repeatedly with a 1 ml sterile syringe, until most of the cells in the bone marrow blow away. Total bone marrow cells were passed into a 40-μm cell strainer and centrifuged (800 × g for 5 min). Bone marrow cell pellets were resuspended in red blood cell lysis buffer for 5 min and washed in PBS. Next, bone marrow cells were differentiated for 7 days in the presence of M-CSF (50 ng/ml, Z03275, GenScript) in RPMI-1640 medium containing 10% FBS and 100 U/ml penicillin-streptomycin. Total bone marrow cells were seeded in 10-cm culture plates at a density of 5 × 10^6 cells per plate and cultured at 37 °C in a 5% CO[2] incubator. On day 7, macrophage purity of the culture was routinely tested by the expression of F4/80 and CD11b by flow cytometry. Around 95–99% of the cells expressed high levels of F4/80 and CD11b after differentiation. BMDMs were cultured at a density of 5 × 10^5 cells per 12-well in 1 ml complete medium, and stimulated for 48 h with IL-10 (20 ng/ml, Z03161, GenScript), IFNγ (25 ng/ml, Z02916, GenScript), LPS (100 ng/ml, Sigma), IL-4 (20 ng/ml, Z02996, GenScript), TGF-β1 (20 ng/ml, Z03431, GenScript) or dexamethasone (100 nM, Sigma) to achieve an M(IL10), M(LPS + IFNγ),Μ(IL4), Μ(TGF-β1) or M(Dex) polarization state, respectively. In some experiments, BMDMs were treated for 48 h (unless stated otherwise) with C75 (20 μM, HY-12364, MCE), Stattic (10 μM, HY-13818, MCE), rosiglitazone (1 μM, HY-17386, MCE), sodium oleate (100 μM, O7501, Sigma-Aldrich), sodium palmitate (100 μM, P9767, Sigma-Aldrich), FFAs (100 μM, the ratio of sodium oleate and sodium palmitate concentrations at 2:1). Skin-cell isolation Mice were anesthetized with ketamine and killed. The dorsal skin was removed and cut into tiny pieces with scissors, followed by incubation with complete medium (RPMI 1640 and 10% FBS) containing 0.1 mg/ml DNase I (Roche), 0.5 mg/ml hyaluronidase (HY-107910, MCE) and 2 mg/ml collagenase IV (HY-E70005D, MCE) in a shaking water bath (37 °C) for 2 h. The cell suspension was filtered with a 40-μm cell nylon mesh to obtain a single-cell suspension. The remaining red blood cells were lysed using RBC lysis buffer. The acquired single-cell suspensions were processed for flow cytometry or MACS. Magnetic-activated cell sorting The skin digest was resuspended in MACS buffer (PBS, 2 mM EDTA, 0.5% BSA) and incubated with anti-F4/80 microbeads (130-110-443, Miltenyi Biotech) for 15 min at 4 °C. Then, the F4/80^+ fraction was isolated using MACS MS columns according to the manufacturer’s instructions (Miltenyi Biotech). Samples were further processed for immunoblotting, flow cytometry analysis or gene expression analysis as described below. Flow cytometry staining and analysis After collection, skin digest or BMDMs were kept in FACS buffer on ice. Non-specific binding was blocked with FcR blocking reagent (130-092-575, Miltenyi Biotech) for 20 min at 4 °C. Antibodies used for flow cytometry analysis were as follows: anti-CD45 (1:100, 157214, Biolegend), anti-F4/80 (1:100, 17-4801-82, Invitrogen), anti-CD11b (1:100, 12-0112-81, Invitrogen), anti-Gr1(1:100, 552093, BD Biosciences), anti-CD3e (1:100, 12-0031-83, Invitrogen), anti-CD206 (1:100, 141705, Biolegend), anti-CD163 (1:100, 17-1631-82, Invitrogen), anti-CD80 (1:100, 12-0801-81, Invitrogen), anti-CD86 (1:100, 553692, BD Biosciences), anti-MHC class II (1:100, 12-5321-82, Invitrogen). For intracellular staining, cells were fixed and permeabilized using the Intracellular Fixation & Permeabilization Buffer Set (#88-8824-00, eBioscience) and intracellularly stained for 1 h at 4 °C with anti-CD206 or anti-CD163 antibody. The lipid content was assessed with 2 μM BODIPY 493/503 (D3922, Invitrogen), following 20 min of incubation at 37 °C. Flow cytometric analysis was performed on a NovoCyte Flow Cytometer (ACEA Biosciences Inc). The gating strategies for Fig. [232]4e and Supplementary Fig. [233]10a are presented in Supplementary Fig. [234]16. Cell fractionation and immunoblot analysis The nuclear and cytoplasmic proteins were extracted according to the manufacturer’s instructions of the kit (P0027, Beyotime). After collection, cells were washed 3 times with PBS, and lysed in RIPA buffer supplemented with protease and phosphatase inhibitor cocktail. Protein concentration was quantified by BCA reagent. Next, 50 μg of protein was then separated by electrophoresis using SDS-PAGE gels and then transferred to PVDF membranes. Membranes were blocked for 1 h in 5% BSA at room temperature. Membranes were incubated overnight with the following primary antibodies: mouse anti-PPARγ (1:200, sc-7273, Santa Cruz), rabbit anti-PPARγ (1:1000, #2443S, Cell Signalling Technology), anti-PPARγ p-Ser273 (1:1000, bs-4888R, Bioss), anti-PPARγ p-Ser112 (1:500, ab195925, Abcam), p-STAT3 (Y705) (1:2000, #9145, Cell Signalling Technology), anti-STAT3 (1:1000, BS1336, Bioworld), anti-ACC1 (1:4000, 21923-1-AP, proteintech), anti-FASN (1:5000, 10624-2-AP, proteintech), Lamin A/C (1:1000, BS1446, Bioworld), anti-β-actin (1:10,000, 700068, Zen-Bioscience) and GAPDH (1:5000, 390035, Zen-Bioscience). Mouse monoclonal antibody against p-T166 (1:500) was generated as described previously^[235]30. Signals were detected by a chemiluminescence System (Tanon-5200S) using Enhanced ECL Chemiluminescent Substrate (36222, YEASEN). The intensities of the bands were quantified by ImageJ software. The expression of proteins was normalized to a housekeeping protein (β-actin), and the phosphorylation status was determined by normalizing to a respective total protein. Quantitative real-time PCR Total RNA was isolated from cells or skin tissue samples with TRIzol reagent (R401-01, Vazyme). Reverse transcription was performed with HiScript III RT SuperMix (R323-01, Vazyme) kit. Quantitative RT-PCR assays were performed on a real-time PCR System (C1000, Bio-Rad) using qPCR SYBR Green Master Mix (Q111-02, Vazyme). The expression of genes of interest was normalized to housekeeping gene Rplp0, and data were normalized to the average of the control group. Primer sequences are described in Supplementary Data [236]1. Single-cell transcriptomic atlas of wound macrophages We obtained scRNA-seq datasets containing wound macrophages from 9 mouse wound healing models, including an anterior cruciate ligament injury (ACLI) model using a noninvasive single dynamic tibial compressive overload method, a polidocanol-mediated airway injury (AWI) model, a kidney injury (KDI) model by short ischemic reperfusion, an acetaminophen-induced acute liver damage (LIVER) model, a bleomycin-mediated acute lung injury (LUNG) model, a myocardial infarction (MCIN) model by left anterior descending coronary artery ligation, a spinal cord injury (SCI) model through laminectomy, a full-thickness excisional skin injury (SKIN) model, and a BaCl[2]-induced skeletal muscle injury (SMI) model (Supplementary Data [237]2). Samples were categorized into inflammatory phase (phase I) and reparative phase (phase II) according to the description provided in the source articles related to these datasets. The scRNA-seq data analysis was performed using Seurat (v4.3.0) R package^[238]58. In brief, quality control was applied to cells based on the total UMI counts and number of detected genes. Cells with less than 200 detected genes were filtered. To remove potential doublets, cells with UMI count above 40,000 and detected genes above 6000 are filtered out. The gene counts were normalized using SCTransform. For clustering, top 2000 variable genes were used, and ribosomal and mitochondrial genes were removed from the variable gene set before running PCA and identifying clusters using a shared nearest neighbor (SNN) modularity optimization-based clustering algorithm. An annotation tool, SingleR (v2.2.0), was employed to assist cluster identification^[239]59. The gene list provided by the source article were also considered as marker genes in macrophage identification (Supplementary Data [240]2). Subclusters which were identified by macrophages were isolated from the original datasets and re-analyzed using above method. The ALRA method^[241]60 was used for imputing dropout genes when necessary. Trajectory analysis Trajectory analysis was performed by monocle (v2.30.1) R package. Monocle introduced the strategy of ordering RNA-sequencing samples in pseudotime, placing them along a trajectory according to their transcriptional similarity. To visualize the macrophage phenotype transition along the cell trajectory, dataset ([242]GSE131364) containing M0, M(IFNγ), M(LPS), M(IL4) and M(IL10) samples was employed as the reference dataset. To train the SingleR classifier, simulated samples were synthesized from the reference dataset by randomly integrating a certain percentage of samples from the same group. The final dataset comprised of 20 samples within each group. Then, the reference data were wrapped into a SummarizedExperiment object for SingleR classification. SCENIC analysis Activated regulons in SKIN macrophages were analyzed using SCENIC (v1.3.1) R package^[243]61. Raw count matrix containing 2953 cells and 15,563 genes was used as input. In brief, the co-expression network was calculated by GRNBoost2 and the regulons were identified by RcisTarget. Next, the regulon activity for each cell was scored by AUCell. The BH method was used to correct multiple hypothesis. RNA sequencing and data analysis Full-thickness excisional skin wounds were generated on the PPARγ^+/+ and PPARγ^A/A mice. The mice were sacrificed and wound tissues were digested into single-cell suspensions at 6 dpi. Then, skin macrophages were isolated using MACS. BMDMs were isolated from PPARγ^+/+ and PPARγ^A/A mice and cultured at a density of 5 × 10^6 cells/10-cm dish. Total RNA was extracted from BMDMs and skin macrophages using TRIzol reagent according to the manufacturer’s protocol. Next generation sequencing library preparations were conducted according to the manufacturer’s protocol. The libraries were sequenced on an Illumina Novaseq 6000 Sequencer. Read quality was assessed with FastQC (v0.12.1) and adapters trimmed using Cutadapt (v4.4)^[244]62. Reads were then aligned to the mm10 genome (UCSC) using HISAT2 (v2.2.1)^[245]63 and reads in exons were counted using featureCounts function in Subread (v2.0.3) package^[246]64. Analysis of differentially expressed genes was conducted using DESeq2 (v1.40.1) R package^[247]65. Count data was transformed to log2 scale and normalized with respect to library size using rlog transformation produces. Absolute log2 fold change (log[2]FC) > 1.5 and p value < 0.05 was considered statistically significant. Enrichment analysis Gene set enrichment analysis (GSEA) was conducted and visualized using fgsea (v1.22.0) R package^[248]66. Single sample gene set enrichment analysis (ssGSEA) was performed using GSVA (v1.44.2) R package^[249]67 with custom gene lists. GO categories (biological process) enrichment analysis was performed and visualized using clusterProfiler (v4.8.3) R package^[250]68. The p values were adjusted using the FDR by BH method. Gene list of six metabolic pathway was obtained from Kyoto Encyclopedia of Genes and Genomes (KEGG) database, including glycolysis (mmu00010), TCA (mmu00020), OXPHOS (mmu00030), PPP (mmu00061), FAO (mmu00071), and FAS (mmu00061). Chromatin immunoprecipitation (ChIP) assays ChIP experiments were carried out using the Enzymatic Chromatin IP Kit (#9003, Cell Signalling Technology) according to the manufacturer’s instructions. In brief, BMDMs were isolated from PPARγ^+/+ and PPARγ^A/A mice and cultured at a density of 5 × 10^6 cells per 10-cm culture dish as described above. 1 × 10^7 cells were fixed with 1% formaldehyde in culture medium at room temperature for 10 min. Subsequently, the crosslinked cells were collected, lysed and digested with MNase. The sheared chromatin was subjected to immunoprecipitation with a rabbit anti-PPARγ antibody (1:100, #2443S, Cell Signalling Technology) or rabbit IgG control overnight at 4 °C with constant rotation. The isolated complexes were collected with Protein G magnetic beads and eluted by incubating at 65 °C in a buffer with high salt concentration. The qPCR following ChIP was performed with the primer sequences listed in Supplementary Data [251]1. For ChIP-seq, libraries were prepared from purified immunoprecipitation DNA by blunting, A-tailing, adapter ligation. Library quality was assessed on the Agilent Bioanalyzer 2100 system. The ChIP-seq libraries were sequenced on an Illumina Novaseq 6000 Sequencer. Then, sequenced reads were aligned to reference mouse genome (mm10 assembly, UCSC) using Bowtie2 (v2.5.1)^[252]69, and clonal reads were removed from further analysis. After mapping reads to the reference genome, we used the MACS2 (v2.1.0)^[253]70 peak calling software to identify regions of IP enrichment over background. p value threshold of 0.0005 was used for all datasets. ChIPseeker (v1.36.0) R package^[254]71 was used to retrieve the nearest genes around the peak and annotate genomic region of the peak. Peak-related genes were used to perform enrichment analysis to identify the function enrichment results. The binding peaks were visualized using Gviz (v1.46.1) R package. CUT&RUN assays CUT&RUN assays were carried out using the Hyperactive pG-MNase CUT&RUN Assay Kit (HD101-01, Vazyme) according to the manufacturer’s instructions. 500,000 BMDMs from C57BL/6J mice were harvested by centrifugation and washed in wash buffer. Concanavalin A-coated beads are activated and stabilized by Ca^2+ and Mn^2+ cations. Cells were bound to activated Concanavalin A-coated magnetic beads, then permeabilized with Wash buffer containing 0.05% Digitonin. The bead-cell slurry was incubated with antibody buffer in a 100 μl volume at 4 °C overnight on a rotator. anti-PPARγ (1:100, #2443S, Cell Signalling Technology) or anti-STAT3 (1:100, BS1336, Bioworld) antibodies were used to IP PPARγ or STAT3 in BMDMs. After 2–3 washes in 1 ml Dig-wash, beads were resuspended in 100 μl pG/MNase and incubated for 1 h at 4 °C. After two washes in Dig-wash, beads were resuspended in ice-cold calcium incubation buffer. Tubes were replaced on ice for 1.5 h, followed by immediate addition of EGTA-STOP buffer. Beads were incubated at 37 °C for 30 min, replaced on a magnet stand and the liquid was removed to a fresh tube and DNA was extracted. The qPCR following CUT&RUN was performed with the primer sequences listed in Supplementary Data [255]1. For CUT&RUN-seq, libraries were prepared from purified immunoprecipitation DNA by blunting, A-tailing, adapter ligation. Library quality was assessed on the Agilent Bioanalyzer 2100 system. The libraries were sequenced on an Illumina Novaseq 6000 Sequencer. Short reads (<20 bp) and adapter sequences were removed using Cutadapt (v4.4). Trimmed reads were aligned to the reference mouse genome (mm10 assembly, UCSC) using Bowtie2 (v2.5.1) with the “--very-sensitive --no-mixed --no-discordant” parameters. Then, the reads were filtered, retaining only the mapped read pairs that are located on the same chromosome and have a fragment length of less than 1000 bp. Subsequently, the data were converted into bedgraph file format using Bedtools (v2.31.0)^[256]72. The bedgraph files were visualized using Gviz (v1.46.1) R package. Luciferase reporter assay Genes encoding mouse PPARγ1 (WT or T166A) were cloned separately into the pcDNA3.1(-) vector. PPRE X3-TK-luc was a gift from Bruce Spiegelman (Addgene plasmid # 1015). HEK293T cells were transfected with PPARγ (WT or T166A), PPRE-Luc, and control Renilla plasmid pRL at a molar ratio of 50:50:1. Following an 8 h transfection period, the culture medium was replaced with fresh culture medium supplemented with 10% FBS. Subsequently, small molecules were added to the culture wells and incubated for a duration of 24 h. Cells were harvested and lysed in luciferase lysis buffer, and the luciferase activity of the cells was carried out following the dual-luciferase reporter assay system (Promega, E1910). Firefly luciferase activity obtained from each sample was normalized to the Renilla luciferase activity from the same sample. For determination of STAT3 binding region on Vegfb, Pdgfb and Tgfb1 promoters, there promoters were amplified and cloned into pGL6 reporter plasmid (D2105, Beyotime). The constructed pGL6 vector was transfected into HEK293T cells with/without STAT3 expression plasmid and luciferase activity was measured and analyzed as above. Tracing with [U-^13C] glucose BMDMs were left unstimulated or stimulated for 48 h with IL-10. Medium was changed to RPMI 1640 containing [U-^13C] glucose (2 g/l) or ^12C-glucose (2 g/l), and 48 h later cell culture plates were washed with PBS and snap-frozen in liquid nitrogen and subjected to fatty acid flux analysis (LipidALL Technologies). Lipids were extracted from BMDMs using 1 ml of ice-cold methanol. Samples were incubated for 30 min at 4 °C and centrifuged for 10 min at 14,000 × g. Supernatant was dried in a SpeedVac under OH mode. Total protein content was determined from the dried pellet using the BCA protein assay. The dried extract was reconstituted in LCMS grade methanol for fatty acid flux analysis on a system comprising an Agilent 1290 II UPLC coupled to Sciex 5600+ quadrupole-TOF MS. Fatty acids were separated on a Waters ACQUITY HSS-T3 column (3.0 × 100 mm, 1.8 μm). Information-dependent acquisition mode was used for MS/MS analyses of the metabolites. Collision energy was set at −30 ± 15 eV. Data acquisition and processing were performed using Analyst TF 1.7.1 Software (AB Sciex). All detected ions were extracted using MarkerView 1.3 (AB Sciex). PeakView 2.2 (AB Sciex) was applied to extract MS/MS data and perform comparisons with the Metabolites database (AB Sciex), HMDB and standard references to