|
|
||||||||
Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI 48109.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|

T lymphocytes. T cell generation
through this thymic-dependent pathway is essential because it
establishes and helps to maintain the peripheral T cell pool. This
process also provides T cell repertoire diversity necessary for
potential immune responses to a vast number of Ags (1).
With increasing age, the lymphatic thymic mass decreases, and thymocyte
production correspondingly declines (thymic involution)
(2). Recent data suggest, however, that the adult thymus
can remain active even late in life supplying functional T cells to the
periphery (3, 4). CD4+ and CD8+ T cells progress through several stages during their long lifespan. Both mature CD4+ and CD8+ thymocytes first emigrate from the thymus into the periphery as recent thymic emigrants (RTE).3 After RTE further mature, they are classified as naive T cells that continually circulate through blood and lymphoid tissues. Measuring the number of RTE produced per day by the thymus allows for assessment of the thymic contribution to the peripheral T cell pool. It is especially important to monitor thymic production under conditions that influence T cell depletion and reconstitution, such as HIV-1 infection, bone marrow transplantation, and immunosuppressive therapy.
In chicken, mouse, and rat models, RTE have been identified using various markers. In chickens, chT1+ T cells are classified as RTE, and the level of chT1+ T cells has been shown to decline in parallel with thymic involution. These cells are distributed in peripheral blood and lymphoid tissues with a half-life of 3 days (5). In mice, FITC labeling of thymocytes and grafting of thymic lobes from congenic Ly5.2+ mice have been used to study RTE (6, 7). Most mouse RTE survive in the periphery for 23 wk after release from the thymus. It is not clear whether these cells die or undergo further maturation. In rats, Thy1+CD45RC-RT6- are markers for RTE (8). These RTE were observed to differentiate into more mature T cells within 1 wk after emigrating from the thymus.
Due to the lack of phenotypic markers for RTE in humans, it is difficult to distinguish them from long-lived naive T cells (5, 9). Whether human RTE are fully functional or require maturation in the periphery needs further study. Therefore, direct quantification and characterization of RTE in humans are not yet possible. Currently, the term RTE in humans generally refers to T cells that have undergone only a few cellular divisions after leaving the thymus (10).
Recently, the concentration of T cell receptor excision circles (TREC)
has been suggested as a measure for quantifying thymic output in humans
(9, 10). During TCR gene rearrangement, excised DNA
fragments are maintained in cells episomally as TREC. Two TREC species,
signal-joint TREC (sjTREC) and coding-joint TREC (cj TREC), which are
products of the deletion of the TCR
locus from the TCR
locus
during the first TCR
gene rearrangement, are produced sequentially
by all 
T cells. In most of them, these TREC are generated
without prior
rearrangement and are hence identical and can be
detected in
70% of 
T cells (9). A maximum of
two sjTREC and two cjTREC can be present within one 
T cell if
the corresponding rearrangement event occurs in both alleles. TREC are
exported from the thymus to the periphery within RTE. Thus, TREC levels
in the periphery could reflect RTE numbers. TREC are stable and are not
duplicated during mitosis; therefore, TREC concentration is diluted out
with each cell division. This explains why thymocytes have higher TREC
concentrations as compared with peripheral blood T cells (11, 12), and why naive T cells have higher TREC concentrations than
memory T cells (10).
TREC concentration has been widely used as a measurement of the number of RTE during HIV-1 infection and treatment, hemopoietic stem cell transplantation, and thymectomy (9, 10, 13, 14, 15, 16). However, controversy exists as to whether TREC concentrations are a good marker for RTE, because TREC concentrations are also affected by peripheral T cell turnover events, such as T cell division and death (9, 10, 17). Different techniques are used to measure TREC concentrations, such as real-time PCR (10, 17), quantitative-competitive-PCR (9, 13), and PCR-ELISA (18). Measurement units of TREC concentration vary as well, including TREC per 106 PBMCs (10), TREC per CD45RA+ T cells (17), TREC per microgram DNA of T cells (9, 14), and TREC per 105 CD4+ T cells (13). Together, these factors make it difficult to interpret and compare TREC data between studies.
It is not yet possible to experimentally differentiate the effects of thymic output and T cell turnover on TREC concentration in a quantitative way, because thymic output cannot be measured. To this end, we use mathematical modeling to predict whether TREC concentrations are a good marker of thymic function, as represented by the number of RTE. Our model captures age-dependent changes in thymopoiesis, RTE, and TREC levels. Using uncertainty and sensitivity analyses, we quantify the potential roles of thymic output, T cell division, T cell death, and intracellular TREC degradation in affecting TREC concentration at each specific age. We further define elements that contribute most to changes in TREC concentration that occur during aging and HIV-1 infection.
| Human Thymopoiesis Model |
|---|
|
|
|---|
|
Once model equations are derived that describe Fig. 1
(see
Appendix), we then solve the system numerically. We present
graphs for the corresponding simulations of this virtual model of human
thymopoiesis. To build the model, we require a representation for
thymic epithelial space (TES), which is the functional thymic tissue
containing all thymocytes. Fig. 2
A shows data for the maximal
number of thymocytes calculated from the volume of the TES region
(solid dots) and our best fit function for these data (solid line)
during an 80-year lifespan (see Appendix).
|
5% per year (Fig. 2
The model predicts the observed ratio of SP4 to SP8 cells as 2,
consistent with human data (19, 20). This ratio in the
thymus determines the CD4+:CD8+ T cell ratio in
the blood, which ranges from 1.5 to 2.5 (21). Our model
predicts that the maximum number of thymic emigrants at the age of one
is 1.1 x 109 cells/day (Fig. 2
D), which is
comparable with another estimate of human thymic emigrants of
109/day (22). Our virtual model also predicts
that the ratio of thymic emigrants per day to total thymocyte numbers
is 1.5%. This is consistent with experimental data ranging between 1
and 5% in mice (23). Thus, this model likely reflects the
process of human thymopoiesis. More importantly, we quantify the number
of RTE entering the periphery per day during the normal aging process,
which cannot be measured experimentally (Fig. 2
C).
| Human TREC Model |
|---|
|
|
|---|
|
The dynamic changes of CD4+ and CD8+ T
cells in peripheral blood/LT based on our simulations are shown in Fig. 4
A. The initial increase of T
cells is due to the large output of newly generated T cells from the
thymus in the early years of life (see Fig. 2
C). Then, T
cell numbers continuously increase until age 30 years. After age 30
years, T cell numbers remain at a relatively stable level with a slow
decline over time, as suggested by data (26, 27). Our
simulation results of T cell dynamics are within the normal range
calculated from human data on cell numbers per mm3 of
blood, blood volume, and body weight (error bars in Fig. 4
A)
(26, 28, 29, 30, 31). Increased peripheral T cell division
maintains the 30-year increase and subsequent relative equilibrium in T
cell numbers (data not shown). This is consistent with data from two
animal models of older sooty mangabeys and mice that lack thymic
function, indicating that most T cell production occurred in the
periphery (27, 32).
|
12 months. Then total
sjTREC level decreases concomitant with the decline in number of RTE
produced per day.
We simulate sjTREC concentration within CD4+ T cells
by dividing total levels of sjTREC within CD4+ T cells with
total CD4+ T cells. In experimental measurements
(9), results were presented in units of sjTREC/1.5 x
105 T cells. Hence we multiply our sjTREC/T cell results by
1.5 x 105 (Fig. 4
C, solid line).
Similarly, we calculate the sjTREC concentration within
CD8+ T cells (Fig. 4
D, solid line). The sjTREC
concentrations for both CD4+ and CD8+ T cells
decline
2 logs during an 80-year period.
Douek et al. (9) reported that sjTREC concentrations
decline by 11.5 logs in both blood CD4+ and
CD8+ T cells during the lifetime of an individual (in units
of sjTREC/1.5 x 105 T cells). Zhang et al.
(10) measured sjTREC level (
1 circles) per
106 PBMCs and found a total decrease of 2 logs over
seven decades. To convert data from Ref. 10 into units of
sjTREC/1.5 x 105 T cells, we divide sjTREC/PBMCs with
the ratios: total lymphocytes/PBMCs and T cells/total lymphocytes over
time (26, 33, 34), and then multiply by 1.5 x
105 T cells. sjTREC data from Douek et al. and Zhang et al.
were both measured using blood samples. The number of T cells in blood
represents
2% of total T cells in the body (35, 36, 37).
We scale the data sets appropriately by assuming that T cells in humans
continuously circulate between blood and lymph, exchanging in the blood
48 times daily in healthy individuals (35), and that the
levels of sjTREC per T cell in blood are similar to sjTREC per T cell
in LT. Our model results of average sjTREC concentrations in blood/LT
during aging are comparable with two experimental data sets from Douek
et al. and Zhang et al. (Fig. 4
, C and D). The
slight differences are likely due to different PCR techniques, target
cell populations, and unit conversions used for data sets.
Elucidating principal factors contributing to TREC concentration in healthy people
The key to interpreting TREC data is to quantitatively distinguish between the effects of thymic output and T cell turnover on sjTREC concentration (38). Although T cell division and T cell death can be measured experimentally, RTE levels cannot. To explore how both thymic and peripheral events affect sjTREC concentration, we study the parameters that correspond to each of four events that determines the sjTREC concentration: thymic output; T cell division; T cell death; and sjTREC degradation. In this section, we address relative contributions of these four events to sjTREC concentration in healthy people, which will help interpret sjTREC data during T cell depletion conditions. To do this, we first examine how these four events affect sjTREC concentration at each specific age, and then we test how these four events combined induce sjTREC concentration decline during 80 years in healthy people.
sjTREC concentration dynamics at a given age. To address relative contributions to sjTREC concentration at a given age, we implement uncertainty and sensitivity analyses for four parameters: SP4 emigration rate; T cell division rate; T cell death rate; and sjTREC degradation rate (cf Ref. 39). Detailed methods for these analyses are discussed in the Appendix. We present results for one of the outcome variables, sjTREC concentration within CD4+ T cells.
We first perform the uncertainty analysis by varying each parameter
independently. Shown in Fig. 5
are the
changes in sjTREC concentration in response to individual variations of
each parameter. The solid line represents the median level of sjTREC
concentration when the parameter is given its baseline value. The range
between upper and lower bounds represents the 95% confidence interval
for the simulated median of sjTREC concentration. According to the 95%
confidence interval, thymic emigration, T cell division, and sjTREC
degradation affect sjTREC concentration to a greater extent than does T
cell death, and this is true at any age over 80 years of life
(p < 0.005). No significant differences are
found between the effects of thymic output and T cell division
(p > 0.05), and we assume the sjTREC half-life
(degradation rate) is likely constant. Therefore, our results indicate
that, at any age, when these events have similar magnitude changes,
thymic output and peripheral T cell division can have an equal impact
on sjTREC concentration, whereas T cell death affects it less.
|
The above experiments were performed by varying each of four events,
thymic output, T cell division, T cell death, and sjTREC degradation,
independently. What remains to be determined is the effect on sjTREC
concentration when all four events are allowed to vary in combination.
This sensitivity analysis will assess the importance of these four
events relative to each other for inducing variability in sjTREC
outcomes. Model simulations of sjTREC concentrations for nine virtual
patients show large variations and different shape curves (Fig. 6
) arising from uncertainties in the four
governing parameters. Differences in the values of four parameters may
explain the considerable discrepancy observed from patient to patient
and study to study (Fig. 4
, C and D) (9, 10, 24).
|
|
|
Our human thymopoiesis and TREC model provides a systematic way to quantitatively analyze how the thymus and peripheral T cell turnover contribute to changes in sjTREC concentration during T cell depletion situations, such as HIV-1 infection, hemopoietic stem cell transplantation preceded by chemoradiotherapy, thymectomy, and chemotherapy.
To apply our TREC model as a predictor of thymic function, we examine an HIV-1 infection scenario. HIV-1 disease is characterized by a gradual decline of CD4+ T cells, with an initial rise followed by a fall in CD8+ T cells during a 10-year period typically (40). sjTREC concentrations within both CD4+ and CD8+ T cells decrease during early HIV-1 infection (9, 13). Some studies suggest that the decline in sjTREC concentration results from a decrease in thymic output (9, 41). Evidence shows that HIV-1 interrupts thymopoiesis by directly infecting CD4+ thymocytes and stromal cells, inducing the decline of thymic output (42, 43, 44). Other studies indicate that global immune activation induced by HIV-1 leads to increased cell division, causing the dilution of sjTREC concentration (17, 25).
Experimental studies suggest that the division rate of CD4+
T cells in patients with HIV-1 infection increases
1- to 2-fold as
compared with uninfected controls, whereas the division rate for
CD8+ T cells increases 7- to 8-fold as compared with
uninfected controls (25, 45, 46, 47). Death rates of both
CD4+ and CD8+ T cells increase 3- to 4-fold
during HIV-1 infection (45, 46). No quantitative data have
been reported as to how RTE levels change during HIV-1 infection, due
to the lack of a phenotypic marker for human RTE. This is the basic
question that hinders distinguishing the relative role of the thymus in
HIV-1 infection. However, our model provides a method to address this
question.
As shown in Fig. 8
, we simulate changes
in sjTREC concentrations and T cell levels during early HIV-1 infection
in a virtual group of people age 30 years. To achieve dynamics similar
to those reported by clinical data for both T cells and sjTREC
concentrations, our model predicts that thymic output of
CD4+ RTEs must decline 10- to 15-fold, and thymic output of
CD8+ RTEs must decline 1- to 6-fold. Experimental data show
that HIV-1 directly infects CD4+ thymocytes and causes the
SP4:SP8 ratio to invert from 2 to 0.5 (20, 42), suggesting
that the decrease in CD4+ RTE is greater than that
occurring in CD8+ RTE, which is consistent with our model
predictions.
|
| Discussion |
|---|
|
|
|---|
Characterization of human RTE has been hampered by the lack of a marker that distinguishes them from long-lived naive T cells. A number of phenotypic markers have been used for the evaluation of naive T cells, with no clear understanding of which naive phenotype is best enriched for by RTE, including CD45RA+CD45RO-, CD45RA+CD62L+, CD45RO-CD27+CD95low (48). One study using a phenotypic approach suggested that CD103-bearing naive CD8+ T cells are a subpopulation of human RTE (12). Another study indicated that cord blood naive T cells are human RTE (49). These identified RTE have rapid rates of apoptosis like thymocytes, and the presence of IL-7 maintains their survival and expansion in an Ag-independent manner.
Our first objective was to quantify age-dependent changes in total
numbers of human RTE produced per day, which cannot be experimentally
measured. We develop a model of human thymopoiesis to capture thymic
involution. This is the first model for studying human thymopoiesis and
quantifying the number of RTE exported per day, and our model is
developed based mainly on human experimental studies. Our results
indicate that thymic output per day decreases exponentially by two
orders of magnitude during 80 years (Fig. 2
D). These are
useful data for monitoring thymic function in healthy individuals, as
well as providing a baseline comparison of thymic output during disease
states.
We quantify thymic output as a rate measurement of the number of RTE produced per day. Currently, the longevity and migration patterns of human RTE are not known, rendering the calculation of total numbers of RTE in the peripheral blood/LT difficult. However, total RTE levels are mainly affected by thymic output per day. Thus, monitoring thymic output per day is a more direct and accurate estimation of thymic function as compared with total RTE numbers.
Recently, TREC concentration has been widely used to measure thymic output. Data suggest that TREC levels change during aging and disease (9, 10, 13, 14, 15, 17, 24, 41). However, it is still an open question as to whether TREC concentrations are mainly affected by thymic output or peripheral T cell events. Different techniques have been used to measure TREC levels, and the unit for measuring TREC levels varies from study to study, which can affect interpretation of TREC data (9, 10, 17, 18).
Our second aim was to verify whether TREC concentration is an accurate predictor of thymic function, represented by the number of RTE produced per day during aging and HIV-1 infection. We first simulate sjTREC concentrations as they change with age and validate our model with two experiment data sets. Using uncertainty and sensitivity analyses, we next quantify how key processes potentially control the dynamics of sjTREC concentration. Our results strongly suggest that thymic output and peripheral T cell division could equally affect sjTREC concentration at any age. T cell death contributes less to changes in sjTREC concentration than do thymic output and T cell division. We further demonstrate that sjTREC concentration is a good measurement for both CD4+ and CD8+ RTE in healthy people and for CD4+ RTE during HIV-1 infection. Our findings have implications for interpreting TREC data during other T cell depletion situations, such as hemopoietic stem cell transplantation, thymectomy, and chemotherapy. We propose that peripheral T cell division and death should be examined before considering TREC concentration as a representation of thymic function.
A previous modeling study (17) suggested that only T cell
division determines changes in TREC concentrations during aging and
HIV-1 infection. This result is based on the assumption that neither
division of naive T cells nor intracellular degradation of TREC occurs.
Although division rates of naive cells are relatively low, they are
unlikely to be zero. Their results later show that naive T cell
division does occur. Similarly, although TREC are stably maintained
within T cells, they have a finite half-life, likely shorter than or
close to the half-life of T cells; TREC in chickens have been shown to
have a half-life of
2 wk (50). Our model predicts that
sjTREC degradation rate is 0.002/day, similar to the T cell death rate,
for example. The previous study (17) further assumes that
naive T cell numbers are proportional to thymic output, and thus TREC
concentration is not a measurement of thymic output but of T cell
division. However, data suggest that naive T cell numbers and thymic
output are not proportional, given that a subset of naive T cells have
a very long lifespan and a subset of memory T cells can revert back to
naive phenotype (22).
Recent clinical data support the thymus as a key contributor to TREC concentration decline during aging and HIV-1 infection (38, 41). One study suggests that changes in peripheral T cell division do not adequately explain many observations of TREC concentrations, including similar ratios of CD4+ TREC to CD8+ TREC in HIV-1-infected patients and healthy controls, similar TREC responses for discordant responders and good responders to antiretroviral therapy, and increased TREC concentration despite increased proliferation of T cells immediately after antiretroviral therapy (41). Another study indicates that the TREC concentration indeed reflects thymic output during aging (38). This study further suggests that thymic output most likely affects TREC concentration within CD4+ T cells and that a combination of increased T cell division and decreased thymic output reduces TREC concentration within CD8+ T cells after HIV-1 infection (38). These results are consistent with our model predications for TREC concentrations during both aging and HIV-1 infection scenarios.
Almost all TREC studies in humans are performed using peripheral blood.
However, the majority of T cells (
98%) exist in the peripheral
lymphoid system (35, 36, 37). Thus, quantification of TREC
levels in different LT compartments would contribute greatly to our
understanding of TREC dynamics. One study using a primate animal model
observed that CD4+ and CD8+ T cells present in
the lymph nodes contain more sjTREC than do peripheral blood T cells,
suggesting that RTE can home into lymphoid tissues (11).
Another study using a rat model observed that the percentages of RTE
among T cell population were comparable in blood and all other LT
compartments, indicating that rat RTE continuously migrate through
blood and lymphoid organs as naive T cells (51). Our TREC
model presents average sjTREC concentrations in peripheral blood/LT
assuming that the sjTREC concentration in T cells within LT is similar
to those in peripheral blood T cells. If the results from the primate
animal model is applicable to humans, then our model simulations of the
sjTREC concentration represent a lower boundary approximation.
Quantification of TREC concentration in two distinct pools of memory and naive T cells is an important next step. TREC concentration has been measured in the naive T cell pool (17). The differences in turnover rates of these T cell subclasses can further elucidate TREC dynamics. Currently, most experiments on TREC concentration are measured using CD4+ and CD8+ T cells (9, 11, 12, 13, 18, 38, 41, 48). Therefore, we can accurately build a model for TREC dynamics in these two compartments. Our model can be further adapted to explore TREC dynamics in memory and naive T cell pools when more data are available.
Overall, this study characterizes the relative roles of thymic output and periphery T cell turnover in regulating TREC concentrations. We propose that peripheral T cell division and death should be examined together with TREC concentration. Our model can be used as an integrated system for testing whether TREC concentration is an accurate marker for RTE levels in all situations by incorporating both T cell dynamics and TREC concentrations that are experimentally measured.
|
|
| Appendix 1 |
|---|
|
|
|---|
Defining mathematical notation for the terms shown in Fig. 1
, T(t) represents TN cells at time t, I(t)
represents ITTP cells, D(t) represents DP cells,
S4(t) represents SP4 cells and
S8(t) represents SP8 cells. The mathematical
representation describing the interaction of these five thymocyte
subsets is based on a previous model of mouse thymocyte subsets
(52 ). However, the thymopoiesis process is different in
humans and mice (compare Refs. 19 and 23 ),
and our model reflects these differences. We suppress the time
dependence, (t), in the variables for ease of readability,
except where needed for emphasis:
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
Human TREC model
TREC concentrations are determined by total TREC dynamics and
total T cell dynamics. Define T4 and
T8 to represent total CD4+ and
CD8+ T cells in both blood/LT, respectively, and
4 and
8 to represent total sjTREC levels
within CD4+ and CD8+ T cells in both blood/LT,
respectively. Then a model for sjTREC and T cell dynamics based on Fig. 3
is given by
![]() | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
1
[k4(t)T4/(T4 +
S1)],
2
[k8(t)T8/(T8 +
S2)]) helps to maintain peripheral T cell homeostasis
during human growth, as a compensation to the involution of thymic
function (27 32 ). Death of T cells
(
1
[T4/k4(t)]T4,
2
[T8/k8(t)]T8) occurs in a
density-dependent fashion. In division and death terms,
k4(t) and k8(t) represent
the average numbers of total CD4+ and total
CD8+ T cells in the peripheral blood/LT, respectively. Here
we model T cell dynamics for healthy individuals during aging; thus T
cell growth and death are not further induced by Ag and inflammatory
signals.
Total sjTREC levels are mainly affected by three processes: input from
the thymus (c1e1S4,
c2e2S8), loss due to death of
T cells (
1
[T4/k4(t)]
4,
2[T8/k8(t)]
8),
and loss due to degradation of sjTREC within T cells
(
1
4,
2
8).
The loss of TREC due to T cell death is proportional to T cell death.
We do not include T cell division and differentiation terms, because
these two processes do not affect total sjTREC levels. With dynamics of
both total sjTREC and total T cells, TREC concentrations can be
simulated by dividing sjTREC by T cells.
Parameter estimation
The variables and parameters used in the human thymopoiesis
model Equations 15![]()
![]()
![]()
![]()
are defined in Table II
. The parameters are the rate
coefficients for each process described. Values for variables and rates
were estimated from the literature. We give weight to human data,
followed by data generated by murine and in vitro systems. For rates
for which poor or no data exist, we perform uncertainty and sensitivity
analyses to explore the change in outcomes over a full range of
possible values. We outline below how we estimated the parameters.
Although the thymus undergoes involution during aging, the distribution
of thymocyte subsets that is seen in adults is similar to the
distribution observed in fetuses (3 ). TREC levels per
thymocyte in adults also appear to be the same in fetuses and newborns
(3 ). Thymic stromal cells, which are associated with
thymocyte development, remain present and active throughout life
(53 ). These data suggest that the effect of aging on the
thymus appears to be quantitative rather than qualitative. To describe
thymic involution dynamics, we represent the maximal number of
thymocytes in TES region, Tes(t), and the bone
marrow cell source, s(t), as time-dependent functions:
![]() |
![]() | (10) |
![]() | (11) |
The lymphatic thymic tissue is composed of TES and perivascular
space (2 ). The TES region, where thymopoiesis occurs,
contains all of the thymocytes. A precise correlation has been shown
between the number of RTE and functional thymic tissue in chickens
(5 ). Here we assume that the volume of the TES region
represents the maximum allowable number of thymocytes, i.e., the
carrying capacity of the thymus for total thymocytes. The TES region
reaches its highest level of 21.8 cm3 by the age of one
(2 ), containing maximally 1011 thymocytes
(54 ). We use the ratio 1011 thymocytes/21.8
cm3 TES volume to calculate the maximal number of
thymocytes from the TES volume (Fig. 2
A). The values for the
TES volume are obtained by subtracting the mean volume of perivascular
space from the mean volume of lymphatic thymic tissue by corresponding
age groups from a stereological evaluation of 136 human thymuses
(2 ). The maximal number of thymocytes,
Tes(t), can be obtained by a best fit function
using nonlinear least squares (Equation 10
for 1 year <
t). The thymic involution rate is 5.66% per year after the
first year of life. During the first year of life, the function
defining the maximal number of thymocytes is derived from the
calculation of thymus weight (Equation 10
for 0
t
1 year) (55 ).
Bone marrow shows no loss of function with age (54 );
however, the recovery of normal T cell subsets after bone marrow
transplantation depends on thymic function (56 57 ). On
the basis of this evidence, we assume that the number of source cells
from bone marrow utilized by thethymus is determined by the number of
available thymocytes. Thus, we track this change in cell source,
s(t) (Equation 11
), as a function of
Tes(t). The seeding rate of progenitor cells
into TN cells is 5 x 104 per day in young mice
(23 ). Because the human thymus is 3 to 4 orders of
magnitude larger than the mouse thymus, 5 x 108 per
day is used as the cell source for TN cells at age 1 year.
The initial value of total thymocytes is 75% of maximal number of thymocytes, implying that one-third of total thymocytes are produced each day. This is derived from a calculation of growth rates, which is consistent with experimental reports (23 58 ). The percentages of TN, ITTP, DP, SP4, and SP8 cells are 3.5, 2.5, 73, 14, and 7%, respectively, from studies of human thymic specimens (3 19 ). Thus, the initial values for the five thymocyte subsets are calculated by multiplying 75% of maximal number of thymocytes at birth by the corresponding subset percentage. This allows for up to a 25% increase in thymocyte number.
Growth rates of thymocytes have been measured experimentally in the range of 0.51.5/day in mice (23 ). Generally, immature thymocytes have greater potential for expansion, whereas mature SP cells grow relatively slowly (23 ). We estimate growth rates of TN, ITTP, DP, SP4, and SP8 cells to be 1.5/day, 1.0/day, 1.5/day, 0.5/day, and 0.5/day, respectively.
It takes 3 days on average for TN cells to differentiate into ITTP cells (19 ); thus, the differentiation rate for TN cells is 0.3/day. Similarly, the differentiation rate for ITTP cells is calculated as 0.5/day (19 ). Approximately 5% of DP cells survive positive and negative selection to become SP cells (59 ), and the ratio of SP4/SP8 equals 2 (3 ); therefore, we estimate the differentiation rates of SP4 and SP8 cells to be 0.04/day and 0.02/day, respectively. Newly formed SP cells spend an average of 14 days in the thymus before emigration into the periphery as observed in a murine thymus system (7 ). Hence emigration rates of SP4 and SP8 cells are estimated as 0.07/day.
Although no available data exist for death rates of subsets of human
thymocytes, the average lifespan of mouse thymocytes has been estimated
to be 3 days (23 ). If we assume that the thymus is in
steady state over a short time, then the maximal number of thymocytes
and cell source are constant and Equations 15![]()
![]()
![]()
![]()
can be set to 0 (i.e.,
no change in the cell rates). Based on that assumption, the death rates
for the five subsets of thymocytes can be calculated to be 0.27/day,
0.17/day, 0.33/day, 0.26/day, and 0.26/day, for TN, ITTP, DP, SP4, and
SP8 cells, respectively.
The variables and parameters used in the human TREC model Equations 69![]()
![]()
![]()
are defined in Table III
. The
average numbers of CD4+ and CD8+ T cells in the
blood linearly increase during the first 30 years of life and then
reach a relative steady state, considering age-dependent changes in
cell number per mm3 blood, blood volume, and body weight
(26 28 29 30 31 ). We assume that the number of T cells in the
LT have a pattern similar to that in the blood during human growth.
Thus, the functions for the average numbers of total CD4+
and total CD8+ T cells in the peripheral blood/LT are as
follows:
![]() | (12) |
![]() | (13) |
There were approximately 200,000 cjTREC per 150,000 CD4+ and CD8+ T cells (100% naive type) observed in cord blood (9 ). Three to four cell divisions occur between rearrangements that produce sjTREC and cjTREC (9 ). Thus, on average 200,000/23.5 = 17677.67 sjTREC exist within 150,000 CD4+ and CD8+ T cells in cord blood. The concentration of sjTREC within cord blood T cells is calculated to be 17677.67/150,000 = 0.118/cell. The sjTREC concentration in cord blood CD4+ T cells is similar to that in CD8+ T cells (9 ). Another study observed that cord blood T cells have high levels of TREC compared with adult naive T cells, indicating that cord blood T cells do represent RTE (49 ). Thus, we use 0.118 sjTREC/cell for concentrations of sjTREC within CD4+ and CD8+ RTE.
The degradation rate of TREC in chickens is 0.05/day, calculated from a half-life of 2 wk (50 ). It is not known how long it takes for sjTREC to degrade in human T cells. We hypothesize that degradation of sjTREC in humans is slower than that in chickens, and the lifespan of sjTREC is less than or equal to that of T cells. Thus, we implement an uncertainty analysis by varying the sjTREC degradation rates for CD4+ and CD8+ T cells. Finally, we use the value of 0.002/day as the sjTREC degradation rate.
The initial conditions (at birth) for total CD4+ and CD8+ T cells are 4.33 x 1010 and 1.80 x 1010, respectively (28 29 30 ). The initial condition for total sjTREC levels within CD4+ T cells is 5.109 x 109, calculated by multiplying total CD4+ T cell number at birth by 0.118 sjTREC/cell (TREC concentration in RTE). Similarly, 2.124 x 109 has been calculated as the initial condition for total sjTREC level within CD8+ T cells.
Once the model equations are derived and the parameter values are estimated, we then solve the system of nonlinear ordinary differential equations using an appropriate numerical method. We solve the system in two separate computational programs to ensure accuracy.
Uncertainty and sensitivity analysis
We use a Latin Hypercube Sampling scheme and PRCC (cf Ref. 39 ) for uncertainty and sensitivity analyses, respectively. Briefly, for the uncertainty analysis, a random sample of each of the parameters to be tested is generated from a list of parameter ranges and distributions. In our case, this sampling is done 90 times. The values in each column of the matrix generated from the sampling are randomly chosen values for a given parameter from a widely defined, evenly partitioned range of values. The sets of values in each row of the matrix are then used in 90 independent simulations. These values provide a range of results for each solution curve generated from the original parameter set. The sensitivity analysis predicts which parameter(s) have the greatest influence on the outcomes resulting from the uncertainty analysis. In our case, by comparing the values of PRCC, the relative importance of thymic and peripheral T cell events to TREC concentrations can be quantitatively evaluated.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Address correspondence and reprint requests to Dr. Denise E. Kirschner, Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI 48109-0620. E-mail address: kirschne{at}umich.edu ![]()
3 Abbreviations used in this paper: RTE, recent thymic emigrants; TREC, TCR excision circles; sjTREC, signal-joint TREC; cjTREC, coding-joint TREC; TN, CD3-CD4-CD8- triple-negative; ITTP, CD3-CD4+CD8- intrathymic T progenitor; DP, CD3+CD4+CD8+ double-positive; SP4, CD3+CD4+CD8- single-positive; SP8, CD3+CD4-CD8+ single positive; TES, thymic epithelial space; LT, lymphoid tissues; PRCC, partial rank correlation coefficients. ![]()
Received for publication December 18, 2001. Accepted for publication March 11, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Ghattas, D. L Wallace, J. A Solon, S. M Henson, Y. Zhang, P. T Ngom, R. Aspinall, G. Morgan, G. E Griffin, A. M Prentice, et al. Long-term effects of perinatal nutrition on T lymphocyte kinetics in young Gambian men Am. J. Clinical Nutrition, February 1, 2007; 85(2): 480 - 487. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. van den Dool and R. J. de Boer The Effects of Age, Thymectomy, and HIV Infection on {alpha} and beta TCR Excision Circles in Naive T Cells J. Immunol., October 1, 2006; 177(7): 4391 - 4401. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Thiebot, B. Vaslin, S. Derdouch, J.-M. Bertho, F. Mouthon, S. Prost, G. Gras, P. Ducouret, D. Dormont, and R. Le Grand Impact of bone marrow hematopoiesis failure on T-cell generation during pathogenic simian immunodeficiency virus infection in macaques Blood, March 15, 2005; 105(6): 2403 - 2409. [Abstract] [Full Text] [PDF] |
||||
![]() |
C Kayser, F L Alberto, N P da Silva, and L E. Andrade Decreased number of T cells bearing TCR rearrangement excision circles (TREC) in active recent onset systemic lupus erythematosus Lupus, December 1, 2004; 13(12): 906 - 911. [Abstract] [PDF] |
||||
![]() |
Y.-W. Chu, S. A. Memon, S. O. Sharrow, F. T. Hakim, M. Eckhaus, P. J. Lucas, and R. E. Gress Exogenous IL-7 increases recent thymic emigrants in peripheral lymphoid tissue without enhanced thymic function Blood, August 15, 2004; 104(4): 1110 - 1119. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. L. Staton, B. Johnston, E. C. Butcher, and D. J. Campbell Murine CD8+ Recent Thymic Emigrants are {alpha}E Integrin-Positive and CC Chemokine Ligand 25 Responsive J. Immunol., June 15, 2004; 172(12): 7282 - 7288. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Permar, W. J. Moss, J. J. Ryon, D. C. Douek, M. Monze, and D. E. Griffin Increased Thymic Output during Acute Measles Virus Infection J. Virol., July 15, 2003; 77(14): 7872 - 7879. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |