Tenofovir-Associated Bone Adverse Outcomes among a US National Historical Cohort of HIV-Infected Veterans: Risk Modification by Concomitant Antiretrovirals.

INTRODUCTION: Tenofovir disoproxil fumarate (TDF) has been associated with greater incidences of bone complications, which might be modified by some concomitantly administered antiretrovirals, possibly by their effect on tenofovir concentrations. We compared bone adverse outcomes among treatment-naïve HIV-infected US veterans initiating efavirenz (EFV)-containing TDF/emtricitabine (FTC) regimens versus those initiating non-EFV-containing TDF/FTC regimens.
METHODS: Using national Veterans Health Administration clinical and administrative data sets, we identified a cohort of treatment-naïve HIV-infected veterans without bone disease who initiated therapy with TDF/FTC plus EFV, rilpivirine, elvitegravir/cobicistat, or ritonavir-boosted protease inhibitors in 2003-2015. The primary composite adverse bone outcome was the unadjusted incidence rate (IR) of osteoporosis, osteopenia, or fragility fracture (any hip, wrist, or spine fracture). To account for selection bias and confounding, we used inverse probability of treatment-weighted Cox proportional hazards regression models to calculate adjusted hazard ratios (HRs) for each outcome associated with EFV + TDF/FTC versus each non-EFV-containing TDF/FTC regimen.
RESULTS: Of 33,048 HIV-positive veterans, 7161 initiated a TDF/FTC-containing regimen (mean age, 50 years; baseline CD4 < 200 cells/mm3, 33.3%; HIV-1 RNA > 100,000 copies/ml, 22.3%; mean follow-up, 13.0 months). Of these, 4137 initiated EFV- and 3024 non-EFV-containing regimens. Veterans initiating EFV- versus non-EFV-containing TDF/FTC regimens had a lower IR of the composite bone outcome (29.3 vs. 41.4 per 1000 patient-years), with significant risk reductions for this outcome [HR, 0.69; 95% confidence interval (CI), 0.58-0.83] and fragility fracture (HR, 0.59; 95% CI, 0.44-0.78).
CONCLUSION: EFV + TDF/FTC is associated with a lower risk of adverse bone outcomes compared with other TDF-containing regimens in the VHA. FUNDING: Bristol-Myers Squibb.

Introduction

Patients with HIV infection have higher rates of osteoporosis and fragility fractures than uninfected individuals, which are not completely explained by differences in traditional risk factors such as age and body mass index [1]. Randomized controlled trials and observational studies suggest that HIV infection itself and the initiation of some antiretroviral therapies (ART) may independently increase the risk of bone adverse outcomes [1-8]. For these reasons, the safety profiles of antiretroviral agents are subject to increased scrutiny in an effort to reduce treatment-related side effects in the aging HIV population.

Tenofovir disoproxil fumarate (TDF)/emtricitabine (FTC) has been the main nucleos(t)ide reverse transcriptase inhibitor backbone of combination ART for more than a decade; TDF/FTC continues to be included in four of the six recommended regimens for treatment-naïve patients with HIV infection in the current Department of Health and Human Services guidelines [9]. TDF has been noted to increase the risk of fragility fractures [7, 10]. In a Veterans Health Administration (VHA) population, cumulative exposure to TDF was independently predictive of fragility fracture [11]; however, no evaluation of TDF-associated fracture risk based on concomitantly administered ART was undertaken. Very few epidemiologic studies in US populations have compared risks of bone adverse outcomes with TDF across differing TDF-containing antiretroviral regimens [12].

ART initiation is associated with increases in bone turnover and modest decreases in bone mineral density (BMD) [13], which are greater with TDF-containing regimens. BMD loss with TDF is amplified when coadministered with boosted protease inhibitors (PIs), but effects are less certain when coadministered with other antiretroviral agents [7, 10, 14–25]. For example, in the AIDS Clinical Trials Group (ACTG) A5224s substudy of ACTG A5202, spine BMD loss was significantly greater in those receiving ritonavir-boosted atazanavir (ATV) + TDF/FTC (−3.1%) than in those receiving efavirenz (EFV) + TDF/FTC (−1.7%; P = 0.035) [7]. However, the clinical relevance of these observations (increased fracture risk) remains uncertain.

