Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 08 November 2022
Sec. Inflammation

Translating the observed differences in interleukin-6 levels between some antiretroviral regimens into potential long-term risk of serious non-AIDS events: A modeling study

Updated
Sergio Serrano-Villar,*&#x;Sergio Serrano-Villar1,2*†Calvin CohenCalvin Cohen3Jason V. Baker,Jason V. Baker4,5Maria Joo Janeiro&#x;Maria João Janeiro6†Filipa Arago,,&#x;Filipa Aragão6,7,8†Kathleen MelbourneKathleen Melbourne3Jose Luis GonzalezJose Luis Gonzalez9Laura LaraLaura Lara9Connie KimConnie Kim3Santiago Moreno,,&#x;Santiago Moreno1,2,10†
  • 1Hospital Universitario Ramón y Cajal, Infectious Diseases, Instituto de Investigación Sanitaria Ramón y Cajal (IRYCIS), Madrid, Spain
  • 2CIBERInfec, Instituto de Salud Carlos III, Madrid, Spain
  • 3HIV Medical Affairs, Gilead Sciences Inc., Foster City, CA, United States
  • 4Division of Infectious Diseases, Hennepin Healthcare Research Institute, Minneapolis, MN, United States
  • 5Department of Medicine, University of Minnesota, Minneapolis, MN, United States
  • 6Maple Health Group, New York, NY, United States
  • 7Incremental Action Consulting Lda, Lisbon, Portugal
  • 8NOVA National School of Public Health, Public Health Research Centre, Universidade NOVA de Lisboa, Lisbon, Portugal
  • 9HIV Medical Affairs, Gilead Sciences Inc., Madrid, Spain
  • 10Department of Medicine, Alcalá University, Madrid, Spain

Introduction: Variable levels of systemic inflammation are observed in people with HIV (PWH), but the clinical significance of differences among antiretroviral therapy (ART) regimens on associated levels of inflammatory markers is unclear. Based on data from previous epidemiologic studies that defined the predicted change in risk of serious non-AIDS events (SNAEs)/death by changes in interleukin-6 (IL-6) and D-dimer, we modeled the effects of differences in these markers between specific ART regimens on the long-term risk of clinical outcomes.

Methods: We used a Markov model to compare the risk of SNAEs/death with differences in IL-6 and D-dimer levels associated with remaining on specific three-drug regimens versus switching to specific two-drug ART regimens over 5 years of treatment. We used IL-6 and D-dimer data based on trajectories over time from the randomized TANGO and observational AIR studies. Age at model entry was set at 39 years. The primary endpoint was the number needed to treat for one additional SNAE/death.

Results: Over 3 years, PWH on one of the three-drug regimens studied were predicted to spend 13% more time in the low IL-6 quartile and 11% less time in the high IL-6 quartile compared with those on one of the two-drug regimens. Over 3 years, the predicted mean number of SNAEs/deaths per 100 PWH was 6.58 for a three-drug regimen associated with lower IL-6 levels versus 6.90 for a two-drug regimen associated with higher IL-6 levels. The number needed to treat for one additional SNAE/death among PWH receiving a two-drug versus three-drug regimen for 3 years was 81. Approximately 7,500 participants would be required for a 5-year clinical study to evaluate the accuracy of the model.

Conclusions: Our Markov model suggests that higher IL-6 levels associated with switching from specific three- to two- drug ART regimens may be associated with an increase in the risk of SNAEs/death. Clinical studies are warranted to confirm or refute these results.

Introduction

People with HIV (PWH) typically have higher levels of immune activation and systemic inflammation than those without (1, 2). Such elevations are observed despite effective antiretroviral therapy (ART) and preserved CD4+ counts, and may persist for life. Chronic immune activation and inflammation are associated with increased risks of serious non-AIDS events (SNAEs) and adverse outcomes (1), with impacts on life expectancy and quality of life (3, 4).

In 2014, UNAIDS established the ‘90-90-90’ target, aiming to achieve ≥90% rates of diagnosis, ART and viral suppression among PWH (5). The addition of good health-related quality of life in PWH with viral load suppression was proposed as a ‘fourth 90’ in 2016 (6). In 2020, the targets were updated to include 95% diagnosed, 95% on treatment, and 95% viral suppression by 2025, with at least 90% of PWH linked to services they need for their overall health and well-being (7). Among PWH on continuous ART, the spectrum of morbidity and mortality commonly reflects non-AIDS comorbidities (8, 9), and despite 95% targets, long-term elevations of inflammatory markers may have a detrimental clinical impact. In an analysis of data from PWH on suppressive oral ART in the continuous ART arms (without IL-2) of three large randomized studies – Strategies for Management of Antiretroviral Therapy (SMART) (10, 11), Evaluation of Subcutaneous Proleukin in a Randomized International Trial (ESPRIT), and Subcutaneous Recombinant, Human Interleukin-2 in HIV-Infected Patients with Low CD4+ Counts under Active Antiretroviral Therapy (SILCAAT) (12) – elevations in the inflammatory markers interleukin-6 (IL-6) and D-dimer were associated with higher risks of SNAEs (cardiovascular, hepatic or renal events, or malignancy) and death (8).

