Impact of loss of Ikaros function on peripheral Treg phenotype.

Example histograms (a) and quantified expression (b) of Ikaros protein by peripheral Treg and Tconv from Foxp3-YFP-Cre (green) and Ikzf1-fl-Foxp3-YFP-Cre (red) mice (N=6 mice per group). (c) Expression of CD25, ICOS, GITR, and PD1 by WT (green) and Ikzf1-cko (red) Treg (N=6 mice per group). (d) Flow cytometric measurement of Foxp3, Aiolos, Eos, and Helios protein expression in WT and Ikzf1-cko Treg (N=6 mice per group). (e) Transcriptomic analysis of WT vs. Ikzf1-cko Treg gene expression. Top differentially expressed genes (FDR<0.05) organized into 8 clusters in ex vivo or in vitro stimulated WT and Ikzf1-cko Treg. The heatmap represents scaled transcripts per million (tpm, N=3 replicates per group).

Survey of relevant differentially expressed genes in WT vs. Ikzf1-cko Treg.

(a) Known negative regulators of Treg function up-regulated in ex vivo Ikzf1-cko Treg, (b) positive regulators of Treg function down-regulated in ex vivo Ikzf1-cko Treg, (c) negative regulators of Treg function up-regulated in in vitro stimulated Ikzf1-cko Treg, (d) positive regulators of Treg function down-regulated in in vitro stimulated Ikzf1-cko Treg. See Supplementary Table 1 for all differential genes and Supplementary Table 2 for a larger list of functionally relevant differentially expressed genes.

In vitro analysis of WT vs. Ikzf1-cko Treg.

(a) Secretion of IL-2 and IFNg protein by WT and Ikzf1-cko Treg measured by ELISA (N=3). (b) IL-2, IFNg, and TNFa production by Treg from WT mice vs. mice with a dominant-negative form of Ikaros (IkDN) measured by ELISA (N=3). (c) IL-2-induced phosphorylation of STAT5 in WT and Ikzf1-cko Treg measured by flow cytometry in vitro. Mean fluorescence intensity (MFI) and individual histograms (inset, N=3) are depicted. (d) Activation-induced proliferation of WT (closed) and Ikzf1-cko (open) Treg measured by dye dilution (N=3). (e) Gene ontology analysis of genes down-regulated (top panel) and up-regulated (bottom panel) in Ikzf1-cko compared to WT Treg. The x-axis is fold enrichment and numbers to the right are unique genes in each pathway. (f) Differential expression of core Treg, nTreg, aTreg, cTreg, and eTreg genes19 in ex vivo (left panel) and in vitro stimulated (right panel) WT and Ikzf1-cko Treg. The heatmap represents z-score and genes significantly differentially expressed are shown in red.

Ectopic activation of the Wnt-catenin pathway in Ikzf1-cko Treg.

(a) Expression of B-catenin by WT (left histogram) vs. Ikzf1-cko (right histogram) Treg measured by flow cytometry (plot depicts B-catenin MFI from N=3 experiments). (b) IFNg secretion by WT (top panels) vs. Ikzf1-cko (bottom panels) Treg activated with (right panels) or without (left panels) the Wnt pathway inhibitor PKF. (c) IFNg secretion (measured by ELISA) by Ikzf1-cko Treg, but not by conventional T cells, is inhibited in a dose-dependent manner by PKF (N=3).

Genome-wide assessment of open chromatin, histone acetylation, and Ikaros occupancy in WT and Ikzf1-cko Treg.

Differential analysis of open chromatin (a) and H3K27ac (b) in WT vs. Ikzf1-cko Treg (FDR<0.05, N=3). Peaks with FC>2 are depicted in red. (c) Ikaros ChIP-seq signal (read density) at genomic regions shared (panel 2) in Tconv (blue) vs. Treg (red), unique to Tconv (panel 3), or unique to Treg (panel 4). Panel 1 depicts read densities at the same regions in control input libraries. (d) Unique vs. shared Ikaros binding sites in Treg (red) vs. Tconv (green) open chromatin. (e) Enrichment of Ikaros ChIP-seq signal (footprint, solid lines) or input background (transparent lines) at accessible Ikaros motifs (AGGAA) in WT Treg (blue) and Tconv (red). (f) Expression (tpm) of genes with open chromatin enriched for the Ikaros consensus binding motif in WT vs. Ikzf1-cko Treg. (g) Enrichment of TF consensus binding motifs in genomic regions that are more accessible in Ikzf1-cko Treg. Inset depicts motifs for Ikaros/Ikzf1, Stat, Nfat, Rbpj, Rela, and AP-1. (h) Open chromatin (top tracks) and H3K27ac (bottom tracks) in WT (blue) and Ikzf1-cko (red) Treg, Ikaros binding sites (purple marks) and Foxp3 binding sites (red marks) in WT Treg, and Foxp3 binding sites in Ikzf1-cko Treg (green marks) at the Bcl6 locus. (i) Enrichment of TF consensus binding motifs in genomic regions that are less accessible in Ikzf1-cko Treg. Inset depicts motifs for forkhead family members. In (g) and (i), factors with roles in Treg function are colored green and purple, forkhead family members are colored red, and factors differentially expressed in Ikzf1-cko Treg are indicated with an asterisk. ATAC-seq was performed on Treg purified directly ex vivo, while ChIP-seq analyses were performed on Treg expanded in vivo using IL-2/anti-IL-2 complexes.

Ikaros-dependent Foxp3 function in Treg.

(a) Input vs. Foxp3 ChIP-seq at genomic regions shared or unique in WT (blue) vs. Ikzf1-cko (red) Treg (N=3 per group). (b) Foxp3- (green), Ikaros- (red) and Foxp3-Ikaros co-bound (orange) OCR. Inset depicts Foxp3-Ikaros co-bound regions in WT (green) vs. Ikzf1-cko (red) Treg. (c) Foxp3 binding sites in WT (green) vs. Ikzf1-cko (red) Treg. (d) Expression (tpm) of genes enriched for accessible Foxp3 consensus motifs in WT vs. Ikzf1-cko Treg. (e) 293T cells transfected with FLAG-tagged full-length Ikaros (Ik1), DNA-binding mutant Ikaros (Ik7), or Runx1 alone (lanes 1-4) or in combination with untagged Foxp3 (lanes 5-7). Whole extracts (panels 1-2) or Foxp3-immunoprecipitated extracts (panels 3-4) probed for Foxp3 or FLAG. (f) IL-2 or IFNg production (left panel) and Foxp3 ChIP-qPCR at Il2 promoter (right panel) in Tconv transduced with vector, Foxp3, Ik7/DN, or Foxp3+Ik7/DN. (g) Concordant vs. discordant genes co-regulated in Ikzf1-cko vs. IkDN Treg (odds ratio=5.19, P value=2.2×10-16). (h,i) Open chromatin (top) and H3K27ac (bottom) in WT (blue) and Ikzf1-cko (red) Treg, Ikaros binding (purple marks) and Foxp3 binding (red marks) in WT Treg, and Foxp3 binding in Ikzf1-cko Treg (green marks) at Il2 (h) and Tcf7 (i).

Immunophenotyping of Tconv from Ikzf1-fl-Foxp3-YFP-Cre and Foxp3-YFP-Cre mice.

Frequencies of total (a), memory (b,c), and naïve (d) phenotype Tconv in secondary lymphoid tissues of 6-8 week-old Ikzf1-fl-Foxp3-YFP-Cre (purple) and Foxp3-YFP-Cre (green) mice (N=7). (e) In vitro suppressive activity of WT (blue) vs. Ikzf1-deficient (orange) Treg against Tconv from WT mice. (f) In vitro suppressive activity of WT Treg against Tconv from WT mice (blue) vs. Ikzf1-deficient Treg against Tconv from Ikzf1-cko mice (orange). (g) In vitro suppressive activity of WT Treg against Tconv from WT (blue) vs. Ikzf1-deficient mice (orange). Tconv proliferation was measured by dye dilution in all cultures (N=3). (h) Frequencies of Foxp3+PD1hiCXCR5hi follicular regulatory T cells (Tfr) and PD1hiCXCR5hi follicular helper T cells (Tfh) in secondary lymphoid tissues of 6 month-old Ikzf1-fl-Foxp3-YFP-Cre (purple) and Foxp3-YFP-Cre (green) mice (N=7). (i) Total serum levels of IgM, IgG, and IgA from Ikzf1-fl-Foxp3-YFP-Cre vs. Foxp3-YFP-Cre mice (N=4). Frequency of IgA-positive B cells and surface density (MFI) of IgA-positive B cells in spleen (j) and mesenteric lymph nodes (k) from 9 month-old WT (green) vs. Ikzf1-deficient (red) mice (N=3). P values are indicated for significant differences.