Regarding fracture risk, available data from randomized controlled trials are limited by the small sample sizes and short observation periods of individual trials and the lack of specific reporting of fragility fractures [26]. Thus, despite differences in BMD across regimens in ACTG A5224s, fractures (all traumatic) were uncommon, and rates were similar across regimens [7]. In contrast, one recent large cohort study of commercial claims data showed lower fracture incidence rates (IRs) for EFV + TDF/FTC compared with elvitegravir/cobicistat (EVG/c)/TDF/FTC, rilpivirine (RPV) + TDF/FTC, and the overall HIV population in the database [12].

The lower BMD loss and fracture risk with EFV-containing TDF/FTC regimens compared with other third agents combined with TDF/FTC may relate to differences in drug-drug and drug-food interactions across these regimens. Tenofovir plasma concentrations are increased by 22–37% when TDF is taken with PIs [27, 28], cobicistat [29], or RPV [30]. The bioavailability of tenofovir is also increased by up to 40% with concomitant food intake [28, 31]. Certain TDF regimens, including RPV, EVG/c, and boosted PIs, are taken with food [32-34]. Conversely, EFV + TDF/FTC is taken in the fasted state, and no clinically relevant drug-drug interactions have been reported between TDF and either EFV or FTC. Therefore, we hypothesized that EFV + TDF/FTC could result in a lower incidence of bone adverse outcomes compared with these other TDF-containing regimens.

The purpose of this analysis was to compare the incidence of bone adverse outcomes among treatment-naïve HIV-infected US veterans without bone disease who initiated different TDF-containing regimens.

Methods

Study Design and Data Sets

Using a national cohort of US veterans, we conducted a population-based historical cohort study using VHA databases containing clinical, pharmacy, and administrative data from more than 150 VHA hospitals and 850 outpatient clinics nationwide [35]. We obtained demographic, laboratory, diagnosis, and utilization data from the Veterans Affairs (VA) Corporate Data Warehouse (CDW), including medical SAS data sets for in- and outpatient encounter data, CDW’s raw pharmacy data, and Decision Support Systems and CDW’s laboratory data. To conduct the analyses, data sets were housed in the VA Informatics and Computing Infrastructure environment, which enables access to data and tools for reporting and analysis in a secure workspace to ensure veterans’ privacy and data security. This article does not contain any studies with human participants or animals performed by any of the authors. The University of Utah Institutional Review Board and the Salt Lake City VA Health Care System Office of Research and Development approved this study.

Patients

The cohort included all HIV-infected antiretroviral-naïve veterans without bone disease who initiated TDF/FTC plus a third antiretroviral agent of interest (see Exposures section) during the period 2003–2015. A validated algorithm [36] with a sensitivity of 86% and positive predictive value of 87% [37] was used to exclude patients with evidence of prior ART received outside of the VHA, which included the following criteria: exposure to any antiretroviral agents during a 1-year period before the index date (the pre-index period), patients whose index ART regimen was a “salvage” regimen (i.e., composed of both a PI and an NNRTI or composed of 5 or more agents), and patients whose HIV RNA levels before the index date were low enough (< 500 copies/ml) to suggest prior antiretroviral exposure. Veterans were identified for inclusion by an available index date for the first pharmacy fill for one of the third agents of interest and if they fulfilled the following criteria at the index date: (1) aged  ≥ 18 years; (2) at least 6 months of pre-index date VHA activity including in- or outpatient services; (3) no evidence of prior treatment including fills for antiretrovirals in the 6 months before the index date; (4) no evidence of prior bone disease (defined as diagnosis of osteoporosis/osteopenia by International Classification of Disease diagnosis codes, Current Procedural Terminology codes, or classification by bone mineral density test result). Codes used to identify patients with HIV infection are shown in Supplemental Digital Content, Table 1.

