Foxp1 is critical for the maintenance of regulatory T-cell homeostasis and suppressive function.

Regulatory T (Treg) cells play central roles in maintaining immune homeostasis and self-tolerance. However, the molecular mechanisms underlying Treg cell homeostasis and suppressive function are still not fully understood. Here, we report that the deletion of another P subfamily members of the forkhead box (Foxp) subfamily member Foxp1 in Treg cells led to increased numbers of activated Treg (aTreg) cells at the expense of quiescent Treg cells, and also resulted in impaired Treg suppressive function. Mice with Foxp1-deficient Treg cells developed spontaneous inflammatory disease with age; they also had more severe inflammatory disease in colitis and experimental autoimmune encephalomyelitis (EAE) models. Mechanistically, we found that Foxp1 bound to the conserved noncoding sequence 2 (CNS2) element of the Foxp3 locus and helped maintain Treg suppressive function by stabilizing the Foxp3 expression. Furthermore, we found that Foxp1 and Foxp3 coordinated the regulation of cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) expression levels. Taken together, our study demonstrates that Foxp1 plays critical roles in both maintaining Treg cell quiescence during homeostasis and regulating Treg suppressive function.

Introduction

Regulatory T (Treg) cells are a subset of CD4+ T cells of vital importance in the maintenance of immunological self-tolerance and T-cell homeostasis [1,2]. Like conventional CD4+ T cells, Treg cells in the periphery have subpopulations with different functions. Based on the expression levels of cell surface markers CD44 and CD62L, Treg cells can be divided into CD44lowCD62Lhigh resting Treg (rTreg) and CD44highCD62Llow activated Treg (aTreg) cells [3,4]. rTreg cells are more quiescent Treg cells in the secondary lymphoid tissues, where paracrine interleukin 2 (IL-2) is required for their homeostasis, and rTreg cells have less spontaneous proliferation [3]. In contrast, aTreg cells are also present in the secondary lymphoid tissues, but they are the predominant Treg cell population in the nonlymphoid tissues [4,5]. Compared with rTreg cells, aTreg cells appear to be more potent suppressors in anti-inflammatory responses [5]. Studies have shown that transcription factors, including forkhead box O1 (Foxo1), myeloblastosis oncogene (Myb), interferon regulatory factor 4 (IRF4), the c-Rel subunit of activation of nuclear factor κB (NF-κB c-Rel), are important for aTreg differentiation and migration [4,6-8]. The transcriptional regulation of rTreg and aTreg cell homeostasis, however, is poorly understood.

Transcription factor Foxp3 has been considered the master regulator of Treg cells, indispensable for Treg cell development, maintenance, and suppressive function [9]. Loss-of-function mutations of Foxp3 are sufficient to drive lethal systemic autoimmunity both in mice and in humans due to defects in Treg cell development [10-12]. The deletion of Foxp3 in mature Treg cells leads to abnormal target gene expression and their impaired suppressive function [13]. It has also been reported that attenuation or unstable Foxp3 expression impairs the function of Treg cells [14-16], suggesting that the stringent regulation of Foxp3 expression levels is important for the maintenance of Treg cell function and identity. Studies have shown that transcription factors, such as GATA binding protein 3 (Gata3), E26 avian leukemia oncogene 1, 5′ domain (Ets-1), cyclic-AMP response element binding protein (CREB), forkhead box O (Foxo), SMAD family member 3 (Smad3), nuclear factor of activated T cells (NFAT), signal transducer and activator of transcription 5 (STAT5), runt related transcription factor 1 (Runx1), reticuloendotheliosis oncogene (c-Rel), and special AT-rich sequence binding protein 1 (Satb1), have the ability to transactivate Foxp3 expression by directly binding to the conserved noncoding sequences (CNSs) of the Foxp3 locus [17-26]. CNSs include a promoter, CNS1 and CNS2 in the first intron, CNS3 downstream of exon 1, and CNS0 in an intron of the neighboring gene 5′ of the Foxp3 locus [26,27]. The functions of CNSs have been intensively studied: CNS1 is important for the generation of Treg cells generated in the periphery (pTreg); CNS3 acts as a pioneer element, indispensable for both the differentiation of thymus-derived Treg cells (tTreg) and pTreg; and CNS2 does not affect thymic generation of Treg cells, but rather it helps maintain the stable expression of Foxp3 and protects Treg cell identity under destabilizing conditions [15,16,27]. Despite all this work, the transcriptional regulation of Foxp3 in Treg cells is not fully understood.

Studies have suggested that Foxp3 alone is not sufficient to determine the Treg cell differentiation; instead, Treg differentiation and functional maintenance need a well-orchestrated transcriptional network, including the interactions between Foxp3 and its partners [9]. A large number of Foxp3 partners have been identified as transcription-related proteins, including transcription factors: nuclear factor of activated T cells, cytoplasmic, calcineurin-dependent 2 (NFATc2), Runx1, B cell leukemia/lymphoma 11B (Bcl11b), forkhead box P1 (Foxp1), forkhead box P4 (Foxp4), Gata3, signal transducer and activator of transcription 3 (STAT3), IKAROS family zinc finger 1 (Ikaros/Ikzf1), IKAROS family zinc finger 3 (Aiolos/Ikzf3), E26 avian leukemia oncogene (Ets), and CCR4-NOT transcription complex, subunit 3 (Cnot3) [28]. Among these identified partners, Foxp1 is one of the P subfamily members of the forkhead box (FOX) transcription factor family, as is Foxp3. Our studies have revealed that Foxp1 is vital for both T-cell and B cell development and function [29-32]. In T cells, Foxp1 actively maintains naive T-cell quiescence by restraining IL-7Rα expression levels [31]. Furthermore, Foxp1 negatively regulates the differentiation and function of follicular helper T cells (TFH) [32]. Although it has been reported that Foxp1 can form a heterodimer with Foxp3 by the conserved leucine zipper domain [28,33], whether Foxp1 plays any role in Treg cell development and function is not known.

Here, we report that Foxp1 regulates Treg cell homeostasis and suppressive function. Loss of Foxp1 specifically in the Treg cells results in increased numbers of aTreg cells at the expense of rTreg cells. The mice with Treg-specific depletion of Foxp1 develop spontaneous inflammatory disease and are more susceptible to dextran sulfate sodium (DSS)-induced colitis and experimental autoimmune encephalomyelitis (EAE). Mechanistically, we found that Foxp1 binds to the CNS2 element of the Foxp3 locus and helps maintain the stable expression of Foxp3. The restoration of Foxp3 expression partially rescues the suppressive function of Foxp1-deficient Treg cells. Moreover, we found that Foxp1 coordinates with Foxp3 in regulating cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) expression levels. Our findings demonstrate that Foxp1 is critical for Treg cell homeostasis, stability, and suppressive function.

Results

Foxp1 is differentially expressed in rTreg and aTreg cells

Previously, we have shown that Foxp1 has no effect on the generation of Treg cells in the thymus of Foxp1CD4 mice [30], in which Foxp1 is deleted at the double-positive (DP) thymocyte stage. Nevertheless, in Foxp1CD4 mice, Foxp1 deletion affects the activation and homeostasis of conventional CD4+ T cells [30], and it is not clear how this may obscure the impact of Foxp1 on Treg cells. To better define the role of Foxp1 in Treg cells, we crossed mice bearing loxP-flanked Foxp1 alleles (Foxp1) [30,31] with Foxp3-YFP-Cre knock-in mice (Foxp3) [34] to deplete Foxp1 specifically in Treg cells (Foxp1Foxp3).

To understand Foxp1 functions in Treg cells, we first examined the Foxp1 expression in wild-type (WT) and Foxp1Foxp3 Treg cells. Given that Foxp1 is differentially expressed in conventional naive CD4+ T cells and activated CD4+ T cells [32], we reason that Foxp1 may also have different expression patterns in CD44lowCD62Lhigh rTreg and CD44highCD62Llow aTreg (). We found that rTreg cells expressed high levels of the full-length Foxp1 isoform, Foxp1A, and low levels of the short isoform Foxp1D (), which resembles the Foxp1 expression in conventional naive CD4+ T cells, suggesting that Foxp1 may also play an important role in maintaining the quiescent state of rTreg cells, like its role in naive CD4+ T cells [31]. Surprisingly, however, we found that compared with rTreg cells, aTreg cells had lower Foxp1A expression at both the mRNA level and the protein level (). This is quite different from the Foxp1 expression pattern in activated conventional CD4+ T cells, in which Foxp1A is expressed at a high level constitutively, and Foxp1D is greatly induced [32]. And as expected, Foxp1 was completely deleted in Treg cells from Foxp1Foxp3 mice (). We further examined Foxp1 expression in Treg cells across various tissues, including lung, liver, and small intestine, in WT mice. Studies have shown a high proportion of aTreg cells in the small intestine [3], and, consistently, we found that among all the tissues examined, the Foxp1 expression levels were the lowest in Treg cells in the small intestine (). These results together suggest that the amount of Foxp1 protein may be important for Treg cell differentiation or function.

Foxp1 is differentially expressed in rTreg and aTreg cells.

(A) Immunoblot analysis of Foxp1 expression in sorted YFP−CD4+CD44lowCD62Lhigh naive CD4+ T cells (WT naive), YFP+CD4+CD44lowCD62Lhigh rTreg cells (WT rTreg), YFP+CD4+CD44highCD62Llow aTreg cells (WT aTreg) from Foxp3 mice, and YFP+ Treg cells from Foxp1Foxp3 mice. (B) Relative Foxp1 mRNA expression in sorted WT naive CD4+ T cells, rTreg cells, aTreg cells from Foxp3 mice, and YFP+ Treg cells from Foxp1Foxp3 mice, n = 3–4. (C) Intracellular staining of Foxp1 in splenic WT naive CD4+ T cells, rTreg cells, aTreg cells in Foxp3 mice and splenic Treg cells in Foxp1Foxp3 mice (left panel), and the corresponding mean fluorescence intensity (MFI) of Foxp1 (right panel) (n = 2). Data in (A, C) represent at least two independent experiments. Data in (B) represent four independent experiments. Data in (B, C: right panel) are mean ± SEM. *P < 0.05 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; CD, cluster of differentiation; Foxp1, forkhead box P1; Foxp1A, Foxp1 isoform A; Foxp1D, Foxp1 isoform D; MFI, mean fluorescence intensity; rTreg, resting Treg; WT, wild-type; YFP, yellow fluorescent protein.

Foxp1 enforces the quiescent state of Treg cells and inhibits aTreg cell differentiation

To address whether Foxp1 is important for Treg cell homeostasis, we first analyzed the frequencies and cell numbers of Treg cells in the thymuses and the peripheral lymphoid and nonlymphoid organs in Foxp1Foxp3 mice and control Foxp3 mice. We found that the frequencies of Treg cells in CD4+ T cells in a majority of the organs examined were very similar between Foxp3 and Foxp1Foxp3 mice, except those in the spleen, the liver, and the lung, which were higher in Foxp1Foxp3 mice (). We found that the numbers of Treg cells in the spleens and the lymph nodes, but not the thymuses, of the Foxp1Foxp3 mice were significantly increased (), suggesting that Foxp1 deficiency has no effect on the generation or maintenance of Treg cells in the thymus, but it impairs the Treg cell homeostasis in the periphery.

Foxp1 enforces the quiescent state of Treg cells and inhibits aTreg cell differentiation.

(A) Flow cytometry analysis of YFP+ Treg cells in CD4+ T cells in indicated organs from 6–8-week-old sex- and age-matched Foxp3 and Foxp1Foxp3 mice. Numbers adjacent to the outlined areas indicate the percentages of YFP+ Treg cells. (B) The frequencies of YFP+ Treg cells in CD4+ T cells in indicated organs from mice as in (A) (n = 3–7). (C) Quantification of total Treg cells in the thymuses, the spleens, and the lymph nodes from the mice as in (A) (n = 3–6). (D) Flow cytometry analysis of CD44 and CD62L expression in Treg cells from 6–8-week-old sex- and age-matched Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined areas indicate the percentages of rTreg cells (CD44lowCD62Lhigh) or aTreg cells (CD44highCD62Llow). (E) Kinetics of frequencies and numbers of rTreg and aTreg cells in Foxp3 and Foxp1Foxp3 mice at different ages (n = 3–6). (F) Flow cytometry analysis of Ki-67 expression in rTreg and aTreg cells in the spleens from 8-week-old Foxp3 and Foxp1Foxp3 mice; numbers above bracketed lines indicate percentages of Ki-67+ Treg cells. (G) Proliferation of rTreg cells stimulated by different concentrations of plate-bound anti-CD3/CD28 antibodies in the presence of 200 U/mL IL-2 for 3 days, as shown by the dilution of CellTrace; numbers above bracketed lines indicate percentages of Treg cells with two or more divisions. Data in (A-E) are representative of at least three independent experiments. Data in (F-G) are representative of two independent experiments. Data in (B, C, E) are mean ± SEM, *P < 0.05, ** P < 0.01, *** P < 0.001 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; CD, cluster of differentiation; Foxp1, forkhead box P1; FSC-H, forward scatter-height; IL, interleukin; Ki-67, antigen identified by monoclonal antibody Ki 67; rTreg, resting Treg; SI, small intestine; Treg, regulatory T; YFP, yellow fluorescent protein.

Because Foxp1 is a key transcription factor that actively maintains conventional naive T-cell quiescence [31], we reasoned that Foxp1 might play important roles in maintaining rTreg cell quiescence and restricting aTreg cell differentiation. To address this question, we examined the frequencies and cell numbers of CD44lowCD62Lhigh rTreg and CD44high CD62Llow aTreg cells in Foxp3 and Foxp1Foxp3 mice at different ages. We found that loss of Foxp1 led to significantly decreased percentages and numbers of rTreg cells, but with elevated percentages and numbers of aTreg cells in both the spleens and the lymph nodes of Foxp1Foxp3 mice (). Such changes were observed in Foxp1Foxp3 mice as early as when they were 3 weeks old (). When the mice were older than 16 weeks, Foxp1Foxp3 mice almost lost the rTreg cell population and compensatorily had many more aTreg cells ().

One of the pronounced characteristics of nonquiescent T cells is their enhanced proliferation capacity. To further determine whether Foxp1-deficient Treg cells lose quiescence, we examined the cell proliferation of these cells. The antigen identified by monoclonal antibody Ki 67 (Ki-67) protein is a cellular marker strictly associated with cell proliferation [35]. We found that in 3-week-old () and 8-week-old mice (), the percentages of Ki-67+ cells in both rTreg and aTreg cells from Foxp1Foxp3 mice were increased compared with their WT counterparts, suggesting that Foxp1-deficient rTreg and aTreg cells proliferated more in vivo. We also sorted out CD44lowCD62Lhigh rTreg cells and CD44highCD62Llow aTreg cells from Foxp3 and Foxp1Foxp3 mice and labeled them with CellTrace, which tracks cell proliferation. We found that upon stimulation by different concentrations of anti-CD3/CD28 antibodies in vitro, Foxp1-deficient rTreg cells proliferated more than WT rTreg cells did (), although we did not observe obvious differences in cell proliferation between WT and Foxp1-deficient aTreg cells (). Of note, the apoptosis of rTreg and aTreg cells in Foxp3 and Foxp1Foxp3 mice showed no difference (), suggesting that the higher numbers of Treg cells in Foxp1Foxp3 mice result from the increased Treg cell proliferation. Taken together, our results suggest that Foxp1 plays a critical role in maintaining the rTreg cell pool and negatively regulates aTreg cell differentiation.

Foxp1 regulates rTreg cell quiescence in a cell-intrinsic manner