Role of Ikaros in Treg-mediated control of inflammatory colitis.

WT CD4+CD25-Tconv were transferred alone (red), or together with WT (green) or Ikzf1-cko (blue) CD4+CD25+ Treg into RAG1ko mice (N=5). Animal weight was monitored for 40 days (a), and intestines were scored for pathology at the gross and histologic levels (b). Example histopathology of colons from RAG1ko recipients of WT Tconv (c), WT Tconv + WT Treg (d), and WT Tconv + Ikzf1-cko Treg (e). H&E (top row), CD4 (middle row), and Foxp3 (bottom row) staining are shown at 200x. Mean Foxp3+ cells per 200X field from N=3 animals is 3.4 in (d) and 22 in (e), P<0.05.

Role of Ikaros in Treg-dependent acquired cardiac transplant tolerance.

(a) B6 Ikzf1-fl-Foxp3-YFP-Cre (red) or Foxp3-YFP-Cre (black) mice (N=5) received BALB/c cardiac allografts under combined CD28+CD40 costimulatory blockade and graft survival was monitored for 100 days. (b) Analysis of intra-graft transcript levels of the indicated genes from grafts harvested at day 19 post-transplant. Histopathological analysis of cardiac grafts harvested at day 19 post-transplant from Foxp3-YFP-Cre and Ikzf1-fl-Foxp3-YFP-Cre recipients by H&E (c) and CD4 and Foxp3 (d) staining (scale=200x).

Comparison of Ikaros vs. Helios binding sites and regulated genes.

Venn diagram in (a) depicts comparison of published ChIP-seq Helios binding sites12 (one replicate) with at least one bp of overlap with Ikaros binding sites from this study. Venn diagram in (b) depicts overlap of differentially expressed genes (FDR < 0.5 & abs(log2FC) > 1) in Ikaros or Helios conditional knockout37.

Immunophenotyping of thymi from Ikzf1-fl-Foxp3-YFP-Cre or Foxp3-YFP-Cre control mice (a). Bisulfite DNA methylation analysis of the Foxp3 CNS2/TSDR region in peripheral Treg and Tconv from Ikzf1-fl-Foxp3-YFP-Cre or Foxp3-YFP-Cre control mice (b). Quantification of total peripheral Treg, CD44lo cTreg, CD44hi eTreg, and PD1hiCXCR5hi Tfreg in Ikzf1-fl-Foxp3-YFP-Cre and Foxp3-YFP-Cre mice (c). Bisulfite DNA methylation analysis of the IFNg promoter and IFNg intronic enhancer (d) regions in peripheral Treg and Tconv from Ikzf1-fl-Foxp3-YFP-Cre or Foxp3-YFP-Cre control mice.

ATAC-seq and H3K27ac ChIP-seq library metrics. Assessment of WT and Ikzf1-cko ex vivo and stimulated Treg library reproducibility by principal component (a) and spearman correlation (b) analysis. (c) Hierarchical clustering of 12,030 differentially accessible regions (DAR, FDR<0.05, abs[logFC]>1) into 7 groups from the root. Red vs. blue indicate increased vs. decreased accessibility. (d) Enrichment of GO terms for genes nearest to DAR (FDR<0.05, >5 genes). (e) Association of DAR accessibility changes with changes in expression of the nearest gene. Enrichment is over the mean increase or decrease in expression, indicated by red vs. blue. Assessment of WT and Ikzf1-cko H3K27ac ChIP-seq library reproducibility by principal component (f) and pairwise spearman correlation (g) analysis. Overlap (h) and Pearson correlation between ATAC-seq and H3K27ac ChIP-seq peaks (i) and and H3K27ac and gene expression (j). Identification of super-enhancer regions based on H3K27ac density (k). The Venn diagram (inset) depicts the number of super-enhancers shared between WT and Ikzf1-cko Treg (orange), unique to WT Treg (green), and unique to Ikzf1-cko Treg (red).