Table 1

Unadjusted baseline characteristics among HIV-infected veterans receiving initial ART with differing TDF/FTC-containing regimens

	Initial ART Containing TDF/FTC Plus
	EFV (n = 4137)	Non-EFVa (n = 3024)	EVG/c (n = 232)	RPV (n = 171)	RTV-boosted PI (n = 2621)
Demographics and physical
Age, years, mean ± SD	50 ± 10	49 ± 9.8	49 ± 13	47 ± 13	50 ± 9.3
Male	4002 (96.7)	2903 (96.0)	221 (95.3)	161 (94.2)	2521 (96.2)
Married	337 (8.1)	221 (7.3)	26 (11.2)	20 (11.7)	175 (6.7)
Race
White	1240 (30.0)	911 (30.1)	73 (31.5)	48 (28.1)	790 (30.1)
Black	2492 (60.2)	1762 (58.3)	124 (53.4)	104 (60.8)	1534 (58.5)
Hispanic	260 (6.3)	207 (6.8)	19 (8.2)	10 (5.8)	178 (6.8)
Asian	34 (0.8)	28 (0.9)	5 (2.2)	2 (1.2)	21 (0.8)
Other	22 (0.5)	28 (0.9)	2 (0.9)	2 (1.2)	24 (0.9)
Missing	89 (2.2)	88 (2.9)	9 (3.9)	5 (2.9)	74 (2.8)
BMI (kg/m2, mean ± SD)	26 ± 5.0	25 ± 5.1	26 ± 5.6	27 ± 5.3	25 ± 5.1
Pre-indexb† prognostic indices
CD4+ count, cells/mm3
< 200	1286 (31.1)	1103 (36.5)	60 (25.9)	23 (13.5)	1020 (38.9)
200–299	698 (16.9)	477 (15.8)	33 (14.2)	30 (17.5)	414 (15.8)
300–399	634 (15.3)	342 (11.3)	34 (14.7)	11 (6.4)	297 (11.3)
400–499	401 (9.7)	222 (7.3)	24 (10.3)	33 (19.3)	165 (6.3)
≥ 500	588 (14.2)	426 (14.1)	61 (26.3)	63 (36.8)	302 (11.5)
Missing	530 (12.8)	454 (15.0)	20 (8.6)	11 (6.4)	423 (16.1)
HIV viral load, copies/ml
< 10,000	1184 (28.6)	950 (31.4)	65 (28.0)	68 (39.8)	817 (31.2)
10,000–100,000	1434 (34.7)	940 (31.1)	78 (33.6)	67 (39.2)	795 (30.3)
> 100,000	913 (22.1)	686 (22.7)	61 (26.3)	15 (8.8)	610 (23.3)
Missing	606 (14.6)	448 (14.8)	28 (12.1)	21 (12.3)	399 (15.2)
Pre-indexb renal function
eGFR, ml/min/1.73 m2
≥ 90	2193 (53.0)	1552 (51.3)	143 (61.6)	108 (63.2)	1301 (49.6)
60–89	1360 (32.9)	985 (32.6)	73 (31.5)	54 (31.6)	858 (32.7)
30–59	164 (4.0)	123 (4.1)	3 (1.3)	4 (2.3)	116 (4.4)
15–29	5 (0.1)	7 (0.2)	0 (0.0)	0 (0.0)	7 (0.3)
< 15	39 (0.9)	53 (1.8)	0 (0.0)	0 (0.0)	53 (2.0)
Missing	376 (9.1)	304 (10.1)	13 (5.6)	5 (2.9)	286 (10.9)
Pre-indexb comorbid diagnoses
CAD/CVD	447 (10.8)	280 (9.3)	28 (12.1)	13 (7.6)	239 (9.1)
Heart failure	132 (3.2)	62 (2.1)	5 (2.2)	2 (1.2)	55 (2.1)
Dyslipidemia	659 (15.9)	412 (13.6)	50 (21.6)	29 (17.0)	333 (12.7)
Hypertension	1521 (36.8)	910 (30.1)	83 (35.8)	64 (37.4)	763 (29.1)
Diabetes mellitus	561 (13.6)	358 (11.8)	30 (12.9)	19 (11.1)	309 (11.8)
Chronic kidney diseasec	281 (6.8)	169 (5.6)	16 (6.9)	9 (5.3)	144 (5.5)
End-stage renal diseased	15 (0.4)	12 (0.4)	1 (0.4)	0 (0.0)	11 (0.4)
Fracture	30 (0.7)	20 (0.7)	4 (1.7)	1 (0.6)	15 (0.6)
Viral hepatitis	1122 (27.1)	902 (29.8)	59 (25.4)	42 (24.6)	801 (30.6)
Tuberculosis	52 (1.3)	46 (1.5)	1 (0.4)	2 (1.2)	43 (1.6)
Psychiatric disorder	1495 (36.1)	1384 (45.8)	119 (51.3)	104 (60.8)	1161 (44.3)
Depression	962 (23.3)	878 (29.0)	83 (35.8)	68 (39.8)	727 (27.7)
Schizophrenia	140 (3.4)	194 (6.4)	12 (5.2)	9 (5.3)	173 (6.6)
Bipolar disorder	635 (15.3)	648 (21.4)	61 (26.3)	63 (36.8)	524 (20.0)
Psychosis	231 (5.6)	308 (10.2)	26 (11.2)	18 (10.5)	264 (10.1)
Posttraumatic stress disorder	366 (8.8)	379 (12.5)	39 (16.8)	52 (30.4)	288 (11.0)
Tobacco usee	1279 (30.9)	918 (30.4)	70 (30.2)	77 (45.0)	771 (29.4)
Alcohol abuse	946 (22.9)	725 (24.0)	59 (25.4)	51 (29.8)	615 (23.5)
Medications
Methadone	49 (1.2)	70 (2.3)	2 (0.9)	2 (1.2)	66 (2.5)
Proton pump inhibitors	1179 (28.5)	713 (23.6)	56 (24.1)	23 (13.5)	634 (24.2)
Bisphosphonates	8 (0.2)	9 (0.3)	0 (0.0)	0 (0.0)	9 (0.3)
Testosterone	66 (1.6)	51 (1.7)	2 (0.9)	3 (1.8)	46 (1.8)

Data are n (%) unless otherwise indicated

aNon-EFV includes the EVG/c, RPV, and RTV-boosted PI groups

bPre-index comorbidities and clinical characteristics were identified in the 6–12-month pre-index period

cDefined as either a chronic kidney disease diagnosis or two consecutive measures of eGFR < 60 ml/min/1.73 m2 occurring at least 30 days apart

dDefined as either a diagnosis of end-stage renal disease, kidney transplant, or dialysis

eIncludes abuse, dependence, rehabilitation, and toxicity related to tobacco

ART antiretroviral therapy, BMI body mass index, CAD coronary artery disease, CVD cerebrovascular disease, EFV efavirenz, eGFR estimated glomerular filtration rate, EVG/c elvitegravir/cobicistat, FTC emtricitabine, PI protease inhibitor (atazanavir, darunavir, or lopinavir), RPV rilpivirine, RTV ritonavir, SD standard deviation, TDF tenofovir disoproxil fumarate

Exposures

Exposures of interest included TDF/FTC (either as a fixed-dose combination or as separate agents) plus one of the following agents: EFV, EVG/c, RPV, or any one of three ritonavir-boosted PIs (i.e., ATV, lopinavir, or darunavir). For regimens with separate dosage forms, the third agent must have overlapped with the backbone within 30 days. For boosted or enhanced regimens (EVG/c and RTV-boosted PIs), the third agent must have also overlapped with the booster/enhancer for the patient to be classified as taking the regimen. Discontinuation of the regimen was defined as having a gap of at least 30 days for either the third agent or the backbone; patients who discontinued their regimen were censored on the first day of the first 30-day gap following the end of the prior days’ supply received by the patient.

