Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Research Article
Transparent Process
Open Access

Dynamic changes in the regulatory T-cell heterogeneity and function by murine IL-2 mutein

View ORCID ProfileDaniel R Lu, Hao Wu, View ORCID ProfileIan Driver, Sarah Ingersoll, View ORCID ProfileSue Sohn, Songli Wang, View ORCID ProfileChi-Ming Li  Correspondence email, View ORCID ProfileHyewon Phee  Correspondence email
Daniel R Lu
2Genome Analysis Unit, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Daniel R Lu
Hao Wu
1Department of Oncology and Inflammation, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ian Driver
2Genome Analysis Unit, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ian Driver
Sarah Ingersoll
1Department of Oncology and Inflammation, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sue Sohn
1Department of Oncology and Inflammation, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sue Sohn
Songli Wang
2Genome Analysis Unit, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chi-Ming Li
2Genome Analysis Unit, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Chi-Ming Li
  • For correspondence: chimingl@amgen.com
Hyewon Phee
1Department of Oncology and Inflammation, Amgen Research, Amgen Inc, South San Francisco, CA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Hyewon Phee
  • For correspondence: hyewonp@amgen.com
Published 8 April 2020. DOI: 10.26508/lsa.201900520
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

The therapeutic expansion of Foxp3+ regulatory T cells (Tregs) shows promise for treating autoimmune and inflammatory disorders. Yet, how this treatment affects the heterogeneity and function of Tregs is not clear. Using single-cell RNA-seq analysis, we characterized 31,908 Tregs from the mice treated with a half-life extended mutant form of murine IL-2 (IL-2 mutein, IL-2M) that preferentially expanded Tregs, or mouse IgG Fc as a control. Cell clustering analysis revealed that IL-2M specifically expands multiple sub-states of Tregs with distinct expression profiles. TCR profiling with single-cell analysis uncovered Treg migration across tissues and transcriptional changes between clonally related Tregs after IL-2M treatment. Finally, we identified IL-2M–expanded Tnfrsf9+Il1rl1+ Tregs with superior suppressive function, highlighting the potential of IL-2M to expand highly suppressive Foxp3+ Tregs.

Introduction

Foxp3+ regulatory T cells (Tregs) play a fundamental role in immunosuppression and immune tolerance, and there is great interest in harnessing Treg populations to treat autoimmune and inflammatory disorders. The differential expression of transcription factors, costimulatory receptors, chemokine receptors, and secreted effectors in quiescent and activated Tregs suggests that the heterogeneous Treg states exist and perform distinct functions (Zheng et al, 2006; Menning et al, 2007; Schiering et al, 2014). Moreover, nonlymphoid tissue Tregs acquire unique phenotypes different from lymphoid-tissue Tregs, suggesting that the anatomical location of Tregs contributes to their heterogeneity (Sather et al, 2007; Miragaia et al, 2019).

Recently, low-dose Interleukin-2 (IL-2) therapies have been tested to induce tolerance in patients with autoimmunity and inflammatory disorders (Koreth et al, 2011; Saadoun et al, 2011; Hartemann et al, 2013; Matsuoka et al, 2013; Klatzmann & Abbas, 2015; Yu et al, 2015). Although the low-dose IL-2 therapies expand Tregs, their effect has been limited by concomitant increases in conventional effector T cells and natural killer cells. To improve selectivity and pharmacokinetics of low-dose IL-2, alternative modalities have been considered (Peterson et al, 2018). However, it is not clear how IL-2–based therapies impact Treg heterogeneity in diverse tissues. Because the goal of Treg-targeted therapies is to expand Treg-mediated tolerance at proper anatomical locations, it is critical to understand how IL-2–mediated expansion impacts the phenotypic and functional heterogeneity of Tregs in lymphoid and nonlymphoid tissues.

Thymic-derived Foxp3+ Tregs undergo TCR-dependent antigen priming and activation-induced expansion in lymphoid organs followed by extravasation into peripheral tissues, where they acquire tissue-specific tolerogenic phenotypes. Given the complex migratory patterns of Tregs, it is unclear how IL-2–mediated therapy affects Tregs within and across tissues. TCR sequencing combined with single-cell profiling provides an opportunity to measure IL-2–induced Treg differentiation and movement by tracing the transcriptional conversions and trafficking patterns of clonal lineages.

To better understand the impact of the IL-2–mediated Treg expansion therapy on Foxp3+ Treg heterogeneity in lymphoid and nonlymphoid tissues, we profiled mouse spleen, lung, and gut Tregs using single-cell RNA-seq (scRNA-seq) with TCR sequencing under murine IL-2 mutein (IL-2M) stimulation or homeostatic (mouse IgG Fc isotype control–treated) conditions. Comparison of resting, primed/activated, and activated Treg states from different tissues revealed unique gene signatures shared between spleen and lung Tregs, as well as distinct activation profiles of gut Tregs. Administration of murine IL-2M dramatically changed the landscape of Tregs in the spleen and the lung, although maintaining tissue-specific identity in the gut. TCR profiling coupled with scRNA-seq revealed gene expression dynamics governing Treg differentiation after IL-2M treatment and uncovered a migratory axis across tissues. In addition, we identified a population of activated Tnfrsf9+Il1rl1+ Tregs in mice that expands after IL-2M and suppresses convention T cells robustly in vitro. Overall, our experiments provide new insights into the relationships between Foxp3+ Treg activation states and their phenotypic heterogeneity in different tissues during homeostasis and after murine IL-2M stimulation.

Results

A half-life–extended mutant form of murine IL2 expands CD25+Foxp3+ Tregs in vivo

To determine the specific role of mouse IL-2 in Foxp3+ Tregs in mice, a half-life–extended mutant form of murine IL-2 (IL-2 mutein, IL-2M) was generated as a mouse IgG2a Fc fusion protein (Fig S1A). Previously, a human form of long-lived IL2 mutein (human IgG-(human IL-2N88D)2) was reported (Peterson et al, 2018). In this human IL-2 mutein, an effector-silent human IgG1 was fused to a mutant form of human IL2 to increase the half-life. Moreover, the N88D mutation was introduced to human IL2 to decrease its binding to the intermediate affinity IL2 receptor, IL2βγ, whereas maintaining its binding to the high-affinity IL2 receptor, IL2αβγ. For the mouse IL-2 mutein, an effector-silent mouse IgG2a Fc (N297G) (Tao & Morrison, 1989) was fused to a mutant form of IL2 to increase the half-life. Furthermore, D34S and N103D mutations were introduced to the mouse IL2 because both amino acids were described to be critical for IL2’s binding to IL2Rβ, whereas minimally affecting interaction with IL2Rα (Zurawski & Zurawski, 1989; Zurawski et al, 1993). The N103 residue of mouse IL2 corresponds to the N88 of human IL2. In addition to D34S and N103D mutations, two additional mutations (C140A and P51T) were introduced to mouse IL2 to facilitate manufacturability (Klein et al, 2017).

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Characterization of murine IL-2 mutein (IL-2M).

(A) Schematic diagram of murine IL-2M. (B) Intracellular staining of phosphorylated STAT5 from CD25+Foxp3+ CD4 Tregs, CD25−Foxp3− CD4 Tconv, and CD25+Foxp3− CD4 Teff cells from mouse spleens stimulated by increasing the dose of wild-type recombinant mouse IL-2 or IL-2M in vitro. Data are representative of two independent experiments. (C) Gating strategy of CD25+Foxp3−, CD25+Foxp3+, CD25−Foxp3+, and CD25−Foxp3− populations from CD4 T cells. (D) Fold increase of cell numbers shown in Fig 1C. Results are representative of at least two independent experiment using three mice per each group.

Intracellular staining of phosphorylated STAT5 as a response to IL2 signaling was performed to determine activity of mouse IL-2M in vitro (Fig S1B and C). As expected, IL-2M induced phopsho-STAT5 in CD25+Foxp3+ Tregs that express the high-affinity IL2Rαβγ in a dose-dependent manner, whereas it did not induce phospho-STAT5 in CD25-Foxp3− Tconv cells that express the intermediate IL2R, IL2Rβγ (Fig S1B, top panel). In contrast, wild-type recombinant mouse IL2 (rmIL2) induced phospho-STAT5 in CD25-Foxp3− Tconv cells as well as CD25+Foxp3+ Tregs (Fig S1B, bottom left panel). Because the mutations introduced were intended to reduce the interaction of mouse IL2 with IL2Rβ but maintain interaction with IL2Rα (CD25) to allow IL-2M’s binding to the high-affinity IL2R (ILRαβγ) on Tregs, the in vitro activity of IL-2M was attenuated compared with wild-type rmIL2. Furthermore, IL-2M slighted activated phospho-STAT5 in CD25+Foxp3− effector Tconv cells, although the degree of activation was markedly reduced compared with wild-type rmIL2 (Fig S1B, bottom right panel).

Administration of murine IL-2M to C57BL/6 mice with different doses revealed that 0.33 mg/kg of murine IL-2M specifically increased Foxp3+ Tregs in the spleen, lymph nodes, blood, and lung, whereas it did not affect CD25−Foxp3− T conventional cells (Tconv) (Figs 1A–D and S1D). Murine IL-2M increased percentages and cell numbers of CD25+Foxp3− activated CD4 T cells slightly at a higher dose (Fig 1A), likely because of the presence of high-affinity IL2R (IL2Rαβγ) in CD25-expressing cells (Malek, 2008; Letourneau et al, 2009). Expansion of Foxp3+ Tregs by IL-2M was comparable with IL-2/anti-IL2 antibody (JES6-1) conjugate (IL-2C), which was previously reported to expand Tregs (Boyman et al, 2006) (Fig S2A–D).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. A half-life–extended mutant form of the murine IL-2M expands CD25+Foxp3+ Tregs in vivo.

(A) Dose-dependent effects of murine IL-2M on the percentages of CD25−Foxp3+ and CD25+Foxp3+ Tregs and CD25+Foxp3− activated Tconv and CD25−Foxp3− Tconv within CD4 T cells in the blood, spleen, lymph nodes, and lung. Experiments were performed using three mice per each group. Control was untreated mice. (B) Effects of murine IL-2M (0.33 mg/kg) on the percentages of each population within CD4 T cells from the blood, spleen, lymph nodes, and lung. (C) Effects of murine IL-2M (0.33 mg/kg) on the cell numbers of each population from the blood, spleen, lymph nodes, and lung. (D) Effects of murine IL-2M (0.33 mg/kg) on the proliferation of each cell population assessed by expression of Ki67 (mean fluorescence intensity). (B, C, D) results are representative of at least two independent experiment. Experiments were performed using three mice per each group. Control: PBS-treated mice. Statistics: two-way ANOVA for multiple comparisons. *0.01 < P < 0.05, **0.001 < P < 0.01, ***0.0001 < P < 0.001, ****P < 0.00001.

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Comparing murine IL-2M and IL-2/anti-IL-2 antibody (JES6-1) complex (IL-2C).

(A, B, C, D) Effects of murine IL-2M and IL-2/anti-IL-2 antibody (JES6-1) complex (IL-2C) on the percentages and cell numbers of CD4 (A); CD8 (B); CD4+CD25+Foxp3+ (C); and CD4+CD25+Foxp3− T cells (D). Results shown are representative of two experiments, each using three mice for conditions. Statistics: one-way ANOVA for multiple comparisons. ***0.0001 < P < 0.001, **0.001 < P < 0.01, *0.01 < P < 0.05.

In vivo–expanded Tregs by IL-2M displayed superiority in suppressing the proliferation of Th1 effector cells (Fig 2A) and generation of IFNγ from Th1 effector cells in vitro compared with Foxp3+ Tregs isolated from spleens of PBS-treated mice (Fig 2B). Although these results revealed that highly suppressive Tregs expanded after IL-2M treatment in vivo, it was not clear how IL-2M impacted on the heterogeneity of Tregs toward functional suppression in lymphoid and nonlymphoid tissues. We, therefore, sought to profile the molecular phenotypes expressed by Tregs to understand how IL-2M elevates Treg expansion and suppression in mice.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. Superior suppressive activity of CD25+Foxp3+ Tregs expanded by murine IL-2M.

(A) In vitro suppression assays using sorted Foxp3-EGFP+ Tregs from the spleen of either PBS or murine IL-2M–treated Foxp3-EGFP mice. Naïve T cells were MACS sorted using Miltenyi naïve T-cell isolation kit, labeled with cell tracer violet (CTV), and then mixed with irradiated APC in the absence (No Tregs) or presence of Tregs from PBS- or IL-2M–treated mice. Tregs were cocultured then with naïve T cells and APCs under Th1 skewing condition. Shown are representative images of proliferation of Tconv indicated as dilution of CTV (left panel). Suppression of proliferation of Tconv was shown by diluted CTV MFI (middle panel) or percent divided (right panel). (B) Suppression of IFNγ generation from Th1 effector cells by Tregs. (B) In vitro Treg suppression assay was performed as described in (B), and expression of IFNγ from Tconv was measured by intracellular staining. Results are representative of two independent experiments. Statistics: one-way ANOVA for multiple comparisons. **0.001 < P < 0.01, ***0.0001 < P < 0.001, ****P < 0.00001.

Isolation and scRNA-seq of Tregs from the spleen, lung, and gut

