# Mathematical model of naive T cell division and survival IL-7 thresholds

^{1}Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds, UK^{2}Centre for Immunology and Infection, University of York, York, UK

We develop a mathematical model of the peripheral naive T cell population to study the change in human naive T cell numbers from birth to adulthood, incorporating thymic output and the availability of interleukin-7 (IL-7). The model is formulated as three ordinary differential equations: two describe T cell numbers, in a resting state and progressing through the cell cycle. The third is introduced to describe changes in IL-7 availability. Thymic output is a decreasing function of time, representative of the thymic atrophy observed in aging humans. Each T cell is assumed to possess two interleukin-7 receptor (IL-7R) signaling thresholds: a survival threshold and a second, higher, proliferation threshold. If the IL-7R signaling strength is below its survival threshold, a cell may undergo apoptosis. When the signaling strength is above the survival threshold, but below the proliferation threshold, the cell survives but does not divide. Signaling strength above the proliferation threshold enables entry into cell cycle. Assuming that individual cell thresholds are log-normally distributed, we derive population-average rates for apoptosis and entry into cell cycle. We have analyzed the adiabatic change in homeostasis as thymic output decreases. With a parameter set representative of a healthy individual, the model predicts a unique equilibrium number of T cells. In a parameter range representative of persistent viral or bacterial infection, where naive T cell cycle progression is impaired, a decrease in thymic output may result in the collapse of the naive T cell repertoire.

## 1. Introduction

The number of naive T cells in the periphery is determined by a balance between cell loss (death or differentiation) and cell renewal due to cell division and thymic export (1, 2). In humans, at least, the decline in thymic export occurs mainly in childhood, from about a year of age until 20 years of age, when the number of naive T cells is increasing (3). In adults, the decline in thymic export is much less pronounced but the number of naive T cells is, more or less, constant (4). Survival of the naive T cell population in the periphery depends on both common gamma chain cytokines and weak “tonic” signals induced by recognition of self-peptides by the T cell receptor (TCR) (5, 6). IL-7 is required for the homeostatic expansion of naive CD8^{+} and CD4^{+} T cells in lymphopenic hosts, while naive T cells disappear over a 1-month period upon adoptive transfer into IL-7 deficient (IL-7^{−}) hosts (7–9).

Signals from recognition of self-peptides bound to major histocompatibility complex (sp-MHC), and IL-7, promote cell survival. Naive T cell survival is impaired when removing access to one of these signals (10–14). Of interest are the mechanisms by which these signals are regulated, and that result in a stable number of naive T cells throughout the lifetimes of mice and humans. In this paper, we focus on IL-7 as a master regulator of naive T cell survival (15). IL-7 is produced by stromal cells in tissues, including fibroblastic reticular cells, marginal reticular cells, and lymphatic endothelial cells (16). These cells produce very small amounts of IL-7 messenger RNA, consistent with IL-7 protein levels limiting T cell expansion. IL-7 is a heparin-sulfate binding protein, and as such, it will bind extra-cellular matrix surrounding stromal cells. Thus, the interaction between naive T cells and stroma controls their homeostasis (17). Recognition of higher affinity, non-self-peptides by the T cell receptor induces naive T cells to undergo an alternative, IL-7 independent, survival program dependent on IL-2 (18).

Naive CD8^{+} T cell responses depend on the amount of IL-7 cells are exposed to (19). At low IL-7 concentrations (<10^{−2} ng ml^{−1}), cell viability was impaired; at higher concentrations (>1 ng ml^{−1}) cells were observed to proliferate in response to IL-7. This difference might arise from changes in the strength of the IL-7R induced signal the cell receives. For an individual cell, IL-7R induced signaling must be greater than some threshold to prevent the accumulation of pro-apoptotic proteins. Similarly, IL-7R signaling must be greater than a second, higher, threshold to induce cell division. Heterogeneity at the single cell level in IL-7 signaling thresholds (a property reported to depend on expression of IL-7R), resulted in differential survival and division (19). Although these observations are based on two different CD8^{+} T cell receptor transgenic mice, it is assumed that the key principles regarding T cell survival will be found in the repertoire of naive CD4^{+} and CD8^{+} T cells.

We introduce a deterministic mathematical model of the naive T cell population to study the change in human naive T cell numbers from birth to adulthood. We will assume cell survival depends on the availability of IL-7. We do not include availability of sp-MHC as a variable within the model, but assume sp-MHC availability is sufficient to allow cell survival and proliferation, in conjunction with sufficient IL-7 stimulus. We also make the approximation that heterogeneity is constant with changes in age. For a mathematical study of the impact sp-MHC availability has on clonal diversity, the reader is referred to Stirk et al. (20, 21). Our model is a mathematical description of the homeostasis of the naive T cell repertoire, but does not consider stimulation by foreign antigens.

## 2. Materials and Methods

### 2.1. A Mathematical Description of the Size of the Peripheral Naive T Cell Population

Stochastic processes provide a method of treating each cell as a distinct, countable object, and permit a more realistic model than a deterministic characterization. Fluctuations in the number of cells can be considered but, in a non-linear stochastic model, approximations are often made to facilitate the analysis. In the linear noise approximation (22), for example, fluctuations are assumed to be of order ${\mathrm{\Omega}}^{\frac{1}{2}}$ for a system of size Ω. The human peripheral T cell compartment is estimated to contain of the order of 10^{11} T cells (3). Letting the system size be the average number of naive T cells in humans, we find 𝒪(Ω) = 10^{11} cells, and correspondingly, fluctuations are expected to be typically 10^{5} − 10^{6} cells in magnitude. That is, we expect fluctuations of approximately 0.001% in the size of the human naive T cell pool due to stochasticity in the per cell division and death rates. Based on these considerations, adopting a deterministic approach to describe the total human peripheral naive T cell population is reasonable.

We assume peripheral naive T cells are either in a resting state, or proceeding through the cell cycle. The deterministic variables *R*(*t*) and *C*(*t*) are introduced to model the total number of T cells in the resting and cycling states, respectively. The variable *I*(*t*) is introduced to model the concentration of IL-7. The deterministic approach we take does not consider any notion of space. Indeed, this approach is tantamount to assuming the resource, IL-7, is shared equally amongst all cells. Competition for the resource is introduced only so far as each cell acts to reduce the global concentration of the resource. Resting cells may receive a signal which induces them to proceed through one round of division. Upon completion of the cell cycle, a cycling cell produces two daughter cells in the resting compartment. Resting cells are assumed to die if the IL-7 induced survival signal is insufficient; cells may also die during cell cycle. The input of cells from the thymus into the resting compartment, in keeping with observations in humans, is a decreasing function of time (23, 24). Production of IL-7 is related to the size of the lymphatic system architecture, which we estimate from the body mass of an individual. In the absence of T cells, IL-7 is assumed to be degraded and/or consumed by other cell types at a constant rate. Upon signal induction through the IL-7 receptor, IL-7 is assumed to be consumed by the T cell. A diagrammatic representation of the model is given in Figure 1.

