Inflammatory biomarkers and abacavir use in the Women's
Interagency HIV Study and the Multicenter AIDS Cohort Study
17 July 2010 - Volume 24 - Issue 11 - p 1657-1665
Palella, Frank J Jr; Gange, Stephen J; Benning, Lorie; Jacobson, Lisa; Kaplan, Robert C; Landay, Alan L; Tracy, Russell P; Elion, Richard
aNorthwestern University Feinberg School of Medicine, Chicago, Illinois, USA
bJohns Hopkins University, Baltimore, USA
cAlbert Einstein College of Medicine, New York, New York, USA
dRush University Medical Center, Chicago, Illinois, USA
eUniversity of Vermont, Burlington, Vermont, USA
fGeorge Washington University, Washington, District of Columbia, USA.
Objective: To assess associations between abacavir (ABC) use and systemic inflammation.
Design: Nested case-control study.
Methods: The Multicenter AIDS Cohort Study (MACS) and Women's Interagency HIV Study (WIHS) cohort participants who initiated ABC were matched, using propensity score methods, to ABC-unexposed persons. Levels of high-sensitivity C-reactive protein (hsCRP) (µg/ml), interleukin-6 (IL-6) (pg/ml), and D-dimer (µg/ml) were measured from pre-HAART and on-HAART plasma. Random-effects models compared markers by ABC exposure and by changes from pre-HAART levels.
Results: Biomarkers were measured in N = 508 matched pairs (328 women; 180 men). Pre-HAART levels did not differ by exposure group except that hsCRP levels were higher among WIHS women who subsequently used ABC (P = 0.04). Regardless of ABC use, mean hsCRP increases and D-dimer reductions were seen when comparing pre-HAART to on-HAART levels, in the overall group (28 and -27%), for MACS men (28 and -31%) and for WIHS women [29 and -24%, P < 0.01 for all]; IL-6 levels declined in MACS men (P = 0.02). No adjusted biomarker level differences existed by ABC exposure at the on-HAART visit. HIV RNA reductions correlated with D-dimer (r = 0.14, P < 0.01) and IL-6 (r = 0.12, P < 0.01) reductions. Associations between ABC use and mean biomarker levels were modified by pre-HAART antiretroviral therapy experience. Renal dysfunction was equally likely among non-ABC and ABC recipients.
Discussion: ABC use was not associated with plasma elevations in hsCRP, IL-6, and D-dimer. Mechanisms other than increased systemic inflammation may account for ABC's reported association with increased cardiovascular disease. HAART-associated reductions in D-dimer and IL-6 were apparent regardless of ABC use and were correlated with HIV RNA reductions.
Reports from several large observational cohorts of HIV-infected persons have linked the recent use (within 6 months) of certain nucleoside reverse transcriptase inhibitor antiretrovirals (NRTIs), most notably abacavir (ABC), with increased risk for coronary events [1-8]. In these groups, recent ABC use was strongly and independently associated with increased risk for the occurrence of myocardial infarctions.
Subsequent analyses of data from the Strategies for Management of Anti-Retroviral Therapy (SMART) study suggested that ABC use was independently associated with increased plasma levels of high-sensitivity C-reactive protein (hsCRP) and interleukin 6 (IL-6), biomarkers that are associated with systemic inflammation [2-4]. These observations have led to hypotheses that ABC's unique association with increased cardiovascular-related death is mediated pathophysiologically, at least in part, through its pro-inflammatory effects. However, data from several other studies have failed to demonstrate associations between ABC use and either increases in cardiovascular disease (CVD) or increases in systemic levels of inflammatory biomarkers [9-12].
To investigate whether use of ABC independently increases an inflammatory response, we comparatively and longitudinally evaluated plasma levels of inflammatory biomarkers among HIV-infected persons who did vs. did not receive ABC-containing HAART regimens. Stored specimens collected prospectively from before and after HAART initiation permitted the examination of within-person changes. These within-cohort and cross-cohort comparisons were undertaken among women participating in the Women's Interagency HIV Study (WIHS) and the Multicenter AIDS Cohort Study (MACS).
In the WIHS and the MACS, we observed changes in the levels of pro-inflammatory biomarkers hsCRP, IL-6 and D-dimer between the baseline (pre-HAART) and index (on HAART) visits, a mean elapsed time of 4.2 years. D-Dimer and IL-6 levels decreased, whereas hsCRP levels increased. These changes were comparable among persons who initiated ABC-containing vs. non-ABC-containing HAART and reflected, at least in part, the effect of HAART in reducing plasma HIV viral burden.
This work should not be interpreted as either refuting or supporting hypotheses suggesting associations between recent use of abacavir (or any other ART for that matter) and CVD in general or specific CVD endpoints. We did not undertake evaluations of use of ABC and CVD nor the assessment of plasma biomarker levels in relationship to specific clinical events. Such evaluations necessarily require evaluations of large numbers of patients followed over long periods, during which sufficiently substantial numbers of clinical endpoints occur so as to permit meaningful comparisons based on exposure. However, our work does suggest that, if recent ABC use is indeed associated with increased risk for adverse cardiovascular events, systemic inflammation is not likely the sole or primary means by which its effects are mediated. It remains entirely possible, as preliminary work from other groups has suggested, that ABC's effects upon vasculature may be mediated by other mechanisms such as reduced endothelial function, that is, flow-mediated arterial dilation .
Changes in biomarker levels observed from the pre-HAART to the on-HAART visit suggest that there may be an effect upon biomarkers exerted by HAART use itself that, at least in terms of the observed reductions in D-dimer levels (as well as the reductions in IL-6 levels seen among MACS men, a group in whom better virologic control was achieved overall compared with WIHS women), corroborates observations that HIV itself, or more specifically HIV plasma viremia, is associated with plasma elevations in inflammatory biomarkers . Indeed, when we examined the association of HIV RNA changes with biomarker level changes between the two visits, we verified that greater HIV RNA reductions were significantly correlated with reductions in D-dimer and IL-6.
Observed increases in plasma levels of hsCRP from the pre-HAART to on-HAART visit are more difficult to explain and may reflect the influence of multiple other non-HIV-related cofactors and co-morbidities that themselves can increase biomarker levels. For example, among WIHS participants, the extent to which higher prevalences of chronic co-infection with HCV (more than one-third of participants overall), smoking (nearly one-half of participants), obesity, or recreational drug use may have accounted for observed differences in plasma hsCRP levels is unclear . Similarly, among MACS men, who are on average older than WIHS women, non-HIV cardiovascular comorbidities such as hypertension (seen in one-quarter of participants) or even aging itself may have exerted an effect on hsCRP levels. Data exist from the MACS demonstrating a positive, although modest, association between plasma hsCRP levels and plasma lipids, including non-HDL cholesterol and high-density lipoprotein cholesterol , and age-associated hsCRP increases .
In analyses limited to persons whose first HAART was also the first exposure to ART, a group representative of contemporary HAART initiators, there were no significant differences in pre-HAART to on-HAART biomarker level changes when comparing ABC and non-ABC HAART users, although adjusted biomarker level changes differences were more pronounced in this 'ART naive prior to HAART' group than in the overall group.
Of note, despite biomarker level changes seen from pre-HAART to on-HAART visits and the cohort-based differences observed, mean absolute plasma levels of each biomarker were within generally accepted 'normal' ranges for healthy adults, both at baseline and index visits. In this regard, it becomes somewhat more difficult to assign a clinical relevance to HAART-associated changes or cohort-associated differences in observed biomarker levels, particularly because we did not (and could not because of a low number of observed events) correlate our findings with risks for cardiovascular events.
There were cohort-based differences in biomarker levels and changes in levels (in the case of IL-6) over time. WIHS women had significantly lower IL-6 and CRP levels at the baseline visit and MACS men experienced greater IL-6 decreases. It is difficult to surmise upon the extent to which these differences are consistent, predictable, and based on sex vs. reflective of differences intrinsic to these specific cohorts. MACS men and WIHS women differed in important ways other than sex. WIHS women were younger, more likely to be nonwhite, more likely to be chronically co-infected with HCV, had higher BMIs, lower CD4 cell count, and higher HIV RNA levels at both baseline and index visits than MACS men. Future work will need to further explore whether sex-based biomarker level differences reliably exist and the extent to which cohort-based differences in other non-HIV factors modify these differences.
There were other limitations to our analysis. 'official' reference ranges for the plasma biomarkers measured do not exist and values obtained may vary depending on several factors including assay reagents and biological samples used. Furthermore, the statistical models that we employed do not account for measurements that are outside assay ranges; however, given the small proportion of observations that were outside of 'normal' ranges, we believe it unlikely that they would have influenced our interpretations.
Obviously, we did not evaluate mechanisms other than increases in systemic inflammation (or, in the case of D-dimer, coagulation, and fibrinolysis) by which ABC use may enhance CVD risk. Several recent reports [29,30] have even suggested an independent role of HIV infection as a risk for atherosclerosis, with the magnitude of increased risk unclear but possibly similar to that which has been attributed to smoking, diabetes mellitus, or advancing age .
In summary, changes in plasma levels of the proinflammatory biomarkers hsCRP, D-dimer, and IL-6 were not observed in association with ABC-based HAART use. Pathophysiologic mechanisms other than increases in systemic inflammation should be explored to explain reported increased rates of myocardial infarctions observed in association with ABC use. We did note that HAART use in general was associated with decreases in D-dimer and IL-6, which may reflect reduced systemic inflammation consequent to HAART-induced HIV suppression, and increases in hsCRP, which are more difficult to explain. Our observed cohort-based differences in biomarker levels beg for further exploration of the effects of sex and of other non-HIV-related co-morbidities upon systemic inflammation in HIV-infected persons.
We measured hsCRP, D-dimer, and IL-6 levels in 508 propensity score matched participant pairs (ABC-exposed vs. non-ABC exposed persons); 180 pairs contributed by MACS men who first reported HAART between November 1995 and March 2008 and 328 pairs contributed by WIHS women who first reported HAART between September 1996 and February 2008. The mean elapsed time between the pre-HAART and on-HAART visits was 4.2 years. Characteristics included in the propensity score models as well as additional patient features are shown in Table 1 for MACS men and in Table 2 for WIHS women. In these tables, characteristics of HAART users initiating ABC are compared to eligible person-visits from all HAART users not reporting ABC use (middle columns), and to those persons selected using propensity-score matching methods (right columns). Prior to matching, MACS men who received ABC-containing HAART had significantly more pre-HAART ART experience, a higher proportion of visits not on HAART, lower median CD4 cell count, higher median plasma HIV RNA levels, and a lower proportion of visits with HIV RNA less than or equal to 500 copies/ml when compared to MACS men who received non-ABC-containing HAART. After matching, the two MACS exposure groups (ABC and non-ABC HAART recipients) were comparable on all matching characteristics except the duration of HAART use (time from HAART initiation to on-HAART visit). MACS ABC users and nonusers also had comparable nadir CD4 cell count, evidence of hypertension, history of diabetes mellitus, and report of CVD.
Prior to ABC initiation, WIHS women who received ABC-containing HAART had lower median serum HDL cholesterol levels, a higher proportion of visits during which they were not receiving HAART, lower median CD4 cell count, higher median plasma HIV RNA levels, and a lower proportion of visits with HIV RNA 500 copies/ml or less when compared with non-ABC HAART person-visits. After matching, WIHS exposure groups were comparable on all matching characteristics, nadir CD4 cell count, evidence of hypertension, and report of CVD. However, a higher proportion of ABC users had a history of diabetes mellitus (16 vs. 10%, P = 0.03).
Comparisons of MACS men and WIHS women illustrated cohort differences other than biological sex: WIHS women were more likely to be nonwhite, HCV co-infected, younger, have lower CD4 cell count, greater median BMI, and a history of diabetes mellitus (P < 0.01 for all comparisons). For these reasons, although the main analysis was conducted on the combined data and adjusted for cohort, we also undertook analyses stratified by cohort.
Figure 2 displays the geometric mean marker levels for the combined and individual cohorts. Pre-HAART biomarker levels in the ABC-exposed and ABC-unexposed groups are displayed as a single mean, as they were similar and nearly all were statistically indistinguishable; the only exception was observed among WIHS participants for whom hsCRP levels were slightly higher among women who subsequently used ABC (P = 0.04). Substantial mean increases in hsCRP and reductions in D-dimer levels were seen when comparing pre-HAART and on-HAART time points, for both the overall group (28 and -27%) and for MACS men (28 and -31%) and WIHS women (29 and -24%) separately (P < 0.01 for all). IL-6 levels, however, declined in MACS men but not in WIHS women, with the latter demonstrating low IL-6 levels at both visits (geometric mean was 1.78 pg/ml pre-HAART and 1.89 pg/ml on-HAART).
These trends were seen in both ABC and non-ABC HAART users. As described above, we conducted two types of analyses to compare these groups. First, we compared mean biomarker levels of persons receiving ABC HAART with levels from the corresponding matched visit of persons receiving non-ABC HAART. Figure 2 shows P values from univariate (PU) and multivariate (PM) random-effects models for comparisons of biomarker levels at the on-HAART visit between ABC and non-ABC users. The only significant adjusted difference was for D-dimer among MACS men, for whom ABC users had lower levels than non-ABC users (PM = 0.05).
Our second analysis evaluated individual biomarker level changes between pre-HAART and on-HAART time points. Table 3 depicts changes from the pre-HAART to on-HAART visits as percentage differences - [(on-HAART level - pre-HAART level)/pre-HAART level] x 100, that is the percentage increase or decrease in marker levels relative to the pre-HAART level as well as adjusted mean differences between non-ABC and ABC HAART users from random-effects models. Table 3a shows results from all 508 matched pairs and Table 3b shows results among a subgroup of 184 individuals who were ART-naive at HAART initiation (not matched). Among all pairs, there were no significant differences in pre-HAART to on-HAART biomarker changes between exposure groups after adjustment. Results were similar for each cohort separately, except that within the MACS there was a marginally greater reduction in D-dimer levels among ABC HAART initiators compared with non-ABC HAART initiators (adjusted difference in change of -24%, P = 0.06).
We evaluated four potential modifiers of the association between ABC use and biomarker levels: HIV RNA levels, ART experience at HAART initiation, HCV antibody status, and BMI. Only ART experience status at the time of HAART initiation was a significant modifier. Although this group is a minority in the MACS and WIHS, the vast majority of patients in the current treatment era are ART naive at HAART initiation. The first and second columns of Table 3b show how observed changes from the pre-HAART visit differ for non-ABC and ABC HAART users in the ART-naive subgroup compared with changes for the complete matched sample shown in the first and second columns of Table 3a. As in the complete matched sample, the estimated adjusted differences of change do not differ significantly between non-ABC and ABC HAART users in the ART-naive subgroup. However, in this subgroup, there was a trend toward greater reductions in biomarkers at the on-HAART visit among ABC HAART recipients (Table 3b column 3 compared with Table 3a column 3). We evaluated whether ABC recipients were more likely to ever have had a history of renal dysfunction [estimated glomerular filtration rate (GFR) <60 ml/min/1.73 m2] using the Modified Diet in Renal Disease equation. Among WIHS women there was no difference in history of renal dysfunction when stratifying by ABC exposure (data not shown). Among MACS men, there was a higher proportion of ABC users with a history of GFR less than 60 ml/min/1.73 m2 than that of non-ABC users (9% for ABC users vs. 4% for non-ABC users, P = 0.04). However, if persons who had only one GFR estimate less than 60 ml/min/1.73 m2 (with all other GFR estimates >60 ml/min/1.73 m2) were re-categorized as not having renal dysfunction, differences in the proportion of MACS men with renal dysfunction in the ABC vs. non-ABC groups disappeared (3 vs. 2%, respectively, P = 0.46).
Prior to HAART initiation, WIHS participants had lower hsCRP and IL-6 levels than MACS participants (1.39 vs. 1.88, P < 0.01, 1.88 vs. 3.75, P < 0.01) but similar D-dimer levels (0.46 vs. 0.41, P = 0.10). At the pre-HAART visit, there were six WIHS women who were pregnant, one (0.3%) in the non-ABC group and five (2%) in the ABC group (P = 0.10). At the on-HAART visit, there were 16 women who were pregnant, seven (2%) in the non-ABC group and nine (3%) in the ABC group (P = 0.63). No differences in any biomarkers levels were seen at either visit when comparing pregnant and nonpregnant women.
Finally, we examined the association of change in plasma HIV RNA levels between the two visits (mean change -1.4 log10 copies/ml) and changes in biomarker levels using Spearman's correlation coefficients. Greater reductions in plasma HIV RNA were significantly correlated with reductions in D-dimer (r = 0.14, P < 0.01) and IL-6 (r = 0.12, P < 0.01) but not with changes in hsCRP (r = -0.01, P = 0.70). Results were similar when restricted to persons with HIV RNA values above the limit of detection (80 copies/ml) at the pre-HAART visit.
Selection of study participants
The MACS and WIHS cohorts have been previously described [13,14]. Briefly, both studies prospectively schedule participants for semi-annual study visits where medical and therapy histories, behavioral assessments, and biological specimens are collected as part of standardized protocols. Individuals eligible for the current study were HIV-seropositive participants who initiated HAART during study follow-up. HAART initiators were restricted to persons for whom 1 year or less had elapsed between the last date at which they reported no HAART use and the date at which first HAART use was reported. Persons who reported use of any non-HAART ABC regimen prior to first HAART use, and person-visits without specimens available, were excluded. HAART was defined in accordance with US Department of Health and Human Services treatment guidelines (www.aidsinfo.gov) as use of at least two NRTIs in combination with at least one protease inhibitor or one nonnucleoside reverse transcriptase inhibitor (NNRTI); one NRTI in combination with at least one protease inhibitor and at least one NNRTI; or an ABC-containing or tenofovir-containing regimen with at least three NRTIs in the absence of both protease inhibitor and NNRTI, except for the three NRTI regimens consisting of ABC/tenofovir/lamivudine or didanosine/tenofovir/lamivudine. Antiretroviral (ART) regimens reported at each semi-annual visit for HAART initiators were categorized as either HAART including ABC (ABC HAART), HAART not including ABC (non-ABC HAART), or non-HAART.
Figure 1 depicts the study design. There were 1508 individuals (567 MACS, 941 WIHS) who initiated HAART during study follow-up and contributed 11 703 HAART person-visits (5136 MACS, 6567 WIHS). From these, we identified 508 'ABC-exposed' individuals (180 MACS, 328 WIHS) who reported ABC use.
To identify an appropriate comparison group of ABC-unexposed individuals, we utilized propensity score methodology . Propensity score methods extend traditional matching methods to eliminate imbalances, and thus confounding, between exposed and unexposed groups in a large set of measured variables. Briefly, we took into account all HAART study visits from HAART initiation to the first report of ABC use or to the last study visit if ABC use was not reported. As visits occurred semi-annually, ABC initiation could have taken place less than 1 to 6 months before the visit at which ABC use was first reported. Propensity scores for each person-visit at which HAART use was reported were obtained from logistic regression models where the outcome was ABC HAART initiation. Predictors included age, calendar year, cigarette use, body mass index (BMI), high-density lipoprotein (HDL) cholesterol level, race/ethnicity, hepatitis C virus (HCV) antibody status, ART experience at HAART initiation, duration of HAART use, the proportion of visits between HAART initiation and the last eligible visit where HAART was not reported, CD4+ cell count, and plasma HIV RNA levels. Propensity score matching was undertaken separately within each cohort (MACS men and WIHS women). For each ABC initiator, we identified a non-ABC HAART user by selecting the participant visit with a propensity score value closest to and no more than 3% different than the propensity score at the visit at which ABC initiation was first reported. Whereas the visits of ABC users prior to ABC initiation were included in model construction so the same individual could contribute to both exposure groups, each matched pair consisted of two individuals, one exposed and one unexposed. The last available visit within 1 year prior to first report of HAART was defined as the baseline (pre-HAART) visit. In summary, our selection identified two visits for each ABC initiator (pre-HAART and on-HAART) matched to two visits of a non-ABC HAART user within the same cohort.
Plasma levels of hsCRP, D-dimer, and IL-6 were measured for each propensity score-matched patient pair at the pre-HAART and on-HAART visits using stored specimens collected in CPT heparin tubes for MACS men and CPT citrate tubes for WIHS women.
hsCRP, IL-6, and D-dimer were chosen as biomarkers because they showed important associations with mortality in the SMART study and appeared to explain much of the association of ABC with CVD events [2-4,15,16]. hsCRP is a well established CVD risk factor in general populations, adding risk information to traditional risk factors such as lipids and blood pressure [15,17]. Levels were measured at the University of Vermont Laboratory for Clinical Biochemistry Research (LCBR) using a nephelometric immunoassay (Dade-Behring BN II, Deerfield, Illinois, USA), with an analytical coefficient of variance of les s than 6%. Consensus cut points for CVD risk are less than 1.0 mg/l, low risk; 1.0-3.0 mg/l, intermediate risk; more than 3.0 mg/l, high risk [15,17]. IL-6 is a proinflammatory cytokine with many important activities related to atherosclerosis and CVD . Levels have been used to explore CVD risk, and IL-6 is a particularly strong risk factor for mortality [18,19]. It was measured by chemiluminescent ELISA (R&D Systems, Minneapolis, MN) with a coefficient of variance of less than 14%. The LCBR multiethnic reference range is 0.5-3.0 pg/ml. D-Dimer, which results from the plasmin-mediated lysis of fibrin formed during ongoing coagulation, is an integrated marker of ongoing coagulation and fibrinolysis and has also been shown to be a CVD risk factor [19,20]. In addition, D-dimer is a risk factor for venous thrombosis . It was measured using an immunoturbidometric assay (Stago STA-R, Parsippany, NJ), with a coefficient of variance of less than 15%. The LCBR multiethnic reference range is 0.07-0.7 µg/ml.
Summary statistics for markers measured at pre-HAART and on-HAART (ABC and non-ABC) visits were computed for all matched pairs and also separately for matched pairs within each cohort. Marker measurements showed considerable skewness, necessitating the use of log tranformations. Random-effects models were used to compare log-transformed plasma biomarker levels at the pre-HAART and on-HAART visits and to evaluate changes in these levels among ABC initiators as compared with non-ABC users . Capitalizing on the repeated measurements for each individual, change was defined as the log-transformed ratio of the biomarker level at the on-HAART visit to the level at the pre-HAART visit - this is equivalent to the difference in log-transformed marker levels at the two visits. Models included random-effect terms for each pair, were weighted by the inverse of the propensity scores [22,23], and were adjusted for age, calendar year of the index visit, time elapsed between baseline and index visits, smoking status, BMI, HDL cholesterol levels, race/ethnicity, HCV antibody status, ART experience at HAART initiation, CD4 cell count, HIV RNA level, cohort, and ABC exposure. Likelihood ratio tests (LRTs) were used to evaluate whether or not each individual factor contributed significantly to the model fit or not. The LRT method was also used to evaluate specific prespecified interactions between ABC exposure and HIV RNA levels, ART experience at HAART initiation, HCV antibody status, and BMI. Lastly, we also examined the association between the change in biomarker levels and the change in log10 HIV RNA levels between the pre-HAART and on-HAART visits using Spearman's correlation coefficients. Analyses were conducted using all data and repeated among those with detectable (>80 copies/ml) values at the pre-HAART visit. All analyses were conducted in SAS version 9 (SAS, Cary, North Carolina, USA).
The MACS and WIHS protocols are reviewed annually by each participating institution's institutional review board.