Levels of IL-6 have been analyzed in many studies of ART. Elevated IL-6 levels have been linked with adverse outcomes in the general population and PWH (13, 14). Different inflammatory responses in PWH receiving treatment may be due to the choice of ART regimen; differences in some inflammatory markers have been observed in randomized studies with ART switches (1517). While IL-6 levels are widely accepted as a robust predictor of adverse outcomes, the clinical significance of differences among ART regimens in their levels of IL-6 and other inflammatory markers is unclear.

Guidelines for the treatment of PWH currently recommend a variety of two- or three-drug ART regimens (18, 19). The choice of regimen is influenced by a variety of factors such as tolerability, short- and long-term safety, cost, patient convenience, genetic barriers to resistance, and baseline laboratory results. However, our understanding of which ART regimens produce the best overall outcomes among PWH, including long-term risk for SNAEs, is incomplete. Levels of immune activation and systemic inflammation may differ among regimens, which can contribute to long-term clinical outcomes. These could be overlooked in clinical trials with limited follow-up time that are powered to detect differences in viral load. The open-label, randomized, multicenter, phase 3 TANGO study compared continuing treatment with a three-drug, tenofovir alafenamide (TAF)-based regimen versus switching to a two-drug regimen (dolutegravir/lamivudine [DTG/3TC]) (20). The three-drug regimen was associated with significantly lower levels of IL-6 (P=0.006 and P=0.039 at Week 48 and 144, respectively) (2022). The AIR study, a nested study in the Spanish AIDS Research Network involving a prospective, multicenter cohort of PWH, reported similar results. IL-6 decreased over time among recipients of a three-drug boosted protease inhibitor (PI)-based or integrase strand transfer inhibitor (INSTI)-based regimen. After an inflection point 3 years from virologic suppression, IL-6 levels in recipients of the three-drug regimens were lower than those among individuals who had switched to a two-drug regimen (DTG + 3TC, DTG + rilpivirine [RPV], or a boosted PI + 3TC; P=0.012 for the difference between IL-6 trajectories). The results were unchanged in a subanalysis that excluded the PI group. Significant differences in the trajectories of high-sensitivity C-reactive protein (hsCRP) and D-dimer were also observed (23).

In the current analysis, we used data on IL-6 and D-dimer changes in the TANGO and AIR studies to mathematically model the degree to which observed differences in inflammatory marker levels between a three-drug ART regimen associated with lower IL-6 levels and a two-drug ART regimen associated with higher IL-6 levels could affect clinical outcomes in virologically suppressed PWH.

Materials and methods

Markov model

We created a Markov model to compare continued treatment of virologically suppressed PWH with a three-drug ART regimen for which two data sources found was associated with lower IL-6 levels versus switching to specific two-drug ART regimens that were associated with higher IL-6 levels. The model is based on observed differences in IL-6 and D-dimer levels between the three- and two-drug ART regimens studied in the TANGO and AIR studies. As such, the terms ‘lower IL-6’ and ‘higher IL-6’ regimens are used to refer to the specific respective three- and two-drug ART regimens used in these two studies and in this model. Participants in the TANGO study were screened for eligibility between January 18, 2018 and May 18, 2018; and participants in the AIR study initiated follow-up between January 02, 2004 and December 28, 2018.

Markov modeling (24, 25) provides a simple framework for linking biomarkers to outcomes over an extended time horizon when it is sufficient to work with average values within health states (in this case IL-6 quartiles), and when differences within those health states are not of key interest. The Markov framework assumes that clinical history is not relevant to inform the future (predicting clinical outcomes), i.e., patients’ future risk of an event depends on their current rather than previous IL-6 quartile. In our model, for each year, an individual could experience: i) no events, ii) an SNAE/death, or iii) death.

Health states at baseline were defined based on the distribution of IL-6 levels observed in the AIR study (23) (median 1.26; standard deviation [SD] 2.27; Serrano-Villar, data on file). Given the asymmetric nature of IL-6 distributions, four health states were created, adopting the quartiles of a gamma distribution as boundaries (‘low’, ‘low–medium’, ‘medium–high’ and ‘high’). The model adopts 48-week long cycles to match the annual visits with IL-6 data in the TANGO study. The risk of an event is influenced by IL-6 levels, and these differed between three-drug and two-drug ART regimens (8).

Our model was based on IL-6 data, including trajectories over time, from the TANGO and AIR studies (20, 21, 23). IL-6 data from TANGO were available from baseline to Week 144 and were used for the first 3 years, with IL-6 data from Weeks 144 to 288 of the AIR study used for subsequent time points. Quartiles for baseline IL-6 levels were defined according to the distribution in the AIR study cohort. Given the different nature of data from the TANGO and AIR studies, transition probabilities were calculated differently for the first and subsequent 3 years. TANGO provided ratios in mean IL-6 values for different time points/visits. These ratios were used to project the IL-6 distribution from one time point to the subsequent, by defining updated cut-offs of IL-6 quartiles. The ratios between previous and subsequent distributions (i.e., the proportion of individuals in each quartile in subsequent time points), were adopted as transition probabilities. In the AIR study, the slope of a regression of IL-6 values over time was used to define transitions across health states. Here, the distribution of IL-6 values was updated every cycle based on a slope adjusted for 48 weeks, therefore redistributing patients among health states. In accordance with the TANGO study, no changes over time in D-dimer levels were assumed. Thus, the model results were based both on the change in IL-6, as well as the lack of change in D-dimer. The TANGO data were from 741 PWH (three-drug regimen, n=372; two-drug regimen, n=369). The AIR study included 148 PWH from a large cohort who met the inclusion criterion and who had at least three plasma samples after HIV RNA suppression that were available for analysis (three-drug regimen, n=90; two-drug regimen, n=58). The three- and two-drug regimens included in the studies are shown in Table 1. In the TANGO study, the median age was 39 years for the three-drug regimen and 40 years for the two-drug regimen, and men accounted for 91% and 93% of the two cohorts, respectively (Table 1). In AIR, the mean ages were 37 and 40 years, and 87% and 86% of the patients were male, respectively. The age at entry in the model was set at 39 years, representing a weighted average from the main publications reporting IL-6 levels among PWH receiving different ART regimens (8, 26, 27).

TABLE 1
www.frontiersin.org

Table 1 Baseline demographics and clinical characteristics of the source data cohortsa (8, 20, 23).

We derived the risk of SNAE/death from analyses of data from the INSIGHT trials network, which enabled the risk to be defined according to changes in IL-6 and D-dimer concentrations. The INSIGHT trials analyses included data from 3,766 virologically suppressed individuals from the control arms of the randomized SMART, ESPRIT, and SILCAAT studies (i.e., on continuous ART without IL-2), in which SNAEs and all-cause death were the primary endpoints (8). Demographic characteristics of PWH in the INSIGHT trials were similar to TANGO and AIR, with a high proportion of males (79% versus 91% and 93% in TANGO, and 87% and 86% in AIR in PWH on three- and two-drug ART regimens, respectively) and median age of 42 years (versus 39 and 40 years in TANGO and 37 and 40 years in AIR in PWH on three- and two-drug ART regimens, respectively) (Table 1) (8). The cumulative risk of events per IL-6 quartile was calculated assuming a lack of change over time in the D-dimer level, based on the results of the TANGO study (20, 21). This was a conservative assumption for longer-term changes in inflammation, considering the differences of higher D-dimer on two ART regimens as noted in AIR and a higher rate of events in the INSIGHT analysis associated with increases in D-dimer (23).

Endpoints

The primary endpoint of the current analysis was the number needed to treat (NNT) for one additional SNAE/death at 144 and 240 weeks, based on the observed changes to IL-6. Predictions of time spent in each IL-6 quartile and numbers of SNAEs/deaths were considered as secondary endpoints, although these were required for calculation of the primary endpoint.

Size of future clinical cohort

We estimated the sizes of clinical cohorts needed for a randomized clinical trial or an observational study to investigate the accuracy of the NNT results using the Markov model. A Poisson distribution was used to calculate sample sizes needed to detect differences in the number of SNAEs/deaths between regimens. A 2-sided test with a significance level (α) of 0.05 and a power of 0.8 was adopted (28).

Ethics/participant consent

All clinical data used in this study were previously published, aggregated, and de-identified. Therefore, neither ethical approval nor patient consent was required or sought.

Results

Baseline characteristics

Baseline demographic and clinical characteristics of the source data cohorts are presented in Table 1. The median age and proportion of males were similar regardless of treatment regimen or study, but median CD4+ cell counts were higher in TANGO and INSIGHT studies than in the AIR study.

IL-6 trajectories from the Markov model

