1. Computational and Systems Biology
  2. Immunology and Inflammation

Adequate immune response ensured by binary IL-2 and graded CD25 expression in a murine transfer model

  1. Franziska Fuhrmann
  2. Timo Lischke
  3. Fridolin Gross
  4. Tobias Scheel
  5. Laura Bauer
  6. Khalid Wasim Kalim
  7. Andreas Radbruch
  8. Hanspeter Herzel
  9. Andreas Hutloff
  10. Ria Baumgrass  Is a corresponding author
  1. Robert Koch Institute, Germany
  2. A Leibniz Institute, Germany
  3. Charité University Medicine, Germany
https://doi.org/10.7554/eLife.20616
Share this article
Cite this article
  1. Franziska Fuhrmann
  2. Timo Lischke
  3. Fridolin Gross
  4. Tobias Scheel
  5. Laura Bauer
  6. Khalid Wasim Kalim
  7. Andreas Radbruch
  8. Hanspeter Herzel
  9. Andreas Hutloff
  10. Ria Baumgrass
(2016)
Adequate immune response ensured by binary IL-2 and graded CD25 expression in a murine transfer model
eLife 5:e20616.
https://doi.org/10.7554/eLife.20616
  2. Download BibTeX
  3. Download .RIS

Abstract

The IL-2/IL-2Ralpha (CD25) axis is of central importance for the interplay of effector and regulatory T cells. Nevertheless, the question how different antigen loads are translated into appropriate IL-2 production to ensure adequate responses against pathogens remains largely unexplored. Here we find that at single cell level, IL-2 is binary (digital) and CD25 is graded expressed whereas at population level both parameters show graded expression correlating with the antigen amount. Combining in vivo data with a mathematical model we demonstrate that only this binary IL-2 expression ensures a wide linear antigen response range for Teff and Treg cells under real spatiotemporal conditions. Furthermore, at low antigen concentrations binary IL-2 expression safeguards by its spatial distribution selective STAT5 activation only of closely adjacent Treg cells regardless of their antigen specificity. These data show that the mode of IL-2 secretion is critical to tailor the adaptive immune response to the antigen amount.

https://doi.org/10.7554/eLife.20616.001

Introduction

The adaptive immune response of healthy vertebrates can be adjusted to match the overall severity of a pathogen. At the population level, T cells translate a total pathogen input into an adequate cytokine response over a large dynamic range. Individual T cells sense both antigen affinity and frequency and tailor the collective antigen-specific T cell activation to the antigen amount (Allison et al., 2016; Gottschalk et al., 2012; Kellogg et al., 2015; Miskov-Zivanov et al., 2013). The cytokine IL-2 and its signaling are of central importance for these processes. IL-2 is mainly produced by activated Foxp3− CD4+ T helper cells (hereinafter referred to as 'Th cells') and mediates diverse pleiotropic actions such as promotion of Th1 and Th2 differentiation, expansion of multiple Th cell populations, as well as generation and maintenance of Foxp3+ regulatory T cells (Treg cells) (Liao et al., 2013). Thus, IL-2 positively acts on Th cells through synergistic effects between antigen and cytokines and negatively through preferred binding to and activation of Foxp3+ Treg cells.

There are a number of in vitro studies, some of them even combined with mathematical modeling, to dissect the dichotomic effect of IL-2 signaling in effector Th cells and Foxp3+ Treg cells (Benary et al., 2008; Bendfeldt et al., 2012; Busse et al., 2010; Feinerman et al., 2010; Tkach et al., 2014; Voisinne et al., 2015). Observations from various knock-out mice (e.g. IL-2Rα, STAT5) and other in vivo studies, however, have challenged some conclusions of these in vitro data and underlined the importance of local antigen concentrations, temporal aspects, cytokine compositions, cell heterogeneity in tissues and much more for the biological outcome orchestrated by IL-2 during an immune response (reviewed in Liao et al., 2013).

There are several articles dealing with certain molecular features of the IL-2 pathway in vivo, such as: First, Sojka et al. studied the kinetics of IL-2 secretion by naive and memory Th cells (Sojka et al., 2004) and some studies observed a paracrine IL-2 action on STAT5 activation (Long and Adler, 2006; O'Gorman et al., 2009; Sabatos et al., 2008). Second, Amado et al. proposed that IL-2 is sensed by both the activated Th cell pool and by Treg cells, which reciprocally regulate cells of the IL-2-producing cell subset in IL-2 reporter mice (Amado et al., 2013). Third, O’Gorman et al. showed that in the initial phase of antigen-specific activation IL-2-dependent STAT5 phosphorylation occurred primarily in Treg cells (O'Gorman et al., 2009). Fourth, Liu et al. underscored the importance of a fine-tuned local IL-2 availability to limit otherwise constant damage from activated auto-responsive Th cells under steady state conditions (Liu et al., 2015).

However, in all of these studies the antigen load-dependent IL-2 expression and IL-2 signaling via CD25/STAT5 in both antigen-specific as well as endogenous Th and Treg cells was not continuously monitored in vivo. Furthermore, there is the open question whether a graded antigen load is translated into binary (also known as digital, all-or-none, switch-like, or bimodal) or graded expression of TCR-dependent genes coding for important cytokines and receptors such as IL-2 and CD25 or CD69, respectively, in vivo. Even in vitro there were conflicting data about binary or graded expression of CD25 and CD69 (Allison et al., 2016; Busse et al., 2010; Feinerman et al., 2010; Huang et al., 2013; Podtschaske et al., 2007; Tkach et al., 2014). This is an interesting issue because binary responses are able to filter out noisy signals and can integrate different information to control cell responses. Thereby, binary responses allow parallel and independent control of single cell activation probability and population heterogeneity (Kellogg et al., 2015; Köck et al., 2014).

Here, we studied the antigen amount-driven regulation of IL-2 and CD25 expression at the single cell and population level as well as IL-2 signaling under the real spatial-temporal conditions of secondary lymphoid organs in vivo. Using a murine adoptive transfer model enabled us to study an antigen-specific T cell response in syngeneic recipient mice and cytokine secretion assays to measure the rate and magnitude of cytokine production of individual cells over time. In addition, mathematical modeling allowed us to compare the effects of the hypothesized graded and observed binary IL-2 expression for the first time. The model revealed advantages of the latter for a wide linear scaling of STAT5 activation in Treg and/or Th cells to antigen concentration. Finally, our results demonstrate that the combination of binary IL-2 secretion and graded CD25 expression at the cellular level enables the parallel control of single cell activation probability and population heterogeneity by the respective antigen amount.