To address whether Foxp1 regulates rTreg cell quiescence in an intrinsic manner, we took advantage of Foxp1Foxp3 female mice, theoretically, in which around half of the Treg cells are Foxp1 sufficient, while the other half of the Treg cells are Foxp1 deficient due to the random X chromosome inactivation. Consistent with the results in Foxp1Foxp3 male mice, we found that the percentage of Ki-67+ Treg cells in Foxp1Foxp3 female mice was also higher than that in control Foxp3 female mice (, indicating that Foxp1 intrinsically regulates Treg cell proliferation. However, in Foxp1Foxp3 females, we found that Foxp1-deficient Treg cells had compromised competitive fitness in the spleens or the lymph nodes (, with aTreg cells being almost outcompeted by WT Treg cells (. Therefore, Foxp1Foxp3 female mice are not suitable to address whether Foxp1 intrinsically controls the balance between rTreg and aTreg. To resolve this issue, we used Foxp1Cre-ERT2 mice, in which Treg cells have normal development and function before tamoxifen treatment. We treated mice with tamoxifen and analyzed the phenotype at day (d)8. At this time point, the conventional T cells in Foxp1Cre-ERT2 mice had no obvious phenotypic alteration based on the expression of CD44 and CD62L (). However, there were more aTreg cells within the Treg cell population in Foxp1Cre-ERT2 mice (. And consistently, acute deletion of Foxp1 led to increased Ki-67+ Treg cells in Foxp1Cre-ERT2 mice (. Taken together, our data suggest that Foxp1 enforces the quiescent state of Treg cells in a cell-intrinsic manner.

Foxp1 regulates rTreg cell quiescence in a cell-intrinsic manner.

(A) Flow cytometry analysis of Ki-67 expression in Foxp3+YFP+ rTreg and aTreg cells in the spleens from Foxp3 and Foxp1Foxp3 female mice; numbers above bracketed lines indicate the percentages of Ki-67+ Treg cells. (B) Flow cytometry analysis of Foxp3+YFP+ Treg cells from spleens and lymph nodes of Foxp3 and Foxp1Foxp3 female mice (left panel), and frequencies of Foxp3+YFP+ Treg cells (right panel) (n = 5). (C) Frequencies of Foxp3+YFP+ aTreg cells in CD4+ T cells, n = 5. (D-F) Foxp1 and Foxp1Cre-ERT2 mice were treated with tamoxifen, and phenotypes were analyzed at d8 after treatment. (D) Flow cytometry analysis of CD44 and CD62L expression in CD4+Foxp3− conventional T cells from Foxp1 and Foxp1Cre-ERT2 mice; numbers adjacent to the outlined areas indicate the percentages of naive cells (CD44lowCD62Lhigh) or activated cells (CD44highCD62Llow). (E) Flow cytometry analysis of CD44 and CD62L expression in Treg cells from Foxp1 and Foxp1Cre-ERT2 mice; numbers adjacent to the outlined areas indicate the percentages of rTreg cells (CD44lowCD62Lhigh) or aTreg cells (CD44highCD62Llow). (F) Flow cytometry analysis of Ki-67 expression in rTreg and aTreg cells from Foxp1 and Foxp1Cre-ERT2 mice. Numbers above bracketed lines indicate percentages of Ki-67+ Treg cells. Data in (A-F) are representative of at least two independent experiments. Data in (B: right panel, C) are mean ± SEM, *P < 0.05 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; CD, cluster of differentiation; d, day; Foxp1, forkhead box P1; Ki-67, antigen identified by monoclonal antibody Ki 67; ns, no significance; rTreg, resting Treg; Treg, regulatory T; YFP, yellow fluorescent protein.

Foxp1 controls Treg suppressive function

We next asked whether Foxp1 plays any role in Treg suppressive function. We first performed phenotype analysis in both 6–8-week-old and >6-month-old mice. Compared with control Foxp3 mice, we found that Foxp1Foxp3 mice had splenomegaly and lymphadenopathy at 8 weeks of age (). Flow cytometry analysis showed that the conventional T cells in 8-week-old Foxp1Foxp3 mice exhibited the memory/effector phenotype (CD44highCD62Lhigh and CD44highCD62Llow) (), with increased percentages of interferon gamma (IFNγ)-producing T cells in the spleens and nonlymphoid organs (). Foxp1Foxp3 mice older than 6 months had even higher percentages of CD44highCD62Llow activated T cells, and more IFNγ or IL-4–producing T cells (). And strikingly, those older Foxp1Foxp3 mice developed a multi-organ inflammatory disease (), strongly suggesting that Foxp1 plays an important role in maintaining the suppressive function of Treg cells under the homeostatic condition.

Foxp1 regulates Treg suppressive function.

(A) Flow cytometry analysis of CD44 and CD62L expression in conventional CD4+ T cells (YFP−CD4+) and CD8+ T cells in the spleens of 28-week-old Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined areas represent the percentages of gated cells. (B) Intracellular staining of cytokines in conventional CD4+ and CD8+ T cells in the spleens of 28-week-old Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined areas represent the percentages of gated cells. (C) Hematoxylin and eosin staining of salivary gland, pancreas, lung, small intestine, and colon sections. The magnification is ×100. Black arrows indicate the areas of immune cell infiltration. Scale bars, 100 μm. (D-G) Foxp3 and Foxp1Foxp3 mice (n = 6) were treated with 2.5% DSS, and the severity of colitis in mice was evaluated by loss of body weight (D), survival curve at indicated time points (E), the representative hematoxylin and eosin staining of colon sections (×40 and ×200, respectively; scale bars, 200 μm and 50 μm, respectively) (F), and the intracellular staining of cytokines in CD4+ T cells in the spleens and the colons (G). (H-K) Sex- and age-matched Foxp3 and Foxp1Foxp3 mice (n = 9–10) were induced with EAE by immunization with MOG peptide and pertussis toxin, and the severity of EAE was evaluated by the clinical score (H), loss of body weight (I), Luxol fast blue-hematoxylin and eosin staining of spinal cord sections (×40, scale bars, 200 μm) (J), and the intracellular staining of cytokines in T cells from the spinal cords and the brains (K). Data in (A-C, H-K) represent at least three independent experiments. Data in (D-G) represent two independent experiments. Data in (D, H, I) are mean ± SEM, *P < 0.05 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. CD, cluster of differentiation; DSS, dextran sulfate sodium; EAE, experimental autoimmune encephalomyelitis; Foxp1, forkhead box P1; IFNγ, interferon gamma; IL, interleukin; MOG, myelin oligodendrocyte glycoprotein; Treg, regulatory T; YFP, yellow fluorescent protein.

Treg cells exert an important function in suppressing various inflammatory responses [9]. Thus, we further addressed whether Foxp1 is indispensable for the restriction of inflammatory responses by Treg cells. In the DSS-induced colitis model, we found that with DSS treatment, WT mice survived well and had no weight loss; but in stark contrast, Foxp1Foxp3 mice began to lose weight at d5 post DSS treatment and started to die at d8, with the mortality rate up to 80% at d11 (). Consistent with the severe weight loss and low survival rate, we found that the colons of Foxp1Foxp3 mice treated with DSS exhibited more immune cell infiltration than those from WT mice (). We also found that, compared with those from WT mice, the conventional T cells in the spleens and the colons from Foxp1Foxp3 mice had higher percentages of inflammatory cytokines IFNγ- and IL-17A–producing cells (). To further confirm these results, we took advantage of the EAE disease model, in which MOG peptide immunization with pertussis toxin (PT) elicits strong TH1 and TH17 inflammatory responses in vivo [36]. We found that Foxp1Foxp3 mice had obviously higher clinical scores and more body weight loss than did WT mice during the development of EAE disease (). In line with more severe EAE symptoms in Foxp1Foxp3 mice, we observed more demyelination, indicated by the loss of Luxol fast blue staining in spinal cords (), more lymphocytes infiltration (), and more IFNγ- or IL-17A–producing T cells in the brains and the spinal cords of Foxp1Foxp3 mice (), suggesting that Foxp1-deficient Treg cells have impaired functions in suppressing inflammation in the EAE model. Thus, by using two different inflammatory models, we have demonstrated that Foxp1 is crucial for Treg suppressive function during anti-inflammatory responses.

Foxp1 regulates iTreg cell differentiation and maintenance

Studies have shown that Foxp3+ Treg cells can also be induced from CD4+CD25− naive T cells in the periphery or in vitro, which are named pTreg and iTreg cells, respectively [37]. As we have demonstrated that Foxp1 is important for the homeostasis and the suppressive function of tTreg cells (), we went further to ask whether Foxp1 plays any role in the differentiation and function of pTreg and iTreg cells. To circumvent the abnormal development of naive CD4+ T cells caused by Foxp1 deletion before T cell maturation, we used the inducible Foxp1 conditional knockout mouse line, the Foxp1Cre-ERT2Rosa26YFP mouse [31]. We treated the Foxp1Cre-ERT2+Rosa26YFP mice with tamoxifen, isolated naive CD4+ T cells, and cultured the CD4+ T cells under iTreg condition with different concentrations of IL-2. We found that the loss of Foxp1 led to dramatically reduced generation of CD25highFoxp3+ iTreg cells under all the IL-2 culture conditions (), suggesting impaired iTreg differentiation in the absence of Foxp1.

Foxp1 regulates iTreg cell differentiation and maintenance.

(A) Foxp1Rosa26 and Foxp1Cre-ERT2+Rosa26 mice were treated with tamoxifen for 6 days, and CD4+ T cells were isolated and cultured in the plate coated with 1 μg/mL anti-CD3/CD28 antibodies in the presence of 0.3 μM 4-hydroxytamoxifen, 5 ng/mL TGF-β, and various concentrations of recombinant human IL-2 for 2 days. Cells were detached and cultured for an additional day. Foxp3 and CD25 expression were analyzed by flow cytometry; numbers adjacent to the outlined areas represent the percentages of gated cells. (B-D) Sorted YFP+ Foxp1Foxp3 or Foxp3 iTreg cells differentiated in vitro were cotransferred with sorted CD44lowCD62LhighCD25− WT naive CD4+ T cells into Rag1− mice (n = 3). Mice were analyzed at 6 weeks post T-cell transfer. (B) Flow cytometry analysis of Foxp3+YFP+ Treg cells in the spleens and the colons of the recipient mice (left panel) and quantification of YFP+ Treg cell frequencies (right panel); the numbers adjacent to the outlined area indicate percentage of YFP+ Treg cells. (C) The severity of colitis was evaluated by loss of body weight of the recipient mice and (D) intracellular staining of cytokines in T cells from the colons; the numbers adjacent to the outlined area indicates percentage of gated cells. Data represent at least two independent experiments. Data in (B: right panel, C) are mean ± SEM, *P < 0.05 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. CD, cluster of differentiation; Foxp1, forkhead box P1; Foxp3, forkhead box P3; IFNγ, interferon gamma; IL, interleukin; iTreg, induced Treg cells; Rag1, recombination activating 1; TGF-β, transforming growth factor, beta; Treg, regulatory T; WT, wild-type; YFP, yellow fluorescent protein.

To elucidate the function and maintenance of iTreg cells in vivo, we cotransferred the iTreg cells differentiated in vitro with naive CD4+ T cells from the WT mice into the Rag1− recipient mice. Six weeks after cells transfer, we found that the frequency of Foxp1-deficient iTreg cells was significantly lower than that of WT iTreg cells in the spleens and the colons in the recipient mice (), suggesting that the maintenance of Foxp1-deficient iTreg cells in vivo is impaired. Consistently, we found that the recipient mice transferred with Foxp1-deficient iTreg cells developed more severe colitis, indicated by the greater body weight loss () and more IFNγ- or IL-17A–producing T cells in the colons (). Taken together, our data demonstrated that Foxp1 is important for both iTreg generation and maintenance.

Foxp1 helps maintain Treg suppressive function through controlling stable Foxp3 expression

To understand the mechanism by which Foxp1 regulates Treg suppressive function, we first examined the expression levels of Foxp3 in Treg cells. We found that the expression of Foxp3 was reduced in Foxp1-deficient rTreg cells at both the mRNA level () and the protein level (), consistent with an earlier observation that the Foxp1-deficient rTreg cells expressed lower levels of YFP, the indicator of Foxp3 expression in mice with a Foxp3 allele (). These results suggest that Foxp1 likely regulates Foxp3 expression in Treg cells at the transcriptional level. Unexpectedly, we found that there was no obvious difference in Foxp3 mRNA and protein levels between WT and Foxp1-deficient aTreg cells (), suggesting that somehow the defective Foxp3 expression in the absence of Foxp1 is compensated in the steady-state aTreg cells with a yet unknown mechanism.

To determine whether Foxp1 may directly regulate Foxp3 expression in rTreg cells, we analyzed the CNSs of the Foxp3 locus by bioinformatics. And the conserved forkhead-binding sites were found in the promoter, CNS2, and CNS0 regions of the Foxp3 locus, respectively (). Using the chromatin immunoprecipitation (ChIP) approach, we found that Foxp1 bound to these three regions of the Foxp3 locus (), suggesting that Foxp3 is a direct target of Foxp1.

It has been reported that CNS2 is required for the maintenance of heritable Foxp3 expression, and it protects Treg lineage identity in inflammatory cytokine environments [15,16]. Because Foxp1 bound to the promoter region and the CNS2 region of Foxp3, we reasoned that Foxp1 might be critical for stable Foxp3 expression, and thereby for Treg stability. Indeed, we found that during cell division upon T-cell receptor (TCR) stimulation, Foxp3 expression was greatly compromised in the Foxp1-deficient rTreg (). Surprisingly, although the Foxp3 levels were normal in Foxp1-deficient aTreg cells in the steady state (), upon TCR stimulation, the heritable Foxp3 expression was not sustained in Foxp1-deficient aTreg cells (). CNS2 has been shown to help maintain Foxp3 expression under a TH2 and TH17 cytokine environment [15,16]. We found that the addition of IL-4 to TCR-stimulated Treg cells greatly down-regulated Foxp3 in both Foxp1-deficient rTreg cells and aTreg cells, but not in WT Treg cells (. Similarly, we observed the unstable Foxp3 expression in Foxp1-deficient iTreg cells upon TCR stimulation (). Taken together, these results suggest that Foxp1, through its binding to the CNS2 and the promoter regions of the Foxp3 locus, plays a vital role in stabilizing the Foxp3 expression.

Fig 6

Foxp1 helps maintain Treg cell suppressive function through controlling stable Foxp3 expression.

(A) Relative Foxp3 mRNA levels in highly purified rTreg (n = 4) and aTreg cells (n = 6) from Foxp3 and Foxp1Foxp3 mice. (B) Immunoblot analysis of Foxp3 in rTreg and aTreg cells from Foxp3 and Foxp1Foxp3 mice. (C) Flow cytometry analysis of the expression of Foxp3 in rTreg and aTreg cells from the spleens of Foxp3 and Foxp1Foxp3 mice (left panel), and the corresponding mean fluorescence intensity (MFI) of Foxp3 (right panel) (n = 4). (D) Foxp1 ChIP was performed in Foxp3 and Foxp1Foxp3 iTreg cells to detect the Foxp1 binding to Foxp3 CNS0, CNS2, and promoter regions; a region in the Gmpr locus was used as a negative control, n = 3–4. (E, F) Sorted rTreg and aTreg cells from Foxp3 and Foxp1Foxp3 mice were labeled with CellTrace and cultured in the plate coated with 1 μg/mL anti-CD3 and -CD28 antibodies in the presence of 200 U/mL IL-2 for 3 days. (E) Foxp3 expression during Treg cell division in rTreg and aTreg cells. The numbers below the red line represent the percentages of Foxp3low Treg cells. (F) Foxp3 expression in rTreg and aTreg cells during Treg cell division; numbers adjacent to the histograms represent Foxp3+ Treg frequencies. (G) Foxp3 expression in Foxp3 and Foxp1Foxp3 rTreg and aTreg cells cultured with plate-bound anti-CD3/CD28 antibodies in the presence or absence of IL-4 for 3 days. (H-I) EAE was induced in Foxp1Foxp3 mice by immunization with MOG peptide and pertussis toxin (n = 3–8), followed by transfer of sorted Foxp3 and Foxp1Foxp3 Treg cells infected with control retroviruses (RV-ctrl) or retroviruses expressing Foxp3 (RV-Foxp3) at d4. The severity of EAE was evaluated by the clinical score (H) and intracellular staining of cytokines in the T cells from the spinal cords and the brains (I); the Foxp1Foxp3 mice without Treg cell transfer were controls. The numbers adjacent to the outlined area indicate the percentage of gated cells. Data in (A-I) represent at least three independent experiments. Data in (A, C: right panel, D, H) are mean ± SEM, *P < 0.05, **P < 0.01, ***P < 0.001 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; CD, cluster of differentiation; ChIP, chromatin immunoprecipitation; CNS0, conserved noncoding sequence 0; CNS2, conserved noncoding sequence 2; d, day; EAE, experimental autoimmune encephalomyelitis; Foxp1, forkhead box P1; Foxp3, forkhead box P3; Gmpr, guanosine monophosphate reductase; IL, interleukin; iTreg, induced Treg cells; max, maximum; MFI, mean fluorescence intensity; MOG, myelin oligodendrocyte glycoprotein; rTreg, resting Treg; RV-ctrl, control retrovirus; RV-Foxp3, retrovirus expressing Foxp3; TCR, T-cell receptor; Treg, regulatory T.

Studies have shown that unstable Foxp3 expression affects Treg cell functions [15,16]. We reasoned that the unstable Foxp3 expression in the absence of Foxp1 would contribute to the defective suppressive function of Foxp1-deficient Treg cells. To test this, we constructed a Foxp3 retroviral expression vector (RV-Foxp3) that restored the Foxp3 expression levels in Foxp1-deficient Treg cells to the levels of WT Treg cells (). In the EAE model, we found that Foxp1Foxp3 mice, which were transferred with Foxp1-deficient Treg cells overexpressing Foxp3, had significantly remitted symptoms from d15 post EAE induction (). Consistently, these mice had fewer IL-17A– and IFNγ-producing T cells in the spinal cords and the brains than Foxp1Foxp3 mice transferred with control Foxp1-deficient Treg cells (). All these results suggest that Foxp1-mediated regulation of Foxp3 expression is important for Treg suppressive function.

Foxp1-dependent transcriptional programs in Treg cells

To further gain insight into the mechanism by which Foxp1 regulates Treg cell homeostasis and function, we analyzed the global gene expression profiles of rTreg and aTreg cells from Foxp1Foxp3 and control Foxp3 (WT) mice by RNA sequencing (RNA-Seq). We identified 444 and 400 genes up-regulated or down-regulated, respectively, by more than 1.5-fold in Foxp1-deficient rTreg cells compared with WT rTreg cells, and 218 and 134 genes up-regulated or down-regulated, respectively, by more than 1.5-fold in Foxp1-deficient aTreg cells compared with WT aTreg cells ( and and ).

Foxp1-dependent transcriptional programs in Treg cells.

(A) Heatmap of the representative genes differentially expressed between highly purified Foxp3 and Foxp1Foxp3 rTreg cells (fold change ≥1.5). (B) Enriched GO categories of the genes differentially expressed between Foxp3 and Foxp1Foxp3 rTreg cells (Benjamin <0.05). (C) Enrichment of gene signatures related to cell cycle and oxidative phosphorylation in Foxp3 versus Foxp1Foxp3 rTreg cells through gene set enrichment analysis (GSEA). Representative enriched gene sets are shown (FDR q value <0.05). (D) Enrichment of gene signatures related to Treg signature gene cluster 6, TGF-β signaling pathway, and Foxo signaling pathway in Foxp3 versus Foxp1Foxp3 rTreg cells through GSEA. Representative enriched gene sets are shown (FDR q value <0.05). Two rounds of RNA sequencing were performed. Data associated with this figure can be found in the supplemental data file (. FDR, false discovery rate; Foxo, forkhead box O; Foxp1, forkhead box P1; Foxp3, forkhead box P3; GO, Gene Ontology; GSEA, gene set enrichment analysis; rTreg, resting Treg; TGF-β, transforming growth factor, beta; Treg, regulatory T.

In WT aTreg cells, there was a set of up-regulated genes compared with WT rTreg cells (2-fold as threshold). We defined these genes as the putative aTreg cell signature genes. Strikingly, we found that 137 up-regulated genes (31% of total up-regulated genes) in Foxp1-deficient rTreg cells were the aTreg cell signature genes, accounting for approximately 11% of the aTreg signature genes ( and ), suggesting that Foxp1 is a key transcription factor that normally restrains the differentiation of aTreg cells from rTreg cells.

The differentially expressed genes between WT and Foxp1-deficient rTreg cells were analyzed and found to be clustered into different groups that include cell cycle–regulating genes, transcription factors, genes involved in signaling transduction, and metabolism-related genes (), suggesting that Foxp1 regulates many cellular programs in Treg cells. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the up-regulated genes in Foxp1-deficient rTreg cells were greatly enriched in the cell cycle–related pathway (). The gene set enrichment analysis (GSEA) showed similar results, that the up-regulated genes in Foxp1-deficient rTreg cells were significantly enriched in cell cycle and oxidative phosphorylation pathways (), suggesting that Foxp1 enforces Treg cell quiescence through regulating cell cycle–and metabolism-related pathways in Treg cells, which is quite similar to our observation in Foxp1-deficient CD8+ T cells [38]. To gain more insights into the mechanism by which Foxp1 regulates Treg cell quiescence, we checked the expression of E2 promoter binding factor (E2f) family genes and phosphoinositide-3-kinase interacting protein 1 (Pik3ip1), which have been shown to regulate CD8+ T-cell quiescence [38]. We found that Foxp1 deletion in rTreg cells led to significantly increased expression levels of E2f1 and its target, cyclin-dependent kinase 1 (Cdk1), but decreased expression levels of Pik3ip1. However, we did not observe obvious changes in the expression levels of E2f2 and E2f3 (. These results suggest that Foxp1 enforces Treg cell quiescence by a similar mechanism involved in regulating CD8+ T-cell quiescence.

In contrast to the up-regulated genes associated with Treg cell activation, the down-regulated genes in Foxp1-deficient rTreg cells were greatly enriched in the Treg signature gene cluster 6 [39], the transforming growth factor, beta (TGF-β) signaling pathway, and the Foxo signaling pathway (), which are important for Treg identity and function. The Treg signature gene cluster 6 has been shown to contain the signature genes not affected by Foxp3, TGF-β signaling, or TCR activation [39]. Interestingly, our data showed that many genes in this cluster may be regulated by Foxp1 (). As TGF-β signaling is important for iTreg differentiation, one of the reasons for the impaired iTreg development of Foxp1-deficient T cells could be the down-regulation of those genes related to the TGF-β signaling pathway in the absence of Foxp1. Collectively, our results suggest that Foxp1 negatively regulates the expression of genes involved in Treg cell activation, and at the same time promotes the expression of a set of genes important for Treg cell identity and function.

Foxp1 and Foxp3 coordinate the regulation of CTLA-4 expression in Treg cells

It has been reported that Foxp subfamily proteins can form homodimers and heterodimers [33,40]. We confirmed the interaction between Foxp1 and Foxp3 in iTreg cells by immunoprecipitation (). We reasoned that Foxp1 and Foxp3 might co-regulate some important genes in Treg cells. Studies of Foxp3 ChIP have revealed many Foxp3 target genes in murine Treg cells [41]. By cross-examining the Foxp3 ChIP data set with our RNA-seq data set, we found that about 5% of up- or down-regulated genes in Foxp1-deficient rTreg cells (41 of 844 differential genes) were Foxp3-bound genes, accounting for about 8% of Foxp3 targets ( and ). These results suggest that Foxp3 and Foxp1 function largely independently in Treg cells.

Nevertheless, CTLA-4, which is critical for Treg suppressive function, was one of the down-regulated genes in Foxp1-deficient Treg cells and among the potential genes co-regulated by Foxp1 and Foxp3 ( and . We recently have shown that Foxp1 directly regulates the expression of CTLA-4 in activated conventional CD4+ T cells [42]. Whether CTLA-4 is also a direct target of Foxp1 in Treg cells is not clear. To examine this, we examined the mRNA and protein levels of CTLA-4 in Treg cells. We found that the expression levels of CTLA-4 were decreased in both rTreg and aTreg cells in the absence of Foxp1 (). Furthermore, Foxp1 ChIP assay in Treg cells showed that Foxp1 bound to the promoter of Ctla4 (), suggesting that CTLA-4 is a direct target of Foxp1 in Treg cells, as well.

Fig 8

Foxp1 and Foxp3 coordinate the regulation of CTLA-4 expression in Treg cells.

(A) Heatmap of the representative common target genes of Foxp1 and Foxp3 in rTreg cells. (B) The relative mRNA levels of CTLA-4 in purified rTreg and aTreg cells from Foxp3 and Foxp1Foxp3 mice, n = 4. (C) Flow cytometry analysis of CTLA-4 expression in splenic rTreg and aTreg cells from Foxp3 and Foxp1Foxp3 mice. (D) Predicted conserved forkhead-binding site in the Ctla4 promoter (left panel); Foxp1 ChIP analysis was performed in Foxp3 and Foxp1Foxp3 iTreg cells (right panel) (n = 3); −2.4 kb upstream of TSS at the Ctla4 locus was used as a negative control. (E) Foxp3 ChIP was performed in Foxp3 and Foxp1Foxp3 iTreg cells to detect Foxp3 binding to the Ctla4 promoter region in the presence or absence of Foxp1, n = 2. A region in the Gmpr locus was used as a negative control. (F) The relative mRNA levels of CTLA-4 in Foxp1Foxp3 iTreg cells infected with control retroviruses (RV-ctrl), retroviruses expressing Foxp1A (RV-Foxp1A), retroviruses expressing Foxp3 (RV-Foxp3), or retroviruses expressing both RV-Foxp1A and RV-Foxp3 (RV-Foxp1A+Foxp3), n = 3. WT iTreg cells infected with control retroviruses were used as a control. Data in (A, D-F) are representative of at least two independent experiments. Data in (B, C) are representative of at least three independent experiments. Data in (B, D, E, F) are mean ± SEM, *P < 0.05, **P < 0.01, ***P < 0.001 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; ChIP, chromatin immunoprecipitation; CTLA-4, cytotoxic T-lymphocyte-associated protein 4; Foxp1, forkhead box P1; Foxp3, forkhead box P3; Gmpr, guanosine monophosphate reductase; iTreg, induced Treg cells; rTreg, resting Treg; Treg, regulatory T; TSS, transcription start site; WT, wild-type.

It has been reported that Foxp3 can bind to the Ctla4 promoter. Thus, it is intriguing to ask whether Foxp1 works together with Foxp3 to regulate the expression levels of CTLA-4 in Treg cells. We performed Foxp3 ChIP assay in Foxp1-sufficient and -deficient iTreg cells, and found that the loss of Foxp1 significantly reduced the binding of Foxp3 to the Ctla4 promoter region (). Of note, the protein levels of Foxp3 had only a very minor decrease in Foxp1-deficient iTreg cells compared with those of WT iTreg cells (), suggesting that the decreased binding of Foxp3 to Ctla4 promoter was mainly attributed to the absence of Foxp1. Consistently, although retroviral overexpression of Foxp1 or Foxp3 alone partially rescued the CTLA-4 expression in Foxp1-deficient iTreg cells, overexpression of Foxp1 together with Foxp3 completely rescued the expression of CTLA-4 in Foxp1-deficient iTreg cells. (). Thus, our results demonstrate that Foxp1 and Foxp3 coordinate to regulate the expression levels of CTLA-4.

Discussion

Foxp1 is an evolutionally conserved transcription factor that is required for the proper development of cardiovascular, intestinal, neural, pulmonary, and lymphoid tissues [29,30,43-47]. In this study, we provided evidence demonstrating that Foxp1 is also an important regulator of Treg cell homeostasis and suppressive function. We found that Foxp1 helps maintain the rTreg cell pool and restrain the differentiation of rTreg cells into aTreg cells. Interestingly, we also found that the ablation of Foxp1 in Treg cells dampens the Treg suppressive function in both homeostatic and pathological situations. Mice with Foxp1 deletion in Treg cells spontaneously developed inflammatory disease, and DSS-induced colitis and EAE were more severe in Foxp1Foxp3 mice. Mechanistically, we found that Foxp1 stabilizes Foxp3 expression to maintain both the Treg suppressive function in anti-inflammatory responses and the Treg cell stability under TCR stimulation or destabilizing conditions. In addition, we found that CTLA-4 is a direct target of Foxp1 in Treg cells, and the absence of Foxp1 attenuates the Foxp3 binding to the Ctla4 promoter. Our study identifies Foxp1 as a critical regulator of Treg cell quiescence, suppressive function, and stability in homeostatic and inflammatory responses.

Akin to conventional CD4+ T cells, tTreg cells can be divided into resting and activated subpopulations with relatively different functions. While rTreg cells that reside in lymphoid tissues play crucial roles in preventing lymphoproliferative disease, aTreg cells are the dominant Treg cell population in nonlymphoid tissues as potent suppressors in anti-inflammatory responses [4,5]. The balance between rTreg and aTreg is important for immune homeostasis. In our study, we found that Foxp1 deficiency leads to a gradual loss of rTreg cells. Given that rTreg cells function as suppressors of lymphoproliferation, the reduced rTreg cell number may explain the early splenomegaly and lymphadenopathy in young Foxp1-deficient mice. Unexpectedly, despite the increased numbers of aTreg cells, Foxp1Foxp3 mice gradually develop spontaneous inflammatory diseases with age, and have more severe disease in DSS-induced colitis and EAE models, indicating that, although Foxp1 helps restrain the differentiation of aTreg cells from rTreg cells, it is also important for the suppressive function of aTreg cells.

In addition to tTreg cells, studies have shown that Foxp3+ Treg cells can also develop in the periphery from mature conventional CD4+ T cells (iTreg cells) [37]. It is an intriguing question whether Foxp1 is important for both tTreg cells and iTreg cells. Recently, Ghosh and colleagues showed that Foxp1 is vital for the iTreg cell differentiation and functions but is dispensable for the function of tTreg cells [48]. However, by using the same Foxp1 mouse line, Konopacki and colleagues showed that Foxp1 has an essential function in regulating tTreg cells [49]. The fact that Ghosh and colleagues did not observe functional defects in Foxp1-deficient tTreg cells could be due to some Foxp1+ Treg cells present in their Foxp1Foxp3 mice [48]. In our study, we have used a different Foxp1 mouse line: Foxp1 is efficiently deleted in Treg cells in our Foxp1Foxp3 mice () [30,31], and we demonstrated that Foxp1 is important for the differentiation and suppressive function of tTreg cells. Thus, we conclude that Foxp1 plays vital roles in both tTreg and iTreg cells.

In conventional naive T cells, Foxp1 actively maintains T-cell quiescence largely by restraining IL-7Rα expression [31]. However, in this study we did not find that Foxp1 regulates the expression of IL-7Rα in rTreg cells, suggesting that the mechanisms underlying the transcriptional regulation of IL-7Rα are different between conventional CD4+ T and Treg cells. Our earlier study has also shown that Foxp1 does not control IL-7Rα expression in pre-B cells that express high levels of Foxp1 and IL-7Rα [29]. Thus, for the same target gene but in different types of cells, Foxp1 seems to engage different transcriptional partners or complexes to exert its regulation. Although Foxp1 does not control IL-7Rα expression in Treg cells, we found that Foxp1 regulates the expression levels of E2f1 and Pik3ip1 in Treg cells, two genes involved in the regulation of CD8+ T cell quiescence [38], suggesting that Foxp1 enforces Treg quiescence in a similar way as in CD8+ T cells: by simultaneously repressing key pathways in cellular metabolism and cell cycle progression.

The work of Ghosh and colleagues shows that Foxp1 binds to the promoter and the CNS2 region of the Foxp3 locus, which retains permissive histone modifications [48]. In our study, we have confirmed that Foxp1 binds to the promoter and the CNS2 region of the Foxp3 locus and, additionally, the CNS0 region. We find that Foxp1 sustains de novo expression of Foxp3 in rTreg cells and maintains stable Foxp3 expression in both rTreg and aTreg cells upon TCR stimulation. Previous studies have reported that Foxp3 CNS2 functions as a sensor of TCR signal to regulate and maintain Foxp3 expression [15,16]. Thus, it is possible that in quiescent rTreg cells, the control of Foxp3 expression by Foxp1 is mainly via promoter or CNS0, but may not involve CNS2; however, upon TCR stimulation, Foxp1 may engage in the formation of the loop between the promoter and the CNS2 of Foxp3 locus, which is required for the maintenance of high levels of Foxp3 expression. Several proteins such as Gata3, Ets-1, CREB, Stat5, Foxo1, Foxp3, Runx1, and CBFβ have been shown to bind to the promoter or the CNS2 region of Foxp3 gene [17-21,23,24,27,50]. This raises interesting questions about whether Foxp1 interacts with those proteins, especially, with Foxp3 in the regulation of Foxp3 expression and whether Foxp1 plays a nonredundant role in the regulation of Foxp3 expression. Others’ and our data show that Foxp1 and Foxp3 can form heterodimers. Thus, it is interesting to propose that the formation of the regulatory loop between the promoter and the CNS2 of the Foxp3 locus requires Foxp1 and Foxp3 heterodimers, which warrants further investigation.

A variety of molecules are involved in the suppressive function of Treg cells, including CTLA-4, a key “checkpoint” in immune tolerance. CTLA-4 deficiency in Treg cells impairs their suppressive function both in vivo and in vitro [51]. CTLA-4 is also known as a Foxp3 direct target. In our study, we showed that Foxp1 and Foxp3 coordinate the regulation of CTLA-4 expression levels. In the absence of Foxp1, Foxp3 binds less to the promoter of Ctla4, suggesting that Foxp3-mediated CTLA-4 regulation is dependent on the presence of Foxp1. As Foxp1 can form heterodimers with Foxp3, it is also interesting to ask how many Foxp3 targets are regulated by Foxp1/Foxp3 dimers. In our RNA-Seq data, although the potential common targets of Foxp1 and Foxp3 are not too many, these genes, such as CTLA-4 and Foxp3, seem to be important for Treg function and identity. Thus, how Foxp1 homodimers or heterodimers with other proteins contribute to Treg cell differentiation and Treg suppressive function is worth further investigation.

In summary, we have defined the critical roles of another Foxp subfamily transcription factor, Foxp1, in regulating Treg cell homeostasis and function. In addition, the interplay between Foxp1 and Foxp3 in regulating key Treg genes also provides new understanding of the transcriptional programs in Treg cells.

Materials and methods

Ethics statement

All procedures were carried out in accordance with guidelines of the Animal Facility of IPS. The experimental procedures were approved by the Institutional Animal Care and Use Committee of Institut Pasteur of Shanghai, China (A2019014). All efforts were made to minimize suffering.

Mice

Foxp1, Cre-ERT2, and Rosa transgenic mice have been previously described [30,31]. Foxp1 were bred with Foxp3-YFP-Cre knock-in mice from Jackson Laboratories to generate Foxp1Foxp3 mice. Control animals were age-/sex-matched control Foxp3 mice (WT) or Foxp1 mice. Unless otherwise stated, all experiments were performed with mice at 6–8 weeks of age. All mice were maintained under specific pathogen-free conditions at Institute Pasteur of Shanghai, China.

Flow cytometry, cell sorting, and intracellular staining

These procedures were done as described previously [31]. The sorted cell populations were >98% pure. Antibodies used in FACS analysis and cell sorting were as follows: APC-Cy7-anti-CD4 (GK1.5, Thermo Fisher Scientific, Waltham, MA), PE-Cy7-anti-CD4 (GK1.5, Thermo Fisher Scientific), PE-Cy7-anti-CD8 (53–6.7, Thermo Fisher Scientific), PerCP-Cy5.5-anti-CD8 (53–6.7, BioLegend, San Diego, CA), PE-anti-CD44 (IM7, Thermo Fisher Scientific), APC-anti-CD62L (MEL-14, Thermo Fisher Scientific), PE-Cy7-anti-CD25 (PC61.5, Thermo Fisher Scientific), Pacific Blue-anti-Foxp3 (FJK-16s, Thermo Fisher Scientific), APC-anti-Foxp3 (FJK-16s, Thermo Fisher Scientific), PerCP-Cy5.5-anti-interferon-γ (XMG1.2, Thermo Fisher Scientific), Alexa Fluor 700-anti-IL-17A (TC11-18H10, BD Pharmingen), Pacific Blue-anti-Granzyme B (NGZB, Thermo Fisher Scientific), PE-anti-CTLA-4 (UC10-4B9, BioLegend), BV605-anti-CD44 (IM7, BioLegend), and Alexa Fluor 647-goat-anti-rabbit IgG (A-21244, Thermo Fisher Scientific). Dead cells were excluded through live/dead staining by using Live/Dead Fixable Aqua Dead Cell staining kit (Thermo Fisher Scientific).

For intracellular staining of Foxp1, the staining procedures were performed as described previously [32]. For intracellular staining of Foxp3 and CTLA-4, cells were fixed with 3.7% formaldehyde following cell surface staining (antibodies mentioned above) and were permeabilized with 0.2% Triton X-100. Then, cells were stained with Pacific Blue-anti-Foxp3 or PE-anti-CTLA-4 at 4°C for 90 minutes. For APC-anti-Foxp3 staining, we used the Foxp3 staining kit (Thermo Fisher Scientific) and performed it as per the manufacturer’s instructions.

For intracellular staining of cytokines, cells were stimulated for 6 hours with 10 ng/mL PMA (phorbol 12-myristate 13-acetate, Sigma, St. Louis, MO) and 1 μg/mL ionomycin (Sigma). Before staining, cells were treated for 2 hours with 10 μg/mL Brefeldin A (BFA; Sigma). The intracellular staining was done as described [31].

iTreg cell differentiation and transfer

Foxp1Rosa26 and Foxp1 Cre-ERT2Rosa26 mice were treated daily with tamoxifen (Sigma) (2 mg per mouse) for 5 days and rested for 1 day. CD4+ T cells were isolated and cultured with 1 μg/mL plate-bound anti-CD3 (BioXcell) and anti-CD28 (BioXcell), 5 ng/mL TGF-β, and indicated recombinant human IL-2 in the presence of 0.3 μM 4-hydroxytamoxifen (Sigma) for 2 days, followed by detaching and culturing in T-cell medium (Dulbecco’s Modified Eagle Medium [DMEM] containing 10% heat-inactivated FCS, 2 mM L-glutamine, penicillin-streptomycin, nonessential amino acids, sodium pyruvate, vitamins, 10 mM HEPES and 50 μM 2-mercaptoethanol) in the presence of 5 ng/mL TGF-β and the indicated recombinant human IL-2.

For iTreg cell transfer, naive (YFP–CD4+CD44lowCD62Lhigh) T cells from the spleens and the lymph nodes of Foxp3 and Foxp1Foxp3 mice were sorted by FACS Aria II device (BD Bioscience, San Jose, CA). Naive T cells were stimulated for 48 hours with plate-bound anti-CD3 (5 μg/mL, 145-2C11, BioXcell) and anti-CD28 (5 μg/mL, 37.51, BioXcell) in complete T-cell medium supplemented with TGF-β (5 ng/mL, R&D Systems, Minneapolis, MN) and IL-2 (200 U/mL, Beijing Four Rings Biopharmaceutical Co., Beijing, China). Then, the cells were detached and cultured in T-cell medium containing 5 ng/mL TGF-β and 200 U/mL recombinant human IL-2 for an additional 3 days. Afterwards, 4 × 105 sorted YFP+ iTreg cells were cotransferred with 1 × 106 purified WT naive CD4+ T cells (CD4+CD44lowCD62LhighCD25–) into Rag1− mice to perform the colitis assay.

Quantitative RT-PCR

rTreg and aTreg cells from Foxp3 or Foxp1Foxp3 mice were sorted out by a FACS Aria II device (BD Bioscience). Then, total RNA was isolated from whole cells using the TRIZOL reagent (Life technologies, Carlsbad, CA) according to the manufacturer’s instructions. Complementary DNA was reverse transcribed using the PrimeScript RT reagent kit (Takara, Tokyo, Japan). Quantitative PCR was performed with SYBR Green mix (Takara) on ABI Prism 7500 sequence detection system (Applied Biosystems). Relative mRNA expression of genes was normalized to GAPDH. The primers are listed as follows: Foxp1, forward: CTGAATCTGGTATCAAGTGTCACCCTCT, reverse: GATTCGAGAATGGCCTGCCTGA; Foxp3, forward: CAGAGTTCTTCCACAACATGGAC, reverse: CAGGGATTGGAGCACTTGTTG; Ctla4, forward: AGAACCATGCCCGGATTCTG, reverse: GCTCTGTTGGGGGCATTTTC; GAPDH, forward: TTCACCACCATGGAGAAGGC, reverse: GGCATGGACTGTGGTCATGA; E2f1, forward: GCATCCAGCTCATTGCCAAGAAGT, reverse: GAAAGCAGTTGCAGCTGTGTGGTA; E2f2, forward: TCAGAGTTGCTCCCTGAGCTTCAA, reverse: TTGAAGTTGCCTACGGCACGGATA; E2f3, forward: AGGGCCCATTGAGGTTTACTT, reverse: GTAGAAACCGAGCAGTCACTA; Cdk1, forward: CATGGACCTCAAGAAGTACCTGG, reverse: CAAGTCTCTGTGAAGAACTCGCC; Mcm5, forward: GGAGGCATTGAGACTGTTCCAG, reverse: AGACACCTGAGAGCCAATGGCA; and Pik3ip1, forward: ATCCTCGCTATTGGAGCTGGCATT, reverse: TCTGGTTGCTGTGCACAATGATGG.

Immunoblot and immunoprecipitation

Cells were washed with ice-cold PBS twice and lysed in RIPA buffer (20 mM Tris-HCl, pH 7.5; 150 mM NaCl; 1% NP-40; 0.1% SDS; 0.5% sodium deoxycholate) containing 1× Protease Inhibitor Cocktail (Roche, Mannheim, Germany). Cell extracts were separated by SDS-PAGE and transferred to a polyvinylidene difluoride membrane (Bio-Rad, Hercules, CA). The membrane was analyzed by immunoblot with rabbit anti-Foxp1 [31], mouse anti-Foxp3 (eBio7979, Thermo Fisher Scientific), or mouse anti-β-actin (AC74, Sigma). For co-immunoprecipitation, cell lysates were immunoprecipitated with appropriate antibodies (1 μg/mL) using protein A sepharose beads (GE Healthcare, Uppsala, Sweden) at 4 °C. Samples were used for immunoblot analysis with the indicated antibodies.

Treg cell proliferation

YFP+ CD44lowCD62LhighCD4+ T cells (rTregs) and YFP+ CD44highCD62LlowCD4+ T cells (aTregs) sorted from Foxp3 or Foxp1Foxp3 mice were labeled with CellTrace following the manufacturer’s instruction (CellTrace Violet Cell Proliferation Kit, Thermo Fisher Scientific). To assess the proliferation of purified Treg cells to TCR stimulation, labeled Treg cells were stimulated with the indicated concentration of anti-CD3 and anti-CD28 in the presence of 200 U/mL IL-2 for 3 days.

Retroviral transduction

The open reading frames of Foxp3 and Foxp1A were subcloned into the retroviral vector MSCV-IRES-VEX and MSCV-IRES-Tomato. Retroviruses were packaged in the HEK293T cell line as described previously [32]. Briefly, retroviral expression plasmids were transfected into HEK293T cells using Lipofectamine 2000 (Thermo Fisher Scientific). For transduction of retrovirus, purified Treg cells were stimulated with 10 μg/mL anti-CD3 and 10 μg/mL anti-CD28 in the presence of 500 U/mL IL-2 for 24 hours. Activated cells were spin transduced with 1 mL virus-containing medium with polybrene (6 μg/mL) at 650g, 30°C for 2 hours; transduced cells were cultured with complete T-cell medium supplemented with recombinant human IL-2 (500 U/mL). The cells were harvested on d4 for staining or adoptive transfer.

Colitis induction and lamina propria mononuclear cell isolation

For the induction of acute colitis, age-matched female Foxp3 or Foxp1Foxp3 mice were treated with 2.5% DSS (molecular weight 36,000–50,000 kDa; MP Biomedicals, Aurora, OH) dissolved in drinking water for 14 days. Mice were weighed every day and killed at d14. Cytokine production of CD4+ T cells from spleens and colons were analyzed by flow cytometry. The colons were also processed for tissue sectioning and HE staining.

For isolation of lamina propria (LP) mononuclear cells, dissected colons were washed with PBS followed by removing fat tissues and Peyer’s patches. Colon tissues were pretreated in 5% FBS DMEM containing 1 mM DTT (Sigma) and 30 mM EDTA (AMResco) at 37°C for 30 minutes in a shaker. Then, pretreated tissues were digested in 5% FBS DMEM containing Liberase (25 μg/mL, Roche) and DNase I (20 μg/mL, Sigma) at 37°C for 30 minutes in a shaker. Digested tissues were meshed in a 70-μm cell strainer, followed by percoll gradient (80%/40%) centrifugation. Mononuclear cells present in the interphase were collected, washed with PBS, and resuspended in T-cell medium.

EAE induction, Treg administration, and CNS mononuclear cell isolation

For EAE induction, 2-month-old Foxp3 or Foxp1Foxp3 mice were immunized subcutaneously with 200 μg of MOG35–55 peptide (GL Biochem, Shanghai, China) emulsified in Complete Freund’s Adjuvant (CFA; Sigma). Additionally, mice were injected intravenously with 200 ng of PT (List biological laboratories, Campbell, CA) on d0 and d2 following immunization. For Treg cell transfer, experimental mice were intravenously injected with 2 × 106 Foxp3-transduced Foxp1Foxp3 Treg cells in 200 μL PBS at d4 post EAE induction. Control groups were either injected with the same number of Foxp3 or Foxp1Foxp3 Treg cells infected with control retroviruses (RV-ctrl) or retroviruses expressing Foxp3 (RV-Foxp3) at d4 post EAE induction. The clinical scores and body weight of mice was evaluated daily after 6 days following immunization. Severity of EAE was assessed daily with a 0–5 scoring system [52]: 0, no clinical signs; 0.5, partially limp tail; 1, paralyzed tail; 2, loss in coordinated movement, hind limb paresis; 2.5, one hind limb paralyzed; 3, bilateral hind limb paralysis; 3.5 hind limbs paralyzed, weakness in forelimbs; 4, forelimb paralysis; 5, moribund.

For isolation of CNS mononuclear cells, mice were perfused with PBS to remove blood cells within tissues. CNS tissues (brains and spinal cords) were cut into small pieces and digested with Liberase (25 μg/mL, Roche) and DNaseI (20 μg/mL, Sigma) at 37°C for 30 minutes. Digested tissues were meshed in a 70-μm cell strainer, followed by percoll gradient (80%/40%) centrifugation. Mononuclear cells present in the interphase were collected, washed by PBS, and resuspended in culture medium.

Histopathology

For histological analysis, mouse tissues were fixed in 3.7% (vol/vol) formaldehyde, and sections were stained with hematoxylin and eosin, or Luxol fast blue-hematoxylin and eosin according to standard procedures.

ChIP

Naive CD4+ T cells from Foxp3 or Foxp1Foxp3 mice were cultured under iTreg-polarizing conditions (5 ng/mL TGF-β, 200 U/mL IL-2). ChIP of Foxp1, Foxp3 was performed as described [32]. Immunoprecipitated DNA and input DNA were analyzed by quantitative PCR (Applied Biosystems). The primers were as follows: Foxp3 CNS2, forward: ATCTGGCCAAGTTCAGGTTG, reverse: GGCGTTCCTGTTTGACTGTT; Foxp3 promoter, forward: CCCTGCAATTATCAGCACAC, reverse: TGTGGGAAACTGCCACATTA; Ctla4 promoter region, forward: GACTCCACGTCTCCAGGTCCTCAG, reverse: GGAAGCCGTGGGTTTAGCTGTTAC; Gmpr, forward: CAGCTGGAACAGCCTTGGAA, reverse: AAATGTCAAGGCCCCTGTGA. Control for CTLA4 promoter region in Foxp1 ChIP assay: forward: GTCAGAGTATTTTATCACAGCCAC, reverse: AGTCTACTGCAAAACCCCAAG.

RNA-Seq analysis

RNA sequencing

Total RNA was extracted using the TRIZOL reagent (Life technologies) according to the manufacturer’s instructions. RNA quality was assessed and sequenced using illumina Hiseq X platform in Novogene (Beijing, China).

RNAseq data analysis

Raw data from the sequencer first underwent quality control using FastQC. Next, reads were treated using trim_galore, in which illumina adapter sequence and low-quality reads(phred score <20) at the 3′ end were trimmed, and the paired reads were removed if any of the two reads did not meet the minimum length (20 bp). After that, duplicates were removed using Picard. For each sample, we counted the reads hits on each transcript using htseq-count.

Differential analysis

Differential analysis between Foxp3 or Foxp1Foxp3 Treg cell groups was done by a count-based method, limma, which is implemented in R, and voom is involved for normalization [53,54]. Paired samples analysis was employed to compensate batch effect. Significantly expressed genes were first screened by BH-adjusted P value 0.05 and further filtered with 1.5 fold-change and 2 fold-change for Foxp1Foxp3 Treg cells versus Foxp3 Treg cells, and WT aTreg cells versus WT rTreg cells, respectively.

Functional analysis

DAVID

Significant genes identified through comparisons between Foxp1Foxp3/Foxp3 or aTreg/rTreg were further used for enrichment analysis using DAVID. Biological processes (BPs) in GO and pathways in KEGG were chosen as significantly enriched terms with a Benjamini value less than 0.05.

GSEA

In parallel, we used GSEA v3.0 (Broad Institute, PreRanked mode) for enrichment analysis. To be consistent with the way for identifying significant genes, we used the t statistic output from limma as the metrics for ranking. A total of 1,000 gene set permutations were set as default, and gene sets were obtained through collecting pathways from KEGG and Biocarta, as well as BPs from GO. A gene set with an FDR q value less than 0.05 was considered significantly enriched.

Heatmap

For the heatmap, the log2 scale of the fragments per kilobase of transcripts per million fragments mapped (FPKM) values of the significant genes were used as input for heatmap generation.

Statistics

Statistical analyses were performed with Prism 6.0 (GraphPad). Two-tailed Student t tests and paired t tests were used for the calculation of P values, except for those in RNA-Seq analysis. *P < 0.05, **P < 0.01, and ***P < 0.001.

The gating strategy for flow cytometry analysis and cell sorting.

The gating strategy of CD44lowCD62Lhigh rTreg and CD44highCD62Llow aTreg (upper panel), and the purity check after cell sorting of CD44lowCD62Lhigh rTreg and CD44highCD62Llow aTreg in Foxp3 and Foxp1Foxp3 mice (lower panel). aTreg, activated Treg; rTreg, resting Treg.

(TIF)

Click here for additional data file.

The expression of Foxp1 in Treg cells from different tissues.

Intracellular staining of Foxp1 in Treg cells in spleens, lungs, livers, small intestines (SI) from Foxp3 mice (WT) (upper panel), and the corresponding MFI of Foxp1 (lower panel). Foxp1-deficient Treg cells were used as staining controls. Data are representative of two independent experiments. Data associated with this figure can be found in the supplemental data file (. Foxp1, forkhead box P1; MFI, mean fluorescence intensity; SI, small intestine; Treg, regulatory T; WT, wild-type.

(TIF)

Click here for additional data file.

Cell proliferation and apoptosis in Foxp3 and Foxp1 Treg cells.

(A) Flow cytometry analysis of Ki-67 expression in rTreg and aTreg cells in spleens from 3-week-old Foxp3 and Foxp1Foxp3 mice; numbers above bracketed lines indicate percent Ki-67+ Treg cells. (B) Proliferation of aTreg cells cultured with different concentrations of anti-CD3 and anti-CD28 antibodies in the presence of 200 U/mL IL-2 for 3 days, as shown by dilution of CellTrace; numbers above bracketed lines indicate percent proliferating cells. (C) Annexin V and 7-AAD staining in rTreg and aTreg cells in spleens from 6-week-old Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined area represent the percentage of gated cells. Data are representative of two independent experiments. aTreg, activated Treg; IL, interleukin; Ki-67, antigen identified by monoclonal antibody Ki 67; rTreg, resting Treg.

(TIF)

Click here for additional data file.

Lymphadenopathy and activated T-cell phenotype in Foxp1 mice.

(A) Representative picture of spleens and lymph nodes from 8-week-old sex- and age-matched Foxp3 and Foxp1Foxp3 mice. (B) Flow cytometry analysis of CD44 and CD62L expression in conventional CD4+ T cells (YFP−CD4+) and CD8+ T cells in spleens and lymph nodes from 8-week-old Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined area represent the percentage of gated cells. (C) Intracellular staining of cytokines in conventional CD4+ and CD8+ T cells in spleens, livers, and small intestines (SI) from 8-week-old Foxp3 and Foxp1Foxp3 mice; numbers adjacent to the outlined area represent the percentage of gated cells. (D) Intracellular staining of IL-4 in YFP−CD4+ T cells in the spleens from 28-week-old Foxp3 and Foxp1Foxp3 mice (left panel) and the corresponding statistics (right panel) (n = 3). Numbers adjacent to the outlined area indicate the percentage of IL-4+ cells in CD4+ T cells. (E) Hematoxylin and eosin staining of spinal cord sections (×100, scale bars, 100 μm) from Foxp3 and Foxp1Foxp3 mice induced with EAE by immunization with MOG peptide and PT. (F) Quantification of CD4+ and CD8+ T-cell frequencies in alive cells from brains and spinal cords of Foxp3 or Foxp1Foxp3 mice (n = 12–17). Data in (A-C, F) are representative of at least three independent experiments. Data in (D and E) are representative of two independent experiments. Data in (D: right panel, F) are mean ± SEM, **P < 0.01, ***P < 0.001 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. CD, cluster of differentiation; EAE, experimental autoimmune encephalomyelitis; IL, interleukin; MOG, myelin oligodendrocyte glycoprotein; PT, pertussis toxin; SI, small intestine.

(TIF)

Click here for additional data file.

Foxp1 sustains stable expression of Foxp3.

(A) Flow cytometry analysis of YFP expression in rTreg and aTreg cells from the spleens of Foxp3 and Foxp1Foxp3 mice. (B) Predicted conserved forkhead-binding sites (highlighted in red) in the CNS0, promoter, and CNS2 region of Foxp3, respectively. (C) Sorted YFP+ Foxp3 and Foxp1Foxp3 iTreg cells were stimulated with anti-CD3 and -CD28 beads for 3 days, and intracellular staining of Foxp3 in Foxp3 and Foxp1Foxp3 iTreg cells were analyzed by flow cytometry. (D) Intracellular staining of Foxp3 in Foxp3 and Foxp1Foxp3 iTreg cells infected with control retroviruses (RV-ctrl) or retroviruses expressing Foxp3 (RV-Foxp3) (left panel), and the corresponding MFI of Foxp3 (right panel). Data in (A, C) are representative of at least three independent experiments. Data in (D) are representative of two independent experiments. Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; CNS, conserved noncoding sequence; Foxp1, forkhead box P1; Foxp3, forkhead box P3; iTreg, induced Treg cells; MFI, mean fluorescence intensity; rTreg, resting Treg; YFP, yellow fluorescent protein.

(TIF)

Click here for additional data file.

Transcriptional programs regulated by Foxp1.

(A) Representative genes differentially expressed between highly purified Foxp3 and Foxp1Foxp3 aTreg cells. (B) Relationship between aTreg signature genes and Foxp1 repressed genes in rTreg cells. Two rounds of RNA sequencing were performed. (C) Relative mRNA levels of E2f1, E2f2, E2f3, Cdk1, Mcm5, and Pik3ip1 in purified rTreg cells and aTreg cells from 6–8-week-old Foxp3 and Foxp1Foxp3 mice, n = 3. Data in (C) represent at least three independent experiments. Data in (C) are mean ± SEM, *P < 0.05, **P < 0.01, ***P < 0.001 (two-tailed Student t test). Data associated with this figure can be found in the supplemental data file (. aTreg, activated Treg; Cdk1, cyclin-dependent kinase 1; E2f1, E2F transcription factor 1; E2f2, E2F transcription factor 2; E2f3, E2F transcription factor 3; Foxp1, forkhead box P1; Mcm5, minichromosome maintenance complex component 5; Pik3ip1, phosphoinositide-3-kinase interacting protein 1; rTreg, resting Treg.

(TIF)

Click here for additional data file.

Foxp1 interacts with Foxp3 in iTreg cells and Foxp1 coordinates with Foxp3 in the regulation of CTLA-4 expression in Treg cells.

(A) Co-immunoprecipitation of Foxp1 and Foxp3 was analyzed in iTreg cells, and isotype matched antibody (IgG) and Foxp1Foxp3 cells were used as controls. (B) Flow cytometry analysis of co-regulation of CTLA-4 by Foxp1 and Foxp3. The MFI of CTLA-4 in Foxp1Foxp3 iTreg cells infected with control retroviruses (RV-ctrl), retroviruses expressing Foxp1A (RV-Foxp1A), retroviruses expressing Foxp3 (RV-Foxp3), or retroviruses expressing both RV-Foxp1A and RV-Foxp3 (RV-Foxp1A+Foxp3) are shown. WT iTreg cells infected with control retroviruses were used as a control. Data in (A and B) are representative of two independent experiments. Data associated with this figure can be found in the supplemental data file (. CTLA-4, cytotoxic T-lymphocyte-associated protein 4; Foxp1, forkhead box P1; Foxp3, forkhead box P3; iTreg, induced Treg cells; MFI, mean fluorescence intensity; Treg, regulatory T; WT, wild-type.

(TIF)

Click here for additional data file.

Related to Fig 7.

Differentially expressed genes between Foxp3 and Foxp1Foxp3 rTreg cells. rTreg, resting Treg.

(XLSX)

Click here for additional data file.

Differentially expressed genes between Foxp3 and Foxp1Foxp3 aTreg cells. aTreg, activated Treg.

(XLSX)

Click here for additional data file.

aTreg signature genes. aTreg, activated Treg.

(XLSX)

Click here for additional data file.

Related to Fig 8.

Common target genes of Foxp1 and Foxp3 in rTreg cells. Foxp1, forkhead box P1; Foxp3, forkhead box P3; rTreg, resting Treg.

(XLSX)

Click here for additional data file.

Numerical values of presented diagrams.

(XLS)

Click here for additional data file.