Class-modeling analysis reveals T-cell homeostasis disturbances involved in loss of immune control in elite controllers
BMC Medicine volume 16, Article number: 30 (2018)
Despite long-lasting HIV replication control, a significant proportion of elite controller (EC) patients may experience CD4 T-cell loss. Discovering perturbations in immunological parameters could help our understanding of the mechanisms that may be operating in those patients experiencing loss of immunological control.
A case–control study was performed to evaluate if alterations in different T-cell homeostatic parameters can predict CD4 T-cell loss in ECs by comparing data from EC patients showing significant CD4 decline (cases) and EC patients showing stable CD4 counts (controls). The partial least-squares–class modeling (PLS-CM) statistical methodology was employed to discriminate between the two groups of patients, and as a predictive model.
Herein, we show that among T-cell homeostatic alterations, lower levels of naïve and recent thymic emigrant subsets of CD8 cells and higher levels of effector and senescent subsets of CD8 cells as well as higher levels of exhaustion of CD4 cells, measured prior to CD4 T-cell loss, predict the loss of immunological control.
These data indicate that the parameters of T-cell homeostasis may identify those EC patients with a higher proclivity to CD4 T-cell loss. Our results may open new avenues for understanding the mechanisms underlying immunological progression despite HIV replication control, and eventually, for finding a functional cure through immune-based clinical trials.
Human immunodeficiency virus-1 (HIV-1) infection is mostly characterized by ongoing viral replication leading to a progressive decline of CD4+ T cells and eventually progression to clinical AIDS. Nevertheless, there is a small subset of HIV-1+ subjects who are able to spontaneously control HIV-RNA replication in plasma to levels below the limit of detection in the absence of antiretroviral therapy. They are called elite controllers (ECs) and represent less than 1% of the total of HIV-infected patients . ECs have been considered as a model of spontaneous functional cure and the role of host genetic factors, innate and adaptive immune responses, viral characteristics, and their interactions have been analyzed in several studies [1,2,3,4,5].
However, ECs constitute a heterogeneous group in terms of genetic and immunologic characteristics , and no single mechanism has been described as being responsible for controlling viral replication [5,6,7]. These patients are also heterogeneous from the clinical point of view, since some of them show either virologic and/or immunologic progression [8,9,10,11]. Studies with large cohorts of EC patients have reported that as many as 25% of them experience virological progression (losing the ability to maintain HIV replication control) , and around 20% experience immunological progression (a decline in CD4 counts) even in the absence of virological progression [8, 9], as we have also previously demonstrated in our cohort of HIV controllers of the Spanish AIDS Research Network (the ECRIS cohort) . However, the mechanisms accounting for this CD4 T-cell loss in ECs remain largely unknown.
Given the relevance of different alterations of T-cell homeostasis in the CD4 depletion observed in HIV patients with uncontrolled viral replication [12,13,14], some of these same mechanisms might be operating in EC patients experiencing CD4 T-cell loss. In this regard, previous studies have shown that even though EC patients are able to suppress HIV replication to an undetectable level, they present different T-cell homeostasis disturbances compared to uninfected controls and/or to HIV non-controllers with suppression of HIV replication mediated by combination antiretroviral therapy (cART) [15,16,17,18,19,20], suggesting that T-cell homeostasis disturbances are still operating in EC patients and could potentially be involved in the immunological progression.
To test this hypothesis, and with the final aim of finding predictors of immunological progression, herein we have comprehensively compared several important parameters of T-cell homeostasis in EC patients with long-lasting HIV replication control but experiencing immunological progression with those for EC patients showing no progression. The partial least-squares–class modeling (PLS-CM) methodology was applied as a statistical modeling technique to discriminate between the two groups of patients, and as a predictive model.
This is a case–control study including adult patients with chronic HIV infection and with long-term spontaneous control of HIV replication (EC patients) during the whole study follow-up period, recruited from the cohort of HIV controllers of the Spanish AIDS Research Network (ECRIS cohort) launched in 2013. ECRIS is an open multicenter cohort of HIV controller patients. The data are from the Spanish Long-Term Non-Progressors (LTNP) cohort and the Spanish AIDS Research Network (CoRIS) cohort , and from different clinical centers in Spain. To be included in the ECRIS cohort, EC patients were defined as asymptomatic individuals with at least three consecutive plasma HIV viral load (pVL) determinations below the detection limit (pVL < 50 copies/mL) during at least 12 months of follow-up, in the absence of any cART. The characteristics of the ECRIS cohort have been described in detail elsewhere . Only those EC subjects maintaining long-term control of HIV replication (more than 3 years) with regular immunovirological follow-ups (CD4 counts and pVL) and with available baseline cryopreserved cellular samples were included in the study. Baseline was the beginning of the follow-up period as an EC patient.
Selection of EC study groups
The relevant CD4 data for each patient were measured during the follow-up period if the patient maintained their EC status. All of these data were used in a linear regression analysis to estimate the CD4 slope for the patient. According to the CD4 slope, the study population was categorized into two groups: EC controls, patients with stable CD4 counts during the follow-up period (CD4 count slope not significantly different from zero) and EC cases, patients whose CD4 counts declined during the follow-up period (a statistically significant and negative CD4 count slope). A total of 36 EC patients were included: 22 being EC controls and 14 being EC cases (Additional file 1 and Fig. 1). Moreover, 14 age- and sex-matched healthy controls (HIV-uninfected individuals) were included as a healthy control group (HCs).
Anonymized samples from both groups of patients were kindly provided by the HIV BioBank, which is integrated into the Spanish AIDS Research Network (RIS) . Anonymized data were taken from the ECRIS cohort database (see Additional file 1), according to the sample assignment agreement reference 19/15 dated 25 August 2015. All analyses were done with cryopreserved peripheral blood mononuclear cells (PBMCs), and were performed when both groups of EC patients (cases and controls) presented similar levels of CD4 counts, that is prior to CD4 T-cell loss in EC cases. EDTA-anticoagulated blood was obtained by venipuncture. PBMCs were immediately isolated by density gradient centrifugation using Ficoll-Hypaque (Sigma Chemical Co., St. Louis, MO), frozen in fetal calf serum plus 10% dimethyl sulfoxide, and stored in liquid nitrogen in the Spanish AIDS Research Network HIV BioBank. The viability of thawed PBMCs was always greater than 85%.
We used a comprehensive approach of simultaneously measuring different T-cell parameters involved in HIV pathogenesis: differentiation stage; recent thymic emigrants (RTEs); activation, exhaustion, and apoptosis of T memory stem cells (Tscm cells); T regulatory cells (Treg cells); replicative senescence; and response to homeostatic cytokines. We selected a total of 16 different monoclonal antibodies distributed in two different panels of 11 and 13 antibodies (plus a cell viability marker), which in combination defined 350 different unique T-cell subsets that were analyzed by multiparametric flow cytometry. A complete list of all monoclonal antibodies and fluorochromes used in the study as well as antibodies included in each of the two staining panels is shown in Additional file 2. A minimum of 2 million of PBMCs from each patient were stained with each antibody panel, and a detailed description of staining conditions and definitions of T-cell subsets analyzed is given in Additional file 3.
The main characteristics of the study population and the different parameters evaluated are expressed as median and interquartile range. Comparisons between groups were done using a Mann–Whitney U test. The two-tailed p values were considered as significant only when they were lower than 0.05. All statistical analyses were performed using the SPSS software version 15 (SPSS Inc., Chicago, IL, USA). PLS-CM was employed to discriminate between the different groups of subjects included in the study (control patients, case patients, and HCs). Two PLS models were generated: one to discriminate between EC patients and HCs (PLS-1) and another to discriminate between cases and control groups of EC patients (PLS-2). In these models, 350 immunological parameters were used as predictor (explicative) variables of the model, and the subjects’ group was the predicted (response) binary variable. For PLS-1, a value of 0 was assigned to EC patients and a value of 1 was assigned to HCs. For PLS-2, a value of 0 was assigned to EC controls and a value of 1 to EC cases.
In the first step, PLS was used to select a minimum set of predictor variable to reduce the original 350 immunological variables to simplify the analysis without losing any information. The parameter employed to make this selection, called variable influence on projection (VIP), summarizes the importance of the predictor variable on the model taking into account the amount of explained variance of the predicted (response) variable.
In the second step, a PLS model was built using the selected set of predictor variables grouped together to generate a minimum set of latent variables (LVs) that are not directly observed (measured) but rather inferred (through a mathematical model) from other variables that are observed (directly measured). Each LV is a linear combination of the original set of observable variables. A clear advantage of using LVs is that it reduces the dimensionality of data. The number of LVs included in the PLS model is optimized to maximize the explained variance of the X (predictor variables) and Y (response variable) components of the model. Applying this PLS model, a predicted value for the response variable was estimated for each sample and for these predicted values, a probability distribution was generated. For this distribution, the probabilities α and ß for type I and II errors of the hypothesis test, respectively, were calculated. The sensitivity [(1- α) × 100] and specificity [(1- ß) × 100] of the PLS model were used to assign each sample to its correct group. A detailed explanation of PLS-CM is given in Additional file 4.
The CD4 slope during the follow-up period for controls and cases is shown in Table 1. Of the 36 EC patients included in the study, 22 were considered controls (stable CD4 counts during the follow-up period) and 14 were considered cases (significant CD4 decline during the follow-up period). All immunological parameters were measured prior to CD4 T-cell loss in EC cases. The evolution of CD4 counts during the follow-up period for controls and cases is shown in Fig. 1. The main characteristics of the controls and cases are shown in Table 1. There were no significant differences between control and case groups in terms of age, median time since diagnosis, CD4 count at the beginning of follow-up, CD4 count at the commencement of the study, percentage of patients infected with hepatitis C virus, percentage of patients presenting with non-AIDS defining events during the follow-up period, and number of HIV-RNA blips during the follow-up period.
PLS class modeling of EC patients versus HC subjects (PLS-1)
First a PLS model was generated to discriminate between EC patients and HIV-seronegative HC subjects. A data matrix with 50 samples (36 EC patients and 14 HCs) and 350 variables was used to fit a first PLS model without any data preprocessing. A binary variable took the value 0 if the sample was from an EC patient and 1 if it was from a HC. This model let us choose 93 variables with VIP values higher than 1 so that we could simplify the analysis without losing any information. A new PLS model (PLS-1) with the same criteria and a reduced data matrix with 93 variables was built. Three LVs were included in the final PLS model, and the percentage of explained variance of the response variable (y) and predictor variables (X) is shown in Table 2 for PLS-1. A normal distribution for each class (EC patients and HCs) was fitted with the values of the response variable calculated by the PLS model because the p values of the normality tests were greater than 0.1. These density functions for each class are plotted in Fig. 2. The mean and standard deviation of the calculated response variable for the PLS-1 model for the EC patient and HC classes are, respectively, 0.075 ± 0.219 and 0.810 ± 0.146 (Fig. 2).
By choosing a threshold to differentiate between EC patients and HC subjects, we can see in Fig. 2 that the small area under the density function of the class of EC patients on the right of the threshold is the probability α of rejecting the null hypothesis (H0: the sample belongs to the EC patient class) when it is true. On the other hand, the small area under the density function of the class of HCs on the left of the threshold is the probability β of accepting the null hypothesis (H0) when it is false. By changing the threshold, a risk curve of type α error versus type ß error can be generated and from this we can estimate the pair of equal values of sensitivity [(1-α) × 100] and specificity [(1- ß) × 100] of the model (Fig. 3). For the PLS-1 model, these values were both 98%. To visualize the ability of the PLS model to discriminate between EC patients and HCs, Fig. 4 (left) shows the projection in a plane of the scores for LV2 and LV3 (rotated 25° counterclockwise) for each of the 50 individuals of the study (36 EC patients and 14 HCs). This plot has a vertical line (a rotated LV2 value) that perfectly separates EC patients (on the left) and HC subjects (on the right).
PLS class modeling of EC cases versus EC controls (PLS-2)
PLS class modeling was applied to discriminate between EC cases (EC patients showing a CD4 decline during the follow-up period) and EC controls (EC patients with stable CD4 counts during the follow-up period). Firstly, by fitting a PLS model with samples from 36 patients (22 EC controls and 14 EC cases), 77 variables from the 350 original ones with a VIP value higher than 1 were selected. No data preprocessing was used and a binary response variable took the value 0 if the sample was from an EC control and the value 1 if it was from an EC case. The same procedure as used for PLS-1 was followed.
A PLS regression (PLS-2) with the same previous criteria and the reduced data matrix with 77 variables was built with four LVs. The percentages of explained variance by the model of the X and y variables are shown in Table 2 (under PLS-2). We checked that the p values of the normality tests were greater than 0.10 and therefore, a normal distribution for each class (EC cases and EC controls) was fitted with the values of the response variable calculated by the PLS model. The normal density functions for each class were plotted (Fig. 2). The mean values and standard deviation of the response variable calculated by the PLS model for the EC case and EC control classes were respectively 0.096 ± 0.247 and 0.846 ± 0.162 (Fig. 2).
By changing the threshold, a risk curve of type α versus type ß errors was generated and from this we estimated the pair of equal values of sensitivity [(1-α) × 100] and specificity [(1- ß) × 100] of the model (Fig. 3). For the PLS-2 model, these values were both 97%. To visualize the ability of the PLS-2 model to discriminate between EC cases and EC controls, Fig. 4 shows the projection in a plane of the scores for LV2 and LV3 rotated 49° counterclockwise. This plot shows has a vertical line (a rotated LV2 value) that perfectly separates EC controls (on the left) and EC cases (on the right).
The PLS model shows that the immunological variables have information with sufficient specificity and sensitivity to discriminate between the EC controls and the EC cases (this is also true for discriminating between EC patients and HC subjects). There is an underlying structure responsible for this capacity of discrimination that is reducible to two LVs. However, the immunological variables are all on the same scale (percentage) and it is not appropriate to normalize them since the variables that are noise in the discrimination acquire the same magnitude as those truly responsible for it and may be incorrectly included in the model. The consequence of not normalizing the data is twofold: (1) one of the LVs is linked to the distance between the vector formed by the mean values of the predictor variables and the origin of coordinates and (2) the immunological variables with a smaller range may not be selected due to a numerical effect. In addition, PLS tends to eliminate highly correlated variables because they do not carry independent information, although their immunological meaning may be relevant and different. Consequently, choosing the immunological variables for which the difference between classes is perceived is done by evaluating the individual distribution of their values, as in the Mann–Whitney U test.
T-cell parameters significantly altered in EC patients with respect to healthy subjects
Of the 93 immunological variables selected by the PLS-1 model to differentiate EC patients from healthy subjects, 43 presented statistically significant differences when EC patients and healthy subjects were compared with the Mann–Whitney test. Since the immunological variables are all on the same scale and it is not appropriate to normalize them, the PLS-CM model did not normalize the variables, and for this reason the immunological variables with a smaller range may not have been selected. Nevertheless, those variables could have some relevant information. Thus, using the Mann–Whitney U test, it was possible to evaluate those variables significantly different between EC patients and healthy subjects but that were not selected by the PLS model. In fact, there were another 36 immunological variables significantly different between EC patients and healthy subjects that were not selected by the PLS model due to a numerical effect. A complete list of immunological variables with significant differences is shown in Additional file 5.
Figure 5 shows levels for the most relevant variables in EC patients and healthy subjects. Among CD4 T-cell subsets, EC patients showed lower levels of Tscm and Treg subsets and higher levels of the CD45RA-CD27-CCR7- subset compared to healthy subjects. Interestingly, expression of activation marker CD38 (but not of HLADR) was also lower in total CD4 T cells and in naïve CD4 T cells of EC patients. In contrast, co-expression of activation markers CD38 and HLADR was higher in Tscm cells of EC patients. The proportion of naïve CD4 cells lacking expression of CD28 was also higher in EC patients.
Among CD8 T-cell subsets, EC patients showed lower levels of naïve RTE (naïve cells expressing CD31) and Tscm subsets compared to healthy subjects. However, co-expression of activation markers CD38 and HLADR in total CD8+ cells and in effector (CD45RA + CD27-CCR7-) subset was higher in EC patients. Expression of CD95 in activated (either CD38+ or HLADR+) subsets of CD8 cells was also higher in EC patients. Lastly, the proportion of CD8 cells lacking CD28 expression was also higher in EC patients (Fig. 5).
T-cell parameters associated with loss of immunological control in EC patients
Of the 77 immunological variables selected by the PLS model to differentiate between EC cases (EC patients showing a significant CD4 count decline) and EC controls (EC patients with a stable CD4 count), 12 presented statistically significant differences between EC cases and EC controls using the Mann–Whitney test. Interestingly, all of them were variables with the highest coefficients (loadings) in the PLS model. The absolute value of the loading of each immunological variable in the PLS model is a measure of the weight of the variable in the model and thus, the higher the loading the higher the probability of significant differences between the two groups (EC cases and EC controls). In fact, for the entire set of 77 immunological variables selected by the PLS model, there was a significant inverse correlation between the absolute value of the loadings and the p values obtained by the Mann–Whitney test (Additional file 6), meaning that the higher the absolute value of the loading, the higher the probability of a statistically significant (p < 0.05) difference.
Moreover, using the Mann–Whitney U test, it was possible to evaluate those variables significantly different between EC cases and EC controls that were not selected by the PLS model due to a numerical effect. This approach allowed us to detect another nine immunological variables with levels significantly different between EC cases and EC controls. A complete list of immunological variables significantly different by the Mann–Whitney test between EC cases and EC controls is shown in Additional file 7.
Figure 6 shows levels for the most relevant variables for EC cases and EC controls. Among CD4 subsets, EC cases presented lower levels of CD38 expression on the effector (CD45RA + CD27-CCR7-) subset of CD4 cells. In contrast, the level of PD1 expression on total and on central memory (CD45RA-CD27 + CCR7+) CD4 cells was higher in EC cases compared to EC controls. Among CD8 T-cell subsets, compared to EC controls, EC cases showed lower levels of naïve (CD45RA + CD27 + CCR7+) and of RTE (naïve cells expressing CD31) subsets. In contrast, the level of the effector (CD45RA + CD27-CCR7-) subset was higher in EC cases as was the level of senescence (CD28-CD57+) in total CD8 cells and in naïve CD8 cells. Lastly, expression of CD95 in the effector subset of CD8 cells was lower in EC cases.
Overall, for those T-cell subsets that were significantly different in EC patients with respect to HCs, the levels for EC controls were intermediate between those found in EC cases and HCs (Fig. 7). Thus, levels in EC controls were more similar to levels in HCs than were levels in EC cases, and for some of these subsets, the levels for EC controls and HCs were not significantly different. This was the case for levels of the CD45RA-CD27-CCR7- subset of CD4 cells and for the levels of CD28 + CD57- subset of CD8 cells (Fig. 7).
The present study was designed to test the hypothesis that perturbations in immunological parameters related to T-cell homeostasis may be associated with CD4 T-cell loss in EC patients maintaining long-term undetectable levels of plasma HIV viremia. The main findings of our study are: (1) there is an important perturbation of T-cell homeostasis in EC patients despite long-term undetectable viremia; (2) based on the T-cell parameters analyzed, a PLS-CM statistical approach was able to discriminate and to predict, with high sensitivity and specificity, EC patients experiencing immunological progression from those EC patients with a stable CD4 count; and (3) immunological parameters associated with the loss of CD4 T cells in EC patients included mainly those related to CD8 T-cell homeostasis.
Despite the persistent control of HIV replication to undetectable levels, several alterations related to a state of systemic inflammation have been recognized in ECs compared to seronegative subjects and to HIV patients with cART-mediated control of HIV replication [23,24,25], and this has been associated with an increased rate of morbidity and mortality in EC patients [24,25,26]. Our study demonstrates that immunological alterations in EC patients are not limited to markers of systemic inflammation and extend to several cellular markers of T-cell homeostasis. Although the EC patients included in our study had maintained long-term undetectable plasma HIV viremia, they clearly departed from seronegative healthy subjects (HCs) regarding the distribution of T-cell subsets. As many as 93 immunological variables were necessary to discriminate between EC patients and HCs and many of these presented levels significantly different between ECs and HCs, suggesting there is an important perturbation of T-cell homeostasis in EC patients despite long-term undetectable viremia.
Different processes involved in T-cell homeostasis were altered in ECs compared to healthy subjects. Among them, naïve and RTE subsets, stem cell memory subsets, regulatory T cells, and expression of the costimulatory receptor CD28 were significantly decreased, whereas activation and expression of CD95 were significantly increased. Some of these perturbations have already been reported in different cohorts of ECs, such as alterations in naïve/memory subsets [17, 19], increased activation [15, 18, 27], increased senescence , and decreased regulatory T cells . However, EC cohorts employed in these previous studies had varying degrees of HIV-associated immune suppression (as reflected in CD4 counts) and different cutoffs of plasma viremia, meaning that some alterations reported in these studies could be influenced by CD4 counts  or low-level plasma HIV viremia . Since all EC patients included in our study had high CD4 counts at the time of the immunophenotypic analysis, our results demonstrate that perturbation of T-cell homeostasis is already present in ECs from the early stages of HIV disease. Moreover, our study is the first reporting alterations in ECs of some important parameters of T-cell homeostasis such as Tscm cells [29, 30]. Interestingly Tscm cells were diminished in both CD4 and CD8 T cells and thus, the potential to repopulate other subsets of memory cells  may be compromised in ECs at early stages of infection.
The most relevant findings of our study are those related to T-cell parameters associated with loss of immunological control in ECs. Using a set of 77 different immunological variables, the PLS model was able to discriminate, with high sensitivity and specificity, between EC patients showing a significant CD4 decline (cases) and those with stable CD4 counts (controls). The immunological variables were assessed when both groups of EC patients presented similar high levels of CD4 counts (prior to CD4 T-cell loss in EC cases) and thus, the differences found cannot be attributed to different degrees of HIV-associated immune suppression, but rather they are characteristic of some EC subjects and likely involved in the decline of CD4 T cells observed in EC cases. Since no functional methods were employed in this study, we can speculate only about the potential mechanisms involved in the differences observed in the extent of T-cell homeostasis disturbance between EC cases and EC controls. Besides viral factors, such as the higher levels of residual viral replication in EC cases , several host-related factors could contribute to the CD4 T-cell loss and T-cell homeostasis disturbance observed in EC cases. Among them are defective thymic function , lymphoid tissue fibrosis , genetic variation in different genes related to immune response [33, 34], and a cytokine-driven inflammatory milieu .
Interestingly, the two-sample non-parametric test (Mann–Whitney U test) was not able to detect subtle but important differences between cases and controls for many of the immunological variables that were detected by the PLS model. PLS-CM is an approach that is best suited to extracting these differences employing latent factors and building a predictive model in situations with a small sample size and many predictive variables. In agreement with this higher sensitivity of PLS modeling to detect differences among a large panel of variables, the immune parameters that had significant differences by the two-sample test were those showing the highest loadings in the PLS model.
The factors involved in the immunological progression observed in some ECs that maintain virological control are largely unknown and studies addressing this issue are scarce, with only a few investigating potential associated immunological factors [10, 15, 36, 37]. In some of these studies, T-cell activation was shown to be associated with disease progression in ECs using a cross-sectional design [15, 36]. The other two studies addressing the association of T-cell activation with CD4 evolution in EC patients and using a longitudinal follow-up period similar to our study [10, 37] found discordant results. One reported increased activation associated with immunologic progression  but the other did not . Differences in the criteria used to define immunological progression may explain the discordant results. Moreover, both studies analyzed T-cell activation when CD4 counts were significantly different between EC patients showing immunological progression and those maintaining stable CD4 counts, which is in clear contrast with our study.
Our results show that several T-cell subsets showed significant differences between EC cases and controls and interestingly they were mainly related to CD8 cell homeostasis. Naïve and RTE subsets of CD8 cells were diminished in EC cases, suggesting either a defect in thymic output or an increased rate of transition of naïve cells to more differentiated subsets. Since a deficit in thymic output would be reflected in the RTE subset of both CD8 and CD4 compartments, our findings more strongly support an increased rate of differentiation. This, together with the finding of increased levels of effector CD8 cells in EC cases, supports the higher levels of residual HIV replication in EC cases. Unfortunately, we did not measure plasma HIV viremia with ultrasensitive assays and thus, we can only speculate on this possibility. However, in support of this hypothesis, previous studies have found higher levels of residual viremia in EC patients showing CD4 depletion [9, 10, 38]. Although levels of effector CD8 cells were higher in EC cases, CD8 T cells showed increased levels of the CD28-CD57+ subset, a phenotype that has been associated with replicative senescence . CD8 T cells expressing this senescent phenotype have lower cytotoxic potential [40, 41] and produce high amounts of inflammatory cytokines such as IL6 and TNFα , promoting an inflammatory state that could contribute to the CD4 depletion observed in EC cases. An interesting finding was that naïve CD8 cells of EC cases also have increased levels of senescence, suggesting that they have undergone an antigen-independent expansion without differentiation to memory cells in an attempt to maintain CD8 T-cell homeostasis. Lastly, regarding alterations of CD4 T-cell subsets, EC cases had increased levels of expression of the cell exhaustion-associated marker PD1 in total CD4 cells and in central memory subset. Since overexpression of PD1 has been associated with a lower ability to reconstitute CD4 T cells in HIV patients with cART-mediated suppression of viral replication [42, 43], this increased expression of PD1 in EC cases could be another factor contributing to CD4 depletion.
Finally, there are some caveats that deserve comment. First, studies employing EC patients are very limited in sample size given the very low frequency of this exceptional population of individuals among the whole population of HIV-infected patients. However, in our study, we addressed the potential bias of a small sample size by studying a very homogeneous study population. Also, for the data analyses we employed the PLS regression, which is a statistical technique for performing multivariate calibration, especially when there are more variables than samples. Second, although the objective of our study had no such scope, it could have been interesting to include a group of HIV-normal progressors to assess whether the underlying mechanisms of the loss of immune control in EC patients are comparable to those underlying the non-immunological control in normal HIV+ progressors patients, as well as to evaluate the power of the PLS model to discriminate between both phenomena.
In summary, the present study is the first to perform an in-depth analysis of T-cell homeostasis disturbances in a well-characterized cohort of ECs with long-term control of HIV replication and showing different evolution of CD4 counts. We used a PLS-CM statistical approach that enabled us to detect subtle but important immunophenotype differences associated with the immunological progression. Overall, our results demonstrate the existence of an important perturbation of T-cell homeostasis in EC patients despite long-term undetectable viremia. For those EC patients experiencing immunological progression, those disturbances are mainly related to CD8 T-cell homeostasis. Our results may open new avenues for understanding the mechanisms underlying immunological progression despite HIV replication control, suggesting that active pathogenic mechanisms are still present in some EC patients. These findings should be considered to improve the clinical management of this exceptional group of patients, and eventually in seeking a functional cure through immune-based clinical trials.
Acquired immune deficiency syndrome
Combination antiretroviral therapy
Ethylenediamine tetraacetic acid
Human immunodeficiency virus-1
Peripheral blood mononuclear cell
Partial least-squares–class modeling
Plasma viral load
Recent thymic emigrant
Statistical Package for the Social Sciences
- Treg cell:
T regulatory cell
- Tscm cell:
T memory stem cell
Variable influence on projection
Sajadi MM, Constantine NT, Mann DL, Charurat M, Dadzan E, Kadlecik P, et al. Epidemiologic characteristics and natural history of HIV-1 natural viral suppressors. J Acquir Immune Defic Syndr. 2009;50:403–8.
Autran B, Descours B, Avettand-Fenoel V, Rouzioux C. Elite controllers as a model of functional cure. Curr Opin HIV and AIDS. 2011;6:181–7.
Cockerham LR, Hatano H. Elite control of HIV: is this the right model for a functional cure? Trends Microbiol. 2015;23:71–5.
Shasha D, Walker BD. Lessons to be learned from natural control of HIV—future directions, therapeutic, and preventive implications. Front Immunol. 2013;4:1–8.
Deeks SG, Walker BD. Human immunodeficiency virus controllers: mechanisms of durable virus control in the absence of antiretroviral therapy. Immunity. 2007;27:406–16.
Pereyra F, Addo M, Kaufmann DE, Liu Y, Miura T, Rathod A, et al. Genetic and immunologic heterogeneity among persons who control HIV infection in the absence of therapy. J Infect Dis. 2008;197:563–71.
Walker BD, Yu XG. Unravelling the mechanisms of durable control of HIV-1. Nat Rev Immunol. 2013;13:487–98.
Madec Y, Boufassa F, Porter K, Meyer L. CASCADE collaboration. Spontaneous control of viral load and CD4 cell count progression among HIV-1 seroconverters. AIDS. 2005;19:2001–7.
Boufassa F, Saez-Cirión A, Lechenadec J, Zucman D, Avettand-Fenoel V, Venet A, et al. CD4 dynamics over a 15 year-period among HIV controllers enrolled in the ANRS French observatory. PLoS One. 2011;6:e18726.
Noel N, Lerolle N, Lecuroux C, Goujard C, Venet A, Saez-Cirión A, et al. Immunologic and virologic progression in HIV controllers: the role of viral blips and immune activation in the ANRS CO21 CODEX study. PLoS One. 2015;10:e0131922.
Leon A, Perez I, Ruiz-Mateos E, Benito JM, Leal M, Lopez-Galindez C, et al. Rate and predictors of progression in elite and viremic HIV-1 controllers. AIDS. 2016;30:1209–20.
Okoye AA, Picker LJ. CD4(+) T-cell depletion in HIV infection: mechanisms of immunological failure. Immunol Rev. 2013;254:54–64.
Paiardini M, Muller-Trutwin M. HIV-associated chronic immune activation. Immunol Rev. 2013;254:78–101.
Moir S, Chun TW, Fauci AS. Pathogenic mechanisms of HIV disease. Annu Rev Pathol. 2011;6:223–48.
Hunt PW, Brenchey J, Sinclair E, McCune JM, Roland M, Page-Shafer K, et al. Relationship between T cell activation and CD4 T cell count in VIV-seropositive individuals with undetectable plasma HIV RNA levels in the absence of therapy. J Infect Dis. 2008;197:126–33.
López M, Soriano V, Peris-Pertusa A, Rallón N, Benito JM. Elite controllers display higher activation on central memory CD8 T cells than HIV patients successfully on HAART. AIDS Res Hum Retrovir. 2011;27:157–65.
Yang Y, Al-Mozaini M, Buzón MJ, Beamon J, Ferrando-Martínez S, Ruiz-Mateos E, et al. CD4 T cell regeneration in HIV-1 elite controllers. AIDS. 2012;26:701–6.
Sajadi MM, Pulijala R, Redfield RR, Talwani R. Chronic immune activation and decreased CD4 counts associated with hepatitis C infection in HIV-1 natural viral suppressors. AIDS. 2012;26:1879–84.
Sajadi MM, Redfield RR, Talwani R. Altered T-cell subsets in HIV-1 natural viral suppressors (elite controllers) with hepatitis C infection. AIDS. 2013;27:1989–95.
Ramirez CM, Sinclair E, Epling L, Lee SA, Jain V, Hsue PY, et al. Immunologic profiles distinguish aviremic HIV-infected adults. AIDS. 2016;30:1553–62.
Sobrino-Vegas P, Gutiérrez F, Berenguer J, Labarga P, et al. [The Cohort of the Spanish HIV Research Network (CoRIS) and its associated biobank; organizational issues, main findings and losses to follow-up.] In Spanish. Enferm Infecc Microbiol Clin. 2011;29:645–53.
García-Merino I, de Las CN, Jiménez JL, Gallego J, Gómez C, Prieto C, et al. Spanish HIV BioBank. The Spanish HIV BioBank: a model of cooperative HIV research. Retrovirology. 2009;6:27.
Krishnan S, Wilson EM, Sheikh V, Rupert A, Mendoza D, Yang J, et al. Evidence for innate immune system activation in HIV type-1 elite controllers. J Infect Dis. 2014;209:931–9.
Pereyra F, Lo J, Triant V, Wei A, Buzon M, Fitch K, et al. Increased coronary atherosclerosis and immune activation in HIV-1 elite controllers. AIDS. 2012;26:2409–12.
Hsue PY, Hunt PW, Schnell A, Kalapus SC, Hoh R, Ganz P, et al. Role of viral replication, antiretroviral therapy and immunodeficiency in HIV-associated atherosclerosis. AIDS. 2009;23:1059–67.
Crowell TA, Gebo KA, Blankson JN, Korthuis PT, Yehia BR, Rutsein RM, et al. Hospitalization rates and reasons among HIV elite controllers and persons with medically controlled HIV infection. J Infect Dis. 2012;211:1692–702.
Ruiz-Mateos E, Ferrando-Martínez S, Machmach K, Viciana P, Pacheco YM, Nogales N, et al. High levels of CD57+CD28- T-cells, low T-cell proliferation and preferential expansion of terminally differentiated CD4+ T-cells in HIV-elite controllers. Curr HIV Res. 2010;8:471–81.
Hunt PW, Landay AL, Sinclair E, Martinson JA, Hatano H, Emu B, et al. A low regulatory cell response may contribute to both viral control and generalized immune activation in HIV controllers. PLoS One. 2011;6:e15924.
Gattinoni L, Speiser DE, Lichterfeld M, Bonini C. T memory cells in health and disease. Nat Med. 2017;23:18–27.
Chahroudi A, Silvestri G, Lichterfeld M. T memory stem cells and HIV: a long term relationship. Curr HIV/AIDS Rep. 2015;12:33–40.
Gattinoni L, Lugli E, Ji Y, Pos Z, Paulos CM, Quigkey MF, et al. A human memory T cell subset with stem cell-like properties. Nat Med. 2011;17:1290–7.
Sanchez JL, Hunt PW, Reilly CS, Hatano H, Beilman GJ, Khoruts A, et al. Lymphoid fibrosis occurs in long-term nonprogressors and persists with antiretroviral therapy but may be reversible with curative interventions. J Infect Dis. 2015;211:1068–75.
Dominguez-Molina B, Tarancón-Diez L, Hua S, Abad-Molina C, Rodriguez-Gallego E, Machmach K, et al. HLA-B*57 and IFNL4-related polymorphisms are associated with protection against HIV-1 disease progression in controllers. Clin Infect Dis. 2017;64:321–8.
Westrop SJ, Cocker ATH, Boasso A, Sullivan AK, Nelson MR, Imami N. Enrichment of HLA types and single-nucleotide polymorphism associated with non-progression in a strictly defined cohort of HIV-1 controllers. Front Immunol. 2017;8:746.
Li JZ, Arnold KB, Lo J, Dugast AS, Plants J, Ribaudo HJ, et al. Differential levels of soluble inflammatory markers by human immunodeficiency virus controller status and demographics. Open Forum Infect Dis. 2015;2:ofu117.
Bansal A, Sterret S, Erdman N, Westfall AO, Dionne-Odom J, Overton ET, et al. Normal T-cell activation in elite controllers with preserved CD4+ T-cell counts. AIDS. 2015;29:2245–54.
Kamya P, Tsoukas CM, Boulet S, Routy JP, Thomas R, Côte P, et al. T cell activation does not drive CD4 decline in longitudinally followed HIV-infected elite controllers. AIDS Res Ther. 2011;8:20–6.
Pereyra F, Palmer S, Miura T, Block BL, Wiegand A, Rothchild AC, et al. Persistent low-level viremia in HIV-1 elite controllers and relationship to immunologic parameters. J Infect Dis. 2009;200:984–90.
Dock J, Effros RB. Role of CD8 T cell replicative senescence in human aging and in HIV-mediated immunosenescence. Aging Dis. 2011;2:382–97.
Dagarag M, Ng H, Lubong R, Effros RB, Yang OO. Differential impairment of lytic and cytokine functions in senescent human immunodeficiency virus type 1-specific cytotoxic T lymphocytes. J Virol. 2003;77:3077–83.
Yang OO, Lin H, Dagarag M, Ng HL, Effros RB, Uittenbogaart CH. Decreased perforin and granzyme B expression in senescent HIV-1-specific cytotoxic T lymphocytes. Virology. 2005;332:16–9.
Grabmeier-Pfistershammer K, Steinberger P, Rieger A, Leitner J, Kohrgruber N. Identification of PD-1 as a unique marker for failing immune reconstitution in HIV-1-infected patients on treatment. J Acquir Immune Defic Syndr. 2011;56:118–24.
Nakanjako D, Ssewanyana I, Mayanja-Kizza H, Kiragga A, Colebunders R, Manabe YC, et al. High T-cell immune activation and immune exhaustion among individuals with suboptimal CD4 recovery after 4 years of antiretroviral therapy in an African cohort. BMC Infect Dis. 2011;11:43.
We especially thank Noelia Rodriguez for her technical assistance. We want to particularly acknowledge the patients in this study for their participation and to the HIV BioBank, which is integrated into the Spanish AIDS Research Network (RIS), and collaborating centers for the generous gifts of clinical samples used in this work. This study would not have been possible without the collaboration of the medical, nursing, staff and data managers who have taken part in the project (see Additional file 8).
This work has been partially funded by projects RD12/0017/0031, RD16/0025/0013, and SAF2015-66193-R as part of the Health Research and Development Strategy, State Plan for Scientific and Technical Research and Innovation (2008–2011 and 2013–2016) and cofinanced by the Institute of Health Carlos III (ISCIII), Sub-Directorate General for Research Assessment and Promotion and European Regional Development Fund. NR is a Miguel Servet investigator from the ISCIII (CP14/00198), Madrid, Spain. C Restrepo was funded by project RD12/0017/0031 and is currently funded by project RD16/0025/0013. M García is a predoctoral student co-funded by grant CP14/00198 and an Intramural Research Scholarship from Instituto de Investigación Sanitaria-Fundación Jiménez Díaz (IIS-FJD). The funders had no role in study design, data collection and analysis, the decision to publish, or the preparation of the manuscript.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article (and its additional files).
Ethics approval and consent to participate
To participate in the study, written informed consent was obtained from all individuals. The study protocol was evaluated and approved by the institutional ethical committee of Carlos III Hospital with reference number HCIII-P3/14.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flow diagram showing the inclusion criteria and the sequential strategy for selecting the patients included in the study. Numbers inside the boxes indicate the number of patients selected after each step of the selection process. (DOC 200 kb)
Monoclonal antibodies and fluorochromes used in the study. (DOC 49 kb)
Staining conditions for immunophenotypic analysis and definitions of T-cell subsets analyzed. (DOC 29 kb)
Partial least-squares–class modeling (PLS-CM). (DOC 33 kb)
Immunological variables (in CD4 and CD8 T-cell subsets) with significant differences between EC patients and healthy subjects (HCs) (Mann–Whitney U test). (DOC 119 kb)
Scatterplot showing loading values (coefficients) for the 77 immunological variables selected by the PLS model versus p values obtained when comparing levels of these variables between EC cases and EC controls with the Mann–Whitney U test. The dotted line marks the threshold for statistical significance (p < 0.05). (DOC 67 kb)
Immunological variables (in CD4 and CD8 T-cell subsets) with significant differences between EC cases and EC controls (Mann–Whitney U test). (DOC 52 kb)
Clinical centers and research groups that contribute to ECRIS. (DOC 33 kb)
About this article
Cite this article
Benito, J.M., Ortiz, M.C., León, A. et al. Class-modeling analysis reveals T-cell homeostasis disturbances involved in loss of immune control in elite controllers. BMC Med 16, 30 (2018). https://doi.org/10.1186/s12916-018-1026-6