**Figure 1. Diagrammatic illustration of the deterministic model**. T cells leaving the thymus enter the resting naive peripheral pool. Cells in either a resting or cycling state may die. The rate of death from the resting state depends on the availability of the resource (IL-7), whereas the death rate for cycling cells is constant. Resting cells enter the cell cycle at a rate that depends on the availability of the resource (IL-7). Cycling cells produce two daughter cells in the resting state upon completion of the cell cycle.

### 2.2. IL-7 Signaling and Heterogeneity in IL-7 Thresholds

In the model, IL-7 signaling is assumed to be uniform across the population. Yet, we introduce heterogeneity in the signaling thresholds for survival and proliferation. Let *S*(*t*) be the average signaling strength across the naive T cell population. Each cell experiences the same strength of signaling for a given concentration of IL-7, *I*(*t*). We relate the internal signaling to the concentration of IL-7 by the equation

The functional form and constants of this relationship are derived from the study of IL-7 receptor dynamics summarized in the Appendix. We assume each individual cell possesses a unique pair of thresholds for survival and proliferation. Furthermore, we assume, in the continuous limit, that these thresholds are distributed log-normally across the entire population of cells. The use of log-normal distributions guarantees, first of all, that all signaling thresholds are positive real numbers. Secondly, the log-normal distribution ensures that no cell can survive or divide independently of IL-7.

Let the random variable Θ* _{s}* represent the survival threshold, and let Θ

*represent the proliferation threshold. We write*

_{p}The respective probability density functions are

Estimates for θ* _{s}* and θ

*are obtained from the model summarized in the Appendix. Lauffenburger et al. found a significant change in cell viability for both OT-1 and F-5 T cells at around 10*

_{p}^{−2.5}ng ml

^{−1}IL-7 (19). Proliferation of OT-1 cells occurred above 1 ng ml

^{−1}IL-7. We estimate the equilibrium signaling at these concentrations as 60 and 600 units, respectively, setting θ

*= 60 and θ*

_{s}*= 600. Modeling heterogeneity in IL-7 responses by assuming heterogeneity in the IL-7 signaling thresholds, allows us to avoid modeling the naive population using either: (i) a PDE approach, where heterogeneity is continuous across the population of cells, or (ii) describing each subset of cells sharing common thresholds with its unique governing set of ODEs. However, these approaches present an obvious avenue for further research beyond the scope of this paper.*

_{p}#### 2.2.1. Death rate of resting cells

Suppose each T cell in the naive population is a distinct member possessing a unique signaling threshold for survival. As in the previous section, we assume these signaling thresholds are distributed log-normally (in the continuous limit). We suppose the death rate of an individual cell is Boolean, in the sense that if the global signaling strength is greater than the cell’s individual survival threshold, then the cell can survive indefinitely. Similarly, if the signaling strength is below the survival threshold, the cell will undergo apoptosis at a rate μ* _{R}*. The death rate for cell

*i*with survival threshold ${\mathrm{\theta}}_{s}^{\left(i\right)}$, is given by

In the continuous limit (assuming signaling thresholds are distributed log-normally), the average death rate for the population of naive T cells is given by

where ${p}_{{\mathrm{\Theta}}_{s}}\left(\mathrm{\theta}\right)$ is the probability density function of the random variable Θ* _{s}*, defined by equation (3) with

*x*=

*s*.

#### 2.2.2. Rate of entry into cell cycle

Analogous to Section 2.2.1, we assume each T cell in the naive population is a distinct member possessing a unique signaling threshold for proliferation. We let the individual rate of entry into cell cycle be given by

Assume, in the continuous limit, the signaling threshold for entry into the cell cycle is represented by the random variable Θ* _{p}*, defined in equation (2), with probability density function ${p}_{{\mathrm{\Theta}}_{p}}\left(\mathrm{\theta}\right)$ (equation (3),

*x*=

*p*). The average rate of entry into cell cycle is given by

#### 2.2.3. Cell cycle progression

Cycling cells take on average λ^{−1} days to complete the cell cycle. After a cell divides, both daughter cells are produced in the resting state and require a second signal before they can progress through another round of cell division. Cell cycle may be interrupted resulting in the death of the cell. Such death events occur at a rate μ* _{C}*.

#### 2.2.4. Thymic export

We assume thymic output to be a decreasing function of time. In particular, we use the functional form given by Bains et al. (3). Let us introduce the thymic output function, *v*(*t*), as follows

where *t* corresponds to the age of the individual, measured in days. A plot of this function is shown in the left panel of Figure 2. The function was chosen by Bains et al. to describe the rate of thymic export of CD4^{+} T cells. We use the same function to describe the export rate of all naive T cells (CD4^{+} or CD8^{+} T cells). This approximation is justified since we require the absolute cell count to roughly approximate the cell count observed in humans (indeed, such an observation is likely subject to large differences). Of interest later in the paper is the relative variation of cell numbers with different choices of parameter values. For our purposes, the important feature of *v*(*t*) is that it is a decreasing function of time.