Within a time horizon of 3 years, PWH maintained on the three-drug ART regimens that were associated with lower IL6 levels (‘lower IL-6’ regimen) were predicted to spend 13% more time in the low IL-6 quartile and 11% less time in the high IL-6 quartile than those switching to a twodrug ART regimen associated with higher IL-6 levels (‘higher IL-6’ regimen) (Figure 1A). Differences between the two regimens in time spent in the middle two IL-6 quartiles were small. We observed a similar pattern over 240 weeks, although differences between the two regimens in the middle two quartiles were more pronounced (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1 Average time spent in each IL-6 quartile according to ART regimen at (A) 3 years and (B) 5 years, according to the model, which used the IL-6 trajectories seen in TANGO (n=741, up to Week 144) and AIR (n=148, Weeks 144–240).

Model-predicted SNAE/death rates and NNT

Overall, 260 participants in the INSIGHT studies experienced SNAE/death (8). The most common SNAEs/deaths were non-AIDS related cancer (37% of participants), cardiovascular disease (32%; including myocardial infarction [17%], stroke [7%], and deaths due to cardiovascular disease [7%]), decompensated liver cirrhosis (7%), end-stage renal disease (3%), and death due to other causes (21%; including AIDS [5%] and violence, accident or suicide [3%] (8).

Our model analysis of the 741 participants in the TANGO study showed that the predicted mean number of SNAEs/deaths per 100 PWH over 144 weeks was lower for continuous treatment with specific three-ART regimens versus switching to a two-ART regimen in three out of the four IL-6 quartiles (Figure 2). When all four IL-6 quartiles were combined, the predicted mean number of SNAEs/deaths per 100 PWH was 6.58 for continuous treatment with lower IL-6 compared with 6.90 with higher IL-6 (Figure 2). Results for the 5-year timeframe based on adding the 148 AIR study participants followed the same pattern, and the predicted mean numbers of SNAEs/deaths per 100 PWH were 14.87 and 17.49 for the three- and two-drug regimens, respectively. Based on these results, we calculated that for every 81 PWH treated for 5 years after switching to the two-drug regimen, there could be one additional SNAE/death as a consequence of higher IL-6, based on the association of IL-6 values and clinical outcomes from the INSIGHT studies. As shown in Table 2, the NNT decreased with increasing duration of treatment.

FIGURE 2
www.frontiersin.org

Figure 2 Incidence of serious non-AIDS events (cardiovascular, hepatic or renal event, or malignancy) or all-cause death according to IL-6 quartile and ART regimen over (A) 3 years and (B) 5 years, according to the model, which used the IL-6 trajectories seen in TANGO (n=741, up to Week 144) and AIR (n=148, Weeks 144–240).

TABLE 2
www.frontiersin.org

Table 2 Number needed to treat to observe one additional serious non-AIDS event (cardiovascular, hepatic or renal event, or malignancy) or all-cause death, with the two-drug ART regimen versus the three-drug regimen, by time on ART.

Clinical trial size required to validate model-predicted results

To evaluate the accuracy of this model in a clinical setting, we first assumed that a two-arm study (either a randomized study or an observational cohort) could be performed with study participants receiving treatment with either a three- versus two-drug regimen (i.e., a regimen associated with a ‘lower versus ‘higher IL-6) in a 1:1 ratio. On this basis, 7,591 participants would be required for a study with a 5-year treatment period. Variations in the sample size according to the study duration are shown in Table 3. Different sample sizes would be required for unequal group sizes (i.e., allocation ratios other than 1:1). For example, a total of almost 8,500 patients would be needed for a 5-year study with a three- versus two-drug regimen ratio of 2:1.

TABLE 3
www.frontiersin.org

Table 3 Estimated sample size required for a clinical study designed to support or refute the results of the Markov modeling.

Discussion

Inflammatory markers such as IL-6, D-dimer and hsCRP are elevated in some PWH, even after virologic suppression has been achieved (2). In some studies, the levels of such markers appear to differ among ART classes and, depending on the specific regimens involved, they may be higher during treatment with some two-drug ART regimens after switching from some three-drug ART regimens (2023). A range of factors may influence the level of systemic inflammation associated with an ART regimen. First, low-level HIV transcription in tissue reservoirs could trigger immune activation despite suppression of viral RNA levels in peripheral blood (29, 30). The extent of this phenomenon could be affected by the pharmacokinetics of the ART regimen (31). Second, the toxicity of certain ART regimens could increase the level of systemic inflammation to a greater extent than others, independent of the number of drugs being taken. Third, the ramifications of suboptimal adherence may lead to higher or more prolonged periods of low-level viremia in plasma, and this may differ according to the pharmacokinetic characteristics of different regimens. Suboptimal adherence has been associated with increased inflammation, even among individuals maintaining virologic suppression (3234). Over the long term, high levels of inflammation have been associated with adverse clinical outcomes among PWH; this motivated the current study modeling potential clinical differences in suppressive regimens and based on observed differences in inflammatory markers. Higher levels of IL-6 and/or D-dimer have been reported to confer increased risk of SNAEs and death (8).

We constructed a Markov model to predict SNAE/death rates using data from the TANGO and AIR studies. These two studies were selected for their long follow-up durations and because they both provided data on IL-6 and D-dimer levels in people switching between ART regimens. Both studies reported higher IL-6 levels in people switching to the specific two-drug regimen(s) in use versus continuing treatment with several three-drug regimens. Our model suggested that switching from some three-drug ART regimens associated with lower IL-6 levels to some two-drug ART regimens associated with higher IL-6 levels may increase the risk of SNAEs and/or death in PWH. Specifically, one additional SNAE/death was predicted to occur for every 43 people treated for 240 weeks or for every 106 people treated for 144 weeks.

The results of our model are specific to the treatment regimens used in the TANGO trial and AIR study and cannot be extrapolated to other ART combinations. More recent two-drug regimens (e.g., long-acting injectable drugs with longer half-lives) may be associated with different levels of systemic inflammation. In other studies, switches from PI- to INSTI-based ART regimens were associated with differences in inflammatory markers, including decreases in not only IL-6 but also soluble CD14 (sCD14), intestinal fatty acid-binding protein, hsCRP, monocyte chemoattractant protein 1, osteoprotegerin, tumor necrosis factor α, and D-dimer (1517), suggesting that differences in rates of SNAEs/deaths could also be observed among the regimens of these studies. Conversely, other clinical trial data found no differences in IL-6 levels over 48 weeks between patients switching from multiple three-drug regimens to DTG + 3TC versus those continuing three-drug regimens (35). While this suggests that no differences would be seen in the SNAE/death rates, the possibility remains that differences could emerge over a longer duration of follow-up or in subsets of the regimens used at entry.

IL-6 and other inflammatory markers are not routinely assessed in PWH, and the impact of inflammation on clinical outcomes is not clearly understood. International guidelines acknowledge the importance of inflammation in PWH but do not provide specific recommendations for the measurement, prevention, or treatment of inflammation (18, 19). This is partly because there is not enough evidence to show that interventions intended to change heightened inflammation levels affect clinical outcomes. This may also be in part due to the relatively small number of studies that examined differences in inflammatory markers with different regimens. Our study helps provide insight into the potential clinical consequences of increased inflammation associated with certain ART regimens, with guidance toward sample sizes of studies, including cohorts, that can be estimated to address this concern.

To contextualize the present results, we draw a comparison with statin therapy, an intervention proven to be beneficial when used at the population level. A review of data from 19 randomized trials showed that the NNT for statins to prevent composite cardiovascular outcomes over a treatment period of 1 to 6 years was 72 (36). In a Cochrane review, the adjusted 5-year NNT was 96 for prevention of death and 56 for prevention of a coronary heart disease event (either fatal or non-fatal) (37). Although the statin data are clearly defined by clinical outcomes in large numbers of patients, whereas the data in the present study are based on our model, the available evidence suggests that preferential use of ART regimens with reduced systemic inflammation could provide a magnitude of overall benefit consistent with other accepted medical interventions, such as statin therapy.

While the current study provides predictions based on the available data, additional data are warranted to refute or confirm the observed differences in inflammatory markers between regimens. In addition, clinical studies are necessary to refute or confirm the predicted effects on clinical outcomes of lower levels versus higher levels of inflammation while on ART. Our model indicates that randomized controlled or observational cohort trials would need to include large numbers of patients and a long treatment period (e.g., almost 3,000 patients for a two-arm, 5-year study). Considering the practical and financial challenges of conducting studies of this scale, a cohort-based approach appears more likely to provide additional clinical evidence to verify or refute these findings.

The principal strength of this study was the use of long-term results from large, parallel-group studies conducted in PWH. Thus, the risks of SNAEs/death were predicted from a robust clinical dataset. Another strength was the assumption of no changes in D-dimer levels (as seen in the shorter TANGO study); given the differences in D-dimer levels observed in the longer AIR study, and since both IL-6 and D-dimer have been found to independently contribute to the risk of severe non-AIDS events (38), this was a conservative approach for the longer-term consequences of changes in inflammation. In addition, as there were large differences in CD4+ counts between the source data cohorts, the results from the present model can be extrapolated to patients across a wide CD4+ count range.

Due to differences in study population and design, data from the TANGO and AIR studies from baseline to Week 144 were not combined. Other limitations include the theoretical nature of the model and consideration of only IL-6 (change) and D-dimer (lack of change) as inflammatory markers. We note that a significantly greater reduction in sCD14 was observed at two time points with two-drug versus three-drug ART regimens at 144 weeks in the TANGO study (20, 21, 39). The clinical relevance of this difference is uncertain, since an analysis of sCD14 in HIV infection showed that high sCD14 levels were associated with an increased risk of mortality but not with SNAEs, and only when results for both sCD14 and IL-6 changed in the same direction, in contrast to what was observed in TANGO (40). We are unaware of data that would have allowed us to model contradictory effects of these markers, so additional work would be needed to understand and incorporate that finding. However, there is evidence to show that IL-6 is the most robust inflammatory marker for predicting clinical events, irrespective of variations in other biomarkers (13). Epidemiologic data estimating their impact on incident SNAEs are less robust than the data for IL-6. An intrinsic limitation of the Markov model is that it does not account for the fact that inflammation over time could affect the disease outcomes. However, the START, SMART, ESPRIT and SILCAAT studies have demonstrated that a single measurement of inflammatory markers can predict long-term outcomes (up to 60 months and 9 years, respectively) (8, 26). We acknowledge that inflammation in PWH is multifactorial, and this study was focused only on the inflammatory effects of several specific ART regimens. The IL-6 data (and therefore the clinical outcomes) presented here may have limited applicability to female PWH due to the predominance of men in the TANGO and AIR studies, and may not be generalizable to older PWH, in whom the metabolism and elimination of ART may be slower than in younger PWH. Indeed, participant demographic differences may have affected levels of inflammatory markers, as may lifestyle factors; these were not considered in our model. Finally, our results should not be generalized to all two- and three-drug regimens; our results are only applicable to the two- and three-drug regimens associated with higher and lower IL-6, respectively, and in particular might not be generalizable to newer ART combinations.

In conclusion, our Markov model suggests that an increase in systemic inflammation associated with switching from a three-drug ART regimen associated with lower IL-6 levels to a two-drug ART regimen associated with higher IL-6 levels may increase the risk of SNAEs and/or death in PWH. Based on observed differences between regimens in IL-6 levels and no change to D-dimer as observed in the TANGO study, one additional SNAE/death was predicted to occur for every 43 PWH aged 30–50 years treated for 240 weeks with a regimen associated with higher IL-6 levels rather than maintaining a regimen associated with lower IL-6 levels. These predictions were based on the assumptions that an individual’s future risk of an event depends on their current rather than previous IL-6 quartile, and that variation in IL-6 levels in TANGO and AIR studies were due to differences in the ART regimens rather than other differences in, for example, ART adherence, comedication, classical risk factors or previous medical history. As TANGO was a large, randomized study, most of these factors are expected to be balanced between the arms. Nevertheless, it was an open-label study, and there may be other factors that have an impact on the observed findings on IL-6. Additional data sets are needed on these and other regimens to confirm or refute the inflammatory marker findings in TANGO and AIR. Finally, as this is a model, clinical studies would be needed to confirm our results. However, if verified, the selection of ART regimens associated with lower levels of systemic inflammation may become an important factor when optimizing long-term care of PWH.

Data availability statement

The datasets presented in this article are not readily available because Gilead Sciences shares data upon request or as required by law or regulation with qualified external researchers based on submitted curriculum vitae and reflecting non-conflict of interest. The request proposal must also include a statistician. Approval of such requests is at Gilead Sciences’ discretion and is dependent on the nature of the request, the merit of the research proposed, the availability of the data, and the intended use of the data. Requests to access the datasets should be directed to Gilead Sciences, ZGF0YXJlcXVlc3RAZ2lsZWFkLmNvbQ==.

Ethics statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

SS-V, CC, and SM were involved in the study design and data collection. All authors contributed to the article and approved the submitted version.

Funding

This is a collaborative work, funded by Gilead Sciences, Inc (Foster City, CA, USA).

Acknowledgments

The authors thank the investigators and participants who provided data for the original studies. Medical writing support, including the development of a draft outline and subsequent drafts in consultation with the authors, assembling tables and figures, collating author comments, copyediting, fact checking, and referencing was provided by Emma McConnell, PhD and Heather Davies, PhD CMPP at Aspire Scientific (Bollington, UK), and funded by Gilead Sciences, Inc (Foster City, CA, USA).

Conflict of interest

SS-V received research grants from Gilead Sciences and MSD; non-financial support or fees from consulting or preparation of educational materials from Gilead Sciences, Janssen Cilag, MSD and ViiV Healthcare. CC, KM, JLG, LL and CK are employed by and own stocks/shares in Gilead Sciences. FA was employed by Maple Health Group. JVB received a research grant from Gilead for a study of remdesivir, 2020. SM has been involved in speaking activities and received grants for research from Gilead, Janssen Cilag, Merck, Sharp & Dohme and ViiV Healthcare.

The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The authors declare that this study received funding from Gilead Sciences. The funder had the following involvement with the study: the study design, collection, analysis, interpretation of data, the writing of this article and the decision to submit it for publication.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Zicari S, Sessa L, Cotugno N, Ruggiero A, Morrocchi E, Concato C, et al. Immune activation, inflammation, and non-AIDS co-morbidities in HIV-infected patients under long-term ART. Viruses (2019) 11:200. doi: 10.3390/v11030200

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Neuhaus J, Jacobs DR Jr., Baker JV, Calmy A, Duprez D, La Rosa A, et al. Markers of inflammation, coagulation, and renal function are elevated in adults with HIV infection. J Infect Dis (2010) 201:1788–95. doi: 10.1086/652749

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Portilla-Tamarit J, Reus S, Portilla I, Fuster Ruiz-de-Apodaca MJ, Portilla J. Impact of advanced HIV disease on quality of life and mortality in the era of combined antiretroviral treatment. J Clin Med (2021) 10:716. doi: 10.3390/jcm10040716

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Kuller LH, Tracy R, Belloso W, De Wit S, Drummond F, Lane HC, et al. Inflammatory and coagulation biomarkers and mortality in patients with HIV infection. PloS Med (2008) 5:e203. doi: 10.1371/journal.pmed.0050203

PubMed Abstract | CrossRef Full Text | Google Scholar

5. UNAIDS. An ambitious treatment target to help end the AIDS epidemic (2014). Available at: https://www.unaids.org/en/resources/909090 (Accessed May 11, 2022).

Google Scholar

6. Lazarus JV, Safreed-Harmon K, Barton SE, Costagliola D, Dedes N, del Amo Valero J, et al. Beyond viral suppression of HIV – the new quality of life frontier. BMC Med (2016) 14:94. doi: 10.1186/s12916-016-0640-4

PubMed Abstract | CrossRef Full Text | Google Scholar

7. UNAIDS. World AIDS day report (2020). Available at: https://www.unaids.org/en/resources/documents/2020/prevailing-against-pandemics (Accessed May 11, 2022).

Google Scholar

8. Grund B, Baker JV, Deeks SG, Wolfson J, Wentworth D, Cozzi-Lepri A, et al. Relevance of interleukin-6 and d-dimer for serious non-AIDS morbidity and death among HIV-positive adults on suppressive antiretroviral therapy. PloS One (2016) 11:e0155100. doi: 10.1371/journal.pone.0155100

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Mocroft A, Reiss P, Gasiorowski J, Ledergerber B, Kowalska J, Chiesi A, et al. Serious fatal and nonfatal non-AIDS-defining illnesses in Europe. J Acquir Immune Defic Syndr (2010) 55:262–70. doi: 10.1097/QAI.0b013e3181e9be6b

PubMed Abstract | CrossRef Full Text | Google Scholar

10. El-Sadr WM, Lundgren J, Neaton JD, Gordin F, Abrams D, Arduino RC, et al. CD4+ count-guided interruption of antiretroviral treatment. N Engl J Med (2006) 355:2283–96. doi: 10.1056/NEJMoa062360

PubMed Abstract | CrossRef Full Text | Google Scholar

11. El-Sadr WM, Grund B, Neuhaus J, Babiker A, Cohen CJ, Darbyshire J, et al. Risk for opportunistic disease and death after reinitiating continuous antiretroviral therapy in patients with HIV previously receiving episodic therapy: a randomized trial. Ann Intern Med (2008) 149:289–99. doi: 10.7326/0003-4819-149-5-200809020-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Abrams D, Lévy Y, Losso MH, Babiker A, Collins G, Cooper DA, et al. Interleukin-2 therapy in patients with HIV infection. N Engl J Med (2009) 361:1548–59. doi: 10.1056/NEJMoa0903175

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Borges AH, O’Connor JL, Phillips AN, Neaton JD, Grund B, Neuhaus J, et al. Interleukin 6 is a stronger predictor of clinical events than high-sensitivity c-reactive protein or d-dimer during HIV infection. J Infect Dis (2016) 214:408–16. doi: 10.1093/infdis/jiw173

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Tanaka T, Narazaki M, Kishimoto T. IL-6 in inflammation, immunity, and disease. Cold Spring Harb Perspect Biol (2014) 6(10):a016295. doi: 10.1101/cshperspect.a016295

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lake J, Currier J, Koteff J, Brennan C, Granier C, Gartland M, et al. Cardiovascular biomarkers after switching to ABC/DTG/3TC: The STRIIVING. Poster 660 presented at the 23rd Conference on Retroviruses and Opportunistic Infections (CROI);(2016) Boston, MA, USA.

Google Scholar

16. Martínez E, D’Albuquerque PM, Llibre JM, Gutierrez F, Podzamczer D, Antela A, et al. Changes in cardiovascular biomarkers in HIV-infected patients switching from ritonavir-boosted protease inhibitors to raltegravir. AIDS (2012) 26:2315–26. doi: 10.1097/QAD.0b013e328359f29c

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Trottier B, Lake JE, Logue K, Brinson C, Santiago L, Brennan C, et al. Dolutegravir/abacavir/lamivudine versus current ART in virally suppressed patients (STRIIVING): a 48-week, randomized, non-inferiority, open-label, phase IIIb study. Antivir Ther (2017) 22:295–305. doi: 10.3851/imp3166

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Department of Health and Human Services Panel on Antiretroviral Guidelines for Adults and Adolescents. Guidelines for the use of antiretroviral agents in adults and adolescents with HIV (2021). Available at: https://clinicalinfo.hiv.gov/sites/default/files/guidelines/archive/AdultandAdolescentGL_2021_08_16.pdf (Accessed October 22, 2021).

Google Scholar

19. European AIDS Clinical Society. Guidelines version 10.1 (2020). Available at: https://www.eacsociety.org/media/guidelines-10.1_30032021_1.pdf (Accessed May 11, 2022).

Google Scholar

20. van Wyk J, Ajana F, Bisshop F, De Wit S, Osiyemi O, Portilla Sogorb J, et al. Efficacy and safety of switching to dolutegravir/lamivudine fixed-dose 2-drug regimen vs continuing a tenofovir alafenamide-based 3- or 4-drug regimen for maintenance of virologic suppression in adults living with human immunodeficiency virus type 1: Phase 3, randomized, noninferiority TANGO study. Clin Infect Dis (2020) 71:1920–9. doi: 10.1093/cid/ciz1243

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Osiyemi O, Ajana F, Bisshop F, De Wit S, Portilla J, Routy J-P, et al. Switching to DTG/3TC fixed-dose combination (FDC) is non-inferior to continuing a TAF-based regimen in maintaining virologic suppression through 144 weeks (TANGO study). Abstract 900 presented at IDWeek 2021; September 29 to October 03, 2021; Virtual Meeting. (2021).

Google Scholar

22. Wang R, Wright J, George N, Ait-Khaled M, Lutz T, Osiyemi O, et al. Comparison of viral replication for the 2-drug regimen (2DR) of dolutegravir/lamivudine (DTG/3TC) versus a 3/4 drug tenofovir alafenamide-based regimen (TBR) in the TANGO study through week, 96 Abstract presented at the 11th IAS Conference on HIV Science; July 18–21, 2021; Virtual Meeting.

Google Scholar

23. Serrano-Villar S, Lopez-Huertas MR, Jimenez D, Galera C, Martinez-Sanz J, Moreno E, et al. Long-term changes of inflammatory biomarkers in individuals on suppressive three-drug or two-drug antiretroviral regimens. Front Immunol (2022) 13:848630. doi: 10.3389/fimmu.2022.848630

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Cowles MK. Bayesian Estimation of the proportion of treatment effect captured by a surrogate marker. Stat Med (2002) 21:811–34. doi: 10.1002/sim.1057

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Sonnenberg FA, Beck JR. Markov Models in medical decision making: a practical guide. Med Decis Making (1993) 13:322–38. doi: 10.1177/0272989X9301300409

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Baker JV, Sharma S, Grund B, Rupert A, Metcalf JA, Schechter M, et al. Systemic inflammation, coagulation, and clinical risk in the START trial. Open Forum Infect Dis (2017) 4:ofx262. doi: 10.1093/ofid/ofx262

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mussini C, Lorenzini P, Cozzi-Lepri A, Lapadula G, Marchetti G, Nicastri E, et al. CD4/CD8 ratio normalisation and non-AIDS-related events in individuals with HIV who achieve viral load suppression with antiretroviral therapy: an observational cohort study. Lancet HIV (2015) 2:e98–106. doi: 10.1016/s2352-3018(15)00006-5

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gu K, Ng HK, Tang ML, Schucany WR. Testing the ratio of two poisson rates. Biom J (2008) 50:283–98. doi: 10.1002/bimj.200710403

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Klatt NR, Chomont N, Douek DC, Deeks SG. Immune activation and HIV persistence: implications for curative approaches to HIV infection. Immunol Rev (2013) 254:326–42. doi: 10.1111/imr.12065

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Lv T, Cao W, Li T. HIV-Related immune activation and inflammation: Current understanding and strategies. J Immunol Res (2021) 2021:7316456. doi: 10.1155/2021/7316456

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Fletcher CV, Staskus K, Wietgrefe SW, Rothenberger M, Reilly C, Chipman JG, et al. Persistent HIV-1 replication is associated with lower antiretroviral drug concentrations in lymphatic tissues. Proc Natl Acad Sci USA (2014) 111:2307–12. doi: 10.1073/pnas.1318249111

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Castillo-Mancilla JR, Brown TT, Erlandson KM, Palella FJ Jr., Gardner EM, Macatangay BJ, et al. Suboptimal adherence to combination antiretroviral therapy is associated with higher levels of inflammation despite HIV suppression. Clin Infect Dis (2016) 63:1661–7. doi: 10.1093/cid/ciw650

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Castillo-Mancilla JR, Phillips AN, Neaton JD, Neuhaus J, Collins S, Mannheimer S, et al. Association of suboptimal antiretroviral therapy adherence with inflammation in virologically suppressed individuals enrolled in the SMART study. Open Forum Infect Dis (2018) 5:ofx275. doi: 10.1093/ofid/ofx275

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Castillo-Mancilla JR, Phillips AN, Neaton JD, Neuhaus J, Sharma S, Baker JV, et al. Incomplete ART adherence is associated with higher inflammation in individuals who achieved virologic suppression in the START study. J Int AIDS Soc (2019) 22:e25297. doi: 10.1002/jia2.25297

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Llibre JM, Brites CA, Cheng C-Y, Osiyemi O, Galera C, Hocquelou xL, et al. Switching to the 2-drug regimen of dovato (dolutegravir/lamivudine) fixed-dose combination is non-inferior to continuing a 3-drug regimen through 48 weeks in a randomized clinical trial (SALSA). Abstract presented at the 11th IAS Conference on HIV Science; July 18–21, 2021; Virtual Meeting.

Google Scholar

36. Chou R, Dana T, Blazina I, Daeges M, Bougatsos C, Grusing S, et al. Statin use for the prevention of cardiovascular disease in adults: a systematic review for the U.S. preventive services task force. In: US Preventive services task force evidence syntheses, formerly systematic evidence reviews. US: Agency for Healthcare Research and Quality (2016). Evidence Syntheses, No. 139. Rockville (MD).

Google Scholar

37. Taylor F, Huffman MD, Macedo AF, Moore TH, Burke M, Davey Smith G, et al. Statins for the primary prevention of cardiovascular disease. Cochrane Database Syst Rev (2013) 2013:Cd004816. doi: 10.1002/14651858.CD004816.pub5

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Tenorio AR, Zheng Y, Bosch RJ, Krishnan S, Rodriguez B, Hunt PW, et al. Soluble markers of inflammation and coagulation but not T-cell activation predict non-AIDS-defining morbid events during suppressive antiretroviral treatment. J Infect Dis (2014) 210:1248–59. doi: 10.1093/infdis/jiu254

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Llibre JM, Cahn PE, Lo J, Barber TJ, Mussini C, van Welzen BJ, et al. Changes in inflammatory and atherogenesis biomarkers with the 2-drug regimen dolutegravir plus lamivudine in antiretroviral therapy-experienced, virologically suppressed people with HIV-1: a systematic literature review. Open Forum Infect Dis (2022) 9:ofac068. doi: 10.1093/ofid/ofac068

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Sandler NG, Wand H, Roque A, Law M, Nason MC, Nixon DE, et al. Plasma levels of soluble CD14 independently predict mortality in HIV infection. J Infect Dis (2011) 203:780–90. doi: 10.1093/infdis/jiq118

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: antiretroviral, HIV, inflammation, interleukin-6, Markov, three-drug regimen, two-drug regimen

Citation: Serrano-Villar S, Cohen C, Baker JV, Janeiro MJ, Aragão F, Melbourne K, Gonzalez JL, Lara L, Kim C and Moreno S (2022) Translating the observed differences in interleukin-6 levels between some antiretroviral regimens into potential long-term risk of serious non-AIDS events: A modeling study. Front. Immunol. 13:976564. doi: 10.3389/fimmu.2022.976564

Received: 24 June 2022; Accepted: 19 October 2022;
Published: 08 November 2022.

Edited by:

Jean-Pierre Routy, McGill University, Canada

Reviewed by:

Pierre Corbeau, Université de Montpellier, France
Scott Letendre, University of California, San Diego, United States

Copyright © 2022 Serrano-Villar, Cohen, Baker, Janeiro, Aragão, Melbourne, Gonzalez, Lara, Kim and Moreno. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sergio Serrano-Villar, c2VyZ2lvLnNlcnJhbm9Ac2FsdWQubWFkcmlkLm9yZw==

ORCID: Sergio Serrano-Villar, orcid.org/0000-0002-5447-3554
Maria João Janeiro, orcid.org/0000-0003-4794-866X
Filipa Aragão, orcid.org/0000-0003-3446-6725
Santiago Moreno, orcid.org/0000-0002-2843-1094

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.