Sustained store-operated calcium entry utilizing activated chromatin state leads to instability in iTregs
Abstract
Thymus-originated tTregs and in vitro induced iTregs are subsets of regulatory T cells. While they share the capacity of immune suppression, their stabilities are different, with iTregs losing their phenotype upon stimulation or under inflammatory milieu. Epigenetic differences, particularly methylation state of Foxp3 CNS2 region, provide an explanation for this shift. Whether additional regulations, including cellular signaling, could directly lead phenotypical instability requires further analysis. Here, we show that upon TCR (T cell receptor) triggering, SOCE (store-operated calcium entry) and NFAT (nuclear factor of activated T cells) nuclear translocation are blunted in tTregs, yet fully operational in iTregs, similar to Tconvs. On the other hand, tTregs show minimal changes in their chromatin accessibility upon activation, in contrast to iTregs that demonstrate an activated chromatin state with highly accessible T cell activation and inflammation related genes. Assisted by several cofactors, NFAT driven by strong SOCE signaling in iTregs preferentially binds to primed-opened T helper (TH) genes, resulting in their activation normally observed only in Tconv activation, ultimately leads to instability. Conversely, suppression of SOCE in iTregs can partially rescue their phenotype. Thus, our study adds two new layers, cellular signaling and chromatin accessibility, of understanding in Treg stability, and may provide a path for better clinical applications of Treg cell therapy.
eLife assessment
This study presents the valuable finding that sustained calcium signaling in induced-Treg (iTreg) cells can lead to the loss of Foxp3 expression and iTreg identity by altering the chromatin landscape. The evidence supporting the claims of the authors is convincing. The work will be of interest to immunologists working on Treg cell therapy.
https://doi.org/10.7554/eLife.88874.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- 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
Introduction
The regulatory T cells (Tregs) are key to maintain peripheral immune tolerance and to block excessive inflammatory reactions. There are two main populations of Tregs. The prototypical one originates from the thymus (tTreg), sharing a developmental program similar to conventional T cells (Rudensky, 2011; Sakaguchi et al., 2008). Their TCRs are believed to have higher basal affinity to MHC (major histocompatibility antigen), and express Foxp3 during CD4 single positive period (Fontenot et al., 2005). The other population originates from the periphery (pTregs) via conversion of Tconvs, arising during inflammation with proper antigenic stimulation and inductive cytokine environment (Curotto de Lafaille et al., 2004; Apostolou and von Boehmer, 2004; Cobbold et al., 2004). iTregs, mimicking the induction of pTregs, is a type of Tregs induced in vitro in the presence of TGF-β (tumor growth factor ) (Chen et al., 2003; Fantini et al., 2004). While iTregs share similarity in suppressive capacities to Tregs in vivo (Shevach and Thornton, 2014; Fu et al., 2004; Singh et al., 2009), they are different from tTregs in several aspects, particularly regarding their stabilities. Although tTregs in some occasions show certain levels of plasticity, their phenotypes are mostly stable, both in vivo and in vitro (Rubtsov et al., 2010). iTregs induced in vitro are more labile with marked loss of Foxp3 expression after activation, limiting the duration of their suppression (Floess et al., 2007; Chen et al., 2011). The short life of iTregs is also a roadblocker of the clinical use of Treg cell therapy, as they cannot be used as an alternative source to replace numerically limited tTregs (Koenecke et al., 2009); and in cases where iTregs are raised with targeted specificities, they carry concerns of potential self-reactivity following the loss of Foxp3.
The stability of tTregs is established prior to their exit from the thymus (Ohkura et al., 2013). In the thymus, Foxp3 works with other transcription factors, such as NFAT and c-Rel to maintain a steady transcriptional program (Fu et al., 2012; Hill et al., 2007; Samstein et al., 2012). At the same time, the extensive demethylation of Foxp3 CNS2 is also critical to their stable phenotype, an event likely taking place before Foxp3 expression (Samstein et al., 2012; Ohkura et al., 2012). While epigenetic differences provide a possible explanation of Treg stability, it is more probable that multiple factors collectively contribute to the difference. Unlike Tconvs that undergo expansion, attrition, and memory states after activation, tTregs lose signature changes and return to their original state following activation (van der Veeken et al., 2016). This persistence in phenotype of tTregs is certainly beneficial to the host, yet mechanisms for this maintenance, in contrast to the lack thereof in iTregs are incompletely understood. Previous work suggests that TCR signal as well as ensuing Ca2+ activities are weaker in tTregs (Sumpter et al., 2008; Vaeth et al., 2012; Yan et al., 2015), whether the reduced TCR signaling strength is the root of tTreg stability remains an intrigue, particularly strong TCR engagement facilitates iTregs’ loss of phenotype.
Here, we report a surprising finding that may contribute to the iTregs’ phenotypical change after their activation. We found that ligation of TCR induces long vibration of Ca2+ signals in Tconvs that lasts for hours, as typically expected in the store-operated calcium entry (SOCE) activation. Initially similar, this vibration is quickly diminished in tTregs, correlating with their reduced response to TCR stimulation. In contrast, in iTregs, this Ca2+ vibration is indistinguishable from Tconvs. The sustained Ca2+ signal leads to a strong NFAT activation in iTregs. Regarding chromatin accessibility, with the aid of our initial observation of iTregs from the perspective of activated Tconv, we found iTregs also show an ‘activated’ chromatin state with T cell activation and inflammation related genes highly open. We demonstrate that the NFAT nuclear translocation driven by sustained calcium signal allows its direct access to prime-opened TH genes and upregulate their expression, destabilizing iTregs. As expected, blocking CRAC channel reverses the loss of Foxp3 in iTregs in vitro and in vivo. Our work therefore suggests that Ca2+ signaling and accessible chromatin state in iTregs are two key causes in their loss of stability following activation, which may be of important value to understand how these two populations of Tregs operate in the disease conditions and lead to better applications of iTregs in the clinic.
Results
iTregs use a suppression mechanism similar to that of tTregs
In our previous work, we found that Tregs show strong binding to DCs, which is a key mechanism for blocking DC’s ability to engage other Tconvs (Chen et al., 2017). The strong binding is due to low basal calcium oscillation of tTregs in the steady state (Wang et al., 2022), which leads to diminished mCalpain activation and commensurate reduced internalization of LFA-1. Furthermore, our previous study showed the low calcium oscillation is maintained by reduced expression of endosomal Ca2+ channel Ryr2 by Foxp3 (Wang et al., 2022). To ascertain that this regulation is also central to iTregs, we produced iTregs from Tconvs. With plate-bound anti-CD3/CD28 and in the presence of TGF-β, iTregs were induced and the percentage was maximized on day 2. We found that in vitro iTregs show stronger suppression than tTregs on OT-II proliferation in response to epitope peptide loaded DCs (Figure 1A). In accordance with the stronger suppression, single-cell force spectroscopy analysis showed that iTregs possessed apparent binding to DCs, even stronger than tTregs (Figure 1B). In addition, iTregs also showed increased binding force compared with activated Tconvs (Figure 1C), suggesting the strong binding force is a common feature of both types of Tregs, but not Tconvs. We measured the expression of Foxp3 more precisely, and noticed that the protein-level expression of Foxp3 reached its peak on day 2 after the induction (Figure 1D). With the accumulation of Foxp3 protein, basal Ca2+ oscillation in the treated Tconvs was gradually diminished, with the former leading the latter roughly by 1 day (Figure 1E). As expected, the transcription level of Ryr2 was suppressed both in iTregs and tTregs (Figure 1—figure supplement 1A), suggesting that the strong binding force with DC may come from reduced calcium oscillation as a result of Foxp3-mediated suppression of Ryr2 transcription. While Ryr2 expression may be regulated by Foxp3, the induction efficiency of Foxp3 was not affected by Ryr2 knock out (Figure 1—figure supplement 1B). Therefore, we show iTregs can apply a suppression mechanism similar to tTregs, which operates via the indirect inhibition via DC blockage. We do not know why iTregs show a stronger suppression than tTregs, an event maybe related to their activation in the induction process.
Sustainment of SOCE and NFAT nuclear translocation are reduced in tTregs, but not in iTregs
It was reported previously that TCR signaling led to a dramatic reduction of Foxp3+ cells in iTregs, but not in tTregs (Chen et al., 2011; Selvaraj and Geiger, 2007). We reproduced these data in our laboratory (Figure 2A). tTregs have been reported to show an attenuated response to TCR triggering (Sumpter et al., 2008; Vaeth et al., 2012; Yan et al., 2015), the exact point of signal attenuation in tTregs has not been elucidated. In immune cells, there are mainly three types of Ca2+ signals, the steady-state low-level Ryr2-based membrane proximal Ca2+ ‘puff’; the TCR ligation-triggered, inositol triphosphate-mediated IP3R opening on the ER (endoplasmic reticulum) endothelia membrane; and ER Ca2+ depletion-induced, sustained Ca2+ influx upon SOCE activation (Guse, 1998). As the extracellular Ca2+ is roughly 4 logs higher than the intracellular, SOCE activation is by far the strongest Ca2+ event in T cells (Figure 2—figure supplement 1A), and is also distinct from the other Ca2+ signaling patterns for its persistent vibration lasting hours to days (Christo et al., 2015; Wülfing et al., 1997).
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/88874/elife-88874-fig2-data1-v1.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/88874/elife-88874-fig2-data2-v1.zip
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/88874/elife-88874-fig2-data3-v1.zip
We wondered if the reported TCR signal truncation in tTregs was manifested as reduced SOCE activities. We titrated anti CD3/CD28 concentrations, and iTregs lost Foxp3 at the lowest dose required to activate, yet tTregs were resistant to an eightfold increase in both (Figure 2—figure supplement 1B). To study the molecular basis for the weak response in tTregs, we first monitored the early SOCE intensity in iTregs and tTregs, we found the first wave of SOCE in iTregs was slightly higher than Tconvs and tTregs, but there was no discernable difference between tTregs and Tconvs (Figure 2B). In addition, the phosphorylation of CD3ζ was similar in Tconv and nTreg immediately after TCR crosslink (Figure 2—figure supplement 1C). As full T cell activation requires sustained SOCE activation (Dolmetsch et al., 1998), we extended the observation time window. As SOCE signal is synchronized upon initial TCR trigger, each individual T cell SOCE oscillation hours into activation is no longer in unison (Christo et al., 2015; Le Borgne et al., 2016), we therefore recorded Ca2+ waves in individual cells. In Tconvs and iTregs, Ca2+ remained cyclic waves expected in SOCE pulsation. In contrast, this Ca2+ signal was suddenly lost in tTregs 40–50 min after the activation (Figure 2C; Videos 1–3).
In contrast to other transcriptional factors, such as NF-κB, activated rapidly in response to short Ca2+ elevation, NFATs are known to be triggered only by sustained cytoplasmic Ca2+ (Dolmetsch et al., 1998; Dolmetsch et al., 1997; Tomida et al., 2003). Indeed, in contrast to iTregs and Tconvs, tTregs’ NFATc1 and c2 failed to translocate into the nucleus (Figure 2D). Upon detailed analysis, in iTregs, anti-CD3/CD28 stimulation led a quick drop in cytoplasmic NFAT in 2 hr, with moderate accumulation of NFAT in the nucleus (Figure 2—figure supplement 1D). This was followed by continuous synthesis of NFAT in the cytoplasm and accumulation of nuclear NFAT (Serfling et al., 2006; Zhou et al., 2002). While tTregs are known to have reduced NFAT activation in response to TCR ligation (Sumpter et al., 2008; Nayak et al., 2009), our results pointed to a truncated Ca2+ response likely to be the molecular basis for this phenotype. With this association potentially established, we wondered if the length of SOCE itself is sufficient to explain the stability difference between iTregs and tTregs. Non-specific Ca2+ membrane permeability modifier such as ionomycin resulted in large-scale cell death (Figure 2—figure supplement 1E), we therefore used SERCA inhibitor Thapsigargin (TG). Blocking SERCA leads to endosome Ca2+ reloading blockage that leads to continuous translocation of STIMs to the plasma membrane–ER junction and Orai channel opening. In contrast to indiscriminating Ca2+ overload, pure SOCE activation was not followed by significant cell death (Figure 2—figure supplement 1E). As expected, some tTregs treated with TG lost their Foxp3 (Figure 3E). These results collectively argue that the truncated Ca2+ in tTregs is essential to protect tTregs from phenotypical shift.
iTregs display activated chromatin state with highly open activation region
As shown, the behaviors of SOCE signaling and NFAT of iTregs were closer to their precursor Tconvs rather than tTregs, we were curious about the systemic similarity among iTregs, Tconvs, and tTregs. Therefore, we performed RNA-seq and Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) of iTregs, together with resting/activated tTregs and Tconvs.
In RNA-seq data, principal component analysis (PCA) and heatmap of DEG (differentially expressed gene analysis) indicated that Tconvs were quite different before and after activation, but the changes of tTregs before and after TCR activation were more limited compared with Tconvs (Figure 3A and Figure 3—figure supplement 1C), verifying that tTregs have low TCR responsiveness. Comparing the transcriptome characteristics with other cell types, iTregs acquired partial Treg features after induction, shown as an intermediate between Tconvs and tTregs on PC2 (Figure 3A and Figure 3—figure supplement 1B). While among all cell states, iTregs had the highest PC1 scores, represent higher proliferation state, probably due to the 4 days of induction (Figure 3A and Figure 3—figure supplement 1A).
Next, we analyzed the chromatin accessibility of these cells using ATAC-seq data. In the PCA plot, tTregs showed fewer changes after activation compared to Tconvs, while iTregs were closer to activated Tconv cells (Figure 3B). In differential peak heatmap clustered by cell type, Tconvs became more accessible in ‘Activation Region’ after activation, and with a dramatic increase in the overall accessibility (Figure 3C). However, the changes of tTregs after activation are limited, the majority of tTreg-specific accessible region, named ‘Treg Region’, were not changed. Also the activation region opening seen in Tconvs were not more accessible in tTregs (Figure 3C). All of these characteristics indicated that tTreg phenotype is highly ensured by their low TCR responsiveness which potentially led to their refusal to change chromatin accessibility seen in Tconvs following activation.
When comparing iTregs with the others, we found iTregs were highly accessible in activation region found in activated Tconvs (Figure 3C, D). However, only a few Treg-specific peaks shown in Treg region were opened in iTregs (Figure 3C, E). For further analysis, we isolated Treg-specific peaks. It was clear that most of them were opened in tTregs but showed little accessibility in iTregs (Figure 3—figure supplement 1D). For example, Foxp3 and Ctla4, expressed in both tTregs and iTregs, show high accessibility in resting and activated tTregs, but were closed or only partially opened in iTregs (Figure 3G). In contrast, genes related to inflammation cytokine in T helper cells, such as Ifng, Il4, Il17ra, and Il21, were highly opened in iTregs (Figure 3F). In further elaboration of Activation Region in Figure 3D, we have selected some genes including typical peaks that are within this region (Figure 3—figure supplement 1E). These genes encompass some T cell activation-associated transcription factors, such as Irf4, Atf3, as well as multiple members of the Tnf family including Lta, Tnfsf4, Tnfsf8, and Tnfsf14. Additionally, genes related to inflammation cytokine and function such as Il12rb2, Il9, and Gzmc are included. These genes display elevated accessibility upon T cell activation, partially open in activated tTreg cells. However, all of them exhibit high accessibility in iTreg cells.
Overall, we showed that at the gene expression level, iTregs have acquired some Treg features, but the chromatin accessibility was mostly trapped in activated Tconv state, not only with regard to the closed Treg genes, but also with genes associated with T activation and inflammation highly opened. The latter is expected to be mainly closed in genuine Tregs in accordance with their lower responsiveness to TCR triggering. We speculate that this ‘unconverted’ chromatin accessibility, especially the highly opened activated genes will contribute to iTreg instability upon TCR stimulation. It is therefore interesting whether NFAT accumulated in nuclei driven by sustained calcium oscillation will indeed target those genes for an activation state unfit for their phenotypical stability.
SOCE and NFAT directly leads to phenotypical instability in iTregs
To further validate the role of calcium and NFAT in iTreg instability, we used inhibitors cyclosporine A (CsA) and CM-4620 to suppress calcium signal and NFAT, respectively. We found that inhibition of calcium or NFAT rescued the Foxp3 expression in iTregs upon restimulation (Figure 4A). Another SOCE inhibitor BTP2 was also verified (Figure 4—figure supplement 1A). Accordingly, ionomycin, which can directly induce calcium release independent of TCR activation, also lead to iTreg instability (Figure 4—figure supplement 1B). Thus, NFAT and calcium signaling blockage were correlated to iTreg stability. As control, inhibition of NF-κB (BAY 11-7082), c-Jun (SP600125), or a c-Jun/c-Fos complex (T5224) had no discernable effect, or in one case, possibly further reduction in stability (Figure 4—figure supplement 1C). These results may indicate that NFAT plays a crucial and special role in TCR activation, which leads to iTreg instability.
To investigate how NFAT and calcium signaling destabilize iTregs, we performed RNA-seq on these cells. From the differential gene expression and PCA, we found that iTregs upon restimulation showed significant differences from the resting state at the transcriptome level, and this difference could be substantially reduced with inhibition of calcium or NFAT (Figure 4B). As shown in the PCA plot, calcium/NFAT-blocked iTregs were significantly closer to the resting state than the unstable ones (Figure 4C). Thus, inhibition of calcium or NFAT can significantly improve the stability of iTregs as reflected by transcriptomic identity. In addition, the large overlap (about 80%) of genes that were rescued by CM-4620 and by CsA further supported that calcium signaling change the fate of iTregs mainly through NFAT (Figure 4D).
Among the genes rescued by CM-4620 and CsA, they can be further divided into two groups. Members of the first group were downregulated upon restimulation and the downregulation was blocked by inhibition of calcium and NFAT. This result indicating that these genes are downregulated by NFAT. Among these genes, several of them are critical for functions and stability of Tregs, such as Foxp3, Id3, Il10rb, Sema4a, Tgfbr1/2, Ms4a4b, and Nlrp6 (Figure 4B and Figure 4—figure supplement 1D). In the functional analysis of these genes, we found that Rap-1 pathway, PD-1/PD-L1 pathway and Foxo pathway were significantly enriched (Figure 4—figure supplement 1E). Loss of these genes was therefore consistent with the phenotype of iTreg instability.
The other group had genes that were otherwise upregulated after restimulation, but failed to do so with CM-4620 or CsA, indicating that these genes are upregulated by NFAT. In these genes, many of them are associated with T helper differentiation and T cell activation, such as Il21, Il12rb2, Tbx21, Gzma, Stat2, Stat3, Atf3, and Fasl (Figure 4B, F). Among them, Il21 encodes a key cytokine in Th17 differentiation, Tbx21, encoding T-bet, and Il12rb2 are key transcriptional factors and cytokine receptor for Th1 differentiation. Functional enrichment analysis of these upregulated genes pointed to Th17 related pathway and Th1/Th2 related pathways (Figure 4E). Further validation of these upregulation T helper related genes (TH genes) by NFAT was perform by QPCR (Figure 4G). These results in all indicated the iTreg restimulation was associated with a strong tendency of differentiation toward non-Treg polarization, and NFAT was likely a key regulator in this process.
NFAT upregulate prime-opened TH genes to destabilize iTregs
To further study the role of NFAT in iTreg instability, we performed Cut&Tag assay to acquire the NFAT-binding sites in resting iTregs and those after restimulation (Figure 5A). NFAT-flag and anti-Flag antibody were used to capture peaks that were bound by NFAT (Figure 5—figure supplement 1A). The binding peaks were of high quality as indicated by the lower number of peaks in the mock control compared to both resting and activated states (Figure 5B). The number of NFAT peaks in restimulated iTregs was significantly higher than that in resting state, which was consistent with the massive nuclear translocation of NFAT indicated by western blot (Figure 5B and Figure 2D). For peak distribution at gene body, NFAT peaks both before and after restimulation were enriched at promoters/transcription start site (TSS) and the first introns of genes (Figure 5—figure supplement 1C).
For those upregulated TH genes, we speculated that this may be associated with the highly ‘activated’ chromatin state of iTregs mentioned earlier. Therefore, we first checked the chromatin accessibility around the TSS in those genes and found that this group of genes was significantly more accessible in iTregs than in tTregs and Tconvs (Figure 5C). As shown in the genome browser view, Il21, Il12rb2, and Tbx21 were more open in iTregs, even more so than the activated Tconvs (Figure 5—figure supplement 1B). Further analysis revealed that NFAT significantly increased their binding to the TSS regions on upregulated Th genes after restimulation (Figure 5D), while showing little change on downregulated Treg genes (Figure 5—figure supplement 1D). As shown by Il21, Tbx21, and Il12rb2, after activation, NFAT accumulated at the sites that were prime-opened resulting in their increased expression (Figure 5E). Therefore, the enhanced chromatin accessibility revealed in Figure 4 was indeed approached by the large amount of NFAT, leading to T cell polarization toward other directions with a sole reciprocal loss of polarization toward Tregs. This tendency may to some extent explain the rapid loss of iTreg phenotype.
As one potential counter argument, NFAT is ostensibly essential to the development of tTregs (Li et al., 2014; Zheng et al., 2010). In the induction of iTregs, NFAT in collaboration with SMAD is also critical (Vaeth et al., 2012; Zheng et al., 2010; Tone et al., 2008). It is therefore intriguing why NFAT re-activation after iTreg formation may serve as a key factor in compromising iTreg stability. To resolve those seemingly contradictory regulations, we performed NFAT motif enrichment analysis. The results showed that putative BATF, AP-1(Fos/Jun), and RORg-binding sites were highly enriched in restimulation state specific regions (Figure 5F). For the preference of different cofactors to different gene subsets, SMAD family tended to bind to genes that were downregulated by NFAT, for example, Foxp3, while AP-1 and RORg preferred binding to genes that were upregulated by NFAT, such as TH genes (Figure 5G). Taken together, our results appear to suggest that upon restimulation of iTregs, NFAT prefers to cooperate with AP-1/RORg to promote the expression of prime-opened TH gene expression, and at the same time denies the SMAD programming that may be critical to maintain Treg functions and stabilities. Certainly, how much this selective cooperation between NFAT and other transcriptional factors contributes to iTreg instability remains speculative until large-scale ChIP-seq analyses are performed on the players involved.
According to the Cut&Tag data, the secretion levels of IL-21 in tTreg and iTreg were also examined by ELISA. As shown in Figure 5H, tTreg did not secrete IL-21 regardless of activation status (undetectable), while iTreg did not secrete IL-21 at resting state but they did so after 48 hr of restimulation (Figure 5H). The secretion of IL-21 was inhibited by CsA and CM-4620 treatment. This observation aligns with our findings where nuclear binding of NFAT to gene loci of these cytokines enhanced their expression, pushing iTreg unstable under inflammatory conditions. These findings further underscore the likelihood that the inhibition of calcium and NFAT signaling might contribute to the stabilization of iTreg by suppressing the secretion of inflammatory cytokines.
Reduction of SOCE enhances iTreg stability and suppressive potential
Given that sustained calcium signal and NFAT translocation can directly disrupt iTreg stability, we asked whether manipulation of calcium signaling or NFAT can conversely be used to stabilize those cells. As Treg suppression assays require the participation of other T cells which are sensitive to inhibitory agents added indiscriminately, after regular iTreg induction, we overexpressed a dominant negative ORAI (DN-ORAI), whose single monomer incorporation will lead to a functionally disabled ORAI hexamer, to control the SOCE signal upon iTreg restimulation (Mignen et al., 2008). As expected, the mutated ORAI reduced the SOCE signaling by about a half (Figure 6A). When we restimulated this DN-ORAI expressing iTregs with anti-CD3/CD28, they showed stronger stability with higher retention of Foxp3 than WT-ORAI iTregs (Figure 6B). Consistently, DN-ORAI iTregs displayed more potent suppressive activity than control iTreg in in vitro suppression assay (Figure 6C).
Finally, the stability of DN-ORAI iTregs in vivo was examined. In brief, OT-II iTregs were infected with DN-ORAI retrovirus and then purified by GFP fluorescence, then the CD45.2+ GFP+ Foxp3-RFP+ iTregs were adoptively transferred into CD45.1 recipient, and OVA/Alum were administrated to restimulate DN-ORAI iTregs in vivo. Analysis of the percentage of Foxp3+ cells in mLN 5 days after adoptive transfer showed that DN-ORAI iTregs were more resistant to the loss of Foxp3 and thus were more stable than control WT-ORAI iTregs (Figure 6D). These results collectively showed that weakened calcium signal and inhibition of NFAT can quantitatively prevent iTregs from losing more Foxp3 both in vitro and in vivo, thus permitting a stronger suppressive capacity and a longer stability.
Based on extensive opening of activation regions approached by NFAT in the process iTreg stimulation, we postulate that openness of activation region could be another crucial feature of Treg stability, in addition to closed Treg region by low TSDR (Treg-specific demythylated region). Beyond the basal regimens of anti-CD3, anti-CD28, IL-2, and TGF-β, several additional factors have been reported to enhance the stability of iTregs, including the use of retinoic acid, rapamycin, vitamin C, and AS2863619 as well as removal of CD28 signaling (Akamatsu et al., 2019; Hill et al., 2008; Mikami et al., 2020; Yue et al., 2016; Zhang et al., 2013). We observed if those new additions in some ways altered the state of TSDR as well as the the openness of activation region. Our findings revealed that rapamycin and AS tended to close the activation region (Figure 6—figure supplement 1A), while Vc, removal of CD28, and AS increased the accessibility of Treg genes (Figure 6—figure supplement 1B). These results indicate that different treatments have varied impacts on the activation region and Treg region. Hence, the selective and combinatorial use of those factors is an interesting proposal in the generation of more stable iTregs.
Discussion
tTregs are different from other subsets of CD4 T cells in the thymic determination of their phenotype. To some extent, this centric determination is closer to CD4/CD8 demarcation, rather than other polarizations in the peripheral environment (Rudensky, 2011; Sakaguchi et al., 2008). The stable phenotype, even after activation, is critical to provide the baseline immune suppression (Rubtsov et al., 2010). pTregs/iTregs are driven by the environment. It is logical that those cells should be deprived of their suppressive ability when are no longer needed to avoid unnecessary accumulation of those cells. The current understanding is that during thymic development, Foxp3 CNS2 region is selectively demethylated, allowing several transcriptional factors, such as Ets1 and Foxp3 itself to binding to this location (Samstein et al., 2012; Ohkura et al., 2012). This is one of several epigenetic features that maintain the stability, and TCR signaling in tTregs does not change the pre-established hypomethylation state. Induced Tregs, on the other hand, use a different program that relies on the binding of SMADs to the CNS1 region to initiate Foxp3 expression (Tone et al., 2008). This SMAD binding mediated transcription is dependent on Foxp3 H3K4me3 modification. In the absence of SMAD, H3K4me3 is partially lost, resulting in iTreg phenotypical shift (Tone et al., 2008). Withstanding this paradigm, iTregs do not immediately lose Foxp3 expression upon TGF-β withdrawal, the loss of Foxp3 is evident hours or days after TCR signaling (Selvaraj and Geiger, 2007). This indicates that while the epigenetic alterations are necessary, an activation event is required to drive the loss of phenotype. This notion is subtly echoed by the limited TCR signaling in tTregs, implying that the hyporesponsiveness to TCR ligation may be a critical feature to maintain their stability.
We observed differences between tTregs and iTregs in calcium signal and Th gene accessibility. tTregs utilize both truncated calcium signal and conserved chromatin accessibility to systematically reduce their TCR responsiveness, ensuring their stability. In contrast, iTregs retain the full strength of the TCR signaling, particularly the sustained vibration of SOCE signal, the key factor that drives NFAT synthesis and nuclear translocation. With the help of several cofactors, such as AP-1, NFAT caused the high expression of these TH genes, which subsequently subverts the original differentiation of iTregs, leading to their instability. Direct inhibition of SOCE signaling and NFAT can enhance iTreg stability in vitro and in vivo (Figure 6E).
Previous observations also showed that the calcium signal in tTreg is weaker than that in Tconvs, and NFAT activity is suppressed in tTregs. The weaker Ca2+ signal in previous reports was seen as early as 10 min into activation, this contradiction likely reflects the experimental settings and stimulation strength (Yan et al., 2015). Nevertheless, in line with our results, inhibition of NFAT by low-dose CsA increased the number of tTreg in atopic dermatitis patients (Brandt et al., 2009). It seems reasonable to establish an inverse association between SOCE signal and Treg stability. In our analysis of the RNA-seq and chromatin accessibility data, we did not identify any particular mechanism that may account for the distinct Ca2+ signal patterns in tTregs and iTregs. However, in contrast to Tconvs and iTregs, there must be a negative signal loop that blunts SOCE vibration in tTregs. We tried to sustain the calcium signaling in tTregs forcibly by thapsigargin, and indeed observed loss of Foxp3. However, as we at the moment cannot appreciate the finesse of this regulation, the underlying molecular details will require much more delicate approaches rather the brutal force of ER store depletion. In our view, future investigation in this direction may lead to better understanding of Treg biology as well as innovations in clinical Treg cell therapy.
Our results are consistent with the understanding that Treg-specific hypomethylation is crucial for tTregs to stably express their Treg-specific genes (Samstein et al., 2012; Ohkura et al., 2012). Those epigenetic changes are mostly absent in iTregs, resulting in their instability. The differences of accessibility in Treg genes were also observed in our ATAC-seq data. Therefore, current treatments such as vitamin C and CD28 blockage applied to achieve demethylation and accessibility of Treg genes in iTregs are worthy attempts (Mikami et al., 2020; Yue et al., 2016). However, our work indicates that at least in ATAC-seq analyses, there are a cluster of activation and differentiation-related peaks present in iTregs that remain silent in tTregs upon TCR stimulation. Critically, these peaks were targeted by abundant NFAT in the iTreg nuclei following activation, leading to the upregulation of Th genes. This tendency of deviation away from Tregs may be another key player to iTreg instability. Our effort to produce highly stable iTregs would likely require a combined approach to demethylate Treg-specific genes and at the same time to block the activation and differentiation-related genes. In our work, targeting NFAT signaling during the restimulation phase appeared to be a potential method. Other ways that can directly close these genes will also likely be beneficial.
In this work, we stress the importance of NFAT signal in the instability of iTregs. Previous research has also unveiled the adverse effects of other TCR signaling pathways on Treg stability, such as NF-κB and mTOR. As all those factors may contribute to the shift in phenotype, it is out of scope of this report to directly compare all those factors. However, the value of this report is the identification of a key difference in Ca2+ signaling under TCR ligation that may explain the stability differences between iTregs and tTregs, which may open the door for additional regulations given the key role of SOCE activity in T cell activation.
Methods
Mice
All mice were in CD45.2 C57BL/6 background unless noted otherwise. OT-II-transgenic mice and CD45.1 mice were purchased from Jackson Laboratories. Foxp3-IRES-GFP transgenic mice (B6.Cg-Foxp3tm2Tch/J, Strain #:006772 from Jackson Laboratory) were a gift from Dr. Hai Qi, School of Medicine, Tsinghua University. Foxp3-RFP mice were a gift from Dr. Xuyu Zhou, Institute of Microbiology, Chinese Academy of Sciences (CAS). Wild-type mice were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd. All mice were bred and housed at Tsinghua University Animal Facilities and maintained under specific pathogen-free conditions. This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved Institutional Animal Care and Use Committee (IACUC) protocols of the Tsinghua University. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Tsinghua University (22-SY6). All surgery was performed under sodium pentobarbital anesthesia, and every effort was made to minimize suffering.
Antibodies and reagents
Recombinant human IL-2 and human TGF-β were from R&D systems. Anti-mouse CD3 and biotin-anti-CD3 (145-2C11) mAb, anti-mouse CD28, and biotin-anti-CD28 (37.51) mAb purified streptavidin were all from biolegend. LIVE/DEAD Fixable Aqua were from Invitrogen. For flow cytometric analysis, in addition to CD45.1 (A20) monoclonal antibody from Invitrogen, others were purchased from eBioscience or BD Pharmingen. The following clones were used: Foxp3 (150D/E4), CD45.2 (104), CD4 (GK1.5), CD25 (PC61.5), CD44 (IM7), and CD62L (MEL-14). Antibody used in immunoblot including NFATc1 Antibody (7A6) and NFATc2 Antibody (4G6-G5) were from Santa Cruz, Lamin A/C (4C11), β-Actin (13E5), and GAPDH (D4C6R) mAb were from CST. Monoclonal ANTI-FLAG M2 antibody used in Cut&Tag assay was from Sigma. Mouse CD4 T Cell Isolation Kit used for the enrichment of CD4+ T cells was from Stem cell. Foxp3/Transcription Factor Staining Buffer Set for Foxp3 staining was from eBioscience. NE-PER Nuclear and Cytoplasmic Extraction Reagents used for immunoblot was from Thermo. All sequencing kit, including TruePrep DNA Library Prep Kit V2, Hyperactive Universal CUT&Tag Assay Kit, and TruePrep Index Kit V2, was from Vazyme. For calcium mobilization, Ionomycin was from Beyotime and thapsigargin from Invitrogen. For inhibitors, cyclosporin A, CM-4620, and BTP2 were all from MCE. OVA323–339 peptides from MBL. For calcium imaging, Fluo-4 AM was from Thermo and Indo-1 AM from BD.
Flow cytometry sorting and FACS
For surface marker detection, cells were incubated with Fc blocker (CD16/32 antibody; 2.4G2) for 5 min, and then incubated with surface antibody for 15 min at room temperature avoiding light. Foxp3 intracellular staining was performed by cell fixation and permeabilization using the Foxp3/Transcription Factor Staining Buffer Set. Samples were analyzed on Fortessa cytometers (BD) with FACS Diva software and data were analyzed with Flowjo.
To prepare Tconv, tTreg, and iTreg for experiments, spleen was isolated from 2- to 4-month-old Foxp3-GFP mice for Tconv and tTreg sorting, and 6-week-old mice for iTreg induction. Subsequently, CD4+ T cells were enriched using CD4 Isolation kit and stained with indicated surface markers. FACSAriaII (BD) was used for collecting particular population. The cell populations used are as follows: naive T cells, CD4+ Foxp3-GFP+ CD44low CD62Lhigh; tTreg cells, CD4+ Foxp3-GFP+; iTreg, Foxp3-RFP+ or Foxp3-GFP+. Infected iTreg, Foxp3-RFP+ MSCV-GFP+.
In vitro iTreg induction and restimulation
CD4+Foxp3-GFP− CD62LhiCD44lo (or Foxp3-RFP−) naive T cells were sorted using 6-week-old mice and stimulated with plate-bound 0.5 μg/ml anti-CD3 and 1 μg/ml anti-CD28 in the presence of 2 ng/ml human TGF-β and 200 U/ml rhIL-2. iTreg was induced for 4 days and sorted using Foxp3-GFP/RFP fluorescence, then cells were rested in fresh RPMI complete medium containing 200 U/ml rhIL2 for 1 day.
For iTreg restimulation, iTreg was harvested plated at 0.4 × 106 cells for restimulation with plate-bound 0.5 μg/ml anti-CD3 and 1 μg/ml anti-CD28 in 96 plate. After 2-day restimulation, Foxp3 was analyzed by intracellular staining.
Calcium measurement
For calcium influx within 10 min upon TCR activation, sorted Foxp3-RFP− Tconv, Foxp3-RFP+ tTreg and Foxp3-GFP+ iTreg were resting for 1 day, then mixed and loaded with 1 μM Indo-1 AM at 37°C for 30 min in RPMI medium without fetal bovine serum (FBS). Cells were then washed for twice in cold RPMI medium containing 1% FBS, followed by staining with biotinylated anti-CD3 and biotinylated anti-CD28 on ice for 1 hr. For calcium measurement, cells were heated in 37°C water bath for 5 min, then the baseline was recorded for 60 s on LSR II(BD). Calcium influx was triggered by the addition of 100 μg/ml streptavidin to crosslink TCR. Data were analyzed using kenetics in FlowJo.
For calcium oscillation within 1 hr, Foxp3-RFP− Tconv, Foxp3-RFP+ tTreg, and Foxp3-RFP+ iTreg were sorted and rested for 1 day, then cells were loaded with Fluo-4 AM at 37°C for 30 min. After washing twice with phosphate-buffered saline and further incubating in RPMI medium for 10 min, cells were activated by plate-bound 0.5 μg/ml anti-CD3 and 1 μg/ml anti-CD28 in poly-L-lysine-coated confocal dish. Excess non-adherent cells were removed by flushing with buffer solution after 10 min. Fluorescence were recorded in the indicated time after stimulation as a time lapse for 20 min with an interval of 6 s. The mean fluorescence intensity changes over time in single cells were analyzed using Surface in Imaris and normalized to resting fluorescence F0 (Fluo-4 F/F0), further the standard deviation of F/F0 in single cells was calculated and presented using GraphPad.
Immunoblot
iTreg, tTreg, and Tconv were rested for 1 day in RPMI medium containing rhIL-2 and treated with CsA for 30 min before activated. Then these cells were stimulated by plate-bound 0.5 μg/ml anti-CD3 and 1 μg/ml anti-CD28 for indicated time, in the absence or presence of TGF-β. After stimulation, cells were immediately lysed using CERI, then nuclear and cytoplasmic component were separated using NE-PER Nuclear and Cytoplasmic Extraction Reagents. After being mixed with sodium dodecyl sulfate loading buffer and boiled for 10 min, the cytoplasmic and nuclear proteins equal to 5 × 105 iTreg and 2 × 10 6 tTreg/Tconv were loaded onto 7.5% PAGE Gels. Then the protein was transferred onto a NC membrane and immunoblotted with NFATc1 and NFATc2, Actin and GAPDH were used as loading control of cytoplasmic protein, and LaminA/C as nuclear. Finally, the immunostained bands were detected by Super ECL Detection Reagent.
Treg cell suppression assay in vitro
104 purified DCs from splenocytes were pulsed with 2 μg/ml OVA323–339 peptide, then iTreg and tTreg (2 × 104) were added onto DCs to occupy the latter for 30 min. CD25− OTII Tconv cells were sorted and stained by CFSE (Thermo Fisher Scientifc). 2 × 104 OTII Tconv cells were mixed in DC-Treg culture in a 96-well U bottom plate. The proliferation of OTII T cells was assessed after 2–4 days by CFSE dilution on flow cytometry.
DN-ORAI and NFAT-flag construction and retrovirus infection
For DN-ORAI, ORAI CDS sequences were amplified from cDNA library and mutated with primer. Then WT-ORAI and DN-ORAI were linked with MSCV-GFP vector. For NFAT-Flag, Nfatc1-flag CDS was amplified from TetO-FUW-flag-Nfatc1 plasmid, and further linked with MSCV-GFP vector. For retrovirus production, platE cells were cultured with Dulbecco’s modified Eagle medium containing 10% FBS to 80–90% confluence. For each 10 cm dish, 10 µg MSCV plasmid and 10 µl Neofect transfection reagent were mixed in 500 µl Opti-MEM medium. After 30-min incubation, the DNA mixture were added into the well and gently mixed. After 48 and 72 hr, supernatant containing viruses was collected, filtered with 0.45 µm membrane. 8 µg/ml polybrene was added to virus supernatant, then Foxp3-RFP+iTreg was infected on the days 2 and 3 of induction by centrifugation at 1500 × g, 37°C 2 hr. On the day 5 of induction, the infected cells with GFP fluorescence and Foxp3-RFP iTreg was sorted and used.
In vivo iTreg stability
CD45.2+ Foxp3-RFP+ iTreg was sorted, and washed twice. 5 × 106 cells were adoptively transferred into CD45.1 mice via tail vein. For DN-ORAI iTreg, CD45.2+ Foxp3-RFP+ iTreg was infected and GFP+RFP+ iTreg was sorted to transfer to CD45.1 recipient. Six hours after cell transfer, mice were immunized with OVA protein mixed with Alum adjuvant in each hind flank. After 4 days of immunization, mLN were collected and the percentage of Foxp3+ cells in CD45.2+ was analyzed by intracellular staining.
AFM-based single-cell force spectroscopy
The experiments were performed as previously described using a JPK CellHesion unit (Chen et al., 2017; Flach et al., 2011). In brief, to measure Treg-DC adhesion forces in bicellular system, DC2.4 cells were cultured on untreated glass disks. iTreg/tTreg/Tconv cells were sorted and rested with rhIL2 overnight. The disks were moved into an AFM-compatible chamber and mounted on to the machine stage. The incubator chamber in which the machine was housed was conditioned at 37°C and at 5% CO2. A clean cantilever was coated with CellTak (BD), and then used to glue individual iTreg/tTreg/Tconv cells added to the disk. Only round and robust cells were selected for AFM gluing. The AFM cantilever carrying a single iTreg/tTreg/Tconv cell was lowered to allow T cell contact with an individual DC and to interact for 15 s before being moved upwards, until two cells were separated completely. The force curves were acquired. The process was then repeated. For each SCFS experiment, a pair of T–DC was used to generate force readings from each up and down cycle over a period of several minutes, a minimum of 14 force curves were collected in each T–DC pair. In all experiments, at least three such pairs were used for each condition. The force curves were further processed using the JPK image processing software.
RNA-seq
The RNA-seq library preparation was performed by a modified SMART-seq2 protocol (Habib et al., 2016; Picelli et al., 2014). Briefly, RNA was added into 4 µl elution mix made of 1 µl RT primer (10 µM), 1 µl dNTP mix (10 mM each), 1 µl RNase inhibitor (4 U/µl), and 1 µl H2O. Eluted samples were incubated at 72°C for 3 min and immediately placed on ice. Each sample was added with 7 µl reverse transcription (RT) mix made of 0.75 µl H2O, 0.1 µl Maxima RNase-minus RT (Thermo Fisher Scientific, #EP0741), 2 µl 5× Maxima RT buffer, 2 µl Betaine (5 M, Sigma-Aldrich, B0300), 0.9 µl MgCl2 (100 mM), 1 µl TSO primer (10 µM), and 0.25 µl RNase inhibitor (40 U/µl). The RT reaction was incubated at 42°C for 90 min and followed by 10 cycles of (50°C for 2 min, 42°C for 2 min), then heat inactivated at 70°C for 15 min. Samples were then amplified with an addition of 14 µl PCR mix made of 1 µl H2O, 0.5 µl ISPCR primer (10 µM), 12.5 µl KAPA HiFi HotStart ReadyMix (KAPA Biosystems, KK2602). The PCR reaction was performed as follows: 98°C for 3 min, 22 cycles of (98°C for 15 s, 67°C for 20 s, 72°C for 6 min), and final extension at 72°C for 5 min.
The amplified cDNA product was purified using VAHTS DNA Clean Beads. Sequencing libraries were prepared using TruePrep DNA Library Prep Kit V2 for Illumina (Vazyme, TD503-02).
Primer sequences: RT primer (Sangon), 5′Biotin-AAGCAGTGGTATCAACGCAGAGTACTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTVN; TSO primer (Sangon), 5′Biotin-AAGCAGTGGTATCAAC CAGAGTACAT/rG//rG//iXNA_G/; ISPCR primer (Sangon), 5′Biotin-AAGCAGTGGTATCAACGC AGA*G*T.
RNA-seq analysis
Raw fastq reads were trimmed by Cutadapt (version 1.18) to trim adapter and low-quality sequence.
Reads were aligned to the mouse genome (mm10) using STAR (version 2.5.3). The number of reads within each gene in each single cell was counted using RSEM (version 1.3.0) with gene annotation file from GENCODE (GRCm38.m23).
Differential expression was estimated by using DESeq2 package (version 1.26) with absolute log2 Fold change >0.5 and adjusted p-value <0.05.
Functional enrichment analysis was performed by using the clusterProfiler R package (v3.14.3) with default parameters (p-value <0.05) and the functional annotations terms in Gene Ontology (GO) (Ashburner et al., 2000).
ATAC-sequencing
ATAC-seq were performed using TruePrep DNA Library Prep Kit V2 for Illumina according to the manufacturer’s recommendation. Briefly, live cells were sorted using PI staining, then 5 × 104 cells were transferred to a low-binding centrifuge tube and resuspended in 50 μl of pre-chilled Lysis Buffer. After 20 times pipetting, cells were incubated on ice for 10 min for lysis. Subsequently, cells were pelleted by centrifugation at 4°C 500 × g for 5 min, and suspended in a mixture consisting of 10 μl of 5× TTBL, 5 μl of TTE Mix V50, and 35 μl of ddH2O. The mixture was gently mixed by pipetting 20 times, and incubated at 37°C for 30 min for fragmentation. The products were purified using 2× VAHTS DNA Clean Beads, eluted with 26 μl of water, and collected for the subsequent amplification. 24 μl of the eluted samples were amplified with different combination of i7 index and i5 index primers. QPCR were performed to determine the cycle number for PCR amplification.0.55× VAHTS DNA Clean Beads was employed for size selection of the PCR product. Finally, 30 ng of each sample were pooled for subsequent second-generation sequencing, on the Illumina NovaSeq 6000 with read length of 150 bp.
ATAC-seq analysis
Raw fastq reads were trimmed by Cutadapt (version 1.18) to trim adapter and low-quality sequence.
Reads were aligned to the mouse genome (mm10) using Bowtie2 (version 2.3.3.1), the parameters are -t -q -N 1 -L 25 -X 2000 --no-mixed --no-discordant.
Reads with alignment quality <Q30, mapped to chrM, overlapping with ENCODE blacklisted regions were discarded. Duplicates were removed using Picard (version 2.20.4) and Samtools (version 1.13). Open chromatin region peaks were called on using MACS2 peak caller (version 2.2.7.1) with the following parameters: -nomodel -nolambda -call-summits. Peaks from all samples were merged and peaks falls in ENCODE blacklisted regions were filtered out.
Differential peaks were estimated by using ‘edgeR’ method in the R package DiffBind (version 2.14.0) with absolute log2 Fold change >1 and adjusted p-value <0.05.
The peaks annotations were performed with the ‘annotatePeak’ function in the R package ChIPseeker (Yu et al., 2015).
The plot of ATAC-seq signals over a set of genomic regions were calculated by using ‘computeMatrix’ function in deepTools (version 2.0) and plotted by using ‘plotHeatmap’ and ‘plotProfile’ functions in deepTools.
Cut&Tag assay
Foxp3-RFP+iTreg were induced, and infected with NFAT-Flag-GFP lentiviruses on days 1.5 and 2.5 of induction. On day 5, Foxp3-RFP+MSCV-NFAT-Flag-GFP+iTregs were sorted, then rested and restimulated for 2 days. Cut&Tag assay was performed using Hyperactive Universal CUT&Tag Assay Kit for Illumina (Vazyme Biotech Co., Ltd) according to the manufacturer’s recommendation. Briefly, 1 × 105 iTregs were collected to low-binding centrifuge tubes and washed. Nucleus isolation were performed by NE buffer for 10 min on ice, and crosslinked by 0.1% formaldehyde for 2 min, then stopped by glycine. Con A Beads were pre-treated with binding buffer, and then incubated with nuclei of iTreg for 10 min. Further, samples were incubated with anti-Flag M2 antibody (1:50) or without primary antibody (Mock) overnight at 4°C. Goat Anti-Mouse IgG secondary antibody (1:50) was added and incubated for 1 hr. Cells were washed three times, and then incubated with pA/G-Tnp enzyme for 1 hr. After incubation, samples were mixed with 5× TTBL and Dig-300 Buffer and incubated at 37°C for 60 min for fragmentation. DNA was extracted and purified, then amplified with different combination of i7 index and i5 index primers. QPCR was performed to determine the cycle number for PCR amplification.VAHTS DNA Clean Beads were employed for size selection of the PCR product. Finally, 30 ng of each sample were pooled for subsequent second-generation sequencing, on the Illumina NovaSeq 6000 with read length of 150 bp.
Cut&Tag data analysis
Raw fastq reads were trimmed by Cutadapt (version 1.18) to trim adapter and low-quality sequence.
Reads were aligned to the mouse genome (mm10) using Bowtie2 (version 2.3.3.1), the parameters are -t -q -N 1 -L 25 -X 2000 --no-mixed --no-discordant.
Reads with alignment quality <Q30, mapped to chrM, overlapping with ENCODE blacklisted regions (https://sites.google.com/site/anshulkundaje/projects/blacklists) were discarded. Duplicates were removed using Picard (version 2.20.4) and Samtools (version 1.13). Open chromatin region peaks were called using MACS2 peak caller (version 2.2.7.1) with the following parameters: -nomodel -nolambda -call-summits. Peaks from all samples were merged and peaks falls in ENCODE blacklisted regions were filtered out.
The peaks annotations were performed with the ‘annotatePeak’ function in the R package ChIPseeker (Yu et al., 2015).
The plot of Cut&Tag signals over a set of genomic regions were calculated by using ‘computeMatrix’ function in deepTools and plotted by using ‘plotHeatmap’ and ‘plotProfile’ functions in deepTools.
The motif enrichment analysis was performed by using the ‘findMotifsGenome.pl’ command in HOMER with default parameters.
The motif occurrences in each peak were identified by using FIMO (MEME suite v5.0.4) with the following settings: a first-order Markov background model, a p-value cutoff of 10−4, and PWMs (Position weight matrix) from the mouse HOCOMOCO motif database (v11).
Statistical analysis
Calculation and graphing were done with Prism (GraphPad). Data were presented as mean ± standard error of the mean, unless indicated otherwise. Student’s t-test was used for comparing endpoint means of different groups. p-value <0.05 was considered significant. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; N.S., not significant. Numbers of independently performed experiment repeats are shown as N, biological replicates of each experiment as n.
Appendix 1
Data availability
RNA/ATAC/CUT&TAG associated sequencing data have been deposited in GEO under accession code GSE246431. All data excepted from the sequencing data, generated or analyzed during this study are included in the Source Data files.
-
NCBI Gene Expression OmnibusID GSE246431. Sustained store-operated calcium entry utilizing activated chromatin state leads to instability in iTregs.
References
-
In vivo instruction of suppressor commitment in naive T cellsThe Journal of Experimental Medicine 199:1401–1408.https://doi.org/10.1084/jem.20040249
-
Gene Ontology: tool for the unification of biologyNature Genetics 25:25–29.https://doi.org/10.1038/75556
-
Conversion of peripheral CD4+CD25- naive T cells to CD4+CD25+ regulatory T cells by TGF-beta induction of transcription factor Foxp3The Journal of Experimental Medicine 198:1875–1886.https://doi.org/10.1084/jem.20030152
-
IL-2 Controls the Stability of Foxp3 Expression in TGF-β–Induced Foxp3+ T Cells In VivoThe Journal of Immunology 186:6329–6337.https://doi.org/10.4049/jimmunol.1100061
-
Strong adhesion by regulatory T cells induces dendritic cell cytoskeletal polarization and contact-dependent lethargyThe Journal of Experimental Medicine 214:327–338.https://doi.org/10.1084/jem.20160620
-
The functional contribution of calcium ion flux heterogeneity in T cellsImmunology and Cell Biology 93:694–704.https://doi.org/10.1038/icb.2015.34
-
CD25- T cells generate CD25+Foxp3+ regulatory T cells by peripheral expansionJournal of Immunology (Baltimore, Md 173:7259–7268.https://doi.org/10.4049/jimmunol.173.12.7259
-
Developmental regulation of Foxp3 expression during ontogenyThe Journal of Experimental Medicine 202:901–906.https://doi.org/10.1084/jem.20050784
-
TGF-beta induces Foxp3 + T-regulatory cells from CD4 + CD25 - precursorsAmerican Journal of Transplantation 4:1614–1627.https://doi.org/10.1111/j.1600-6143.2004.00566.x
-
Ca2+ signaling in T-lymphocytesCritical Reviews in Immunology 18:419–448.https://doi.org/10.1615/critrevimmunol.v18.i5.20
-
Alloantigen-specific de novo-induced Foxp3+ Treg revert in vivo and do not protect from experimental GVHDEuropean Journal of Immunology 39:3091–3096.https://doi.org/10.1002/eji.200939432
-
Orai1 subunit stoichiometry of the mammalian CRAC channel pore: CRAC channel subunit stoichiometryThe Journal of Physiology 586:419–425.https://doi.org/10.1113/jphysiol.2007.147249
-
Sumoylation of the transcription factor NFATc1 leads to its subnuclear relocalization and interleukin-2 repression by histone deacetylaseThe Journal of Biological Chemistry 284:10935–10946.https://doi.org/10.1074/jbc.M900465200
-
Full-length RNA-seq from single cells using Smart-seq2Nature Protocols 9:171–181.https://doi.org/10.1038/nprot.2014.006
-
Regulatory T cells and Foxp3Immunological Reviews 241:260–268.https://doi.org/10.1111/j.1600-065X.2011.01018.x
-
A kinetic and dynamic analysis of Foxp3 induced in T cells by TGF-betaJournal of Immunology 179:7667.https://doi.org/10.4049/jimmunol.178.12.7667
-
NFATc1 autoregulation: a crucial step for cell-fate determinationTrends in Immunology 27:461–469.https://doi.org/10.1016/j.it.2006.08.005
-
tTregs, pTregs, and iTregs: similarities and differencesImmunological Reviews 259:88–102.https://doi.org/10.1111/imr.12160
-
Regulation of the NFAT pathway discriminates CD4+CD25+ regulatory T cells from CD4+CD25- helper T cellsJournal of Leukocyte Biology 83:708–717.https://doi.org/10.1189/jlb.0507321
-
NFAT functions as a working memory of Ca2+ signals in decoding Ca2+ oscillationThe EMBO Journal 22:3825–3832.https://doi.org/10.1093/emboj/cdg381
-
Smad3 and NFAT cooperate to induce Foxp3 expression through its enhancerNature Immunology 9:194–202.https://doi.org/10.1038/ni1549
-
Kinetics and extent of T cell activation as measured with the calcium signalThe Journal of Experimental Medicine 185:1815–1825.https://doi.org/10.1084/jem.185.10.1815
-
Control of Foxp3 stability through modulation of TET activityThe Journal of Experimental Medicine 213:377–397.https://doi.org/10.1084/jem.20151438
-
Induced regulatory T cells promote tolerance when stabilized by rapamycin and IL-2 In VivoThe Journal of Immunology 191:5291–5303.https://doi.org/10.4049/jimmunol.1301181
-
Regulation of the murine Nfatc1 gene by NFATc2The Journal of Biological Chemistry 277:10704–10711.https://doi.org/10.1074/jbc.M107068200
Article and author information
Author details
Funding
National Natural Science Foundation of China (31370878)
- Yan Shi
National Natural Science Foundation of China (31630023)
- Yan Shi
National Natural Science Foundation of China (81621002)
- Yan Shi
Canadian Institutes of Health Research (PJT-166155)
- Yan Shi
Canadian Institutes of Health Research (PJT-156334)
- Yan Shi
Natural Sciences and Engineering Research Council of Canada (RGPIN/03748-2018)
- Yan Shi
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Hai Qi for providing Foxp3-GFP mice, platE cell line and guidance of calcium-related experimental design. We thank Dr. Saiyong Zhu for providing TetO-FUW-flag-Nfatc1 plasmid. We also thank Imaging Core Facility, Technology Center for Protein Sciences of Tsinghua University for assistance. Y.S. is supported by the joint Peking-Tsinghua Center for Life Sciences, the National Natural Science Foundation of China General Program (31370878), State Key Program (31630023), and Innovative Research Group Program (81621002), CIHR (PJT-156334 and PJT-166155) and NSERC (RGPIN/03748-2018).
Ethics
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved Institutional Animal Care and Use Committee (IACUC) protocols of the Tsinghua University. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Tsinghua University (22-SY6). All surgery was performed under sodium pentobarbital anesthesia, and every effort was made to minimize suffering.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.88874. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Lyu, Yuan et al.
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
-
- 556
- views
-
- 107
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Genetics and Genomics
Solitary fibrous tumors have gene expression signatures similar to those of neuroendocrine tumors.
-
- Genetics and Genomics
Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling due to steric similarities to endogenous hormones. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) and differentially expressed genes (DEGs) that predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many EDCs persists over multiple generations. Many studies have described direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDC-induced epimutations are introduced or subsequently propagated, whether there are cell type-specific susceptibilities to the same EDC, and whether this correlates with differential expression of relevant hormone receptors. We exposed cultured pluripotent (iPS), somatic (Sertoli and granulosa), and primordial germ cell-like (PGCLC) cells to BPS and found that differential incidences of BPS-induced epimutations and DEGs correlated with differential expression of relevant hormone receptors inducing epimutations near relevant hormone response elements in somatic and pluripotent, but not germ cell types. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and DEGs was largely retained, however, >90% of the specific epimutations and DEGs were replaced by novel epimutations and DEGs. These results suggest a unique mechanism by which an EDC-induced epimutated state may be propagated transgenerationally.