Results

Binary IL-2 production of Th cells in vivo after immunization

Previously, we showed in vitro that IL-2 is produced in a binary (‘digital’ or ‘all-or-none’) manner by Th cells (Podtschaske et al., 2007), because variation of strength of the stimulation changed the number of IL-2 expressing cells rather than the intensity of IL-2 per cell. It is challenging to analyze whether this holds true in vivo, too.

To address this question, we adoptively transferred naïve transgenic Foxp3− CD4+ T helper cells from DO11.10 mice into recipient mice and analyzed IL-2 production in the spleen using a cytokine secretion assay (Assenmacher et al., 1998) 14 hr after immunization with ovalbumin (OVA). Different amounts of antigen used for immunization resulted in changes in the number of IL-2-producing cells (Figure 1A, left), but almost no change in the mean fluorescence intensity (MFI) of IL-2 expression per Th cell (Figure 1A, right). Thus, the amount of antigen correlates with the proportion of Th cells that produced IL-2, but does not change the amount of IL-2 per cell, which confirms that IL-2 production is regulated in a binary fashion in vivo, too. The IL-2 expression peaked between 6 and 18 hr (Figure 1D). The IL-2 secretion assay detected up to 75% IL-2-producing Th cells at the highest amount of antigen (2 mg) used. This number reflects the in vivo situation better than direct measurement of IL-2 by intracellular staining, which only detected maximally ~9% IL-2-producing Th cells, considered to be too low (Long and Adler, 2006). Moreover, we excluded an underestimation of cytokine production by confirming that the cytokine capture matrix is not saturated in our experiments (Figure 1—figure supplement 1). The data shown were generated with transfer of OVA-specific DO11.10 T cells (Figure 1), but were also reproduced using T cells from OT-II mice (Figure 1—figure supplement 2). Our results were confirmed in a more thorough statistical analysis in which we explicitly tested for bimodality in IL-2 expression and used mixed gaussian models to extract expression values and percentages of cell subpopulations (Figure 1—figure supplement 3).

Figure 1 with 6 supplements see all
Download asset Open asset
Binary IL-2 secretion of adoptively transferred Th cells in vivo after immunization with increasing amounts of antigen and kinetics of IL-2 secretion.

BALB/c mice, adoptively transferred with OVA-specific T cells from DO11.10 mice, were immunized intravenously with the indicated amounts of OVA and 10 μg LPS as adjuvant and analyzed 14 hr later (A–C) or with 2 mg OVA + 10 µg LPS and analyzed over time (D). IL-2 secretion of OVA-TCR+ CD4+ T cells was analyzed using the IL-2 secretion assay. Data from gated transgenic T cells (live B220− CD4+ OVA-TCR+ Foxp3−) were concatenated (3 mice per antigen dose). (A) IL-2 production was plotted against the amount of OVA used for immunization: Very left dot plot: fluorescence intensities of IL-2 per cell; left graph: frequencies of IL-2 producing cells; right histogram overlay: IL-2 fluorescence intensities (FI) for all gated OVA-TCR+ CD4+ T cells; very right histogram overlay: IL-2 producing OVA-TCR+ CD4+ T cells. (B) Similar analysis for the expression of CD25 on gated OVA-TCR+ CD4+ T cells. In the very right graph the mean fluorescence intensity (MFI) of CD25 is plotted against the amount of OVA. (C) Similar analysis for the expression of CD69 on gated OVA-TCR+ CD4+ T cells. CD69 expressing T cells were distinguished in antigen-specific activated CD69high and bystander-activated CD69intermediate expressing T cells. The gates for CD69intermediate and CD69high T cells were set according to unimmunized and non-transfer immunized controls. (D) IL-2 production (top row) and CD25 expression (bottom row) in OVA-TCR+ CD4+ T cells are plotted against time after immunization. Dot plots (left) depict fluorescence intensities per cell; the middle graphs depict frequencies of positive cells; the right graphs depict MFIs. Statistics: mean and standard deviation were plotted in all graphs. Data are representative of three independent experiments.

https://doi.org/10.7554/eLife.20616.002

In parallel we studied whether the amount of the transcription factors c-Fos and NFATc2 limit the IL-2 production per mouse cell in vivo, as already shown in vitro for human peripheral blood memory Th cells (Podtschaske et al., 2007). To this end we used adoptively transferred naïve transgenic CD4+ T cells 7 hr after immunization of mice and co-stained c-Fos, NFATc2, and IL-2 using intracellular staining and cytokine secretion assays, respectively. Interestingly, only the amount of c-Fos but not of NFATc2 per Th cell correlates with the probability of IL-2 expression in vivo. Quartiles with the same cell number but different fluorescence intensities of c-Fos or NFATc2 (Figure 1—figure supplement 4, left) showed that the probability of IL-2 production correlates with the c-Fos amount per cell but is independent of NFATc2 amounts within their physiological range (Figure 1—figure supplement 4, right). This is in line with the observation of Allison et al. that the MEK/ERK/AP1 axis plays an important role in translating TCR signaling into proportional activation of genes essential for T cell function (Allison et al., 2016).

Graded IL-2Rα expression by Th cells in vivo after immunization

Next we analyzed how CD25, the high affinity IL-2Rα chain, is expressed in vivo. So far, in vitro quantitative single cell studies have produced conflicting results proposing opposite manners of CD25 expression, namely binary (Busse et al., 2010) and graded (Feinerman et al., 2010). Given that CD25 expression peaks between 12 and 24 hr after immunization (Lischke et al., 2012; Long and Adler, 2006), we observed graded CD25 and CD69 expression of adoptively transferred transgenic CD4+ T cells 14 hr after immunization of mice. Different antigen concentrations used for immunization changed the mean fluorescence intensity (MFI) of CD25 expression per Th cell (Figure 1B, right). The amount of antigen correlated with the amount of CD25 expressed per cell (Figure 1B, left), and thus its ability to respond to IL-2. Notably, CD25 expression was up-regulated only by transferred antigen-specific Th and Treg cells (Figure 1—figure supplements 5A and 6), but not by endogenous T cells. This fact is pointing to preferential up-regulation of CD25 expression via TCR-stimulated T cell activation. The same holds true for high expression of CD69 (Figure 1—figure supplement 5B). However, a small bystander expression of CD69 (CD69intermediate cells) was observed on endogenous Th and Treg cells after immunization of mice with OVA + LPS.

Next, we analyzed the kinetics of IL-2 secretion and CD25 expression in adoptively transferred transgenic Th cells after immunization of mice with 2 mg OVA + 10 µg LPS (Figure 1D). Both processes reached their maximum value between 3 to 6 hr and obviously CD25 protein expression persisted longer than IL-2 secretion.

CD25 upregulation is similar in IL-2-producing and IL-2-nonproducing transgenic T cells in vivo

The amount of expression of CD25 is critical for Th cells in order to respond to IL-2, in particular if the IL-2 concentration is low (Busse et al., 2010; Feinerman et al., 2010; Popmihajlov and Smith, 2008; Tkach et al., 2014). In contrast to Treg cells, Th cells show high expression of CD25 only transiently, due to the short duration of both, antigen-induced TCR signaling and IL-2/IL-2R signaling (Malek and Castro, 2010; Tkach et al., 2014).

Interestingly, co-staining of CD25 and IL-2 in adoptively transferred transgenic Th cells (Figure 2A, top row) revealed that the amount and variation of CD25 expression is similar in IL-2-producing and IL-2-nonproducing cells 14 hr after immunization with different antigen concentrations (Figure 2A, bottom row). Accordingly, the kinetics of IL-2 secretion and CD25 expression in adoptively transferred transgenic Th cells after immunization of mice with 2 mg OVA + 10 µg LPS was similar in IL-2-producing and IL-2-nonproducing cells during the time of expression of both molecules (Figure 2B). These data suggest that mainly the amount of antigen and a paracrine, rather than an autocrine, action of IL-2 regulates the expression of CD25 in vivo.

Figure 2
Download asset Open asset
In vivo kinetics of IL-2 secretion and CD25 expression in adoptively transferred Th cells after immunization and their correlation upon titration of antigen amount as well as over time.

BALB/c mice, adoptively transferred with OVA-specific T cells from DO11.10 mice, were immunized intravenously with increasing amounts of OVA and 10 μg LPS as adjuvant and analyzed at 14 hr (A) or 2 mg OVA + 10 μg LPS and analyzed over time (B). IL-2 secretion of OVA-TCR+ CD4+ T cells was analyzed using the IL-2 secretion assay. Data from gated OVA-TCR+ CD4+ T cells (live B220− CD4+ OVA-TCR+ Foxp3−) were concatenated (3 mice per antigen dose or time point). (A) IL-2 production and CD25 expression (fluorescence intensities per cell) are depicted for increasing amounts of OVA used for immunization (top row, dot plots). A histogram overlay shows the IL-2 fluorescence intensities of gated OVA-TCR+ CD4+ T cells for different antigen amounts for direct comparison (bottom row, left). Histogram overlays of CD25 fluorescence intensities in IL-2 producing (red) and non-producing (black) OVA-TCR+ CD4+ T cells for individual OVA amounts are shown for direct comparison of CD25 expression (bottom row, right). (B) IL-2 production and CD25 expression (fluorescence intensities per cell) are depicted over time (top row, dot plots). A histogram overlay shows the IL-2 fluorescence intensities of gated OVA-TCR+ CD4+ T cells over time for direct comparison (bottom row, left). Histogram overlays of CD25 fluorescence intensities in IL-2 producing (red) and non-producing (black) OVA-TCR+ CD4+ T cells for individual time points are shown for direct comparison of CD25 expression (bottom row, right). Statistics: mean and standard deviation were plotted in all graphs. Data are representative of three independent experiments.

https://doi.org/10.7554/eLife.20616.009

Computational modeling suggests physiological advantages of binary IL-2 secretion

To get an idea about the physiological importance of binary IL-2 expression for activation of antigen-specific Th and Treg cells, we compared the observed binary IL-2 pattern to the hypothetical scenario of graded IL-2 secretion. The physiological consequences of different secretion patterns are not intuitively obvious and might depend on many factors such as the spatial distribution of cells in the spleen and on diffusion and secretion rates. Computational modeling allowed us to simulate the consequences of different hypothetical scenarios taking into account quantitative empirical information about the real system. All parameters of the model are listed in Supplementary file 1 and the spatial distribution of T cells is depicted in Figure 3—figure supplement 1.

Our computational model describes a population of T cells that can secrete IL-2 and/or express IL-2R at varying levels corresponding to the antigen stimulus. Cells are considered ‘activated’ when the amount of receptor-bound IL-2 exceeds a particular threshold. Figure 3A schematically represents the differences between the compared scenarios at low antigen concentrations. Intuitively, one might expect that activation of Th cells requires higher levels of antigen stimulus in the graded scenario because for IL-2 diffusion it is harder to overcome the critical distances between cells. It is not obvious, however, whether there are differences in the activation of Treg cells as they are able to sense very low concentrations of IL-2. The model can assist in making sense of or in correcting our intuitions.

Figure 3 with 1 supplement see all
Download asset Open asset
Model and experimental validation of binary versus graded IL-2 production.

(A) Schematic representation for the effects of different modes of IL-2 secretion (red) by Th (blue) on surrounding Th (blue) and Treg (green) cells. Graded IL-2 secretion at low antigen dose (left) is compared to binary IL-2 secretion (right). IL-2 mediated T cell activation is depicted via conversion of quiescent CD25 (black Y) to activated CD25 (red Y). (B) Model and comparison to experiment: Diffusion of IL-2 and interaction of IL-2 and IL-2R is modeled for a population of T cells in the spleen. The figure depicts the fraction of T cells that have more than a critical amount IL-2 bound to IL-2R on their surface (and thus reach their activation threshold) and compares it to the fraction of p-STAT5+ T cells in the experiment. (C) BALB/c mice, adoptively transferred with OVA-specific T cells from DO11.10 mice, were immunized intravenously with increasing amounts of OVA and 10 µg LPS as adjuvant. CD25 expression and STAT5 phosphorylation were analyzed after 14 hr on gated CD4+ T cells (live B220− CD4+ OVA-TCR+/− Foxp3+/−). Left figure: Dot plots of p-STAT5 versus CD25 expression for gated transgenic Th and Treg cells without immunization or after immunization with 2 mg OVA + 10 µg LPS. Right graph: Percentage of p-STAT5+ signal in different T cell subpopulations versus amount of OVA used for immunization (3 mice per antigen dose). Statistics: mean and standard deviation were plotted. Data are representative of two independent experiments.

https://doi.org/10.7554/eLife.20616.010

