CXXC-finger protein 1 associates with FOXP3 to stabilize homeostasis and suppressive functions of regulatory T cells

Xiaoyu Meng; Yezhang Zhu; Kuai Liu; Yuxi Wang; Xiaoqian Liu; Chenxin Liu; Yan Zeng; Shuai Wang; Xianzhi Gao; Xin Shen; Jing Chen; Sijue Tao; Qianying Xu; Linjia Dong; Li Shen; Lie Wang

doi:10.7554/eLife.103417.1

eLife Assessment

This study presents important findings on the role of CXXC-finger protein 1 in regulatory T cell gene regulation and function. The evidence supporting the authors' claims is solid, with mostly state-of-the-art technology, although the inclusion of more mechanistic insights would have strengthened the study. The work will be of relevance to immunologists interested in regulatory T cell biology and autoimmunity.

https://doi.org/10.7554/eLife.103417.1.sa3

Significance of findings

important: Findings that have theoretical or practical implications beyond a single subfield

landmark
fundamental
important
valuable
useful

Strength of evidence

solid: Methods, data and analyses broadly support the claims with only minor weaknesses

exceptional
compelling
convincing
solid
incomplete
inadequate

During the peer-review process the editor and reviewers write an eLife assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife assessments

Abstract

FOXP3-expressing regulatory T (T_reg) cells play a pivotal role in maintaining immune homeostasis and tolerance, with their activation being crucial for preventing various inflammatory responses. However, the mechanisms governing the epigenetic program in T_reg cells during their dynamic activation remain unclear. In this study, we demonstrate that CXXC finger protein 1 (CXXC1) interacts with the transcription factor FOXP3 and facilitates the regulation of target genes by modulating H3K4me3 deposition. Cxxc1 deletion in T_reg cells leads to severe inflammatory disease and spontaneous T-cell activation, with impaired immunosuppressive function. As a transcriptional regulator, CXXC1 promotes the expression of key T_reg functional markers under steady-state conditions, which are essential for the maintenance of T_reg cell homeostasis and their suppressive functions. Epigenetically, CXXC1 binds to the genomic regulatory regions of T_reg program genes in mouse T_reg cells, overlapping with FOXP3 binding sites. Given its critical role in T_reg cell homeostasis, CXXC1 presents itself as a promising therapeutic target for autoimmune diseases.

Introduction

Regulatory T (T_reg) cells are a distinct subset of CD4⁺ T cells that play a critical role in maintaining immune homeostasis and self-tolerance by suppressing excessive or aberrant immune responses to foreign or self-antigens ^1–3. These cells can be further categorized into thymus-derived regulatory T cells (tT_reg cells), periphery-derived T_reg cells (pT_reg cells), and induced T_reg cells (iT_reg cells) ⁴. They uniquely express the transcription factor FOXP3, a member of the forkhead winged-helix family, which is essential for T_reg cell lineage commitment and suppressive function ^5,6. Deletion or mutation of the Foxp3 gene leads to a range of immunological disorders, including allergies, immunopathology, and autoimmune diseases in both mice and humans ^7,8.

It is well established that FOXP3 recruits various cofactors to form complexes that either promote or repress the expression of downstream genes, with histone and DNA modifications playing pivotal roles in this process. FOXP3 can activate or repress the transcription of key regulators of T_reg cell activation and function by recruiting the histone acetyltransferases (HATs) or histone deacetylases (HDACs)^9,10. Notably, FOXP3-bound sites exhibit significant enrichment of H3K27me3, a modification essential for FOXP3-mediated repressive chromatin remodeling under inflammatory conditions¹¹. However, the direct role of FOXP3 as a transcriptional activator through interactions with epigenetic regulators, particularly via modulation of H3K4 trimethylation, remains poorly documented.

In mammals, six proteins have been identified that catalyze H3K4 methylation. These proteins contain the SET domain and include MLL1 (KMT2A), MLL2 (KMT2B), MLL3 (KMT2C), MLL4 (KMT2D), SETD1A, and SETD1B ^12,13. For example, MLL1 plays a critical role as an epigenetic regulator in T_reg cell activation and functional specialization¹⁴. Additionally, Placek et al. demonstrated that MLL4 is essential for T_reg cell development by catalyzing H3K4me1 at distant unbound enhancers through chromatin looping ¹⁵. H3K4me3, which is enriched at the transcription start site (TSS) and the CpG island (CGI), converts chromatin into active euchromatin by recruiting activating factors ¹⁶. CXXC finger protein 1(CXXC1, also known as CFP1), which contains a SET1 interaction domain (SID), is required for binding to the histone H3K4 methyltransferases SETD1A and SETD1B ¹⁷. Previous studies have shown that CXXC1 regulates promoter patterns during T cell maturation, GM-CSF-derived macrophage phagocytosis, TH17 cell differentiation, and the function of ILC3 cells with aging by modulating H3K4me3 modification^18–21. Despite the well-documented role of CXXC1 in various immune effector cells, its role in T_reg cells remains unclear.

Here, we demonstrate that CXXC1 interacts with FOXP3 and enhances the expression of FOXP3 target genes. T_reg cell-specific deletion of Cxxc1 triggers systemic autoimmunity, accompanied by multiorgan inflammation, ultimately resulting in early-onset fatal inflammatory disease in mice. Cxxc1-deficient T_reg cells exhibit a significant disadvantage in proliferation and homeostasis, even in non-inflammatory mice where coexisting wild-type (WT) T_reg cells were present. Moreover, Cxxc1-deficient T_reg cells display intrinsic defects in the expression of key suppression molecules, including CTLA-4, CD25, ICOS, and GITR. Consistent with these findings, mechanistic studies strongly suggest that CXXC1 functions as an essential cofactor of FOXP3, playing a crucial role in sustaining H3K4me3 modifications at key T_reg cell genes.

Materials and methods

Mice