Outcomes

The primary composite outcome was a bone adverse event defined as a diagnosis of osteoporosis; a BMD T-score in the osteoporotic or osteopenic ranges for the femoral neck, total spine, distal radius, or total hip; or a diagnosis or procedure code for likely fragility fracture (any hip, wrist/forearm, or spine fracture). BMD T-scores were extracted from patient radiology dual-energy X-ray absorptiometry (DEXA) reports and clinical notes using a previously developed natural language processing tool, with accuracy in the range of 90.4–92.8% [38, 39]. All codes used to identify outcomes of interest are provided in Supplemental digital content, Table 2.

Covariates

To control for confounding and selection bias, we measured baseline covariates that were selected on the basis of potential associations with treatment and/or outcomes, as found in published literature and based on prior clinical knowledge of ART and HIV. These included baseline demographics, baseline HIV laboratory measures, baseline BMD measures and related diagnoses, lifestyle exposures, other comorbidities, medication exposures, and calendar year of the index regimen. All covariates were identified over a 12-month look-back period preceding the index date. Specific definitions for all covariates are provided in the Supplemental digital content, Table 3.

Statistical Analysis

We calculated baseline characteristics overall and by treatment group and used standardized mean differences (SMDs) to compare differences between groups, with SMDs outside the bounds of ± 0.1 indicating meaningful differences [40]. We calculated crude IRs of bone adverse outcomes per 1000 patient-years of exposure and associated exact 95% confidence intervals (CIs) in the unweighted cohort assuming a Poisson distribution. To control for confounding by indication and selection bias, we used inverse probability of treatment weighting (IPTW) for each patient [41]. Weighted Cox proportional hazards regression models were used to estimate adjusted hazard ratios (HRs) for bone outcomes associated with EFV + TDF/FTC compared with EVG/c + TDF/FTC, RPV + TDF/FTC, ritonavir-boosted PIs + TDF/FTC, and all non-EFV regimens combined. To reduce risk of bias, model variability, and unreliable CIs, we excluded any group with < 5 events from analysis comparisons. We used IPTW in the primary analysis; however, if covariate balance was not achieved with IPTW, we conducted a sensitivity analysis using matching weights, which are less sensitive to residual and unmeasured confounding. All analyses were done in SAS version 9.2 (SAS Institute, Cary, NC).

Results

Patient Characteristics

A total of 7161 patients met all eligibility criteria (Fig. 1). Of these, 4137 initiated EFV and 3024 initiated non-EFV regimens, including 232 with EVG/c, 171 with RPV, and 2621 with ritonavir-boosted PIs. At baseline, 26 (0.4%) had received a bisphosphonate (Table 1). For all unadjusted comparisons, psychiatric disorders were overrepresented in those receiving non-EFV-containing regimens versus EFV-containing regimens (Fig. 2 and Supplemental digital content Table 4). After IPTW, covariate balance was achieved among all variables for EFV versus all non-EFV combined (Fig. 2), versus EVG/c (Supplemental digital content, Fig. 1A), and versus ritonavir-boosted PIs (Supplemental digital content, Fig. 1B). For EFV versus RPV, the covariate balance was achieved for 83.0% of variables (Supplemental digital content, Fig. 1C). Mean (standard deviation) follow-up times were 13.0 (19.2) months overall and 15.0 (21.2), 7.6 (7.6), 10.0 (10.5), and 10.5 (16.7) months for EFV, EVG/c, RPV, and ritonavir-boosted PIs, respectively. Corresponding median (interquartile range; min–max) values were 5.0 (2.0–14.2; 0.1–127.7), 5.7 (2.1–17.9; 0.1–110.3), 4.6 (2.1–10.4; 0.2–34.6), 5.9 (2.8–13.4; 0.5–48.3), and 4.2 (1.8–10.7; 0.1–127.7) months.