The therapeutic success of Treg expansion depends on the efficacy of Treg immunosuppression at specific anatomical sites. To characterize the heterogeneity of Tregs and assess how IL-2M reshapes this diversity, we performed scRNA-seq on Foxp3+ Tregs from multiple tissues at steady state and after IL-2M treatment. CD4+eGFP+ single cells were sorted from the spleen, lung, and gut of Foxp3-eGFP mice and single cell libraries were prepared using the 10× Chromium platform (Fig 3A and Table S1). After filtering and cross-sample normalization using Seurat (Butler et al, 2018), we recovered 17,097 spleen Tregs, 10,353 lung Tregs, and 4,458 gut Tregs across three replicates (except gut Tregs treated with IL-2M [n = 2]) with roughly equivalent Tregs in mouse IgG Fc isotype control (Iso)–treated and IL-2M–treated (IL-2M) conditions (16,152 and 15,756 cells, respectively). Single-cell molecule and gene recovery rates showed good concordance between Iso- and IL-2M–treated cells (Fig S3A and B). To establish that we had accurately sorted and sequenced Tregs, we additionally performed scRNA-seq on CD4 conventional T cells sorted using FACS. Comparison of Tregs and Tconvs confirmed that Tregs and Tconvs expressed distinct transcriptional profiles and clustered into distinct groups (Fig S4A). In addition, Tregs expressed higher transcript levels of established Treg genes such as Foxp3, Il2ra, Ctla4, Ikzf2, and Nrp1, whereas both cell types expressed similar levels of Cd4 (Fig S4B and C).

Table S1 CD4+ regulatory and conventional T-cell single-cell RNA-seq samples across treatments and tissues.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Classification of heterogeneous Treg resting, primed, and activation states across tissues.

(A) Experimental workflow for interrogation of Tregs in the spleen, lung, and gut after administration of IL-2M or mouse IgG Fc control (Iso) by scRNA-seq. (B) t-SNE projection of all Tregs grouped by cell subset (top) or by cell state (bottom) using K-nearest neighbor graph–based clustering. Individual cells are colored by cell state cluster assignments. (C) Heat map depicting the mean expression (Z-score) of most representative signature genes (rows) for each cluster (columns). Genes shown in blue text with an asterisk are genes used to characterize Treg states. (D) Violin plots depicting gradual down-regulation of classical resting Treg markers (left) and gradual up-regulation of classical activated Treg markers (right) in progressively higher Treg activation states. (E) Violin plots showing marker genes associated with sub-populations of activated Tregs. Statistics: Wilcoxon rank sum test. **P < 0.01, ***P < 0.001, n.s. = P > 0.05.

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. Single-cell RNA-seq (scRNA-seq) variation across samples and treatment.

(A) Number of genes and unique mRNA molecules for all single-cell RNA-seq libraries after isotype control or IL-2M treatment. Single cell libraries are colored by sample. (B) Histogram of the distribution of unique mRNA molecules across all single-cell RNA-seq libraries. Blue dotted line indicates the mean number of molecules per cell for isotype control– and IL-2M–treated cells.

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Treg expression of established Treg genes relative to Tconv cells.

(A) t-SNE projection of single regulatory T cells (Treg) and conventional CD4 T cells (Tconv). (B, C) t-SNE projections (B) and violin plots (C) show the single-cell expression levels of established Treg marker genes Foxp3, Il2ra, Ctla4, Ikzf2, and Nrp1, as well as even Cd4 expression between Treg and Tconv cells. n.s. = P > 0.05 (not significant), *** = P < 0.001.

Molecular characterization of resting, primed, and activated Treg cell states across all tissues

Previous scRNA-seq profiling efforts have identified resting, early-activated (primed/activated), and terminal effector subsets along the Treg cell differentiation continuum (Cheng et al, 2012; Gaublomme et al, 2015; Guo et al, 2018; Zemmour et al, 2018; Miragaia et al, 2019), with each subset containing multiple transcriptionally defined cell states. We sought to organize all Tregs from our dataset into different cellular states to generate a single-cell Treg classification schema, which could then define cellular variation across tissues and after IL-2M stimulation at a detailed resolution. To do this, we first impartially defined cell states using unsupervised clustering, and then identified highly representative signature genes from each cell state to anchor our cell state assignments in the context of known Treg biology.

Cell clustering of all Tregs using the graph-based Louvain algorithm from Seurat (Satija et al, 2015) identified three established Treg subsets (Fig 3B, top), which could be divided further into 10 Treg cell states of activation (Fig 3B, bottom). These clusters were evenly represented in all replicates, indicating that they represent biological heterogeneity and not batch-specific processing artifacts (Fig S5A and B).

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. Variation of single-cell RNA-seq cell state classification across experimental replicates.

(A) t-SNE projection consisting of the cells for each experimental replicate (also known as one isotype-treated mouse and one IL2-mutein treated mouse) from which Tregs were isolated and profiled using single-cell RNA-seq. (B) Frequency of Treg single cell libraries per experimental replicate assigned to each to Treg cell state classification. Statistical testing for batch effects on the observed counts of Tregs in each cluster was performed for each pair of batches using Fisher’s exact test.

Signature genes were then defined for each cell state to understand transcriptional differences between Treg populations (see the Materials and Methods section) (Figs 3C and S6A). When grouped by resting, primed/activated, and activated Treg subsets, cell states had differential expression of known lymphoid-associated and activated genes. Resting Tregs (C1) have high expression of lymphoid-tissue homing receptors (Ccr7 and S1pr1) and Treg-establishing transcriptional regulators (Lef1 [Xing et al, 2019], Satb1 [Kitagawa et al, 2017], and Klf2 [Pabbisetty et al, 2016]) and low-activation marker expression (Ctla4, Icos, Tnfrsf1b, and Tnfrsf4) (Figs 3D and S7A), suggesting that they either have not yet undergone TCR-mediated activation or have a central memory-like phenotype.

Figure S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S6. Distinct gene expression profiles between Treg states.

(A) Heat map showing top 5–7 overlapping genes for each Treg cell state after differential expression using Wilcoxon rank sum test. Vertical lines depict individual cells belonging to each cell state. (B, C, D, E) Select heat maps showing the top 10 most differentially expressed genes from all Tregs (spleen, lung, and gut) of both isotype control– and IL-2M–treated mice when comparing between two Treg cell states using MAST. Genes are colored by the difference in log2-fold change expression between cell states.

Figure S7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S7. Cellular heterogeneity within primed-early–activated Treg sub-populations.

(A) Violin plots showing intermediate C2 and C9 expression of classical resting and activated Treg marker genes. (B) Volcano plot showing gene expression profile between C1-resting and C2-primed/activated Treg states. Labeled genes show most significant genes based on P-value and fold change thresholds. (C) Volcano plot showing gene expression profile between C1 resting and C9 primed/activated Treg states. Labeled genes show most significant genes based on P-value and fold-change thresholds. (D) Violin plot showing high Ms4a4b expression primed-early–activated Treg states. Clusters with significantly higher gene expression relative to all other cell states are shown with asterisks. ** = P < 0.01, *** = P < 0.001.

Two primed/activated Treg states (C2 and C9) express intermediate levels of both resting and activation gene sets, suggesting they may be transitioning toward an activated phenotype (Figs 3D and S7A). In support of this observation, primed/activated Tregs are distinguishable from C1 resting Tregs with higher expression of Treg-specific genes that stabilize inhibitory activity of Treg (C2: Ikzf2 [Kim et al, 2015] and S100 family proteins) or inflammatory response mediators (C9: Stat1, Cxcl10, and interferon-response genes) (Fig S7B and C). Both C2 and C9 have high Ms4a4b expression (Fig S7D), which modulates the activation by regulating proliferation and augmenting co-stimulation through GITR (Howie et al, 2009; Xu et al, 2010). Furthermore, C2- and C9-Tregs could be distinguished from each other, as C2-Tregs express more Nrp1 (which is expressed highly in natural Tregs [Yadav et al, 2012] and promotes Treg survival by interacting with Sema4a on Tconv cells [Delgoffe et al, 2013]), whereas C9-Tregs uniquely express interferon-response genes (Epeldegui et al, 2015) (Bst2, Tgtp2, Ifit1, and Isg15) (Fig S8A and B).

Figure S8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S8. Distinct interferon and Treg stability signatures differentiate primed/activated Treg sub-populations.

(A) Violin plots depicting selected signature genes differentially expressed between C2 and C9 primed–activated Treg states. Clusters with significantly higher gene expression relative to all other cell states are shown with asterisks. P-values showing differential expression between C2 and C9 were calculated using MAST. (B) Heat map showing the top 10 most differentially expressed genes between C2 and C9 Tregs using MAST, colored by log2-fold change expression. ** = P < 0.01, *** = P < 0.001.

Six Treg clusters (C3-C5, C7, C8, and C10) showed an activated phenotype with the lowest expression of lymphoid-tissue homing markers and highest expression of activation genes (Fig 3D). Whereas activated cell states express high Tnfrsf9 (encoding for costimulatory activator 4-1BB), each population also expresses distinct genes (Fig 3E). Differential gene expression of these activated clusters compared with C1-resting cluster or C2-primed/activated cluster or each other was shown (Figs S6B–E, S9A–D, and S10A–C). Both C4 and C8 Tregs express high Il1rl1 (ST2) compared with all other states (Liao et al, 2010; Bandala-Sanchez et al, 2013; Mahmud et al, 2014; Schiering et al, 2014; Geiger et al, 2016; Siede et al, 2016), whereas they differentially express T-cell activation regulators and effector molecules (C4: Ass1, Klrg1, Gzmb, and Tagln2 versus C8: Itgb7, Cd52, and GITR) (Figs 3E, S6B–E, S9B, and S10B). Similarly, C4 and C8 Tregs express high CXCR3, a chemokine receptor regulated by T-bet and associated with tissue Tregs that suppress Th1 responses (Koch et al, 2009) (Fig 3E). On the other hand, C3 Tregs express genes required for follicular regulatory phenotype (C3: Maf [Wheaton et al, 2017]) and genes that promote Treg survival and persistence (C3: Nr4a2, Cst7) (Figs 3E, 4C, S6B, and S9A), whereas C5 Tregs express effectors that promotes suppression (C5: Il10 [Chaudhry et al, 2011], Gzma [Velaga et al, 2015], and Gzmb [Cao et al, 2007; Loebbermann et al, 2012]) or tissue protection (C5: Areg [Arpaia et al, 2015]) (Figs 3E, 4C, S6B, and S9C). In addition, we identified one proliferative cell state (C6) (Mki67 and Top2a) that co-expresses activation markers (Fig S9D). C7 Tregs are a small population that expresses stress response genes (Hspa1a and Hspa1b) (Fig S10A). C10 Tregs is a small population as well, expressing genes associated with NK cells (NKG7, CCl5, CXCR6, and Gzmb) and cytokine (Ifng), which is expressed in ex-Tregs (Daniel et al, 2015) (Fig S10C).

Figure S9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S9. Different gene expression profiles between C1 Tregs and C3, C4, C5, and C6 Tregs.

(A, B, C, D) Volcano plots showing gene expression profiles between C1 Tregs and C3, C4, C5, and C6 Tregs. Labeled genes show most significant genes based on P-value and fold change thresholds.

Figure S10.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S10. Different gene expression profiles between C1 Tregs and C7, C8, and C10 Tregs.

(A, B, C) Volcano plots showing gene expression profiles between C1 Tregs and C7, C8, and C10 Tregs. Labeled genes show most significant genes based on P-value and fold change thresholds.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Heterogeneous landscape of Treg cell states defined between the spleen, lung, and gut.

(A) t-SNE projection of Treg cell states in the spleen, lung, gut, and combined at steady state (isotype control–treated). Individual cells are colored by Treg cell state classification from Fig 3. (B) Comparison of Treg signature genes between all tissue-specific cells in the spleen versus lung, spleen versus gut, and lung versus gut at steady state. Significant genes are shown in red (adjusted P-value < 0.05, MAST). (C, D, E, F) Select heat maps showing the top 10 most differentially expressed genes from all Tregs (spleen, lung, and gut) of isotype control (Iso)–treated mice when comparing between two Treg cell states using MAST. Genes are colored by the difference in log2-fold change expression between cell states.

Analysis of Treg clusters at steady-state identified tissue adaptations of Tregs

We next sought to characterize the Treg landscape in the spleen, lung, and gut from isotype control–treated mice using the classification of Treg states described in Fig 3 as a reference criterion. At the tissue level, the spleen and lung share a >40% frequency of C1-resting and minor frequencies of primed/activated and activated Tregs. Conversely, >80% of gut Tregs are activated and the majority are C5-Tregs (Fig 5B–D). Although lymphoid organs, such as the spleen, are known to contain large reservoirs of resting Tregs that can be recruited during immune challenge, the prominence of resting Tregs in the lung is noteworthy and suggests that the lung plays an active role in steady-state immune surveillance.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. IL-2 mutein induces a convergent expansion toward distinct Treg states in the spleen, lung, and gut.

(A) t-SNE projection of all Tregs shown by cell state from Iso-treated mice (left) or Tregs from IL-2M–treated mice (right) using K-nearest neighbor graph–based clustering. Individual cells are colored by cell state cluster assignments. (B) t-SNE projection of Treg cell states in the spleen, lung, and gut from isotype control–treated (Iso) and IL-2 mutein treated mice. Treg clusters that are enriched at least twofold in isotype control–treated or in IL-2 mutein–treated mice relative to the other treatment condition are outlined in with a black dotted line. Prominent populations (>10% of total Tregs in a given tissue) whose frequencies are consistent across treatment are highlighted in gray. (C) The tissue-specific frequency of Treg cell states for each experimental replicate. (D) Bar plots showing shift in Treg population in the spleen, lung, and gut. N = 3 for all isotype and IL-2 mutein treatments except for the IL-2 mutein–treated gut Tregs (N = 2). Statistics: Welch’s t test. *P < 0.05, **P < 0.01, ***P < 0.001.

