|
|
||||||||
Theoretical Biology, Utrecht University, Utrecht, The Netherlands
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Cell division of memory cells is also the dominant maintenance mechanism for the T cell repertoire in human adults (8). Division rates of human naive CD45RA+ and memory CD45RO+ T cells have been estimated by the rates at which patients treated with radiotherapy lose lymphocytes with chromosome damage (9, 10). It was thus estimated that memory CD45RO+ T cells divide once every 22 wk (10). The CD45RA+ naive T cells were estimated to divide every 3.5 yr (10). Having an order of magnitude of 1011 naive T cells in either the CD4+ or the CD8+ compartment, dividing once every 1000 days amounts to a production of order of magnitude of 108 naive T cells/day in each subclass of T cells. Similarly, an average division rate of memory T cells of once every 100 days amounts to a production of about 109 memory T cells/day. Note that these numbers are probably upper estimates because these patients are recovering from radiotherapy and have not yet attained normal steady state peripheral T cell counts.
Cancer patients treated with chemotherapy typically fail to reconstitute a normal T cell repertoire, and the reconstitution rates correlate with the enlargement of the thymus and the appearance of CD45RA+ CD4+ naive T cells (11). For example, the three oldest patients in this study (11) are 24 yr old, and have daily recovery rates of 0.54, 0.26, and 0.63 CD4+ T cells/µl blood. For adults having approximately 5 liters of blood, and having about 2% of the CD4+ T cells in the circulation (12), these figures indicate an average total body production of about 108 CD4+ T cells per day during this recovery phase. This is similar to the naive T cell production estimated above by chromosome damage (9). Similarly, multiple sclerosis patients treated with anti-CD4 mAb have markedly depleted peripheral CD4+ T cell counts (13). Upon withdrawal of the treatment the CD4+ T cell counts again recover very slowly, with similar recovery rates of CD45RO+ memory and CD45RA+ naive T cells (13).
| Materials and Methods |
|---|
|
|
|---|
More recently, it became clear that activated T cells can express telomerase activity (22, 23, 24, 25, 26, 27), and avoid telomere shortening during the first 10 to 12 divisions in culture (23, 27). Thus, the CD4+ T cell telomere shortening of 33-bp year-1 (14) and the 1400-bp difference between the T cell subpopulations (14, 15) paradoxically coexist with apparent telomerase activity. This implies that there should be cell division that is not fully compensated by telomerase activity. Developing a mathematical model in which naive and memory T cells are maintained both by renewal and clonal expansion, we are able to resolve this paradox by showing that one can obtain full consistency with the data even if telomerase activity fully compensates during clonal expansion. If memory T cells do divide more frequently than naive T cells during renewal (10), we obtain similar telomere shortening rates of both subpopulations as a steady state solution of our model. This accounts for the fixed difference of average telomere length between T cell populations, irrespective of division rate of the memory T cells, and irrespective of telomerase activity during clonal expansion.
| Results |
|---|
|
|
|---|
|
Let ni be the number of T cells of telomere
loss index, i, in which i = 0 denotes cells
that have an initial maximal telomere length. Cells may die, with rate
, or divide, with rate
. Division of a cell of index i
clearly removes it from the cell count ni, while
its two daughter cells are added to cell count
ni+1. The basic dynamics of our model can thus
be expressed by the following system of differential equations
![]() | (1) |
0. The
cellular division rates
and/or death rates
are expected to be
regulated homeostatically. For deriving our results, it is not required
to specify these functions, however.
The dynamics of the total number of T cells, N, is found
simply by summing the ni equations
![]() | (2) |
One advantage of this simplification is that one can derive a very
simple differential equation describing the growth of the "average
telomere loss index," µ. This index, µ, should be proportional to
the average loss of telomere length (i.e., L =
L0 - bµ, where L is the measured average
telomere length, L0 is its starting value, and
b is the number of base pairs lost per cell division. The
average division index, µ, is simply the mean value of the telomere
loss index, i.e.,
![]() | (3) |
![]() | (4) |
. According to this model,
the telomere-shortening rate is independent of the average telomere
length. In tumor cells, additional factors have recently been
described, such as the telomeric repeat-binding factor TRF1
(31), which are involved in repairing telomere shortening
dependent on the telomere length. We have not incorporated such factors
in our model because it seems unlikely that similar factors play a
significant role in human T cells, which are well known to suffer from
telomere shortening, and, ultimately, senescence, by cell division
(14, 15, 22, 23, 24, 25, 26, 27). Telomere loss per cell division
Experimental measurements of the average rate of telomere loss are
based upon a longitudinal in vitro expansion of cell lines
(17, 18, 19, 20, 21). Implicitly assuming synchronous cell divisions,
and ignoring cell death, one typically computes the number of bases
lost per mean population doubling (MPD)3
as
![]() | (5) |
L is the total telomere loss over
the experiment, N0 is the number of cells the
culture started with, and Nt is the final
population size. These experiments have yielded 50
bMPD
100.
Because cell division is not expected to be synchronous, this
bMPD is not equal to the number of base pairs
lost per cell division (i.e., to our parameter b). We here
calculate the difference between the two. Solving Equation 2
, one
obtains that (
-
)t = ln
[Nt/N0]. Employing the
solution of Equation 4
, i.e., µ = 2
t, and our
definition that b =
L/µ, and also
ignoring cell death, i.e., assuming
<<
, we obtain the
expression
![]() | (6) |
![]() | (7) |
Thus, the true average rate of telomere loss per cell division is
72% of the previous estimates of telomere loss per MPD. The existing
data (17, 18, 19, 20, 21) then yield an estimate of 35
b
70. For definiteness we use b = 50
bp/cell division from now on. The rate at which a population of
dividing T cells loses its telomere base pairs (i.e., 2b
)
can now be equated to the experimental value of 33 base pairs/yr
(14). This yields a division rate of
= 1/3 per year.
Fortuitously, this is identical to the previous estimate
(14).
Finally, as total T cell numbers tend to remain at steady state,
we may put dN/dt = 0 in Equation 2
and
obtain
=
, as expected. Hence, the steady state average cellular
life span (1/
) equals 1/
= 3 yr.
Two-compartment model
In Figure 2
we generalize the simple
one-compartment model to a model describing two coupled compartments,
for "naive" and "memory" T cells, respectively. We assume
separate homeostatic control of division and/or death rates for the
naive and the memory cell types, since one would otherwise expect one
of the two types to outcompete the other (32, 33).
Additionally, mouse experiments do suggest that the naive and memory
subsets are regulated independently (34). The human
telomere data (14) demonstrate that naive CD4+
T cells (or at least their precursors) are, indeed, slowly dividing
cells. Experiments with mice also argue for slow division rates of
CD45RA+ naive T cells (2, 35). Although in a
rat model it has recently been suggested that, in the absence of
antigenic stimulation, memory CD4+ T cells revert to a
naive phenotype (36), it remains unclear whether one
should allow for a similar reversal in the human situation. In the
present model, we do not allow for such a reversion of the memory
marker (see alternative model under Discussion). Thus, we
devise a two-compartment model in which naive T cells cannot only
divide or die, but can also be primed by external Ags to expand into a
clone of memory T cells.
|
![]() | (8) |
0 as before. Extending our
previous notation, the
V and
V denote the division and death rates, which
will probably be homeostatically controlled by the total count of naive
cells V
ivi. The new
parameter,
, is the probabilistic rate at which naive T cells are
externally primed to become a proliferating clone of activated T cells.
Clearly, all results of the one-compartment model carry over to the
naive cells of this two-compartment model. We need only repeat the
analysis of the one-compartment model to obtain for the total number of
naive cells
![]() | (9) |
iivi/V
![]() | (10) |
As before, the naive telomere lengths decrease with twice the
naive T cell division rate. The only difference is that at a steady
state of Equation 9
we obtain
V =
V
+
, which is simply the "net turnover rate" of naive cells.
The "priming" rate,
, thus increases the naive cell
telomere-shortening rate, but this is estimated to be a negligible
contribution (see below).
Turning to the memory cell compartment, we have to consider clonal
expansion. The 1400-base pair difference in average telomere lengths
between human memory and naive T cells was interpreted to reflect a
clonal expansion of 14 to 28 cell divisions (14, 15). In
our model, we consider naive T cells that may be primed with external
Ags at a probabilistic rate
per day, to ultimately yield a clone of
C memory T cells. Due to the telomerase activity evoked by
the activation with Ag and/or the costimulatory factors
(22, 23, 24, 25, 26, 27), the degree of clonal expansion is unlikely to be
reflected in a proportional degree of telomere loss (23, 27). Thus, we write that right after clonal expansion a memory T
cell has shifted its telomere loss index by an amount K.
Note that perfect compensation by telomerase activity during clonal
expansion would correspond to K = 0.
In Appendix B, we derive that the total number of memory cells
M satisfies
![]() | (11) |
![]() | (12) |
M, but
depends on the clonal expansion, the naive:memory ratio, and their
difference in average telomere lengths. Thus, memory telomere lengths
are much less directly related to the division or turnover rates than
the very simple proportionality Equation 10
The difference
in naive and memory telomere lengths
For memory T cells, it is apparent from Equation 12
that the
shortening rate slows down when µM exceeds
µV + K, which is the mean index of cells that
have just been produced by clonal expansion of a primed (former naive)
cell. To analyze this effect, it is helpful to write a differential
equation for the distance between the naive and memory telomere loss
index, which we define as
µM -
µV. Subtracting Equation 10
from Equation 12
gives
![]() | (13) |
Inspecting this result, one notices a "time constant"
= M/(
CV), and a "source" term
=
2(
M -
V) + K/
. Since we only
consider data from adults (14), in whom memory and naive T
cell counts tend to remain in steady state,
and
can be treated
as finite constants allowing us to simplify Equation 13
into
![]() | (14) |
![]() | (15) |
. Comparing these results to experimental
data, we can draw several interesting conclusions.
Importantly, our results imply that the steady state
telomere-shortening rate of memory cells should approach that of the
naive cells, after a time period of a few times
. The empirical
results on CD4+ T cells (14) demonstrate
precisely this similarity of shortening rates. Hence, in steady state,
the rate of telomere shortening in memory T cells should approach
2
V, i.e., that of the naive T cells, and
should not have any relation to the parameters of the memory T cells.
There is, indeed, good agreement between both telomere-shortening rates
and the average division rates
V
0.0008
day-1 of human naive T cells (9, 10, 14).
This resolves what appeared to be a contradiction between measured
memory cell telomere-shortening rates and data suggesting that memory T
cells are dividing considerably more frequently than naive T cells
(1, 2, 9, 10, 35, 37), seem to be cross-reactive
(5), and have less stringent maintenance requirements
(38).
Intuitively, this result can be understood in terms of the influx of
primed naive cells with long telomeres into the memory compartment. If
memory T cells are dividing more frequently than naive T cells their
average telomere length will initially shorten more rapidly than that
of the naive T cells. The average memory telomere length cannot "run
freely," however, because it is bound to the average naive telomere
length by the influx of primed naive cells. This bond ultimately leads
to a steady distance
between the two average telomere lengths. It
is also natural therefore that this distance
depends on the
difference in naive and memory T cell division rates, on the clonal
expansion, and on the ratio of total numbers of naive and memory T
cells.
The effect of telomerase activity during clonal expansion on the
naive-memory distance
can be understood completely in terms of the
parameter K. One extreme case is "perfect telomerase
compensation" corresponding to K = 0. By Equation 15
one sees that K = 0 does allow for
> 0, i.e., a
fixed telomere length difference
between naive and memory T cells
remains a natural result even if there is no telomere shortening at all
during clonal expansion. This does require that memory cells divide
more frequently than naive T cells, i.e.,
M >
V, which is indeed supported by data (1, 2, 9, 10, 35, 37, 38). This resolves the apparent contradiction
between transient telomerase activity during clonal expansion
(23, 27) and the measured telomere shortening: the latter
could, in principle, all be due to telomere loss during renewal.
Obviously, the parameter K need not be zero, allowing the
1400-base pair difference (14, 15) to reflect both clonal
expansion and differences in renewal rates.
Parameter estimation
The time constant,
, and the priming rate of naive T cells,
, can both be estimated from Equation 11
. Arguing from the viewpoint
that in adults memory T cells are maintained largely by renewal, and
only marginally by the priming of naive T cells (1, 2, 6, 8, 35, 37), we obtain from Equation 11
that
CV <<
MM. Because in most adults the ratio of memory over
naive T cells is approximately one (39, 40), i.e.,
M/V
1, this simplifies into
C <<
M. By the same reasoning the time constant
simplifies to
1/(
C). Employing the estimated
memory division frequency of once every 22 wk, i.e.,
M
0.006 (10), we therefore
obtain that
> 22 wk. Approaching the steady state telomere
distance
in Equation 15
may therefore take several years. Since the
experimental data were all from human adults (14, 15) such
a time scale could indeed be consistent with our steady state
approximation. Assuming that T cell precursor frequencies generally
increase about a 1000-fold following a typical immune reaction
(41), i.e., assuming C
1000, we obtain
for the priming rate that
<< 6 x 10-6 per
day. In other words, on average, a naive T cell would run a chance of
less than 10-6 day-1 of being primed by a
foreign Ag to clonally expand.
| Discussion |
|---|
|
|
|---|
The naive and memory T cell repertoires in our model are maintained by proliferative self-renewal, and we have shown that this is consistent with the available data. Alternatively, one could argue that both naive and memory T cells are long lived, and are only stimulated to divide by clonal expansion upon stimulation with specific Ag. The latter view is, in fact, supported by the slow recovery rate of memory T cells in adults, and the relationship of this recovery rate with the rate of naive T cell recovery (11, 13). Furthermore, it has recently been suggested that the so-called "memory" isoforms of the CD45R T cell markers reflect an activation stage rather then a true memory stage (36). Thus, in the absence of Ag, memory cells would revert to a quiescent state with a naive CD45R phenotype (9, 36). These revertant T cells could, in principle, account for the observed shortening of the average telomere length of naive CD4+ T cells (14) in the absence of any cell division in the naive T cell compartment. The shortening of the average memory T cell telomere length could, indeed, be due to incomplete telomerase compensation during clonal expansion. Although it remains unclear whether this alternative scenario is realistic, we are currently developing alternative mathematic models to investigate its consistency.
In our model, we have not considered a possible source of naive T cells
from thymic or extrathymic lymphoid compartment. It is, indeed,
believed that such a source is small in human adults (8, 11). Our approach of a two-compartment model can, however, be
reiterated for the progenitor T cells providing novel naive T cells. If
the compartment of mature immunocompetent T cells would have a source
term from a self-renewing stem progenitor compartment
(20), we could reinterpret our two-compartment model for
progenitor T cell and mature T cell renewal. Calling progenitor T cells
vi and mature T cells mi,
C would represent the differentiation and proliferation of a
progenitor cell into an immunocompetent mature T cell. Thus, the rate
of mature T cell telomere length shortening should reflect twice the
division rate of the progenitor cells, and should be quite independent
of the mature T cell division rate.
Other groups (R. Antia and S. Frost, personal communication) are analyzing the same telomere data (14) with mathematic models incorporating a Hayflick limit (28). These authors interpret the data as showing a skew in the distribution of memory T cell telomere lengths in older humans, which would be due to the death of memory T cells when their telomeres become too short. Thus, memory T cells can be dividing more frequently than the rate estimated from the average shortening rate of their telomeres. Although the recent data on T cell telomerase activity cast some doubt on the relevance of the Hayflick limit for T cells, this remains a valid alternative interpretation of the data. By choosing to ignore the Hayflick limit, we have been able to demonstrate that the average rate at which memory T cells lose their telomere ends is expected to be independent of their division rate anyway. Our model can in principle be extended to include the effect of a Hayflick limit. The major effect of this is that the difference between the memory and naive telomere lengths need not be in perfect equilibrium, and may become smaller in elderly individuals. Thus, qualitatively, this hardly affects our results.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Address correspondence and reprint requests to Rob J. De Boer, Theoretical Biology U.U., Padualaan 8, 3584 CH Utrecht, The Netherlands. E-mail address: ![]()
3 Abbreviation used in this paper: MPD, mean population doubling. ![]()
4 K. C. Wolthers, A. J. Noest, S. A. Otto, F. Miedema, and R. J. De Boer. Mathematical modelling of CD4+ T cell subset TRF length in HIV-1 infection allows for only a marginally increased T cell turnover. Submitted for publication. ![]()
Received for publication September 8, 1997. Accepted for publication February 18, 1998.
| References |
|---|
|
|
|---|
secreting TH1 cells in humans. J. Clin. Invest. 99:2225.[Medline]
This article has been cited by other articles:
![]() |
N. Vrisekoop, R. van Gent, A. B. de Boer, S. A. Otto, J. C. C. Borleffs, R. Steingrover, J. M. Prins, T. W. Kuijpers, T. F. W. Wolfs, S. P. M. Geelen, et al. Restoration of the CD4 T Cell Compartment after Long-Term Highly Active Antiretroviral Therapy without Phenotypical Signs of Accelerated Immunological Aging J. Immunol., July 15, 2008; 181(2): 1573 - 1581. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lundegaard, O. Lund, C. Kesmir, S. Brunak, and M. Nielsen Modeling the adaptive immune system: predictions and simulations Bioinformatics, December 15, 2007; 23(24): 3265 - 3275. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Ke, J. A. Borghans, R. J. de Boer;, T. P. Arstila, A. Casrouge, V. Baron, J. Even, J. Kanellopoulos, and P. Kourilsky; Diversity of Human T Cell Receptors Science, May 19, 2000; 288(5469): 1135a - 1135. [Full Text] |
||||
![]() |
S. Bonhoeffer, H. Mohri, D. Ho, and A. S. Perelson Quantification of Cell Turnover Kinetics Using 5-Bromo-2'-deoxyuridine1 J. Immunol., May 15, 2000; 164(10): 5049 - 5054. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Weipoltshammer, C. Schofer, M. Almeder, V. V. Philimonenko, K. Frei, F. Wachtler, and P. Hozak Intranuclear Anchoring of Repetitive DNA Sequences: Centromeres, Telomeres, and Ribosomal DNA J. Cell Biol., December 27, 1999; 147(7): 1409 - 1418. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |