- Research article
- Open Open Peer Review
Risk of cardiovascular events from current, recent, and cumulative exposure to abacavir among persons living with HIV who were receiving antiretroviral therapy in the United States: a cohort study
BMC Infectious Diseasesvolume 17, Article number: 708 (2017)
There is ongoing controversy regarding abacavir use in the treatment of HIV infection and the risk of subsequent development of cardiovascular disease. It is unclear how the risk varies as exposure accumulates.
Using an administrative health-plan dataset, risk of cardiovascular disease events (CVDe), defined as the first episode of an acute myocardial infarction or a coronary intervention procedure, associated with abacavir exposure was assessed among HIV-infected individuals receiving antiretroviral therapy across the U.S. from October 2009 through December 2014. The data were longitudinal, and analyzed using marginal structural models.
Over 114,470 person-years (n = 72,733) of ART exposure, 714 CVDe occurred at an incidence rate (IR) (95% CI) of 6·23 (5·80, 6·71)/1000 person-years. Individuals exposed to abacavir had a higher IR of CVDe of 9·74 (8·24, 11·52)/1000 person-years as compared to 5·75 (5·30, 6·24)/1000 person-years for those exposed to other antiretroviral agents. The hazard (HR; 95% CI) of CVDe was increased for current (1·43; 1·18, 1·73), recent (1·41; 1·16, 1·70), and cumulative [(1·18; 1·06, 1·31) per year] exposure to abacavir. The risk for cumulative exposure followed a bell-shaped dose-response curve peaking at 24-months of exposure. Risk was similarly elevated among participants free of pre-existing heart disease or history of illicit substance use at baseline.
Current, recent, and cumulative use of abacavir was associated with an increased risk of CVDe. The findings were consistent irrespective of underlying cardiovascular risk factors.
Cardiovascular disease (CVD) accounts for approximately 16% of deaths among persons living with HIV (PLWH) . Risk factors for CVD are more prevalent among PLWH , and use of various antiretroviral (ARV) drugs has been shown to be associated with an increased risk of CVD . With rapid expansion of antiretroviral therapy (ART) coverage both domestically and abroad, researchers and clinicians have become increasingly aware of potential ARV drug-related adverse events. Whether the commonly used ARV drug abacavir is associated with an increased risk of CVD has been intensely debated. Abacavir sulfate is a guanosine analog nucleoside reverse transcriptase inhibitor that possesses retroviral suppressive properties similar to tenofovir , and is a commonly prescribed backbone ARV agent. However, the writing of prescriptions of abacavir declined after the Data Collection on Adverse Events of Anti-HIV Drugs (D:A:D) study group reported in 2008 an increased risk of acute myocardial infarction (AMI) among PLWH exposed to abacavir [5,6,7]. Independent investigations that were subsequently carried out have both supported [7,8,9,10,11,12,13,14,15,16,17] and refuted [18,19,20,21,22,23] the D:A:D study group’s findings.
While studies conducted more recently have generally suggested an increased risk of CVD from abacavir exposure [8, 10, 14, 17], they were limited by few outcomes, with results occasionally underpowered [8, 17]. Failure to identify a clear underlying biological mechanism to explain the epidemiologic findings has added to the deliberation . Furthermore, there has also been a lack of consensus regarding whether the risk of CVD from exposure to abacavir reverses within a few months of stopping the drug [5, 16] and a lack of understanding on how the risk varies as exposure accumulates. In this study, we have sought to address these limitations by investigating the risk of CVD events (CVDe) from current, recent, and cumulative exposure to abacavir among PLWH using conventional and causal statistical methods.
Study design, sample collection and participants
The risk of CVDe was assessed among PLWH who started ARV drugs in the U.S. between October 1, 2009 and December 31, 2014. Data were obtained from medical and prescription claims data included in the IMS’ PharMetrics Plus database. October 1, 2009 was the earliest possible date for complete availability of relevant data; ART prescription history prior to this date was not available. PharMetrics Plus is a large health plan insurance claims database in the U.S., and is comprised of adjudicated claims for more than 150 million unique enrollees from across four regions of the U.S. . The data undergo a series of quality checks to minimize errors. This study used a pre-defined algorithm (Fig. 1) to extract and define the study population of PLWH exposed to any ART in the database. The study population was restricted to those ≥18 years of age. The baseline time point was defined as the date of ART initiation in the database and individual follow up time was censored at the first of three events after baseline: 1) first occurrence of CVDe, 2) last recorded date of ART receipt, 3) December 31, 2014. The study was approved by the Committee for Protection of Human Subjects at the University of California, Berkeley.
Exposure, covariate, and outcome definitions
Exposures to specific ARV agents were identified by generic product identifier (GPI) codes. Person-time of exposure to abacavir was compared to exposure to ARV agents other than abacavir. Any two prescriptions for an ARV agent separated by <30 days were combined to represent a single continuous exposure; gaps ≥30 days were not combined and this person-time was not included in the analysis. These data are longitudinal, and each subject’s follow up time was divided into consecutive one-month periods during which treatment was allowed to vary. The values of covariates were updated at the start of each month. The outcome of CVDe for an individual was defined as the first occurrence of an AMI or receipt of a coronary artery intervention procedure (i.e. percutaneous coronary intervention or coronary artery bypass graft) after baseline. AMI and coronary artery intervention procedures were ascertained using the International Classification of Disease, 9th Revision, Clinical Modification (ICD-9-CM) or Current Procedural Terminology (CPT) codes, respectively (Additional file 1: Table S1). The ICD-9 code used for AMI (410.xx) has been previously validated in another claims database .
The temporal ordering of covariate, treatment, and outcome allowed for a time-varying analysis, and the opportunity for causal interpretations. The first observation of a time-dependent covariate corresponded to its baseline value and once a health condition developed, an individual was assumed to have the condition for the remainder of the study. Current exposure to abacavir was defined as exposure (yes/no) during each one-month observation period. Recent exposure was defined as exposure (yes/no) in the previous six months (inclusive of the current month). Cumulative exposure was defined as the total duration of exposure an individual had received at a particular time point in one-month increments, updated monthly. Duration of exposure ceased to accumulate upon discontinuation of the drug but resumed if the drug was restarted. HIV-infection status and covariates were ascertained using the ICD-9-CM or CPT codes (Additional file 1: Table S1).
The risk of CVDe from a current, recent, and cumulative exposure to abacavir was estimated by the parameters of pooled logistic regression marginal structural models using stabilized inverse probability of treatment weights (sIPTW) . The sIPTW was generated from four treatment models – two each for the numerator and the denominator of the weight . For the denominator, the time point specific probability of exposure initiation was first estimated by fitting a main term pooled logistic regression to data up to the individual’s first month of receiving the exposure or end of follow up for those who were never exposed. The probability of exposure continuation was then estimated by fitting the model to data after the first month of starting the exposure. The denominator was modelled as a function of baseline covariates: gender, tobacco use/smoking (ever), substance or alcohol abuse (ever), serologic evidence of hepatitis B and C infections, history of stroke, cancer or old myocardial infarction, and time-dependent covariates: age, year of ART initiation, body weight, receipt of hypoglycemic agents (i.e. sulfonylureas, biguanides, insulin, thiazolidinedione) or medications for CVDe (i.e. aspirin, beta-blocker, angiotensin converting enzyme inhibitor, angiotensin receptor blocker, calcium channel blocker, statins) or diagnoses of: chronic kidney disease (CKD), dyslipidemia, heart failure, cardiac dysrhythmia, atherosclerosis, diabetes mellitus, and hypertension. The exposure continuation model additionally contained a variable for past month’s exposure status. The probabilities for the numerator of the sIPTW were similarly modelled but as a function of baseline covariates only. The follow-up time was modeled as a function of natural cubic splines with three internal knots placed at 25th, 50th and 75th percentiles. The marginal structural model was adjusted for the sIPTW and the baseline covariates. Same treatment weights were used for estimation of CVDe risk from current, recent, and cumulative exposure to abacavir. In order to assess the change in risk over time, the adjusted and marginal models were fit as a function of categories of cumulative exposure, i.e., never exposed, 1–6, 7–12, 13–18, 19–24, and >25 months of exposure. In sensitivity analyses, the study population was restricted to individuals free of CVD at baseline, and to individuals without a history of alcohol and substance abuse at baseline. Sensitivity analyses were additionally carried out to assess if the risk of CVDe from abacavir exposure differed after adjusting for other anti-retroviral agents. Using the same sIPTW models, we tested for interaction to see whether risk of CVDe from current abacavir exposure is modified in the presence of 13 different risk factors (Additional file 1: Table S5). In addition to the marginal structural results, corresponding results from unadjusted and adjusted Cox models were calculated. This study assumes uninformative censoring. Data were extracted and processed from the main claims databases using TERADATA (Dayton, OH), SAS version 9.1 (SAS Institute, Cary, NC), and STATA version 13.1 (StataCorp, College Station, TX). The marginal structural models were implemented in STATA version 13.1, based on Fewell et al. . The rationale, definition, and implementation of the marginal structural models are described in Additional file 1: Appendix 1.
Study population and incidence rates
There were 72,733 participants contributing 114,470 person-years of exposure to antiretroviral agents. On average, participants were exposed to ART for 1.5 years. The mean age of the study population was 46 years and 82% were males. The characteristics of the study population at baseline and summary of exposure to various antiretroviral drugs are described in Tables 1, 2 respectively. Overall, 714 CVDe occurred at an incidence rate of 6.23 (95% CI: 5.80, 6.71)/1000 person-years. Of the 714 outcomes, 137 were observed over 14,060 person-years of current exposure to abacavir at an incidence rate of 9.74 (95% CI: 8.24, 11.52)/1000 person-years, as compared to 577 outcomes over 100,410 person-years with an incidence rate of 5.75 (95% CI: 5.30, 6.24)/1000 person-years for those currently exposed to other ARV drugs. The incidence rate was highest for those exposed to abacavir between 13 and 18 months (11.32/1000 person-years) (Table 3). Of the 714 CVDe, 548 were cases of AMI. The overall incidence rate of AMI was 4.78 (95% CI: 4.39, 5.19)/1000 person-years (Additional file 1: Table S2). We calculated a population attributable risk (PAR) associated with abacavir exposure as:
Factors associated with Abacavir use
At baseline, abacavir recipients had a higher prevalence of essential hypertension, diabetes mellitus, chronic kidney disease (CKD), dyslipidemia, lipodystrophy, heart disease, and use of cardiovascular medications (Table 1). In the pooled logistic regression model, older age, a diagnosis of CKD, symptomatic HIV infection, and presence of lipodystrophy were associated with an increased probability of receiving abacavir (Additional file 1: Table S3).
Predictors of outcome
The sIPTW models showed the risk of CVDe (HR; 95% CI) was increased for current (1.43; 1.18, 1.73), recent (1.41; 1.16, 1.70) and cumulative (1.18; 1.06, 1.31) exposure (per year) to abacavir (Table 4). Separate models were run for each of current, recent, and cumulative exposure. The unadjusted and adjusted Cox models also showed increased risk for these exposures (Table 4). On further assessment of the risk from cumulative exposure, the HR varied with the duration of exposure in an inverted U-shaped pattern (Table 5 and Fig. 2 ); the relative hazard continued to increase up to 24 months of exposure, after which it decreased to non-significant levels but remained elevated compared to those never exposed to abacavir. Older age, male sex, tobacco use, other heart diseases, prior AMI, use of CVD-related medications, diabetes mellitus, and dyslipidemia were each associated with increased hazard of CVDe in the adjusted Cox model (Additional file 1: Table S4). We also assessed whether the risk was reversible after six months of stopping abacavir by comparing those with any abacavir exposure prior to but not in the last six months including the current month to those never exposed and found that the risk (HR; 95% CI) remained elevated (sIPTW model: 1.69; 0.89, 3.20; adjusted Cox model: 2.08; 1.17, 3.71). In tests of interactions, we observed that the risk of CVDe associated with abacavir use was more pronounced for age < 45 years (interaction p-value: 0.028) and for people without prior heart disease (interaction p-value: 0.016) (Additional file 1: Table S5).
In a sensitivity analysis, we observed a 53% higher risk of CVDe (sIPTW model) for current exposure to abacavir among individuals without a prior AMI or heart disease at baseline (Additional file 1: Table S6). This relationship was also assessed by excluding other heart diseases (heart failure, cardiac arrhythmia, atherosclerosis, or receipt of cardiovascular medications) from the adjustment set of covariates for both the marginal and the adjusted Cox model, with similar results. The risk also remained elevated by 41% when the study population was restricted to individuals not using illicit substances or alcohol at baseline (Additional file 1: Table S6). We further tested for CVDe risk from abacavir use after adjusting for cumulative exposure to other antiretroviral agents (tenofovir, emtricitabine, zidovudine, lamivudine, lopinavir, atazanavir, darunavir, efavirenz, nevirapine, rilpivirine, and raltegravir) using sIPTW models and found elevated risk (HR; 95% CI) for current (1.38; 1.12, 1.68), recent (1.34; 1.09, 1.64), and cumulative exposure (1.16; 1.03, 1.31). We then replicated the D:A:D model  for cumulative exposure by including recent exposure in the same model as cumulative exposure and observed that although our hazard ratio estimate for risk from cumulative exposure (per year) remained elevated (HR:1.08; 95% CI: 0.89, 1.30), it decreased to a non-statistically significant level. When we modelled the risk by partitioning the cumulative exposure into various durations, we observed a similar increased risk [HR (95% CI)] pattern as observed in our primary analysis (Table 4): 1–6 months: 1.91 (0.95–3.83); 7–12 months: 2.58 (1.16–5.71); 13–18 months: 2.68 (1.17–6.11); 19–24 months: 2.90 (1.37–6.17); and ≥25 months: 2.13 (0.93–4.88).
In a large database claims-based study, we found an increased risk of CVDe associated with exposure to current, recent, and cumulative exposure to abacavir using both adjusted Cox and marginal structural models estimated with inverse probability treatment weights. The overall incidence rate of AMI in this study was 4.78/1000 person-years, which compares to 3.3/1000 person-years in the 2008 D:A:D study. AMI incidences of 1.41/1000 people and 1.2/1000 people were seen in the general population in Olmstead county in Minnesota in 2006 and in men 35–65 years of age in the Framingham study population, respectively [29, 30]. This relatively higher incidence of AMI in the PLWH could be due to HIV infection [31,32,33], ART use , or both; PLWH have been shown to have more risk factors for CVD as compared to the general population [31,32,33]. The incidence rates of AMI associated with exposure to abacavir in this study (6.9/1000 person-years) and in the D:A:D study (6.1/1000 person-years) were ~4–5 fold higher than the general population estimates and approximately 2-fold higher than in the general population of PLWH [2, 31,32,33,34]. Some of the difference in results between this study and the D:A:D study including higher incidence rate of AMI in this study could be because participants in this study were all exposed to ART whereas the D:A:D study included individuals who had not yet started ART, as well as those who had discontinued ART totally. We calculated a population attributable risk of 8%. This means 8% (n = 57) of the total CVDe risk in the PLWH could be prevented if abacavir was not used, assuming a causal relationship between abacavir use and CVDe risk.
In an attempt to characterize an underlying biological mechanism for the increase in CVDe risk associated with abacavir use, we assessed how the risk varied with duration of exposure. The relative hazard of AMI increased with increasing duration of exposure in an inverted U-shaped pattern, peaking between 13 and 24 months of exposure and leveling off thereafter, suggesting a dose response relationship between cumulative time exposed to abacavir and risk, up to 24 months. This result agrees with earlier finding by Young et al. in which they first showed that the risk of CVD increased with increasing duration of exposure, with greatest risk between 6 and 36 months and exposure beyond 36 months adding little to the risk, suggesting a dose-response pattern. We observed the dose-response relationship for various durations of cumulative exposure after controlling for recent exposure as well in addition to other variables in the model; the D:A:D study group  had reported that the observed risk for cumulative exposure disappeared after adjusting for recent exposure, meaning that the CVD risk existed only up to first 6 months of exposure, after which the risk reversed. In a separate model, we tested the risk reversibility and found a 69% increased risk of CVDe among those who had stopped abacavir prior to last six months, suggesting a risk not reversible within six months of stopping the drug. We did not formally test whether the inverted U-shaped curve described for cumulative exposure provides a better fit to the observed risk estimates than a simple linear association.
Whereas this and Young et al.’s study results do not support an underlying mechanism related to immediate exposure to abacavir, the results are not consistent with an atherogenic mechanism, in which an ongoing or increasing risk would be expected with an increasing duration of exposure, without necessarily reaching a peak effect and leveling off after 24 months. The finding of an early peak in the increased risk of AMI in the course of abacavir treatment is helpful in understanding how risk may change with continuing versus changing therapy. The study results presented here suggest a reversible but more gradual underlying mechanism with a longer lasting impact that regresses more slowly after removal of the exposure.
Prior work has suggested that abacavir-induced platelet hyper-reactivity and aggregation could potentially lead to thrombosis and myocardial infarction [35,36,37]. Specifically, abacavir may induce platelet hyper-reactivity by competitive inhibition of a nitric oxide-induced soluble guanylyl cyclase via its active metabolite, carbovir-triphosphate, leading to a decreased production of cyclic guanosine monophosphate, an inhibitor of platelet aggregation and secretion [24, 35, 36, 38]. It is possible that abacavir may trigger an acute platelet response leading to endothelial injury with a longer lasting impact. It is also unclear whether abacavir may exert its effect on CVD risk through an increase in inflammatory biomarkers. While the SMART/INSIGHT study investigators , Kristoffersen et al. , and Hileman et al.  showed evidence for a possible role of inflammatory biomarkers in causing CVD among abacavir users [e.g. increased levels of high sensitivity c-reactive protein (hsCRP) and interleukin-6], several other studies have shown that levels of inflammatory biomarkers such as hsCRP, interleukin-6, selectin P and E, D-dimer, vascular adhesion molecule-1, intercellular adhesion molecule-1, and tumor necrosis factor alpha are not elevated after exposure to abacavir [41,42,43,44,45,46,47,48,49,50,51,52,53,54]. Future interdisciplinary studies may explore these areas by bridging basic, translational and clinical science to provide additional insights into the mechanisms underlying abacavir-associated cardiovascular risk. We have not established a clear reason for observing a higher risk of CVDe associated with abacavir use among the younger age-group and individuals without a pre-existing cardiac condition in the test of interactions. While we acknowledge the exploratory nature of the analyses for interaction testing with the possibility that the results could be due to chance, the observation of a higher CVDe risk in individuals without prior heart disease may stand to support the finding of an increased risk in younger age people. The increased CVDe risk in younger age people could also reflect a higher prevalence of cocaine and injection drug use among them . It would be important to test in other populations whether CVD risk associated with abacavir use differs by age.
We used the sIPTW approach because individuals with certain risk factors for CVD such as CKD, hypertension, diabetes mellitus, and dyslipidemia, may be preferentially channeled into (or away from) receiving abacavir based on its known toxicity in the presence of these conditions. The sIPTW approach may also be necessary because post-baseline values of these variables may simultaneously serve as confounders and causal intermediates; adjusting for these through traditional methods can lead to biased results . Under such settings, the use of inverse probability weights provides a valuable tool for balancing confounders across exposure groups without conditioning on variables affected by treatment [55, 56]. Some of our results for various durations of cumulative exposure appreciably differed between conventional Cox models and marginal structural models. For example, the hazard ratios (95% CI, p value) for 7–12 months, 13–18 months, and 19–24 months of cumulative exposure were 1.27 (0.87–1.86; p = 0.219), 1.71 (1.10, 2.65; p = 0.016), and 1.62 (0.98, 2.69; p = 0.060), respectively, in adjusted Cox models. The corresponding hazard ratios from marginal structural models were 1.41 (0.97, 2.06; p = 0.073), 1.78 (1.16, 2.72; p = 0.009), and 1.90 (1.16, 3.11; p = 0.011) (Table 4).
A key strength of this study is the application of conventional and robust methods to address key study questions while using a very large U.S. health-plan dataset containing longitudinal information on usage of ART in >70,000 PLWH receiving care across the U.S. The recency of the data is an asset. Most studies that showed an association between abacavir use and CVD risk so far were hospital based [5, 9, 10, 14, 16,17,18,19] and hence may be subject to similar bias, such as channeling bias, that could arise from specific prescription behavior of physicians. Therefore, reproduction of the results in another representative population, such as that enrolled in the claims database, would be relevant and important. The similarity of these results to those from prior studies, the reproducibility of the results in the sensitivity analyses, and the finding of a background incidence rate of AMI comparable to that found in other studies are reassuring. A limitation of the study is that the ICD-9 and CPT diagnostic codes used may be prone to coding errors; however, such errors are likely to affect the exposure groups non-differentially and may not bias the study results. It is possible that information on covariates, such as body-weight, for which re-imbursement may not be sought could be under-reported in the database. Again, we expect this problem to exist non-differentially across exposure groups. This is an observational cohort study and is therefore subject to confounding from unmeasured factors and possible channeling bias; we have attempted to account for the latter by adopting an sIPTW-based analytic approach. Covariates that could be relevant but not available in the claims database and hence missing in our study are race/ethnicity, CD4 cell count, and HIV viral load. Adjustment for CD4 cell count and HIV viral load made little difference to the relative rate of AMI in a prior study . There is potential for bias in the study results from residual confounding that may arise from the binary categorization of certain variables in the study, rather than having a graded continuous response. We assumed uninformative censoring for the study because participants in both the exposure groups, i.e., PLWH receiving abacavir based ART regimen and PLWH receiving non-abacavir based ART regimen, may be at similar risk of adverse HIV–related life events that may cease their continued enrollment into the health-plan and hence representation in the database. We chose AMI and/or coronary artery interventions only to define CVDe so as to be as specific as possible with the study outcome’s representation of ischemic CVD; however, we might use a broader definition including other cardiac conditions or cerebrovascular events for the study outcome.
In summary, exposure to abacavir is associated with an increased risk of CVDe. We recommend a careful consideration of the risks and benefits of abacavir treatment while formulating antiretroviral treatment regimens with patients.
Acute Myocardial Infarction
Chronic Kidney Disease
Current Procedural Terminology
Cardiovascular Disease Events
Data Collection on Adverse Events of Anti-HIV Drugs
Human Immune Deficiency Virus
High Sensitivity C-Reactive Protein
International Classification of Disease, 9th Revision, Clinical Modification
People Living With HIV
Stabilized Inverse Probability Treatment Weights
Antiretroviral Therapy Cohort C. Causes of death in HIV-1-infected patients treated with antiretroviral therapy, 1996-2006: collaborative analysis of 13 HIV cohort studies. Clin Infect Dis. 2010;50(10):1387–96.
Triant VA, Lee H, Hadigan C, Grinspoon SK. Increased acute myocardial infarction rates and cardiovascular risk factors among patients with human immunodeficiency virus disease. J Clin Endocrinol Metab. 2007;92(7):2506–12.
Bavinger C, Bendavid E, Niehaus K, Olshen RA, Olkin I, Sundaram V, Wein N, Holodniy M, Hou N, Owens DK, et al. Risk of cardiovascular disease from antiretroviral therapy for HIV: a systematic review. PLoS One. 2013;8(3)
Cruciani M, Mengoli C, Malena M, Serpelloni G, Parisi SG, Moyle G, Bosco O. Virological efficacy of abacavir: systematic review and meta-analysis. J Antimicrob Chemother. 2014;69(12):3169–80.
Sabin CA, Worm SW, Weber R, Reiss P, El-Sadr W, Dabis F, De Wit S, Law M, Monforte AD, Friis-Moller N, et al. Use of nucleoside reverse transcriptase inhibitors and risk of myocardial infarction in HIV-infected patients enrolled in the D : a : D study: a multi-cohort collaboration. Lancet. 2008;371(9622):1417–26.
Antoniou T, Gillis J, Loutfy MR, Cooper C, Hogg RS, Klein MB, Machouf N, Montaner JSG, Rourke SB, Tsoukas C, et al. Impact of the data collection on adverse events of anti-HIV drugs cohort study on abacavir prescription among treatment-naive, HIV-infected patients in Canada. J Int Assoc Providers AIDS Care. 2014;13(2):153–9.
Sabin CA, Reiss P, Ryom L, Phillips AN, Weber R, Law M, Fontas E, Mocroft A, de Wit S, Smith C, et al. Is there continued evidence for an association between abacavir usage and myocardial infarction risk in individuals with HIV? A cohort collaboration. BMC Med. 2016;14:61.
Brouwer ES, Napravnik S, Eron Jr JJ, Stalzer B, Floris-Moore M, Simpson Jr RJ, Stürmer T. Effects of combination antiretroviral therapies on the risk of myocardial infarction among HIV patients. Epidemiology. 2014;25(3):406–17.
Choi AI, Vittinghoff E, Deeks SG, Weekley CC, Li YM, Shlipak MG. Cardiovascular risks associated with abacavir and tenofovir exposure in HIV-infected persons. AIDS. 2011;25(10):1289–98.
Desai M, Joyce V, Bendavid E, Olshen RA, Hlatky M, Chow A, Holodniy M, Barnett P, Owens DK. Risk of cardiovascular events associated with current exposure to HIV antiretroviral therapies in a US veteran population. Clin Infect Dis. 2015;61(3):445–52.
Durand M, Sheehy O, Lelorier J, Tremblay CL. Association between use of antiretroviral therapy and risk of acute myocardial infarction: a nested case control study using Quebec's public health insurance database (RAMQ). J Popul Ther Clin Pharmacol. 2011;18(2):e178–9.
Martin A, Bloch M, Amin J, Baker D, Cooper DA, Emery S, Carr A, Grp SS. Simplification of antiretroviral therapy with Tenofovir-Emtricitabine or Abacavir-lamivudine: a randomized, 96-week trial. Clin Infect Dis. 2009;49(10):1591–601.
Obel N, Farkas DK, Kronborg G, Larsen CS, Pedersen G, Riis A, Pedersen C, Gerstoft J, Sørensen HT. Abacavir and risk of myocardial infarction in HIV-infected patients on highly active antiretroviral therapy: a population-based nationwide cohort study. HIV Med. 2010;11(2):130–6.
Palella FJ Jr, Althoff K, Moore R, Zhang J, Kitahata M, Gange S, Crane H, Drozd D, Brooks J, Elion R. Abacavir use and risk for myocardial infarction in the NA-ACCORD. In: CROI 2015: February 23–26, 2015 2015; Seattle; 2015.
Strategies for Management of Anti-Retroviral Therapy I, Groups DADS. Use of nucleoside reverse transcriptase inhibitors and risk of myocardial infarction in HIV-infected patients. AIDS (London). England. 2008;22(14):F17–24.
Young J, Xiao Y, Moodie EEM, Abrahamowicz M, Klein MB, Bernasconi E, Schmid P, Calmy A, Cavassini M, Cusini A, et al. Effect of cumulating exposure to Abacavir on the risk of cardiovascular disease events in patients from the Swiss HIV cohort study. J Acquir Immune Defic Syndr. 2015;69(4):413–21.
Marcus JL, Neugebauer RS, Leyden WA, Chao CR, Xu L, Quesenberry CP Jr, Klein DB, Towner WJ, Horberg MA, Silverberg MJ. Use of Abacavir and risk of cardiovascular disease among HIV-infected individuals. J Acquir Immune Defic Syndr. 2016;71(4):413–9.
Bedimo RJ, Westfall AO, Drechsler H, Vidiella G, Tebas P. Abacavir use and risk of acute myocardial infarction and cerebrovascular events in the highly active antiretroviral therapy era. Clin Infect Dis. 2011;53(1):84–91.
Lang S, Mary-Krause M, Cotte L, Gilquin J, Partisani M, Simon A, Boccara F, Costagliola D. Impact of individual antiretroviral drugs on the risk of myocardial infarction in human immunodeficiency virus-infected patients: a case-control study nested within the French hospital database on HIV ANRS cohort CO4. Arch Intern Med. 2010;170(14):1228–38.
Cruciani M, Zanichelli V, Serpelloni G, Bosco O, Malena M, Mazzi R, Mengoli C, Parisi SG, Moyle G. Abacavir use and cardiovascular disease events: a meta-analysis of published and unpublished data. AIDS. 2011;25(16):1993–2004.
Ding X, Andraca-Carrera E, Cooper C, Miele P, Kornegay C, Soukup M, Marcus KA. No association of abacavir use with myocardial infarction: findings of an FDA meta-analysis. J Acquir Immune Defic Syndr. 2012;61(4):441–7.
Ribaudo HJ, Benson CA, Zheng Y, Koletar SL, Collier AC, Lok JJ, Smurzynski M, Bosch RJ, Bastow B, Schouten JT. No risk of myocardial infarction associated with initial antiretroviral treatment containing abacavir: short and long-term results from ACTG A5001/ALLRT. Clin Infect Dis. 2011;52(7):929–40.
Brothers CH, Hernandez JE, Cutrell AG, Curtis L, Ait-Khaled M, Bowlin SJ, Hughes SH, Yeo JM, Lapierre DH. Risk of myocardial infarction and Abacavir therapy: no increased risk across 52 GlaxoSmithKline-sponsored clinical trials in adult subjects. J Acquir Immune Defic Syndr. 2009;51(1):20–8.
Gresele P, Falcinelli E, Momi S, Francisci D, Baldelli F. Highly active antiretroviral therapy-related mechanisms of endothelial and platelet function alterations. Rev Cardiovasc Med. 2014;15(Suppl 1):S9–20.
IMS’ PharMetrics Plus Data Dictionary [http://www.imshealth.com/en/thought-leadership/quintilesims-institute/research-support/research-support%E2%80%93data-and-information].
Kiyota Y, Schneeweiss S, Glynn RJ, Cannuscio CC, Avorn J, Solomon DH. Accuracy of Medicare claims-based diagnosis of acute myocardial infarction: estimating positive predictive value on the basis of review of hospital records. Am Heart J. 2004;148(1):99–104.
Cole SR, Hernan MA. Constructing inverse probability weights for marginal structural models. Am J Epidemiol. 2008;168(6):656–64.
Fewell Z, Hernan MA, Wolfe F, Tilling K, Choi H, Sterne JAC. Controlling for time-dependent confounding using marginal structural models. Stata J. 2004;4(4):402–20.
Parikh NI, Gona P, Larson MG, Fox CS, Benjamin EJ, Murabito JM, O'Donnell CJ, Vasan RS, Levy D. Long-term trends in myocardial infarction incidence and case fatality in the National Heart, Lung, and Blood Institute's Framingham heart study. Circulation. 2009;119(9):1203–10.
Roger VL, Weston SA, Gerber Y, Killian JM, Dunlay SM, Jaffe AS, Bell MR, Kors J, Yawn BP, Jacobsen SJ. Trends in incidence, severity, and outcome of hospitalized myocardial infarction. Circulation. 2010;121(7):863–9.
Currier JS, Taylor A, Boyd F, Dezii CM, Kawabata H, Burtcel B, Maa JF, Hodder S. Coronary heart disease in HIV-infected individuals. J Acquir Immune Defic Syndr. 2003;33(4):506–12.
Data Collection on Adverse Events of Anti HIVdSG, Smith C, Sabin CA, Lundgren JD, Thiebaut R, Weber R, Law M, Monforte A, Kirk O, Friis-Moller N, et al. Factors associated with specific causes of death amongst HIV-positive individuals in the D:a:D study. AIDS. 2010;24(10):1537–48.
Worm SW, De Wit S, Weber R, Sabin CA, Reiss P, El-Sadr W, Monforte AD, Kirk O, Fontas E, Dabis F, et al. Diabetes mellitus, preexisting coronary heart disease, and the risk of subsequent coronary heart disease events in patients infected with human immunodeficiency virus: the data collection on adverse events of anti-HIV drugs (D:a:D study). Circulation. 2009;119(6):805–11.
Freiberg MS, Chang CC, Kuller LH, Skanderson M, Lowy E, Kraemer KL, Butt AA, Bidwell Goetz M, Leaf D, Oursler KA, et al. HIV infection and the risk of acute myocardial infarction. JAMA Intern Med. 2013;173(8):614–22.
Baum PD, Sullam PM, Stoddart CA, McCune JM. Abacavir increases platelet reactivity via competitive inhibition of soluble guanylyl cyclase. AIDS (London, England). 2011;25(18):2243–8.
Falcinelli E, Francisci D, Belfiori B, Petito E, Guglielmini G, Malincarne L, Mezzasoma A, Sebastiano M, Conti V, Giannini S, et al. Vivo platelet activation and platelet hyperreactivity in abacavir-treated HIV-infected patients. Thromb Haemost. 2013;110(2):349–57.
Satchell CS, O'Halloran JA, Cotter AG, Peace AJ, O'Connor EF, Tedesco AF, Feeney ER, Lambert JS, Sheehan GJ, Kenny D, et al. Increased platelet reactivity in HIV-1-infected patients receiving Abacavir-containing antiretroviral therapy. J Infect Dis. 2011;204(8):1202–10.
Gresele P, Francisci D, Falcinelli E, Belfiori B, Petito E, Fierro T, Mezzasoma AM, Baldelli F. Possible role of platelet activation in the cardiovascular complications associated with HIV infection: differential effects of abacavir (ABC) vs. tenofovir (TDF). J Thromb Haemost. 2011;9:799.
Kristoffersen US, Kofoed K, Kronborg G, Benfield T, Kjaer A, Lebech AM. Changes in biomarkers of cardiovascular risk after a switch to abacavir in HIV-1-infected individuals receiving combination antiretroviral therapy. HIV Med. 2009;10(10):627–33.
Hileman CO, Wohl DA, Tisch DJ, Debanne SM, McComsey GA. Short communication initiation of an Abacavir-containing regimen in HIV-infected adults is associated with a smaller decrease in inflammation and endothelial activation markers compared to non-Abacavir-containing regimens. AIDS Res Hum Retrovir. 2012;28(12):1561–4.
Young B, Squires KE, Ross LL, Santiago L, Sloan LM, Zhao HH, Wine BC, Pakes GE, Margolis DA, Shaefer MS, et al. Inflammatory biomarker changes and their correlation with Framingham cardiovascular risk and lipid changes in antiretroviral-naive HIV-infected patients treated for 144 weeks with Abacavir/lamivudine/Atazanavir with or without ritonavir in ARIES. AIDS Res Hum Retrovir. 2013;29(2):350–8.
Martinez E, Larrousse M, Podzamczer D, Perez I, Gutierrez F, Lonca M, Barragan P, Deulofeu R, Casamitjana R, Mallolas J, et al. Abacavir-based therapy does not affect biological mechanisms associated with cardiovascular dysfunction. AIDS. 2010;24(3):F1–9.
Palella FJ, Gange SJ, Benning L, Jacobson L, Kaplan RC, Landay AL, Tracy RP, Elion R. Inflammatory biomarkers and abacavir use in the Women's interagency HIV study and the multicenter AIDS cohort study. AIDS. 2010;24(11):1657–65.
Padilla S, Masia M, Garcia N, Jarrin I, Tormo C, Gutierrez F. Early changes in inflammatory and pro-thrombotic biomarkers in patients initiating antiretroviral therapy with abacavir or tenofovir. BMC Infect Dis. 2011;11
De Luca A, Donati KD, Cozzi-Lepri A, Colafigli M, De Curtis A, Capobianchi MR, Antinori A, Giacometti A, Magnani G, Vullo V, et al. Exposure to Abacavir and biomarkers of cardiovascular disease in HIV-1-infected patients on suppressive antiretroviral therapy: a longitudinal study. Jaids-Journal of Acquired Immune Deficiency Syndromes. 2012;60(3):E98–E101.
Kim C, Gupta SK, Green L, Taylor BM, Deuter-Reinhard M, Desta Z, Clauss M. Abacavir, didanosine and tenofovir do not induce inflammatory, apoptotic or oxidative stress genes in coronary endothelial cells. Antivir Ther. 2011;16(8):1335–9.
Patel P, Bush T, Overton T, Baker J, Hammer J, Kojic E, Conley L, Henry K, Brooks JT. Effect of abacavir on acute changes in biomarkers associated with cardiovascular dysfunction. Antivir Ther. 2012;17(4):755–61.
Rasmussen TA, Tolstrup M, Melchjorsen J, Frederiksen CA, Nielsen US, Langdahl BL, Ostergaard L, Laursen AL. Evaluation of cardiovascular biomarkers in HIV-infected patients switching to abacavir or tenofovir based therapy. BMC Infect Dis. 2011;11
Wohl DA, Arnoczy G, Fichtenbaum CJ, Campbell T, Taiwo B, Hicks C, McComsey GA, Koletar S, Sax P, Tebas P, et al. Comparison of cardiovascular disease risk markers in HIV-infected patients receiving abacavir and tenofovir: the nucleoside inflammation, coagulation and endothelial function (NICE) study. Antivir Ther. 2014;19(2):141–7.
Hammond E, McKinnon E, Mallal S, Nolan D. Longitudinal evaluation of cardiovascular disease-associated biomarkers in relation to abacavir therapy. AIDS. 2008;22(18):2540–3.
O’Halloran J, Dunne E, Tinago W, Denieffe S, Kenny D, Mallon P. Effect of switch from Abacavir to Tenofovir DF on platelet function markers: a SWIFT trial substudy. In: CROI 2014: March 3–6, 2014 2014; Boston. 2014:484.
Martin A, Amin J, Cooper DA, Carr A, Kelleher AD, Bloch M, Baker D, Woolley I, Emery S. Group Ss: Abacavir does not affect circulating levels of inflammatory or coagulopathic biomarkers in suppressed HIV: a randomized clinical trial. AIDS. 2010;24(17):2657–63.
De Luca A, de Gaetano Donati K, Cozzi-Lepri A, Colafigli M, De Curtis A, Capobianchi MR, Antinori A, Giacometti A, Magnani G, Vullo V, et al. Exposure to abacavir and biomarkers of cardiovascular disease in HIV-1-infected patients on suppressive antiretroviral therapy: a longitudinal study. J Acquir Immune Defic Syndr. 2012;60(3):e98–101.
Patel P, Bush T, Overton T, Baker J, Hammer J, Kojic E, Conley L, Henry K, Brooks JT, Investigators SUNS. Effect of abacavir on acute changes in biomarkers associated with cardiovascular dysfunction. Antivir Ther. 2012;17(4):755–61.
Robins JM, Hernan MA, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology. 2000;11(5):550–60.
Petersen ML, van der Laan MJ. Causal models and learning from data: integrating causal modeling and statistical estimation. Epidemiology. 2014;25(3):418–26.
We would like to thankfully acknowledge the study participants and the clinic staff at the study centers.
Availability of data and materials
The datasets used and /or analyzed during the current study are available from the corresponding author on reasonable request and shall be made available at any time, upon request or suggestion by the editorial board, through an online public repository such as DYRAD https://datadryad.org/pages/organization, etc.
K.D. is currently a post-doctoral fellow at the Johns Hopkins School of Medicine Infectious Diseases. K.D. is working with Richard Chaisson, MD, on global control of tuberculosis and HIV-AIDS, especially in the pediatric population. K.D. recently completed his Ph.D. degree in Epidemiology from UC-Berkeley School of Public Health under supervision of Arthur Reingold, MD, and Alan Hubbard, PhD, who have overseen his doctoral dissertation work including the work contained in this manuscript. Dr. Reingold, co-author on this paper, is the head of UC-Berkeley School of Public Health Division of Epidemiology and Dr. Hubbard is the head of UC-Berkeley School of Public Health Division of Biostatistics. Dr. Baxi has been an infectious disease fellow at the University of California, San Francisco and a Ph.D. student at UC-Berkeley at the time of this work.
Ethics approval and consent to participate
The study was approved by the Committee for Protection of Human Subjects at the University of California, Berkeley.
Consent for publication
The authors declare that they have no competing interests.
Disclosures: K.D. was an intern with the Division of Epidemiology at Gilead Sciences (Foster City, CA, USA), which supported the acquisition of the data. Gilead Sciences was not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. This manuscript is a result of K.D.’s PhD dissertation work at UC-Berkeley.
Support: S.M.B. is supported by the UCSF Traineeship in AIDS Prevention Studies (US National Institutes of Health (NIH) T32 MH-19105). The funders were not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ICD-9-CM and CPT codes for defining various covariates and outcomes. Table S2. Age-specific incidence rate (IR) of acute myocardial infarction (AMI) among persons living with HIV receiving antiretroviral therapy. Table S3. Factors associated with initiation of abacavir among persons living with HIV, by pooled logistic regression. Table S4 The influence of various risk factors on the development of CVD among persons living with HIV receiving anti-retroviral therapy. Table S5 Risk of CVD from current exposure to abacavir in sub-groups of variables at baseline (test of interactions). Table S6 Risk of cardiovascular disease from exposure to abacavir among persons living with HIV free of heart diseasea or substance or alcohol abuse at baseline. Appendix 1. Detailed approach to developing marginal structural models. (DOCX 58 kb)