Background. Human immunodeficiency virus (HIV)–infected persons may be at increased risk for developing type 2 diabetes mellitus because of viral coinfection and adverse effects of treatment.
Methods. We studied associations of new-onset diabetes mellitus with hepatitis B virus and hepatitis C virus coinfections and antiretroviral therapy in participants in the Swiss HIV Cohort Study, using Poisson regression.
Results. A total of 123 of 6513 persons experienced diabetes mellitus during 27,798 person-years of follow-up (PYFU), resulting in an incidence of 4.4 cases per 1000 PYFU (95% confidence interval [CI], 3.7–5.3 cases per 1000 PYFU). An increased incidence rate ratio (IRR) was found for male subjects (IRR, 2.5; 95% CI, 1.5–4.2), older age (IRR for subjects >60 years old, 4.3; 95% CI, 2.3–8.2), black (IRR, 2.1; 95% CI, 1.1–4.0) and Asian (IRR, 4.9; 95% CI, 2.2–10.9) ethnicity, Centers for Disease Control and Prevention disease stage C (IRR, 1.6; 95% CI, 1.04–2.4), and obesity (IRR, 4.7; 95% CI, 3.1–7.0), but results for hepatitis C virus infection or active hepatitis B virus infection were inconclusive. Strong associations were found for current treatment with nucleoside reverse-transcriptase inhibitors (IRR, 2.22; 95% CI, 1.11–4.45), nucleoside reverse-transcriptase inhibitors plus protease inhibitors (IRR, 2.48; 95% CI, 1.42–4.31), and nucleoside reverse-transcriptase inhibitors plus protease inhibitors and nonnucleoside reverse-transcriptase inhibitors (IRR, 3.25; 95% CI, 1.59–6.67) but were not found for treatment with nucleoside reverse-transcriptase inhibitors plus nonnucleoside reverse-transcriptase inhibitors (IRR, 1.47; 95% CI, 0.77–2.82).
Conclusions. In addition to traditional risk factors, current treatment with protease inhibitor– and nucleoside reverse-transcriptase inhibitor–containing regimens was associated with the risk of developing type 2 diabetes mellitus. Our study did not find a significant association between viral hepatitis infection and risk of incident diabetes.
In 1997, the US Food and Drug Administration reported the association of hyperglycemia and new-onset type 2 diabetes mellitus (DM) with protease inhibitors (PIs; saquinavir, ritonavir, indinavir, and nelfinavir) [1]. Subsequent studies have confirmed the association of hyperglycemia or DM with PI use [2,3,4,5,6–7]. More recently, nucleoside and nucleotide reverse-transcriptase inhibitors (NRTIs), but not nonnucleoside reverse-transcriptase inhibitors (NNRTIs), were found to contribute to the disturbance of glucose metabolism [8,9,10–11]. Furthermore, associations of hyperglycemia and DM with hepatitis C virus (HCV) infection have been reported both in HIV-negative [12,13–14] and HIV-positive [15,16–17] populations. Potential mechanisms may include HCV-induced insulin resistance mediated by proinflammatory cytokines [18], immune reactions against pancreatic β-cells, or direct infection of β-cells by HCV [19].
Patterns of use of antiretroviral regimens have changed over the years in response to perceived toxicity and to the availability of new drugs. Because of improved life expectancy, increased cumulative exposure to antiretroviral drugs in HIV-infected persons may have resulted in cumulative toxicity in some patients. This study was designed to assess the impact of hepatitis coinfection on the development of new-onset DM in the Swiss HIV Cohort Study (SHCS), taking into account changes in anthropometric risk factors and antiretroviral therapy (ART) during follow-up.
Study population. The SHCS is an ongoing, prospective, clinic-based study that was established in 1988 and that continuously enrolls and observes HIV-1–infected individuals aged ⩾16 years at 5 university outpatient clinics, 2 large district hospitals, affiliated regional hospitals, and private practices [20]. The study was approved by local ethical review boards, and written informed consent was obtained from all participants.
For this analysis, we included SHCS participants if they had at least 2 study visits with at least 1 year of follow-up after 1 March 2000. Patients with prevalent cases of DM were excluded.
Ascertainment of DM. As part of the Data Collection on Adverse Events of Anti-HIV Drugs (DAD) protocol [21], cardiovascular risk factors and events, including plasma glucose levels and diagnosis of DM, have been routinely collected and verified at biannual follow-up visits since 1 March 2000. DM is diagnosed according to the criteria of the Expert Committee on the Diagnosis and Classification of DM [22], with confirmed plasma glucose level cut-off values of >7.0 mmol/L (fasting) and >11.1 mmol/L (nonfasting). Treating physicians were contacted to ascertain the status of patients with glucose levels that exceeded the cut-off value or who were receiving antidiabetic medication without an explicit diagnosis of DM.
Definitions. Body mass index (BMI), calculated as weight in kilograms divided by the square of the height in meters, was stratified into <18.5 (underweight), 18.5–24.9 (normal), 25–29.9 (overweight), and ⩾30 (obese) [23], without modification for Asians [24]. For central obesity, we used the new worldwide definition [25], with sex- and ethnicity-specific waist circumference cut-off values as follows: Europeans and sub-Saharan Africans, ⩾94 cm for males and ⩾80 cm for females; south Asians, Chinese, and south and central Americans, ⩾90 cm for males and ⩾80 cm for females; and Japanese, ⩾85 cm for males and ⩾90 cm for females. Elevated blood pressure was defined as present in patients with diastolic blood pressure ⩾85 mmHg or systolic blood pressure ⩾130 mmHg or who were receiving antihypertensive treatment. HCV infection was defined as present in patients who were seropositive for HCV and who had test results positive for HCV RNA. Active hepatitis B virus (HBV) infection was defined as present in patients who were seropositive for HBV and hepatitis B surface antigen or hepatitis B e antigen or who had test results positive for HBV DNA.
Statistical analysis. We defined the rate of new-onset DM as the number of cases of DM divided by the total number of person-years of follow-up (PYFU). Follow-up was counted from the first visit after 1 March 2000 (baseline) to the date of the first diagnosis of DM, death, or the patient's last cohort visit, whichever occurred first.
Associations between incident DM were analyzed in univariable and multivariable Poisson regression models. Fixed covariables were sex, ethnicity, injection drug use as the presumed mode of HIV acquisition, and HBV and/or HCV infection status. The following characteristics were analyzed in 2 separate models as fixed (baseline) and time-updated covariables, respectively: Age (prespecified strata of 16–39 years old, 40–49 years old, 50–59 years old, and ⩾60 years old), smoking status (never smoked, former smoker, or current smoker), central obesity, elevated blood pressure, and CD4+ cell count, stratified according to Centers for Disease Control and Prevention categories (<200 cells/µL, 200–499 cells/µL, and ⩾500 cells/µL) [26].
Dates of starting and stopping treatment with antiretroviral drugs were documented prospectively by the treating physician and checked for plausibility and completeness by trained data managers prior to entry into the SHCS database. We defined combination ART (cART), as distinct from any ART, as an ART regiment in which at least 3 antiretroviral drugs were given simultaneously. Time-updated antiretroviral treatment information included (1) years of exposure to ART, cART, the different drug classes (NRTI, PI, and NNRTI), or individual drugs and (2) current receipt of a specific drug class or individual drugs. The time-updated cumulative exposure variables are increasing with continuous intake and are kept constant after discontinuation of the specific treatment (drug class or drug), whereas current treatment variables can have alternating values of 0 and 1 over time, depending on whether the specific treatment is currently received (1) or not (0). Because of obvious collinearities, we present separate multivariable models for cumulative exposure and current exposure variables.
Characteristics of different patient groups were compared using Wilcoxon rank-sum tests for continuous variables and χ2 tests or Fisher's exact test for categorical variables. We used Stata software, version 9.2 (StataCorp), for analysis.
Sensitivity analysis. We performed different sensitivity analyses by adjusting for calendar period (2000–2001, 2002–2003, and 2004 and after), looking for trends over time, and using separate models for male and female subjects, individuals aged <50 years and those aged ⩾50 years, participants with African ethnicity, patients receiving ritonavir-boosted PI regimens, and patients who were not receiving ART versus patients receiving their first regimen versus patients receiving later regimens. We considered continuous values for years of cumulative exposure, but we also tested the linear relationship by stratifying exposure into 0 years, >0 to x years, and >x years, with x taking the values 1 year, 2 years, and 3 years.
Because current antiretroviral drug use could be subject to reversed causality problems (e.g., patients who switched to another drug because of increasing glucose levels and whose diagnosis of DM was then attributed to the newly initiated drug regimen), we performed several sensitivity analyses in which we lagged the starting and stopping dates of drugs by 1 month, 2 months, and 3 months, yielding estimates for regimens that the patients had received 1 month, 2 months, or 3 months prior to the diagnosis of DM. Finally, we checked whether findings were influenced by the inclusion of other time-updated anthropometric measures, such as BMI, waist-to-hip ratio, and perceived lipohypertrophy and lipoatrophy, or by the exclusion of hypertension, because this may be on the causal pathway to DM.
Of 8253 SHCS participants seen after 1 March 2000, 6681 (81%) had at least 2 follow-up visits over a period of at least 1 year. of these, we excluded 130 patients (1.9%) with a preexisting diagnosis of DM and 38 individuals for whom waist circumference or BMI were unknown. The present study is thus based on 6513 individuals with follow-up visits between 1 March 2000 and 17 July 2006. Baseline characteristics are listed in table 1. Participants (31% of whom were female) had a median age of 38 years and were well distributed between the major HIV transmission categories. Almost one-quarter of the patients had illness that was classified as Centers for Disease Control and Prevention disease stage C, and 73% were receiving or had received ART for a median duration of 1.9 years. Baseline and nadir CD4+ cell counts were 403 cells/µL and 230 cells/µL, respectively. HCV infection was diagnosed in 27.4% of subjects, and HBV infection was active in 5.3%. The mean duration of follow-up was 4.3 years, and the cumulative follow-up for all subjects was 27,798 PYFU. During the follow-up period, 123 patients developed new-onset DM, resulting in an incidence of 4.42 cases per 1000 PYFU (95% CI, 3.71–5.28 cases per 1000 PYFU). In the course of the study, 341 patients (5.2%) died, and 656 patients (10.1%) were lost to follow-up. Injection drug users and persons with nonwhite ethnicity were more likely to become lost to follow-up (injection drug users vs. non–injection drug users, 15.0% vs. 8.6%; nonwhite ethnicity vs. white ethnicity, 13.5% vs. 7.3%; P < .001, by the χ2 test). In addition, persons lost to follow-up had less advanced HIV infection than did persons who were not lost to follow-up, with a higher median nadir CD4+ cell count (270 cells/µL vs. 226 cells/µL; P < .001) and fewer individuals with Centers for Disease Control and Prevention disease stage C at baseline (18% vs. 23%; P = .004).
Incidence rate ratios (IRRs) for the development of new-onset type 2 diabetes mellitus (DM) based on 123 events among 6513 participants with 27,798 person-years of follow-up. Shown are associations with current receipt of specific drug classes and individual protease inhibitor (PI) and nucleoside or nucleotide reverse-transcriptase inhibitor (NRTI) combinations. Multivariable Poisson models were adjusted for all variables listed in table 2. 3TC, lamivudine; ABC, abacavir; ATV, atazanavir; AZT, zidovudine; D4T, stavudine; DDI, didanosine; IDV, indinavir; LPV, lopinavir; NFV, nelfinavir; NNRTI, nonnucleoside reverse-transcriptase inhibitor; SQV, saquinavir; TDF, tenofovir.
Comparison of sex- and age-specific incidence rates of type 2 diabetes mellitus (DM) between the present study of HIV-infected persons, the MONICA Augsburg Cohort (persons not infected with HIV) [27], and the US population in 2004 [30]. Whiskers indicate standard error. PYFU, person-years of follow-up.
Baseline characteristics of 6513 subjects with and without new-onset type 2 diabetes mellitus (DM) during follow-up.
Multivariable models. Baseline demographic, clinical, and anthropometric covariables known to be potentially associated with an increased risk of developing DM were used to build a baseline multivariable Poisson model (table 2). Male sex, older age, African or Asian ethnicity, clinical AIDS, and central obesity were strong predictors of DM. Simultaneous estimates for BMI and central obesity are potentially inaccurate because of collinearity. Therefore, we removed BMI from the multivariable model.
Univariable and multivariable Poisson regression of baseline demographic, clinical, and anthropometric covariables potentially affecting the risk of developing type 2 diabetes mellitus, based on 6513 individuals with 123 events.
We then built a time-updated model by replacing the baseline variables for age, CD4+ cell count, nadir CD4+ cell count, Centers for Disease Control and Prevention disease stage, smoking status, hypertension, and central obesity with their time-updated counterparts (i.e., the latest data available at each moment in time). Results of this time-updated model were virtually identical (data not shown), except that time-update hypertension also became statistically significant, with an incidence rate ratio (IRR) of 1.65 (95% CI, 1.10–2.48). Both the baseline model and the time-updated model will be used for adjustment in the following analyses.
Association of hepatitis and DM. Next, we analyzed the contribution of HBV and HCV infection to the development of DM in univariable models and together with the baseline and time-updated models (table 3). None of the HBV and HCV infection categories were associated with the incidence of DM (all P values were >.2). These findings remained unchanged in sensitivity analyses involving persons not receiving cART, persons receiving their first cART regimen, and persons receiving subsequent cART regimens (all interaction terms had P values >.6). To unmask potential confounding by ART-induced DM, we then included treatment-associated cofactors.
Univariable and multivariable Poisson regression of hepatitis-associated risk of developing type 2 diabetes mellitus, based on 6513 individuals with 123 events.
Association of antiretrovial treatment and DM. There was no clear effect of cumulative exposure to the different drug classes on the incidence of DM; univariable IRR per year of exposure to NRTI, PI, and NNRTI therapy were 1.04 (95% CI, 0.99–1.10), 1.05 (95% CI, 0.98–1.12) and 1.02 (95% CI, 0.90–1.14), respectively. Also, in the baseline and the time-updated multivariable models, as well as in 3 models with cumulative exposure (stratified into none, 0–1 year, and >1 year; none, 0–2 years, and >2 years; and none, 0–3 years, and >3 years), there was no statistically significant association (all P values were >.5).
However, current exposure to NRTI therapy, NRTI and PI combination therapy, or NRTI, PI, and NNRTI combination therapy increased the risk of developing DM in the univariable model, with IRRs of 2.22 (1.11–4.45), 2.48 (1.42–4.31), and 3.25 (1.59–6.67), respectively; there was no such association with current exposure to NRTI and NNRTI combination therapy, which was associated with an IRR of 1.47 (0.77–2.82). Adjusting with the variables from the baseline model (table 2) and the time-updated model confirmed the findings from the univariable models (figure 1). It should be noted that receiving a PI and receiving an NRTI are highly collinear, because, until recently, PIs were usually combined with pairs of NRTIs. Therefore, we analyzed individual drugs from the PI class and the most frequent combinations of NRTIs in 2 separate models and combined in a single model. Again, these 3 models were repeated multivariably with adjustment for the baseline and time-updated variables. Results were similar (figure 1). In this multivariable analysis, DM was independently associated with current exposure to indinavir (IRR, 2.03 [95% CI, 1.05–3.93]), lamivudine-stavudine (IRR, 2.62 [95% CI, 1.22–5.61]), didanosine-stavudine (IRR, 3.09 [95% CI, 1.29–7.39]), and didanosine-tenofovir (IRR, 3.94 [95% CI, 1.57–9.87]), but other PI and NRTI combinations also showed trends.
There were no apparent interactions between PI and NRTI combinations, and results from sensitivity analyses with ART variables lagged by 1 month, 2 months, and 3 months yielded consistent results (data not shown). Findings also did not change when the model was adjusted for calendar period (2000–2001, 2002–2003, or 2004 and after), and we did not find a statistically significant trend over time (P = .28). Further sensitivity analyses included separate models for male and female subjects, for individuals aged <50 years and those aged ⩾50 years, for participants with African ethnicity, for patients with ritonavir-boosted PI regimens, and for patients receiving no ART versus patients receiving their first ART regimen versus patients receiving later regimens. None of these analyses revealed appreciable alterations of the main findings. Similarly, findings were unchanged in models incorporating BMI or waist-to-hip ratio instead of central obesity, including lipohypertrophy and lipoatrophy, or excluding hypertension. Finally, estimates for the association of HCV and HBV coinfection with DM remained virtually unchanged when adjusting for antiretroviral treatment: HCV infection was associated with an IRR of 1.10 (95% CI, 0.57–2.12), and active HBV infection was associated with an IRR of 1.42 (95% CI, 0.65–3.13).
In this analysis of 6513 participants who were followed-up for 27,798 person-years in the SHCS, we found an incidence of DM of 4.42 cases per 1000 PYFU. Similar to the HIV-seronegative population, factors such as male sex, older age, African or Asian ethnicity, and obesity were strong predictors of DM in this population [27,28–29].
Because there are no data available for the incidence of DM in Switzerland, we compared our results with published results of a large cohort of HIV-negative persons from Germany for the period 1984–1998 [27]. As shown in figure 2, the sex- and age-specific incidence rates are very similar, except for the highest age groups. There are, however, substantial differences in prevalence and incidence of DM between Europe and the United States. The incidence in the United States [30] is approximately twice as high as the incidence in Germany, and it is also higher than the findings from our study in most strata (figure 2). This higher background incidence may, therefore, partly explain the observed higher incidences in US cohorts of HIV-infected persons, such as the Multicenter AIDS Cohort Study [7], which reported incidences of 47 cases per 1000 PYFU (95% CI, 32–71 cases per 1000 PYFU) among individuals receiving cART, 17 cases per 1000 PYFU (95% CI, 6–45 cases per 1000 PYFU) among individuals not receiving cART, and 14 cases per 1000 PYFU (95% CI, 8–26 cases per 1000 PYFU) among HIV-seronegative individuals. Other factors contributing to the high incidences in the Multicenter AIDS Cohort Study [7] may include that the definition of DM was based on a single fasting plasma glucose determination without confirmation, that the subjects had a median BMI that was 3 greater than that in our study, and that their subjects had a median age that was ∼10 years older than that of our subjects. From 1994 through 1998, the Women's Interagency HIV Study found incidences of 28 cases per 1000 PYFU (95% CI, 16–41 cases per 1000 PYFU) among individuals receiving PI-containing regimens, 12 cases per 1000 PYFU (95% CI, 7–18 cases per 1000 PYFU) among individuals receiving NRTIs or NNRTIs only or no ART, and 13 cases per 1000 PYFU (95% CI, 7–22 cases per 1000 PYFU) among HIV-seronegative participants [6]. Although age and BMI data appear to be similar to data in our study, 55% and 26% of HIV-infected Women's Interagency HIV Study participants were black and Hispanic, respectively, compared with only 22.7% and 2.6% of women with black and Hispanic ethnicity in our study. Analysis of the Veterans Affairs administrative database, which includes almost 27,000 male individuals free of DM at baseline, identified substantial increases in the risk of developing DM in the cART period, compared with the pre-cART period, among nonwhite individuals and older individuals [15]. In addition, HCV seropositivity in the cART era was associated with a rather modest, but statistically significant, hazard ratio of 1.39 (95% CI, 1.27–1.53) in the multivariable analyses. Our study may have failed to find a significant association with HCV infection because of our limited sample size.
Several antiretroviral drugs and drug combinations were related to the development of DM; in particular, these include indinavir, lamivudine-stavudine, didanosine-stavudine, and didanosine-tenofovir. Although several studies have reported associations between DM and PI use [1,2,3,4,5,6–7], only limited data are available on the association of DM or hyperglycemia with exposure to NRTIs. Regimens including stavudine plus indinavir were found to increase the risk for DM in an Italian study [8], and hyperglycemia was associated with regimens that included didanosine plus tenofovir in Spain [10]. In another report from the Multicenter AIDS Cohort Study, increased insulin resistance was associated with cumulative exposure to regimens containing indinavir, stavudine, and lamivudine [9].
Several mechanisms have been postulated to explain how antiretroviral drug toxicity may lead to insulin resistance and DM. Insulin sensitivity is reduced by a single dose of indinavir [31], and there are specific effects of PIs on the glucose transporter GLUT4 [32] and of NRTIs on mitochondria. Among the currently used NRTIs, the strongest association with mitochondrial toxicity, measured as inhibition of the mitochondrial DNA polymerase-γ, is found for didanosine and stavudine (known as the d-drugs) [33]; notably, these 2 drugs were strongly associated with DM in our study. Lipodystrophy is a crucial aspect of the association of cART with insulin resistance, leading to relative preponderance of visceral fat, hepatic steatosis, and fat deposition at other "ectopic" sites. In the fatless mouse animal model, diabetes and insulin resistance develop but respond to transplantation of subcutaneous fat [34]. HIV-infected persons with lipodystrophy, compared with those without lipodystrophy, have a reduction in plasma adiponectin and adipose tissue adiponectin mRNA levels of ∼50%, correlating with insulin resistance and with increased cytokine levels [35]. Finally, there may also be a direct effect of HIV on the pancreas, because there have been reports of the spontaneous resolution of DM in patients starting cART [36, 37].
The strengths of this clinic-based study are the large number of patient-years with prospectively collected anthropometric data and plasma glucose and treatment information, as well as the verification of hyperglycemia values in relation to a DM diagnosis by the treating physicians. Unfortunately, we do not have information on family history of DM and adherence to antiretroviral drug regimens.
As the life expectancy of HIV-infected persons with access to cART has dramatically increased over the last decade, the risks for metabolic problems have also increased. Strong predictors for DM are nonmodifiable characteristics, such as age and ethnicity, but, importantly, the strong predictors also include obesity, which should become a major target for prevention. Furthermore, there are risks for DM associated with ART, especially therapy with PIs and some NRTI combinations. Because of their association with other metabolic disorders, regimens containing stavudine and didanosine are avoided as long as possible in developed countries, but they belong to first-line regimens in resource-limited areas. Together with the probably genetically determined elevated DM risk associated with Asian and African ethnicity, this may have an important impact on the long-term tolerability of anti-HIV treatment in the regions that are most affected.
In conclusion, type 2 DM is a serious but partially preventable complication in HIV-infected persons. In this study, we have identified several drugs and drug combinations that are associated with an increased risk of DM. As the number of antiretroviral drugs increases, continued monitoring in large cohort studies and collaborations remains an important prerequisite for the early identification of new toxicities and long-term side effects.
M. Battegay, E. Bernasconi, J. Böni, H. Bucher, Ph. Bürgisser, S. Cattacin, M. Cavassini, R. Dubs, M. Egger, L. Elzi, P. Erb, M. Fischer, M. Flepp, A. Fontana, P. Francioli (President of the SHCS, Centre Hospitalier Universitaire Vaudois, Lausanne), H. Furrer (Chairman of the Clinical and Laboratory Committee), M. Gorgievski, H. Günthard, H. Hirsch, B. Hirschel, I. Hösli, Ch. Kahlert, L. Kaiser, U. Karrer, C. Kind, Th. Klimkait, B. Ledergerber, G. Martinetti, B. Martinez, N. Müller, D. Nadal, M. Opravil, F. Paccaud, G. Pantaleo, M. Rickenbach (Head of Data Center), C. Rudin (Chairman of the Mother and Child Substudy), P. Schmid, D. Schultze, J. Schüpbach, R. Speck, P. Taffé, P. Tarr, A. Telenti, A. Trkola, P. Vernazza (Chairman of the Scientific Board), R. Weber, and S. Yerly
Financial support. The Swiss National Science Foundation.
Potential conflicts of interest. B.L. has received travel grants and honoraria from Abbott, Aventis, Bristol-Myers Squibb, Gilead Sciences, GlaxoSmithKline, Merck Sharp & Dohme, Roche, and Tibotec. H.F. has participated in advisory boards of Abbott, GlaxoSmithKline, Bristol-Myers Squibb, Roche, Gilead, Merck Sharp & Dohme, and Boehringer-Ingelheim. The institution with which H.F. is affiliated has received unrestricted educational grants from Abbott, GlaxoSmithKline, Bristol-Myers Squibb, Roche, Gilead, Merck Sharp & Dohme, and Boehringer-Ingelheim. M.R. has received a travel grant from GlaxoSmithKline. E.B. has received travel grants and honoraria from GlaxoSmithKline, Roche, Abbott, Merck Sharp & Dohme, and Tibotec. R.W. has received travel grants and honoraria for teaching from Abbott, Boehringer Ingelheim, Bristol-Myers Squibb, Gilead, GlaxoSmithKline, Merck, Pfizer, Roche, and TRB Chemedica. All other authors: no conflicts.
↵Members of the Swiss HIV Cohort Study are listed at the end of the text.
IDSA Members: For your free access to this journal, log in via the IDSA members area.
Open access options for authors visit Oxford Open
This journal enables compliance with the NIH Public Access Policy