Our simulations show that in the binary scenario for IL-2 the number of activated Th and Treg cells slowly increases with increasing antigen stimulus, while the behaviour is rather switch-like in the graded case (Figure 3B). This can be understood in the following way: If secretion is graded (i.e. many producers), the distribution of IL-2 is relatively homogeneous across the whole volume. Consequently, all cells are exposed to similar levels of IL-2 and thus are either all activated or not in the graded IL-2 scenario. If secretion is binary (i.e. the number of secreting cells increases with increasing stimulus, the amount is constant), by contrast, then the number of activated cells increases more gradually because cells that are in close proximity to secreting cells are activated first. Binary secretion at the cellular level thus leads to graded activation at the population level and vice versa.

In summary, binary IL-2 and graded CD25 expression as well as higher CD25 expression by Treg cells ensures a fast action and wide linear antigen response range of Treg cells. Interestingly, the model simulations predict that Th cells are activated at lower antigen concentrations in the binary than in the graded scenario. However, Treg cells are always activated at lower antigen concentrations than Th cells. This is obviously important for effective control of autoimmunity ensuring Treg cell activation in spatial proximity even if only rare self-reactive IL-2 secreting Th cells exist (Liu et al., 2015).

Binary IL-2 secretion tailors the number of activated Th and Treg cells to the amount of antigen

To compare the model data with the in vivo situation, we analysed IL-2 signaling by STAT5 activation in antigen-specific T cells. We co-stained phosphorylated STAT5 and CD25 directly ex vivo in adoptively transferred antigen-specific Th and Treg cells 14 hr after immunization of mice. As expected, with increasing amounts of antigen (and hence increasing IL-2 secretion) STAT5 phosphorylation increased proportionally in antigen-specific Th and Treg cells. Moreover, antigen-specific Treg cells showed a much higher STAT5 phosphorylation than antigen-specific Th cells. However, the p-STAT5 signals were quite low. Even at the highest antigen concentration (2 mg OVA) only ~11% of Th cells stained positive for p-STAT5 despite the fact that about 75% of Th cells secrete IL-2 and about 95% express CD25 at this condition (Figure 3C). Within the Treg cell population about 50% were p-STAT5+ (Figure 3C). It is possible that p-STAT5+ frequencies are underestimated due to difficulties in preserving the p-STAT5 signal during the preparation of cells ex vivo and the staining process. O’ Gorman et al. described that the induction of p-STAT5 is lower in naïve Th cells compared to repeatedly activated cells but lasts longer (O'Gorman et al., 2009). In line with this, our endogenous antigen-unspecific Th cells had a very low expression of CD25 and almost no p-STAT5. In contrast, STAT5 activation (Figure 3C) was observed at a similar degree on both Treg cell populations, transgenic (antigen-specific) and endogenous Treg cells.

To analyse whether the low p-STAT5 levels of antigen-specific Th cells are due to their unresponsiveness to IL-2 or because they do not have access to IL-2 in vivo, we treated mice with either recombinant IL-2 (20 μg) or IL-2-complex (2.5 μg IL-2 with 10 μg anti-IL-2 antibody JES6-5H4). Furthermore, we incubated one part of the spleen in vitro with IL-2 (200 ng recombinant IL-2 for 25 min). The data clearly show that neither IL-2 treatment of mice nor in vitro IL-2 treatment of Th cells increased STAT5 activation in transgenic Th cells (Figure 4). Thus it is very probable that antigen-specific Th cells are unresponsive to IL-2 by negative crosstalk between antigen and IL-2 signaling as described before (Tkach et al., 2014; Villarino et al., 2007; Waysbort et al., 2013). Taken together, our model data are consistent with the experimental data (Figure 3C) and reveal the importance of IL-2-activated Treg cells for limiting the immune response of Th cells. First, at low amounts of antigen (<30 µg OVA) preferentially a proportion of Treg cells is activated by IL-2, and thus can limit the immune response at low danger. Second, binary IL-2 expression tailors the fraction of activated Treg cells proportionately to the antigen amount and avoids that all Treg cells are almost simultaneously activated as would be expected in the graded IL-2 expression scenario. Third, at low antigen amounts binary IL-2 expression ensures by its spatial distribution that only those antigen-specific and unspecific Treg cells are activated, which are in close spatial proximity to the few IL-2 producing Th cells. This is in good agreement with data of Liu et al. concerning observations on self tolerance (Liu et al., 2015).

Figure 4
Download asset Open asset
Treatment of mice with recombinant IL-2 does not increase STAT5 activation.

C57BL/6 mice, adoptively transferred with OVA-specific T cells from OT-II mice, were immunized intravenously with 2000 µg OVA and 10 µg LPS or were left unimmunized. After 12 hr mice received either 20 µg recombinant IL-2, 2.5 µg IL-2 complexed with 10 µg anti-IL-2 (JES6-5H4; targets IL-2 to CD122), or PBS as control. Mice were sacrificed 1 hr later. Cells from one part of the spleen were incubated for 25 min in vitro with 200 ng recombinant IL-2. (A) Representative original data (gated on transgenic CD4+ T cells) are shown in the flow cytometric plots. (B) Percentage of p-STAT5+ cells within either endogenous CD4+ T cells (including FoxP3+ cells) or transgenic FoxP3-negative CD4+ T cells are depicted in bar graphs. Pooled data are from two independent experiments (together six mice per group). Error bars represent SEM.

https://doi.org/10.7554/eLife.20616.012

The ratio of antigen-specific Th and Treg cells remains stable over a wide antigen concentration range

Next, we investigated whether the observed preferential Treg cell activation at low antigen doses is translated into a higher expansion or survival rate compared to Th cells. Using the adoptive transfer model, we found that antigen-specific Th and Treg cells did not proliferate after immunization at low antigen concentrations (2–8 µg OVA) in CFSE dilution assays (Figure 5A). In line with this, we observed similar total cell counts and relative frequencies as in the unimmunized controls under these conditions (Figure 5B,C). The percentage of dividing progenitors was approximately 10%, which appeared to be the background level already present in Th and Treg cells without immunization (Figure 5D). Above 8 µg OVA, total cell numbers of transgenic Th and Treg cells increased proportionally to the antigen amount, but with a dramatic different slope. CFSE dilution profiles revealed a stronger proliferative response of OVA-specific Th cells compared to Treg cells at higher concentrations. As an example, at 500 µg OVA 94% of Th cell progenitors were dividing compared to 65% of Treg cell progenitors (Figure 5D).

Figure 5 with 1 supplement see all
Download asset Open asset
Dependency of in vivo proliferation rate of adoptively transferred Th and Treg cells on amount of antigen.

BALB/c mice, adoptively transferred with CFSE-labeled OVA-specific T cells from DO11.10 mice, were immunized intravenously with increasing amounts of OVA and 10 μg LPS as adjuvant and analyzed at 72 hr. Gated OVA-TCR+ CD4+ T cells (live B220− CD4+ OVA-TCR+ Foxp3+/−) were analyzed for CFSE dilution (2 mice per antigen dose). (A) Gating strategy for proliferation analysis: after doublet exclusion in FSC-A/FSC-H and SSC-A/SSC-H plots (not shown), cells were gated according to scatter characteristics, dead cells (Pacific Orange positive) were excluded, followed by gating on B220− CD4+ OVA-TCR+ Foxp3+/− T cells and analysis of CFSE dilution. (B) The total T cell count per spleen was calculated and plotted against amount of antigen for OVA-TCR+ Foxp3− Th cells (blue circles) and OVA-TCR+ Foxp3+ Treg cells (green circles). (C) The relative frequencies of Foxp3− Th cells (blue circles) and Foxp3+ Treg cells (green circles) in all transgenic OVA-TCR+ T cells were plotted against amount of antigen. (D) Furthermore, the percentage of dividing cells of the progenitor cells (generation 0) was calculated and plotted against the amount of OVA used for immunization. Statistics: mean and standard deviation were plotted in all graphs. Data are representative of four independent experiments.

https://doi.org/10.7554/eLife.20616.013

Subsequently, we analyzed the apoptosis rate of transferred antigen-specific Th and Treg cells 72 hr after immunization using annexin V binding (a marker for loss of cell membrane integrity upon apoptosis) (Vermes et al., 1995) and FLICA staining (detecting activated caspases) (Thornberry et al., 1997). In both assays, antigen-specific Treg cells showed less apoptosis than antigen-specific Th cells, in particular at low antigen concentrations (4–5% and 7–9% apoptotic cells, respectively, at 2–8 µg OVA). Apoptosis of Treg cells was proportionally increasing with the antigen amount up to the level of Th cells at 2000 µg OVA (Figure 5—figure supplement 1). In contrast, the apoptosis rate of antigen-specific Th cells was almost constant over the covered antigen concentration range.

Comparing cell counts at low antigen concentrations (2–8 µg OVA) (left parts of Figure 5B and Figure 5—figure supplement 1), it is obvious that apoptosis does not play an important role because T cell numbers were only slightly reduced at 72 hr. However, one must pay attention to the fact that apoptosis data represent only a momentary picture in the course of the immune response whereas proliferation analyses (CFSE dilution assay) are reflecting the whole process. In fact, the dilution of CFSE upon cell divisions reveals the proliferative history of the adoptively transferred T cells and reflects the total amount of proliferation that has occurred in the course of the immune response.

Altogether, the initial response to an antigen is characterized by antigen-specific IL-2 production of Th cells and an early activation of antigen-specific and endogenous Treg cells (Figures 1 and 3). At low antigen concentrations (up to 8 µg OVA) there is: (i) activation of Th and Treg cells, but (ii) almost no proliferation of all cells because neither Treg nor Th cells reach the activation threshold for entry into it, and (iii) a slightly lower apoptosis rate of Treg cells compared to Th cells. At higher concentrations (>8 µg OVA) the cells expand proportionately to the antigen load, but it is obvious that living antigen-specific Th expand more than Treg cells (Figure 5B–D), which is in contrast to a previous report using IL-2 reporter mice (Amado et al., 2013).

Discussion

In order to contribute to the main question, how the adaptive immune system ensures tolerance at very low antigen concentration and adequate immune response at rising antigen amounts in respect of IL-2 expression and IL-2 action, we asked here: Is a graded antigen load translated into a binary or graded expression of the cytokine IL-2, the IL-2 receptor CD25, and the activation marker CD69 shortly after immunization of mice? Do binary processes indeed contribute to an enhanced robustness (Kellogg et al., 2015) of antigen-specific activation of Th cells? Do binary processes support antigen amount-dependent Treg cell-mediated tolerance? Is IL-2 acting in a paracrine as well as autocrine fashion on CD25 expression?

Binding of IL-2 influences CD25 expression in Th cells as shown by several groups using single cell in vitro experiments including mathematical modeling (Busse et al., 2010; Feinerman et al., 2008; Tkach et al., 2014). IL-2 enhances TCR-induced up-regulation of CD25 level at suboptimal stimulation (Feinerman et al., 2008). Our data clearly demonstrate a paracrine and ruled out a sole autocrine in vivo action of IL-2. Co-staining of CD25 and IL-2 revealed that IL-2-producing and IL-2-non-producing transgenic CD4+ T cells have the same expression of CD25 at all the different antigen concentrations used (Figure 2A) and over the whole time frame of IL-2 and CD25 co-expression (Figure 2C). These data are consistent with previous reports showing via STAT5 activation that IL-2 is delivered to neighboring T cells via synapses (Sabatos et al., 2008) and that IL-2 production does not coincide with STAT5 phosphorylation (Long and Adler, 2006). The observed peak times of IL-2 (6–14 hr) and CD25 (12–24 hr) show that IL-2 expression precedes CD25 expression in vivo, which is in good agreement with previous observations (Lischke et al., 2012; Long and Adler, 2006).

Beside the delivery of an important survival and maintenance signal for Treg cells, only a paracrine action of IL-2 allows in addition at least two other most recently discovered processes to occur. First, weaker T cell clones can be activated and proliferate just by integrating antigen and IL-2 response. Thus, IL-2 from strongly activated clones co-opt a fraction of weaker clones into activation and proliferation (Voisinne et al., 2015) and contribute thereby to continuous heterogeneity of the responding population. Second, IL-2 may act as a gatekeeper for T cell homeostasis during an immune response, initiating the down-modulation of excessive or boosting of weak effector T cell responses by the interplay of Th and Treg cells (Almeida et al., 2012, 2006; Amado et al., 2013; Reynolds et al., 2014).

In vitro studies claimed that many readouts downstream of TCR-activation such as IL-2, CD69, CD25, and ppERK are always binary (Feinerman et al., 2010; Tkach et al., 2014). However, there are conflicting data from in vitro experiments concerning CD69 and CD25 expression (Allison et al., 2016; Busse et al., 2010; Feinerman et al., 2010; Kimmig and Baumgrass, 2009; Podtschaske et al., 2007; Tkach et al., 2014). Using the adoptive transfer approach, we could confirm binary IL-2 expression in vivo, because the quantities of this cytokine on the cellular level were independent of the antigen concentrations used (Figure 1A). However, we also found that the overall IL-2 responses were not binary but graded at the population level. Increasing antigen amount resulted in a proportional higher percentage of IL-2-producing T cells. Thus, increasing or decreasing the number of IL-2 producers and not regulating the magnitude of IL-2 output of individual cells calibrates the antigen response. Interestingly, unlike IL-2, and in contrast to some in vitro data, CD25 and CD69 expression is clearly graded in vivo at both, single cell and population level (Figure 1B,C) (Feinerman et al., 2010; Tkach et al., 2014).

Binary responses of T cells play not just a role in filtering out noisy signals but moreover they can integrate information, such as timing or intensity, of the danger signal to adapt cell responses (Kellogg et al., 2015). On the other hand, as discussed recently by Tkach et al., “binary single parameter readouts carry some limitations: they are prone to saturation of input detection range, exhibit narrow output dynamic ranges, are subject to shifts in value by cellular death and migration, and provide little information on the functional capacity of activated cells” (Tkach et al., 2014).

Therefore, we asked whether the interplay of binary IL-2 and graded CD25 expression might have physiological advantages. To this end, we generated a computational model to compare the potential consequences of the observed binary and a hypothetical graded IL-2 secretion for the adaptive immune response. The model allowed us to simulate the outcome of these different scenarios taking into account quantitative spatial information in the spleen (Busse et al., 2010; Chapman et al., 1981; Lischke et al., 2012; Mueller and Germain, 2009; Sabatos et al., 2008; Smith, 2006). Our model confirmed that Treg cells become activated by IL-2 at lower antigen concentrations than Th cells and predicted that graded IL-2 secretion would result in activation of all T cells in a very narrow range of antigen stimulus, which is represented by a very steep slope of the activation curve of Th and Treg cells. Moreover, the binary IL-2 secretion results in a broadened linear range of initial T cell activation. This suits the needs of a fine-tuned immune response for a broad range of possible antigen loads.

Crucially, p-STAT5 measurements could substantiate our modeling results. Binary IL-2 secretion of Th cells led to a graded antigen dose-dependent increase in STAT5 phosphorylation in antigen-specific Th and Treg cells in vivo, which correlated perfectly with the activation curve of antigen-specific Th and Treg cells in our computational model. Very remarkably, not only the number of antigen-specific Treg cells but to a similar degree also endogenous Treg cells are p-STAT5 positive and correlate with the antigen amount. However, it is known that not only IL-2 signaling but also TCR signaling is required for the Treg cells to gain and sustain their full suppressive activity (Levine et al., 2014; Vahl et al., 2014).

Recently, Liu et al. provided information on how the spatial organization of IL-2-producing Th cells and activated Treg cells might help to maintain immune homeostasis (Liu et al., 2015). Using imaging flow cytometry (histo-cytometry) they discovered that STAT5 was phosphorylated in most clustering Treg cells and that an IL-2-producing Th cell was usually at the center of Treg cell clusters. This study convincingly highlighted the role of spatial proximity of few IL-2 producing Th cells and clustering of Treg cells close to them for critical IL-2-driven immune regulatory processes. The reduction of motility of Th and Treg cells within the clusters could secure effective IL-2 delivery for Treg cells even if IL-2 is released in a non-directional manner. Moreover, the number of Treg cells within the cluster might determine how many Treg cells are affected. The point, how IL-2 is shared between Th and Treg cells is extensively studied and discussed in two recent theoretical articles (Labowsky, 2016; Thurley et al., 2015).

Almeida et al. deduced a central role for IL-2 in homeostasis of Th cells by indexing of Treg cells to the number of IL-2 producing cells resulting in a constant ratio of both in IL-2 reporter mice (Almeida et al., 2006; Amado et al., 2013). In contrast, we showed in our transfer model a higher increase in the number of antigen-specific Th cells compared to Treg cells over the antigen concentration range used (Figure 5B–D). Nevertheless, the interdependence of IL-2 producing cells and Treg cells allows immune responses but also feedback control.

Most recently, Chinen et al. proved the essential role of IL-2R-dependent activation of STAT5 in controlling suppressor function of mature Treg cells (Chinen et al., 2016). Using genetic gain and loss of function approaches they demonstrated that augmented STAT5 activation in differentiated Treg cells increases the formation of Treg-DC-cell-conjugates and potentiates suppressor function in a TCR-independent manner.

In summary, we have shown that there are substantial differences between in vitro and in vivo experimental data concerning IL-2 and CD25 expression and regulation. Furthermore, we have demonstrated that computational modeling serves as a valuable and indispensable tool to understand complex regulatory mechanisms in vivo. Our insights contribute to a better understanding of antigen dose-dependent activation of the IL-2/CD25 axis of Th and Treg cells in an in vivo context and therefore under the real spatio-temporal conditions.

Materials and methods

Mice

OVA-TCR transgenic OT-II mice (Jackson Laboratory, Bar Harbor, ME, Stock 004194) were crossed to B6PL mice (Jackson Stock 000406). OVA-TCR transgenic DO11.10 mice (Jackson Stock 003303) were for some experiments additionally crossed with DEREG mice (Lahl et al., 2007) to obtain GFP reporter mice for Foxp3. All these mice including non-transgenic C57BL/6NCrl (Charles River, Sulzfeld, Germany, Strain Code 027) were bred under specific-pathogen-free conditions in the animal facility of the Federal Institute for Risk Assessment (Berlin, Germany). BALB/cAnNCrl mice were purchased from Charles River (Strain Code 028).

Adoptive transfer

Request a detailed protocol

T cell receptor transgenic cells were isolated from spleens according to standard procedures. Naive OT-II cells were enriched by positive MACS sort with L-Selectin (CD62L) beads (Miltenyi Biotec, Bergisch Gladbach, Germany) and transferred into C57BL/6 mice (2.5–3.0 × 106 transgenic T cells per recipient). DO11.10 cells were depleted of CD8+ cells by negative MACS sort with CD8 beads (Miltenyi Biotec) and transferred into BALB/c mice (1.0–1.5 × 106 transgenic T cells per recipient). To analyze proliferation, cells were labeled with CFSE (Invitrogen) before adoptive transfer, according to standard protocols. The mice were immunized intravenously with different amounts of endotoxin-free Ovalbumin (<5 pg of endotoxin/mg of protein by Limulus amebocyte lysate assay) and 10 µg lipopolysaccharide (LPS, from E. coli O55:B5, Sigma-Aldrich, St. Louis, MO, USA) as adjuvant. For IL-2 in vivo treatment, mice received i.p. either 20 µg recombinant IL-2 or 2.5 µg IL-2 complexed with 10 µg anti-IL-2 (clone JES6-5H4 targets IL-2 to CD122). Splenocytes from C57BL/6 or BALB/c recipient mice were isolated and analyzed at indicated times after immunization.

IL-2 secretion assay

Request a detailed protocol

IL-2 secretion assays (Miltenyi Biotech) were performed with 7.5 × 106 splenocytes from C57BL/6 mice as previously described (Assenmacher et al., 1998). In brief, cells were washed, resuspended in ice-cold buffer (0.5% BSA, 2 mM EDTA in PBS) and incubated with the anti-IL-2 antibody capture matrix on ice for 5–10 min. After adding pre-warmed medium (RPMI 1640 + 4% FCS), cells were incubated at 37°C in a shaking water bath for 45 min. The secretion phase was stopped by adding ice cold PBS/BSA/EDTA buffer to the cells. IL-2 secreting cells were stained with anti-IL-2-PE antibody and different surface marker antibodies for 20 min on ice. To control whether the IL-2 capture matrix is saturated by endogenous IL-2 secretion, 0.4 µg/ml recombinant murine IL-2 (eBioscience) was added 5 min before the end of the secretion phase (see Figure 1—figure supplement 1).

Flow cytometric analysis

Request a detailed protocol

Cell suspensions from spleen were counted with a Guava EasyCyte capillary flow cytometer and ViaCount Assay (Merck Millipore, Germany). For reduction of unspecific antibody binding, cells were preincubated with 100 µg/ml 2.4G2 (FcγRII/III; ATCC, Manassas, VA) and 50 μg/ml purified rat Ig (Nordic, Tilburg, The Netherlands) for 5 min. Surface staining was performed on ice for 20 min with monoclonal antibodies conjugated to FITC, PE, PerCP-Cy5.5, PE-Cy7, APC, AF700, and PacB: αCD4 (RM4-5), αCD25 (PC61.5), αB220 (RA3-6B2), αOVA-TCR (KJ1-26), αThy1.1 (OX-7; all purified from hybridoma supernatants), and αCD69 (H1.2F3; eBioscience, San Diego, CA, RRID:AB_465119). Dead cells were excluded from analysis by DAPI (4-,6-diamidino-2-phenylindole) staining.

For detection of early apoptotic cells, either PE-conjugated annexin V (BioLegend, RRID:AB_2561298) or an Alexa Fluor 660-labeled inhibitor of caspases (FLICA; Immunochemistry Technologies) was used in combination with DAPI.

When analyzing CFSE dilution profiles the percentage of dividing progenitors, i.e. the precursor frequency of dividing cells, was calculated as described (Lyons and Parish, 1994).

For intracellular staining of the transcription factor p-STAT5 (Alexa Fluor 647-labeled antibody clone 47/p-STAT5, BD Biosciences) spleens were disintegrated directly into BD Fixation Buffer containing 2% paraformaldehyde and stained in BD Perm Buffer III (BD Biosciences) for 30 min on ice.

For intracellular staining of the transcription factors NFATc2 (Cy5-labeled own polyclonal antibodies (Bendfeldt et al., 2012), c-Fos (Alexa Fluor 488-labeled rabbit polyclonal IgG antibodies from Santa Cruz Biotechnology, RRID:AB_2231996), and Foxp3 (PE-Cy7-labeled antibody clone FJK-16s from eBioscience, RRID:AB_891554) the cells were stained with 1.34 µM Pacific Orange succinimidyl ester to exclude dead cells, fixed with Foxp3 fixation buffer (eBioscience), stained in Foxp3 permeabilization buffer (eBioscience), and analyzed using a LSR II flow cytometer (BD Biosciences). Data were analyzed with FlowJo software (Treestar, Ashland, OR).

Computational model

Basic features of the model

Request a detailed protocol

The model describes a population of T cells on a 2D circular domain (spleen/lymph node) and explicitly includes the secretion and diffusion of IL-2 and the formation of IL-2:IL-2R-complexes in space and time via a set of reaction-diffusion equations. Cells are modeled as circular areas with reactive boundaries that account for either secretion of IL-2 or IL-2:IL-2R-complex formation. Dirichlet boundary conditions were imposed on the outer boundary of the domain. The cells were distributed randomly with a mean cell-to-cell distance derived from considerations about the number of T cells and the size of the T cell zone in the spleen. The dynamics of IL-2:IL-2R-interactions on the cell boundaries was included in terms of a simple Michaelis-Menten type model which accounts for complex formation, dissociation and internalization. We assumed that the receptor dynamics can be approximated by considering only the influence of the expression of the high affinity chain CD25. The influence of unspecific CD25 positive cells was modeled in terms of a homogeneous background degradation of IL-2. The model was written and evaluated in MATLAB making use of the Partial Differential Equation Toolbox.

Implementation of binary/graded scenarios

Request a detailed protocol

Binary (all-or-none) IL-2 secretion was implemented by randomly choosing a subset of producers that all secrete IL-2 at the maximal rate. In the graded scenario, by contrast, the maximum amount of producers secrete IL-2 at a rate that changes with increasing antigen stimulus. The fraction of producers in the binary case and the secretion rate in the graded case were chosen such that the total amount of IL-2 secreted is the same for both scenarios. The quantitative relationship between IL-2 secretion rate/CD25 expression and the antigen stimulus were determined by evaluating the corresponding FACS experiments (Figure 1).

Output of the model

Request a detailed protocol

Simulations were performed for varying antigen stimuli according to the doses used in the corresponding FACS experiments. For each stimulus the process was simulated over a time span of 48 hr which encompasses the initial increase and eventual decline of both IL-2 secretion and CD25 expression. In each simulation the distribution of cells and the random selection of producers/consumers was carried out anew. For each consumer cell the total amount of internalized complex was measured over the course of the simulation and taken as a proxy for downstream STAT5 signaling (not included in the model). We assumed that a cell is activated when the amount of complex exceeds a critical threshold. This threshold was chosen such as to optimize the fit between the fraction of activated cells in the model and the fraction of p-STAT5+ cells in the experiment (Figure 3C).

Equations

The basic scheme of the diffusion equation used in the model is

∂u∂t=D∂2u∂x2−dUR⋅IL2R⋅uu+k,

where u represents the solution for the concentration of IL-2 and is the diffusion coefficient. The second term on the right hand side accounts for the homogeneous background degradation due to unspecific regulatory cells. IL2R is the expression rate of CD25 and the factor dUR stands for the density of unspecific Tregs. The parameter k can be interpreted as the concentration at which the consumption of IL-2 is half its maximum value. This parameter depends on the rates of formation, dissociation, and internalization of complexes in the following way:

k=kdeg+koffkon.

The secretion of IL-2 by producer cells was accounted for in terms of generalized Neumann conditions on the outer cell boundaries:

n⋅(D∇u)=IL2⋅f(A,t)cCell,

where n is the outward unit normal, IL2 the maximal secretion rate, and cCell the cell circumference. The factor f(A,t) accounts for the dependence of the secretion rate on the antigen stimulus and the time after stimulation. Dependence on time was modeled as the composition of a linearly increasing function from 0 to the maximum value between t = 0 and t = 9 hr and a linearly decreasing function from the maximum to 0 between t = 9 hr and t = 18 hr (compare to Figure 3).

The interaction of IL-2 and IL-2R on the surface of consumer cells is also accounted for in terms of generalized Neumann boundary conditions:

n⋅(D∇u)=−IL2R⋅g(A,t)cCell⋅uu+k.

Here, IL2R denotes the maximal expression of CD25, while g(A,t) accounts for the dependence of CD25 expression on stimulus and time. Time dependence in this case was modeled as a linear increase between t = 0 and t = 18 hr and a subsequent decrease between t = 18 hr and t = 30. The dependence of CD25 expression on the antigen stimulus is different for T helper cells and regulatory T cells and was determined by evaluating the corresponding FACS experiments.

References

Article and author information

Author details

  1. Franziska Fuhrmann

    1. Robert Koch Institute, Berlin, Germany
    2. German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    FF, Conception and design, Acquisition of data, Analysis and interpretation of data
    Contributed equally with
    Timo Lischke and Fridolin Gross
    Competing interests
    The authors declare that no competing interests exist.

  2. Timo Lischke

    German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    TL, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Franziska Fuhrmann and Fridolin Gross
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0413-4252

  3. Fridolin Gross

    Institute for Theoretical Biology, Charité University Medicine, Berlin, Germany
    Contribution
    FG, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Franziska Fuhrmann and Timo Lischke
    Competing interests
    The authors declare that no competing interests exist.

  4. Tobias Scheel

    German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    TS, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.

  5. Laura Bauer

    1. Robert Koch Institute, Berlin, Germany
    2. German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    LB, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.

  6. Khalid Wasim Kalim

    German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    KWK, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.

  7. Andreas Radbruch

    1. German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    2. Charité University Medicine, Berlin, Germany
    Contribution
    AR, Conception and design, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.

  8. Hanspeter Herzel

    Institute for Theoretical Biology, Charité University Medicine, Berlin, Germany
    Contribution
    HH, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0414-7889

  9. Andreas Hutloff

    1. Robert Koch Institute, Berlin, Germany
    2. German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    AH, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.

  10. Ria Baumgrass

    German Rheumatism Research Center Berlin (DRFZ), A Leibniz Institute, Berlin, Germany
    Contribution
    RB, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    baumgrass@drfz.de
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3289-1608

Funding

Bundesministerium für Bildung und Forschung (T-Sys (e:Bio))

  • Andreas Radbruch
  • Hanspeter Herzel
  • Ria Baumgrass

Deutsche Forschungsgemeinschaft (SFB650)

  • Ria Baumgrass

Deutsche Forschungsgemeinschaft (TR52)

  • Ria Baumgrass

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Tim Sparwasser, Jochen Hühn, Ute Hoffmann and Alf Hamann for providing DO11.10×DEREG mice. We are grateful to the staff of the animal facility and the labmanagers. We thank Jenny Kirsch and Toralf Kaiser, operators of the flow cytometry core facility (FCCF). We thank Melanie Venzke, Christian Gabriel, Manja Jargosch, Anna Abajyan, Martin Karl, Katharina Hecklau, Stefanie Gryzik and Yen Hoang for lively discussion and helpful advices.

Ethics

Animal experimentation: All animal experiments were performed according to state guidelines and approved by the responsible governmental authority LAGeSo (Landesamt für Gesundheit und Soziales) Berlin in animal experiment licenses T 0187/01 and G 0070/13.

Copyright

© 2016, Fuhrmann et al.

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

Metrics

  • 2,060
    views
  • 307
    downloads
  • 15
    citations

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

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Franziska Fuhrmann
  2. Timo Lischke
  3. Fridolin Gross
  4. Tobias Scheel
  5. Laura Bauer
  6. Khalid Wasim Kalim
  7. Andreas Radbruch
  8. Hanspeter Herzel
  9. Andreas Hutloff
  10. Ria Baumgrass
(2016)
Adequate immune response ensured by binary IL-2 and graded CD25 expression in a murine transfer model
eLife 5:e20616.
https://doi.org/10.7554/eLife.20616

Share this article

https://doi.org/10.7554/eLife.20616

Categories and tags

Research organism

Further reading

    1. Cancer Biology
    2. Computational and Systems Biology

    Luminal epithelial cells integrate variable responses to aging into stereotypical changes that underlie breast cancer susceptibility

    Rosalyn W Sayaman, Masaru Miyano ... Mark LaBarge
    Research Article

    Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger <30 (y)ears old and older >55y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.

    1. Computational and Systems Biology
    2. Microbiology and Infectious Disease

    An antimicrobial drug recommender system using MALDI-TOF MS and dual-branch neural networks

    Gaetan De Waele, Gerben Menschaert, Willem Waegeman
    Research Article

    Timely and effective use of antimicrobial drugs can improve patient outcomes, as well as help safeguard against resistance development. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is currently routinely used in clinical diagnostics for rapid species identification. Mining additional data from said spectra in the form of antimicrobial resistance (AMR) profiles is, therefore, highly promising. Such AMR profiles could serve as a drop-in solution for drastically improving treatment efficiency, effectiveness, and costs. This study endeavors to develop the first machine learning models capable of predicting AMR profiles for the whole repertoire of species and drugs encountered in clinical microbiology. The resulting models can be interpreted as drug recommender systems for infectious diseases. We find that our dual-branch method delivers considerably higher performance compared to previous approaches. In addition, experiments show that the models can be efficiently fine-tuned to data from other clinical laboratories. MALDI-TOF-based AMR recommender systems can, hence, greatly extend the value of MALDI-TOF MS for clinical diagnostics. All code supporting this study is distributed on PyPI and is packaged at https://github.com/gdewael/maldi-nn.

    1. Computational and Systems Biology
    2. Genetics and Genomics

    Functional characteristics and computational model of abundant hyperactive loci in the human genome

    Sanjarbek Hudaiberdiev, Ivan Ovcharenko
    Research Article

    Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.