Activated states of Tregs differ substantially between the tissues, indicative of tissue adaptations acquired by Tregs (Fig 4A and B). The spleen in the unperturbed state uniquely contained more than 30% C3-Nr4a2hi–activated Tregs, which express differentiation-promoting transcription factor Nr4a2 (Sekiya et al, 2011). In contrast, the lung did not contain the C3-Nr4a2hi Tregs as a main activation state. Rather, the lung contains the C4-Ass1hi and C8-Itgb7hi–activated clusters and C6-Stmn1hi proliferative cluster. Comparison of C3 and C4 revealed that C3 cluster differentially express TBC1D4, a transcription factor highly expressed in follicular regulatory Tregs, suggesting C3-Nr4a2+–activated Tregs include Tfr in the spleen (Fig 4E). The co-expression of immunomodulatory genes (Izumo1r [Yamaguchi et al, 2007] and Nt5e (CD73) [Kobie et al, 2006]) and a reduced activation phenotype (Tnfrsf9medCd83med) without expression of effector proteins compared with other activated states suggest that Nr4a2+ Tregs are activated Tregs that are not terminally differentiated but still have suppressive capabilities. Compared with C1-resting Tregs, C4-cluster Tregs express high levels of effector molecules (GzmB, Glrx, and Klrg1) and chemokine receptor (Ccr2), suggesting this cluster may contain terminally differentiated Tregs with tissue-homing receptors (Cheng et al, 2012) (Fig 4C–F).

The gut Tregs contain small proportion of C1-resting Tregs and very little C2 or C9 prime/activated Tregs. Instead, they primarily contain C5- and C7-activated Treg populations. The C5-IL10hi Tregs uniquely express effectors such as IL10 and tissue-repairing Areg (Rubtsov et al, 2008) (Figs 3E and S9C) and proinflammatory transcriptional regulator Rora (Schiering et al, 2014) and elevated levels of gut-homing chemokine receptor Ccr9 (Iwata et al, 2004). The C7-Dnajb1hi–activated Tregs express early response inflammatory (Dnajb1 [Regateiro et al, 2012], Hspa1a [Wachstein et al, 2012; Collins et al, 2013], Tagap [Arshad et al, 2018], and interferonresponse–related) genes, some of which are associated with enhanced Treg function and autoimmunity (Fig S10A). Given that cells in the gut regularly encounter diverse commensal microbes and food antigens, these Tregs may actively promote tolerance even during steady-state conditions.

Reorganization of Treg landscape in multiple tissues after treatment with murine IL-2 mutein

We demonstrated the expansion of Tregs after murine IL-2M in vivo (Fig 1), but whether this treatment will expand all Treg sub-populations equally is not clear. Therefore, we compared the frequency of Treg states in each tissue before and after IL-2M to detect changes in molecular phenotypes and function. Treatment with IL-2M shifted the frequency of Treg clusters, reducing C1-resting and C3-activated Tregs while increasing primed/activated (C2) and activated Treg states (C4 and C8) (Fig 5A [distribution of all Tregs from isotype versus IL-2M–treated mice], Fig 5B and C [distribution of specific tissue Tregs from isotype versus IL-2M–treated mice]). These shifts were consistent in all animals (Fig 5C) and reflected differences at the tissue level, as lung and spleen Tregs shared most of the genes after IL-2M treatment (Fig S11A), but they were distinct when compared with gut Tregs (Fig S11B and C). In the spleen, lung, and gut, C2-primed Tregs and C4-activated Tregs all displayed an more than fourfold increase in all animals (Fig 5D). This expansion of C2- and C4-Tregs in the spleen and lung established these two states as the most prominent Treg populations after IL-2M stimulation, coinciding with a drastic reduction in C1-resting Tregs in both tissues and in C3-activated Tregs in the spleen. A twofold increase was also observed in C8-activated Tregs in the spleen and lungs. The selective expansion of C4- and C8-Tregs, which both co-express Tnfrsf9 and Il1rl1, suggests that IL-2M may skew expansion in eligible tissues toward specific activation phenotypes. Gut Tregs still maintained C5-activated Tregs after IL-2M treatment as the most prominent cell state (Fig 5C and D).

Figure S11.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S11. Transcriptional profiles converge in the spleen and lung after IL2-mutein treatment, whereas gut retains a distinct transcriptional profile from the spleen and lung.

(A, B, C) Comparison of Treg signature genes between all Tregs in the (A) spleen versus lung, (B) spleen versus gut, and (C) lung versus gut after IL-2M treatment. Significant genes are shown in red (adjusted P-value < 0.05, MAST).

IL-2 mutein skews clonal Treg expansion toward specific activation states

During development in the thymus, each T cell generates a unique TCRα and -β rearrangement (or clonotype), which predicates antigen specificity and influences Treg fate. This clonotype is retained in cells that descend from a common progenitor T cell, enabling characterization of Tregs that share a developmental lineage. Because we observed an expansion in Tregs (1) toward distinct cell states (C2, C4, and C6) in all three tissues, (2) toward the C8 state in the spleen and lungs, and (3) away from the C3 state in the spleen, we analyzed TCR information combined with Treg cell state classifications from single cells to elucidate whether these cell states expand independently of one another or whether they are developmentally related.

From the 3,600 sequenced single-cell transcriptomes where TCRα and TCRβ were also recovered (see the Materials and Methods section), 3,405 unique Treg clonotypes (1,428 in Iso- and 1,977 in IL-2M–treated conditions) were identified. Of these clonotypes, 119 CFs (i.e., clonotypes that were shared by at least two Tregs) consisting of 314 total Tregs were identified. Among the 119 CFs, 67 CFs of Tregs existed in different phenotypically distinct activated states (e.g., C3 and C4) or in primed/activated and activated states (e.g., C2 and C4), at both steady state and after IL-2M stimulation (Fig 6), revealing that different Treg cell states can be derived from the same parental Treg. The identification of transcriptional diversity among Tregs from the same clonal family (CF) was an interesting result because we also find that pairs of T cells belonging to the same clonotype tend to be transcriptionally correlated than randomly sampled pairs of Tregs at the population level (P = 2.106 × 10−4 when comparing equal numbers of cells from same and different clonotypes), although this was a modest effect, as previously reported (Zemmour et al, 2018) (Fig S12). CF membership also revealed evidence of Treg migratory patterns between tissues because we observed that Tregs of the same clonotype spanned the spleen, lung, and gut (Fig 6). Taken together, we interpret these findings as such that, although the TCR responses may generally influence the transcription of Tregs toward a programmed transcriptional fate, clonally related Tregs still can differentiate into multiple cell states depending on the tissue.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Tregs from the same progenitor differentiate across multiple states within and across tissues.

Chord diagrams of TCR clonotype sharing involving Tregs within the spleen (left), lung (middle), and across tissues (right) in control and IL-2M–treated condition. Outer ring denotes the tissue from which a Treg was isolated. Inner ring denotes the cell state identified by Seurat cell clustering. Green lines link splenic Tregs of the same clonotype. Blue lines link lung Tregs of the same clonotype. Purple lines link clonotypes shared between different tissues. P-values show Fisher’s exact test for independence between clonal family frequency and IL-2 mutein treatment.

Figure S12.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S12. Tregs belonging to the same clonotype share modest yet significant transcriptional programs.

Cell–cell correlation analysis between the transcriptomes of Treg cells that express the same rearranged TCRαβ (n = 478 cell–cell correlations) or unrelated TCRs (n = 10,000 correlations between any Tregs in the data not sharing clonal family membership, randomly sampled). ***P < 0.0001 (two-tailed t test).

We next characterized cell state classifications of clonotype pairs to identify cell states that may be developmentally linked (which we refer to as “cell state axes”). At the steady state, no cell state axis comprised more than 1% of all clonotypes because the majority of clonotypes (95.6%) did not belong to a CF (Fig 7B, top panel). In contrast, IL-2M–induced clonal proliferation facilitated the detection of prominent cell state axes connecting clonotype pairs by increasing the proportion of Tregs detected in a CF by 105% (211 CF Tregs in IL-2M compared with 103 CF Tregs in isotype), despite only a 41% increase in Tregs with TCRs recovered (2,107 TCRs in IL-2M, 1,493 TCRs in isotype). Multiple axes higher than 1% frequency were identified after IL-2M treatment (Fig 7B, bottom panel). Among those, multiple clonotypic pairs were observed in the C4-Tregs, including proliferation-dependent C4-C4 (3.9%) and C4-C6 (2.7%) clonotype pairs, as well as C2-C4 (1.0%) and C3-C4 (1.1%) pairs, revealing that C2, C3, and C4 are developmentally linked (Fig 7A and B). C4 also shared clonotypes with activated states C7 and C8, albeit at a lower frequency of 0.5% and 0.2%, respectively.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7. Treg lineage tracing and trajectory analyses identify axes of Treg differentiation across cell states.

(A) Network graph representing the percentage of shared clonotype pairs (i.e., pairs of cells belonging to the same clonotype that exist across two cell states) of spleen and lung Tregs under isotype control– or IL-2M–treated conditions. The thickness and darkness of the gray edges connecting two cell states corresponds to the percentage of shared clonotypes between two cell states or within a single state. (B) Top 15 pairs of cell states that share clonal family membership. Numbers above each bar indicate the total number of clonotype pairs in a given cell state axis. P-values show Fisher’s exact test for independence between matched cell state pairs frequency and IL-2 mutein treatment. (C) Comparison of the percentage of inter-tissue and intra-tissue clonotypes between isotype control– and IL-2M–treated conditions. Fisher’s exact test was used to test for the effect of IL-2M treatment on inter-tissue or intra-tissue frequency. (D) Pseudotemporal ordering of Tregs with recovered TCRs (n = 3,600 cells). Individual cells are colored by cell state (left) or by pseudotime or progress along Treg differentiation (right). (E) Histograms showing the distribution of each Treg cell state along the pseudotime axis. (F) Heat map showing the Z-score scaled expression of 30 highly variant genes along pseudotime. Genes are grouped into modules with similar expression by unsupervised hierarchical clustering.

Given that IL-2M increases clonal Treg expansion, we also examined how IL-2M influences the localization of CF Tregs by comparing the frequency of CF Tregs that were shared across tissues versus within the same tissue. Although the percentage of inter-tissue CFs remained the same in both isotype and IL-2M conditions (3.9% versus 4.2%, respectively), the percentage of CFs found only in one tissue was nearly doubled (3.0% versus 5.8%, respectively) (Fig 7C). Although we cannot ascertain whether in situ expansion or migration from other tissues is responsible for this phenomenon, these results indicate that CF Tregs after IL-2M–induced clonal expansion are more frequently observed in the same tissue than across multiple tissues.

C4 and C8 Tregs expand after IL-2M stimulation through C2 intermediate state

Given the gene expression profiles and robust lineage relationship of the C2-primed states and C3/C4/C8 activated states, we used pseudotime analyses using Tregs with recovered TCRs (n = 3,600 cells) to define their developmental relationship (Figs 7D and E and S13A). Treg cell states occupied distinct territories in pseudotime. We defined the node enriched for C1-resting Tregs as the root node (pseudotime t = 0), and pseudotime values were assigned in an unbiased manner to the manifold based on the distance from that root node. As expected, C1-resting Tregs occupied the earliest period of pseudotime. C2-primed/activated and C3-activated Tregs express markers of previous TCR signaling and were dispersed throughout intermediate points in pseudotime. At the latest periods in pseudotime/differentiation, we observed two distinct differentiation branchpoints consisting of C5-Tregs at one terminus and C4/C8 Tregs at the other (Fig S13B). This suggests that Tregs generally follow a differentiation trajectory from C1 into either C5 Il10-producing Tregs or C4/C8 Treg that express effector molecules such as Klrg1, Gzmb, and Glrx. We obtained similar ordering of cell states when applying pseudotime analyses to the other Tregs in our dataset where no TCR information was obtained (Fig S13C and D).

Figure S13.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S13. Trajectory analyses show consistent Treg differentiation patterns across replicates and gene expression dynamics along Treg differentiation.

(A) Two representative clonal families involving Tnfrsf9+Il1rl1+ Tregs (clonal family 836) and Il10-producing Tregs (clonal family 493) are plotted along the trajectory manifold. Individual cells are colored by cell state. (B) Individual Tregs with recovered TCRs are grouped by Treg cell state and plotted separately along the trajectory manifold. (C) Pseudotemporal ordering of Tregs from replicate 1 (n = 12,840 cells) without recovered TCRs. Individual cells are colored by cell state (left) or by pseudotime (right). (D) Pseudotemporal ordering of Tregs from replicate 2 (n = 11,822 cells). Individual cells are colored by cell state (left) or by pseudotime (right).

The analysis of Treg trajectories also revealed gene expression dynamics that advance Tregs toward effector states. Of the 30 most variant genes identified from this analysis, four major gene modules were identified that correspond to cell-state classifications (Fig 7F). Module 1 consisted of Izumo1r and Nr4a2, the key markers of C3-Tregs, which are moderately expressed all throughout pseudotime until terminal differentiation. The expression pattern of Nr4a2 agrees with previous findings that establish its role in tonic TCR signaling to allow Treg persistence and survival (Sekiya et al, 2011). Module 2 consisted of genes most frequently found in resting (C1) and primed/activated (C2/C9) Tregs, which are highly expressed at early cell states and become quickly down-regulated once differentiation initiates. Modules 3 and 4 consisted mainly of genes that were highly expressed in the C4/C8 and C5 terminal states, respectively. Overall, analysis of clonal Treg differentiation trajectories suggests IL-2M promotes differentiation into the terminally differentiated C4 and C8 Treg state by expanding through an intermediate state such as C2 in the spleen and lung.