Fig. 1

Patient selection according to eligibility criteria. ART antiretroviral therapy, ATV atazanavir, BMI body mass index, DRV darunavir, EFV efavirenz, EVG/c elvitegravir/cobicistat, FTC emtricitabine, LPV lopinavir, PI protease inhibitor (atazanavir, darunavir, or lopinavir), RPV rilpivirine, RTV ritonavir, TDF tenofovir disoproxil fumarate, VHA Veterans Health Administration. *Bone disease defined as diagnosis of osteoporosis/osteopenia by international classification of disease diagnosis codes, current procedural terminology codes, or classification by bone mineral density test result

Fig. 2

Baseline demographics: verification that IPT weighting achieves baseline covariate balance between TDF/FTC-containing regimens with EFV versus all non-EFV. CAD/CVD coronary artery disease/cerebrovascular disease, CKD chronic kidney disease, EFV efavirenz, eGFR estimated glomerular filtration rate, FTC emtricitabine, IPT inverse probability of treatment, TDF tenofovir disoproxil fumarate

Incidence and Risk of Any Bone Outcome

The unadjusted (crude) IRs and adjusted HRs for all comparisons are summarized in Fig. 3. The crude IR of the composite bone outcome was lower in the EFV group than in the EVG/c, RPV, and ritonavir-boosted PI groups. In adjusted analyses, EFV was associated with a statistically significant 31% lower risk of the composite bone outcome than non-EFV groups combined, 25% lower than RPV, and 30% lower than the ritonavir-boosted PI group.

Fig. 3

Unadjusted incidence rates and IPTW adjusted hazard ratios of bone adverse outcomes among patients treated with different TDF/FTC-containing regimens with at least five events per group. CI confidence interval, CVD cerebrovascular disease, diag. diagnosis, EFV efavirenz, EVG/c elvitegravir/cobicistat, FTC emtricitabine, IPTW inverse probability of treatment weighting, PI protease inhibitor (atazanavir, darunavir, or lopinavir), PY patient-years, RPV rilpivirine, RTV ritonavir, TDF tenofovir disoproxil fumarate, – insufficient events to calculate hazard ratio. aMeasured at the femoral neck, total spine, distal radius, or total hip

Incidence and Risk of Osteoporosis or Osteopenia

The crude IR for osteoporosis was lower for EFV compared with all other regimens. In adjusted analyses, EFV was associated with a statistically significant 36% lower risk of osteoporosis than non-EFV groups combined, 53% lower than EVG/c, and 35% lower than the ritonavir-boosted PI group. Too few events were observed with RPV to make comparisons. The crude IR of osteopenia was lower for EFV compared with EVG/c, but for all other comparisons IRs were similar. In adjusted analyses, the risk of osteopenia was similar between EFV and the ritonavir-boosted PI group (too few events were observed to report comparisons for EVG/c or RPV). When combining osteoporosis and osteopenia events, EFV was associated with a 21% lower risk compared with the ritonavir-boosted PI group. Frequencies of available pre- and post-index DEXA scans as well as baseline characteristics by DEXA scan status are presented in the Supplemental digital content, Tables 6 and 7, respectively.

Incidence and Risk of Fracture

The crude IR of fragility fractures was lower for EFV compared with EVG/c, RPV, and ritonavir-boosted PIs. In adjusted analyses, risk differences were statistically significant for EFV versus all non-EFV combined (41% lower), RPV (57% lower), and ritonavir-boosted PIs (40% lower).

Of the individual fracture outcomes, EFV had a lower incidence compared with each non-EFV regimen for all fracture sites in all comparisons except for hip fractures with ritonavir-boosted PIs, which had too few events to make comparisons. In adjusted analyses, this corresponded to statistically significant reduced risks for EFV compared with all non-EFV regimens combined for vertebral fracture (51% lower) and wrist/forearm fracture (60% lower). Significantly lower risks were also observed versus ritonavir-boosted PIs for vertebral fracture (48% lower) and wrist/forearm fracture (60% lower). Too few events were observed for individual fracture outcomes to make comparisons with EVG/c or RPV.

For a summary of event numbers, time to events and follow-up times for all bone adverse outcomes, see Supplemental digital content, Table 5.

Sensitivity Analysis

In matching weight analyses, the risks of any bone outcome or fracture for RPV versus EFV showed a similar magnitude, suggesting the IPTW findings were robust regarding residual confounding present for 7/47 (14.9%) of variable comparisons for EFV versus RPV.

Discussion

EFV + TDF/FTC was generally associated with lower risks of bone adverse outcomes compared with TDF-containing regimens with EVG/c, RPV, or ritonavir-boosted PIs in the VHA cohort over an average follow-up time of 13.0 months. These findings suggest that patients with HIV infection receiving EFV + TDF/FTC, which has no known drug-drug interactions with TDF or FTC, or drug-food interactions, may be at lower overall risk of bone adverse outcomes compared with those receiving TDF-containing regimens with EVG/c, RPV, or ritonavir-boosted PIs.

This study is the first to identify an association between EFV + TDF/FTC use among veterans and a reduced risk of bone adverse outcomes versus other TDF-containing regimens. These findings confirm those of Nkhoma et al. [12], who conducted a large claims database analysis of bone fractures associated with TDF/FTC and various third agents (although adjusted analyses were not possible because of insufficient numbers of fracture events). Nkhoma et al. [12] found a lower fracture IR with EFV + TDF/FTC (3.4 per 1000 patient-years) compared with RPV + TDF/FTC (3.6 per 1000 patient-years) or with EVG/c + TDF/FTC (4.4 per 1000 patient-years). The IRs in this study followed the same general trend as in Nkhoma et al. [12], but were higher in magnitude, possibly because the current study was conducted in a higher-risk population in a relatively closed, integrated care network.

Cobicistat, rilpivirine, and ritonavir increase serum creatinine and decrease the estimated glomerular filtration rate (GFR) via inhibition of renal tubular transporters [42, 43], but this does not appear to affect actual GFR as measured by iohexol clearance [44]. Therefore, the bone effects we observed were more likely due to relative differences in tenofovir exposure across the evaluated regimens. Both RPV/FTC/TDF and EVG/COBI/FTC/TDF are recommended to be administered with food, which increases tenofovir exposure compared with fasting administration [33, 34]. In addition, both protease inhibitors (including RTV) [45] and cobicistat [46] increase tenofovir exposure by inhibition of intestinal P-gp-mediated efflux of TDF. Taken together, these effects on tenofovir exposure may have contributed to the increased risk of adverse bone outcomes that we observed with EVG/c/FTC/TDF, RPV/TDF/FTC, and RTV-boosted PIs in this analysis.

For patients with a high risk of bone adverse outcomes, the use of the novel formulation tenofovir alafenamide (TAF) or abacavir has been associated with lower BMD losses at the time of therapy initiation compared with TDF [9]. TAF, a prodrug of tenofovir, is associated with a 91% lower plasma tenofovir concentration than that following TDF administration while maintaining higher intracellular tenofovir concentrations in peripheral blood mononuclear cells for HIV suppression [47]. Biomarkers of bone turnover appear to be less affected by TAF-containing regimens compared with TDF-containing regimens. TAF is now recommended in HIV treatment guidelines [9], and its use is being preferred to TDF by HIV-care providers concerned about bone toxicity of TDF in the aging HIV population.

Given the ongoing widespread use of EFV + TDF/FTC and the absence of an EFV co-formulation with TAF, the results of the current study may reassure physicians and their patients about the bone safety of this combination. Moreover, where alternatives to TDF are limited (such as in resource-limited settings), or where use of generic EFV or TDF as a cost reduction strategy is available, or where EFV + TDF/FTC is available through the Medicines Patent Pool, the choice of a third agent remains critical to long-term safety. The reduced risk of bone adverse outcomes for EFV + TDF/FTC found in this study is of high relevance, especially in resource-limited settings where the cost effectiveness of the fixed-dose combination has achieved widespread use [48]. However, while it appears that EFV + TDF/FTC has a lower risk for bone adverse outcomes than PIs and other boosted regimens in combination with TDF/FTC, our study is not generalizable to the use of unboosted integrase inhibitors, for which additional research is needed.

Other options to counter antiretroviral-associated BMD loss include vitamin D/calcium and/or zoledronic acid supplementation. Vitamin D/calcium supplementation lessened BMD loss among patients receiving EFV + TDF/FTC over 48 weeks [21, 23], and single-dose zoledronic acid administered at the initiation of a PI-containing TDF/FTC regimen can prevent the initial BMD loss [49]. These options, together with careful antiretroviral choice, should be considered between patients and health care providers, taking into account the drug resistance profile, treatment history, other comorbidities, and the risk or tolerability of side effects.

Strengths and limitations of this study are those common to large epidemiologic studies. The main strengths of our study were its large sample size and detailed data from VHA data sets, as well as the use of a natural language-processing tool, which allowed us to extract BMD results from radiology and clinical notes. Our study also has limitations. The VHA population was more than 95% male, which limits generalizability to female populations that could be affected differently by the predictors identified in our study. Follow-up times were short, which may have led to underascertainment of relevant bone adverse outcomes; however, statistically significant between-regimen differences for relevant bone adverse outcomes were noted despite these follow-up times. Moreover, other data have demonstrated that the incidence of fracture is highest during the first and second years after ART initiation, tailing off thereafter [50]; fracture risk among men with HIV infection is higher among older individuals [51] as was the case in this VA population. Thus, we consider that follow-up times in the current study were sufficient to detect differences in BMD and fracture risk. Although IPTW was successfully used to adjust for selection bias and measured confounders for all comparisons with the exception of EFV versus RPV, the potential for unmeasured confounders and incomplete adjustment for measured confounders cannot be ruled out. For example, patients with a range of psychiatric disorders were less likely to be prescribed EFV-containing regimens across all comparisons (Fig. 2 and Supplemental digital content, Fig. 1), likely because of channeling bias consequent to the side effect profile of EFV in patients with severe psychiatric symptoms [52]. Various psychiatric disorders, such as schizophrenia and depression, and psychotropic medications are associated with BMD loss and an increased risk of fracture [53]. However, our method of controlling for confounding effectively balanced the observed differences in these measured characteristics. Specifically, the six psychiatric covariates examined were balanced after IPTW for the EFV versus EVG/c and EFV versus ritonavir-boosted PI comparisons. For the EFV versus RPV comparison, only one remained imbalanced after IPTW, and this was by a small margin (SMD of 0.102). In addition, the more conservative matching weights analyses produced qualitatively similar results to those using IPTW, making it highly unlikely that channeling bias affected our results. Finally, as in any retrospective observational study, causal associations cannot be proven; thus, these findings require confirmation in further prospective studies.

Conclusion

In conclusion, EFV + TDF/FTC was generally associated with a lower incidence of bone adverse outcomes, including osteoporosis, any major fracture, vertebral fracture, and wrist/forearm fracture, compared with other TDF/FTC-containing regimens in the VHA. The third agent in antiretroviral regimens may have a significant effect on the risk of bone adverse events associated with TDF.

Below is the link to the electronic supplementary material.

Supplementary material 1 (PDF 512 kb)