Multi-omic integration of differential ATAC-seq, ChIP-seq and RNA-seq analyses. (a-c) Ikaros-dependent changes in accessibility, enhancer signatures, and Foxp3 occupancy at genes differentially expressed in Ikzf1-cko Treg. Open chromatin (top tracks) and H3K27ac (bottom tracks) in WT (blue) and Ikzf1-cko (red) Treg, Ikaros binding sites (purple marks) and Foxp3 binding sites (red marks) in WT Treg, and Foxp3 binding sites in Ikzf1-cko Treg (green marks) at the Ifng (a) and Notch2 (b) and Irf4 (c) loci. (d) Overlap between Ikaros ChIP-seq peaks and differentially acetylated regions is indicated in red in the pie charts. (e) Differential H3K27ac (z-score) at top Ikaros-bound DEG in ex vivo and stimulated WT and Ikzf1-cko Treg.

Foxp3 ChIP-seq library metrics. Assessment of WT and Ikzf1-cko library reproducibility by principal component (a) and pairwise spearman correlation (b) analysis. (c) Strand cross-correlation analysis of FoxP3 ChIP-seq signal vs. noise corresponding to the broad peak (red) and phantom peak (blue; mappability bias). (d) Differential H3K27ac (z-score) at top Foxp3-bound DEG in ex vivo and stimulated WT and Ikzf1-cko Treg. (e) Ikaros-dependent changes in accessibility, enhancer signatures, and Foxp3 occupancy at the Il2ra, Rbpj, and Maml3 genes differentially expressed in Ikzf1-cko Treg. Open chromatin (top tracks) and H3K27ac (bottom tracks) in WT (blue) and Ikzf1-cko (red) Treg, Ikaros binding sites (purple marks) and Foxp3 binding sites (red marks) in WT Treg, and Foxp3 binding sites in Ikzf1-cko Treg (green marks).

(a) Analysis of H&E-stained paraffin sections of liver, lung, kidney, and colon from 1 year-old Ikzf1-fl-Foxp3-YFP-Cre (N=3) or age-matched Foxp3-YFP-Cre control mice (N=3). Representative images at 100X magnification are shown. All tissues showed normal appearance with no evidence of spontaneous autoimmunity. (b) Immunophenotyping of Tconv from Ikzf1-fl-Foxp3-YFP-Cre and Foxp3-YFP-Cre mice. Frequencies of total (left panel), CD62LhiCD44lo naive (middle panel), and CD62LloCD44hi memory (right panel) phenotype Tconv in secondary lymphoid tissues of 10 month-old Ikzf1-fl-Foxp3-YFP-Cre (purple) and Foxp3-YFP-Cre (green) mice (N=3).

Role of Ikaros in Treg-mediated protection from colitis. Spleens and colons from Rag1ko recipients of WT Tconv (a), WT Tconv plus WT Treg (b) and WT Tconv plus Ikzf1-cko Treg (c) harvested at day 40 post-transfer. Quantitation of total CD4+ T cells (d,e), CD4+CD44+ effector Tconv (f,g), and CD4+Foxp3+ Treg (h,i) in the spleens (d,f,h) and mesenteric lymph nodes (e,h,i) from Rag1ko recipients of WT Tconv (orange), WT Tconv plus WT Treg (blue) and WT Tconv plus Ikzf1-cko Treg (green) harvested at day 40 post-transfer. Role of Ikaros in Treg-dependent acquired cardiac transplant tolerance. (j) B6 Ikzf1-fl-Foxp3-YFP-Cre (red) or Foxp3-YFP-Cre (black) mice received BALB/c cardiac allografts, donor-specific transfusion (DST), and anti-CD40L. Graft survival was monitored for 100 days. (k) Histopathological analysis of cardiac grafts harvested at day 14 post-transplant from Foxp3-YFP-Cre and Ikzf1-fl-Foxp3-YFP-Cre recipients (N=3, scale=200x).