Among IL-2M expanded sub-populations, Tnfrsf9+Il1rl1+ Tregs demonstrated superiority in in vitro suppression functional assay

Pseudotime analysis demonstrated enrichment of genes in four different modules (Fig 7F). Among those genes, Il1rl1 and Tnfrsf9 are highly expressed in Module 3, and they are highly expressed in C4/C8 clusters (Fig 3E). Likewise, scRNA-seq demonstrated increased expression of Il1rl1 and Tnfrsf9 in the spleen and lung after the IL-2M treatment (Fig 8A). Il1rl1 or Tnfrsf9 encodes proteins ST2 or 4-1BB, respectively, and FACS reagents recognizing these proteins were readily available. To verify protein expression of these two genes, we performed flow cytometry analysis. Flow cytometry analysis confirmed that Foxp3+ Tregs expressing both Il1rl1(ST2) and Tnfrsf9 (4-1BB) were increased to ∼10-fold from the lung and spleen of the mice treated with IL-2M at protein level (Fig 8A–C). Pan-Tregs isolated after IL-2M treatment increased the suppressive activity (Fig 2A and B), but it is not clear which sub-state Tregs contribute to this enhanced suppressive activity. Because Tregs with expression of Il1rl1 and Tnfrsf9 after IL-2M treatment represent the end stage of the differentiation, we determined suppressive activity of Tregs with differential expression of ST2 and 4-1BB using in vitro suppression assays examining suppression activity of Tregs on proliferation and IFNγ production. To this end, Tregs were expanded in vivo by murine IL-2M administration. 4 d after treatment, four populations of Foxp3+ Treg cells based on the expression of ST2 and 4-1BB were sorted, and in vitro suppression assay was performed. Proliferation of effector T cells was suppressed by Tregs expressing either ST2 or 4-1BB, but ST2+4-1BB+ Foxp3+ Tregs displayed the most superior suppression (Figs 8D). and S14A and B). Similarly, ST2+4-1BB+ Foxp3+ Tregs were most efficient in suppressing IFNγ production of T effector cells under Th1 skewing condition (Fig 8E). Combined with increased suppression of Foxp3+ Tregs expanded with murine IL-2M in Fig 1, these data demonstrated that murine IL-2 mutein treatment preferentially increased the proportion of ST2+4-1BB+ Foxp3+ Tregs, which possess superior ability to suppress proliferation and IFNγ production of T effector cells. Thus, endogenous expansion of Tregs using in vivo murine IL-2M administration has the potential either to induce or expand Foxp3+ Tregs with the highest capacity to suppress autoreactive T cells.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8. Tnfrsf9+Il1rl1+ Tregs are superior in suppression functional assay of Tregs in vitro.

(A) Expression of Il1rl1(ST2) and Tnfrsf9(4-1BB) in RNA single-cell analysis of Tregs from the spleen (top) or lung (bottom) of isotype control–treated (Iso) or murine IL-2M–treated mice (IL-2M). (B) Expression of ST2 and 4-1BB using FACS analysis from the spleen (top) or lung (bottom) of Iso or IL-2M–treated mice. (C) Quantification of ST2+4-1BB+ cells in Foxp3+ Tregs. (B, C) Results are representative of two independent experiments, using three mice in each experiment. (D) Fopx3-EGFP+ Tregs were stained based on ST2 and 4-1BB expression. Single-cell suspension from two IL-2M–treated mice was combined and four quadrants of Fopx3-EGFP+ Tregs were sorted based on ST2 and 4-1BB. In vitro suppression assays were performed in duplicate using four populations of sorted Tregs or without Tregs. Experiments were repeated using three mice and in triplicate (data not shown). (E) Suppression function of four Treg sub-populations on IFNγ production. Four quadrants of Fopx3-EGFP+ Tregs were sorted based on ST2 and 4-1BB. Tregs were cocultured with effector T cells and APCs differentiated under Th1 skewing condition. IFNγ intracellular staining under the conditions using Tregs from PBS- or IL-2M–treated mice. No Tregs were added in one condition as a control. (B, C, D, E) Results are representative of two independent experiments, using 2–3 mice in each experiment. (C, D) Statistics: (C, D), one-way ANOVA for multiple comparisons. **0.001 < P < 0.01, ***0.0001 < P < 0.001, ****P < 0.00001.

Figure S14.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S14. Tnfrsf9+Il1rl1+ Tregs are superior in suppression functional assay of Tregs in vitro.

(A) Fopx3-EGFP+ Tregs were stained based on ST2 and 4-1BB expression. Single-cell suspensions from three IL-2M–treated mice were combined and four quadrants of Fopx3-EGFP+ Tregs were sorted based on ST2 and 4-1BB. In vitro suppression assays were performed in triplicate using four populations of sorted Tregs or without Tregs. (A, B) In vitro suppression assays were performed in triplicate using the sorted Tregs described in (A) under Th1 skewing condition. One-way ANOVA for multiple comparisons. 0.01 < *P < 0.05, **0.001 < P < 0.01, ***0.0001 < P < 0.001, ****P < 0.00001.

Discussion

IL-2 plays an indispensable role in immune tolerance by governing the proliferation and maintenance of Foxp3+ Tregs. IL-2 GWAS studies identified multiple polymorphisms affecting IL2 pathway genes encoding IL2Rα (CD25) and IL-2 as key genetic risk alleles in autoimmunity, corroborating IL2 as an essential cytokine for maintaining immune homeostasis (Vella et al, 2005; International Multiple Sclerosis Genetics et al, 2007; Yamanouchi et al, 2007; International Multiple Sclerosis Genetics, 2008; Onengut-Gumuscu et al, 2015; Wallace et al, 2015). Together with IL-2 therapy studies and mouse models defective in IL2 signaling, it is now well accepted that IL2 is the key cytokine necessary for generation, function, and survival of Foxp3+ Tregs. Consistent with genetic studies in human and mouse, reduction in numbers and functions of Foxp3+ Tregs had been reported as a contributing factor for a variety of autoimmune diseases (International Multiple Sclerosis Genetics, 2008; Long et al, 2011; Garg et al, 2012; Cerosaletti et al, 2013; Carbone et al, 2014; Moulton & Tsokos, 2015; Yang et al, 2015).

Given its role in Foxp3+ Treg biology, providing IL2 to Foxp3+ Tregs as an intervention for autoimmune diseases had been tested in clinic over the years. Recent clinical trials using low-dose recombinant human IL2 (Proleukin, aldeskeukin) have shown preliminary success in selectively expanding Tregs in chronic graft-versus-host-disease (Koreth et al, 2011; Hartemann et al, 2013; Kennedy-Nasser et al, 2014) and in steroid-refractory moderate-to-severe systemic lupus erythematosus patients (von Spee-Mayer et al, 2016). Additional efforts have focused on novel forms of IL2 with superior pharmacokinetics and Treg selectivity. These include a long-lived bivalent IgG-human IL2 fusion protein with N88D mutation that reduces IL2Rβ binding (Peterson et al, 2018), a polyethylene glycol (PEG)-modified murine IL2 that increases IL2 retention in the body by protection from enzymatic digestion and renal clearance (Wu et al, 2016), and a murine IL2/anti–IL2 antibody (JES6-1) conjugate (Spangler et al, 2018). Despite the effort to develop various forms of IL2 with superior pharmacokinetics and selectivity for treatment of autoimmune and inflammatory diseases, it is not clear how they affect heterogeneity of Tregs, given that Tregs exists in a variety of sub-states at different tissue sites, which are phenotypically and functionally heterogeneous.

To establish an effective Treg-mediated therapy, it is critical to characterize the phenotypic and functional heterogeneity induced by the treatment and define mechanisms that promote functional Treg differentiation. Because of the diverse profiles of proteins expressed in Tregs, it is challenging to define functional Tregs, let alone to identify Tregs with superior effector functions. CD4+CD45RAlowCD25hiCD127low cells have been previously used as human surface Treg markers to identify the pan-Tregs, and various activated markers of Tregs have been adopted to understand the functional and phenotypic diversity of Tregs. However, the detailed molecular characterization of Treg cell states is limited with a mere handful of activation markers. scRNA-seq overcomes this obstacle by providing a comprehensive view of transcriptional heterogeneity of Tregs, thereby making it possible to better design a targeted Treg therapy.

Using scRNA-seq and ex vivo cellular assays, we mapped the landscape of Treg sub-states at the single cell level during steady state and after murine IL-2M treatment in mice. We used an engineered, half-life extended murine IL2 mutant protein (IL-2M) and demonstrated that it expanded Tregs robustly in vivo. Expanded Tregs by murine IL-2M treatment suppressed proliferation of Tconv as well as production of IFNγ from Th1 effector cells. scRNA-seq study uncovered that IL-2M induces a convergent shift of the Treg landscape toward primed (C2) and highly suppressive activated Treg states (C4 and C8) in addition to proliferative Treg state (C6). Of note, because expansion of Tregs in the C2/C4/C8 clusters after IL-2M treatment is the result of action of IL-2M, we cannot distinguish which cluster of Tregs responded to the IL-2M to give rise to the expansion of the C2/C4/C8 clusters. Likewise, it is possible there is a preferential increase in a specific cluster to respond to IL-2M, resulting in the expansion of the C2/C4/C8 clusters. Furthermore, TCR-seq combined with scRNA-seq revealed a migratory network across tissues. IL-2M–mediated Treg expansion increased shared TCR clonality, and this expansion was more prominent in CFs that were identified in one tissue, whereas the percentage of inter-tissue CFs remained the same in both isotype and IL-2M treatment. Interestingly, shared TCR clones were also found in multiple and distinct cell states, providing us with a molecular snapshot of the transcriptional changes after IL-2M stimulation.

The similarities in single cell profiles between the spleen and lung versus the gut at steady state suggest that the spleen and lung in healthy mice possess a large proportion of resting-Tregs (C1-resting Tregs) that serve as a reservoir of Tregs in preparation for potential inflammatory stimuli. The spleen uniquely contains a high frequency of activated Tregs (C3-activated Tregs) expressing Nr4a2, a Foxp3-binding transcription factor (Kamalipour et al), Cst7 (Hamilton et al, 2008), a marker for previous Treg activation and regulator of cytotoxicity, and Izumo1r (Yamaguchi et al, 2007), a marker for natural Tregs in the spleen. Compared with other activated Treg cluster C4, C3-activated cluster also differentially express Tbc1d4 and Maf, transcription factors for follicular regulatory T cells (Wheaton et al, 2017; Wing et al, 2017). Given the role of the spleen in immune surveillance and its lymphatic connection to peripheral tissues as well as germinal center formation (Lim et al, 2004), these splenic Nr4a2+ Tregs may include follicular regulatory T cells (Tfr) and recently activated Tregs licensed to exit the spleen and migrate into sites of inflammation. On the other hand, the gut is composed of uniquely activated Treg clusters, IL10-producing suppressive C5 and inflammatory/stress response-associated C7 clusters, likely because of its constant exposure to commensal microbes, pathogens, and food antigens, it requires many activated Tregs that can maintain tissue homeostasis and promote tolerance (Geuking et al, 2011; Lathrop et al, 2011).

IL-2M treatment reduced the proportion of resting C1 Tregs in the spleen and lungs and the C3 Tregs in the spleen, whereas expanding the numbers of proliferating (C6), primed/activated Tregs (C2) and potently suppressive Tnfrsf9+Il1rl1+ activated Tregs (C4 and C8). Because we observed a more than twofold expansion in the number of Tregs in all tissues after IL-2M stimulation, the decrease in C1 and C3 Treg percentage was likely driven more by the concomitant expansion of various expanding sub-populations as opposed to a significant reduction in the absolute numbers of the prominent Treg populations during steady state. Both expanded activated Treg populations, C4 and C8, express high Il1rl1 (ST2), which is the receptor for tissue alarmin IL-33 (Siede et al, 2016). ST2+ Tregs have been previously shown to prevent excessive tissue damage in organs such as the skin, liver, and gut (Matta et al, 2014; Schiering et al, 2014; Gajardo et al, 2015). C4 and C8 Tregs, although both expressing Tnfrsf9 and Il1rl1, as well as effector molecules such as Klrg1 and Gzmb, differ in expression of regulatory mediators such as Ass1, Cd52, and Tnfrsf18. Trajectory analysis suggests that the C4 and C8 Tregs are terminally differentiated. Previously, it has been shown that Klrg1+ Tregs represent a terminally differentiated Treg subset that is recently activated by antigen and resides in the lamina propria of small intestine (Cheng et al, 2012). Klrg1+ Tregs are highly activated and express enhanced levels of suppressive Treg molecules. Interestingly, the development of Klrg1+ Tregs requires extensive IL-2R signaling (Cheng et al, 2012). Based on gene expression, it is highly likely that C4 and C8 Tregs contain Klrg1+ Tregs. Between the two clusters, the C8 appears to be more terminally differentiated because of its expression of effector cell surface receptors Cd52 and Tnfrsf18-encoded GITR (Bandala-Sanchez et al, 2013; Ephrem et al, 2013), which suppress other immune cells by binding to ITIM-containing Siglec10 and GITR-L, respectively. Identification of these expanded clusters suggests that IL-2M may significantly remodel the Treg landscape toward more terminally differentiated effector Tregs and suppressive Tregs, which mediate potent immunosuppression in these tissues.

