Nave T-Cell Dynamics in Human Immunodeficiency Vir
http://www.100md.com
病菌学杂志 2006年第6期
National Institutes of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, Maryland
SAIC, Frederick, Maryland
Clinical Center, National Institutes of Health, Bethesda, Maryland
National Cancer Institute-Frederick, National Institutes of Health, Frederick, Maryland
ABSTRACT
Both nave CD4+ and nave CD8+ T cells are depleted in individuals with human immunodeficiency virus type 1 (HIV-1) infection by unknown mechanisms. Analysis of their dynamics prior to and after highly active antiretroviral therapy (HAART) could reveal possible mechanisms of depletion. Twenty patients were evaluated with immunophenotyping, intracellular Ki67 staining, T-cell receptor excision circle (TREC) quantitation in sorted CD4 and CD8 cells, and thymic computed tomography scans prior to and 6 and 18 months after initiation of HAART. Nave T-cell proliferation decreased significantly during the first 6 months of therapy (P < 0.01) followed by a slower decline. Thymic indices did not change significantly over time. At baseline, nave CD4+ T-cell numbers were lower than naive CD8+ T-cell numbers; after HAART, a greater increase in nave CD4+ T cells than nave CD8+ T cells was observed. A greater relative change (n-fold) in the number of TREC+ T cells/μl than in nave T-cell counts was observed at 6 months for both CD4+ (median relative change [n-fold] of 2.2 and 1.7, respectively; P < 0.01) and CD8+ T cell pools (1.4 and 1.2; P < 0.01). A more pronounced decrease in the proliferation than the disappearance rate of nave T cells after HAART was observed in a second group of six HIV-1-infected patients studied by in vivo pulse labeling with bromodeoxyuridine. These observations are consistent with a mathematical model where the HIV-1-induced increase in proliferation of nave T cells is mostly explained by a faster recruitment into memory cells.
INTRODUCTION
The mechanisms leading to CD4+ T-cell depletion during human immunodeficiency virus (HIV) infection and CD4+ T-cell restoration during therapy with highly active antiretroviral therapy (HAART) remain undefined. Direct viral cytopathogenicity or redistribution of lymphocytes between lymphoid tissues and the circulation can account only in part for these changes (28, 32, 41). Characterizing lymphocyte turnover in HIV-infected patients can potentially lead to a better understanding of these mechanisms.
A number of studies utilizing indirect methods (2, 21, 31, 40, 48) and, more recently, direct methods using deuterated glucose (25, 36) or 5-bromo-2'-deoxyuridine (BrdU) (29) to measure CD4 cell dynamics have shown that CD4 cell turnover is increased during chronic HIV type 1 (HIV-1) infection. Moreover, while preliminary cross-sectional studies described an increase in CD4 cell turnover in patients following initiation of HAART, suggesting a defect in CD4 cell production secondary to HIV infection (12, 25), longitudinal studies have clearly demonstrated a rapid and persistent decrease in CD4 cell proliferation following initiation of HAART, suggesting that the increase in CD4 cell turnover itself may be an important pathogenic mechanism of CD4 depletion (21, 29, 36). While this increased turnover was initially postulated to represent a homeostatic response to CD4 depletion (11, 12), such a hypothesis is inconsistent with the rapid reduction in proliferation of CD4+ as well as CD8+ T cells after viral suppression with HAART, prior to normalization of CD4 cell numbers (2, 29, 31, 36). These observations led to alternative hypotheses proposing that either HIV-directed or nonspecific immune activation drives increased turnover. Moreover, based in part on studies demonstrating that levels of immune activation in T cells, especially CD8 cells, are independent predictors of CD4 depletion and disease progression, immune activation is currently felt by many investigators to play a direct role in HIV-associated CD4 depletion (17, 34, 43).
An increase in turnover has been demonstrated in nave T cells (24, 26) as well as memory T cells during pathogenic lentiviral infection. Whereas memory CD4+ but not memory CD8+ T cells decrease in number during chronic HIV infection, both nave CD4+ as well as nave CD8+ T cells are depleted during such infection (5, 22, 39). This observation led to the conclusion that nave T-cell depletion is one of the hallmarks of HIV infection. While infection of nave T cells has been documented, this appears to be a relatively rare event that cannot quantitatively explain the loss of nave CD4+ T cells. The observation that both nave CD4+ and nave CD8+ T cells decrease during HIV infection led to the hypothesis that persistent hyperactivation of the immune system leads to erosion of nave T cells by their increased recruitment into memory cells (20), probably through antigen- and nonantigen-specific stimulation, as has been shown in animal models (18, 37).
Because the thymus is the source of new T cells, examining thymic function during HIV infection and therapy is critical to studies of T-cell dynamics (14, 19). Due to the difficulties in directly studying thymic function, quantitation of T-cell receptor excision circles (TRECs) has been utilized as a surrogate of thymic function. Early studies measuring the number of TRECs per nave T cell suggested that HIV infection leads to a decrease in thymic function and that improved thymic function contributes to immune reconstitution following HAART (9). More recently, studies using mathematical modeling showed that the observed changes in TREC content cannot be explained solely by changes in thymic function (21) or by redistribution of T lymphocytes from lymphoid tissues to the blood (32); the observed TREC dynamics were more consistent with changes in peripheral proliferation and disappearance rates of T-lymphocytes.
In the current study we undertook a detailed examination of the relationship between T-cell turnover, thymic function, and immune activation in HIV-1-infected patients to better understand the contribution of these various parameters to the immunologic changes seen during HIV infection and therapy.
MATERIALS AND METHODS
Patient characteristics. Twenty-three HIV-1-infected patients who were initiating HAART (n = 20), starting HAART after a 6- and 24-month treatment interruption (n = 2), or switching from an ineffective to an effective HAART regimen (n = 1) were studied. Twenty previously untreated (or with limited nucleoside analogue exposure) HIV-1-infected patients (group 1) enrolled between 1997 and 2000 in a study to identify viral reservoirs in patients receiving antiretroviral therapy. These patients initiated therapy with indinavir, nevirapine, lamivudine, and zidovudine (or stavudine [one patient]). Drug changes were permitted for intolerance; for 14 patients, zidovudine was replaced by stavudine, and for 6 patients, indinavir was replaced by nelfinavir. After at least 1 year on study, additional changes were permitted at the patient's and referring physician's discretion; seven patients made such changes, usually for simplification of the regimen. All patients ultimately demonstrated a virologic response, with declines in plasma viral loads to levels below branched DNA assay limits (<50 copies/ml, 18 patients; <500 copies/ml, 2 patients) and all but 2 patients had an increase in CD4 cell numbers of at least 100 cells/mm3. Group 1 patients had immunophenotyping analysis and evaluation of thymic function by TREC analysis and thymic computed tomography (CT) scans at the following time points: 0 (pre-HAART), 1 (5 to 8 months, median of 6 months after start of therapy), and 2 (13 to 42 months, median of 18 months after start of therapy). Nave cell labeling kinetics were analyzed for six HIV-1-infected patients (group 2) who had undergone in vivo pulse labeling with BrdU prior to and 3 to 6 months after the initiation of an effective HAART regimen. Three of these patients are also included in group 1; the remaining three received stavudine, lamivudine, indinavir, and nevirapine (one patient); stavudine, lamivudine, and efavirenz (1 patient); and abacavir, amprenavir, nelfinavir, and efavirenz (one patient). Results for CD4 and CD8 cell labeling, but not nave cell labeling, have been previously reported for this group (29). Baseline characteristics for group 1 and pre- and post-HAART parameters for group 2 are shown in Tables 1 and 2.
Viral load. HIV-1 RNA levels were determined using a modification of the Roche Amplicor HIV Monitor assay kit (Indianapolis, IN) (27) or a branched DNA assay (Bayer Diagnostics, Tarrytown, NY) (3). Each assay had a lower limit of detection of 50 copies/ml; for two patients enrolled early in the study, the limit of detection was 500 copies/ml.
Immunophenotyping and intracellular staining for Ki67. Immunophenotypic analysis of cryopreserved peripheral blood mononuclear cells was performed using four-color immunofluorescence as previously described (42). Nave cells were defined as CD45RO– CD27+, central memory were defined as CD45RO+ CD27+, and effector memory were defined as CD45RO+ CD27– for CD4 T cells and as CD27– (CD45RO+ or CD45RO–) for CD8 T cells. Cells were stained intracellularly with Ki67-phycoerythrin (clone B56) or isotype (mouse immunoglobulin G1-phycoerythrin; clone MOPC-21) from BD/Pharmingen. T-cell proliferation was defined as the percentage of cells expressing Ki67 (15).
TREC determination. Signal joint TRECs (Sj-TRECs) in purified cell subsets were quantitated by real-time PCR by the cell lysis method as described previously (38). The consistency of the DNA content of the cell lysate was checked by real-time PCR using a ribosomal protein gene and a TaqMan gene expression assay kit from Applied Biosystems, Inc. (Foster City, CA). Because no more than one Sj-TREC can be produced per cell, the number of TRECs per unit volume of blood also represents the number of TREC+ T cells in the same unit volume.
Thymic CT scans. CT scans of the thymus were obtained prior to and at a median of 6 and 18 months after starting HAART. Scans were graded as previously described on a 0 (no thymic tissue) to 5 (thymic mass) scale by two independent radiologists blinded to clinical and laboratory results (35). In addition, computer-based density and volume analysis of the thymus was performed by transferring CT data to a GE Advantage Windows workstation (versions 2.1 and 4.0; GE Medical Systems, Advanced Windows Workstation Training Program, Milwaukee, WI). Contours of the anterior mediastinum were outlined by a radiologist-trained technician and corroborated by a radiologist (13).
BrdU infusion and flow cytometry. The fractions of BrdU+ CD4+ CD45RO– and BrdU+ CD8+ CD45RO– T cells were analyzed by flow cytometry as previously described (29).
Statistics. Changes in median values for each variable were tested for significance by the permutation test with paired samples computed with an exact method (44). Tests were performed using StatXact software. Association between variables was determined by the Spearman rank correlation test. Adjustment of P values for multiple testing was done by the Bonferroni method. Occasional data points were missing for group 1; however, all paired analyses were tested with n values that were 15.
The relative change (n-fold) of a variable H between time (t) points (e.g., t0 and t1) is defined as the ratio between the value of H at t1 versus t0: H(t1)/H(t0).
Modeling. Differential equations were solved using Labview 7.0 (National Instruments, Austin, TX). The data were fitted to the differential equations using the Levenberg-Marquardt method (16).
Mathematical model for TREC analysis. To help interpret the data obtained in this study, we propose a slight generalization of a mathematical model originally described by Hazenberg et al. (21). The number of nave T cells/μl (T) and the number of TREC+ T cells/μl (T+) are governed by the following equations:
(1)
Nave CD4+ or CD8+ T cells, T, in the periphery receive a constant input from the thymus, , proliferate at rate p (day–1) and disappear at a rate d (day–1). TREC+ T cells, T+, appear at a lower rate of thymic production (f) and disappear from the same compartment as nave T cells at a rate d. We assume that nave T cells can proliferate without losing their nave phenotype as has been reported (45-47). Since TRECs do not replicate during cell mitosis (9), proliferation of TREC+ T cells decreases the fraction of TREC+ T cells per nave T cell (T+/T).
Since changes in T and T+ occur very slowly, the steady-state values obtained by equation 1 can be used to analyze how changes in the parameters (, d, p, and f) would affect changes in the number of nave cells, T, the number of TREC+ cells, T+, and the fraction of TREC+ T cells per nave T cell (T+/T):
(2)
For the generalization in equation 1, we assume that the disappearance rate, d, is a composite of two different factors: d = A + R, where A represents the death rate of nave cells and R is the rate of priming of nave cells into memory cells (7, 23). We also assume that the increase in the proliferation rate of nave T cells during chronic HIV-1 infection (p p + ) is largely explained by the increase in the rate of priming of nave T cells into memory cells (R R + ). Thus, the generalized model has been simulated to predict the changes in nave T-cell counts, TREC+ T cells, and the fraction of TRECs per nave T cell after HAART, when the perturbation of the quasi-steady state induced by the administration of HAART is described mathematically by the reduction of .
Mathematical model for the kinetics of BrdU-labeled nave T cells. To describe the in vivo kinetics of BrdU-labeled nave T cells in the blood after a 30-min BrdU infusion, we used the following semiempirical equation (29):
(3)
The pool of labeled cells in the blood, L, is refilled at a constant rate s until time , and labeled cells disappear from the pool of nave T cells in the blood with a disappearance rate d. In developing this semiempirical model, it is assumed that lymphoid tissue serves as an effective source of labeled cells that are distributed to the blood until equilibration is reached (time ), at which point the effective source ceases to affect changes in the concentration of labeled cells (29). Because BrdU+ chromosomes segregate independently into daughter cells, labeled cells that have divided will still be counted as BrdU+ cells, as long as the intensity of BrdU in each cell is higher than the threshold of flow cytometric detection. In human lymphocytes we estimate that the BrdU intensity decreases below the detection threshold after two to three divisions (data not shown). Thus, for highly proliferating cells equation 3 is approximately similar to a model where d consists only of the disappearance rate of labeled cells (29). For more slowly proliferating cells, such as nave T cells, a possible contribution of proliferation can be taken into account by replacing d with d – p, which is similar to an equation proposed by Debacq and colleagues (4). When d is small (<<1) the solution of equation 3 for t is given by , or the fraction of labeled cells increases approximately linearly over time with a rate s. Thus, the solution of equation 3 used in this analysis is given by
(4)
Moreover, if the peak of labeling is reached at similar times between recipients, the value of the fraction of BrdU-labeled cells at the peak will correlate with s and, thus, with proliferation rates (4).
For BrdU labeling, nave CD4 cells are defined as CD45RO–, as additional markers were not utilized in these analyses. For CD4 cells, this is a good approximation of true nave cells (6, 30).
RESULTS
Effect of HAART on viral load, T-cell counts, and T-cell proliferation. Sixteen of 20 group 1 patients had suppression of viral loads to <50 copies/ml at time point 1, and 18 of 20 had the same result at time point 2. Subjects had a mean CD4 count increase compared to a baseline of 144 cells/μl (range, 33 to 515 cells/μl) at time point 1 (6 months) and 286 cells/μl (range, 17 to 643 cells/μl) at time point 2 (18 months; P < 0.01) (Table 3).
Statistically significant increases in nave CD4+ T cells were observed at both time points 1 and 2 compared to baseline (from 87 to 127 and 164 cells/μl, respectively; P < 0.01), while statistically significant increases in nave CD8+ T cells were only observed at time point 2 (from 114 to 175 cells/μl; P < 0.05) (Table 3). At baseline, nave CD4+ T-cell numbers were lower than nave CD8 T-cell numbers (P < 0.05) but not at time point 1 or time point 2. Thus, the increase in nave T cells from baseline was higher for nave CD4+ T cells than nave CD8+ T cells at time point 1 (mean increase, 55 versus 24 cells/μl; P < 0.05) and at time point 2 (mean increase, 104 versus 68 cells/μl; P < 0.05). The additional changes between time points 1 and 2 in nave T-cell numbers did not differ between nave CD4+ and nave CD8+ T cells.
Significant inverse Spearman rank correlations were observed between baseline nave T-cell counts and the relative change (n-fold) in nave T cells for both CD4+ and CD8+ T cells at time point 1 ( = –0.53, P < 0.05 for CD4; = –0.48, P < 0.05 for CD8) and time point 2 ( = –0.82, P < 0.01 for CD4; = –0.64, P < 0.01 for CD8).
Significant declines in Ki67 expression were observed in both CD4+ and CD8+ nave and central memory T cells at time point 1 (P < 0.01) (Table 3). Before HAART, a median of 3.8% of nave CD4 cells and 5% of nave CD8 cells were Ki67+, with the difference between the means reaching borderline statistical significance (P = 0.05). At time point 1 these numbers decreased to 1.7 and 2.3%, respectively (P < 0.01). The additional observed declines at time point 2 to 1.4% and 1.3%, respectively, were not statistically significant (P > 0.05, between time points 1 and 2). Moreover, the relative change (n-fold) in percent Ki67+ nave CD4+ T cells between time point 0 to time point 1 inversely correlated with the corresponding relative change in CD4+ nave T-cell counts at the same time points (Spearman rank correlation, = –0.54, P = 0.017). No similar correlation was observed between time points 1 and 2 or between any time points for nave CD8+ T cells.
Greater relative change (n-fold) in TRECs/μl than nave T cells/μl in CD4+ and CD8+ T cells after HAART. A statistically significant increase in TRECs/μl from time point 0 to time point 1 (P < 0.01) was seen for both nave CD4+ T cells and nave CD8+ T cells (Fig. 1). A greater relative increase (n-fold) was seen in the number of TRECs/μl than in the nave T-cell counts in both CD4+ (median relative increases of 2.2- and 1.7-fold, respectively, P < 0.01) and CD8+ T-cell pools (1.4- and 1.2-fold, respectively, P < 0.01) at time point 1. This observation is equivalent to an increase in the fraction of TRECs, i.e., TRECs per million nave T cells, after initiation of HAART (Fig. 1), consistent with previously reported data (21, 32). The relative increase was observed to be lower for TRECs/μl than for nave T-cell counts from time point 1 to time point 2, with this difference being statistically significant for CD4+ T cells (P < 0.05) but not for CD8+ T cells. No significant correlations were observed between the relative change in percent Ki67+ nave CD4+ T cells at time point 1 to time point 0 and the relative change in TREC+ cells at the same time points (for CD4+, n = 14, = –0.03, and P = 0.47; for CD8+, n = 16, = –0.05, and P = 0.41). The lack of correlation between changes in the proliferation of nave T cells and the change in the fraction of TREC+ cells per nave T cell suggests that it is not the extent of reduction in proliferation per se that can explain the greater relative change in TRECs/μl versus nave T-cell counts, but the latter does probably depend on changes in both the proliferation and disappearance rate of the nave pool after the initiation of HAART.
Thymic indices did not significantly change during HAART and did not correlate with changes in TRECs. CT scoring was carried out by two radiologists blinded to clinical and laboratory data. The median baseline score was 1.5 (0.5 to 3.0). Median scores remained constant at time point 1 (1.5; range, 0.5 to 4) and at time point 2 (1.5; range, 0.5 to 4), with a median change for both time points of 0 (Table 3). When scores were clustered according to subjects with greater than and less than or equal to the median CD4 increase of 96 cells/μl at time point 1 and 257 cells/μl at time point 2, differences again were not significant. By volumetric analysis, the median score prior to initiating HAART was 1.58 cm3 (0.75 to 11.9). At time point 1 there was a median increase of 0.17 cm3 (–4.24 to 9.33); an additional period of approximately 12 months of HAART resulted in a median change of –0.04 cm3 (range, –5.27 to 6.5). None of these changes was statistically significant (Table 3). Baseline nave CD4+, but not nave CD8+, T-cell counts positively correlate with thymic volume (data not shown). An inverse statistically significant correlation was observed between the relative change (n-fold) in nave T cell counts and baseline CT scores for nave CD4+ T cells ( = –0.52, P < 0.05) but not for nave CD8+ T cells ( = –0.08, P > 0.05). Partial Spearman rank correlation analysis that included age or baseline nave T-cell counts did not qualitatively change the statistical significances of the above correlations. Similar patterns of correlations were observed when CT scores were replaced by thymic volumes. In addition, thymic volume changes did not correlate with changes in any TREC parameters.
Analysis of the model. A possible explanation for the greater relative change in TRECs/μl versus nave T-cell counts can be provided by a model of nave T-cell dynamics similar to the one described in equation 1 with the disappearance rate, d, a composite of two different factors: d = A + R, where A represents the death rate of nave cells and R is the rate of priming of nave cells into memory cells (7, 23). Assuming that the increase in the proliferation rate of nave T cells during chronic HIV-1 infection (p p + ) is mostly explained by the increase in the rate of priming of nave T cells into memory cells (R R + ) (Fig. 2A), nave T-cell counts will not substantially change, since the increased recruitment of nave T cells into memory cells is counterbalanced by the simultaneous increase in the proliferation rate of nave T cells (without losing their phenotypes). Conversely, if the main effect of HAART consists of reducing the rate of priming ( 0, with time), again, nave T-cell counts will be only marginally affected since the consequent decrease in recruitment (and consequent loss) of nave T cells is now counterbalanced by the simultaneous decrease in the proliferation rate. However, the reduction in proliferation will lead to an increase of the nave T-cell (and thus TREC+ T cell) life span or, equivalently, the decrease in disappearance rate d. This leads to an increase in TREC+ T cells/μl, and since the nave T-cell counts are only marginally affected, the fraction of TREC+ T cells per nave T cell is also predicted to increase.
With values of >0 this generalization predicts a greater decrease in the proliferation than the disappearance rate when a new quasi-steady state is reached following initiation of HAART (data not shown). The dynamics of nave CD8+ T cells observed in this study are consistent with this: a significant increase in the number of TREC+ T cells/μl is accompanied by a limited increase in nave T-cell counts after initiation of HAART (Table 3 and Fig. 1 and 2B). This model also predicts that the decrease in the proliferation rate after HAART is not expected to correlate with the increase in nave T-cell counts, since the latter will be only marginally affected regardless of the rate at which nave T-cell proliferation normalizes. Again, our data are consistent with this: there is no significant correlation between the increase of CD8+ nave T-cell counts and the decrease of percent Ki67+ CD8+ nave T cells between time points 0 and 1 or later.
For CD4+ nave T cells, the above model is inadequate to explain the observed dynamics. However, the transient increase in the fraction of TREC+ T cells per nave T cell and the following new steady state reached after 18 months of antiretroviral therapy can be explained by assuming that the death rate, A, is also enhanced (A A + ) during HIV-1 infection due to activation-induced cell death. The effect of HAART is again to normalize the proliferation rate and the disappearance rate of the CD4+ nave T cells. However, the latter now includes both the death rate and the priming rate, which demonstrate differential dynamics related to changes in and , respectively. A HAART-induced normalization of the death rate that is approximately 10-fold slower than the normalization of the proliferation rate adequately accounts for the changes in both TREC fractions and TREC+ T cells/μl that are observed for CD4+ nave T cells (Fig. 2C).
Previous reports have identified a positive correlation between baseline percentages of proliferating and apoptotic (terminal deoxynucleotidyltransferase-mediated UTP nick end labeling-positive cells) T cells (36). If the HIV-1-induced increases in proliferation and death rates of nave T cells ( and ) are similarly proportionally related, then the model also predicts an inverse correlation between the HAART-induced decrease in proliferation of nave T cells and the recovery of nave T-cell counts (data not shown) which we observed for CD4+ nave T cells. Thus, HIV-1-infected patients with higher proliferation and death rates of nave T cells during chronic infection would demonstrate a faster recovery of nave T-cell counts after HAART, with most of the recovery explained by normalization of the death rate rather than by normalization of the proliferation rate, which, as in the earlier model, is still counterbalanced by the simultaneous reduction in the rate of priming.
Kinetics of BrdU-labeled nave T cells prior to and after HAART. To further evaluate the effects of HAART on the proliferation and disappearance rates of nave T cells, we analyzed the kinetics of BrdU-labeled nave T cells prior to and after initiation of HAART in six HIV-1-infected patients. Figure 3 shows the theoretical curves obtained by best fitting equation 3 to the fraction of BrdU-labeled nave CD4+ T cells for each patient, before and after initiation of therapy. Because the percentages of BrdU-labeled cells within the nave T-cell population are relatively small, and thus closer to the detection threshold, this analysis resulted in a weak convergence of the nonlinear best-fitting procedure using equation 3. Thus, equation 4 was used to estimate the slope s from the time of infusion to the peak of labeling, and a single exponential decay function has been used to estimate d from the time of peak labeling. As shown in Table 4, this modeling predicts a decay in s greater than a decay in d for nave CD4+ T cells after initiation of HAART (median relative decrease [n-fold] in s of 2.3 and median relative decrease in d of 1.5; P = 0.03). Since changes in s imply changes in mean proliferation rates, a greater decay in s than d implies a more dramatic decrease in proliferation than disappearance rates of nave T cells after HAART, whether d represents the apoptotic death rate plus the rate of priming into memory T cells or the difference between the disappearance rate in the blood, d, and the proliferation rate in the blood, p (data not shown).
DISCUSSION
In this study we have analyzed the relationship between T-cell turnover, thymic function, and immune activation in HIV-1-infected patients, focusing on nave CD4+ and nave CD8+ T cells, to better understand the contribution of these various parameters to the immunologic changes seen during HIV infection and therapy. We specifically targeted a broad range of baseline viral load and CD4+ T-cell counts in order to highlight the dynamics across the spectrum of HIV infection. At baseline, naive T-cell numbers were lower in the CD4 pool compared to the CD8 pool, but no difference between the two groups was observed in the percentages of proliferating nave T cells or in the number of TRECs. Even though a dramatic decrease in proliferation is observed for both nave CD4+ and CD8+ T cells after the first 6 months of therapy, we found that nave CD4+ T-cell numbers significantly increased after initiation of HAART, but nave CD8+ T-cell numbers were only marginally affected. Baseline nave T-cell counts inversely correlated with the relative changes in nave T cells in both subpopulations of cells, and TRECs/μl increased in both subsets of cells. We also observed an increase in the fraction of TREC+ T cells per nave T cell, which is equivalent to a faster growth of TREC+ T cells versus total nave T cells, for both CD4 and CD8 cells following initiation of HAART. We did not observe significant changes of thymic volumes during time or statistically significant correlations between changes in thymic volumes and changes in nave T-cell numbers or TRECs/μl.
To provide a unified description of the concomitant changes in nave T cells (TRECs/μl and the percentage of proliferating nave T cells during HAART) that can explain the observed differences in the dynamics of these variables in the CD4 and CD8 subpopulations of T lymphocytes, we developed a mathematical model based on a generalization of a model originally described by Hazenberg et al. (21). We have generalized the disappearance rate of nave T cells as the sum of the rates of nave cells priming into memory cells and nave cell death and assumed that the increase in proliferation of nave T cells during chronic infection is primarily explained by the increase in the rate of priming of nave T cells into memory cells. This simple theoretical framework is sufficient to predict the simultaneous increases in the fraction of TREC+ T cells per nave T cell and the number of TREC+ T cells in the periphery after initiation of HAART. The model explains the dynamics of nave CD8+, but not CD4+, T cells after institution of HAART. To explain the dynamics of nave CD4+ T cells, we postulate that there is an increase in the apoptotic death rate of nave T cells during HIV-1 infection related to immune activation or to the increase in proliferation rate and that there is delayed normalization of the apoptotic death rate compared to the proliferation rate, as has been reported for total CD4+ T cells in lymph node samples of HIV-1-infected patients before and after initiation of HAART (48). This model does not require (or exclude) changes in thymic output or redistribution of T lymphocytes from the lymphoid tissue as additional mechanisms contributing to nave T-cell recovery.
The normalization of the death rates can also account for the inverse correlation between baseline nave T-cell counts and the relative change in nave T-cell counts after initiation of HAART. However, the increase in nave T cells appears to be lower for nave CD8 T cells than for the nave CD4+ T cells, which suggests independent mechanisms of peripheral normalization for the different populations. This dichotomy, observed in the dynamics of nave CD4+ T cells compared to nave CD8+ T cells after HAART, as well as the presence of an inverse correlation between baseline thymic scores and the relative change (n-fold) in nave T-cell numbers for CD4+ but not CD8+ T cells is difficult to explain solely as a result of changes in thymic output rates or trafficking effects, since these should not have differential effects on the two populations of nave T cells. Moreover, based on this model, the observed inverse correlation between baseline counts and the relative change in nave T-cell counts after HAART for both nave CD4+ and nave CD8+ T cells suggests that increases in the death rate affect both populations. However, the greater baseline depletion of nave CD4+ T cells compared to nave CD8+ T cells, together with the concomitant increase in TRECs per microliter in both compartments following therapy, suggests, as highlighted by the model, that similar mechanisms drive both nave CD4+ and nave CD8+ T cells to be primed into memory cells, but for unknown reasons the increase in the death rate is more pronounced in the CD4+ than in the CD8+ nave T-cell populations. Interestingly, Li et al. have recently shown that, at least in the settings of acute simian immunodeficiency virus infection, higher levels of Fas- and Fas ligand-mediated apoptosis are observed within CD4+ but not CD8+ T lymphocytes in the lamina propria, which may result from massive exposure of CD4+ T cells to virion gp120 (33). Our data suggest that during chronic infection, the differential ability to tolerate similar increases in proliferation, , is an intrinsic property of each subpopulation of nave T cells. An alternative scenario in which is different between the two subpopulations would require that nave CD4 T cells have a higher proliferation rate than nave CD8 T cells to account for the relative loss in nave CD4 T cells. However, this was not observed in our data when we looked at the baseline fractions of proliferating nave CD4 and CD8 T cells.
It is important to note that this model represents an idealized situation and that deviations from this model, resulting, for instance, from the presence of nonlinear contributions of trafficking of lymphocytes or of replenishment of the peripheral pool by thymic output, might result in a situation that is far from the quasi-steady-state condition assumed in equation 1. In the latter circumstances, TREC content after initiation of HAART can potentially be affected in an unpredictable manner.
The hypothesis of a more pronounced decrease in the proliferation rate than the disappearance rate of nave T cells after HAART is supported by the kinetics of BrdU-labeled nave T cells studied longitudinally (pre- and post-HAART) that we observed in a smaller group of HIV-1-infected patients.
In principle, changes in the fraction of nave T cells carrying TRECs upon exiting the thymus (f in equation 1) might also explain a greater relative change in TREC+ T cells/μl than nave T-cell counts after initiation of HAART (32). Among the four parameters discussed in this analysis (d, p, , and f), f is the least investigated. Dion and colleagues (8) have recently shown an increase of the ratio -TRECs/-TRECs after initiation of HAART, which suggests that the newly produced nave T cells undergo more intrathymic divisions before entering the peripheral pool. Since -TRECs are produced before -TRECs, this would lead to a decrease, not an increase, in f after initiation of HAART, thus excluding changes in f as a major factor affecting the dynamics of the fraction of TREC+ T cells after initiation of HAART.
This analysis provides evidence that changes in peripheral proliferation and disappearance rates of nave T cells, rather than changes in thymic output, explain the observed dynamics of TRECs and the fraction of proliferating T cells during HAART. But what is the mechanism that drives nave T cells to proliferate faster during chronic infection The first pathogenic effect of HIV-1 infection might consist of an increase in the rate of priming of nave T cells due to a generalized state of chronic immune activation. In this scenario, the HIV-1-induced increase in the proliferation rate could serve as a compensatory (homeostatic) mechanism aimed at maintaining the nave T-cell count constant, in response to the loss of nave T cells that have been primed into memory cells. Alternatively, the p(t) expression of the proliferation rate of our model could also be modeled as a function of the number of cells, T(t), for instance, following a density-dependent law (10). Under this scenario, changes in thymic output and consequent (homeostatic) changes in peripheral proliferation of nave T cells could also explain the observed TREC dynamics, as suggested by Dutilh et al. to explain the changes of TREC content during aging (10). However, given the relatively short time frame (a few weeks) in which significant changes of proliferation and disappearance rates are seen (36) and the negligible contribution by the thymus seen in thymectomy studies (as recently described for nonhuman primates by Arron et al. [1]), as well as our thymic CT data, we feel that peripheral mechanisms (including homeostatic mechanisms) leading to changes in proliferation and disappearance rates are the major factors affecting the dynamics of TRECs. The observed dichotomy in the dynamics of nave CD4+ and nave CD8+ T cells is difficult to explain by invoking changes in thymic output as the sole mechanism responsible for changes in peripheral (homeostatic) proliferation rates of nave T cells.
If a homeostatic mechanism governs the proliferation of nave T cells, we would expect, as has been previously suggested for the entire population of T cells (24), an inverse correlation between relative change (n-fold) in nave T-cell counts after HAART and the relative changes in the percentages of proliferating nave T cells. In our study we do observe such an inverse correlation for nave CD4+ but not nave CD8+ T cells. Alternatively, the first effect of activation might consist of inducing nave T cells to proliferate faster without losing their phenotype (45-47), bringing these cells closer to the priming activation threshold which leads to an increased disappearance rate of nave T cells. In this scenario, the increase in the proliferation rate is primarily explained by the increase in the rate of priming during chronic infection. The simultaneous decrease in both the proliferation and priming rates after initiation of HAART should generate a lack of correlation between the recovery of nave T-cell counts (only marginally affected) and the decrease in the percentage of proliferating nave T cells. This latter paradigm would explain the observed lack of such correlation for nave CD8+ T cells. The additional assumption that higher levels of proliferation are associated with higher levels of apoptosis is required to explain the presence of this inverse correlation for CD4+ nave T cells. Both scenarios show consistency with a differential change in the proliferation and disappearance rates of nave T-cells. Based on current available data, it is difficult to make conclusive arguments in support of either hypothesis. Thus, further investigations are required to clarify the mechanisms responsible for the increased proliferation of nave T cells induced by HIV-1.
ACKNOWLEDGMENTS
We thank the patients for their willingness to participate in these studies, the physicians, nurses, and other members of the National Institute of Allergy and Infectious Diseases-Clinical Center HIV/AIDS Program for their support of these studies, and Betty Wise for performing the quantitative thymic measurements.
This research was supported in part by the Intramural Research Program of the National Institutes of Health, National Institute of Allergy and Infectious Diseases, and Warren Grant Magnuson Clinical Center.
REFERENCES
Arron, S. T., R. M. Ribeiro, A. Gettie, R. Bohm, J. Blanchard, J. Yu, A. S. Perelson, D. D. Ho, and L. Zhang. 2005. Impact of thymectomy on the peripheral T cell pool in rhesus macaques before and after infection with simian immunodeficiency virus. Eur. J. Immunol. 35:46-55.
Davey, R. T., Jr., N. Bhat, C. Yoder, T. W. Chun, J. A. Metcalf, R. Dewar, V. Natarajan, R. A. Lempicki, J. W. Adelsberger, K. D. Miller, J. A. Kovacs, M. A. Polis, R. E. Walker, J. Falloon, H. Masur, D. Gee, M. Baseler, D. S. Dimitrov, A. S. Fauci, and H. C. Lane. 1999. HIV-1 and T cell dynamics after interruption of highly active antiretroviral therapy (HAART) in patients with a history of sustained viral suppression. Proc. Natl. Acad. Sci. USA 96:15109-15114.
Davey, R. T., Jr., R. L. Murphy, F. M. Graziano, S. L. Boswell, A. T. Pavia, M. Cancio, J. P. Nadler, D. G. Chaitt, R. L. Dewar, D. K. Sahner, A. M. Duliege, W. B. Capra, W. P. Leong, M. A. Giedlin, H. C. Lane, and J. O. Kahn. 2000. Immunologic and virologic effects of subcutaneous interleukin 2 in combination with antiretroviral therapy: A randomized controlled trial. JAMA 284:183-189.
Debacq, C., B. Asquith, P. Kerkhofs, D. Portetelle, A. Burny, R. Kettmann, and L. Willems. 2002. Increased cell proliferation, but not reduced cell death, induces lymphocytosis in bovine leukemia virus-infected sheep. Proc. Natl. Acad. Sci. USA 99:10048-10053.
De Boer, R. J., H. Mohri, D. D. Ho, and A. S. Perelson. 2003. Turnover rates of B cells, T cells, and NK cells in simian immunodeficiency virus-infected and uninfected rhesus macaques. J. Immunol. 170:2479-2487.
De Rosa, S. C., L. A. Herzenberg, and M. Roederer. 2001. 11-color, 13-parameter flow cytometry: identification of human naive T cells by phenotype, function, and T-cell receptor diversity. Nat. Med. 7:245-248.
Di Mascio, M., I. Sereti, L. Matthews, V. Natarajan, C. Yoder, E. C. Jones, C. Chow, I. A. Sidorov, D. S. Dimitrov, J. A. Metcalf, M. A. Polis, and J. A. Kovacs. 2004. Presented at the XV International AIDS Conference, Bangkok, Thailand, 11 to 16 July 2004.
Dion, M. L., J. F. Poulin, R. Bordi, M. Sylvestre, R. Corsini, N. Kettaf, A. Dalloul, M. R. Boulassel, P. Debre, J. P. Routy, Z. Grossman, R. P. Sekaly, and R. Cheynier. 2004. HIV infection rapidly induces and maintains a substantial suppression of thymocyte proliferation. Immunity 21:757-768.
Douek, D. C., R. D. McFarland, P. H. Keiser, E. A. Gage, J. M. Massey, B. F. Haynes, M. A. Polis, A. T. Haase, M. B. Feinberg, J. L. Sullivan, B. D. Jamieson, J. A. Zack, L. J. Picker, and R. A. Koup. 1998. Changes in thymic function with age and during the treatment of HIV infection. Nature 396:690-695.
Dutilh, B. E., and R. J. de Boer. 2003. Decline in excision circles requires homeostatic renewal or homeostatic death of naive T cells. J. Theor. Biol. 224:351-358.
Fleury, S., R. J. de Boer, G. P. Rizzardi, K. C. Wolthers, S. A. Otto, C. C. Welbon, C. Graziosi, C. Knabenhans, H. Soudeyns, P. A. Bart, S. Gallant, J. M. Corpataux, M. Gillet, P. Meylan, P. Schnyder, J. Y. Meuwly, W. Spreen, M. P. Glauser, F. Miedema, and G. Pantaleo. 1998. Limited CD4+ T-cell renewal in early HIV-1 infection: effect of highly active antiretroviral therapy. Nat. Med. 4:794-801.
Fleury, S., G. P. Rizzardi, A. Chapuis, G. Tambussi, C. Knabenhans, E. Simeoni, J. Y. Meuwly, J. M. Corpataux, A. Lazzarin, F. Miedema, and G. Pantaleo. 2000. Long-term kinetics of T cell production in HIV-infected subjects treated with highly active antiretroviral therapy. Proc. Natl. Acad. Sci. USA 97:5393-5398.
Franco, J. M., A. Rubio, M. Martinez-Moya, M. Leal, E. Merchante, A. Sanchez-Quijano, and E. Lissen. 2002. T-cell repopulation and thymic volume in HIV-1-infected adult patients after highly active antiretroviral therapy. Blood 99:3702-3706.
Gaulton, G. N., J. V. Scobie, and M. Rosenzweig. 1997. HIV-1 and the thymus. AIDS 11:403-414.
Gerdes, J., H. Lemke, H. Baisch, H. H. Wacker, U. Schwab, and H. Stein. 1984. Cell cycle analysis of a cell proliferation-associated human nuclear antigen defined by the monoclonal antibody Ki-67. J. Immunol. 133:1710-1715.
Gill, P. E., W. Murray, and M. H. Wright. 1981. Practical optimization, p. 136-137. Academic Press, London, United Kingdom.
Giorgi, J. V., L. E. Hultin, J. A. McKeating, T. D. Johnson, B. Owens, L. P. Jacobson, R. Shih, J. Lewis, D. J. Wiley, J. P. Phair, S. M. Wolinsky, and R. Detels. 1999. Shorter survival in advanced human immunodeficiency virus type 1 infection is more closely associated with T lymphocyte activation than with plasma virus burden or virus chemokine coreceptor usage. J. Infect. Dis. 179:859-870.
Hazenberg, M., S. Galkina, G. Chkhenkeli, C. Stoddart, and M. McCune. 2004. Presented at the 11th Conference on Retroviruses and Opportunistic Infections, San Francisco, Calif.
Hazenberg, M. D., D. R. Clark, and F. Miedema. 1999. Tilted balance of cell renewal in HIV-1 infection. AIDS Rev. 1:67-73.
Hazenberg, M. D., D. Hamann, H. Schuitemaker, and F. Miedema. 2000. T cell depletion in HIV-1 infection: how CD4+ T cells go out of stock. Nat. Immunol. 1:285-289.
Hazenberg, M. D., S. A. Otto, J. W. Cohen Stuart, M. C. Verschuren, J. C. Borleffs, C. A. Boucher, R. A. Coutinho, J. M. Lange, T. F. Rinke de Wit, A. Tsegaye, J. J. van Dongen, D. Hamann, R. J. de Boer, and F. Miedema. 2000. Increased cell division but not thymic dysfunction rapidly affects the T-cell receptor excision circle content of the naive T cell population in HIV-1 infection. Nat. Med. 6:1036-1042.
Hazenberg, M. D., S. A. Otto, B. H. van Benthem, M. T. Roos, R. A. Coutinho, J. M. Lange, D. Hamann, M. Prins, and F. Miedema. 2003. Persistent immune activation in HIV-1 infection is associated with progression to AIDS. AIDS 17:1881-1888.
Hazenberg, M. D., S. A. Otto, A. M. van Rossum, H. J. Scherpbier, R. de Groot, T. W. Kuijpers, J. M. Lange, D. Hamann, R. J. de Boer, J. A. Borghans, and F. Miedema. 2004. Establishment of the CD4+ T-cell pool in healthy children and untreated children infected with HIV-1. Blood 104:3513-3519.
Hazenberg, M. D., J. W. Stuart, S. A. Otto, J. C. Borleffs, C. A. Boucher, R. J. de Boer, F. Miedema, and D. Hamann. 2000. T-cell division in human immunodeficiency virus (HIV)-1 infection is mainly due to immune activation: a longitudinal analysis in patients before and during highly active antiretroviral therapy (HAART). Blood 95:249-255.
Hellerstein, M., M. B. Hanley, D. Cesar, S. Siler, C. Papageorgopoulos, E. Wieder, D. Schmidt, R. Hoh, R. Neese, D. Macallan, S. Deeks, and J. M. McCune. 1999. Directly measured kinetics of circulating T lymphocytes in normal and HIV-1-infected humans. Nat. Med. 5:83-89.
Hellerstein, M. K., R. A. Hoh, M. B. Hanley, D. Cesar, D. Lee, R. A. Neese, and J. M. McCune. 2003. Subpopulations of long-lived and short-lived T cells in advanced HIV-1 infection. J. Clin. Investig. 112:956-966.
Highbarger, H. C., W. G. Alvord, M. K. Jiang, A. S. Shah, J. A. Metcalf, H. C. Lane, and R. L. Dewar. 1999. Comparison of the Quantiplex version 3.0 assay and a sensitized Amplicor monitor assay for measurement of human immunodeficiency virus type 1 RNA levels in plasma samples. J. Clin. Microbiol. 37:3612-3614.
Kaufmann, G. R., J. Zaunders, J. Murray, A. D. Kelleher, S. R. Lewin, A. Solomon, D. Smith, and D. A. Cooper. 2001. Relative significance of different pathways of immune reconstitution in HIV type 1 infection as estimated by mathematical modeling. AIDS Res. Hum. Retrovir. 17:147-159.
Kovacs, J. A., R. A. Lempicki, I. A. Sidorov, J. W. Adelsberger, B. Herpin, J. A. Metcalf, I. Sereti, M. A. Polis, R. T. Davey, J. Tavel, J. Falloon, R. Stevens, L. Lambert, R. Dewar, D. J. Schwartzentruber, M. R. Anver, M. W. Baseler, H. Masur, D. S. Dimitrov, and H. C. Lane. 2001. Identification of dynamically distinct subpopulations of T lymphocytes that are differentially affected by HIV. J. Exp. Med. 194:1731-1741.
Kovacs, J. A., S. Vogel, J. A. Metcalf, M. Baseler, R. Stevens, J. Adelsberger, R. Lempicki, R. L. Hengel, I. Sereti, L. Lambert, R. L. Dewar, R. T. Davey, Jr., R. E. Walker, J. Falloon, M. A. Polis, H. Masur, and H. C. Lane. 2001. Interleukin-2 induced immune effects in human immunodeficiency virus-infected patients receiving intermittent interleukin-2 immunotherapy. Eur. J. Immunol. 31:1351-1360.
Lempicki, R. A., J. A. Kovacs, M. W. Baseler, J. W. Adelsberger, R. L. Dewar, V. Natarajan, M. C. Bosche, J. A. Metcalf, R. A. Stevens, L. A. Lambert, W. G. Alvord, M. A. Polis, R. T. Davey, D. S. Dimitrov, and H. C. Lane. 2000. Impact of HIV-1 infection and highly active antiretroviral therapy on the kinetics of CD4+ and CD8+ T cell turnover in HIV-infected patients. Proc. Natl. Acad. Sci. USA 97:13778-13783.
Lewin, S. R., R. M. Ribeiro, G. R. Kaufmann, D. Smith, J. Zaunders, M. Law, A. Solomon, P. U. Cameron, D. Cooper, and A. S. Perelson. 2002. Dynamics of T cells and TCR excision circles differ after treatment of acute and chronic HIV infection. J. Immunol. 169:4657-4666.
Li, Q., L. Duan, J. D. Estes, Z. M. Ma, T. Rourke, Y. Wang, C. Reilly, J. Carlis, C. J. Miller, and A. T. Haase. 2005. Peak SIV replication in resting memory CD4+ T cells depletes gut lamina propria CD4+ T cells. Nature 434:1148-1152.
Liu, Z., W. G. Cumberland, L. E. Hultin, H. E. Prince, R. Detels, and J. V. Giorgi. 1997. Elevated CD38 antigen expression on CD8+ T cells is a stronger marker for the risk of chronic HIV disease progression to AIDS and death in the Multicenter AIDS Cohort Study than CD4+ cell count, soluble immune activation markers, or combinations of HLA-DR and CD38 expression. J. Acquir. Immune Defic. Syndr. Hum. Retrovirol. 16:83-92.
Lu, A. C., E. C. Jones, C. Chow, K. D. Miller, B. Herpin, D. Rock-Kress, J. A. Metcalf, H. C. Lane, and J. A. Kovacs. 2003. Increases in CD4+ T lymphocytes occur without increases in thymic size in HIV-infected subjects receiving interleukin-2 therapy. J. Acquir. Immune Defic. Syndr. 34:299-303.
Mohri, H., A. S. Perelson, K. Tung, R. Ribeiro, B. Ramratnam, M. Markowitz, R. Kost, A. Hurley, L. Weinberger, D. Cesar, M. Hellerstein, and D. D. Ho. 2001. Increased turnover of T lymphocytes in HIV-1 infection and its reduction by antiretroviral therapy. J. Exp. Med. 194:1277-1287.
Murali-Krishna, K., and R. Ahmed. 2000. Cutting edge: naive T cells masquerading as memory cells. J. Immunol. 165:1733-1737.
Natarajan, V., R. A. Lempicki, I. Sereti, Y. Badralmaa, J. W. Adelsberger, J. A. Metcalf, D. A. Prieto, R. Stevens, M. W. Baseler, J. A. Kovacs, and H. C. Lane. 2002. Increased peripheral expansion of naive CD4+ T cells in vivo after IL-2 treatment of patients with HIV infection. Proc. Natl. Acad. Sci. USA 99:10712-10717.
Roederer, M., J. G. Dubs, M. T. Anderson, P. A. Raju, L. A. Herzenberg, and L. A. Herzenberg. 1995. CD8 naive T cell counts decrease progressively in HIV-infected adults. J. Clin. Investig. 95:2061-2066.
Sachsenberg, N., A. S. Perelson, S. Yerly, G. A. Schockmel, D. Leduc, B. Hirschel, and L. Perrin. 1998. Turnover of CD4+ and CD8+ T lymphocytes in HIV-1 infection as measured by Ki-67 antigen. J. Exp. Med. 187:1295-1303.
Schacker, T., S. Little, E. Connick, K. Gebhard, Z. Q. Zhang, J. Krieger, J. Pryor, D. Havlir, J. K. Wong, R. T. Schooley, D. Richman, L. Corey, and A. T. Haase. 2001. Productive infection of T cells in lymphoid tissues during primary and early human immunodeficiency virus infection. J. Infect. Dis. 183:555-562.
Sereti, I., K. B. Anthony, H. Martinez-Wilson, R. Lempicki, J. Adelsberger, J. A. Metcalf, C. W. Hallahan, D. Follmann, R. T. Davey, J. A. Kovacs, and H. C. Lane. 2004. IL-2-induced CD4+ T-cell expansion in HIV-infected patients is associated with long-term decreases in T-cell proliferation. Blood 104:775-780.
Sousa, A. E., J. Carneiro, M. Meier-Schellersheim, Z. Grossman, and R. M. Victorino. 2002. CD4 T cell depletion is linked directly to immune activation in the pathogenesis of HIV-1 and HIV-2 but only indirectly to the viral load. J. Immunol. 169:3400-3406.
Sprent, P. 1993. Applied nonparametric statistical methods, 2nd ed. Chapman and Hall, New York, N.Y.
Tough, D. F., and J. Sprent. 1994. Turnover of naive- and memory-phenotype T cells. J. Exp. Med. 179:1127-1135.
Unutmaz, D., F. Baldoni, and S. Abrignani. 1995. Human naive T cells activated by cytokines differentiate into a split phenotype with functional features intermediate between naive and memory T cells. Int. Immunol. 7:1417-1424.
Unutmaz, D., P. Pileri, and S. Abrignani. 1994. Antigen-independent activation of naive and memory resting T cells by a cytokine combination. J. Exp. Med. 180:1159-1164.
Zhang, Z. Q., D. W. Notermans, G. Sedgewick, W. Cavert, S. Wietgrefe, M. Zupancic, K. Gebhard, K. Henry, L. Boies, Z. Chen, M. Jenkins, R. Mills, H. McDade, C. Goodwin, C. M. Schuwirth, S. A. Danner, and A. T. Haase. 1998. Kinetics of CD4+ T cell repopulation of lymphoid tissues after treatment of HIV-1 infection. Proc. Natl. Acad. Sci. USA 95:1154-1159.(Michele Di Mascio, Irini )
SAIC, Frederick, Maryland
Clinical Center, National Institutes of Health, Bethesda, Maryland
National Cancer Institute-Frederick, National Institutes of Health, Frederick, Maryland
ABSTRACT
Both nave CD4+ and nave CD8+ T cells are depleted in individuals with human immunodeficiency virus type 1 (HIV-1) infection by unknown mechanisms. Analysis of their dynamics prior to and after highly active antiretroviral therapy (HAART) could reveal possible mechanisms of depletion. Twenty patients were evaluated with immunophenotyping, intracellular Ki67 staining, T-cell receptor excision circle (TREC) quantitation in sorted CD4 and CD8 cells, and thymic computed tomography scans prior to and 6 and 18 months after initiation of HAART. Nave T-cell proliferation decreased significantly during the first 6 months of therapy (P < 0.01) followed by a slower decline. Thymic indices did not change significantly over time. At baseline, nave CD4+ T-cell numbers were lower than naive CD8+ T-cell numbers; after HAART, a greater increase in nave CD4+ T cells than nave CD8+ T cells was observed. A greater relative change (n-fold) in the number of TREC+ T cells/μl than in nave T-cell counts was observed at 6 months for both CD4+ (median relative change [n-fold] of 2.2 and 1.7, respectively; P < 0.01) and CD8+ T cell pools (1.4 and 1.2; P < 0.01). A more pronounced decrease in the proliferation than the disappearance rate of nave T cells after HAART was observed in a second group of six HIV-1-infected patients studied by in vivo pulse labeling with bromodeoxyuridine. These observations are consistent with a mathematical model where the HIV-1-induced increase in proliferation of nave T cells is mostly explained by a faster recruitment into memory cells.
INTRODUCTION
The mechanisms leading to CD4+ T-cell depletion during human immunodeficiency virus (HIV) infection and CD4+ T-cell restoration during therapy with highly active antiretroviral therapy (HAART) remain undefined. Direct viral cytopathogenicity or redistribution of lymphocytes between lymphoid tissues and the circulation can account only in part for these changes (28, 32, 41). Characterizing lymphocyte turnover in HIV-infected patients can potentially lead to a better understanding of these mechanisms.
A number of studies utilizing indirect methods (2, 21, 31, 40, 48) and, more recently, direct methods using deuterated glucose (25, 36) or 5-bromo-2'-deoxyuridine (BrdU) (29) to measure CD4 cell dynamics have shown that CD4 cell turnover is increased during chronic HIV type 1 (HIV-1) infection. Moreover, while preliminary cross-sectional studies described an increase in CD4 cell turnover in patients following initiation of HAART, suggesting a defect in CD4 cell production secondary to HIV infection (12, 25), longitudinal studies have clearly demonstrated a rapid and persistent decrease in CD4 cell proliferation following initiation of HAART, suggesting that the increase in CD4 cell turnover itself may be an important pathogenic mechanism of CD4 depletion (21, 29, 36). While this increased turnover was initially postulated to represent a homeostatic response to CD4 depletion (11, 12), such a hypothesis is inconsistent with the rapid reduction in proliferation of CD4+ as well as CD8+ T cells after viral suppression with HAART, prior to normalization of CD4 cell numbers (2, 29, 31, 36). These observations led to alternative hypotheses proposing that either HIV-directed or nonspecific immune activation drives increased turnover. Moreover, based in part on studies demonstrating that levels of immune activation in T cells, especially CD8 cells, are independent predictors of CD4 depletion and disease progression, immune activation is currently felt by many investigators to play a direct role in HIV-associated CD4 depletion (17, 34, 43).
An increase in turnover has been demonstrated in nave T cells (24, 26) as well as memory T cells during pathogenic lentiviral infection. Whereas memory CD4+ but not memory CD8+ T cells decrease in number during chronic HIV infection, both nave CD4+ as well as nave CD8+ T cells are depleted during such infection (5, 22, 39). This observation led to the conclusion that nave T-cell depletion is one of the hallmarks of HIV infection. While infection of nave T cells has been documented, this appears to be a relatively rare event that cannot quantitatively explain the loss of nave CD4+ T cells. The observation that both nave CD4+ and nave CD8+ T cells decrease during HIV infection led to the hypothesis that persistent hyperactivation of the immune system leads to erosion of nave T cells by their increased recruitment into memory cells (20), probably through antigen- and nonantigen-specific stimulation, as has been shown in animal models (18, 37).
Because the thymus is the source of new T cells, examining thymic function during HIV infection and therapy is critical to studies of T-cell dynamics (14, 19). Due to the difficulties in directly studying thymic function, quantitation of T-cell receptor excision circles (TRECs) has been utilized as a surrogate of thymic function. Early studies measuring the number of TRECs per nave T cell suggested that HIV infection leads to a decrease in thymic function and that improved thymic function contributes to immune reconstitution following HAART (9). More recently, studies using mathematical modeling showed that the observed changes in TREC content cannot be explained solely by changes in thymic function (21) or by redistribution of T lymphocytes from lymphoid tissues to the blood (32); the observed TREC dynamics were more consistent with changes in peripheral proliferation and disappearance rates of T-lymphocytes.
In the current study we undertook a detailed examination of the relationship between T-cell turnover, thymic function, and immune activation in HIV-1-infected patients to better understand the contribution of these various parameters to the immunologic changes seen during HIV infection and therapy.
MATERIALS AND METHODS
Patient characteristics. Twenty-three HIV-1-infected patients who were initiating HAART (n = 20), starting HAART after a 6- and 24-month treatment interruption (n = 2), or switching from an ineffective to an effective HAART regimen (n = 1) were studied. Twenty previously untreated (or with limited nucleoside analogue exposure) HIV-1-infected patients (group 1) enrolled between 1997 and 2000 in a study to identify viral reservoirs in patients receiving antiretroviral therapy. These patients initiated therapy with indinavir, nevirapine, lamivudine, and zidovudine (or stavudine [one patient]). Drug changes were permitted for intolerance; for 14 patients, zidovudine was replaced by stavudine, and for 6 patients, indinavir was replaced by nelfinavir. After at least 1 year on study, additional changes were permitted at the patient's and referring physician's discretion; seven patients made such changes, usually for simplification of the regimen. All patients ultimately demonstrated a virologic response, with declines in plasma viral loads to levels below branched DNA assay limits (<50 copies/ml, 18 patients; <500 copies/ml, 2 patients) and all but 2 patients had an increase in CD4 cell numbers of at least 100 cells/mm3. Group 1 patients had immunophenotyping analysis and evaluation of thymic function by TREC analysis and thymic computed tomography (CT) scans at the following time points: 0 (pre-HAART), 1 (5 to 8 months, median of 6 months after start of therapy), and 2 (13 to 42 months, median of 18 months after start of therapy). Nave cell labeling kinetics were analyzed for six HIV-1-infected patients (group 2) who had undergone in vivo pulse labeling with BrdU prior to and 3 to 6 months after the initiation of an effective HAART regimen. Three of these patients are also included in group 1; the remaining three received stavudine, lamivudine, indinavir, and nevirapine (one patient); stavudine, lamivudine, and efavirenz (1 patient); and abacavir, amprenavir, nelfinavir, and efavirenz (one patient). Results for CD4 and CD8 cell labeling, but not nave cell labeling, have been previously reported for this group (29). Baseline characteristics for group 1 and pre- and post-HAART parameters for group 2 are shown in Tables 1 and 2.
Viral load. HIV-1 RNA levels were determined using a modification of the Roche Amplicor HIV Monitor assay kit (Indianapolis, IN) (27) or a branched DNA assay (Bayer Diagnostics, Tarrytown, NY) (3). Each assay had a lower limit of detection of 50 copies/ml; for two patients enrolled early in the study, the limit of detection was 500 copies/ml.
Immunophenotyping and intracellular staining for Ki67. Immunophenotypic analysis of cryopreserved peripheral blood mononuclear cells was performed using four-color immunofluorescence as previously described (42). Nave cells were defined as CD45RO– CD27+, central memory were defined as CD45RO+ CD27+, and effector memory were defined as CD45RO+ CD27– for CD4 T cells and as CD27– (CD45RO+ or CD45RO–) for CD8 T cells. Cells were stained intracellularly with Ki67-phycoerythrin (clone B56) or isotype (mouse immunoglobulin G1-phycoerythrin; clone MOPC-21) from BD/Pharmingen. T-cell proliferation was defined as the percentage of cells expressing Ki67 (15).
TREC determination. Signal joint TRECs (Sj-TRECs) in purified cell subsets were quantitated by real-time PCR by the cell lysis method as described previously (38). The consistency of the DNA content of the cell lysate was checked by real-time PCR using a ribosomal protein gene and a TaqMan gene expression assay kit from Applied Biosystems, Inc. (Foster City, CA). Because no more than one Sj-TREC can be produced per cell, the number of TRECs per unit volume of blood also represents the number of TREC+ T cells in the same unit volume.
Thymic CT scans. CT scans of the thymus were obtained prior to and at a median of 6 and 18 months after starting HAART. Scans were graded as previously described on a 0 (no thymic tissue) to 5 (thymic mass) scale by two independent radiologists blinded to clinical and laboratory results (35). In addition, computer-based density and volume analysis of the thymus was performed by transferring CT data to a GE Advantage Windows workstation (versions 2.1 and 4.0; GE Medical Systems, Advanced Windows Workstation Training Program, Milwaukee, WI). Contours of the anterior mediastinum were outlined by a radiologist-trained technician and corroborated by a radiologist (13).
BrdU infusion and flow cytometry. The fractions of BrdU+ CD4+ CD45RO– and BrdU+ CD8+ CD45RO– T cells were analyzed by flow cytometry as previously described (29).
Statistics. Changes in median values for each variable were tested for significance by the permutation test with paired samples computed with an exact method (44). Tests were performed using StatXact software. Association between variables was determined by the Spearman rank correlation test. Adjustment of P values for multiple testing was done by the Bonferroni method. Occasional data points were missing for group 1; however, all paired analyses were tested with n values that were 15.
The relative change (n-fold) of a variable H between time (t) points (e.g., t0 and t1) is defined as the ratio between the value of H at t1 versus t0: H(t1)/H(t0).
Modeling. Differential equations were solved using Labview 7.0 (National Instruments, Austin, TX). The data were fitted to the differential equations using the Levenberg-Marquardt method (16).
Mathematical model for TREC analysis. To help interpret the data obtained in this study, we propose a slight generalization of a mathematical model originally described by Hazenberg et al. (21). The number of nave T cells/μl (T) and the number of TREC+ T cells/μl (T+) are governed by the following equations:
(1)
Nave CD4+ or CD8+ T cells, T, in the periphery receive a constant input from the thymus, , proliferate at rate p (day–1) and disappear at a rate d (day–1). TREC+ T cells, T+, appear at a lower rate of thymic production (f) and disappear from the same compartment as nave T cells at a rate d. We assume that nave T cells can proliferate without losing their nave phenotype as has been reported (45-47). Since TRECs do not replicate during cell mitosis (9), proliferation of TREC+ T cells decreases the fraction of TREC+ T cells per nave T cell (T+/T).
Since changes in T and T+ occur very slowly, the steady-state values obtained by equation 1 can be used to analyze how changes in the parameters (, d, p, and f) would affect changes in the number of nave cells, T, the number of TREC+ cells, T+, and the fraction of TREC+ T cells per nave T cell (T+/T):
(2)
For the generalization in equation 1, we assume that the disappearance rate, d, is a composite of two different factors: d = A + R, where A represents the death rate of nave cells and R is the rate of priming of nave cells into memory cells (7, 23). We also assume that the increase in the proliferation rate of nave T cells during chronic HIV-1 infection (p p + ) is largely explained by the increase in the rate of priming of nave T cells into memory cells (R R + ). Thus, the generalized model has been simulated to predict the changes in nave T-cell counts, TREC+ T cells, and the fraction of TRECs per nave T cell after HAART, when the perturbation of the quasi-steady state induced by the administration of HAART is described mathematically by the reduction of .
Mathematical model for the kinetics of BrdU-labeled nave T cells. To describe the in vivo kinetics of BrdU-labeled nave T cells in the blood after a 30-min BrdU infusion, we used the following semiempirical equation (29):
(3)
The pool of labeled cells in the blood, L, is refilled at a constant rate s until time , and labeled cells disappear from the pool of nave T cells in the blood with a disappearance rate d. In developing this semiempirical model, it is assumed that lymphoid tissue serves as an effective source of labeled cells that are distributed to the blood until equilibration is reached (time ), at which point the effective source ceases to affect changes in the concentration of labeled cells (29). Because BrdU+ chromosomes segregate independently into daughter cells, labeled cells that have divided will still be counted as BrdU+ cells, as long as the intensity of BrdU in each cell is higher than the threshold of flow cytometric detection. In human lymphocytes we estimate that the BrdU intensity decreases below the detection threshold after two to three divisions (data not shown). Thus, for highly proliferating cells equation 3 is approximately similar to a model where d consists only of the disappearance rate of labeled cells (29). For more slowly proliferating cells, such as nave T cells, a possible contribution of proliferation can be taken into account by replacing d with d – p, which is similar to an equation proposed by Debacq and colleagues (4). When d is small (<<1) the solution of equation 3 for t is given by , or the fraction of labeled cells increases approximately linearly over time with a rate s. Thus, the solution of equation 3 used in this analysis is given by
(4)
Moreover, if the peak of labeling is reached at similar times between recipients, the value of the fraction of BrdU-labeled cells at the peak will correlate with s and, thus, with proliferation rates (4).
For BrdU labeling, nave CD4 cells are defined as CD45RO–, as additional markers were not utilized in these analyses. For CD4 cells, this is a good approximation of true nave cells (6, 30).
RESULTS
Effect of HAART on viral load, T-cell counts, and T-cell proliferation. Sixteen of 20 group 1 patients had suppression of viral loads to <50 copies/ml at time point 1, and 18 of 20 had the same result at time point 2. Subjects had a mean CD4 count increase compared to a baseline of 144 cells/μl (range, 33 to 515 cells/μl) at time point 1 (6 months) and 286 cells/μl (range, 17 to 643 cells/μl) at time point 2 (18 months; P < 0.01) (Table 3).
Statistically significant increases in nave CD4+ T cells were observed at both time points 1 and 2 compared to baseline (from 87 to 127 and 164 cells/μl, respectively; P < 0.01), while statistically significant increases in nave CD8+ T cells were only observed at time point 2 (from 114 to 175 cells/μl; P < 0.05) (Table 3). At baseline, nave CD4+ T-cell numbers were lower than nave CD8 T-cell numbers (P < 0.05) but not at time point 1 or time point 2. Thus, the increase in nave T cells from baseline was higher for nave CD4+ T cells than nave CD8+ T cells at time point 1 (mean increase, 55 versus 24 cells/μl; P < 0.05) and at time point 2 (mean increase, 104 versus 68 cells/μl; P < 0.05). The additional changes between time points 1 and 2 in nave T-cell numbers did not differ between nave CD4+ and nave CD8+ T cells.
Significant inverse Spearman rank correlations were observed between baseline nave T-cell counts and the relative change (n-fold) in nave T cells for both CD4+ and CD8+ T cells at time point 1 ( = –0.53, P < 0.05 for CD4; = –0.48, P < 0.05 for CD8) and time point 2 ( = –0.82, P < 0.01 for CD4; = –0.64, P < 0.01 for CD8).
Significant declines in Ki67 expression were observed in both CD4+ and CD8+ nave and central memory T cells at time point 1 (P < 0.01) (Table 3). Before HAART, a median of 3.8% of nave CD4 cells and 5% of nave CD8 cells were Ki67+, with the difference between the means reaching borderline statistical significance (P = 0.05). At time point 1 these numbers decreased to 1.7 and 2.3%, respectively (P < 0.01). The additional observed declines at time point 2 to 1.4% and 1.3%, respectively, were not statistically significant (P > 0.05, between time points 1 and 2). Moreover, the relative change (n-fold) in percent Ki67+ nave CD4+ T cells between time point 0 to time point 1 inversely correlated with the corresponding relative change in CD4+ nave T-cell counts at the same time points (Spearman rank correlation, = –0.54, P = 0.017). No similar correlation was observed between time points 1 and 2 or between any time points for nave CD8+ T cells.
Greater relative change (n-fold) in TRECs/μl than nave T cells/μl in CD4+ and CD8+ T cells after HAART. A statistically significant increase in TRECs/μl from time point 0 to time point 1 (P < 0.01) was seen for both nave CD4+ T cells and nave CD8+ T cells (Fig. 1). A greater relative increase (n-fold) was seen in the number of TRECs/μl than in the nave T-cell counts in both CD4+ (median relative increases of 2.2- and 1.7-fold, respectively, P < 0.01) and CD8+ T-cell pools (1.4- and 1.2-fold, respectively, P < 0.01) at time point 1. This observation is equivalent to an increase in the fraction of TRECs, i.e., TRECs per million nave T cells, after initiation of HAART (Fig. 1), consistent with previously reported data (21, 32). The relative increase was observed to be lower for TRECs/μl than for nave T-cell counts from time point 1 to time point 2, with this difference being statistically significant for CD4+ T cells (P < 0.05) but not for CD8+ T cells. No significant correlations were observed between the relative change in percent Ki67+ nave CD4+ T cells at time point 1 to time point 0 and the relative change in TREC+ cells at the same time points (for CD4+, n = 14, = –0.03, and P = 0.47; for CD8+, n = 16, = –0.05, and P = 0.41). The lack of correlation between changes in the proliferation of nave T cells and the change in the fraction of TREC+ cells per nave T cell suggests that it is not the extent of reduction in proliferation per se that can explain the greater relative change in TRECs/μl versus nave T-cell counts, but the latter does probably depend on changes in both the proliferation and disappearance rate of the nave pool after the initiation of HAART.
Thymic indices did not significantly change during HAART and did not correlate with changes in TRECs. CT scoring was carried out by two radiologists blinded to clinical and laboratory data. The median baseline score was 1.5 (0.5 to 3.0). Median scores remained constant at time point 1 (1.5; range, 0.5 to 4) and at time point 2 (1.5; range, 0.5 to 4), with a median change for both time points of 0 (Table 3). When scores were clustered according to subjects with greater than and less than or equal to the median CD4 increase of 96 cells/μl at time point 1 and 257 cells/μl at time point 2, differences again were not significant. By volumetric analysis, the median score prior to initiating HAART was 1.58 cm3 (0.75 to 11.9). At time point 1 there was a median increase of 0.17 cm3 (–4.24 to 9.33); an additional period of approximately 12 months of HAART resulted in a median change of –0.04 cm3 (range, –5.27 to 6.5). None of these changes was statistically significant (Table 3). Baseline nave CD4+, but not nave CD8+, T-cell counts positively correlate with thymic volume (data not shown). An inverse statistically significant correlation was observed between the relative change (n-fold) in nave T cell counts and baseline CT scores for nave CD4+ T cells ( = –0.52, P < 0.05) but not for nave CD8+ T cells ( = –0.08, P > 0.05). Partial Spearman rank correlation analysis that included age or baseline nave T-cell counts did not qualitatively change the statistical significances of the above correlations. Similar patterns of correlations were observed when CT scores were replaced by thymic volumes. In addition, thymic volume changes did not correlate with changes in any TREC parameters.
Analysis of the model. A possible explanation for the greater relative change in TRECs/μl versus nave T-cell counts can be provided by a model of nave T-cell dynamics similar to the one described in equation 1 with the disappearance rate, d, a composite of two different factors: d = A + R, where A represents the death rate of nave cells and R is the rate of priming of nave cells into memory cells (7, 23). Assuming that the increase in the proliferation rate of nave T cells during chronic HIV-1 infection (p p + ) is mostly explained by the increase in the rate of priming of nave T cells into memory cells (R R + ) (Fig. 2A), nave T-cell counts will not substantially change, since the increased recruitment of nave T cells into memory cells is counterbalanced by the simultaneous increase in the proliferation rate of nave T cells (without losing their phenotypes). Conversely, if the main effect of HAART consists of reducing the rate of priming ( 0, with time), again, nave T-cell counts will be only marginally affected since the consequent decrease in recruitment (and consequent loss) of nave T cells is now counterbalanced by the simultaneous decrease in the proliferation rate. However, the reduction in proliferation will lead to an increase of the nave T-cell (and thus TREC+ T cell) life span or, equivalently, the decrease in disappearance rate d. This leads to an increase in TREC+ T cells/μl, and since the nave T-cell counts are only marginally affected, the fraction of TREC+ T cells per nave T cell is also predicted to increase.
With values of >0 this generalization predicts a greater decrease in the proliferation than the disappearance rate when a new quasi-steady state is reached following initiation of HAART (data not shown). The dynamics of nave CD8+ T cells observed in this study are consistent with this: a significant increase in the number of TREC+ T cells/μl is accompanied by a limited increase in nave T-cell counts after initiation of HAART (Table 3 and Fig. 1 and 2B). This model also predicts that the decrease in the proliferation rate after HAART is not expected to correlate with the increase in nave T-cell counts, since the latter will be only marginally affected regardless of the rate at which nave T-cell proliferation normalizes. Again, our data are consistent with this: there is no significant correlation between the increase of CD8+ nave T-cell counts and the decrease of percent Ki67+ CD8+ nave T cells between time points 0 and 1 or later.
For CD4+ nave T cells, the above model is inadequate to explain the observed dynamics. However, the transient increase in the fraction of TREC+ T cells per nave T cell and the following new steady state reached after 18 months of antiretroviral therapy can be explained by assuming that the death rate, A, is also enhanced (A A + ) during HIV-1 infection due to activation-induced cell death. The effect of HAART is again to normalize the proliferation rate and the disappearance rate of the CD4+ nave T cells. However, the latter now includes both the death rate and the priming rate, which demonstrate differential dynamics related to changes in and , respectively. A HAART-induced normalization of the death rate that is approximately 10-fold slower than the normalization of the proliferation rate adequately accounts for the changes in both TREC fractions and TREC+ T cells/μl that are observed for CD4+ nave T cells (Fig. 2C).
Previous reports have identified a positive correlation between baseline percentages of proliferating and apoptotic (terminal deoxynucleotidyltransferase-mediated UTP nick end labeling-positive cells) T cells (36). If the HIV-1-induced increases in proliferation and death rates of nave T cells ( and ) are similarly proportionally related, then the model also predicts an inverse correlation between the HAART-induced decrease in proliferation of nave T cells and the recovery of nave T-cell counts (data not shown) which we observed for CD4+ nave T cells. Thus, HIV-1-infected patients with higher proliferation and death rates of nave T cells during chronic infection would demonstrate a faster recovery of nave T-cell counts after HAART, with most of the recovery explained by normalization of the death rate rather than by normalization of the proliferation rate, which, as in the earlier model, is still counterbalanced by the simultaneous reduction in the rate of priming.
Kinetics of BrdU-labeled nave T cells prior to and after HAART. To further evaluate the effects of HAART on the proliferation and disappearance rates of nave T cells, we analyzed the kinetics of BrdU-labeled nave T cells prior to and after initiation of HAART in six HIV-1-infected patients. Figure 3 shows the theoretical curves obtained by best fitting equation 3 to the fraction of BrdU-labeled nave CD4+ T cells for each patient, before and after initiation of therapy. Because the percentages of BrdU-labeled cells within the nave T-cell population are relatively small, and thus closer to the detection threshold, this analysis resulted in a weak convergence of the nonlinear best-fitting procedure using equation 3. Thus, equation 4 was used to estimate the slope s from the time of infusion to the peak of labeling, and a single exponential decay function has been used to estimate d from the time of peak labeling. As shown in Table 4, this modeling predicts a decay in s greater than a decay in d for nave CD4+ T cells after initiation of HAART (median relative decrease [n-fold] in s of 2.3 and median relative decrease in d of 1.5; P = 0.03). Since changes in s imply changes in mean proliferation rates, a greater decay in s than d implies a more dramatic decrease in proliferation than disappearance rates of nave T cells after HAART, whether d represents the apoptotic death rate plus the rate of priming into memory T cells or the difference between the disappearance rate in the blood, d, and the proliferation rate in the blood, p (data not shown).
DISCUSSION
In this study we have analyzed the relationship between T-cell turnover, thymic function, and immune activation in HIV-1-infected patients, focusing on nave CD4+ and nave CD8+ T cells, to better understand the contribution of these various parameters to the immunologic changes seen during HIV infection and therapy. We specifically targeted a broad range of baseline viral load and CD4+ T-cell counts in order to highlight the dynamics across the spectrum of HIV infection. At baseline, naive T-cell numbers were lower in the CD4 pool compared to the CD8 pool, but no difference between the two groups was observed in the percentages of proliferating nave T cells or in the number of TRECs. Even though a dramatic decrease in proliferation is observed for both nave CD4+ and CD8+ T cells after the first 6 months of therapy, we found that nave CD4+ T-cell numbers significantly increased after initiation of HAART, but nave CD8+ T-cell numbers were only marginally affected. Baseline nave T-cell counts inversely correlated with the relative changes in nave T cells in both subpopulations of cells, and TRECs/μl increased in both subsets of cells. We also observed an increase in the fraction of TREC+ T cells per nave T cell, which is equivalent to a faster growth of TREC+ T cells versus total nave T cells, for both CD4 and CD8 cells following initiation of HAART. We did not observe significant changes of thymic volumes during time or statistically significant correlations between changes in thymic volumes and changes in nave T-cell numbers or TRECs/μl.
To provide a unified description of the concomitant changes in nave T cells (TRECs/μl and the percentage of proliferating nave T cells during HAART) that can explain the observed differences in the dynamics of these variables in the CD4 and CD8 subpopulations of T lymphocytes, we developed a mathematical model based on a generalization of a model originally described by Hazenberg et al. (21). We have generalized the disappearance rate of nave T cells as the sum of the rates of nave cells priming into memory cells and nave cell death and assumed that the increase in proliferation of nave T cells during chronic infection is primarily explained by the increase in the rate of priming of nave T cells into memory cells. This simple theoretical framework is sufficient to predict the simultaneous increases in the fraction of TREC+ T cells per nave T cell and the number of TREC+ T cells in the periphery after initiation of HAART. The model explains the dynamics of nave CD8+, but not CD4+, T cells after institution of HAART. To explain the dynamics of nave CD4+ T cells, we postulate that there is an increase in the apoptotic death rate of nave T cells during HIV-1 infection related to immune activation or to the increase in proliferation rate and that there is delayed normalization of the apoptotic death rate compared to the proliferation rate, as has been reported for total CD4+ T cells in lymph node samples of HIV-1-infected patients before and after initiation of HAART (48). This model does not require (or exclude) changes in thymic output or redistribution of T lymphocytes from the lymphoid tissue as additional mechanisms contributing to nave T-cell recovery.
The normalization of the death rates can also account for the inverse correlation between baseline nave T-cell counts and the relative change in nave T-cell counts after initiation of HAART. However, the increase in nave T cells appears to be lower for nave CD8 T cells than for the nave CD4+ T cells, which suggests independent mechanisms of peripheral normalization for the different populations. This dichotomy, observed in the dynamics of nave CD4+ T cells compared to nave CD8+ T cells after HAART, as well as the presence of an inverse correlation between baseline thymic scores and the relative change (n-fold) in nave T-cell numbers for CD4+ but not CD8+ T cells is difficult to explain solely as a result of changes in thymic output rates or trafficking effects, since these should not have differential effects on the two populations of nave T cells. Moreover, based on this model, the observed inverse correlation between baseline counts and the relative change in nave T-cell counts after HAART for both nave CD4+ and nave CD8+ T cells suggests that increases in the death rate affect both populations. However, the greater baseline depletion of nave CD4+ T cells compared to nave CD8+ T cells, together with the concomitant increase in TRECs per microliter in both compartments following therapy, suggests, as highlighted by the model, that similar mechanisms drive both nave CD4+ and nave CD8+ T cells to be primed into memory cells, but for unknown reasons the increase in the death rate is more pronounced in the CD4+ than in the CD8+ nave T-cell populations. Interestingly, Li et al. have recently shown that, at least in the settings of acute simian immunodeficiency virus infection, higher levels of Fas- and Fas ligand-mediated apoptosis are observed within CD4+ but not CD8+ T lymphocytes in the lamina propria, which may result from massive exposure of CD4+ T cells to virion gp120 (33). Our data suggest that during chronic infection, the differential ability to tolerate similar increases in proliferation, , is an intrinsic property of each subpopulation of nave T cells. An alternative scenario in which is different between the two subpopulations would require that nave CD4 T cells have a higher proliferation rate than nave CD8 T cells to account for the relative loss in nave CD4 T cells. However, this was not observed in our data when we looked at the baseline fractions of proliferating nave CD4 and CD8 T cells.
It is important to note that this model represents an idealized situation and that deviations from this model, resulting, for instance, from the presence of nonlinear contributions of trafficking of lymphocytes or of replenishment of the peripheral pool by thymic output, might result in a situation that is far from the quasi-steady-state condition assumed in equation 1. In the latter circumstances, TREC content after initiation of HAART can potentially be affected in an unpredictable manner.
The hypothesis of a more pronounced decrease in the proliferation rate than the disappearance rate of nave T cells after HAART is supported by the kinetics of BrdU-labeled nave T cells studied longitudinally (pre- and post-HAART) that we observed in a smaller group of HIV-1-infected patients.
In principle, changes in the fraction of nave T cells carrying TRECs upon exiting the thymus (f in equation 1) might also explain a greater relative change in TREC+ T cells/μl than nave T-cell counts after initiation of HAART (32). Among the four parameters discussed in this analysis (d, p, , and f), f is the least investigated. Dion and colleagues (8) have recently shown an increase of the ratio -TRECs/-TRECs after initiation of HAART, which suggests that the newly produced nave T cells undergo more intrathymic divisions before entering the peripheral pool. Since -TRECs are produced before -TRECs, this would lead to a decrease, not an increase, in f after initiation of HAART, thus excluding changes in f as a major factor affecting the dynamics of the fraction of TREC+ T cells after initiation of HAART.
This analysis provides evidence that changes in peripheral proliferation and disappearance rates of nave T cells, rather than changes in thymic output, explain the observed dynamics of TRECs and the fraction of proliferating T cells during HAART. But what is the mechanism that drives nave T cells to proliferate faster during chronic infection The first pathogenic effect of HIV-1 infection might consist of an increase in the rate of priming of nave T cells due to a generalized state of chronic immune activation. In this scenario, the HIV-1-induced increase in the proliferation rate could serve as a compensatory (homeostatic) mechanism aimed at maintaining the nave T-cell count constant, in response to the loss of nave T cells that have been primed into memory cells. Alternatively, the p(t) expression of the proliferation rate of our model could also be modeled as a function of the number of cells, T(t), for instance, following a density-dependent law (10). Under this scenario, changes in thymic output and consequent (homeostatic) changes in peripheral proliferation of nave T cells could also explain the observed TREC dynamics, as suggested by Dutilh et al. to explain the changes of TREC content during aging (10). However, given the relatively short time frame (a few weeks) in which significant changes of proliferation and disappearance rates are seen (36) and the negligible contribution by the thymus seen in thymectomy studies (as recently described for nonhuman primates by Arron et al. [1]), as well as our thymic CT data, we feel that peripheral mechanisms (including homeostatic mechanisms) leading to changes in proliferation and disappearance rates are the major factors affecting the dynamics of TRECs. The observed dichotomy in the dynamics of nave CD4+ and nave CD8+ T cells is difficult to explain by invoking changes in thymic output as the sole mechanism responsible for changes in peripheral (homeostatic) proliferation rates of nave T cells.
If a homeostatic mechanism governs the proliferation of nave T cells, we would expect, as has been previously suggested for the entire population of T cells (24), an inverse correlation between relative change (n-fold) in nave T-cell counts after HAART and the relative changes in the percentages of proliferating nave T cells. In our study we do observe such an inverse correlation for nave CD4+ but not nave CD8+ T cells. Alternatively, the first effect of activation might consist of inducing nave T cells to proliferate faster without losing their phenotype (45-47), bringing these cells closer to the priming activation threshold which leads to an increased disappearance rate of nave T cells. In this scenario, the increase in the proliferation rate is primarily explained by the increase in the rate of priming during chronic infection. The simultaneous decrease in both the proliferation and priming rates after initiation of HAART should generate a lack of correlation between the recovery of nave T-cell counts (only marginally affected) and the decrease in the percentage of proliferating nave T cells. This latter paradigm would explain the observed lack of such correlation for nave CD8+ T cells. The additional assumption that higher levels of proliferation are associated with higher levels of apoptosis is required to explain the presence of this inverse correlation for CD4+ nave T cells. Both scenarios show consistency with a differential change in the proliferation and disappearance rates of nave T-cells. Based on current available data, it is difficult to make conclusive arguments in support of either hypothesis. Thus, further investigations are required to clarify the mechanisms responsible for the increased proliferation of nave T cells induced by HIV-1.
ACKNOWLEDGMENTS
We thank the patients for their willingness to participate in these studies, the physicians, nurses, and other members of the National Institute of Allergy and Infectious Diseases-Clinical Center HIV/AIDS Program for their support of these studies, and Betty Wise for performing the quantitative thymic measurements.
This research was supported in part by the Intramural Research Program of the National Institutes of Health, National Institute of Allergy and Infectious Diseases, and Warren Grant Magnuson Clinical Center.
REFERENCES
Arron, S. T., R. M. Ribeiro, A. Gettie, R. Bohm, J. Blanchard, J. Yu, A. S. Perelson, D. D. Ho, and L. Zhang. 2005. Impact of thymectomy on the peripheral T cell pool in rhesus macaques before and after infection with simian immunodeficiency virus. Eur. J. Immunol. 35:46-55.
Davey, R. T., Jr., N. Bhat, C. Yoder, T. W. Chun, J. A. Metcalf, R. Dewar, V. Natarajan, R. A. Lempicki, J. W. Adelsberger, K. D. Miller, J. A. Kovacs, M. A. Polis, R. E. Walker, J. Falloon, H. Masur, D. Gee, M. Baseler, D. S. Dimitrov, A. S. Fauci, and H. C. Lane. 1999. HIV-1 and T cell dynamics after interruption of highly active antiretroviral therapy (HAART) in patients with a history of sustained viral suppression. Proc. Natl. Acad. Sci. USA 96:15109-15114.
Davey, R. T., Jr., R. L. Murphy, F. M. Graziano, S. L. Boswell, A. T. Pavia, M. Cancio, J. P. Nadler, D. G. Chaitt, R. L. Dewar, D. K. Sahner, A. M. Duliege, W. B. Capra, W. P. Leong, M. A. Giedlin, H. C. Lane, and J. O. Kahn. 2000. Immunologic and virologic effects of subcutaneous interleukin 2 in combination with antiretroviral therapy: A randomized controlled trial. JAMA 284:183-189.
Debacq, C., B. Asquith, P. Kerkhofs, D. Portetelle, A. Burny, R. Kettmann, and L. Willems. 2002. Increased cell proliferation, but not reduced cell death, induces lymphocytosis in bovine leukemia virus-infected sheep. Proc. Natl. Acad. Sci. USA 99:10048-10053.
De Boer, R. J., H. Mohri, D. D. Ho, and A. S. Perelson. 2003. Turnover rates of B cells, T cells, and NK cells in simian immunodeficiency virus-infected and uninfected rhesus macaques. J. Immunol. 170:2479-2487.
De Rosa, S. C., L. A. Herzenberg, and M. Roederer. 2001. 11-color, 13-parameter flow cytometry: identification of human naive T cells by phenotype, function, and T-cell receptor diversity. Nat. Med. 7:245-248.
Di Mascio, M., I. Sereti, L. Matthews, V. Natarajan, C. Yoder, E. C. Jones, C. Chow, I. A. Sidorov, D. S. Dimitrov, J. A. Metcalf, M. A. Polis, and J. A. Kovacs. 2004. Presented at the XV International AIDS Conference, Bangkok, Thailand, 11 to 16 July 2004.
Dion, M. L., J. F. Poulin, R. Bordi, M. Sylvestre, R. Corsini, N. Kettaf, A. Dalloul, M. R. Boulassel, P. Debre, J. P. Routy, Z. Grossman, R. P. Sekaly, and R. Cheynier. 2004. HIV infection rapidly induces and maintains a substantial suppression of thymocyte proliferation. Immunity 21:757-768.
Douek, D. C., R. D. McFarland, P. H. Keiser, E. A. Gage, J. M. Massey, B. F. Haynes, M. A. Polis, A. T. Haase, M. B. Feinberg, J. L. Sullivan, B. D. Jamieson, J. A. Zack, L. J. Picker, and R. A. Koup. 1998. Changes in thymic function with age and during the treatment of HIV infection. Nature 396:690-695.
Dutilh, B. E., and R. J. de Boer. 2003. Decline in excision circles requires homeostatic renewal or homeostatic death of naive T cells. J. Theor. Biol. 224:351-358.
Fleury, S., R. J. de Boer, G. P. Rizzardi, K. C. Wolthers, S. A. Otto, C. C. Welbon, C. Graziosi, C. Knabenhans, H. Soudeyns, P. A. Bart, S. Gallant, J. M. Corpataux, M. Gillet, P. Meylan, P. Schnyder, J. Y. Meuwly, W. Spreen, M. P. Glauser, F. Miedema, and G. Pantaleo. 1998. Limited CD4+ T-cell renewal in early HIV-1 infection: effect of highly active antiretroviral therapy. Nat. Med. 4:794-801.
Fleury, S., G. P. Rizzardi, A. Chapuis, G. Tambussi, C. Knabenhans, E. Simeoni, J. Y. Meuwly, J. M. Corpataux, A. Lazzarin, F. Miedema, and G. Pantaleo. 2000. Long-term kinetics of T cell production in HIV-infected subjects treated with highly active antiretroviral therapy. Proc. Natl. Acad. Sci. USA 97:5393-5398.
Franco, J. M., A. Rubio, M. Martinez-Moya, M. Leal, E. Merchante, A. Sanchez-Quijano, and E. Lissen. 2002. T-cell repopulation and thymic volume in HIV-1-infected adult patients after highly active antiretroviral therapy. Blood 99:3702-3706.
Gaulton, G. N., J. V. Scobie, and M. Rosenzweig. 1997. HIV-1 and the thymus. AIDS 11:403-414.
Gerdes, J., H. Lemke, H. Baisch, H. H. Wacker, U. Schwab, and H. Stein. 1984. Cell cycle analysis of a cell proliferation-associated human nuclear antigen defined by the monoclonal antibody Ki-67. J. Immunol. 133:1710-1715.
Gill, P. E., W. Murray, and M. H. Wright. 1981. Practical optimization, p. 136-137. Academic Press, London, United Kingdom.
Giorgi, J. V., L. E. Hultin, J. A. McKeating, T. D. Johnson, B. Owens, L. P. Jacobson, R. Shih, J. Lewis, D. J. Wiley, J. P. Phair, S. M. Wolinsky, and R. Detels. 1999. Shorter survival in advanced human immunodeficiency virus type 1 infection is more closely associated with T lymphocyte activation than with plasma virus burden or virus chemokine coreceptor usage. J. Infect. Dis. 179:859-870.
Hazenberg, M., S. Galkina, G. Chkhenkeli, C. Stoddart, and M. McCune. 2004. Presented at the 11th Conference on Retroviruses and Opportunistic Infections, San Francisco, Calif.
Hazenberg, M. D., D. R. Clark, and F. Miedema. 1999. Tilted balance of cell renewal in HIV-1 infection. AIDS Rev. 1:67-73.
Hazenberg, M. D., D. Hamann, H. Schuitemaker, and F. Miedema. 2000. T cell depletion in HIV-1 infection: how CD4+ T cells go out of stock. Nat. Immunol. 1:285-289.
Hazenberg, M. D., S. A. Otto, J. W. Cohen Stuart, M. C. Verschuren, J. C. Borleffs, C. A. Boucher, R. A. Coutinho, J. M. Lange, T. F. Rinke de Wit, A. Tsegaye, J. J. van Dongen, D. Hamann, R. J. de Boer, and F. Miedema. 2000. Increased cell division but not thymic dysfunction rapidly affects the T-cell receptor excision circle content of the naive T cell population in HIV-1 infection. Nat. Med. 6:1036-1042.
Hazenberg, M. D., S. A. Otto, B. H. van Benthem, M. T. Roos, R. A. Coutinho, J. M. Lange, D. Hamann, M. Prins, and F. Miedema. 2003. Persistent immune activation in HIV-1 infection is associated with progression to AIDS. AIDS 17:1881-1888.
Hazenberg, M. D., S. A. Otto, A. M. van Rossum, H. J. Scherpbier, R. de Groot, T. W. Kuijpers, J. M. Lange, D. Hamann, R. J. de Boer, J. A. Borghans, and F. Miedema. 2004. Establishment of the CD4+ T-cell pool in healthy children and untreated children infected with HIV-1. Blood 104:3513-3519.
Hazenberg, M. D., J. W. Stuart, S. A. Otto, J. C. Borleffs, C. A. Boucher, R. J. de Boer, F. Miedema, and D. Hamann. 2000. T-cell division in human immunodeficiency virus (HIV)-1 infection is mainly due to immune activation: a longitudinal analysis in patients before and during highly active antiretroviral therapy (HAART). Blood 95:249-255.
Hellerstein, M., M. B. Hanley, D. Cesar, S. Siler, C. Papageorgopoulos, E. Wieder, D. Schmidt, R. Hoh, R. Neese, D. Macallan, S. Deeks, and J. M. McCune. 1999. Directly measured kinetics of circulating T lymphocytes in normal and HIV-1-infected humans. Nat. Med. 5:83-89.
Hellerstein, M. K., R. A. Hoh, M. B. Hanley, D. Cesar, D. Lee, R. A. Neese, and J. M. McCune. 2003. Subpopulations of long-lived and short-lived T cells in advanced HIV-1 infection. J. Clin. Investig. 112:956-966.
Highbarger, H. C., W. G. Alvord, M. K. Jiang, A. S. Shah, J. A. Metcalf, H. C. Lane, and R. L. Dewar. 1999. Comparison of the Quantiplex version 3.0 assay and a sensitized Amplicor monitor assay for measurement of human immunodeficiency virus type 1 RNA levels in plasma samples. J. Clin. Microbiol. 37:3612-3614.
Kaufmann, G. R., J. Zaunders, J. Murray, A. D. Kelleher, S. R. Lewin, A. Solomon, D. Smith, and D. A. Cooper. 2001. Relative significance of different pathways of immune reconstitution in HIV type 1 infection as estimated by mathematical modeling. AIDS Res. Hum. Retrovir. 17:147-159.
Kovacs, J. A., R. A. Lempicki, I. A. Sidorov, J. W. Adelsberger, B. Herpin, J. A. Metcalf, I. Sereti, M. A. Polis, R. T. Davey, J. Tavel, J. Falloon, R. Stevens, L. Lambert, R. Dewar, D. J. Schwartzentruber, M. R. Anver, M. W. Baseler, H. Masur, D. S. Dimitrov, and H. C. Lane. 2001. Identification of dynamically distinct subpopulations of T lymphocytes that are differentially affected by HIV. J. Exp. Med. 194:1731-1741.
Kovacs, J. A., S. Vogel, J. A. Metcalf, M. Baseler, R. Stevens, J. Adelsberger, R. Lempicki, R. L. Hengel, I. Sereti, L. Lambert, R. L. Dewar, R. T. Davey, Jr., R. E. Walker, J. Falloon, M. A. Polis, H. Masur, and H. C. Lane. 2001. Interleukin-2 induced immune effects in human immunodeficiency virus-infected patients receiving intermittent interleukin-2 immunotherapy. Eur. J. Immunol. 31:1351-1360.
Lempicki, R. A., J. A. Kovacs, M. W. Baseler, J. W. Adelsberger, R. L. Dewar, V. Natarajan, M. C. Bosche, J. A. Metcalf, R. A. Stevens, L. A. Lambert, W. G. Alvord, M. A. Polis, R. T. Davey, D. S. Dimitrov, and H. C. Lane. 2000. Impact of HIV-1 infection and highly active antiretroviral therapy on the kinetics of CD4+ and CD8+ T cell turnover in HIV-infected patients. Proc. Natl. Acad. Sci. USA 97:13778-13783.
Lewin, S. R., R. M. Ribeiro, G. R. Kaufmann, D. Smith, J. Zaunders, M. Law, A. Solomon, P. U. Cameron, D. Cooper, and A. S. Perelson. 2002. Dynamics of T cells and TCR excision circles differ after treatment of acute and chronic HIV infection. J. Immunol. 169:4657-4666.
Li, Q., L. Duan, J. D. Estes, Z. M. Ma, T. Rourke, Y. Wang, C. Reilly, J. Carlis, C. J. Miller, and A. T. Haase. 2005. Peak SIV replication in resting memory CD4+ T cells depletes gut lamina propria CD4+ T cells. Nature 434:1148-1152.
Liu, Z., W. G. Cumberland, L. E. Hultin, H. E. Prince, R. Detels, and J. V. Giorgi. 1997. Elevated CD38 antigen expression on CD8+ T cells is a stronger marker for the risk of chronic HIV disease progression to AIDS and death in the Multicenter AIDS Cohort Study than CD4+ cell count, soluble immune activation markers, or combinations of HLA-DR and CD38 expression. J. Acquir. Immune Defic. Syndr. Hum. Retrovirol. 16:83-92.
Lu, A. C., E. C. Jones, C. Chow, K. D. Miller, B. Herpin, D. Rock-Kress, J. A. Metcalf, H. C. Lane, and J. A. Kovacs. 2003. Increases in CD4+ T lymphocytes occur without increases in thymic size in HIV-infected subjects receiving interleukin-2 therapy. J. Acquir. Immune Defic. Syndr. 34:299-303.
Mohri, H., A. S. Perelson, K. Tung, R. Ribeiro, B. Ramratnam, M. Markowitz, R. Kost, A. Hurley, L. Weinberger, D. Cesar, M. Hellerstein, and D. D. Ho. 2001. Increased turnover of T lymphocytes in HIV-1 infection and its reduction by antiretroviral therapy. J. Exp. Med. 194:1277-1287.
Murali-Krishna, K., and R. Ahmed. 2000. Cutting edge: naive T cells masquerading as memory cells. J. Immunol. 165:1733-1737.
Natarajan, V., R. A. Lempicki, I. Sereti, Y. Badralmaa, J. W. Adelsberger, J. A. Metcalf, D. A. Prieto, R. Stevens, M. W. Baseler, J. A. Kovacs, and H. C. Lane. 2002. Increased peripheral expansion of naive CD4+ T cells in vivo after IL-2 treatment of patients with HIV infection. Proc. Natl. Acad. Sci. USA 99:10712-10717.
Roederer, M., J. G. Dubs, M. T. Anderson, P. A. Raju, L. A. Herzenberg, and L. A. Herzenberg. 1995. CD8 naive T cell counts decrease progressively in HIV-infected adults. J. Clin. Investig. 95:2061-2066.
Sachsenberg, N., A. S. Perelson, S. Yerly, G. A. Schockmel, D. Leduc, B. Hirschel, and L. Perrin. 1998. Turnover of CD4+ and CD8+ T lymphocytes in HIV-1 infection as measured by Ki-67 antigen. J. Exp. Med. 187:1295-1303.
Schacker, T., S. Little, E. Connick, K. Gebhard, Z. Q. Zhang, J. Krieger, J. Pryor, D. Havlir, J. K. Wong, R. T. Schooley, D. Richman, L. Corey, and A. T. Haase. 2001. Productive infection of T cells in lymphoid tissues during primary and early human immunodeficiency virus infection. J. Infect. Dis. 183:555-562.
Sereti, I., K. B. Anthony, H. Martinez-Wilson, R. Lempicki, J. Adelsberger, J. A. Metcalf, C. W. Hallahan, D. Follmann, R. T. Davey, J. A. Kovacs, and H. C. Lane. 2004. IL-2-induced CD4+ T-cell expansion in HIV-infected patients is associated with long-term decreases in T-cell proliferation. Blood 104:775-780.
Sousa, A. E., J. Carneiro, M. Meier-Schellersheim, Z. Grossman, and R. M. Victorino. 2002. CD4 T cell depletion is linked directly to immune activation in the pathogenesis of HIV-1 and HIV-2 but only indirectly to the viral load. J. Immunol. 169:3400-3406.
Sprent, P. 1993. Applied nonparametric statistical methods, 2nd ed. Chapman and Hall, New York, N.Y.
Tough, D. F., and J. Sprent. 1994. Turnover of naive- and memory-phenotype T cells. J. Exp. Med. 179:1127-1135.
Unutmaz, D., F. Baldoni, and S. Abrignani. 1995. Human naive T cells activated by cytokines differentiate into a split phenotype with functional features intermediate between naive and memory T cells. Int. Immunol. 7:1417-1424.
Unutmaz, D., P. Pileri, and S. Abrignani. 1994. Antigen-independent activation of naive and memory resting T cells by a cytokine combination. J. Exp. Med. 180:1159-1164.
Zhang, Z. Q., D. W. Notermans, G. Sedgewick, W. Cavert, S. Wietgrefe, M. Zupancic, K. Gebhard, K. Henry, L. Boies, Z. Chen, M. Jenkins, R. Mills, H. McDade, C. Goodwin, C. M. Schuwirth, S. A. Danner, and A. T. Haase. 1998. Kinetics of CD4+ T cell repopulation of lymphoid tissues after treatment of HIV-1 infection. Proc. Natl. Acad. Sci. USA 95:1154-1159.(Michele Di Mascio, Irini )