All mice used in this study were bred for a minimum of seven generations on a C57BL/6 background. Mouse experiments, or cells from mice of the same genotype, compared littermates or age-matched control animals. The Cxxc1^fl/fl mouse strain has been previously described ²². The Foxp3^cre mice (JAX,016959) were generously provided by Bin Li (Shanghai Jiao Tong University School of Medicine, Shanghai, China). CD45.1 (NM-KI-210226) mice were purchased from the Nanjing Biomedical Research Institute of Nanjing University. Rag1^-/- mice (stock#T004753) were purchased from GemPharmatech. 2D2 (MOG35-55-specific TCR transgenic) mice were graciously supplied by Prof. Linrong Lu (Zhejiang University School of Medicine, Hangzhou, Zhejiang, China). Foxp3^cre (WT) and Foxp3^cre Cxxc1^fl/fl (cKO) were utilized at 3 weeks old unless otherwise stated. All mice were housed in the Zhejiang University Laboratory Animal Center under specific pathogen-free conditions, and all animal experimental procedures were approved by the Zhejiang University Animal Care and Use Committee (approval no.ZJU20230246).

Cell culture

HEK 293T cells were obtained from ATCC, and Plat E cells were kindly provided by Prof. Xiaolong Liu (Shanghai Institutes for Biological Sciences). Both cell lines were cultured in Dulbecco-modified Eagle medium (DMEM) containing 10% (v/v) fetal bovine serum (FBS), supplemented with 1% penicillin/streptomycin.

Immunofluorescence (IF) Staining

As previously described ²³, coverslips were treated with a 0.01% poly-L-lysine solution (P4707; Sigma) for 10 minutes, air-dried, and then coated with CD4⁺ YFP⁺ T_reg cells. The cells were then fixed in 4% formaldehyde for 15 minutes at room temperature, permeabilized with 0.2% Triton X-100, and blocked with 1% BSA. Antibodies against CXXC1 (ab198977; Abcam) and FOXP3 (17-5773-82; Invitrogen) were diluted in Image iT FX signal enhancer (I3693; Invitrogen) and incubated with cells overnight at 4 °C. After washing with phosphate-buffered saline (PBS), the cells were incubated with a goat anti-rabbit antibody Alexa Fluor 594 (1:250;10015289; Invitrogen) secondary antibody and stained with DAPI (200 ng/mL; D523; Dojindo). Slides were washed with PBS and sealed with an antifade solution(P36934; Invitrogen) before imaging with an Olympus FV3000 fluorescence microscope. The images were visualized using the FV31-SW software.

Co-immunoprecipitation and western blot

Harvest the appropriate transfected cell lines and primary cells from culture and wash them with ice-cold PBS. Lyse the cells in NETN300 buffer (300 mM NaCl, 0.5 mM EDTA, 0.5% (v/v) NP-40, 20 mM Tris-HCl pH 8.0) supplemented with a protease inhibitor (1:100, P8340, Sigma-Aldrich) and PMSF (1 mM) on ice for 10 minutes. Take a small portion of the whole-cell lysate as input, and incubate the remaining lysate with either Anti-FLAG M2 Beads (M8823; Sigma) or Anti-HA Beads (HY-K0201; MCE) on a rotator at 4°C overnight. For the endogenous Co-IP assay targeting CXXC1 and FOXP3, incubate the cell lysate with protein G magnetic beads along with anti-CXXC1 (ab198977; Abcam) or anti-FOXP3 (14-4774-82; Invitrogen) antibodies on a rotator at 4°C overnight. Wash the beads three times with IP buffer (100 mM NaCl, 0.5 mM EDTA, 0.5% (v/v) NP-40, 20 mM Tris-HCl pH 8.0) to remove non-specific binding. Boil the washed beads with 1× Laemmli sample buffer (1610747; Bio-Rad) to elute the bound proteins. Separate the denatured proteins by SDS-PAGE. Transfer the separated proteins onto PVDF membranes (Millipore) for immunoblotting.Immunoblot the PVDF membranes (IPVH00010) with the following antibodies: anti-CXXC1 (1:1000; ab198977; Abcam), anti-FOXP3 (1:500; 14-7979-80; Invitrogen), anti-FlAG (1:1000; 14793; Cell Signaling Technology), anti-HA (1:1000; 3724S; Cell Signaling Technology).Detect the immunoblotted proteins using a secondary HRP-conjugated goat anti-rabbit antibody (1:1000; HA1001-100; Huabio) and visualize the bands using an appropriate detection method.

Flow cytometry

Spleen and lymph nodes were harvested from 3-week-old mice. Tissues were placed in a dish containing PBS with 2% fetal bovine serum (FBS). The tissues were minced into small pieces and then filtered through a 70 µM nylon filter to obtain a single-cell suspension. For surface marker analysis, cells were incubated for 15 minutes with purified anti-mouse CD16/32 antibody (101320; BioLegend) to block Fc receptors. After blocking, cells were stained with the indicated antibodies for surface markers. To determine cytokine expression, cells were stimulated for 4 hours at 37℃ with phorbol12-myristate13-acetate (50ng/mL; S1819; Beyotime), ionomycin (1mg/ml; S1672; Beyotime), and brefeldin A (BFA;00-4506-51; Invitrogen). After stimulation, cells were labeled with a fixable viability dye and surface markers. Cells were then fixed and permeabilized according to the manufacturer’s instructions (00-8222-49; Invitrogen). For transcription factor staining, samples were fixed using the FOXP3/Transcription Factor Staining Buffer Set (00-5523-00; Invitrogen). Flow cytometry was performed using a BD Fortessa (BD Biosciences). Flow cytometry data were acquired and analyzed using FlowJo software.

The following antibodies were purchased from Invitrogen or BioLegend: Zombie Violet f-ixable viability(423113), Zombie NIR fixable(423105), CD25(PC61), CD8α(53-6.7), CD62L(MEL-14), PD-1(J43), CD44(IM7), IL-17A(TC1118H10), KL-RG1(2F1), CD4 (GK1.5), TCRβ(H57-597), IFN-γ(XMG1.2), FOXP3(FJK-16s), CTLA-4(UC10-4B9), ICOS(15F9), GITR(DTA-1), LAG3(C9B7W), TCRVβ11(RR3-15), CD45RB(C363-16A) CD69 (H1.2F3), MHCII(M5/114.15.2), CD80(16-10A1), T-bet(eBio4B10), IKZF4(ESB7C2), IL-4(11B11), CCR7(4B12), CD73(TY/11.8).

Real-time PCR

Total RNA was extracted from T_reg cells using the RNAAiso Plus(9109; Takara) reagent according to the manufacturer’s instructions, and cDNA synthesis was performed using the Prime Script RT Reagent Kit (Takara). TB Green Premix Ex Taq (RR420A; Takara) was used for quantitative real-time PCR (qPCR). The expression levels of target mRNA were normalized to the level of β-actin expression. The primers for qPCR are as follows:

Cxxc1 qPCR Forward: ATCCGGGAATGGTACTGTCG

Cxxc1 qPCR Reverse: CTGTGGAGAAGATTTGTGGG

β-actin qPCR Forward: CTGTCCCTGTATGCCTCTG

β-actin qPCR Reverse: ATGTCACGCACGATTTCC

CD4⁺T and YFP⁺ T_reg cells adoptive transfer in Experimental autoimmune encephalomyelitis (EAE)

CD4⁺YFP⁺ T_reg cells from Foxp3^cre and Foxp3^cre Cxxc1^fl/fl mice were enriched using the Mouse CD4 T Cell Isolation Kit (480005; Biolegend) and then sorted using the BD Aria II flow cytometer. Naïve CD4⁺T cells from 2D2 (MOG35-55-specificTCR transgenic) mice were isolated using the Mouse CD4 Naïve T cell Isolation Kit480039 (480039; Biolegend). As previously described²⁴, 2D2 naïve CD4⁺ T cells alone (5 x 10⁵ per mouse), or 2D2 naïve CD4⁺ T cells (5x 10⁵ per mouse) together with WT or Foxp3^creCxxc1^fl/fl T_reg cells (2x10⁵ per mouse), were transferred into Rag1^−/− mice via the tail vein. One day after cell transfer, the recipient mice were inoculated subcutaneously (s.c.) with 200μg MOG35-55 peptide (MEVGWYRSPFSRVVHLYRNGK; GenemeSynthesis) emulsified in complete Freund’s adjuvant (CFA) (F5506; Sigma). Intravenous administration of 200 ng of Pertussis toxin (181; List Biological Laboratories) was performed on days 0 and 2 after peptide inoculation. The severity of EAE was monitored and blindly graded using a clinical score from 0 to 5: 0, no clinical signs; 1, limp tail; 2, paraparesis (weakness, incomplete paralysis of one or two hind limbs); 3, paraplegia (complete paralysis of two hind limbs); 4, paraplegia with forelimb weakness or paralysis; 5, dying or death.

Isolation lymphocytes from the central nervous system (CNS)

On day 14 after EAE induction, mice were perfused with transcardially administered PBS to eliminate contaminating blood cells in the central nervous system (CNS). The forebrain and cerebellum were dissected to expose the spinal cord, which was then carefully removed from the spinal canal. The fresh spinal cord was harvested and cut into 2 mm pieces. The CNS tissue pieces were homogenized using a syringe and passed through a 70 µM cell strainer to obtain a single-cell suspension. The single-cell suspension was digested with collagenase D (2μg/ml;11088858001; Roche) and deoxyribonuclease I (DNase I; 1 μg/ml; DN25; Sigma-Aldrich) at 37°C for 20 minutes under rotation. After digestion, the cell suspension was centrifuged to pellet the cells. The cell pellets were resuspended in 40% Percoll and layered onto a discontinuous Percoll gradient. Centrifugation at 80% Percoll allowed for the separation of cells at the 40–80% Percoll interface, which were collected as CNS mononuclear cells. The collected CNS mononuclear cells were washed with PBS to remove any remaining Percoll. CNS mononuclear cells were stimulated for 4 hours with PMA and ionomycin in the presence of Brefeldin A to induce cytokine production. After stimulation, cells were fixed, rendered permeable, and stained with appropriate antibodies for intracellular cytokine detection.

Histological analyses

The lungs, skin, liver, and colon were excised from three-week-old mice. Prior to histological analysis, the samples were fixed in formalin, embedded in paraffin, and stained with hematoxylin and eosin. For CNS histology, spinal cords were fixed in 4% paraformaldehyde, paraffin-embedded, sectioned, and stained with Luxol Fast Blue and hematoxylin and eosin (H&E). To examine colon histology, colons from Rag1^−/− hosts were similarly processed and stained with H&E.

Adoptive transfer colitis model

Colitis was induced following the protocol described²⁵. In brief, CD4⁺ YFP⁺ T_reg cells were isolated from 3-week-old CD45.2⁺Foxp3^creCxxc1^fl/fl and CD45.2⁺ Foxp3^cre mice. A total of 2x10⁵ T_reg cells from each group were mixed with 4x10⁵ T_eff cells (CD45.1⁺CD4⁺CD45RB^hi) sorted from CD45.1⁺ mice and transferred into the Rag1^-/- mice via intraperitoneal injection. T_eff cells alone were transferred as a control group. Mouse body weight was measured weekly post-adoptive transfer. The percentage change in body weight was calculated by comparing the current weight with the initial weight on day 0. Mice were euthanized when any had reached 80% of their initial body weight. The large intestines were sectioned into 4μm thick slices and stained with hematoxylin.

In vitro T_reg suppression assay

Naïve CD4⁺ T cells isolated from WT mice were labeled with CFSE (C34554; Invitrogen). CD4⁺YFP⁺ T_reg cells from Foxp3^cre and Foxp3^creCxxc1^fl/fl mice were cultured with naïve CD4⁺ T cells (1ⅹ10⁵ cells) at various ratios in the presence of 2μg/mL anti-CD3(16-0031-85; Invitrogen) and 3μg/mL anti-CD28(16-0281-85; Invitrogen). On day 3, cells were analyzed by flow cytometry.

CUT&Tag

CUT&Tag assays of CD4⁺YFP⁺ T_reg cells were conducted as previously described²⁶. Briefly, approximately 1x10⁵ single cells were carefully pipetted into wash buffer twice. The pelleted cells were resuspended in wash buffer, activated concanavalin (BP531; Bangs Laboratories), and incubated for 15 minutes at room temperature. Cells bound to the beads were resuspended in Dig-Wash Buffer and incubated with a 1:50 dilution of primary antibodies (rabbit anti-H3K4me3, Active Motif,39016; rabbit anti-CXXC1, abcam, ab198977; normal IgG, Cell Signaling, 2729) at 4 ℃ overnight. The beads were incubated with a secondary antibody (goat anti-rabbit IgG; SAB3700883; Sigma-Aldrish) diluted 1:100 in Dig-Wash buffer for 60 minutes at room temperature. Cells were treated with Hyperactive pG-Tn5 Transposase (S602; Vazyme) diluted in Dig-300 Buffer for 1 hour at room temperature. The cells were subsequently resuspended in Tagmentation buffer (10mM MgCl2 in Dig-300 Buffer) and incubated at 37℃ for 1 hour. To halt tagmentation, 10μl was spiked with 0.5M EDTA, 3μl with 10%SDS, and 3μl with 20mg/ml Proteinase K and incubated at 55℃ for 1h. DNA library amplification was performed according to the manufacturer’s instructions and purified using VAHTS DNA Clean Beads (N411; Vazyme). Libraries were sequenced on the Illumina NovaSeq platform (Annoroad Gene Technology).

CUT&Tag and ChIP-seq data analysis

FOXP3 ChIP-seq data was obtained from GSE121279. H3K27me3 ChIP-seq data was obtained from GSE14254. CUT&Tag and ChIP-seq reads were trimmed to 50 bp and aligned against the mouse genome build mm9 using Bowtie2(v2.3.4.1) with default parameters. All PCR duplicates and unmapped reads were removed. Peak calling was performed using MACS2(v2.1.1.20160309) and signal tracks for each sample were generated using the ‘wigToBigWig’ utility of UCSC. We classified the H3k4me3 peaks around TSSs into three groups: broad (>5 kb), medium (1-5kb), and narrow (<1 kb). The top 5% of the widest peaks were considered as broad peaks. The average intensity profiles were generated using deepTools (v2.5.4). Motif analysis was performed using the ‘findmotifsGenome.pl’ command inHomer2 package. Epigenetic factors were identified using the Epigenetic Factor Database(https://epifactors.autosome.org/) and then screened for those that exclusively regulate the expression of their target genes by modulating the deposition of H3K4me3.

Genomic distribution was analyzed using the “genomation” R package. GO pathway analysis was performed using the “clusterProfiler” R package. The sequencing information of CUT&Tag data used in this study is summarized in Supplementary Table S1.

Clustering analysis

Promoters were defined as ±2 kb regions flanking the annotated TSS. Reads in promoters were counted using the “coverage” command in bedtools (v2.26.0) and further normalized to RPKM. The k-means clustering of H3K4me3 and H3K27me3 enrichment at promoters was conducted using the “kmeans” function in R.

RNA-seq and data analysis

Total RNA was extracted from sorted CD4⁺YFP⁺ T_reg cells using the RNeasy Plus Mini Kit (Qiagen, #74134), following the manufacturer’s protocol. RNA-Seq libraries were constructed and sequenced by Haplox (Nanchang, China), using an Illumina platform with paired-end reads of 150 bp. RNA-seq data was obtained from GSE82076. Raw reads were trimmed to 50 bp and mapped to the mouse genome (mm9) using TopHat (v2.1.1) with default parameters. Only uniquely mapped reads were kept for downstream analysis. The RNA abundance of each gene was quantified using Cufflinks (v2.2.1).

Single-Cell RNA-Sequencing

A total of 300,000 sort-purified CD4⁺YFP⁺ T_reg cells from Foxp3^cre and Foxp3^cre Cxxc1^fl/fl were resuspended in BD Pharmingen Stain Buffer (FBS) (554656; BD). Single cells were isolated using a chromium controller (BD platform, BD Bioscience) according to the manufacturer’s instructions, as previously described²⁷. The single cells were labeled with sample tags using the BD Mouse Immune Single-Cell Multiplexing Kit (633793; BD). Following standard protocols, cDNA amplification and library construction were performed to generate scRNA-seq libraries.

Targeted scRNA-seq data processing

The raw FASTQ files were processed by BD Rhapsody using the Targeted analysis pipeline. After alignment and filtering, the distribution-based error correction(DBEC)-adjusted molecules were loaded into R studio(version 4.3.2). All subsequent analyses were performed using the package Seurat (version 4.4.0) with default parameters. Specifically, the scRNA-seq data counts were log-normalized. All targeted genes were scaled and then were used for Principal components analysis (PCA). The batch effects were removed by the HarmonyMatrix function in the Harmony package (version 1.1.0). The first 20 principal components were used to calculate nonlinear dimensionality reduction using RunUMAP. Differential gene expression (DEGs) between clusters was assessed using the FindAllMarkers function. The clusters were then annotated based on DEGs. Barplots were generated using ggplot2 (version 3.4.4). Heatmaps were generated using pheatmap (version 1.0.12).

Analysis of the single-cell TCR-seq repertoire

Raw V(D)J fastq reads were processed using BD Rhapsody Pipeline and then were analyzed using scRepertoire (version 1.12.0). The TCR clonotype was called using the nucleotide sequence of the CDR3 region for both TCR alpha and beta chains. For cells with multiple chains, the top two clonotypes with the highest expression were selected for downstream analysis. Clonal overlap between different cell types was calculated using the clonalOverlap function of scRepertoire. A clonotype was defined as expansion if it could be detected in at least two cells.

scRNA-seq trajectory analysis

UMAP embeddings obtained from the Seurat package were projected into the Slingshot(version 2.10.0) package to construct pseudotime Trajectories for T_reg cells. Naïve subsets were set as the root state.

Whole genome bisulfite sequencing (WGBS) and data analysis

Sorted CD4⁺YFP⁺ T_reg cells(3ⅹ10⁶) were lysed in cell lysis buffer to release DNA. The bisulfite-treated DNA was used to prepare the sequencing library. DNA libraries were transferred to the Illumina Platform for sequencing using 150 bp paired-end reads. Raw reads were trimmed using TrimGalore (v0.4.4) with default parameters. Subsequently, the reads were mapped against the mm9 reference genome using Bismark v0.19.0 with parameters ‘--bowtie2’. PCR duplicates were removed and the methylation levels were calculated using ‘bismark_methylation_extractor’. We calculated the mean CpG methylation levels of various genome elements: promoter, 5′-UTR, exon, intron, 3′-UTR, genebody, intergenic, CGIs, and repeats using in-house scripts. The sequencing information of WGBS data used in this study is summarized in Supplementary Table S2.

Statistical analysis

The statistical significance analysis was performed using Prism 8.0 (GraphPad). Error bars are presented as mean ± SD. P values of < 0.05 were deemed statistically significant (*P < 0.05, **P < 0.01, ***P < 0.001 and ****P < 0.000). Statistical analyses were performed with unpaired Student’s t-test or two-way ANOVA and Holm–Sidak post hoc test.

Results

FOXP3 binds regulatory loci primed for activation and repression in T_reg cells

FOXP3-mediated gene expression is well-recognized, with several studies highlighting its dual role as both a transcriptional activator and repressor ^28–32. However, the regulation of T_reg-specific gene transcription by FOXP3 through intermediary mechanisms, such as epigenetic factors, remains incompletely understood. To address this, we employed CUT&Tag to generate global H3K4me3 maps in T_reg cells and compared these with an H3K27me3 ChIP-Seq dataset from Wei et al.³³ at previously identified FOXP3-bound loci(Konopacki et al. 2019)³⁴. This comparison aimed to identify FOXP3-dependent genes associated with H3K4me3 and H3K27me3 deposition. As expected, H3K4me3 was enriched at gene promoters (fig. S1, A and B). A Venn diagram revealed a significant overlap between FOXP3 binding sites and H3K4me3 peaks, with minimal overlap with H3K27me3 peaks (Fig. 1, A and B). The overlapping regions between FOXP3 binding sites and H3K4me3 or H3K27me3 peaks were predominantly located at promoters (Fig. 1, A and B). We then clustered the promoters into four groups based on the enrichment of H3K4me3 and H3K27me3. Clusters 1 and 3 showed strong enrichment of H3K4me3; cluster 2 was enriched with H3K27me3; and cluster 4 showed weak enrichment of both modifications. FOXP3 enrichment analysis revealed that FOXP3 preferentially binds to the promoters of clusters 1 and 3, which display high levels of H3K4me3 (Fig. 1C). Correspondingly, genes in these clusters exhibited high transcription levels, as demonstrated by our reanalysis of previously published RNA-seq data (Oh et al. 2017) shown in the right panels of Figure 1D³⁵. In contrast, genes with H3K27me3 enrichment at their promoters were transcribed at low levels. Gene Ontology (GO) analysis of these four clusters revealed that cluster 1 was significantly enriched in genes involved in mRNA processing, covalent chromatin modification, and histone modification, while cluster 3 was enriched in genes related to DNA repair and mitochondrion organization (Fig. 1E). Cluster 2, on the other hand, was enriched in the pattern specification process, whereas cluster 4 showed no significant correlation with T_reg cells. As expected, signature T_reg cell genes such as Tnfrsf18, Nrp1, Stat5a, Lag3, Icos, and Pdcd1 were enriched in clusters 1 and 3, which exhibited strong H3K4me3 marks (fig. S1C). Conversely, FOXP3-bound sites in cluster 2 were specifically enriched for H3K27me3, including genes like Hic1, Trp73, and Rnf157, which are associated with inflammatory responses(fig. S1C). These observations support the hypothesis that FOXP3 primarily functions as an activator by promoting H3K4me3 deposition in T_reg cells.

H3K4me3 is required for FOXP3-dependent gene activation in T_reg cells.
(A)Venn diagram showing overlap of H3K4me3 enriched regions (this study) and FOXP3 binding sites (Konopacki *et al*.2019) in sorted CD4⁺YFP⁺ T_reg cells (left). Genomic distribution of overlapped peaks(right). Note that the overlapped peaks are predominantly enriched at promoters.
(B)Venn diagram showing overlap of H3K27me3 enriched regions (Wei *et al*.2009) and FOXP3 binding sites in T_reg cells (left). Genomic distribution of overlapped peaks (right). Note that the overlapped peaks are predominantly enriched at promoters.
(C)Heatmap showing enrichment of H3K27me3, H3K4me3, and FOXP3 surrounding the TSS. Unsupervised K-means clustering was conducted on H3K27me3 and H3K4me3 signals.
(D)Heatmap showing gene expression levels in T_reg cells (RNA-seq data was obtained from Oh *et al.*2017) The clusters were consistent as in Fig. 1C.
(E)Gene Ontology (GO) pathway analysis of different clusters.

CXXC1 interacts with FOXP3 and binds H3K4me3-enriched sites in T_reg cells

We conducted an enrichment analysis of known motifs at the overlapping peaks of FOXP3 ChIP-seq and H3K4me3 CUT&Tag in T_reg cells to identify epigenetic factors that directly interact with FOXP3 to mediate chromatin remodeling and transcriptional reprogramming. Motif analysis of the overlapping peaks between FOXP3 binding sites and regions enriched in H3K4me3 revealed that, in addition to transcription factors, the most abundant motif associated with H3K4me3 was the epigenetic factor CXXC1(fig. S2A). To investigate this further, we performed CUT&Tag for endogenous CXXC1 in T_reg cells to examine the genome-wide co-occupancy of CXXC1 and FOXP3. Over half of these CXXC1 binding sites were located at promoter regions (fig. S2B). Additionally, CXXC1 exhibited strong binding at transcription start site (TSS) and CpG islands (CGIs) (Fig. 2A and fig. S2C). As illustrated by the Venn diagram (Fig. 2B), more than half of the FOXP3-bound genes and H3K4me3-enriched genes were also bound by CXXC1. Similarly, more than half of CXXC1 peaks were overlapped with FOXP3 peaks (fig. S2D). Furthermore, the CXXC1-specific and FOXP3-specific binding sites also demonstrated modest binding of FOXP3 and CXXC1, respectively (Fig. 2C). These findings indicate that FOXP3 and CXXC1 share a substantial number of target genes in T_reg cells. To confirm this interaction, we further validated the reciprocal immunoprecipitation of both endogenous and exogenous CXXC1 and FOXP3 (Fig. 2D and fig. S2E). An immunofluorescence assay revealed predominant colocalization of CXXC1 with FOXP3 in the nucleus (Fig. 2E). Overall, these results suggest that CXXC1 primarily functions as a coactivator of FOXP3-driven transcription in T_reg cells.

CXXC1 interacts with FOXP3 in T_reg cell.
(A)Average CXXC1 CUT&Tag signals around promoters in T_reg cells. IgG was used as the control.
(B)Venn diagrams showing the overlap of FOXP3 binding genes, CXXC1 binding genes, and H3K4me3 enriched genes in T_reg cells. Gene promoters covered by FOXP3 binding sites, CXXC1 binding sites, or exhibited high H3K4me3 levels were defined as FOXP3-bound genes, CXXC1-bound genes, or H3K4me3-enriched genes.
(C)Heatmaps showing FOXP3 ChIP-seq and CXXC1 CUT&Tag signals at indicated regions.
(D)Interaction between FOXP3 and CXXC1 was assessed by co-IP (forward and reverse) using T_reg cell lysates.
(E)Immunofluorescence for FOXP3 and CXXC1 colocalization in T_reg cells. Scale bars, 2 μm.

Complete ablation of Cxxc1 in T_reg cells leads to a fatal autoimmune disease

To investigate the role of CXXC1 in T_reg-cell homeostasis and function, we generated Foxp3^CreCxxc1^fl/fl mice (conditional knockout [cKO] mice) by crossing Cxxc1^fl/fl with Foxp3^YFP-Cre ³⁶ mice, thereby specifically deleting Cxxc1 in T_reg cells. The effective depletion of Cxxc1 in T_reg cells was confirmed through quantitative PCR (qPCR) and Western blotting (fig. S3A). Notably, cKO mice appeared normal at birth but later exhibited spontaneous mortality starting around three weeks of age (Fig. 3A). Deletion of Cxxc1 in T_reg cells led to the development of severe inflammatory disease, characterized by reduced body size, stooped posture, crusting of the eyelids, ears, and tail, and skin ulceration, particularly on the head and upper back (Fig. 3, B and C). Additionally, cKO mice developed extensive splenomegaly and lymphadenopathy (Fig. 3D). Histopathological analysis revealed massive lymphocyte and myeloid cell infiltration in the skin, lungs, liver sinusoids, and colon mucosa (Fig. 3E). Moreover, cKO mice displayed a significant increase in CD8⁺ T cell percentages (fig. S3B), along with a marked rise in cells exhibiting an effector/memory phenotype (CD44hi CD62Llo) (Fig. 3F). Furthermore, T cells from cKO mice produced elevated levels of IFN-γ, IL-17, and IL-4 in CD4⁺ T cells, as well as increased IFN-γ production in CD8⁺ T cells (Fig. 3G and fig. S3, C and D). These phenotypes closely resembled those observed in Foxp3-deficient mice ³⁷ or mice with depleted T_reg cells³⁸, suggesting a deficiency in immune suppression.

*Foxp3*^creCxxc1^fl/fl mice spontaneously develop a fatal early-onset inflammatory disorder.
(A)Survival curves of *Foxp3*^cre (black line) and *Foxp3*^creCxxc1^fl/fl (red line) mice (n = 20).
(B)Gross body weight of *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice (n = 6).
(C)A representative image of *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(D)Representative images showing the spleen and peripheral lymph nodes from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(E)Hematoxylin and eosin staining of the skin, lung, liver, and colon from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice (scale bar, 100 μm).
(F)Flow cytometry analysis of CD62L and CD44 expression on peripheral lymph node CD4⁺ and CD8⁺T cells from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice (left). Right, frequency of CD44^lowCD62L^hi and CD44^hiCD62L^low population in CD4⁺ or CD8⁺ T cells (n = 9).
(G)Lymph node cells from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice were stimulated ex vivo with PMA + ionomycin for 4h and analyzed for IFN-γ expressing in CD4⁺ or CD8⁺ T cells using flow cytometry (left). Right, percentages of IFN-γ⁺CD4⁺ or IFN-γ⁺CD8⁺ T cells in the lymph nodes of *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice (n = 7).
All mice analyzed were 18-20 days old unless otherwise specified. Error bars show mean ± SD. The log-rank survival curve was used for survival analysis in A, and multiple unpaired t-test or two-tailed Student’s t-test were used for statistical analyses in B, F, and G (****P<0.0001). The flow cytometry results are representative of three independent experiments.

CXXC1 is necessary for the maintenance of T_reg-cell suppressive activity

Despite the development of severe autoimmune disease, we observed an increase in both the absolute number and percentage of FOXP3⁺ T_reg cells in the lymph nodes (fig. S4A). The expression level of the FOXP3 protein was only slightly altered in Cxxc1-deficient T_reg cells (fig. S4B). In an in vitro suppression assay, T_reg cells from Foxp3^CreCxxc1^fl/fl and WT mice exhibited similar suppressive effects on naïve T (Tn) cell proliferation (fig. S4C). The expression of the hallmark T_reg-cell marker CTLA-4 showed a modest increase in Cxxc1-deficient T_reg cells compared to WT T_reg cells, while the expression of GITR remained unchanged (fig. S4D). To further assess the suppressive capacity of Cxxc1-deficient T_reg cells in vivo, we employed the experimental autoimmune encephalomyelitis (EAE) model. Naïve CD4⁺ T cells from 2D2 mice were co-transferred with T_reg cells from either Foxp3^Cre or Foxp3^Cre Cxxc1^fl/fl mice into Rag1^-/- recipients, and EAE was induced in these recipient mice. Mice that received only naïve CD4⁺ T cells from 2D2 mice developed more severe EAE symptoms (Fig. 4A). The addition of WT T_reg cells from Foxp3^Cre mice slightly mitigated EAE progression and reduced Th17 cells in the spinal cord (Fig. 4, A to D). In contrast, Foxp3^Cre Cxxc1^fl/fl T_reg cells failed to suppress EAE (Fig. 4, A to D), and the cKO mice showed a significant reduction in T_reg cells frequency in CNS tissues(Fig. 4E). Finally, we examined the role of CXXC1 in T_reg cell-mediated suppression using T cell transfer-induced colitis, in which naïve T cells were transferred to Rag1^-/- recipients either alone or together with WT or Foxp3^Cre Cxxc1^fl/fl T_reg cells. The transfer of naïve T cells led to weight loss and intestinal pathology in recipient mice(Fig. 4, F and G). Mice receiving WT T_reg cells continued to gain weight (Fig. 4F), whereas those that received T_reg cells from cKO mice were unable to prevent colitis and exhibited a reduced percentage of T_reg cells (Fig. 4, F to H). These findings underscore the critical role of CXXC1 in maintaining T_reg cell function in vivo.

CXXC1 is essential for T_reg cells to suppress T cell-mediated EAE and colitis.
(A)Mean clinical scores for EAE in *Rag1^−/−* recipients of 2D2 CD4⁺ T cells, either alone or in combination with *Foxp3*^cre or *Foxp3*^creCxxc1^fl/fl mice after immunization with MOG35–55, complete Freund’s adjuvant (CFA), and pertussis toxin (n=5).
(B and C)Representative histology of the spinal cord of *Rag1^−/−* mice after EAE induction. Hematoxylin and eosin(H&E) staining (upper), Luxol fast blue (F&B) staining (lower). Scale bars, 50 μm (400×) and 200 μm (100×).
(D)Representative flow cytometry plots and quantification of the the percentages of IFNγ⁺ or IL-17A⁺ CD4⁺Vβ11⁺ T cells (n=5).
(E)Statistical analysis of the percentage CD4⁺ FOXP3⁺ T_reg cell in CNS tissues 14 days after EAE induction.
(F)Changes in body weight of *Rag1^−/−* mice after colitis induction(n=6).
(G)Haematoxylin and eosin (H&E) staining of colons from T cell-induced colitis mice 6 weeks after T cell transfer. Scale bars, 50 μm (400×) and 200 μm (100×).
(H)Statistical analysis of the percentage CD4⁺ FOXP3⁺ T_reg cell in the spleen, mesenteric lymph nodes, and colon 6 weeks after colitis induction.
Error bars show mean ± SD. P values are determined by a two-tailed Student’s t-test or two-way ANOVA and Holm–Sidak post hoc test (A, D, and E). (*P<0.05, **P<0.01, ***P < 0.001, ****P<0.0001).

T_reg cell lineage stability and proliferation depend upon CXXC1

T_reg cells harbor a diverse T-cell receptor (TCR) repertoire, which likely play a critical role in their immune suppression function ^39–41. To explore the role of CXXC1 in T_reg-mediated suppression, we performed single-cell RNA sequencing (scRNA-seq) combined with TCR sequencing (TCR-seq) on CD4⁺YFP⁺ T_reg cells isolated from mouse lymph nodes. After quality control and removal of doublets, 18,577 cells were retained for further analysis. Through unsupervised clustering and uniform manifold approximation and projection (UMAP) analysis, we identified eight distinct T_reg cell clusters based on the expression of well-characterized markers, with a particular focus on two clusters of activated T_reg cells that exhibited higher expression of markers and gene sets relative to naïve T_reg cells (Fig. 5A and fig. S5, A to C). A comparison between Cxxc1-deficient and WT T_reg cells within each cluster revealed a reduction in Cxxc1-deficient cells in the naïve subsets, while an increase was observed in the Gzmb⁺ and H2-Eb1⁺ subsets (Fig. 5B). To further elucidate the transition of T_reg cells along a dynamic biological timeline, we constructed pseudo-time trajectories using Slingshot ⁴². The pseudo-time gradient depicted a progression from quiescent to activated T_reg cells, ultimately encompassing the Gzmb⁺ and H2-Eb1⁺ subsets (Fig. 5C). Given the antigen-specific suppression capabilities of T_reg cells ^43,44, we examined their clonal expansion. The analysis revealed that expanded WT TCR clonotypes (n ≥ 2) were predominantly distributed among the Nt5e⁺ subsets, while Cxxc1-deficient T_reg cells showed expanded clonotypes primarily within the Gzmb⁺ and H2-Eb1⁺ subsets(Fig. 5D and Fig. S5D). TCR sharing analysis indicated significant clonotype sharing among various clusters of WT T_reg cells, suggesting a degree of homogeneity. However, the reduced TCR sharing in Cxxc1-deficient cells implies that decreased TCR diversity may impair the suppressive activity of T_reg cells (Fig. 5E) ⁴¹. Furthermore, the Cxxc1-deficient group exhibited lower expression of several T_reg-specific genes associated with suppressive functions, such asNt5e, Il10, Pdcd1, Klrg1, as well as genes that inhibit effector T cells differentiation, including Sell and Tcf7 (Fig. 5F and fig. S5E). Conversely, Cxxc1-deficient T_reg cells demonstrated elevated expression of Gzmb, Il2ra, and Cd69 compared to WT, reflecting a profile indicative of increased activation (Fig. 5, F and G). Additionally, we observed increased expression of genes linked to Th1-type inflammation, such as Ifng, Tbx21, and Hif1a, in Cxxc1-deficient T_reg cells, likely due to extreme inflammatory conditions (Fig. 5, F to H). The proportion of FOXP3⁺Ki67⁺ T_reg cells was lower in cKO mice compared to WT mice (Fig. 5I). These findings underscore the crucial role of CXXC1 in maintaining T_reg cell stability.

Single-cell transcriptomics reveals distinct T_reg cell populations.
(A)UMAP plot showing clusters identified based on variable gene expression of sorted CD4⁺YFP⁺ T_reg cells. Each dot represents a cell, and each color corresponds to a different population of cell types. Clustering analysis revealed 8 distinct T_reg cell populations.
(B)Mean fold changes in cluster abundance between *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(C)Pseudotime trajectories of T_reg cells based on Slingshot, color-coded by T_reg cell subpopulations.
(D)Visualization of density and clonotype richness across T_reg clusters from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(E)TCR sharing of expanded clonotypes across all possible combinations of T_reg cells from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(F)Heatmap showing Z scores for the average expression of T_reg-specific genes in each cluster between *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice.
(G-H)Representative flow cytometry plots and quantification of (G) expression of CD25, CD69, ICOS, T-bet, and (H) IFN-γ in CD4⁺YFP⁺ T_reg cells from *Foxp3*^cre and *Foxp3*^creCxxc1^fl/fl mice (n=5 CD25, n=10 CD69, n=5 ICOS, n=6 T-bet, n=5 IFN-γ).
(I)Ki-67 expression (left) and frequency (right) in CD4⁺FOXP3⁺ T_reg cells from *Foxp3*^cre a-nd *Foxp3*^creCxxc1^fl/fl mice (n=6).
Error bars show mean ± SD. P values are determined by a two-tailed Student’s t-test (G-H). (**P<0.01, ***P < 0.001, ****P<0.0001). The flow cytometry results are representative of three independent experiments.

Cxxc1 deletion leads to a reduction in essential T_reg cell markers

To confirm that the deficiency in T_reg cell function in T_reg-specific Cxxc1-deficient animals is due to intrinsic defects caused by Cxxc1 deficiency, rather than severe autoimmune inflammation in Foxp3^Cre Cxxc1^fl/fl mice, we examined Cxxc1-sufficient and Cxxc1-deficient T_reg cell subsets in heterozygous Foxp3^Cre/+ Cxxc1^fl/fl (designated as “het-KO”) and littermate Foxp3^Cre/+ Cxxc1^fl/+ (designated as “het-WT”) female mice (Fig. 6A). Notably, het-KO female mice did not exhibit overt signs of autoimmunity, as random X-chromosome inactivation led to the coexistence of both Cxxc1-cKO and Cxxc1-WT T_reg cells. However, both the frequency and absolute numbers of FOXP3⁺YFP⁺ T_reg cells within the total T_reg population were reduced in het-KO mice compared to their counterparts in het-WT littermates, indicating that Cxxc1 deficiency imposes a competitive disadvantage on T_reg cells (Fig. 6B). Additionally, Cxxc1-deficient YFP⁺ T_reg cells failed to upregulate the proliferation marker Ki-67 (Fig. 6C). Moreover, YFP⁺ T_reg cells in het-KO female mice showed reduced expression of key genes essential for suppressive function, including Icos, Il2ra, Ctla4, and Tnfrsf18, compared to YFP^-T_reg cells from the same mice (Fig. 6D). To investigate the molecular program affected by the deletion of Cxxc1 in T_reg cells, we performed RNA-sequencing (RNA-seq) analysis on CD4⁺YFP⁺ T_reg cells isolated from het-WT and het-KO mice. We then conducted a differential gene expression (DGE) analysis based on the RNA-seq data. Among all expressed genes, 865 were upregulated (“Up” genes) and 761 were downregulated (“Down” genes) by ≥1.5-fold in CD4⁺YFP⁺ T_reg cells from het-KO mice compared to het-WT mice, with a false discovery rate (FDR)–adjusted P-value cutoff of 0.05 (fig. S6A). Gene Ontology (GO) enrichment analysis revealed that the downregulated genes in Cxxc1-deficient T_reg cells were predominantly enriched in pathways related to the negative regulation of immune system process and regulation of cell−cell adhesion (fig. S6B). The Cxxc1-deficient T_reg cells also showed reduced expression of several genes associated with T_reg cell suppressive function, including Il10, Tigit, Lag3, Icos, Nt5e(encoding CD73) and Itgae (encoding CD103) (fig. S6C). Thus, although YFP⁻ WT T_reg cells can prevent autoimmunity in het-KO mice, the absence of Cxxc1 in YFP⁺ T_reg cells leads to the dysregulation of key T_reg cell markers under steady-state conditions.

*Cxxc1*-deficient T_reg cells in het-KO female mice exhibit reduced suppression markers under steady-state conditions.
(A)Schematic representation of wild-type and Cre-positive T_reg cells in female *Foxp3*^Cre/+mice.
(B)Flow cytometry analysis of the YFP^-FOXP3⁺ (WT) and YFP⁺FOXP3⁺ (KO) T_reg cells in *Foxp3*^Cre/+ *Cxxc1*^fl/+ (het-WT) and *Foxp3*^Cre/+ *Cxxc1*^fl/fl (het-KO) female mice (left), along with the frequency and absolute numbers of YFP⁺ cells within the total T_reg population (right)(n=5).
(C)Flow cytometry analysis of Ki-67expression (left) and MFI (right) in YFP⁻ and YFP⁺ cells within the CD4⁺FOXP3⁺ T_reg cells from 6- to 8-week-old het-KO female mice (n=4).
(D)Representative flow cytometry plots and quantification of ICOS, CD25, CTLA4, and GITR expression in YFP⁻ and YFP⁺ cells within CD4⁺FOXP3⁺ T_reg cells from het-KO female mice (n=5).
Error bars show mean ± SD. P values are determined by a two-tailed Student’s t-test (G-H). (**P<0.01, ***P < 0.001, ****P<0.0001). The flow cytometry results are representative of three independent experiments.

The FOXP3-CXXC1 complex regulates the expression of key factors in T_reg cells that are associated with the breadth of H3K4me3

CXXC1 binds to unmethylated CpG DNA via its N-terminal CXXC finger domain, facilitating its interaction with DNA methyltransferase 1 (DNMT1). This binding stabilizes the DNMT1 protein, thereby regulating DNA methylation ^45,46. To investigate whether CXXC1 depletion affects DNA methylation in T_reg cells, we performed whole genome bisulfite sequencing (WGBS) on T_reg cells isolated from both WT and cKO mice. On average, Cxxc1-deficient T_reg cells exhibited no significant changes in DNA methylation at gene loci or across genome-wide CpG sites, irrespective of chromosomal region (fig. S7, A to C). Furthermore, Cxxc1 knockout T_reg cells did not show a significant increase in DNA methylation at key T_reg signature gene loci (fig. S7D). Given the pivotal role of MLL4-mediated H3K4me1 in establishing the enhancer landscape and facilitating long-range chromatin interactions during T_reg cell development¹⁵, we performed CUT&Tag to assess changes in H3K4me1 levels in Cxxc1-deficient T_reg cells. This analysis revealed that H3K4me1 levels were similar in both WT and Cxxc1-deficient T_reg cells (fig. S7, E to G).

While H3K4me3 modifications typically form sharp 1-to 2-kb peaks around promoters, some genes exhibit broader H3K4me3 regions, referred to as broad H3K4me3 domains (H3K4me3-BDs), which can extend to cover part or all of the gene’s coding sequences (up to 20 kb) ^47,48. Broad H3K4me3 domains are preferentially associated with genes essential for the identity or function of specific cell types ^47,49 and have been implicated in enhancing transcriptional elongation and increasing enhancer activity ⁴⁹. To further explore the connection between broad H3K4me3 domains and the expression of immune-regulatory genes, we analyzed genes containing broad H3K4me3 regions. We classified the H3K4me3 domains surrounding transcription start sites (TSSs) into three categories: broad (more than 5 kb), medium (between 1 and 5 kb), and narrow (less than 1 kb) (Fig. 7A). Notably, Cxxc1-deficient T_reg cells exhibited weaker H3K4me3 signals compared to WT cells within the broad H3K4me3 domains where CXXC1 binding is prominent (Fig. 7, A and B). Using the criteria established by Benayoun et al.⁴⁷, which defines the top 5% of the widest H3K4me3 domains as BDs, we observed similar enrichment results (fig. S7, H and I). We then compared three groups of genes: BD-associated genes with reduced H3K4me3 levels following Cxxc1 deletion, genes with direct CXXC1 binding, and genes with direct FOXP3 binding. The Venn diagram revealed that the majority of genes (283 out of 294, 96%) with CXXC1 binding and reduced H3K4me3 levels overlap with FOXP3-bound genes, suggesting that CXXC1 presence enhances FOXP3 binding within the broad H3K4me3 domain (Fig. 7C). Furthermore, GO term analysis indicated that BD-associated genes are enriched in biological processes related to the negative regulation of immune system processes (Fig. 7D). Genome browser views displayed the enrichments of FOXP3, CXXC1, and H3K4me3 at key signature genes in T_reg cells, such as Ctla4, Il2ra, Icos, and Tnfrsf18, with lower H3K4me3 densities observed at these loci in Cxxc1-deficient T_reg cells (Fig. 7E). Similar patterns were observed at core genes involved in T_reg homeostasis and suppressive function (e.g., Lag3, Nt5e, Ikzf4, and Cd28) (fig. S7, J)^50,51. These findings suggest that CXXC1 and FOXP3 collaboratively promote sustained T_reg cell homeostasis and function by preserving the H3K4me3 modification at key T_reg cell genes.

CXXC1 regulates Foxp3-dependent molecule H3K4 trimethylation in T_reg cells.
(A-B)Heatmaps showing H3K4me3(A)and CXXC1(B) signals centered on narrow, medium, and broad domains. The top panel shows the average CUT&Tag signals around indicated domains.
(C)Venn diagrams showing the overlap of FOXP3 binding genes, CXXC1 binding genes, and H3K4me3-BD-associated genes with decreased H3K4me3 levels after CXXC1 depletion in T_reg cells.
(D)Gene Ontology (GO) pathway analysis of the overlap peaks.
(E)Representative genome browser view showing the enrichments of FOXP3, CXXC1, and H3K4me3 in T_reg cells.

Discussion

T_reg cells specifically express the transcription factor FOXP3, which is essential for maintaining T_reg lineage stability and suppressive function⁵². However, FOXP3 alone is insufficient to fully regulate the transcriptional signature and functionality of T_reg cells; its interaction with protein partners is crucial for this regulation. In this study, we identify CXXC1 as a key transcription factor that, by acting as a critical cofactor of FOXP3, orchestrates the T_reg transcriptional program, which in turn plays a vital role in T_reg cell homeostasis and function.

FOXP3 is known to promote histone H3 acetylation at the promoters and enhancers of its target genes, such as Il2ra, Ctla4, and Tnfrsf18, following T_reg cell activation, thereby functioning as a transcriptional activator^53,54. Conversely, FOXP3 can act as a repressor by silencing target genes like Il2 and Ifng through the induction of histone H3 deacetylation, mediated by the recruitment of histone deacetylases and transcriptional co-repressors^10,55. Additionally, FOXP3 exerts this repression by recruiting the Ezh2-containing polycomb repressive complex to target genes during activation, as FOXP3-repressed genes are associated with H3K27me3 deposition and reduced chromatin accessibility^11,56. Recent studies have begun to explore the biological significance of H3K4me3 breadth, revealing a positive correlation between H3K4me3 breadth and gene expression ^49,57,58. These studies also suggest that H3K4me3 breadth contributes to defining specific cell identities during development and disease, including systemic autoimmune diseases like systemic lupus erythematosus and various cancers ^47,59–62. Our study demonstrates that the loss of Cxxc1 leads to reduced H3K4me3 levels, predominantly in genes with broader peaks, such as Il2ra, Tnfrsf18, Ctla4 and Icos, directly impairing their suppressive function in T_reg cells. Although our research primarily focused on the role of CXXC1 in T_reg cells, it is plausible that similar mechanisms may be operative in other cell types.

The T_reg-specific deletion of Cxxc1 leads to a rapid and fatal autoimmune disorder, characterized by systemic inflammation and tissue damage, underscoring the essential role of CXXC1 in maintaining immune self-tolerance within T_reg cells. Interestingly, despite this severe phenotype, our findings show that H3K4me1 levels were comparable between WT and Cxxc1-deficient T_reg cells. In contrast, MLL4 is critical for T_reg cell development in the thymus, primarily by regulating H3K4me1, though it is not required for peripheral T_reg cell function¹⁵. Numerous studies investigating the effectors of T_reg cell-mediated suppression have identified a wide range of molecules and mechanisms. These include the upregulation of the inhibitory co-stimulatory receptor CTLA-4, which initiates inhibitory signaling; the sequestration of the T-cell growth factor IL-2 via CD25; the secretion of inhibitory cytokines, such as interleukin (IL)-10, IL-35 and TGF-β and the activity of ectoenzymes CD39 and CD73 on the T_reg cells surface, which convert extracellular ATP, a potent pro-inflammatory mediator, into its anti-inflammatory counterpart, adenosine^41,63. Mechanistically, we demonstrated that CXXC1 interacts with FOXP3 to regulate T_reg cell function by trimethylating H3K4 at broad H3K4me3 domains of multiple genes involved in suppressive functions such as Il2ra, Nt5e, and Ctla4. Additionally, MLL1, another KMT, controls T_reg cell activation and function by specifically regulating H3K4 trimethylation at genes encoding key T_reg-related molecules such as Tigit, Klrg1, Tbx21, Cxcr3, and serves as a crucial epigenetic regulator in establishing a stable Th1-T_reg lineage¹⁴.

Our results offer provide novel insights into the suppressive functions, heterogeneity, and regulatory mechanisms of T_reg cells. Maintaining T_reg cell homeostasis and function remains a significant challenge in harnessing T_reg cells for the treatment of autoimmune diseases and the prevention of graft rejection.

Acknowledgements

We thank B. Li (Shanghai Jiao Tong University School of Medicine, Shanghai, China) for providing Foxp3^YFP-Cre mice, and X. L. Liu (Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences) for his generous gifts of cell lines. We thank Y. Huang, Y. Li, J. Wan, C.Guo, and Z. Lin from the Core Facilities, Zhejiang University School of Medicine for their technical support and S. Hong, Y. Ding, H. Jin, Q. Wang, and X. Zhang from Animal Facilities, Zhejiang University, for feeding the mice.

Additional information

Competing interests

The authors declare that they have no competing interests.

Footnotes

This work was supported by grants from the National Natural Science Foundation of China (32341002, 32030035, 32321002, 32100693 and 32270839), the National Key R & D Program of China (2023YFA1800202), the Zhejiang Provincial Natural Science Foundation of China (LZ21C080001, LZ23C070003), Science and Technology Innovation 2030-Major Project (2021ZD0200405), Key project of the Experimental Technology Program of Zhejiang University (SZD202203).

Author contributions

L.W. and L.S.: supervision, conceptualization, project administration, and writing—review and editing. L.W., L.S. and X.M: funding acquisition. X.M: investigation, methodology, project administration, and writing— original draft. Y.Z.and K.L.: bioinformatics analysis. Y. W., X.L., J.C., C.L., Y.Z., S.W., S.T., Q.X., L.D., X.S., and X.G.: investigation and methodology.

References

1
1. Sakaguchi S.
2. Yamaguchi T.
3. Nomura T.
4. Ono M.
2008Regulatory T cells and immune toleranceCell 133:775–787https://doi.org/10.1016/j.cell.2008.05.009
2
1. Josefowicz S. Z.
2. Lu L. F.
3. Rudensky A. Y.
2012Annual Review of ImmunologyAnnual Review of Immunology 30:531–564
3
1. Singer B. D.
2. King L. S.
3. D’Alessio F. R
2014Regulatory T cells as immunotherapyFrontiers in Immunology 5https://doi.org/10.3389/fimmu.2014.00046
4
1. Ohkura N.
2. Sakaguchi S
2020Transcriptional and epigenetic basis of Treg cell development and function: its genetic anomalies or variations in autoimmune diseasesCell Research 30:465–474https://doi.org/10.1038/s41422-020-0324-7
5
1. Littman D. R.
2. Rudensky A. Y
2010Th17 and Regulatory T Cells in Mediating and Restraining InflammationCell 140:845–858https://doi.org/10.1016/j.cell.2010.02.021
6
1. Sakaguchi S.
2. Vignali D. A. A.
3. Rudensky A. Y.
4. Niec R. E.
5. Waldmann H
2013The plasticity and stability of regulatory T cellsNature Reviews Immunology 13:461–467https://doi.org/10.1038/nri3464
7
1. Fontenot J. D.
2. Gavin M. A.
3. Rudensky A. Y
2003Foxp3 programs the development and function of CD4+CD25+ regulatory T cellsNat Immunol 4:330–336https://doi.org/10.1038/ni904
8
1. Van Gool F.
2. Nguyen M. L. T.
3. Mumbach M. R.
4. Satpathy A. T.
5. Rosenthal W. L.
6. Giacometti S.
7. Le D. T.
8. Liu W.
9. Brusko T. M.
10. Anderson M. S.
11. Rudensky A. Y.
12. Marson A.
13. Chang H. Y.
14. Bluestone J. A
2019A Mutation in the Transcription Factor Foxp3 Drives T Helper 2 Effector Function in Regulatory T CellsImmunity 50https://doi.org/10.1016/j.immuni.2018.12.016
9
1. Katoh H.
2. Qin Z. S.
3. Liu R.
4. Wang L.
5. Li W.
6. Li X.
7. Wu L.
8. Du Z.
9. Lyons R.
10. Liu C.-G.
11. Liu X.
12. Dou Y.
13. Zheng P.
14. Liu Y
2011FOXP3 Orchestrates H4K16 Acetylation and H3K4 Trimethylation for Activation of Multiple Genes by Recruiting MOF and Causing Displacement of PLU-1Molecular Cell 44:770–784https://doi.org/10.1016/j.molcel.2011.10.012
10
1. Wang L.
2. Liu Y.
3. Han R.
4. Beier U. H.
5. Bhatti T. R.
6. Akimova T.
7. Greene M. I.
8. Hiebert S. W.
9. Hancock W. W
2015FOXP3(+) regulatory T cell development and function require histone/protein deacetylase 3Journal of Clinical Investigation 125:1111–1123https://doi.org/10.1172/jci77088
11
1. Arvey A.
2. van der Veeken J.
3. Samstein R. M.
4. Feng Y.
5. Stamatoyannopoulos J. A.
6. Rudensky A. Y
2014Inflammation-induced repression of chromatin bound by the transcription factor Foxp3 in regulatory T cellsNature Immunology 15:580–587https://doi.org/10.1038/ni.2868
12
1. Ruthenburg A. J.
2. Allis C. D.
3. Wysocka J
2007Methylation of lysine 4 on histone H3: Intricacy of writing and reading a single epigenetic markMolecular Cell 25:15–30https://doi.org/10.1016/j.molcel.2006.12.014
13
1. Takahashi Y.-h.
2. Westfield G. H.
3. Oleskie A. N.
4. Trievel R. C.
5. Shilatifard A.
6. Skiniotis G.
2011Structural analysis of the core COMPASS family of histone H3K4 methylases from yeast to humanProceedings of the National Academy of Sciences of the United States of America 108:20526–20531https://doi.org/10.1073/pnas.1109360108
14
1. Wang T.
2. Guo J.
3. Li L.
4. Jin Q.
5. Zhang F.
6. Hou B.
7. Zhang Y.
8. Zhou X
2024The histone lysine methyltransferase MLL1 regulates the activation and functional specialization of regulatory T cellsCell Reports 43https://doi.org/10.1016/j.celrep.2024.114222
15
1. Placek K.
2. Hu G.
3. Cui K.
4. Zhang D.
5. Ding Y.
6. Lee J.-E.
7. Jang Y.
8. Wang C.
9. Konkel J. E.
10. Song J.
11. Liu C.
12. Ge K.
13. Chen W.
14. Zhao K
2017MLL4 prepares the enhancer landscape for Foxp3 induction via chromatin loopingNature Immunology 18https://doi.org/10.1038/ni.3812
16
1. Yang Y.
2. Yang Y.
3. Chan K.
4. Couture J.-F
2021Analyzing the impact of CFP1 mutational landscape on epigenetic signalingFaseb Journal 35https://doi.org/10.1096/fj.202100427R
17
1. Tate C. M.
2. Lee J.-H.
3. Skalnik D. G
2010CXXC finger protein 1 restricts the Setd1A histone H3K4 methyltransferase complex to euchromatinFebs Journal 277:210–223https://doi.org/10.1111/j.1742-4658.2009.07475.x
18
1. Cao W.
2. Guo J.
3. Wen X.
4. Miao L.
5. Lin F.
6. Xu G.
7. Ma R.
8. Yin S.
9. Hui Z.
10. Chen T.
11. Guo S.
12. Chen W.
13. Huang Y.
14. Liu Y.
15. Wang J.
16. Wei L.
17. Wang L
2016CXXC finger protein 1 is critical for T-cell intrathymic development through regulating H3K4 trimethylationNat Commun 7https://doi.org/10.1038/ncomms11687
19
1. Hui Z.
2. Zhou L.
3. Xue Z.
4. Zhou L.
5. Luo Y.
6. Lin F.
7. Liu X.
8. Hong S.
9. Li W.
10. Wang D.
11. Lu L.
12. Wang J.
13. Wang L
2018Cxxc Finger Protein 1 Positively Regulates GM-CSF-Derived Macrophage Phagocytosis Through Csf2rα-Mediated SignalingFrontiers in Immunology 9https://doi.org/10.3389/fimmu.2018.01885
20
1. Lin F.
2. Meng X.
3. Guo Y.
4. Cao W.
5. Liu W.
6. Xia Q.
7. Hui Z.
8. Chen J.
9. Hong S.
10. Zhang X.
11. Wu C.
12. Wang D.
13. Wang J.
14. Lu L.
15. Qian W.
16. Wei L.
17. Wang L
2019Epigenetic initiation of the T(H)17 differentiation program is promoted by Cxxc finger protein 1Science Advances 5https://doi.org/10.1126/sciadv.aax1608
21
1. Shen X.
2. Gao X.
3. Luo Y.
4. Xu Q.
5. Fan Y.
6. Hong S.
7. Huang Z.
8. Liu X.
9. Wang Q.
10. Chen Z.
11. Wang D.
12. Lu L.
13. Wu C.
14. Liang H.
15. Wang L
2023Cxxc finger protein 1 maintains homeostasis and function of intestinal group 3 innate lymphoid cells with agingNature Aging 3https://doi.org/10.1038/s43587-023-00453-7
22
1. Cao W.
2. Guo J.
3. Wen X.
4. Miao L.
5. Lin F.
6. Xu G.
7. Ma R.
8. Yin S.
9. Hui Z.
10. Chen T.
11. Guo S.
12. Chen W.
13. Huang Y.
14. Liu Y.
15. Wang J.
16. Wei L.
17. Wang L
2016CXXC finger protein 1 is critical for T-cell intrathymic development through regulating H3K4 trimethylationNature Communications 7https://doi.org/10.1038/ncomms11687
23
1. Zhou L.
2. Wang S.
3. Hu W.
4. Liu X.
5. Xu L.
6. Tong B.
7. Zhang T.
8. Xue Z.
9. Guo Y.
10. Zhao J.
11. Lu L.
12. Fan H.
13. Qian W.
14. Chen J.
15. Chen W.
16. Wang L
2023T cell proliferation requires ribosomal maturation in nucleolar condensates dependent on DCAF13Journal of Cell Biology 222https://doi.org/10.1083/jcb.202201096
24
1. Chou W.-C.
2. Guo Z.
3. Guo H.
4. Chen L.
5. Zhang G.
6. Liang K.
7. Xie L.
8. Tan X.
9. Gibson S. A.
10. Rampanelli E.
11. Wang Y.
12. Montgomery S. A.
13. Brickey W. J.
14. Deng M.
15. Freeman L.
16. Zhang S.
17. Su M. A.
18. Chen X.
19. Wan Y. Y.
20. Ting J. P. Y
2021AIM2 in regulatory T cells restrains autoimmune diseasesNature 591https://doi.org/10.1038/s41586-021-03231-w
25
1. Zeng H.
2. Yang K.
3. Cloer C.
4. Neale G.
5. Vogel P.
6. Chi H
2013mTORC1 couples immune signals and metabolic programming to establish T(reg)-cell functionNature 499:485–490https://doi.org/10.1038/nature12297
26
1. Lu D.
2. Liu L.
3. Sun Y.
4. Song J.
5. Yin Q.
6. Zhang G.
7. Qi F.
8. Hu Z.
9. Yang Z.
10. Zhou Z.
11. Hu Y.
12. Zhang L.
13. Ji J.
14. Zhao X.
15. Jin Y.
16. McNutt M. A.
17. Yin Y
2020The phosphatase PAC1 acts as a T cell suppressor and attenuates host antitumor immunityNature Immunology 21https://doi.org/10.1038/s41590-019-0577-9
27
1. Chen J.
2. Liu K.
3. Luo Y.
4. Kang M.
5. Wang J.
6. Chen G.
7. Qi J.
8. Wu W.
9. Wang B.
10. Han Y.
11. Shi L.
12. Wang K.
13. Han X.
14. Ma X.
15. Liu W.
16. Ding Y.
17. Wang L.
18. Liang H.
19. Wang L.
20. Chen J
2023Single-Cell Profiling of Tumor Immune Microenvironment Reveals Immune Irresponsiveness in Gastric Signet-Ring Cell CarcinomaGastroenterology 165:88–103https://doi.org/10.1053/j.gastro.2023.03.008
28
1. Marson A.
2. Kretschmer K.
3. Frampton G. M.
4. Jacobsen E. S.
5. Polansky J. K.
6. Maclsaac K. D.
7. Levine S. S.
8. Fraenkel E.
9. von Boehmer H.
10. Young R. A
2007Foxp3 occupancy and regulation of key target genes during T-cell stimulationNature 445:931–935https://doi.org/10.1038/nature05478
29
1. Zheng Y.
2. Josefowicz S. Z.
3. Kas A.
4. Chu T.-T.
5. Gavin M. A.
6. Rudensky A. Y
2007Genome-wide analysis of Foxp3 target genes in developing and mature regulatory T cellsNature 445:936–940https://doi.org/10.1038/nature05563
30
1. Ono M.
2. Yaguchi H.
3. Ohkura N.
4. Kitabayashi I.
5. Nagamura Y.
6. Nomura T.
7. Miyachi Y.
8. Tsukada T.
9. Sakaguchi S
2007Foxp3 controls regulatory T-cell function by interacting with AML1/Runx1Nature 446:685–689https://doi.org/10.1038/nature05673
31
1. Fu W.
2. Ergun A.
3. Lu T.
4. Hill J. A.
5. Haxhinasto S.
6. Fassett M. S.
7. Gazit R.
8. Adoro S.
9. Glimcher L.
10. Chan S.
11. Kastner P.
12. Rossi D.
13. Collins J. J.
14. Mathis D.
15. Benoist C
2012A multiply redundant genetic switch ‘locks in’ the transcriptional signature of regulatory T cellsNature Immunology 13:972–980https://doi.org/10.1038/ni.2420
32
1. Ohkura N.
2. Hamaguchi M.
3. Morikawa H.
4. Sugimura K.
5. Tanaka A.
6. Ito Y.
7. Osaki M.
8. Tanaka Y.
9. Yamashita R.
10. Nakano N.
11. Huehn J.
12. Fehling H. J.
13. Sparwasser T.
14. Nakai K.
15. Sakaguchi S
2012T Cell Receptor Stimulation-Induced Epigenetic Changes and Foxp3 Expression Are Independent and Complementary Events Required for Treg Cell DevelopmentImmunity 37:785–799https://doi.org/10.1016/j.immuni.2012.09.010
33
1. Wei G.
2. Wei L.
3. Zhu J.
4. Zang C.
5. Hu-Li J.
6. Yao Z.
7. Cui K.
8. Kanno Y.
9. Roh T.-Y.
10. Watford W. T.
11. Schones D. E.
12. Peng W.
13. Sun H.-w.
14. Paul W. E.
15. O’Shea J. J.
16. Zhao K.
2009Global Mapping of H3K4me3 and H3K27me3 Reveals Specificity and Plasticity in Lineage Fate Determination of Differentiating CD4(+) T CellsImmunity 30:155–167https://doi.org/10.1016/j.immuni.2008.12.009
34
1. Konopacki C.
2. Pritykin Y.
3. Rubtsov Y.
4. Leslie C. S.
5. Rudensky A. Y
2019Transcription factor Foxp1 regulates Foxp3 chromatin binding and coordinates regulatory T cell functionNature Immunology 20https://doi.org/10.1038/s41590-018-0291-z
35
1. Oh H.
2. Grinberg-Bleyer Y.
3. Liao W.
4. Maloney D.
5. Wang P.
6. Wu Z.
7. Wang J.
8. Bhatt D. M.
9. Heise N.
10. Schmid R. M.
11. Hayden M. S.
12. Klein U.
13. Rabadan R.
14. Ghosh S
2017An NF-kappaB Transcription-Factor-Dependent Lineage-Specific Transcriptional Program Promotes Regulatory T Cell Identity and FunctionImmunity 47:450–465https://doi.org/10.1016/j.immuni.2017.08.010
36
1. Rubtsov Y. P.
2. Rasmussen J. P.
3. Chi E. Y.
4. Fontenot J.
5. Castelli L.
6. Ye X.
7. Treuting P.
8. Siewe L.
9. Roers A.
10. Henderson W. R.
11. Muller W.
12. Rudensky A. Y
2008Regulatory T cell-derived interleukin-10 limits inflammation at environmental interfacesImmunity 28:546–558https://doi.org/10.1016/j.immuni.2008.02.017
37
1. Fontenot J. D.
2. Gavin M. A.
3. Rudensky A. Y
2003Foxp3 programs the development and function of CD4(+)CD25(+) regulatory T cellsNature Immunology 4:330–336https://doi.org/10.1038/ni904
38
1. Kim J. M.
2. Rasmussen J. P.
3. Rudensky A. Y
2007Regulatory T cells prevent catastrophic autoimmunity throughout the lifespan of miceNature Immunology 8:191–197https://doi.org/10.1038/ni1428
39
1. Shevyrev D.
2. Tereshchenko V.
2020Treg Heterogeneity, Function, and HomeostasisFrontiers in Immunology 10https://doi.org/10.3389/fimmu.2019.03100
40
1. Zagorulya M.
2. Yim L.
3. Morgan D. M.
4. Edwards A.
5. Torres-Mejia E.
6. Momin N.
7. McCreery C. V.
8. Zamora I. L.
9. Horton B. L.
10. Fox J. G.
11. Wittrup K. D.
12. Love J. C.
13. Spranger S
2023Tissue-specific abundance of interferon-gamma drives regulatory T cells to restrain DC1-mediated priming of cytotoxic T cells against lung cancerImmunity 56https://doi.org/10.1016/j.immuni.2023.01.010
41
1. Dikiy S.
2. Rudensky A. Y
2023Principles of regulatory T cell functionImmunity 56:240–255https://doi.org/10.1016/j.immuni.2023.01.004
42
1. Street K.
2. Risso D.
3. Fletcher R. B.
4. Das D.
5. Ngai J.
6. Yosef N.
7. Purdom E.
8. Dudoit S
2018Slingshot: cell lineage and pseudotime inference for single-cell transcriptomicsBmc Genomics 19https://doi.org/10.1186/s12864-018-4772-0
43
1. Tarbell K. V.
2. Yamazaki S.
3. Olson K.
4. Toy P.
5. Steinman R. M
2004CD25+ CD4+ T cells, expanded with dendritic cells presenting a single autoantigenic peptide, suppress autoimmune diabetesJournal of Experimental Medicine 199:1467–1477https://doi.org/10.1084/jem.20040180
44
1. Hori S.
2. Haury M.
3. Coutinho A.
4. Demengeot J
2002Specificity requirements for selection and effector functions of CD25+4+ regulatory T cells in anti-myelin basic protein T cell receptor transgenic miceProceedings of the National Academy of Sciences of the United States of America 99:8213–8218https://doi.org/10.1073/pnas.122224799
45
1. Butler J. S.
2. Lee J.-H.
3. Skalnik D. G
2008CFP1 Interacts with DNMT1 Independently of Association with the Setd1 Histone H3K4 Methyltransferase ComplexesDNA and Cell Biology 27:533–543https://doi.org/10.1089/dna.2007.0714
46
1. Butler J. S.
2. Palam L. R.
3. Tate C. M.
4. Sanford J. R.
5. Wek R. C.
6. Skalnik D. G
2009DNA Methyltransferase Protein Synthesis Is Reduced in CXXC Finger Protein 1-Deficient Embryonic Stem CellsDNA and Cell Biology 28:223–231https://doi.org/10.1089/dna.2009.0854
47
1. Benayoun B. A.
2. Pollina E. A.
3. Ucar D.
4. Mahmoudi S.
5. Karra K.
6. Wong E. D.
7. Devarajan K.
8. Daugherty A. C.
9. Kundaje A. B.
10. Mancini E.
11. Hitz B. C.
12. Gupta R.
13. Rando T. A.
14. Baker J. C.
15. Snyder M. P.
16. Cherry J. M.
17. Brunet A
2014H3K4me3 Breadth Is Linked to Cell Identity and Transcriptional ConsistencyCell 158:673–688https://doi.org/10.1016/j.cell.2014.06.027
48
1. Zacarias-Cabeza J.
2. Belhocine M.
3. Vanhille L.
4. Cauchy P.
5. Koch F.
6. Pekowska A.
7. Fenouil R.
8. Bergon A.
9. Gut M.
10. Gut I.
11. Eick D.
12. Imbert J.
13. Ferrier P.
14. Andrau J.-C.
15. Spicuglia S
2015Transcription-Dependent Generation of a Specialized Chromatin Structure at the TCRβ LocusJournal of Immunology 194:3432–3443https://doi.org/10.4049/jimmunol.1400789
49
1. Chen K.
2. Chen Z.
3. Wu D.
4. Zhang L.
5. Lin X.
6. Su J.
7. Rodriguez B.
8. Xi Y.
9. Xia Z.
10. Chen X.
11. Shi X.
12. Wang Q.
13. Li W
2015Broad H3K4me3 is associated with increased transcription elongation and enhancer activity at tumor-suppressor genesNature Genetics 47https://doi.org/10.1038/ng.3385
50
1. Zhang R.
2. Huynh A.
3. Whitcher G.
4. Chang J.
5. Maltzman J. S.
6. Turka L. A
2013An obligate cell-intrinsic function for CD28 in TregsJournal of Clinical Investigation 123:580–593https://doi.org/10.1172/jci65013
51
1. Gokhale A. S.
2. Gangaplara A.
3. Lopez-Occasio M.
4. Thornton A. M.
5. Shevach E. M
2019Selective deletion of Eos (Ikzf4) in T-regulatory cells leads to loss of suppressive function and development of systemic autoimmunityJournal of Autoimmunity 105https://doi.org/10.1016/j.jaut.2019.06.011
52
1. Wan Y. Y.
2. Flavell R. A
2007Regulatory T-cell functions are subverted and converted owing to attenuated Foxp3 expressionNature 445:766–770https://doi.org/10.1038/nature05479
53
1. Chen C.
2. Rowell E. A.
3. Thomas R. M.
4. Hancock W. W.
5. Wells A. D
2006Transcriptional regulation by Foxp3 is associated with direct promoter occupancy and modulation of histone acetylationJournal of Biological Chemistry 281:36828–36834https://doi.org/10.1074/jbc.M608848200
54
1. Morikawa H.
2. Ohkura N.
3. Vandenbon A.
4. Itoh M.
5. Nagao-Sato S.
6. Kawaji H.
7. Lassmann T.
8. Carninci P.
9. Hayashizaki Y.
10. Forrest A. R. R.
11. Standley D. M.
12. Date H.
13. Sakaguchi S.
14. Consortium F
2014Differential roles of epigenetic changes and Foxp3 expression in regulatory T cell-specific transcriptional regulationProceedings of the National Academy of Sciences of the United States of America 111:5289–5294https://doi.org/10.1073/pnas.1312717110
55
1. Pan F.
2. Yu H.
3. Dang E. V.
4. Barbi J.
5. Pan X.
6. Grosso J. F.
7. Jinasena D.
8. Sharma S. M.
9. McCadden E. M.
10. Getnet D.
11. Drake C. G.
12. Liu J. O.
13. Ostrowski M. C.
14. Pardoll D. M
2009Eos Mediates Foxp3-Dependent Gene Silencing in CD4+ Regulatory T CellsScience 325:1142–1146https://doi.org/10.1126/science.1176077
56
1. DuPage M.
2. Chopra G.
3. Quiros J.
4. Rosenthal W. L.
5. Morar M. M.
6. Holohan D.
7. Zhang R.
8. Turka L.
9. Marson A.
10. Bluestone J. A
2015The Chromatin-Modifying Enzyme Ezh2 Is Critical for the Maintenance of Regulatory T Cell Identity after ActivationImmunity 42:227–238https://doi.org/10.1016/j.immuni.2015.01.007
57
1. Liu X.
2. Wang C.
3. Liu W.
4. Li J.
5. Li C.
6. Kou X.
7. Chen J.
8. Zhao Y.
9. Gao H.
10. Wang H.
11. Zhang Y.
12. Gao Y.
13. Gao S
2016Distinct features of H3K4me3 and H3K27me3 chromatin domains in pre-implantation embryosNature 537https://doi.org/10.1038/nature19362
58
1. Sze C. C.
2. Ozark P. A.
3. Cao K.
4. Ugarenko M.
5. Das S.
6. Wang L.
7. Marshall S. A.
8. Rendleman E. J.
9. Ryan C. A.
10. Zha D.
11. Douillet D.
12. Chen F. X.
13. Shilatifard A
2020Coordinated regulation of cellular identity-associated H3K4me3 breadth by the COMPASS familyScience Advances 6https://doi.org/10.1126/sciadv.aaz4764
59
1. Dahl J. A.
2. Jung I.
3. Aanes H.
4. Greggains G. D.
5. Manaf A.
6. Lerdrup M.
7. Li G.
8. Kuan S.
9. Li B.
10. Lee A. Y.
11. Preissl S.
12. Jermstad I.
13. Haugen M. H.
14. Suganthan R.
15. Bjoras M.
16. Hansen K.
17. Dalen K. T.
18. Fedorcsak P.
19. Ren B.
20. Klungland A
2016Broad histone H3K4me3 domains in mouse oocytes modulate maternal-to-zygotic transitionNature 537https://doi.org/10.1038/nature19360
60
1. Zhang B.
2. Zheng H.
3. Huang B.
4. Li W.
5. Xiang Y.
6. Peng X.
7. Ming J.
8. Wu X.
9. Zhang Y.
10. Xu Q.
11. Liu W.
12. Kou X.
13. Zhao Y.
14. He W.
15. Li C.
16. Chen B.
17. Li Y.
18. Wang Q.
19. Ma J.
20. Yin Q.
21. Kee K.
22. Meng A.
23. Gao S.
24. Xu F.
25. Na J.
26. Xie W
2016Allelic reprogramming of the histone modification H3K4me3 in early mammalian developmentNature 537https://doi.org/10.1038/nature19361
61
1. Zhang Z.
2. Shi L.
3. Dawany N.
4. Kelsen J.
5. Petri M. A.
6. Sullivan K. E
2016H3K4 tri-methylation breadth at transcription start sites impacts the transcriptome of systemic lupus erythematosusClinical Epigenetics 8https://doi.org/10.1186/s13148-016-0179-4
62
1. Wong S. H. K.
2. Goode D. L.
3. Iwasaki M.
4. Wei M. C.
5. Kuo H.-P.
6. Zhu L.
7. Schneidawind D.
8. Duque-Afonso J.
9. Weng Z.
10. Cleary M. L
2015The H3K4-Methyl Epigenome Regulates Leukemia Stem Cell Oncogenic PotentialCancer Cell 28:198–209https://doi.org/10.1016/j.ccell.2015.06.003
63
1. Deaglio S.
2. Dwyer K. M.
3. Gao W.
4. Friedman D.
5. Usheva A.
6. Erat A.
7. Chen J.-F.
8. Enjyoji K.
9. Linden J.
10. Oukka M.
11. Kuchroo V. K.
12. Strom T. B.
13. Robson S. C
2007Adenosine generation catalyzed by CD39 and CD73 expressed on regulatory T cells mediates immune suppressionJournal of Experimental Medicine 204:1257–1265https://doi.org/10.1084/jem.20062512

Article and author information

Author information

Xiaoyu Meng
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
- These authors contributed equally to this work.
Yezhang Zhu
Department of Hematology, Tongji Hospital, School of Medicine, Tongji University, Shanghai, China, Shanghai Key Laboratory of Signaling and Disease Research, Frontier Science Center of Stem Cell Research, National Stem Cell Translational Resource Center, School of Life Sciences and Technology, Tongji University, Shanghai, China
- These authors contributed equally to this work.
Kuai Liu
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
- These authors contributed equally to this work.
Yuxi Wang
Laboratory Animal Center, Zhejiang University, Hangzhou, China
Xiaoqian Liu
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
Chenxin Liu
Zhejiang University School of Medicine, Hangzhou, China
Yan Zeng
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
Shuai Wang
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
Xianzhi Gao
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China
Xin Shen
Co-Facility Center, Zhejiang University School of Medicine, Hangzhou, China
Jing Chen
Department of Gastrointestinal Surgery, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China
Sijue Tao
Laboratory Animal Center, Zhejiang University, Hangzhou, China
Qianying Xu
Zhejiang University School of Medicine, Hangzhou, China
Linjia Dong
School of Basic Medical Sciences and Forensic Medicine, Hangzhou Medical College, Hangzhou, China
Li Shen
MOE Key Laboratory of Biosystems Homeostasis & Protection and Zhejiang Provincial Key Laboratory for Cancer Molecular Cell Biology, Life Sciences Institute, Zhejiang University, Hangzhou, China, Department of Orthopedics Surgery, the Second Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou, China
ORCID iD: 0000-0002-5696-2191
- For correspondence: li_shen@zju.edu.cn
Lie Wang
Institute of Immunology and Bone Marrow Transplantation Center, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou, China, Zhejiang University School of Medicine, Hangzhou, China, Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China, Laboratory Animal Center, Zhejiang University, Hangzhou, China, Future Health Laboratory, Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing, China
ORCID iD: 0000-0001-5094-012X
- For correspondence: wanglie@zju.edu.cn

Version history

Sent for peer review: September 25, 2024
Preprint posted: October 4, 2024
Reviewed Preprint version 1: November 21, 2024

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

views: 115
downloads: 10
citations: 0

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Significance of findings

Strength of evidence

Abstract

Introduction

Materials and methods

Mice

Cell culture

Immunofluorescence (IF) Staining

Co-immunoprecipitation and western blot

Flow cytometry

Real-time PCR

CD4+T and YFP+ Treg cells adoptive transfer in Experimental autoimmune encephalomyelitis (EAE)

Isolation lymphocytes from the central nervous system (CNS)

Histological analyses

Adoptive transfer colitis model

In vitro Treg suppression assay

CUT&Tag

CUT&Tag and ChIP-seq data analysis

Clustering analysis

RNA-seq and data analysis

Single-Cell RNA-Sequencing

Targeted scRNA-seq data processing

Analysis of the single-cell TCR-seq repertoire

scRNA-seq trajectory analysis

Whole genome bisulfite sequencing (WGBS) and data analysis

Statistical analysis

Results

FOXP3 binds regulatory loci primed for activation and repression in Treg cells

H3K4me3 is required for FOXP3-dependent gene activation in Treg cells.

CXXC1 interacts with FOXP3 and binds H3K4me3-enriched sites in Treg cells

CXXC1 interacts with FOXP3 in Treg cell.

Complete ablation of Cxxc1 in Treg cells leads to a fatal autoimmune disease

Foxp3creCxxc1fl/fl mice spontaneously develop a fatal early-onset inflammatory disorder.

CXXC1 is necessary for the maintenance of Treg-cell suppressive activity

CXXC1 is essential for Treg cells to suppress T cell-mediated EAE and colitis.

Treg cell lineage stability and proliferation depend upon CXXC1

Single-cell transcriptomics reveals distinct Treg cell populations.

Cxxc1 deletion leads to a reduction in essential Treg cell markers

Cxxc1-deficient Treg cells in het-KO female mice exhibit reduced suppression markers under steady-state conditions.

The FOXP3-CXXC1 complex regulates the expression of key factors in Treg cells that are associated with the breadth of H3K4me3

CXXC1 regulates Foxp3-dependent molecule H3K4 trimethylation in Treg cells.

Discussion

Acknowledgements

Additional information

Competing interests

Footnotes

Author contributions

References

Article and author information

Author information

Xiaoyu Meng†

Yezhang Zhu†

Kuai Liu†

Yuxi Wang

Xiaoqian Liu

Chenxin Liu

Yan Zeng

Shuai Wang

Xianzhi Gao

Xin Shen

Jing Chen

Sijue Tao

Qianying Xu

Linjia Dong

Li Shen

Lie Wang

Version history

Copyright

Metrics

CD4⁺T and YFP⁺ T_reg cells adoptive transfer in Experimental autoimmune encephalomyelitis (EAE)

In vitro T_reg suppression assay

FOXP3 binds regulatory loci primed for activation and repression in T_reg cells

H3K4me3 is required for FOXP3-dependent gene activation in T_reg cells.

CXXC1 interacts with FOXP3 and binds H3K4me3-enriched sites in T_reg cells

CXXC1 interacts with FOXP3 in T_reg cell.

Complete ablation of Cxxc1 in T_reg cells leads to a fatal autoimmune disease

Foxp3^creCxxc1^fl/fl mice spontaneously develop a fatal early-onset inflammatory disorder.

CXXC1 is necessary for the maintenance of T_reg-cell suppressive activity

CXXC1 is essential for T_reg cells to suppress T cell-mediated EAE and colitis.

T_reg cell lineage stability and proliferation depend upon CXXC1

Single-cell transcriptomics reveals distinct T_reg cell populations.

Cxxc1 deletion leads to a reduction in essential T_reg cell markers

Cxxc1-deficient T_reg cells in het-KO female mice exhibit reduced suppression markers under steady-state conditions.

The FOXP3-CXXC1 complex regulates the expression of key factors in T_reg cells that are associated with the breadth of H3K4me3

CXXC1 regulates Foxp3-dependent molecule H3K4 trimethylation in T_reg cells.

Xiaoyu Meng

Yezhang Zhu

Kuai Liu