We provide evidence using TCR analyses and trajectory analyses to suggest that C4 and C8 Tregs differentiate from C2/C3 primed/activated states. TCR analyses identify that C2, C3, C4, and C8 Tregs in the spleen and lung share clonotypes, demonstrating that they can all derive from a common progenitor. Cells from these populations that share the same clonotype are also present across the spleen and the lung, revealing an immune-trafficking axis that exists between the two tissues. Trajectory analysis also identifies a bifurcation in Treg differentiation, which either differentiate into suppressive Il10+Rora+ C5 Tregs, which are most prevalent in the gut, or into C4/C8 Tregs that are prominent in the spleen and lungs. However, the trajectory analysis does not inform about whether IL-2M actually directs this bifurcation.

In this study, we characterize Tregs in the spleen, lung, and gut using scRNA-seq and analyze how the Treg landscape in these tissues is influenced by IL-2M stimulation. We identified an immune axis of clonal TCRs between the spleen and the lung and found an activated, Tnfrsf9+Il1rl1+ Treg population that expand highly upon IL-2M stimulation. Supporting analyses indicate that this population differentiates through multiple primed and recently activated intermediate cell states and it is a highly suppressive Treg state that inhibits Th1 effector cell activity ex vivo. The mechanism of increased suppressive activity by IL-2M remains to be elucidated further. We demonstrated that IL-2M expands Tnfrsf9+Il1rl1+ Tregs with superior suppressive function. Another non-mutually exclusive possibility is that IL-2M induces transcriptional changes improving suppressor capacity within Tnfrsf9+Il1rl1+ Tregs. Overall, this study reveals a potential mechanism by which IL-2M induces tissue-specific Treg immune modulation and classifies Treg state signatures that may serve as biomarkers for studying Treg responses. Although our analysis focused on genes with previously characterized immune cell functions, previously uncharacterized genes that were associated with Treg states of differentiation, such as certain ribosomal genes, warrant further investigation to elucidate their role in Treg biology.

Materials and Methods

Mice and IL-2 mutein treatment

