A mechanistic model for long-term immunological outcomes in South African HIV-infected children and adults receiving ART
Abstract
Long-term effects of the growing population of HIV-treated people in Southern Africa on individuals and the public health sector at large are not yet understood. This study proposes a novel ‘ratio’ model that relates CD4+ T-cell counts of HIV-infected individuals to the CD4+ count reference values from healthy populations. We use mixed-effects regression to fit the model to data from 1616 children (median age 4.3 years at ART initiation) and 14,542 adults (median age 36 years at ART initiation). We found that the scaled carrying capacity, maximum CD4+ count relative to an HIV-negative individual of similar age, and baseline scaled CD4+ counts were closer to healthy values in children than in adults. Post-ART initiation, CD4+ growth rate was inversely correlated with baseline CD4+ T-cell counts, and consequently higher in adults than children. Our results highlight the impacts of age on dynamics of the immune system of healthy and HIV-infected individuals.
eLife digest
The human immunodeficiency virus (HIV) remains an ongoing global pandemic. There is currently no cure for HIV, but antiretroviral therapies can keep the virus in check and allow individuals with HIV to live longer, healthier lives. These drugs work in two ways. They block the ability of the virus to multiply and they allow numbers of an important type of infection-fighting cell called CD4+ T cells to rebound.
As more patients with HIV survive and transition from one life stage to the next, it is critical to understand how long-term antiretroviral therapies will affect normal age-related changes in their immune systems. The health of an immune system can be evaluated by looking at the number of CD4+ T cells an individual has, though this will vary by age and location. Clinicians use the same metrics to assess the immune health of individuals with HIV, however, as they age, it becomes a challenge to identify if a patient’s immune system recovers normally or insufficiently. Thus, learning more about age-related differences in CD4+ T cells in people living with HIV may help improve their care.
Using data from 1,616 children and 14,542 adults from South Africa, Ujeneza et al. created a simple mathematical model that can compare the immune system of person with HIV with the immune system of a similarly aged healthy individual. The model shows that among individuals with HIV receiving antiretroviral therapies, children have CD4+ T-cell numbers that are closest to the numbers seen in healthy individuals of the same age. This suggests that children may be more able to recover immune system function than adults after beginning treatment. Children also start antiretroviral therapies before their immune system has been severely damaged, while adults tend to start treatment much later when they have fewer CD4+ T cells left.
Ujeneza et al. show that the fewer CD4+ T cells a person has when they start treatment, the faster the number of these cells grows after starting treatment. This suggests that the more damaged the immune system is, the harder it works to recover. This reinforces the need to identify people infected with HIV as soon as possible through testing and to begin treatment promptly. The new model may help clinicians and policy makers develop screening and treatment protocols tailored to the specific needs of children and adults living with HIV.
Introduction
The efficacy with which antiretroviral therapy (ART) suppresses HIV viral load and restores lymphocyte responses to pathogen-derived antigens, normally lost due to viral replication, is well established (Cooney, 2002; Autran et al., 1997). Following ART initiation, for the majority of HIV-infected individuals, the CD4+ T-cell counts rapidly increase for approximately 4–6 months (Lawn et al., 2006), followed by a slower increase in the next 2–4 years, after which cell numbers plateau (Pinzone et al., 2012). Despite the successes of ART, several studies have demonstrated that large proportions of patients experience incomplete immune restoration following treatment initiation. In Sub-Saharan Africa, approximately 10–16% of children on ART (Mutwa et al., 2014) and variable proportions of adults (Barth et al., 2011; Fatti et al., 2014; Mee et al., 2008) fail to suppress viral load within 12 months of ART initiation. Further, among those that do suppress viral load, many are suboptimal immunological responders (35–40%), that is, they do not reach a CD4+ T-cell count greater than 500 cells/µL within 5 years (Nakanjako, 2016; Swiss HIV Cohort Study et al., 2005). Some patients demonstrate no improvement at all (Lawn et al., 2006).
A variety of statistical methods have been used to model the recovery of CD4+ T-cell counts after ART initiation. A prior review of sub-Saharan African studies of this type established that generalized linear mixed models and generalized estimating equations were most commonly used (Sempa et al., 2017). However, such models are not ‘mechanistic’, in the sense that they make no assumptions regarding the underlying biological processes involved in CD4+ T-cell reconstitution. For this reason, these methods do not allow inferences regarding the dynamics of CD4+ recovery. Recent studies have employed an asymptotic ‘semi-mechanistic’ mixed model to describe CD4+ count recovery in children (Lewis et al., 2012; ARROW Trial Team et al., 2013; Lewis et al., 2017; De Beaudrap et al., 2008) and adults (ANRS 1215/90 Study Group et al., 2009; Means et al., 2016). These models simply assume that CD4+ T-cell counts follow an asymptotic recovery process in all patients. However, it has been demonstrated that in patients with suboptimal or no immune recovery, response trajectories do not follow this profile (ARROW Trial Team et al., 2013). Complicating this picture, diverse CD4+ T-cell count variable transformations were applied, different demographic groups were studied and there is large natural variability of such cell counts across different population groups (Schaberg et al., 1997; Abuye et al., 2005; Maini et al., 1996). As a result, these models do not allow for comparisons across heterogeneous groups of individuals.
In this study, we present a novel mechanistic model that describes the dynamics of CD4+ T-cell count responses in HIV-infected individuals relative to those of healthy individuals. We use nonlinear mixed modelling methods and data from large cohorts of adults and children in South Africa. We assess the generalizability of this model from children to adults. We compare it to the previously published asymptotic model and investigate the impact of demographic and clinical characteristics on long-term immune outcomes.
Methodology
ART patient data
De-identified longitudinal data from HIV-infected patients receiving ART in 13 South African cohorts was provided by the International Epidemiologic Databases to Evaluate AIDS in Southern Africa (IeDEA-SA) (https://www.iedea-sa.org/) (ART-LINC Collaboration of IeDEA et al., 2008; IeDEA and COHERE Cohort Collaborations and Anderegg, 2018). The initial data set was composed of 2,858,743 CD4+ count observations for 223,688 patients of whom 202,108 were adults, 21,267 were children and 313 had no recorded date of birth or date of ART initiation. Children were defined as those aged 17 or younger at treatment initiation. After removing observations that were unrealistic (outside plausible biological ranges) or missing date of measurement, 189,647 adults and 19,060 children remained (Figure 1).
‘Ratio’ model construction
The ratio model relates the CD4+ T-cell dynamics from HIV+ individuals to those of healthy individuals. We first defined, then merged, separate models for each group, using the logistic growth model. This model assumes that a given quantity grows exponentially until it approaches a constant carrying capacity k (Figure 2), or limit, to represent the dynamics of the immune system. This concept has been used in prior studies to describe the proliferation rate of CD4+ T-cells for individuals on therapy (Kaufmann et al., 2001; Di Mascio et al., 2006). We denote the CD4+ T-cell count ‘growth’ (or regeneration) rate of an individual on ART by , and the ‘environmental’ (or physiological) carrying capacity by , which can be thought of as the maximal number of CD4+ T-cells that can be sustained by all available biological resources.
For HIV-infected people, the rate of change of CD4+ T-cells per per unit time, is expressed by the following ordinary differential equation:
where starts at ART initiation. The theoretical growth rate is actually an effective rate, a constant obtained by taking the average of the growth rate at different time points. This ‘instant’ growth rate increases for values of smaller than the inflexion point and starts decreasing thereafter, as the value of approaches the carrying capacity .
The solution to Equation 1 gives the expression for CD4+ T-cell counts per at time as:
where is the CD4+ T-cell count of the individual, at ART initiation.
A study on healthy South African population showed that their CD4+ counts increased from young adulthood until about 70 years of age (Malaza et al., 2013). Thus, we can also describe the CD4+ T-cell counts in a healthy individual at time , as:
where is the CD4+ T-cell growth rate, the carrying capacity, and time has an origin equivalent to ART initiation in an age-matched HIV infected patient. Thus, the equation describes the immune system dynamics of the healthy individual, starting from when they have a similar age as their corresponding HIV-infected counterpart.
The ‘Ratio’ model is then defined by dividing Equation 2 by Equation 3, substituting , and , and simplifying to obtain, , the scaled CD4+ T-cell count at time :
where and are the CD4+ T-cell scaled carrying capacities of an HIV-infected and a healthy age-matched individual, respectively; and the ratio of the baseline CD4+ T-cell count of the HIV-infected individuals to that of the healthy individual, in order words, the baseline scaled CD4+ T-cell counts. An implication of this model is that when the value of is close to unity, the patient’s long-term carrying capacity is approaching the particular patient’s ‘homeostatic optimum’ for CD4+ T-cell counts, that is, the amount of CD4 cells that can be sustained by the available biological resources. In this case a subsequent increase in CD4+ T-cell counts would be unlikely. Further, associated with an ‘immune set-point’ achieved in healthy adults their value of would tend to unity. In healthy children, this is not so owing to growth, that is, their CD4+ T-cell count changes with the increase of their total blood volume (Bains et al., 2009). Thus, their value of is expected to be less than one.
Additionally, a baseline scaled CD4+ T-cell count of 1 means that the individual started therapy with a normal CD4+ T-cell count for their age, while a value smaller than one indicates that an individual started therapy with a lower CD4+ count compared with those of healthy individuals of the same age. Note that both growth rates for HIV-infected and healthy individuals are the same (unchanged), from the original logistic growth models to the ‘ratio’ model. A table indicating the ranges for each parameter is available in Supplementary file 1.
Model fitting techniques
We used nonlinear mixed effects techniques (NLMM) to fit all models to data. NLMMs are characterized by two main components: a fixed effect part that describes the population mean, and the random effects that describe individual subjects' deviations from the mean. If we denote by the vector of population estimates, and by the matrix of the random effects’ estimates, a general mathematical description of NLMM for continuous variables is given by,
where is the scaled CD4+ T-cell counts of patient at time function is the nonlinear model of interest and is the structure of the error model; is the total number of subjects and the number of observations for each subject . Both and depend on the predictor variables , the population parameters , and the random effects . In addition, the function depends on the form chosen for the error model, which is governed by a set of parameters In this study, we used a proportional error model, such that and (single constant parameter ), and Equation 4 as our basis function :
For purposes of comparison, we also applied the same NLMM methodology to an ‘asymptotic’ model, similar to that previously described (Lewis et al., 2012), such that,
where is the value of the asymptote for patient , their intercept, and their logarithm of rate of increase of scaled CD4+ T-cell counts. The term represents the scaled increase of CD4+ T-cells following ART initiation.
The NLMM method assumes that the random effects are normally distributed, with mean zero and variance-covariance matrix :; and the errors are also normally distributed with mean zero and variance 1: . These represent noise and errors in the data. The individual parameters (denoted by a matrix ) are given by , with the matrix of individual covariates. A demo code is available through a github repository (https://github.com/EvaLiliane/RM_Code_eLife [copy archived at swh:1:rev:624ff31c5fc969885f29b7291ee06886d24c64f7]; Ujeneza, 2020).
Variable scaling
The key outcome variable is scaled CD4+ T-cell count. In both infected adults and children, the cell counts post-ART initiation were scaled by reference values from healthy populations, to obtain the outcome variable.
For HIV-infected children these reference values were calculated from the cross-sectional data (see description in Appendix 1) of healthy children at specific ages, due to the large variability in CD4+ T-cell counts in the early years of life. For the reference values by age, a single exponential model was fitted to the healthy children’s cross-sectional data and continuous population estimates were simulated (see Appendix 1—table 1 and Appendix 1—Figure 1). These were within the normal CD4+ T-cell counts ranges published in South Africa (Lawrie et al., 2009; Lawrie et al., 2015) and elsewhere (Idigbe et al., 2010; Pediatric AIDS Clinical Trials Group et al., 2003). We then scaled all HIV-infected children’s CD4+ T-cell counts as follows:
where is the scaled CD4+ T-cell counts of patient of age (in months) at time (measured as time since ART initiation); is the CD4+ T-cell counts of an HIV-infected child of age and is the CD4+ T-cell counts of a healthy child of similar age as patient . Scaling CD4+ T-cell counts of HIV-infected children by that of healthy children of similar ages enabled the comparison of CD4+ T-cell counts responses across ages, while simultaneously accounting for the child’s growth and immune system development.
For adults, normal reference values were estimated using values obtained from the literature (Malaza et al., 2013; Lawrie et al., 2009, Ngowi et al., 2009; Lugada et al., 2004; Institute of Human Virology/Plateau State Specialist Hospital AIDS Prevention in Nigeria Study Team et al., 2005; Oladepo et al., 2009) (see details in Appendix 1—table 2 and Appendix 1—figure 2). Although healthy adult CD4+ T-cell counts are known to vary with age, this variability is small compared to that observed in children, with an average difference of 100 cells between 25 and 60 years old healthy adults [Oladepo et al., 2009] vs 1500 cells difference between 3 months and 15 years old healthy children [Lawrie et al., 2015]. Similarly, average CD4+ T-cell counts differences between male and female adults in South Africa were in the range of 80–150 cells (Malaza et al., 2013; Lawrie et al., 2009). We evaluated results obtained from scaling HIV-infected CD4+ counts by simulated age and sex-dependent reference CD4+ counts values, and found minor or no difference with those obtained from CD4+ T-cell counts scaled by a single value. We obtained exactly the same population estimates and very similar individual parameter estimates. These were evaluated and we found no significant differences (see Supplementary file 2). Thus, due to the unavailability of individual age-specific data for a South African healthy adult population, a single normal reference value, that is, a constant of 800 CD4+ T-cells per µL, was used to scale all CD4+ T-cell count observations for adults on ART ().
Baseline scaled CD4+ T-cell counts
Our analysis was restricted to patients with available baseline (i.e. at antiretroviral treatment initiation) CD4+ T-cell counts and sufficient observations to estimate all model parameters (Figure 1). We categorized as baseline any CD4+ count measurement that was taken within 15 days before or after the ART initiation date. Two model-fitting scenarios were defined based on the availability of enough observations to estimate the model parameters:
Scenario 1: In which baseline scaled CD4+ T-cell counts were estimated, including 1312 children and 12,238 adults, with a minimum of five CD4+ T-cell counts measurements and no missing values for our variables of interest (specified below).
Scenario 2: Where baseline scaled CD4+ T-cell counts were used as a predictor, including 1616 children and 14,542 adults, with a minimum of four CD4+ T-cell counts measurements and no missing values for our variables of interest. These variables were CD4+ T-cell counts since ART initiation; viral load, age and body mass index at ART initiation; sex and suppression (or not) of viral load within 12 months of ART initiation.
In this paper we will mostly discuss results of scenario 1, where all parameters are estimated. The methodology used and most of the obtained results are applicable for scenario 2, when baseline CD4+ T-cells counts data are available. We also present results of scenario 2.
Structure of the random effects
Under both scenarios, in the initial model fitting, all parameters were assumed to vary per individual (i.e. there were random effects on all parameters) and their distributions was set as log-normal. Thus, for each parameter , the distribution of individual values was , with mean μ and variance ω2. Random effects were assumed to be independent and the identity diagonal matrix was used for the variance–covariance structure.
This assumption was later relaxed and different variance–covariance structures were subsequently evaluated. To compare models, the Akaike and Schwarz information criteria (AIC and BIC respectively), computed by importance sampling, were used (see Comets et al., 2017). Models with a full variance–covariance matrix were retained as this structure gave the lowest AIC and BIC values.
Covariates
Our choice of covariates was based on literature reviews (Pinzone et al., 2012; Sempa et al., 2017), data availability, and biological plausibility. We considered baseline characteristics that were measured within 15 days before or after ART initiation, namely age, z-score body-mass-index (BMI, for children), and viral load. Other covariates included sex and viral suppression (or not) within 12 months of ART initiation.
Patient age at ART initiation was expressed in months. Z-score BMI at baseline for children was calculated using WHO-Igrowup’s package (WHO Multicentre Growth Reference Study Group,, 2006). ‘Suppress’, a binary, was defined as reaching an undetectable viral load (<1000 copies/mL) or not, within 12 months of ART initiation. Baseline viral load was log transformed to simplify model fitting. Baseline age and sex effects were included for all parameters, while viral load suppression, baseline log viral load and BMI z-scores for children, were only included for the parameters describing the longitudinal immune responses of individuals on ART, namely: scaled carrying capacity post-ART, CD4+ T-cell growth rate post-ART and baseline scaled CD4+ T-cell count .
Additional considerations
Final adjusted models were obtained using a backward–forward stepwise approach using p-value criterion for covariate inclusion and AIC, BIC for the overall models. To evaluate the robustness of our results, we compared the estimated parameters of the full analysis described above using only three covariates: sex, baseline age, and baseline viral load, with those from models with all covariates.
All graphs and computations were produced in the statistical environment R (R Development Core Team, 2017). Nonlinear mixed model fitting employed the saemix R package (Comets et al., 2005), which uses a stochastic approximation, expectation-maximization algorithm for parameter estimation. All final adult and children models (Figures 3 and 4) converged relatively rapidly towards their estimated values, as shown by the log-likelihood graphs (Figure 3—figure supplement 1 and Figure 4—figure supplement 1). Even when initial values were marginally varied, the final models converged to similar estimates.
Results
Clinical characteristics
Cohort characteristics are summarized in Table 1. For scenario 1, the median number of clinical visits was 8 per patient, in both adults IQR (6,9) and children IQR (6,10), with a median follow-up time of 3 years for adults IQR (2,5) and 4 years for children IQR (3,5). At baseline, that is, ART initiation, the median age for children was 4.5 years IQR (1.4, 7.9) and 36 years IQR (30.7, 43) for adults. A clinical WHO stage III was most common in all patients at ART initiation. Among adults, there were more females than males (56.7% versus 43.2% respectively), but only slightly more females among children (51.3% versus 48.6% males). All patients in our data sets initiated therapy between 1997 and 2013, with 95% of them initiating in the period, 2003–2012. About 53.4% of children reached a CD4+ T-cell count of 500 per µL within 12 months of ART initiation, while only 4.4% of adults reached the same threshold. This number increases to 10.1% after 2 years of therapy for adults.
The median baseline CD4+ T-cell counts for children and adults were low at 404 per µL IQR (159.7,706.2) and 128 per µL IQR (62,196) respectively (Table 1). Higher counts at ART initiation were more common in younger (0–5 year’s age group) versus older children (>5 years). In adults, no differences were found in CD4+ T-cell counts at baseline, for different age groups. Both children and adults that presented with higher counts at baseline had lower viral loads versus those with low baseline counts.
Comparison of models and random effects structures
The performance of both asymptotic and ratio models, in terms of AIC and BIC criteria, was found to be strongly dependent on the structure selected for the variance–covariance matrices (see Supplementary file 3). Initially, both models were defined with random effects on all parameters and a diagonal matrix structure for the random effects, that is, individual random effects were assumed independent. Different matrix structures were then compared and those with full matrices, that is, in which random effects are correlated, were found to be the best.
For both adults and children, the AIC and BIC for the Ratio model were smaller than for the Asymptotic model when baseline scaled CD4+ T-cell counts were estimated (Scenario 1). When the baseline scaled CD4+ T-cell counts were used as a predictor (Scenario 2) this reversed. However, the number of parameters estimated by the two models was different. For the Asymptotic model there were three for scenario one vs two for scenario 2, a 33% change, and for the ratio model there five for scenario one vs four for scenario 2, a 20% change (Table 2). Given the different numbers of parameters, the Ratio performed well in comparison to the Asymptotic model.
In terms of parameter estimations, both models estimated baseline scaled CD4+ of similar magnitudes, for HIV-infected adults and children (Table 3 and Supplementary file 4 – Tables 1 and 2). The asymptotic model predicts a sharp increase of CD4+ counts, which rapidly reaches an asymptote, while the ratio model predicts a sharp increase that is followed by a slower but still gradual increase of scaled CD4+, up to more than 10 years after ART initiation (Figures 3 and 4).
We note that the parameter estimates obtained using a baseline CD4 scaling constant of 800 cells/µL for healthy adults were comparable to those obtained using age-dependent healthy adults CD4 values (Supplementary file 2). Similarly, the estimates obtained using a baseline CD4 scaling constant of 800 cells/µl for healthy adults were also comparable to those obtained without scaling (i.e. using a scaling constant of 1 cell/µL) for all parameters except z0, whose estimate increased by about 800-fold (Supplementary file 2) as expected.
Final ratio models
We found no major differences in the estimated fixed effects when adjusting the models for age and sex only, versus adjusting for all our covariates. The minor differences between both fits were the sizes of estimated parameters (Supplementary file 4 – Table 2), a few inconsistences in the effects of age on particular parameters (adults and children) and sex effects on post-ART scaled carrying capacity (adults only). Therefore, we only describe results for the models adjusted for all covariates. These results are presented in Table 3. Note that a sample of individual fits is shown in Figure 3—figure supplement 2 and Figure 4—figure supplement 2.
Population estimates/fixed effects
Both adults and children, on average, start ART with an impaired immune system. We found that baseline scaled CD4+ T-cell count was 0.20 (95% CI: 0.18, 0.22) for adults and 0.74 (95% CI: 0.49, 1.00) for children (Scenario 1, Table 3).
The growth rates of CD4+ T-cell counts were more than tenfold higher in HIV infected versus healthy individuals. We found that in adults, the growth rate for HIV-infected individuals was 0.49 cell counts per µL per day, while it was 0.02 for healthy individuals. In infected children, a growth rate of 0.27 cell counts per µL per day was found, while it was much lower (0.01) in healthy children. This means, for example, that an infected child starting ART with 100 cells will cross the 500 cells threshold in 1470 days, whereas an infected adult with the same number of cells at ART initiation will cross the threshold much faster, after about 816 days.
For HIV-infected individuals the scaled carrying capacity was higher than that of healthy individuals in both adults and children. We found that the scaled carrying capacity for HIV-infected adults was 1.75, which is three times higher than that of healthy adults (0.49). We found that the scaled carrying capacity for HIV-infected children was above one (1.03). The scaled carrying capacity post-ART of 1.75 in adults implied that an average 36.5 year old adult on ART would have to experience a 75% increase of their initial CD4+ T-cells to reach their ultimate CD4+ T-cell counts levels. For an average 5-year-old child, this rebound needs to be much lower: 3% increase of their initial CD4+ T-cell counts, to reach normal ranges.
Individual estimates and random effects
For both adults and children, the coefficient of variation for the random effects for all parameter estimates was lower for the parameters of HIV-infected patients (k, r, and z0) and higher for the healthy individual parameters (q and s). This is understandable given the lack of individual matching in the data, that is, a person could not be HIV-infected and healthy simultaneously.
As can be seen in Table 4, the correlation between the random effects was similar in adults and children, with the exception of the correlation between the scaled carrying capacity and growth rate of CD4+ T-cells for healthy individuals, that was positive and strong (0.61) in adults, but weak and negative (−0.13) in children. The strongest correlation was between the scaled carrying capacity for HIV-infected individuals and baseline scaled CD4+ T-cell counts, in both adults (−0.67) and children (−0.85).
Covariate estimates in children
We found that being 1 month older resulted in a decrease in baseline scaled CD4+ T-cell count of −0.007 (95% CI: −0.009, −0.006). Similarly, being 1 month older was associated with a higher CD4+ T-cell growth rate in healthy children (0.008) and a significant decreasing effect on post-ART CD4+ T-cell growth rate (−0.001). This means that, given a threshold of 200 CD4+ cell counts, younger children will reach it earlier than older children that had a similar CD4+ T-cell count at baseline. Furthermore, we found that age had a decreasing effect on the scaled carrying capacity of healthy children (−0.005).
Both scenarios show similar effects of baseline viral load on children’s parameters. We found that one unit increase in baseline log viral load was associated with a −0.07 decrease in the baseline scaled CD4+ T-cell count, and a +0.06 increase in the CD4+ T-cell scaled carrying capacity post-ART. It also led to a +0.03 increase in the growth rate of CD4+ T-cells post-ART which indicates that those with higher viral loads at ART initiation will reach a defined threshold faster than those with lower viral loads, if they started with the same initial number of cells.
Baseline BMI z-score was associated with children’s baseline scaled CD4 T-cell count, post-ART CD4+ T-cell scaled carrying capacity, and growth rate. One unit increase in baseline BMI z-score resulted in an increase of +0.15 in baseline scaled CD4+ T-cell count, a decrease in the growth rate of CD4+ T-cells post-ART (−0.08) and in the post-ART scaled carrying capacity (−0.12). An untreated child with advanced HIV-disease is likely to have a low BMI z-score and will consequently have a low scaled CD4+ T-cell count at ART initiation.
Covariate estimates in adults
In adults, being 1 month older was associated with an increase of +0.0006 (95% CI: 0.0004, 0.0008) in baseline scaled CD4+ T-cell counts. We found no effect of age on post-ART CD4+ T-cell growth rate, but it was associated with a decline in CD4+ T-cell growth rate in healthy adults (−0.0012 per month older). Our results show that being older was associated with a decrease in the scaled carrying capacity in adults on ART (−0.0007 per month older), though it had no effect on the scaled carrying capacity of healthy adults.
Sex and baseline log viral load had an effect on CD4+ T-cell dynamics in adults. We found that adult females started ART therapy at higher scaled CD4+ T-cell values compared to males (+ 0.08) and one unit increase in baseline log viral load was associated with a decrease in the baseline scaled CD4+ T-cell counts (−0.08). We also observed a +0.07 increase in the CD4+ T-cell scaled carrying capacity post-ART and a +0.11 increase in the growth rate of CD4+ T-cells post-ART, per unit of increase in baseline log viral load.
Discussion
CD4+ T-cell count remains the main proxy used to evaluate the long-term effects of ART on the immune system. Understanding its dynamics is crucial, to ensure that the best care is delivered to HIV-infected patients. Although previous studies described CD4+ T-cell count dynamics in HIV-treated patients, using non-mechanistic (Sempa et al., 2017) and mechanistic models (Lewis et al., 2012; Means et al., 2016; Di Mascio et al., 2006), none provided a method that allows for direct comparison with age-matched healthy controls. This study proposes a mechanistic model for the immune system reconstitution, which relates CD4+ T-cells count of an individual on ART with that of a healthy individual with similar characteristics. We demonstrate that CD4+ T-cell growth rates are higher in HIV-treated patients than in healthy individuals and that age has opposite effects on CD4+ counts dynamics in HIV-treated children, compared to healthy children. This model was compared to the asymptotic model, previously used to model CD4+ T-cell trajectories after ART initiation, under two different scenarios: estimating baseline scaled CD4+ T-cell count (scenario 1) and using it as a predictor (scenario 2). Moreover, this study is the first of its kind in evaluating large samples of children and adults data in a comparative way.
Corbeau and Reynes, 2011 have described the three phases of CD4+ counts recovery in HIV-treated patients: a sharp increase in the first 1–6 months, a still high increase up to 2 years, and a slow gradual increase that goes on beyond 4 years of therapy. Our results show that though both the ratio and asymptotic model predict similar scaled baseline CD4+ counts, only the ratio model is able to reproduce these three phases of immune reconstitution. Thus, it is likely to be better in accurately predicting the long-term behaviour of CD4+ T-cell counts trajectory. Furthermore, the ratio model is derived directly from the logistic growth model which describes a known biological process, and as such, CD4+ T-cells growth rate is directly relatable to the average of the actual increase of CD4+ T-cells per month over the observation period.
Baseline scaled CD4+ T-cell counts (z0)
Fitting our model to the data demonstrated that larger values of baseline scaled CD4+ T-cell counts are associated with larger values of scaled CD4+ T-cell counts at any subsequent time. This is the consequence of equation 4, which is in line with the findings of prior studies in which higher baseline CD4+ T-cell counts were associated with a higher final or plateauing value for CD4+ T-cell counts (Moore and Keruly, 2007). Similarly, patients with very low baseline CD4+ T-cell counts, <350 cells/µL, often fail to reach normal values even after long durations of therapy (Nakanjako, 2016; Swiss HIV Cohort Study et al., 2005; Moore and Keruly, 2007; Kelley et al., 2009). We also found that the scaled carrying capacity of CD4+ T-cells in HIV-infected individuals was negatively correlated with the baseline scaled CD4+ T-cell counts in both adults and children, respectively, −0.64 and −0.81. This seems reasonable in that the closer an individual is to their normal or optimal CD4+ T-cell value at ART initiation, the less cell population expansion is required to reach normal levels.
Scaled carrying capacity of CD4+ T-cells for HIV-infected and healthy individuals
Our parameter estimates demonstrated that the scaled cellular carrying capacity was higher in individuals on ART than in those who were healthy. In both HIV-infected adults and children the scaled carrying capacity post-ART was greater than 1, meaning that baseline CD4+ T-cell counts were lower than their corresponding long-term homeostatic optimum. That is, individual CD4+ T-cell counts had to grow to reach normal levels. This is consistent with the notion that an impaired immune system usually experiences repair following treatment initiation.
In contrast, the scaled capacity for healthy adults and children were both lower than unity. In children the value was 0.68, which is consistent with a mechanistic understanding of the dynamics of healthy CD4+ T-cell counts in early life. In particular, there is an increase until the age of approximately 1 year and then a decrease thereafter (Bains et al., 2009). Total blood and body volume increase throughout childhood associated with the shrinkage of the thymus, which is accompanied by a reduction of naïve T-cell production (Hapuarachchi et al., 2013; Hazenberg et al., 2004). In combination, this results in decreasing CD4+ T-cell counts per volume with age. In healthy adults, we found a scaled carrying capacity of 0.49, suggesting a decreasing CD4+ T-cell trend with age. This is in agreement with other studies from elsewhere which found that CD4+ T-cell counts decreased from young adulthood to middle age (Lugada et al., 2004; Zeh et al., 2011). However, one study in healthy South African adults has described CD4+ T-cell counts increasing slightly until the age of 64 years old (Malaza et al., 2013).
Post-ART CD4+ T-cell count growth rate depends on scaled baseline CD4+ T-cell counts
We found that post-ART CD4+ T-cell growth rate depends inversely on the cell count at treatment initiation, that is, they are negatively correlated (Table 4). This implies that the higher the scaled CD4+ T-cell counts of an individual at ART initiation, that is, the closer it is to its healthy/normal true value, the lower the rate of recovery will be. This is reasonable given a decrease in the need to achieve normal levels and is in agreement with the findings of prior studies (Lawn et al., 2006; Sachsenberg et al., 1998).
Post-ART CD4+ T-cell growth rates in HIV-infected versus healthy individuals
We found that CD4+ T-cell growth rates were higher in individuals on ART than in healthy individuals. This agrees with prior studies in which cell growth rates in adults on ART were estimated to be sixfold to tenfold greater than in healthy adults (Sachsenberg et al., 1998; Hazenberg et al., 2000). This is understandable bearing in mind the immune system’s effort after treatment initiation to re-fill the void of peripheral CD4+ T-cells destroyed during HIV infection. Biological studies have described this as an initial redistribution of memory T-cells from the lymph nodes into the blood stream, followed by homeostatic proliferation and production of naïve cells by the thymus (Tsukamoto et al., 2009; Autran, 1999). Thus, given that HIV-infected individuals have much lower CD4+ T-cell counts they require higher CD4+ T-cell growth rates than healthy individuals. This behaviour was captured by our model.
Interestingly, the differences in our estimates for CD4+ T-cell growth rate in healthy individuals versus patients on ART were slightly lower than those found elsewhere (Sachsenberg et al., 1998). This may be due to: the differences in the total time of follow-up of our patients, selection bias in the original population data, differences in the demographics of the populations studied, the types of treatments administered, and changes in the WHO guidelines for ART initiation over time. The WHO guidelines for minimal CD4+ T-cell counts at initiation changed from 200, to 350, to 500 and later to initiation at diagnosis. Thus, post-ART CD4+ T-cell growth rate in HIV-infected individual might be smaller now than it may have been on average in the past, as CD4+ T-cell counts at ART initiation are now higher.
The effects of age
Baseline scaled CD4+ T-cell count was lower in adults compared to children (0.2 vs 0.8), that is, adults started treatment when their immune systems were more compromised compared to children. This could be due to the fact that children (<17 years old) are more likely vertically infected and, thus, they are more likely to be diagnosed early, owing to early testing and follow-up in the South African program for prevention of mother-to-child transmission. Adults, on the other hand, might be infected for an extended period and are consequently more likely to be more highly immune-compromised prior to diagnosis, compared to children.
We found that older age is associated with a lower value of baseline scaled CD4+ T-cell count in children and a higher value in adults. Vertical HIV- infection and an extended duration without treatment may lead to greater immune compromise in older children. In contrast, higher baseline scaled CD4+ T-cell counts in older adults suggest that they may have been more health conscious than younger adults (Prohaska et al., 1985), that is, younger adults tend to seek treatment later than older adults.
Such findings might also be explained by the negative correlation between post-ART CD4+ T-cell growth rates and baseline CD4+ T-cell counts: younger adults have a higher growth rate as they start at a lower CD4+ T-cell count value (Cornell et al., 2012). Thus, their CD4+ T-cell rebound value may be higher than that of older adults over a similar period of time (Means et al., 2016; International epidemiological Database to Evaluate AIDS (IeDEA) West Africa Collaboration et al., 2012). Our results make no inferences regarding the period an individual takes to reach their normal CD4+ T-cell count value.
Sex effects
We found that adult females have lower post-ART scaled carrying capacities than men, meaning they need a lower CD4+ T-cell count rebound to reach normality. This suggests that female adults tend to initiate treatment earlier than men (Cornell and Myer, 2013). This was also validated by the fact that baseline scaled CD4+ T-cell counts were higher in females than in males. A prior study has also found that women tend to have higher plateauing CD4+ T-cell counts than men (Means et al., 2016). In our study, women had lower rates of recovery than males which also supports our finding that higher baseline scaled CD4+ T-cell count is associated with a lower CD4+ T-cell growth rate. As they have higher CD4+ T-cell counts at ART initiation, and given the inverse correlation between the post-ART CD4+ T-cell growth rate and baseline scaled CD4+ T-cell count, women demonstrate lower growth rates.
Previous studies have found that adult males and females spend a similar time on therapy before reaching their rebound set-point CD4+ T-cell count (EuroSIDA group et al., 2003; Patterson et al., 2007). We do not consider this a contradiction of our results as an individual with a higher baseline CD4+ T-cell count will have a lower cell growth rate compared to an individual that started ART with a low CD4+ count. Thus, after an equivalent period following ART initiation similar increases in CD4+ T-cell counts may have occurred. We believe that the diverging opinions regarding the effects of sex on immune outcomes in adults post-ART initiation might be explained by variations in the definitions of ‘immunological outcomes’ by different authors and variations in analyses conducted.
This study has strengths and limitations. To ensure that patients had sufficient data to enable parameterization of the model, only subsets of the full data set were used in the analysis. This prevented over-fitting, but it might have introduced selection bias. However, with the exception of the percentage of people that suppressed viral load within 12 months of starting ART (Table 1), comparison of the summary statistics of the subset versus the full data set demonstrated that they were similar. Further, the results obtained from adjusting and not adjusting the model parameters with the variable, viral load suppression within 12 months, were also in agreement (see Supplementary file 4). We consequently believe that the subset used was representative of the full data set. Our model does not distinguish between naïve and memory sub-types of CD4+ T-cells. Prior studies have shown that these have different dynamics (Di Mascio et al., 2006; Bajaria et al., 2002). However, in routine HIV monitoring, for which we had data, the different subtypes are not measured. As ART treatment is for life, separating analyses for children from that for adults does not account for children growing into adulthood. Future studies might express the scaled carrying capacity as a function of age rather than as a variable parameter. Ignoring the thymus’s contribution to CD4+ T-cell recovery, particularly in the earlier years of life, might have resulted in an overestimation of and parameters which are aggregate growth rates. Lastly, scaling the CD4+ T-cell counts of adults by a single average healthy CD4+ count value for adults might have introduced bias in the estimates of the scaled carrying capacity of healthy adults, due to variations of CD4+ T-cell counts across different age and demographic groups.
Our study does provide insight into the effect of ageing on immune system dynamics in adults and children on ART compared to healthy individuals. The ratio model provides a more accurate estimation of CD4+ counts reconstitution than the asymptotic model as well as the ability to compare different immune system outcomes, for both healthy and HIV-treated individuals. Using scaled CD4+ T-cell counts allows for the evaluation of CD4+ counts trajectories, which is not possible with unscaled CD4+ T-cell counts. We found large variations in CD4+ T-cell growth rates and scaled carrying capacities between individuals, highlighting the need to evaluate ART outcomes on an individual level. This calls for improved patient monitoring strategies. The strong inverse correlation between baseline scaled CD4+ T-cell count and the scaled carrying capacity emphasizes the importance of early ART initiation, regardless of age or state of disease progression. We found that post-ART CD4+ T-cell growth rate is not associated with a patient’s age, but it is associated with higher baseline viral load. With the expansion of an aging population on ART, understanding long-term effects of the treatment on the immune system is critical to ensure that the best care is delivered to HIV-infected patients.
Appendix 1
A1: Simulated references CD4+ T-cells values per age, for healthy children
Cross-sectional CD4+ T-cell counts data of 381 healthy children were collected as part of a study that sought to evaluate the use of international paediatric reference intervals by NHLS for South African populations (Lawrie et al., 2015). The study was conducted at a child wellness clinic, in Wesbank, a semi-informal settlement located in Cape Town. Healthiness was defined as not being on any prescribed medication, no chronic illness, a full clinical history, and attendance of the child to the clinic, with his/her biological mother. These children had a median age of 18.9 months IQR (6.5, 45.6) and a median CD4+ T-cell counts of 2007 cells/µL IQR (1433, 2749).
We simulated age-dependent normal CD4+ counts for healthy children, using a single exponential model as suggested in previous (unpublished) work. We model healthy CD4+ counts of a healthy child at age as follow:
,
where the constants and (unit: cells per µL) give an indication of the CD4+ count of a child at birth, and indicates the rate of decline of CD4+ counts (unit: cells per µL per month), thereafter. We used generalized nonlinear least squares method to estimate the model’s parameters. We found that , and were equal to 1070, 2.174, and 0.031 respectively. Using these values, we simulated CD4+ counts references values for each month of age, from birth until 17 years. As an indication, a sample of the simulated values are shown in the table below. The graph (see Appendix 1—figure 1) shows the raw data with the simulated references values as a red line.
-
Appendix 1—figure 1—source data 1
- https://cdn.elifesciences.org/articles/42390/elife-42390-app1-fig1-data1-v2.csv
A2: Age and sex-dependent CD4+ T-cell counts reference values for healthy adults
In the study, we scaled all adults CD4+ counts by a single constant value, matching the studies below.
We also fitted cubic splines to age-group published median CD4+ counts data (Malaza et al., 2013), for men and women separately. We then simulated age-dependent reference values as shown in the graph below (Appendix 1—figure 2). We found no difference between our initial results and those obtained from scaling data of HIV-infected adults by age and sex-dependent reference values.
-
Appendix 1—figure 2—source data 1
- https://cdn.elifesciences.org/articles/42390/elife-42390-app1-fig2-data1-v2.csv
Data availability
Data used are from the International Epidemiologic Databases to Evaluate AIDS Southern Africa collaboration. They maintain a database of routinely collected data from various clinics, mostly located in South Africa. We recommend that interested readers contact Dr Morna Cornell, Project Manager IeDEA-SA in Cape Town (morna.cornell@uct.ac.za) to establish a data-sharing agreement. A research proposal highlighting how the data will be used is required. Source data for figures and figure supplements are provided, and the source code is available at https://github.com/EvaLiliane/RM_Code_eLife copy archived at https://archive.softwareheritage.org/swh:1:rev:624ff31c5fc969885f29b7291ee06886d24c64f7/.
References
-
Modeling CD4+ cell count increase over a six-year period in HIV-1-infected patients on highly active antiretroviral therapy in SenegalThe American Journal of Tropical Medicine and Hygiene 80:1047–1053.https://doi.org/10.4269/ajtmh.2009.80.1047
-
Restoration of the immune system with anti-retroviral therapyImmunology Letters 66:207–211.https://doi.org/10.1016/S0165-2478(98)00159-X
-
Dynamics of naive and memory CD4+ T lymphocytes in HIV-1 disease progressionJAIDS Journal of Acquired Immune Deficiency Syndromes 30:41–58.https://doi.org/10.1097/00042560-200205010-00006
-
Long-Term outcome of an HIV-Treatment programme in rural africa: viral suppression despite early mortalityAIDS Research and Treatment 2011:1–5.https://doi.org/10.1155/2011/434375
-
SoftwareSAEMIX, an R Version of the SAEM Algorithm 20th MeetSAEMIX, an R Version of the SAEM Algorithm 20th Meet.
-
Parameter estimation in nonlinear mixed effect models using saemix, an R implementation of the SAEM algorithmJournal of Statistical Software 80:i03.https://doi.org/10.18637/jss.v080.i03
-
Clinical indicators of immune restoration following highly active antiretroviral therapyClinical Infectious Diseases 34:224–233.https://doi.org/10.1086/323898
-
CD4 cell response before and after HAART initiation according to viral load and growth indicators in HIV-1-infected children in Abidjan, côte d'ivoireJAIDS Journal of Acquired Immune Deficiency Syndromes 49:70–76.https://doi.org/10.1097/QAI.0b013e3181831847
-
Virologic, immunologic, and clinical response to highly active antiretroviral therapy: the gender issue revisitedJAIDS Journal of Acquired Immune Deficiency Syndromes 32:452–461.https://doi.org/10.1097/00126334-200304010-00017
-
T cell depletion in HIV-1 infection: how CD4+ T cells go out of stockNature Immunology 1:285–289.https://doi.org/10.1038/79724
-
T-Lymphocyte subsets in apparently healthy nigerian childrenInternational Journal of Pediatrics 2010:1–7.https://doi.org/10.1155/2010/474380
-
Global Trends in CD4 Cell Count at the Start of Antiretroviral Therapy: Collaborative Study of Treatment ProgramsClinical Infectious Diseases : An Official Publication of the Infectious Diseases Society of America 66:893–903.https://doi.org/10.1093/cid/cix915
-
Reference values of CD4 T lymphocytes in human immunodeficiency virus-negative adult nigeriansClinical Diagnostic Laboratory Immunology 12:525–530.https://doi.org/10.1128/CDLI.12.4.525-530.2005
-
Relative significance of different pathways of immune reconstitution in HIV type 1 infection as estimated by mathematical modelingAIDS Research and Human Retroviruses 17:147–159.https://doi.org/10.1089/08892220150217238
-
Incomplete peripheral CD4+ cell count restoration in HIV-infected patients receiving long-term antiretroviral treatmentClinical Infectious Diseases : An Official Publication of the Infectious Diseases Society of America 48:787–794.https://doi.org/10.1086/597093
-
Local reference ranges for full blood count and CD4 lymphocyte count testingSouth African Med. J 99:243–248.https://doi.org/10.7196/SAMJ.2745
-
Age and CD4 count at initiation of antiretroviral therapy in HIV-infected children: effects on long-term T-cell reconstitutionJournal of Infectious Diseases 205:548–556.https://doi.org/10.1093/infdis/jir787
-
Population-based hematologic and immunologic reference values for a healthy ugandan populationClinical Diagnostic Laboratory Immunology 11:29–34.https://doi.org/10.1128/CDLI.11.1.29-34.2004
-
Reference ranges and sources of variability of CD4 counts in HIV-seronegative women and menSexually Transmitted Infections 72:27–31.https://doi.org/10.1136/sti.72.1.27
-
Establishment of reference values of CD4 and CD8 lymphocyte subsets in healthy nigerian adultsClinical and Vaccine Immunology 16:1374–1377.https://doi.org/10.1128/CVI.00378-08
-
Lymphocyte subsets in healthy children from birth through 18 years of age: the pediatric AIDS clinical trials group P1009 studyThe Journal of Allergy and Clinical Immunology 112:973–980.https://doi.org/10.1016/j.jaci.2003.07.003
-
HIV RNA suppression and immune restoration: can we do better?Clinical and Developmental Immunology 2012:1–12.https://doi.org/10.1155/2012/515962
-
Health practices and illness cognition in young, middle aged, and elderly adultsJournal of Gerontology 40:569–578.https://doi.org/10.1093/geronj/40.5.569
-
BookR: A Language and Environment for Statistical ComputingVienna, Austria: R Foundation for Statistical Computing.
-
Turnover of CD4+ and CD8+ T lymphocytes in HIV-1 infection as measured by Ki-67 antigenJournal of Experimental Medicine 187:1295–1303.https://doi.org/10.1084/jem.187.8.1295
-
SoftwareGithub “RM_Code_eLife”, version swh:1:rev:624ff31c5fc969885f29b7291ee06886d24c64f7Software Heritage.
-
Book“WHO Child Growth Standards: Length/height-for-Age, Weight-for-Age, Weight-for-Length, Weight-for-Height and Body Mass Indexfor-Age: Methods and Development.”, 312Geneva, Switizerland: World Health Organization.
Article and author information
Author details
Funding
Schlumberger Stichting Fund, under the context of the Future Fellowship Program
- Eva Liliane Ujeneza
South African Department of Science and Technology and the National Research Foundation's Center of Excellence for Modelling and Analysis of Epidemiological Data
- Eva Liliane Ujeneza
National Institute Of Allergy And Infectious Diseases of the National Institutes of Health (U01AI069924)
- Mary-Ann Davies
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors gratefully acknowledge the support from Dr Joanna Lewis in the initial conceptualization of the study and the Centre for High Performance computing in Cape Town (http://www.chpc.ac.za) for providing part of the computational facilities.
Ethics
Human subjects: This study was approved as part of the IeDEA Southern African collaboration's protocol, by the Human Research Ethics Committee of the University of Cape Town, with a reference number N1810119 RECIP UCT 084/2006. Informed consent was obtained from all participants by the clinics collecting the data according to IeDEA protocols.
Copyright
© 2021, Ujeneza et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 731
- views
-
- 122
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Microbiology and Infectious Disease
Zika virus (ZIKV) infection causes significant human disease that, with no approved treatment or vaccine, constitutes a major public health concern. Its life cycle entirely relies on the cytoplasmic fate of the viral RNA genome (vRNA) through a fine-tuned equilibrium between vRNA translation, replication, and packaging into new virions, all within virus-induced replication organelles (vROs). In this study, with an RNA interference (RNAi) mini-screening and subsequent functional characterization, we have identified insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) as a new host dependency factor that regulates vRNA synthesis. In infected cells, IGF2BP2 associates with viral NS5 polymerase and redistributes to the perinuclear viral replication compartment. Combined fluorescence in situ hybridization-based confocal imaging, in vitro binding assays, and immunoprecipitation coupled to RT-qPCR showed that IGF2BP2 directly interacts with ZIKV vRNA 3’ nontranslated region. Using ZIKV sub-genomic replicons and a replication-independent vRO induction system, we demonstrated that IGF2BP2 knockdown impairs de novo vRO biogenesis and, consistently, vRNA synthesis. Finally, the analysis of immunopurified IGF2BP2 complex using quantitative mass spectrometry and RT-qPCR revealed that ZIKV infection alters the protein and RNA interactomes of IGF2BP2. Altogether, our data support that ZIKV hijacks and remodels the IGF2BP2 ribonucleoprotein complex to regulate vRO biogenesis and vRNA neosynthesis.
-
- Microbiology and Infectious Disease
Bacillus velezensis is a species of Bacillus that has been widely investigated because of its broad-spectrum antimicrobial activity. However, most studies on B. velezensis have focused on the biocontrol of plant diseases, with few reports on antagonizing Salmonella Typhimurium infections. In this investigation, it was discovered that B. velezensis HBXN2020, which was isolated from healthy black pigs, possessed strong anti-stress and broad-spectrum antibacterial activity. Importantly, B. velezensis HBXN2020 did not cause any adverse side effects in mice when administered at various doses (1×107, 1×108, and 1×109 CFU) for 14 days. Supplementing B. velezensis HBXN2020 spores, either as a curative or preventive measure, dramatically reduced the levels of S. Typhimurium ATCC14028 in the mice’s feces, ileum, cecum, and colon, as well as the disease activity index (DAI), in a model of infection caused by this pathogen in mice. Additionally, supplementing B. velezensis HBXN2020 spores significantly regulated cytokine levels (Tnfa, Il1b, Il6, and Il10) and maintained the expression of tight junction proteins and mucin protein. Most importantly, adding B. velezensis HBXN2020 spores to the colonic microbiota improved its stability and increased the amount of beneficial bacteria (Lactobacillus and Akkermansia). All together, B. velezensis HBXN2020 can improve intestinal microbiota stability and barrier integrity and reduce inflammation to help treat infection by S. Typhimurium.