**Figure 2. Left panel: rate of thymic export**. Plot of equation (8). Middle panel: IL-7 internalization rate. Plot of the function γ(*I*(*t*). Right panel: plot of the function *m*(*s*) (equation (11)), a model for the average body mass of male individuals.

#### 2.2.5. Internalization of IL-7

We use the model summarized in the Appendix to estimate the rate of IL-7 internalization. The total number of IL-7 molecules internalized by a single T cell in 1 day, exposed to IL-7 at concentration *I* ng ml^{−1}, is described by the function

It is reported IL-7 has a molecular mass of 17 kDa (≈2.8 × 10^{−11} ng) (25). Based on this, we define the per cell rate of IL-7 internalization to be

In order to convert the rate of change of mass to the rate of change of concentration, we must choose a volume for the system. Naive T cells are typically found in the lymph nodes, spleen, and gut of the human body. We shall make the rough estimation that the total volume is of the order of 1 l. This implies the rate of IL-7 internalization by all naive T cells in the population is given by

where *R*(*t*) is the number of resting naive T cells. The IL-7 internalization rate, γ(*I*(*t*), is shown in the middle panel of Figure 2.

#### 2.2.6. Production of IL-7

We assume the rate of IL-7 production is proportional (with proportionality constant $\widehat{\mathrm{\beta}}$) to the body mass of an individual. To estimate average body mass we use the model given by Burmaster and Crouch (26). The explicit relationship between mass and age is given by the function *m*(*s*) as follows

where *s* is age measured in years. A plot of this function is given in the right panel of Figure 2. The rate of IL-7 production is given by the function

where *t* denotes age as measured in days.

#### 2.2.7. Intra-cellular degradation of IL-7

We assume IL-7 is degraded and internalized by other cell types at a constant rate. We let the degradation rate of IL-7 be denoted by the parameter δ.

#### 2.2.8. Deterministic mathematical model

From the above assumptions, the system of differential equations governing the behavior of the naive T cells (resting and cycling) and the concentration of IL-7 is given by

This system is subject to the initial conditions *I*_{0}, *R*_{0}, and *C*_{0}.

## 3. Results

### 3.1. Parameter Estimates

The function describing the rate of internalization of IL-7, equation (10), and the signaling relation, equation (1), were both estimated from our studies of IL-7 receptor dynamics in naive T cells (summarized in the Appendix). In the same study we have estimated the average signaling thresholds to be, respectively, θ* _{s}* = 60 and θ

*= 600. Our estimate for θ*

_{p}*was found to be close (θ*

_{p}*≈ 585, which we have rounded to 600) to the limit of the signaling function equation (1), as*

_{p}*I*(

*t*) → +∞. These estimates imply that, on average, approximately half of the naive T cell population does not proliferate in response to excess IL-7 when receptor dynamics is in equilibrium. Such a finding is consistent with the observations made by Lauffenburger et al. in the experiments described in Ref. (19). In these experiments, F-5 T cells did not proliferate (whereas OT-1 cells did proliferate), even in excess amounts of IL-7. Estimates for the IL-7 average signaling thresholds, θ

*and θ*

_{s}*, have been based on OT-1 cells. However, we note that changes to these thresholds only result in quantitative differences, provided θ*

_{p}*< θ*

_{s}*.*

_{p}We assume cycling naive T cells take 12 h to complete the cell cycle and produce two daughter cells (λ^{−1} = 0.5 day). Resting naive T cells are assumed to die after 2 days of IL-7 starvation $\left({\mathrm{\mu}}_{R}^{-1}=2\phantom{\rule{0.5em}{0ex}}\mathrm{day}\right)$. We, furthermore, set the rate of entry into cell cycle to be ρ = 5 day^{−1}. Cycling cells are assumed to die at rate μ* _{C}* = 1 day

^{−1}in healthy individuals, whereas we choose μ

*= 3 day*

_{C}^{−1}(>λ) to be representative of cell cycle impairment. These parameters have been estimated from the literature. We note that the model behavior we discuss in the following sections was found to be robust to changes in these parameters. Robustness was concluded since a 10-fold change in each or any combination thereof, of these four parameters did not change the qualitative behavior of the model, provided the relation λ > μ

*(or λ < μ*

_{C}*) was maintained. We examine the changes in model behavior arising from altering the relative values of λ and μ*

_{C}*in the following sections.*

_{C}We choose δ such that δ*I*(*t*) is similar in magnitude to γ(*I*(*t*)*R*(*t*), when *R*(*t*) ≈ 10^{11}, and *I*(*t*) ≈ 10^{−2} ng ml^{−1}. The proportionality constant $\widehat{\mathrm{\beta}}$ is chosen such that we observe 𝒪(10^{11}) naive T cells in equilibrium at 20 years of age (for *I*(*t*) ≈ 10^{−2} ng ml^{−1}). We had the least inclination when choosing the parameter α. This parameter describes the spread in the individual signaling thresholds across the naive T cell population. We choose α = 2, however this choice has no justification from the literature. The parameter set is summarized in Table 1.

### 3.2. There Exists a Unique and Steady State When λ > μ_{C}

_{C}

Let us suppose changes in thymic output and IL-7 production occur in time scales slower than those of the changes in the number of naive T cells. Under this assumption, we look for adiabatic solutions of the system as follows:

For the parameter set studied, the relative error between this solution and the exact solution is within 3% for the resting naive T cell population and the concentration of IL-7, and within 14% for the cycling T cell population (see Figure 3).

**Figure 3. Relative error between the adiabatic solution and the exact solution, where initial conditions were chosen to be equal to the adiabatic solution at time t = 0**.

All numerical results presented in the paper have been obtained with a Python code^{1}: differential equations (13), (14), and (15) have been solved using a fourth-order Runge-Kutta scheme. Quasi-stationary solutions were found using the scipy.optimize package. Bifurcation plots were computed using a bisection scheme to search for multiple solutions of equation (21) in the interval [0,1]. Corresponding T cell numbers were calculated using equations (19) and (20).

Notice that the adiabatic solution for both cell types is uniquely defined for a given value of cytokine concentration, $\widehat{I}\left(t\right)$, namely

The problem of finding adiabatic solutions can then be reduced to finding solutions, $\widehat{I}\left(t\right)$, to the one-dimensional equation

By construction, $\mathrm{\gamma}\left(\widehat{I}\left(t\right)\right)$ is a monotonically increasing function of $\widehat{I}\left(t\right)$. Furthermore, the existence of positive adiabatic solutions requires $\mathrm{\beta}\left(t\right)>\mathrm{\delta}\widehat{I}\left(t\right)$. Therefore, for a fixed time *t*, the left-hand side of equation (21) is an increasing function of $\widehat{I}\left(t\right)$. For λ > μ* _{C}*, the right-hand side of equation (21) is a monotonically decreasing function of $\widehat{I}\left(t\right)$. Lastly, $\mathrm{\gamma}\left(\widehat{I}\left(t\right)\right)=0$ for $\widehat{I}\left(t\right)=0$, and the limit as $\widehat{I}\left(t\right)\to 0$ of the right-hand side is equal to μ

*. It follows that the intersection of the left and right sides must be unique and positive. We deduce that for λ > μ*

_{R}*there exists a unique adiabatic solution to the system (see left plot of Figure 4). This solution is stable for the parameter set given in Table 1. We have also numerically explored parameter space, but have not found a parameter set for which this solution is unstable. In Figure 5 we present numerical solutions to equations (13–15), computed with a fourth-order Runge-Kutta method implemented in Python. Initial conditions were chosen to be the adiabatic solutions of equations (16–18) at*

_{C}*t*= 0.

**Figure 4. Left panel: for λ > μ_{C} (t = 25 years), the right-hand side of equation (21) is a decreasing function of I(t) (red line)**. There exists a unique, asymptotically stable solution $\widehat{I}\left(t\right)$ found at the intersection of the solid black and red lines. Right panel: there exist three intersections between the red and solid black lines when λ < μ

*for*

_{C}*t*= 25 years. We require $\mathrm{\beta}\left(t\right)>\mathrm{\gamma}\left(\widehat{I}\left(t\right)\right)$ for existence of stable solutions, therefore we neglect all intersections with the dashed black line.

**Figure 5. Numerical solutions to equations (13–15) for λ > μ_{C}**. Initial conditions are chosen to be the adiabatic solutions of equations (16–18) at

*t*= 0. We observe a large increase in the naive T cell population from birth to adulthood. During this period we see a 10

^{7}-fold increase in the number of cycling cells, resulting from increased IL-7 availability and reduced thymic output. IL-7 availability increases continually as the individual ages.

### 3.3. There Exist Two Steady States When λ < μ_{C}

_{C}

Let us now consider the case λ < μ* _{C}*. The right-hand side of equation (21) is no longer a decreasing function of $\widehat{I}\left(t\right)$. We find there may exist up to three solutions $\widehat{I}\left(t\right)$, two of which may be stable simultaneously, whilst the third is unstable (see right plot of Figure 4). Further examination (by numerically finding all solutions using the bisection method) of the model reveals a saddle-node bifurcation as thymic output changes with age (see left panel, Figure 6). We assume that individuals with (around) 10

^{11}naive T cells are healthy, in the sense that they have a sufficient number of T cells to provide protection against immune challenges. When a cell in cycle is more likely to die than to produce two daughter cells, we define cell cycle progression to be impaired. Mathematically, this corresponds to λ < μ

*. When cell cycle progression is impaired, the model predicts an individual may possess a healthy number of T cells, thereby being immuno-competent, up until the age at which the model bifurcates. For the parameter set we have investigated, this bifurcation is inevitable given the estimated decline in thymic output established by Bains et al. (3). Indeed, for a given parameter set, from the known rate at which thymic output declines, one can estimate the age at which the model bifurcates. The bifurcation results in a decrease in the naive T cell population size of approximately two orders of magnitude, that is, following the bifurcation we expect roughly 99 out of every 100 naive T cells to be lost.*

_{C}**Figure 6. Left panel: bifurcation diagram for varying the parameter v(t) for λ < μ_{C}**. Right panel: bifurcation plot as a function of varying the rate of IL-7 production. Thymic export is fixed at

*v*(60 years).

## 4. Discussion

In a healthy individual, it is reasonable to expect that naive T cells entering the cell cycle are more likely to complete division and produce two daughter cells, than to die during the division process. Therefore, we suppose the parameter relation λ < μ* _{C}* represents a healthy individual. We have shown, under this hypothesis, that the model allows a single asymptotically stable adiabatic solution. The total number of naive T cells per kilogram of body mass was found to increase over the first 18 years of life and decrease thereafter. The decline in this ratio for adults is seemingly a consequence of the reduction in thymic output.

The number of cycling naive T cells in adiabatic conditions was found to increase from practically 0 (approximately 0.75 cells) cells at birth to 10^{7} cells by adulthood. Whilst the number of cycling cells increased thereafter, the increase was at a markedly slower rate. The increase in the number of cycling cells is probably due to the decline in thymic output, wherein competition for IL-7 decreases. Despite the increase in the cycling population, the number of cells in cycle is several orders of magnitude smaller than the resting population. For adults, this proportion is approximately [0.01, 0.04]% of the total naive population.

When cell cycle is impaired (λ < μ* _{C}*), the model exhibits a saddle-node bifurcation as thymic output declines. For the parameter set investigated, the two adiabatic solutions are generally separated by two orders of magnitude. This feature motivates the following theoretical scenario: suppose a healthy individual experiences some event which causes naive cell cycle progression to become impaired at 18 years of age. The model predicts the naive T cell repertoire will experience a dramatic reduction in cell numbers at roughly 33 years of age. See Figure 7 for the full solution of the model illustrating this scenario. From an immunological perspective, this scenario is interesting because the noticeable effect of the event (the dramatic loss of naive T cells) occurs roughly 15 years later than the event itself (naive cell cycle progression to become impaired at 18 years of age). Indeed, it is the reduction in thymic output that triggers the loss rather than the event itself. If the thymus was not atrophic, the loss in naive T cells would not occur.

**Figure 7. Theoretical scenario in which we set μ_{C} = 3 at 18 years of age**. At roughly 33 years of age the model bifurcates resulting in a dramatic loss of naive T cells.

Suppose now we fix thymic output at a constant rate and let the cell cycle be impaired. The model undergoes a saddle-node bifurcation as the rate of IL-7 production changes. More specifically, the model bifurcates, resulting in a decrease in the total number of T cells, as a consequence of increasing IL-7 production. Intuitively, this can be understood as follows: for increased IL-7 production, the rate of entry into cell cycle is enhanced, however, since cells are more likely to die in cell cycle than to produce two daughter cells, the enhanced rate of entry into cell cycle actually serves to decrease the total amount of naive T cells. The bifurcation diagram for this scenario is shown in the right panel of Figure 6. We found the critical value of this bifurcation decreases with thymic output. In the limiting case, corresponding to thymic output at 60 years of age, we found this critical value to be approximately 9.1 ng ml^{−1}. Consider again the theoretical scenario in which a healthy individual undergoes an event resulting in cell cycle impairment at 18 years of age. Suppose now at age 25 the same individual undergoes some treatment to limit IL-7 production to 8 ng ml^{−1}, corresponding to thymic production at approximately 11 years of age. The T cell count is reduced by approximately 60%. This reduction is a significant improvement on the 99% T cell loss in the untreated individual at 33 years of age. For this theoretical scenario, the model predicts limiting IL-7 availability will partially avoid the dramatic T cell loss arising from reducing thymic export when the cell cycle is impaired. See Figure 8 for the full model solution in the treated individual.

**Figure 8. At age 18 the individual experiences some event resulting in impaired cell cycle progression**. At age 25 the individual undergoes treatment to limit IL-7 production to levels comparable to those at 11 years of age. Whilst the T cell count is reduced, this attrition is a significant improvement on the reduction observed when the individual has not received treatment (see Figure 7).

In the model presented here we have neglected the fact naive T cells become activated in response to recognition of ligand specific to their unique TCR. Research by Koenen et al. has shown T cell survival is IL-7 independent following T cell activation (18). Suppose now we include a term in the governing ODEs to represent differentiation of naive T cells into cells with a different phenotype, such as activated T cells. Assuming no reversion back to the naive phenotype, such a term would appear in the model as a loss term equivalent to the death rate ${\widehat{\mathrm{\mu}}}_{R}\left(S\left(t\right)\right)$. The simplest approach to including differentiation (due to activation) would be to assume naive T cells differentiate into activated T cells at a constant rate proportional to the rate of antigenic challenge. In this case, we would replace the term ${\widehat{\mathrm{\mu}}}_{R}\left(S\left(t\right)\right)$ by ${\widehat{\mathrm{\mu}}}_{R}\left(S\left(t\right)\right)+{\mathrm{\mu}}_{D}$, where μ* _{D}* is constant. Consider again the red curves in Figure 4. The differentiation term will cause a translation in the red curve of length μ

*up the vertical axis. For λ > μ*

_{D}*, there still exists a unique solution, however there will be a quantitative change in comparison to the case when μ*

_{C}*= 0. When λ < μ*

_{D}*, there still exists the possibility of a saddle-node bifurcation for a general parameter set. For the parameter set we have investigated, there exists a maximum value ${\mathrm{\mu}}_{D}^{*}$, for a given time*

_{C}*t*, such that we can find more than one solution. For differentiation rates ${\mathrm{\mu}}_{D}>{\mathrm{\mu}}_{D}^{*}$, we can only find a unique solution corresponding to the stable adiabatic one, in which we have reduced T cell numbers.

In this paper we have developed and analyzed a deterministic mathematical model of a population of naive T cells, whose survival depends on the availability of the cytokine IL-7. We have shown this model predicts a stable population of cells when cell cycle progression is healthy. More interestingly, when cell cycle progression is impaired, our results indicate declining thymic output may result in a dramatic loss in the number of naive T cells. Furthermore, we have been able to establish that limiting IL-7 production partially rescues this decline. However, our study has been restricted to naive T cells and we have not taken into account the accumulation of memory T cells as an individual ages. Memory T cells are generated in response to antigen or homeostatic cytokines (27), and during repeated homeostasis-driven divisions of naive T cells (28). Previous studies have shown the percentage of memory T cells increases with age, yet the percentage of naive T cells decreases with age (29). CD4^{+} memory T cells also require IL-7 for survival and proliferation in the periphery (30). In the case of CD8^{+} memory T cells, IL-15 is largely responsible for their peripheral survival, but IL-7 may also be required. It is, then, reasonable to assume that memory T cells compete for the IL-7 required for naive T cell survival. At least, that is, for those memory T cells which access IL-7 in the same tissues as naive T cells, such as the lymph nodes. As mentioned before, memory T cells can be derived from repeated divisions of naive T cells (31). In our model, the naive T cells that will acquire a memory-like phenotype first are those with low IL-7 survival and proliferation thresholds. We would expect, then, that over time the ratio of naive to memory T cells would decrease. Furthermore, the naive population would lose those cells with lowest survival and division thresholds. Given the distribution of signaling thresholds in our model, we would expect to see a shift to the right for the distribution of both survival and signaling thresholds. This shift, which implies that on average naive T cells require a higher concentration of IL-7 to enter cell cycle with age, together with the additional competition from an increasing memory population, is expected to cause a decrease in the total number of naive T cells and the percentage of naive T cells in cell cycle. These considerations might help explain the discrepancy between our model and data showing a reduction in the percentage of naive T cells expressing Ki67^{+} during childhood^{2}. There is support for the fact that childhood might be the phase in which most memory T cells are acquired (3). Indeed, it might be that the model presented here, better describes the total number of T cells that require IL-7, naive and memory, rather than just the naive T cell pool alone. This, however, would require a more detailed analysis, out of the scope of this paper.

Many different receptor-mediated signals are integrated by T cells in their micro-environment, such as cytokines, adhesion molecules, T cell receptors, and co-receptors (CD28, CTLA-4). However, the survival of naive T cells, the focus of this paper, has been shown to depend on IL-7 and low level TCR stimulation (8, 9, 15, 32). Other cytokines, such as IL-15, play a significant role in the homeostasis of memory CD8^{+} T cells (33). Given the experimental support for the hypothesis that IL-7 is critical for the homeostatic proliferation and the survival of naive T cells, in this first model we have neglected other signals (7).

Previous mathematical models of naive T cell homeostasis have focused on the relative contribution of thymic export and cell division in the periphery (3, 34, 35). These models conclude that in humans, thymic export makes an important contribution to the size of the naive T cell population in early life (<20 years of age), whereas later in life the number of naive T cells is maintained by homeostatic proliferation in the periphery. Some recent studies have measured the relative contribution of thymic export by examining the average number of TRECs in naive T cells (3, 35, 36). The use of mathematical models has allowed these groups to infer the relative (young versus old) kinetics for naive T cells, based on experimental estimates of the total number of naive T cells and recent thymic emigrants (3, 34, 35). In the model introduced here, we have aimed to provide a mechanistic perspective by investigating at the molecular and cellular levels, the role IL-7 plays in regulating the homeostasis of naive T cells (37). The increase in the proportion of cycling cells in our model is in agreement with previous experimental studies (38). Our study suggests the increase in peripheral division rates can be explained by the availability of IL-7, which is a consequence of a combined effect: (i) an increased net IL-7 production as an individual ages, and (ii) a reduction of recently exported thymocytes competing for this resource.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

We thank Robin Callard, Megan Palmer, and Douglas Lauffenburger for sharing their advice on the approaches taken in this study. This work has been partially supported by the research grant BB/G023395/1 (BBSRC, Carmen Molina-París).

## Footnotes

## References

1. Almeida A, Rocha B, Freitas A, Tanchot C. Homeostasis of T cell numbers: from thymus production to peripheral compartmentalization and the indexation of regulatory T cells. *Semin Immunol* (2005) **17**(3):239–49. doi: 10.1016/j.smim.2005.02.002

2. Surh CD, Sprent J. Homeostasis of naive and memory T cells. *Immunity* (2008) **29**(6):848–62. doi:10.1016/j.immuni.2008.11.002

3. Bains I, Thiébaut R, Yates A, Callard R. Quantifying thymic export: combining models of naive T cell proliferation and TCR excision circle dynamics gives an explicit measure of thymic output. *J Immunol* (2009) **183**(7):4329–36. doi:10.4049/jimmunol.0900743

4. den Braber I, Mugwagwa T, Vrisekoop N, Westera L, Mögling R, Bregje de Boer A, et al. Maintenance of peripheral naive T cells is sustained by thymus output in mice but not humans. *Immunity* (2012) **36**(2):288–97. doi:10.1016/j.immuni.2012.02.006

5. Takada K, Jameson S. Naive T cell homeostasis: from awareness of space to a sense of place. *Nat Rev Immunol* (2009) **9**(12):823–32. doi:10.1038/nri2657

6. Sprent J, Surh C. Normal T cell homeostasis: the conversion of naive cells into memory-phenotype cells. *Nat Immunol* (2011) **12**(6):478–84. doi:10.1038/ni.2018

7. Tan J, Dudl E, LeRoy E, Murray R, Sprent J, Weinberg K, et al. IL-7 is critical for homeostatic proliferation and survival of naive T cells. *Proc Natl Acad Sci U S A* (2001) **98**(15):8732–7. doi:10.1073/pnas.161126098

8. Schluns K, Kieper W, Jameson S, Lefrançois L. Interleukin-7 mediates the homeostasis of naive and memory CD8 T cells in vivo. *Nat Immunol* (2000) **1**(5):426–32. doi:10.1038/80868

9. Fry T, Mackall C. The many faces of IL-7: from lymphopoiesis to peripheral T cell maintenance. *J Immunol* (2005) **174**(11):6571–6.

10. Martin B, Bécourt C, Bienvenu B, Lucas B. Self-recognition is crucial for maintaining the peripheral CD4^{+} T-cell pool in a nonlymphopenic environment. *Blood* (2006) **108**(1):270–7. doi:10.1182/blood-2006-01-0017

11. Brocker T. Survival of mature CD4 T lymphocytes is dependent on major histocompatibility complex class II-expressing dendritic cells. *J Exp Med* (1997) **186**(8):1223–32. doi:10.1084/jem.186.8.1223

12. Markiewicz MA, Brown I, Gajewski TF. Death of peripheral CD8^{+} T cells in the absence of MHC class I is Fas-dependent and not blocked by Bcl-xL. *Eur J Immunol* (2003) **33**(10):2917–26. doi:10.1002/eji.200324273

13. Vivien L, Benoist C, Mathis D. T lymphocytes need IL-7 but not IL-4 or IL-6 to survive in vivo. *Int Immunol* (2001) **13**(6):763–8. doi:10.1093/intimm/13.6.763

14. Kondrack R, Harbertson J, Tan J, McBreen M, Surh C, Bradley L. Interleukin 7 regulates the survival and generation of memory CD4 cells. *J Exp Med* (2003) **198**(12):1797–806. doi:10.1084/jem.20030735

15. Fry T, Mackall C. Interleukin-7: master regulator of peripheral T-cell homeostasis? *Trends Immunol* (2001) **22**(10):564–71. doi:10.1016/S1471-4906(01)02028-2

16. Mueller SN, Germain RN. Stromal cell contributions to the homeostasis and functionality of the immune system. *Nat Rev Immunol* (2009) **9**(9):618–29. doi:10.1038/nri2588

17. Kang J, Coles M. IL-7: the global builder of the innate lymphoid network and beyond, one niche at a time. *Semin Immunol* (2012) **24**(3):190–7. doi:10.1016/j.smim.2012.02.003

18. Koenen P, Heinzel S, Carrington EM, Happo L, Alexander WS, Zhang J-G, et al. Mutually exclusive regulation of T cell survival by IL-7R and antigen receptor-induced signals. *Nat Commun* (2013) **4**:1735. doi:10.1038/ncomms2719

19. Palmer M, Mahajan V, Chen J, Irvine D, Lauffenburger D. Signaling thresholds govern heterogeneity in IL-7-receptor-mediated responses of naïve CD8^{+}; T cells. *Immunol Cell Biol* (2011) **89**(5):581–94. doi:10.1038/icb.2011.5

20. Stirk ER, Lythe G, van den Berg HA, Molina-París C. Stochastic competitive exclusion in the maintenance of the naïve T cell repertoire. *J Theor Biol* (2010) **265**(3):396–410. doi:10.1016/j.jtbi.2010.05.004

21. Stirk ER, Lythe G, van den Berg HA, Hurst GA, Molina-París C. The limiting conditional probability distribution in a stochastic model of T cell repertoire maintenance. *Math Biosci* (2010) **224**(2):74–86. doi:10.1016/j.mbs.2009.12.004

22. Wallace EW. A simplified derivation of the linear noise approximation. *arXiv Preprint arXiv:1004.4280* (2010).

23. Mitchell WA, Meng I, Nicholson SA, Aspinall R. Thymic output, ageing and zinc. *Biogerontology* (2006) **7**(5–6):461–70. doi:10.1007/s10522-006-9061-7

24. Steinmann GG. Changes in the human thymus during aging. In: Müller-Hermelink HK, editor. *The Human Thymus*. Berlin: Springer (1986). p. 43–88.

25. Haugen F, Norheim F, Lian H, Wensaas AJ, Dueland S, Berg O, et al. IL-7 is expressed and secreted by human skeletal muscle cells. *Am J Physiol Cell Physiol* (2010) **298**(4):C807–16. doi:10.1152/ajpcell.00094.2009

26. Burmaster DE, Crouch EA. Lognormal distributions for body weight as a function of age for males and females in the united states, 1976–1980. *Risk Anal* (1997) **17**(4):499–505. doi:10.1111/j.1539-6924.1997.tb00890.x

27. Geginat J, Lanzavecchia A, Sallusto F. Proliferation and differentiation potential of human CD8+ memory T-cell subsets in response to antigen or homeostatic cytokines. *Blood* (2003) **101**(11):4260–6. doi:10.1182/blood-2002-11-3577

28. Goldrath AW, Bogatzki LY, Bevan MJ. Naive T cells transiently acquire a memory-like phenotype during homeostasis-driven proliferation. *J Exp Med* (2000) **192**(4):557–64. doi:10.1084/jem.192.4.557

29. Saule P, Trauet J, Dutriez V, Lekeux V, Dessaint J-P, Labalette M. Accumulation of memory T cells from childhood to old age: central and effector memory cells in CD4(+) versus effector memory and terminally differentiated memory cells in CD8(+) compartment. *Mech Ageing Dev* (2006) **127**(3):274–81. doi:10.1016/j.mad.2005.11.001

30. Seddon B, Tomlinson P, Zamoyska R. Interleukin 7 and T cell receptor signals regulate homeostasis of CD4 memory cells. *Nat Immunol* (2003) **4**(7):680–6. doi:10.1038/ni946

31. Cho B, Rao V, Ge Q, Eisen H, Chen J. Homeostasis-stimulated proliferation drives naive T cells to differentiate directly into memory T cells. *J Exp Med* (2000) **192**(4):549. doi:10.1084/jem.192.4.549

32. Fry T, Mackall C. Interleukin-7: from bench to clinic. *Blood* (2002) **99**(11):3892–904. doi:10.1182/blood.V99.11.3892

33. Tan JT, Ernst B, Kieper WC, LeRoy E, Sprent J, Surh CD. Interleukin (IL)-15 and IL-7 jointly regulate homeostatic proliferation of memory phenotype CD8+ cells but are not required for memory phenotype CD4+ cells. *J Exp Med* (2002) **195**(12):1523–32. doi:10.1084/jem.20020066

34. Thomas-Vaslin V, Altes HK, de Boer RJ, Klatzmann D. Comprehensive assessment and mathematical modeling of T cell population dynamics and homeostasis. *J Immunol* (2008) **180**(4):2240–50.

35. Murray JM, Kaufmann GR, Hodgkin PD, Lewin SR, Kelleher AD, Davenport MP, et al. Naive T cells are maintained by thymic output in early ages but by proliferation without phenotypic change after age twenty. *Immunol Cell Biol* (2003) **81**(6):487–95. doi:10.1046/j.1440-1711.2003.01191.x

36. Ribeiro RM, Perelson AS. Determining thymic output quantitatively: using models to interpret experimental T-cell receptor excision circle (TREC) data. *Immunol Rev* (2007) **216**(1):21–34. doi:10.1111/j.1600-065X.2006.00493.x

37. Fallon EM, Lauffenburger DA. Computational model for effects of ligand/receptor binding properties on interleukin-2 trafficking dynamics and T cell proliferation response. *Biotechnol Prog* (2000) **16**(5):905–16. doi:10.1021/bp000097t

38. Prelog M, Keller M, Geiger R, Brandstätter A, Würzner R, Schweigmann U, et al. Thymectomy in early childhood: significant alterations of the CD4(+)CD45RA(+)CD62L(+) T cell compartment in later life. *Clin Immunol* (2009) **130**(2):123–32. doi:10.1016/j.clim.2008.08.023

39. Park J, Yu Q, Erman B, Appelbaum J, Montoya-Durango D, Grimes H, et al. Suppression of IL7Rα transcription by IL-7 and other prosurvival cytokines: a novel mechanism for maximizing IL-7-dependent T cell survival. *Immunity* (2004) **21**(2):289–302. doi:10.1016/j.immuni.2004.07.016

40. Henriques CM, Rino J, Nibbs RJ, Graham GJ, Barata JT. IL-7 induces rapid clathrin-mediated internalization and JAK3-dependent degradation of IL-7Rα in T cells. *Blood* (2010) **115**(16):3269–77. doi:10.1182/blood-2009-10-246876

41. Palmer MJ, Mahajan VS, Trajman LC, Irvine DJ, Lauffenburger DA, Chen J, et al. Interleukin-7 receptor signaling network: an integrated systems perspective. *Cell Mol Immunol* (2008) **5**(2):79–89. doi:10.1038/cmi.2008.10

## Appendix

### A.1. A Model of IL-7R Dynamics

In this section we present a summary of a stochastic model in which we consider the number of IL-7 receptors on the surface of naive T cells. This study is used to derive equations (1) and (9) in the main text. The model has been formulated as a continuous time Markov process. We consider the number of receptors on a single naive T cell. For our purposes, we present the mean field approximation to the model with which we estimate the parameters. We introduce the following four variables:

• *m*_{1}(*t*) – the total number of IL-7 receptors on the surface of a naive T cell,

• *m*_{2}(*t*) – the number of unbound internalized receptors,

• *m*_{3}(*t*) – the number of bound internalized receptors,

• *m*_{4}(*t*) – the quantity of IL-7 induced signal.

We assume, in equilibrium, the fraction of surface receptors bound to IL-7 is given by *f _{I}*. The system of ODEs denoting the mean field approximation to the stochastic model is given by:

### A.2. Reduced Model in the Case When *I* = 0

Consider a T cell in an IL-7 free medium with initial conditions such that the IL-7 induced signaling vanishes and the number of IL-7:IL-7R internal complexes is zero. Then the mean field model can be reduced to the following set of ODEs:

This reduced system is governed by four parameters ϕ, ξ* _{U}*, σ

*, and ζ*

_{U}*, and possesses the following stable steady state:*

_{U}We assume in the reduced model 10% of the total number of receptors are internalized in equilibrium. We set ${m}_{1}^{*}=9{m}_{2}^{*}$. Based on the measurements of Singer et al. (39), we shall assume 4 × 10^{4} receptors in total when the reduced model is in steady state. Therefore, we let

In Ref. (40), cells were cultured with the translation inhibitor cycloheximide (CHX) to prevent transcription of the IL-7 receptor. Total expression of the IL-7 receptor was measured over several time points, from which the authors estimate the half-life of the receptor in an unstimulated cell to be approximately 24 h.

In the reduced model, all receptors are guaranteed to be degraded in a finite amount of time. The expected time for a receptor, that is initially on the cell surface, to be degraded in the lysosome is given by

Assuming exponential decay, the half-life for a receptor to undergo lysosomal degradation, starting on the cell surface, is then given by

Thus, we can write

Combining equations (A5), (A6), and (A7) we find

The value of ξ* _{U}* relative to ζ

*effectively dictates the ratio of receptors which are degraded to those recycled back to the cell surface. We assume the system has evolved to minimize waste of functional proteins and tentatively let ξ*

_{U}*> ζ*

_{U}*. That is, we assume that a greater fraction of receptors are recycled back to the surface of the cell. We somewhat arbitrarily set*

_{U}**Figure A1. Diagrammatic representation of the transition probabilities of the stochastic model for the IL-7 and IL-7 receptor system**.

### A.3. Receptor-Ligand Kinetics

Suppose the number of surface receptors is constant and denoted by *R _{T}*. Let us also assume the extra-cellular concentration of IL-7 is constant and denoted by

*I*. Define

*R*(

_{B}*t*) to be the number of IL-7 receptors bound to IL-7. Note that we assume the time to recruit the common gamma chain, ${\mathrm{\gamma}}_{c}$, is negligible. Then, we can describe changes in the number of bound complexes by the following ODE

where *k*_{+} and *k*_{−}are, respectively, the binding and unbinding rates of the IL-7:IL-7R receptor-ligand system. We assume the timescales for this reaction are faster than the timescales for changes in receptor numbers, such that we can consider these reactions to be in equilibrium. This ODE has a unique stable steady state:

The Appendix for Ref. (41) provides estimates for *k*_{+} and *k*_{−}, from which we set *k*_{+} = 1 nM min^{−1} and *k*_{−} = 0.1 min^{−1}. It is reported IL-7 has a molecular mass of around 17 kDa, from which we estimate the ratio *k*_{−}/*k*_{+} ≈ 1.7 ng ml^{−1}.

### A.4. Early Internalization Events

In Ref. (40), surface receptor expression was assessed in human thymocytes. It is reported, cells in 50 ng ml^{−1} IL-7 culture down-regulated IL-7R expression. A 20% reduction was observed after 10 min. Using the above estimate for *k*_{−}/*k*_{+}, we find *f _{I}*|

_{I}_{=50}≈ 0.97. Based on this, we can neglect internalization of the unbound receptor. In the first 10 min, we shall also neglect recycling and inhibition of receptor transcription. We assume surface receptor expression loss is modeled by the ODE

Given initial surface expression levels equal to *m*_{1}(0), this ODE has solution

We assume the previous estimates obtained from the reduced model for *m*_{1}(0) and ϕ. That is, we let *m*_{1}(0) = 3.6 × 10^{4} receptors and ϕ = 1.2 × 10^{3} receptors h^{−1}. Then using the above expression for *m*_{1}(*t*) equation (A8), we obtain an estimate for σ* _{B}*. We find σ

*≈ 1.4 h*

_{B}^{−1}. The authors of Ref. (40) estimate the half-life of the IL-7 receptor in cells treated with CHX, cultured in 50 ng ml

^{−1}, to be approximately 3 h. The expected time to lysosomal degradation for a surface receptor is given by

Assuming exponential decay, with a half-life of 3 h, we find ξ* _{B}* ≈ 0.2ζ

*+ 0.3 h*

_{B}^{−1}. Again, without a direct measurement, let us set ξ

*= 1 h*

_{B}^{−1}, to obtain an estimate for ζ

*≈ 3.5 h*

_{B}^{−1}.

### A.5. Remaining Parameters

Consider the observation of an approximately 98% reduction in surface receptor expression following overnight culture in 6 ng ml^{−1} IL-7 made in Ref. (39). We let ${\left.{m}_{1}^{*}\right|}_{I=6}=0.02\phantom{\rule{0.5em}{0ex}}{\left.{m}_{1}^{*}\right|}_{I=0}=0.02\frac{\mathrm{\varphi}\left({\mathrm{\xi}}_{U}+{\mathrm{\zeta}}_{U}\right)}{{\mathrm{\sigma}}_{U}{\mathrm{\zeta}}_{U}}\approx 763\phantom{\rule{0.5em}{0ex}}\mathrm{receptors}$. We set the derivatives equal to zero in the system of ODEs equations (A1–A4), and manipulate the resulting system of equations to obtain

From the steady solutions of ODEs equations (A3) and (A4), we find

We use this expression to rearrange equation (A9) in terms of ${m}_{4}^{*}$. This gives

Solving the above expression numerically, we find ${m}_{4}^{*}\u2215\mathrm{\kappa}\approx 0.1$. We estimate a value for χ based on the observation that following culture in 6 ng ml^{−1} IL-7, mRNA levels took approximately 12 h to return to 99% of the control levels. The transcription rate is given by ϕexp(−*m*_{4}(*t*)/κ). We again assume the IL-7 induced signal decays according to the equation *m*_{4}(*t*) = *m*_{4}(0)e^{−}* ^{χt}*, where

*m*

_{4}(0) = 0.1. Combing these assumptions with the experimental observations, we have

from which we find χ ≈ 0.19 h^{−1} ⇒ φ ≈ 6.1 × 10^{−4} receptors^{−1} h^{−1} κ. The parameter κ was chosen to be 1000. This choice was made from the stochastic model: κ = 1000 is the minimum value (to the nearest power of 10) such that fluctuations in the signaling quantity are greater than zero for low (10^{−2} ng ml^{−1}) concentrations of IL-7. Using this value for κ, we find ϕ ≈ 0.61 h^{−1}. The parameter estimates are summarized in Table A1.

#### A.5.1. Changes in the concentration of IL-7

The functional form

is used to approximate the numerical solution of ${m}_{4}^{*}$ as a function of the concentration of IL-7. Using this function we find *a* ≈ 600 signaling units and *b* ≈ 0.025 ng ml^{−1}. In a similar manner we use the functional form

to approximate the numerical solution of ${m}_{1}^{*}$ as a function of the concentration of IL-7. We find *c* ≈ 600 receptors, *d* ≈ 1000 receptors ng ml^{−1} and *e* ≈ 0.03 ng ml^{−1}. The number of IL-7 molecules internalized by each T cell per day is assumed to be the same as the number of internalized IL-7:IL-7R complexes per day. We let the complex internalization rate, as a function of *I*, be given by

Keywords: IL-7, T cell, homeostasis, threshold, IL-7R, mathematical model, thymic output

Citation: Reynolds J, Coles M, Lythe G and Molina-París C (2013) Mathematical model of naive T cell division and survival IL-7 thresholds. *Front. Immunol.* **4**:434. doi: 10.3389/fimmu.2013.00434

Received: 31 July 2013; Accepted: 22 November 2013;

Published online: 23 December 2013.

Edited by:

Miles Davenport, University of New South Wales, AustraliaReviewed by:

Edward John Collins, The University of North Carolina at Chapel Hill, USAChristian Schönbach, Nazarbayev University, Kazakhstan

Robin Callard, University College London, UK

Copyright: © 2013 Reynolds, Coles, Lythe and Molina-París. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Carmen Molina-París, Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK e-mail: carmen@maths.leeds.ac.uk