All experimental studies were conducted under protocols approved by the Institutional Animal Care and Use Committee of Amgen. Animals were housed at Association for Assessment and Accreditation of Laboratory Animal Care International–accredited facilities at Amgen in ventilated micro-isolator housing on corncob bedding. Animals had access ad libitum to sterile pelleted food and reverse osmosis–purified water and were maintained on a 12:12 h light:dark cycle with access to environmental enrichment opportunities. Foxp3-EGFP (#006772), C57BL/6J, and B6.SJL-PrprcaPepcb/BoyJ (BoyJ) mice (8–12 wk old) were purchased from Jackson Laboratory. Once purchased, mice were housed under specific pathogen-free conditions in the laboratory animal facility and were handled according to protocols approved by Institutional Animal Care and Use Committee at Amgen. To determine the effect of murine IL-2 mutein (IL-2M), different amounts of IL-2M (0.1, 0.33, and 1 mg/kg) were administered into naïve C57BL/6J mice. After 4 d of IL-2M treatment, flow cytometry was performed to determine its effect on the regulatory and conventional T cells. To determine the effect of murine IL-2M on T-cell proliferation, Ki67 staining was performed. Untreated, PBS- or mouse IgG Fc isotype (31205; Thermo Fisher Scientific)–treated mice were compared and used as negative controls. Administration of murine IL-2M via i.p. or s.c. routes resulted in similar effects. For scRNA-seq experiment, the mice were administered with murine IL-2 mutein (0.33 mg/kg) or mouse IgG Fc isotype control (Iso) by i.p. injection and were sacrificed 4 d after murine IL-2M administration.

IL-2/anti-IL2 antibody complex (IL-2C)

IL-2 complex (IL-2C) was prepared as described (Webster et al, 2009) by incubating 1 μg recombinant mouse IL-2 (PeproTech) with 5 μg purified antimouse IL-2 antibody (JES6-1; BioXcell) for 30 min at 37°C (molar ratio for IL2:anti-IL-2 is 2:1). IL-2C (total 6 μg) was administered i.p. in a final volume of 200 μl daily for 3 d. 4 d after the last administration of IL-2C, expansion of Tregs and Tconv was determined by FACS analysis.

Flow cytometry reagents and Treg sorting

Anti-CD11c (N418), anti-F4/80 (BM8), anti-CD19 (6D5), anti-CD11b (M1/70), and anti-CD45.1 (A20) were from BioLegend. Anti-ST2 (U29-93), anti-CD3 (17A2), anti-CD4 (RM4-5), anti-IFNγ (XMG1.2), and anti-CD16/32 (2.4G2) were from BD Biosciences. CellTrace Violet, Fixable Viability Dye eFluor 780, and Foxp3 (FJK-16s) were from Thermo Fisher Scientific. Anti-4-1BB (158332) was from R&D.

Isolation of murine spleen and lung Treg cells for scRNA-seq

The spleen and lung were collected and processed according to the protocol from Miltenyi murine Lung Dissociation kit (130-095-927). Lungs were perfused via the right ventricle before the process. Briefly, tissues were cut and rinsed in PBS. The tissues were transferred into the gentleMACS C Tube containing an enzyme mix of Enzyme D and A. gentleMACS Program m_lung_01 was run twice and the samples were incubated for 30 min at 37°C under continuous shaking using the gentleMACS dissociator (Miltenyi Biotec) at 150 rpm. Then, gentleMACS Program m_lung_02 was run. The tissue disaggregate was filtered through 70-μm strainers, and the cells were washed with PEB (PBS [pH 7.2] 2 mM EDTA, and 0.5% bovine serum albumin) buffer twice. The samples were ready for downstream flow staining and sorting.

Isolation of murine gut Treg cells for scRNA-seq

The small intestine (extending from the end of stomach to the cecum) was collected and processed according to the protocol from Miltenyi murine Lamina Propria Dissociation kit (130-097-410). Briefly, the feces and Peyer’s patches were removed, and the small intestine was filleted and cut into pieces. The samples were washed in predigestion solution, vortexed, and passed through 100-μm filters twice. Then, the samples were washed in HBSS (w/o), vortexed, and passed through 100-μm filters. The samples were transferred to gentleMACS C Tube containing preheated 2.35 ml of digestion solution with Enzyme D, R, and A. The samples were incubated for 30 min at 37°C under continuous shaking using the gentleMACS dissociator (Miltenyi Biotec) at 150 rpm. The C Tube was plugged into gentalMACS and run Program m_intestine_01 twice. The tissue disaggregates were filtered through 100-μm strainers, and the cells were washed with PB (PBS [pH 7.2] and 0.5% bovine serum albumin) buffer twice. The samples were ready for downstream flow staining and sorting.

Treg suppression assay

Various populations of Treg cells were FACS-sorted on Aria II sorter (BD Biosciences) based on their expression of EGFP reporter. Naïve CD4+ T cells were isolated by Miltenyi mouse naive CD4+ T cell isolation kit and were used as effector T cells. APCs were prepared by depleting CD90+ T cells from splenocytes and irradiated for 3,000 rads (Faxitron). Naïve CD4+ T cells and APCs were from BoyJ mice, which is CD45.1+ to differentiate from EGFP+ Treg cells which were in B6 background by flow cytometry. To test the suppression function of Treg cells on proliferation, various populations of Treg cells (1 × 104/well) were cocultured with effector T cells (2 × 104/well) and APCs (3 × 104/well), as indicated, along with anti-CD3 (145-2C11, 1 μg/ml; BioXcell) for 72 h. To test the suppression function of Treg cells on IFNγ production, various populations of Treg cells (2,500/well) were cocultured with effector T cells (2 × 104/well) and APCs (3 × 104/well), along with anti-CD3 (1 μg/ml), anti-IL-4 (11B11, 10 μg/ml; BioXcell), and IL-12 (6 ng/ml; PeproTech) for 96 h. Then, the whole culture will be stained for intracellular cytokine staining. Briefly, the cells were washed once and were stimulated with cell stimulation cocktail (eBioscience) for 5 h in cell culture medium, and then fixed and stained for IFNγ.

Single-cell RNA-seq

scRNA-seq was performed in three replicates, with each replicate containing the spleen, lung, and gut Tregs from one isotype control–treated mouse (Iso) and one IL2-mutein treated mouse (IL-2M). Combining different treatment conditions and tissues in each replicate allowed us to identify and regress out gene expression differences driven by batch processing.

CD4+Foxp3+ single cells were sorted from tissue-dissociated single-cell suspensions after gating on CD14−CD19−CD11c−F4/80− cells. After sorting, the cells were immediately washed and resuspended in chilled PBS + 0.04% BSA at 100–1,000 cells per microliter. Washed cells were then encapsulated in one lane of a 10× Chromium Single Cell Controller, and libraries were constructed using either the 10× Single Cell 3′ Reagent Kit (V2 chemistry) (https://support.10xgenomics.com/single-cell-gene-expression/library-prep/doc/user-guide-chromium-single-cell-3-reagent-kits-user-guide-v2-chemistry) or the 10× Single Cell 5′ Reagent Kit (https://support.10xgenomics.com/single-cell-vdj/library-prep/doc/user-guide-chromium-single-cell-vdj-reagent-kits-v1-chemistry) according to the manufacturer’s protocols. Completed libraries were sequenced on the HiSeq 4000 or NovaSeq 6000 platforms (26/8/0/98 Read1/i7/i5/Read2). Data were analyzed using the Cell Ranger 2.0.0 Pipeline by aligning reads to a genome containing mm10 genome and eGFP (Table S1).

Post-alignment filtering and data aggregation

The Seurat (version 2.4) implementation (Satija et al, 2015) in Scanpy (version 0.94) (Wolf et al, 2018) was used to aggregate data from multiple experimental replicates and cluster Tregs into distinct states. Cells with low-quality transcriptomes (<500 detected genes) and doublets (>8,000 genes) were removed from the analysis. Canonical correlation analysis from the Seurat package was performed to identify common sources of variation and align data from different library preparation batches (each batch consisted of one isotype and one IL-2M treated sample of Tregs from each tissue), thereby mitigating batch effects. After filtering, our dataset detected 18,367 total genes and 31,908 cells across all replicates (12,840 cells in replicate 1, 11,822 cells in replicate 2, and 7,246 cells in replicate 3). We recovered an average of 5,032 unique mRNA molecules after collapsing duplicate unique molecular identifiers and 1,635 genes per cell.

Unsupervised cell clustering and differential expression analysis of clusters

Unsupervised cell clustering of scRNA-seq data was performed using the FindClusters() function from Seurat. The optimal clustering resolution in Seurat was determined by clustering integrated single-cell expression data at 10 different resolutions from 0.1 to 1.0 using the “resolution” parameter in the FindClusters() function. At each increasing resolution, the top marker genes of the cluster containing the fewest cells were evaluated against previously published literature to support inclusion. The determination of 0.6 for the resolution parameter was made because of the identification of an IFNghi Treg cluster because lower resolutions did not identify this population. Additional clusters at higher resolutions did not reveal additional known immune markers.

Differential expression testing of each cluster (i.e., Treg cell state) versus all other clusters was performed using MAST (Finak et al, 2015) to generate a list of differentially expressed genes for each cluster. A list of genes that were significantly up-regulated genes by both tests were then used to define hallmark marker genes for each cell state. Hallmark genes for each of the Treg sub-populations was determined by selecting from among the top 10 genes with an adjusted P-value < 0.01, average log2-fold change expression >0.3, and statistically significant enrichment in less than half of the total cell clusters. The inclusion of genes significantly enriched in multiple clusters was set under the assumption that whereas the combinatorial marker gene expression is unique to each cluster, individual marker genes may be up-regulated by multiple cell states. Mitochondrial and pseudogenes were excluded from consideration as gene markers. T-stochastic neighbor embedding using RunTSNE() was used from Seurat for dimensionality-reduced visualization of single cell data.

Single-cell TCR-seq

For samples prepared using the 5′ Reagent Kit, 2 μl of amplified cDNA was set aside to amplify mouse TCRa and TCRb chains in two sequential target enrichment PCRs using custom reverse primers: 5′-TGAAGATATCTTGGCAGGTG-3′ and 5′-TGCTCAGGCAGTAGCTATAATTGCT-3′ were used for the first PCR and 5′-GATCTTTTAACTGGTACACA-3′ and 5′-TTTGATGGCTCAAACAAGGA-3′ were used for the second PCR. All forward primers, primer concentrations, and PCR conditions and cycles were the same as those specified in the manufacturer’s protocol. Completed libraries were sequenced on the HiSeq 4000 platform (150/8/0/150 Read1/i7/i5/Read2). Data were analyzed using the Cell Ranger 2.0 Pipeline by aligning reads to all IMGT mouse TCR variable and constant region sequences (Lefranc et al, 2009). Germline gene segments and CDR3 TCRa and TCRb sequences were identified and clonotypes were defined by single cells containing a recovered TCRa and TCRb with identical germline gene segments and CDR3 sequences. In total, 3,600 cells contained complete TCR and transcriptomic information out of 7,246 TCRs recovered using the 10× Chromium 5′-VDJ method.

Defining cell state networks using TCR clonotype and transcriptomic data

Custom R scripts were used to define relationships between TCR and gene expression information. Single cell TCR data were merged with transcriptomic data using the common cell nucleotide barcode incorporated during reverse transcription. All clonotype analysis was performed using cells with one single productive TCRa and TCR chain and transcriptomic information (filtered using Scanpy), whereas all other cells were excluded.

To generate chord diagrams to illustrate the clonal relationships and cell states shared by Tregs from the same clonotype, the circlize package in R was used. Tissue of isolation and cell state definitions determined by Seurat were used to construct the outer and inner tracks of the chord diagram. Links were drawn for any two cells sharing a clonotype.

To generate network plots depicting the magnitude of clonotype sharing across cell states, the igraph package in R was used in Fig 5A. The frequency of clonotypes across cell states was first calculated by summing the number of cell pairs shared between two cell states and dividing by the total number of clonotypes. This value was then used to determine edge widths connecting two cell states, or vertices. Only spleen and lung Tregs were included in this analysis, as we recovered low numbers of gut Treg TCRs.

Pseudotime analysis of Treg differentiation

Monocle2 (2.8.0) was used to order cells according to pseudotime. To generate the input data for the analysis, we combined cells from IL-2M and isotype conditions because (1) cell state classifications were determined independently of treatment and (2) we reasoned that the developmental trajectory should be preserved across treatment condition. We removed cells (i.e., C6 Tregs) whose primary marker genes were driven by proliferative status and not by their relative status along Treg-specific differentiation, although an independent analysis including these cells revealed that they are dispersed evenly throughout the trajectory manifold, which is expected because both resting and activated cell states are capable of proliferation. The gene list used for pseudotime analysis was filtered to include genes expressed in greater than 10% of all cells in the data (5,311 genes). We then used Monocle DDRTree to perform dimension reduction and manifold constructing, representing the trajectory of Treg differentiation. The cells were then ordered along the manifold. Genes that varied the most along the pseudotime axis were determined using the Monocle function: differentialGeneTest(fullModelFormulaStr = “∼sm.ns(Pseudotime)”).

Data Availability

The scRNA-seq data from this publication have been deposited to the European Molecular Biology Laboratory’s European Bioinformatics Institute public repository (https://www.ebi.ac.uk/ena/browser/home), and assigned the accession code PRJEB36020 (https://www.ebi.ac.uk/ena/browser/view/PRJEB36020).

Acknowledgements

We thank Olga Pryshchep, Hyun Ra, Nathan Deer, Jennifer Cheung, Laura Dieu, Aaron Fojas, Marisela Killian, Min-Zu Wu, Dev Bhatt, Oliver Homann, and Oh-Kyu Yoon for their technical support and intellectual input.

Author Contributions

  • DR Lu: formal analysis and investigation.

  • H Wu: formal analysis and investigation.

  • I Driver: data curation, formal analysis, and investigation.

  • S Ingersoll: validation and methodology.

  • S Sohn: validation, visualization, and methodology.

  • S Wang: conceptualization and resources.

  • C-M Li: resources, software, investigation, and methodology.

  • H Phee: conceptualization, resources, formal analysis, validation, investigation, and visualization.

Conflict of Interest Statement

DR Lu, S Sohn, S Wang, C-M Li, and H Phee are employees of Amgen. H Wu is currently employed at Pharmacyclics. I Driver is currently employed at Gordian Biotechnology. S Ingersoll is currently employed at Nektar therapeutics.

  • Received August 7, 2019.
  • Revision received March 30, 2020.
  • Accepted March 31, 2020.
  • © 2020 Lu et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Arpaia N,
    2. Green JA,
    3. Moltedo B,
    4. Arvey A,
    5. Hemmers S,
    6. Yuan S,
    7. Treuting PM,
    8. Rudensky AY
    (2015) A distinct function of regulatory T cells in tissue protection. Cell 162: 1078–1089. doi:10.1016/j.cell.2015.08.021
    OpenUrlCrossRefPubMed
  2. ↵
    1. Arshad M,
    2. Bhatti A,
    3. John P,
    4. Jalil F,
    5. Borghese F,
    6. Kawalkowska JZ,
    7. Williams RO,
    8. Clanchy FIL
    (2018) T cell activation Rho GTPase activating protein (TAGAP) is upregulated in clinical and experimental arthritis. Cytokine 104: 130–135. doi:10.1016/j.cyto.2017.10.002
    OpenUrlCrossRef
  3. ↵
    1. Bandala-Sanchez E,
    2. Zhang Y,
    3. Reinwald S,
    4. Dromey JA,
    5. Lee BH,
    6. Qian J,
    7. Bohmer RM,
    8. Harrison LC
    (2013) T cell regulation mediated by interaction of soluble CD52 with the inhibitory receptor Siglec-10. Nat Immunol 14: 741–748. doi:10.1038/ni.2610
    OpenUrlCrossRefPubMed
  4. ↵
    1. Boyman O,
    2. Kovar M,
    3. Rubinstein MP,
    4. Surh CD,
    5. Sprent J
    (2006) Selective stimulation of T cell subsets with antibody-cytokine immune complexes. Science 311: 1924–1927. doi:10.1126/science.1122927
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Butler A,
    2. Hoffman P,
    3. Smibert P,
    4. Papalexi E,
    5. Satija R
    (2018) Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol 36: 411–420. doi:10.1038/nbt.4096
    OpenUrlCrossRefPubMed
  6. ↵
    1. Cao X,
    2. Cai SF,
    3. Fehniger TA,
    4. Song J,
    5. Collins LI,
    6. Piwnica-Worms DR,
    7. Ley TJ
    (2007) Granzyme B and perforin are important for regulatory T cell-mediated suppression of tumor clearance. Immunity 27: 635–646. doi:10.1016/j.immuni.2007.08.014
    OpenUrlCrossRefPubMed
  7. ↵
    1. Carbone F,
    2. De Rosa V,
    3. Carrieri PB,
    4. Montella S,
    5. Bruzzese D,
    6. Porcellini A,
    7. Procaccini C,
    8. La Cava A,
    9. Matarese G
    (2014) Regulatory T cell proliferative potential is impaired in human autoimmune disease. Nat Med 20: 69–74. doi:10.1038/nm.3411
    OpenUrlCrossRefPubMed
  8. ↵
    1. Cerosaletti K,
    2. Schneider A,
    3. Schwedhelm K,
    4. Frank I,
    5. Tatum M,
    6. Wei S,
    7. Whalen E,
    8. Greenbaum C,
    9. Kita M,
    10. Buckner J, et al.
    (2013) Multiple autoimmune-associated variants confer decreased IL-2R signaling in CD4+ CD25(hi) T cells of type 1 diabetic and multiple sclerosis patients. PLoS One 8: e83811. doi:10.1371/journal.pone.0083811
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chaudhry A,
    2. Samstein RM,
    3. Treuting P,
    4. Liang Y,
    5. Pils MC,
    6. Heinrich JM,
    7. Jack RS,
    8. Wunderlich FT,
    9. Bruning JC,
    10. Muller W, et al.
    (2011) Interleukin-10 signaling in regulatory T cells is required for suppression of Th17 cell-mediated inflammation. Immunity 34: 566–578. doi:10.1016/j.immuni.2011.03.018
    OpenUrlCrossRefPubMed
  10. ↵
    1. Cheng G,
    2. Yuan X,
    3. Tsai MS,
    4. Podack ER,
    5. Yu A,
    6. Malek TR
    (2012) IL-2 receptor signaling is essential for the development of Klrg1+ terminally differentiated T regulatory cells. J Immunol 189: 1780–1791. doi:10.4049/jimmunol.1103768
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Collins CB,
    2. Aherne CM,
    3. Yeckes A,
    4. Pound K,
    5. Eltzschig HK,
    6. Jedlicka P,
    7. de Zoeten EF
    (2013) Inhibition of N-terminal ATPase on HSP90 attenuates colitis through enhanced Treg function. Mucosal Immunol 6: 960–971. doi:10.1038/mi.2012.134
    OpenUrlCrossRef
  12. ↵
    1. Daniel V,
    2. Trojan K,
    3. Adamek M,
    4. Opelz G
    (2015) IFNgamma+ Treg in-vivo and in-vitro represent both activated nTreg and peripherally induced aTreg and remain phenotypically stable in-vitro after removal of the stimulus. BMC Immunol 16: 45. doi:10.1186/s12865-015-0111-2
    OpenUrlCrossRefPubMed
  13. ↵
    1. Delgoffe GM,
    2. Woo SR,
    3. Turnis ME,
    4. Gravano DM,
    5. Guy C,
    6. Overacre AE,
    7. Bettini ML,
    8. Vogel P,
    9. Finkelstein D,
    10. Bonnevier J, et al.
    (2013) Stability and function of regulatory T cells is maintained by a neuropilin-1-semaphorin-4a axis. Nature 501: 252–256. doi:10.1038/nature12428
    OpenUrlCrossRefPubMed
  14. ↵
    1. Epeldegui M,
    2. Blom B,
    3. Uittenbogaart CH
    (2015) BST2/tetherin is constitutively expressed on human thymocytes with the phenotype and function of Treg cells. Eur J Immunol 45: 728–737. doi:10.1002/eji.201444787
    OpenUrlCrossRefPubMed
  15. ↵
    1. Ephrem A,
    2. Epstein AL,
    3. Stephens GL,
    4. Thornton AM,
    5. Glass D,
    6. Shevach EM
    (2013) Modulation of Treg cells/T effector function by GITR signaling is context-dependent. Eur J Immunol 43: 2421–2429. doi:10.1002/eji.201343451
    OpenUrlCrossRefPubMed
  16. ↵
    1. Finak G,
    2. McDavid A,
    3. Yajima M,
    4. Deng J,
    5. Gersuk V,
    6. Shalek AK,
    7. Slichter CK,
    8. Miller HW,
    9. McElrath MJ,
    10. Prlic M, et al.
    (2015) MAST: A flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biol 16: 278. doi:10.1186/s13059-015-0844-5
    OpenUrlCrossRefPubMed
  17. ↵
    1. Gajardo T,
    2. Morales RA,
    3. Campos-Mora M,
    4. Campos-Acuna J,
    5. Pino-Lagos K
    (2015) Exogenous interleukin-33 targets myeloid-derived suppressor cells and generates periphery-induced Foxp3(+) regulatory T cells in skin-transplanted mice. Immunology 146: 81–88. doi:10.1111/imm.12483
    OpenUrlCrossRef
  18. ↵
    1. Garg G,
    2. Tyler JR,
    3. Yang JH,
    4. Cutler AJ,
    5. Downes K,
    6. Pekalski M,
    7. Bell GL,
    8. Nutland S,
    9. Peakman M,
    10. Todd JA, et al.
    (2012) Type 1 diabetes-associated IL2RA variation lowers IL-2 signaling and contributes to diminished CD4+CD25+ regulatory T cell function. J Immunol 188: 4644–4653. doi:10.4049/jimmunol.1100272
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Gaublomme JT,
    2. Yosef N,
    3. Lee Y,
    4. Gertner RS,
    5. Yang LV,
    6. Wu C,
    7. Pandolfi PP,
    8. Mak T,
    9. Satija R,
    10. Shalek AK, et al.
    (2015) Single-cell genomics unveils critical regulators of Th17 cell pathogenicity. Cell 163: 1400–1412. doi:10.1016/j.cell.2015.11.009
    OpenUrlCrossRefPubMed
    1. Geiger R,
    2. Rieckmann JC,
    3. Wolf T,
    4. Basso C,
    5. Feng Y,
    6. Fuhrer T,
    7. Kogadeeva M,
    8. Picotti P,
    9. Meissner F,
    10. Mann M, et al.
    (2016) L-arginine modulates T cell metabolism and enhances survival and anti-tumor activity. Cell 167: 829–842.e13. doi:10.1016/j.cell.2016.09.031
    OpenUrlCrossRefPubMed
  20. ↵
    1. Geuking MB,
    2. Cahenzli J,
    3. Lawson MA,
    4. Ng DC,
    5. Slack E,
    6. Hapfelmeier S,
    7. McCoy KD,
    8. Macpherson AJ
    (2011) Intestinal bacterial colonization induces mutualistic regulatory T cell responses. Immunity 34: 794–806. doi:10.1016/j.immuni.2011.03.021
    OpenUrlCrossRefPubMed
  21. ↵
    1. Guo X,
    2. Zhang Y,
    3. Zheng L,
    4. Zheng C,
    5. Song J,
    6. Zhang Q,
    7. Kang B,
    8. Liu Z,
    9. Jin L,
    10. Xing R, et al.
    (2018) Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat Med 24: 978–985. doi:10.1038/s41591-018-0045-3
    OpenUrlCrossRef
  22. ↵
    1. Hamilton G,
    2. Colbert JD,
    3. Schuettelkopf AW,
    4. Watts C
    (2008) Cystatin F is a cathepsin C-directed protease inhibitor regulated by proteolysis. EMBO J 27: 499–508. doi:10.1038/sj.emboj.7601979
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Hartemann A,
    2. Bensimon G,
    3. Payan CA,
    4. Jacqueminet S,
    5. Bourron O,
    6. Nicolas N,
    7. Fonfrede M,
    8. Rosenzwajg M,
    9. Bernard C,
    10. Klatzmann D
    (2013) Low-dose interleukin 2 in patients with type 1 diabetes: A phase 1/2 randomised, double-blind, placebo-controlled trial. Lancet Diabetes Endocrinol 1: 295–305. doi:10.1016/s2213-8587(13)70113-x
    OpenUrlCrossRef
  24. ↵
    1. Howie D,
    2. Nolan KF,
    3. Daley S,
    4. Butterfield E,
    5. Adams E,
    6. Garcia-Rueda H,
    7. Thompson C,
    8. Saunders NJ,
    9. Cobbold SP,
    10. Tone Y, et al.
    (2009) MS4A4B is a GITR-associated membrane adapter, expressed by regulatory T cells, which modulates T cell activation. J Immunol 183: 4197–4204. doi:10.4049/jimmunol.0901070
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. International Multiple Sclerosis Genetics Consortium (IMSGC)
    (2008) Refining genetic associations in multiple sclerosis. Lancet Neurol 7: 567–569. doi:10.1016/S1474-4422(08)70122-4
    OpenUrlCrossRefPubMed
  26. ↵
    1. International Multiple Sclerosis Genetics Consortium,
    2. Hafler DA,
    3. Compston A,
    4. Sawcer S,
    5. Lander ES,
    6. Daly MJ,
    7. De Jager PL,
    8. de Bakker PI,
    9. Gabriel SB,
    10. Mirel DB, et al.
    (2007) Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med 357: 851–862. doi:10.1056/NEJMoa073493
    OpenUrlCrossRefPubMed
  27. ↵
    1. Iwata M,
    2. Hirakiyama A,
    3. Eshima Y,
    4. Kagechika H,
    5. Kato C,
    6. Song SY
    (2004) Retinoic acid imprints gut-homing specificity on T cells. Immunity 21: 527–538. doi:10.1016/j.immuni.2004.08.011
    OpenUrlCrossRefPubMed
  28. ↵
    1. Kamalipour H,
    2. Ahmadi S,
    3. Kamali K,
    4. Moaref A,
    5. Shafa M,
    6. Kamalipour P
    (2016) Ultrasound for localization of central venous catheter: A good alternative to chest X-ray? Anesth Pain Med 6: e38834. doi:10.5812/aapm.38834
    OpenUrlCrossRef
  29. ↵
    1. Kennedy-Nasser AA,
    2. Ku S,
    3. Castillo-Caro P,
    4. Hazrat Y,
    5. Wu MF,
    6. Liu H,
    7. Melenhorst J,
    8. Barrett AJ,
    9. Ito S,
    10. Foster A, et al.
    (2014) Ultra low-dose IL-2 for GVHD prophylaxis after allogeneic hematopoietic stem cell transplantation mediates expansion of regulatory T cells without diminishing antiviral and antileukemic activity. Clin Cancer Res 20: 2215–2225. doi:10.1158/1078-0432.ccr-13-3205
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Kim HJ,
    2. Barnitz RA,
    3. Kreslavsky T,
    4. Brown FD,
    5. Moffett H,
    6. Lemieux ME,
    7. Kaygusuz Y,
    8. Meissner T,
    9. Holderried TA,
    10. Chan S, et al.
    (2015) Stable inhibitory activity of regulatory T cells requires the transcription factor Helios. Science 350: 334–339. doi:10.1126/science.aad0616
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Kitagawa Y,
    2. Ohkura N,
    3. Kidani Y,
    4. Vandenbon A,
    5. Hirota K,
    6. Kawakami R,
    7. Yasuda K,
    8. Motooka D,
    9. Nakamura S,
    10. Kondo M, et al.
    (2017) Guidance of regulatory T cell development by Satb1-dependent super-enhancer establishment. Nat Immunol 18: 173–183. doi:10.1038/ni.3646
    OpenUrlCrossRefPubMed
  32. ↵
    1. Klatzmann D,
    2. Abbas AK
    (2015) The promise of low-dose interleukin-2 therapy for autoimmune and inflammatory diseases. Nat Rev Immunol 15: 283–294. doi:10.1038/nri3823
    OpenUrlCrossRefPubMed
  33. ↵
    1. Klein C,
    2. Waldhauer I,
    3. Nicolini VG,
    4. Freimoser-Grundschober A,
    5. Nayak T,
    6. Vugts DJ,
    7. Dunn C,
    8. Bolijn M,
    9. Benz J,
    10. Stihle M, et al.
    (2017) Cergutuzumab amunaleukin (CEA-IL2v), a CEA-targeted IL-2 variant-based immunocytokine for combination cancer immunotherapy: Overcoming limitations of aldesleukin and conventional IL-2-based immunocytokines. Oncoimmunology 6: e1277306. doi:10.1080/2162402x.2016.1277306
    OpenUrlCrossRef
  34. ↵
    1. Kobie JJ,
    2. Shah PR,
    3. Yang L,
    4. Rebhahn JA,
    5. Fowell DJ,
    6. Mosmann TR
    (2006) T regulatory and primed uncommitted CD4 T cells express CD73, which suppresses effector CD4 T cells by converting 5′-adenosine monophosphate to adenosine. J Immunol 177: 6780–6786. doi:10.4049/jimmunol.177.10.6780
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Koch MA,
    2. Tucker-Heard G,
    3. Perdue NR,
    4. Killebrew JR,
    5. Urdahl KB,
    6. Campbell DJ
    (2009) The transcription factor T-bet controls regulatory T cell homeostasis and function during type 1 inflammation. Nat Immunol 10: 595–602. doi:10.1038/ni.1731
    OpenUrlCrossRefPubMed
  36. ↵
    1. Koreth J,
    2. Matsuoka K,
    3. Kim HT,
    4. McDonough SM,
    5. Bindra B,
    6. Alyea EP 3rd.,
    7. Armand P,
    8. Cutler C,
    9. Ho VT,
    10. Treister NS, et al.
    (2011) Interleukin-2 and regulatory T cells in graft-versus-host disease. N Engl J Med 365: 2055–2066. doi:10.1056/nejmoa1108188
    OpenUrlCrossRefPubMed
  37. ↵
    1. Lathrop SK,
    2. Bloom SM,
    3. Rao SM,
    4. Nutsch K,
    5. Lio CW,
    6. Santacruz N,
    7. Peterson DA,
    8. Stappenbeck TS,
    9. Hsieh CS
    (2011) Peripheral education of the immune system by colonic commensal microbiota. Nature 478: 250–254. doi:10.1038/nature10434
    OpenUrlCrossRefPubMed
  38. ↵
    1. Lefranc MP,
    2. Giudicelli V,
    3. Ginestoux C,
    4. Jabado-Michaloud J,
    5. Folch G,
    6. Bellahcene F,
    7. Wu Y,
    8. Gemrot E,
    9. Brochet X,
    10. Lane J, et al.
    (2009) IMGT, the international ImMunoGeneTics information system. Nucleic Acids Res 37: D1006–D1012. doi:10.1093/nar/gkn838
    OpenUrlCrossRefPubMed
  39. ↵
    1. Letourneau S,
    2. Krieg C,
    3. Pantaleo G,
    4. Boyman O
    (2009) IL-2- and CD25-dependent immunoregulatory mechanisms in the homeostasis of T-cell subsets. J Allergy Clin Immunol 123: 758–762. doi:10.1016/j.jaci.2009.02.011
    OpenUrlCrossRefPubMed
  40. ↵
    1. Liao G,
    2. Nayak S,
    3. Regueiro JR,
    4. Berger SB,
    5. Detre C,
    6. Romero X,
    7. de Waal Malefyt R,
    8. Chatila TA,
    9. Herzog RW,
    10. Terhorst C
    (2010) GITR engagement preferentially enhances proliferation of functionally competent CD4+CD25+FoxP3+ regulatory T cells. Int Immunol 22: 259–270. doi:10.1093/intimm/dxq001
    OpenUrlCrossRefPubMed
  41. ↵
    1. Lim HW,
    2. Hillsamer P,
    3. Kim CH
    (2004) Regulatory T cells can migrate to follicles upon T cell activation and suppress GC-Th cells and GC-Th cell-driven B cell responses. J Clin Invest 114: 1640–1649. doi:10.1172/jci200422325
    OpenUrlCrossRefPubMed
  42. ↵
    1. Loebbermann J,
    2. Thornton H,
    3. Durant L,
    4. Sparwasser T,
    5. Webster KE,
    6. Sprent J,
    7. Culley FJ,
    8. Johansson C,
    9. Openshaw PJ
    (2012) Regulatory T cells expressing granzyme B play a critical role in controlling lung inflammation during acute viral infection. Mucosal Immunol 5: 161–172. doi:10.1038/mi.2011.62
    OpenUrlCrossRefPubMed
  43. ↵
    1. Long SA,
    2. Cerosaletti K,
    3. Wan JY,
    4. Ho JC,
    5. Tatum M,
    6. Wei S,
    7. Shilling HG,
    8. Buckner JH
    (2011) An autoimmune-associated variant in PTPN2 reveals an impairment of IL-2R signaling in CD4(+) T cells. Genes Immun 12: 116–125. doi:10.1038/gene.2010.54
    OpenUrlCrossRefPubMed
  44. ↵
    1. Mahmud SA,
    2. Manlove LS,
    3. Schmitz HM,
    4. Xing Y,
    5. Wang Y,
    6. Owen DL,
    7. Schenkel JM,
    8. Boomer JS,
    9. Green JM,
    10. Yagita H, et al.
    (2014) Costimulation via the tumor-necrosis factor receptor superfamily couples TCR signal strength to the thymic differentiation of regulatory T cells. Nat Immunol 15: 473–481. doi:10.1038/ni.2849
    OpenUrlCrossRefPubMed
  45. ↵
    1. Malek TR
    (2008) The biology of interleukin-2. Annu Rev Immunol 26: 453–479. doi:10.1146/annurev.immunol.26.021607.090357
    OpenUrlCrossRefPubMed
  46. ↵
    1. Matsuoka K,
    2. Koreth J,
    3. Kim HT,
    4. Bascug G,
    5. McDonough S,
    6. Kawano Y,
    7. Murase K,
    8. Cutler C,
    9. Ho VT,
    10. Alyea EP, et al.
    (2013) Low-dose interleukin-2 therapy restores regulatory T cell homeostasis in patients with chronic graft-versus-host disease. Sci Transl Med 5: 179ra143. doi:10.1126/scitranslmed.3005265
    OpenUrlCrossRef
  47. ↵
    1. Matta BM,
    2. Lott JM,
    3. Mathews LR,
    4. Liu Q,
    5. Rosborough BR,
    6. Blazar BR,
    7. Turnquist HR
    (2014) IL-33 is an unconventional alarmin that stimulates IL-2 secretion by dendritic cells to selectively expand IL-33R/ST2+ regulatory T cells. J Immunol 193: 4010–4020. doi:10.4049/jimmunol.1400481
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Menning A,
    2. Hopken UE,
    3. Siegmund K,
    4. Lipp M,
    5. Hamann A,
    6. Huehn J
    (2007) Distinctive role of CCR7 in migration and functional activity of naive- and effector/memory-like Treg subsets. Eur J Immunol 37: 1575–1583. doi:10.1002/eji.200737201
    OpenUrlCrossRefPubMed
  49. ↵
    1. Miragaia RJ,
    2. Gomes T,
    3. Chomka A,
    4. Jardine L,
    5. Riedel A,
    6. Hegazy AN,
    7. Whibley N,
    8. Tucci A,
    9. Chen X,
    10. Lindeman I, et al.
    (2019) Single-cell transcriptomics of regulatory T cells reveals trajectories of tissue adaptation. Immunity 50: 493–504.e7. doi:10.1016/j.immuni.2019.01.001
    OpenUrlCrossRef
  50. ↵
    1. Moulton VR,
    2. Tsokos GC
    (2015) T cell signaling abnormalities contribute to aberrant immune cell function and autoimmunity. J Clin Invest 125: 2220–2227. doi:10.1172/jci78087
    OpenUrlCrossRefPubMed
  51. ↵
    1. Onengut-Gumuscu S,
    2. Chen WM,
    3. Burren O,
    4. Cooper NJ,
    5. Quinlan AR,
    6. Mychaleckyj JC,
    7. Farber E,
    8. Bonnie JK,
    9. Szpak M,
    10. Schofield E, et al.
    (2015) Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nat Genet 47: 381–386. doi:10.1038/ng.3245
    OpenUrlCrossRefPubMed
  52. ↵
    1. Pabbisetty SK,
    2. Rabacal W,
    3. Volanakis EJ,
    4. Parekh VV,
    5. Olivares-Villagomez D,
    6. Cendron D,
    7. Boyd KL,
    8. Van Kaer L,
    9. Sebzda E
    (2016) Peripheral tolerance can be modified by altering KLF2-regulated Treg migration. Proc Natl Acad Sci U S A 113: E4662–E4670. doi:10.1073/pnas.1605849113
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Peterson LB,
    2. Bell CJM,
    3. Howlett SK,
    4. Pekalski ML,
    5. Brady K,
    6. Hinton H,
    7. Sauter D,
    8. Todd JA,
    9. Umana P,
    10. Ast O, et al.
    (2018) A long-lived IL-2 mutein that selectively activates and expands regulatory T cells as a therapy for autoimmune disease. J Autoimmun 95: 1–14. doi:10.1016/j.jaut.2018.10.017
    OpenUrlCrossRef
  54. ↵
    1. Regateiro FS,
    2. Chen Y,
    3. Kendal AR,
    4. Hilbrands R,
    5. Adams E,
    6. Cobbold SP,
    7. Ma J,
    8. Andersen KG,
    9. Betz AG,
    10. Zhang M, et al.
    (2012) Foxp3 expression is required for the induction of therapeutic tissue tolerance. J Immunol 189: 3947–3956. doi:10.4049/jimmunol.1200449
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Rubtsov YP,
    2. Rasmussen JP,
    3. Chi EY,
    4. Fontenot J,
    5. Castelli L,
    6. Ye X,
    7. Treuting P,
    8. Siewe L,
    9. Roers A,
    10. Henderson WR, Jr., et al.
    (2008) Regulatory T cell-derived interleukin-10 limits inflammation at environmental interfaces. Immunity 28: 546–558. doi:10.1016/j.immuni.2008.02.017
    OpenUrlCrossRefPubMed
  56. ↵
    1. Saadoun D,
    2. Rosenzwajg M,
    3. Joly F,
    4. Six A,
    5. Carrat F,
    6. Thibault V,
    7. Sene D,
    8. Cacoub P,
    9. Klatzmann D
    (2011) Regulatory T-cell responses to low-dose interleukin-2 in HCV-induced vasculitis. N Engl J Med 365: 2067–2077. doi:10.1056/nejmoa1105143
    OpenUrlCrossRefPubMed
  57. ↵
    1. Sather BD,
    2. Treuting P,
    3. Perdue N,
    4. Miazgowicz M,
    5. Fontenot JD,
    6. Rudensky AY,
    7. Campbell DJ
    (2007) Altering the distribution of Foxp3(+) regulatory T cells results in tissue-specific inflammatory disease. J Exp Med 204: 1335–1347. doi:10.1084/jem.20070081
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Satija R,
    2. Farrell JA,
    3. Gennert D,
    4. Schier AF,
    5. Regev A
    (2015) Spatial reconstruction of single-cell gene expression data. Nat Biotechnol 33: 495–502. doi:10.1038/nbt.3192
    OpenUrlCrossRefPubMed
  59. ↵
    1. Schiering C,
    2. Krausgruber T,
    3. Chomka A,
    4. Frohlich A,
    5. Adelmann K,
    6. Wohlfert EA,
    7. Pott J,
    8. Griseri T,
    9. Bollrath J,
    10. Hegazy AN, et al.
    (2014) The alarmin IL-33 promotes regulatory T-cell function in the intestine. Nature 513: 564–568. doi:10.1038/nature13577
    OpenUrlCrossRefPubMed
  60. ↵
    1. Sekiya T,
    2. Kashiwagi I,
    3. Inoue N,
    4. Morita R,
    5. Hori S,
    6. Waldmann H,
    7. Rudensky AY,
    8. Ichinose H,
    9. Metzger D,
    10. Chambon P, et al.
    (2011) The nuclear orphan receptor Nr4a2 induces Foxp3 and regulates differentiation of CD4+ T cells. Nat Commun 2: 269. doi:10.1038/ncomms1272
    OpenUrlCrossRefPubMed
  61. ↵
    1. Siede J,
    2. Frohlich A,
    3. Datsi A,
    4. Hegazy AN,
    5. Varga DV,
    6. Holecska V,
    7. Saito H,
    8. Nakae S,
    9. Lohning M
    (2016) IL-33 receptor-expressing regulatory T cells are highly activated, Th2 biased and suppress CD4 T cell proliferation through IL-10 and TGFbeta release. PLoS One 11: e0161507. doi:10.1371/journal.pone.0161507
    OpenUrlCrossRef
  62. ↵
    1. Spangler JB,
    2. Trotta E,
    3. Tomala J,
    4. Peck A,
    5. Young TA,
    6. Savvides CS,
    7. Silveria S,
    8. Votavova P,
    9. Salafsky J,
    10. Pande VS, et al.
    (2018) Engineering a single-agent cytokine/antibody fusion that selectively expands regulatory T cells for autoimmune disease therapy. J Immunol 201: 2094–2106. doi:10.4049/jimmunol.1800578
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Tao MH,
    2. Morrison SL
    (1989) Studies of aglycosylated chimeric mouse-human IgG. Role of carbohydrate in the structure and effector functions mediated by the human IgG constant region. J Immunol 143: 2595–2601.
    OpenUrlAbstract
  64. ↵
    1. Velaga S,
    2. Ukena SN,
    3. Dringenberg U,
    4. Alter C,
    5. Pardo J,
    6. Kershaw O,
    7. Franzke A
    (2015) Granzyme A is required for regulatory T-cell mediated prevention of gastrointestinal graft-versus-host disease. PLoS One 10: e0124927. doi:10.1371/journal.pone.0124927
    OpenUrlCrossRef
  65. ↵
    1. Vella A,
    2. Cooper JD,
    3. Lowe CE,
    4. Walker N,
    5. Nutland S,
    6. Widmer B,
    7. Jones R,
    8. Ring SM,
    9. McArdle W,
    10. Pembrey ME, et al.
    (2005) Localization of a type 1 diabetes locus in the IL2RA/CD25 region by use of tag single-nucleotide polymorphisms. Am J Hum Genet 76: 773–779. doi:10.1086/429843
    OpenUrlCrossRefPubMed
  66. ↵
    1. von Spee-Mayer C,
    2. Siegert E,
    3. Abdirama D,
    4. Rose A,
    5. Klaus A,
    6. Alexander T,
    7. Enghard P,
    8. Sawitzki B,
    9. Hiepe F,
    10. Radbruch A, et al.
    (2016) Low-dose interleukin-2 selectively corrects regulatory T cell defects in patients with systemic lupus erythematosus. Ann Rheum Dis 75: 1407–1415. doi:10.1136/annrheumdis-2015-207776
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Wachstein J,
    2. Tischer S,
    3. Figueiredo C,
    4. Limbourg A,
    5. Falk C,
    6. Immenschuh S,
    7. Blasczyk R,
    8. Eiz-Vesper B
    (2012) HSP70 enhances immunosuppressive function of CD4(+)CD25(+)FoxP3(+) T regulatory cells and cytotoxicity in CD4(+)CD25(−) T cells. PLoS One 7: e51747. doi:10.1371/journal.pone.0051747
    OpenUrlCrossRefPubMed
  68. ↵
    1. Wallace C,
    2. Cutler AJ,
    3. Pontikos N,
    4. Pekalski ML,
    5. Burren OS,
    6. Cooper JD,
    7. Garcia AR,
    8. Ferreira RC,
    9. Guo H,
    10. Walker NM, et al.
    (2015) Dissection of a complex disease susceptibility region using a bayesian stochastic search approach to fine mapping. PLoS Genet 11: e1005272. doi:10.1371/journal.pgen.1005272
    OpenUrlCrossRefPubMed
  69. ↵
    1. Webster KE,
    2. Walters S,
    3. Kohler RE,
    4. Mrkvan T,
    5. Boyman O,
    6. Surh CD,
    7. Grey ST,
    8. Sprent J
    (2009) In vivo expansion of T reg cells with IL-2-mAb complexes: Induction of resistance to EAE and long-term acceptance of islet allografts without immunosuppression. J Exp Med 206: 751–760. doi:10.1084/jem.20082824
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Wheaton JD,
    2. Yeh CH,
    3. Ciofani M
    (2017) Cutting edge: c-Maf is required for regulatory T cells to adopt RORgammat(+) and follicular phenotypes. J Immunol 199: 3931–3936. doi:10.4049/jimmunol.1701134
    OpenUrlAbstract/FREE Full Text
  71. ↵
    1. Wing JB,
    2. Kitagawa Y,
    3. Locci M,
    4. Hume H,
    5. Tay C,
    6. Morita T,
    7. Kidani Y,
    8. Matsuda K,
    9. Inoue T,
    10. Kurosaki T, et al.
    (2017) A distinct subpopulation of CD25(-) T-follicular regulatory cells localizes in the germinal centers. Proc Natl Acad Sci U S A 114: E6400–E6409. doi:10.1073/pnas.1705551114
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Wolf FA,
    2. Angerer P,
    3. Theis FJ
    (2018) SCANPY: Large-scale single-cell gene expression data analysis. Genome Biol 19: 15. doi:10.1186/s13059-017-1382-0
    OpenUrlCrossRefPubMed
  73. ↵
    1. Wu K,
    2. Ma J,
    3. Bai W,
    4. Cui X,
    5. Han T,
    6. Wang S,
    7. Xie Y,
    8. Xie Y
    (2016) Short-term intratracheal use of PEG-modified IL-2 and glucocorticoid persistently alleviates asthma in a mouse model. Sci Rep 6: 31562. doi:10.1038/srep31562
    OpenUrlCrossRef
  74. ↵
    1. Xing S,
    2. Gai K,
    3. Li X,
    4. Shao P,
    5. Zeng Z,
    6. Zhao X,
    7. Zhao X,
    8. Chen X,
    9. Paradee WJ,
    10. Meyerholz DK, et al.
    (2019) Tcf1 and Lef1 are required for the immunosuppressive function of regulatory T cells. J Exp Med 216: 847–866. doi:10.1084/jem.20182010
    OpenUrlAbstract/FREE Full Text
  75. ↵
    1. Xu H,
    2. Yan Y,
    3. Williams MS,
    4. Carey GB,
    5. Yang J,
    6. Li H,
    7. Zhang GX,
    8. Rostami A
    (2010) MS4a4B, a CD20 homologue in T cells, inhibits T cell propagation by modulation of cell cycle. PLoS One 5: e13780. doi:10.1371/journal.pone.0013780
    OpenUrlCrossRefPubMed
  76. ↵
    1. Yadav M,
    2. Louvet C,
    3. Davini D,
    4. Gardner JM,
    5. Martinez-Llordella M,
    6. Bailey-Bucktrout S,
    7. Anthony BA,
    8. Sverdrup FM,
    9. Head R,
    10. Kuster DJ, et al.
    (2012) Neuropilin-1 distinguishes natural and inducible regulatory T cells among regulatory T cell subsets in vivo. J Exp Med 209: 1713–1722, S1711-1719. doi:10.1084/jem.20120822
    OpenUrlAbstract/FREE Full Text
  77. ↵
    1. Yamaguchi T,
    2. Hirota K,
    3. Nagahama K,
    4. Ohkawa K,
    5. Takahashi T,
    6. Nomura T,
    7. Sakaguchi S
    (2007) Control of immune responses by antigen-specific regulatory T cells expressing the folate receptor. Immunity 27: 145–159. doi:10.1016/j.immuni.2007.04.017
    OpenUrlCrossRefPubMed
  78. ↵
    1. Yamanouchi J,
    2. Rainbow D,
    3. Serra P,
    4. Howlett S,
    5. Hunter K,
    6. Garner VE,
    7. Gonzalez-Munoz A,
    8. Clark J,
    9. Veijola R,
    10. Cubbon R, et al.
    (2007) Interleukin-2 gene variation impairs regulatory T cell function and causes autoimmunity. Nat Genet 39: 329–337. doi:10.1038/ng1958
    OpenUrlCrossRefPubMed
  79. ↵
    1. Yang JH,
    2. Cutler AJ,
    3. Ferreira RC,
    4. Reading JL,
    5. Cooper NJ,
    6. Wallace C,
    7. Clarke P,
    8. Smyth DJ,
    9. Boyce CS,
    10. Gao GJ, et al.
    (2015) Natural variation in interleukin-2 sensitivity influences regulatory T-cell frequency and function in individuals with long-standing type 1 diabetes. Diabetes 64: 3891–3902. doi:10.2337/db15-0516
    OpenUrlAbstract/FREE Full Text
  80. ↵
    1. Yu A,
    2. Snowhite I,
    3. Vendrame F,
    4. Rosenzwajg M,
    5. Klatzmann D,
    6. Pugliese A,
    7. Malek TR
    (2015) Selective IL-2 responsiveness of regulatory T cells through multiple intrinsic mechanisms supports the use of low-dose IL-2 therapy in type 1 diabetes. Diabetes 64: 2172–2183. doi:10.2337/db14-1322
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Zemmour D,
    2. Zilionis R,
    3. Kiner E,
    4. Klein AM,
    5. Mathis D,
    6. Benoist C
    (2018) Single-cell gene expression reveals a landscape of regulatory T cell phenotypes shaped by the TCR. Nat Immunol 19: 291–301. doi:10.1038/s41590-018-0051-0
    OpenUrlCrossRef
  82. ↵
    1. Zheng SG,
    2. Wang JH,
    3. Stohl W,
    4. Kim KS,
    5. Gray JD,
    6. Horwitz DA
    (2006) TGF-beta requires CTLA-4 early after T cell activation to induce FoxP3 and generate adaptive CD4+CD25+ regulatory cells. J Immunol 176: 3321–3329. doi:10.4049/jimmunol.176.6.3321
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Zurawski SM,
    2. Vega F, Jr.,
    3. Doyle EL,
    4. Huyghe B,
    5. Flaherty K,
    6. McKay DB,
    7. Zurawski G
    (1993) Definition and spatial location of mouse interleukin-2 residues that interact with its heterotrimeric receptor. EMBO J 12: 5113–5119. doi:10.1002/j.1460-2075.1993.tb06206.x
    OpenUrlCrossRefPubMed
  84. ↵
    1. Zurawski SM,
    2. Zurawski G
    (1989) Mouse interleukin-2 structure-function studies: Substitutions in the first alpha-helix can specifically inactivate p70 receptor binding and mutations in the fifth alpha-helix can specifically inactivate p55 receptor binding. EMBO J 8: 2583–2590. doi:10.1002/j.1460-2075.1989.tb08397.x
    OpenUrlCrossRefPubMed
View Abstract
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Dynamic changes in the regulatory T-cell heterogeneity and function by murine IL-2 mutein
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Single-cell RNA-seq analysis of regulatory T cells
Daniel R Lu, Hao Wu, Ian Driver, Sarah Ingersoll, Sue Sohn, Songli Wang, Chi-Ming Li, Hyewon Phee
Life Science Alliance Apr 2020, 3 (5) e201900520; DOI: 10.26508/lsa.201900520

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Single-cell RNA-seq analysis of regulatory T cells
Daniel R Lu, Hao Wu, Ian Driver, Sarah Ingersoll, Sue Sohn, Songli Wang, Chi-Ming Li, Hyewon Phee
Life Science Alliance Apr 2020, 3 (5) e201900520; DOI: 10.26508/lsa.201900520
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 3, No. 5
May 2020
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Data Availability
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Medical Research
  • Immunology

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2021 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.