Relationship between combination antiretroviral therapy regimens and diabetes mellitus-related comorbidities among HIV patients in Gaborone Botswana

Background

Combination antiretroviral therapy (cARTs) regiments are known to prolong the recipients’ life even though they are risk factors for diabetes mellitus-related comorbidities (DRCs). We sought to: (i) examine cART relationship with DRCs among patients attending HIV clinics in Gaborone, Botswana (which cART regimens are associated with shorter/longer time to the event), (ii) characterize patients’ underlying biomedical and demographic risk factors of DRC and identify the most important, (iii) investigate survival of patients on different cART regimens in the presence of these risk factors.

Methods

Data from two major HIV clinics in Botswana were reviewed. Relationships between different cART regimens and DRCs were investigated among 531 recipients. Recipients’ DRC risk factors were identified. Cox regression model was run. Unadjusted and adjusted hazard ratios were computed, and hazard and survival functions for different cART regimens were plotted.

Results

Major findings were: patients on second- and third-line cART were less likely to develop DRCs earlier than those on first-line cART. Patients with CD4 count ≤ 200 cells/mm³ at cART initiation were more likely to develop DRCs earlier than those who had CD4 count > 200 cells/mm³. Overweight patients at cART initiation had a higher risk of developing DRCs earlier than those who had normal body mass index. Males had a lower risk of developing DRCs earlier than females.

Conclusion

The risk of new onset of DRC among cART recipients is a function of the type of cART regimen, duration of exposure and patients’ underlying biomedical and demographic DRC risk factors. The study has provided a survival model highlighting DRCs’ significant prognostic factors to guide clinical care, policy and management of recipients of cARTs. Further studies in the same direction will likely improve the survival to the development of DRC of every cART recipient in this community.

Human Immunodeficiency Virus (HIV) patients are now living longer than before [1, 2]. This is mainly attributed to access to free combination antiretroviral therapy (cART) [2] or combinations consisting of a minimum of two active drugs from two classes of antiretroviral agents; usually they contain more different active drugs and are referred to as HAART. The improvement in life expectancy and quality of life observed among cART recipients could have been even better had these drugs not been associated with illnesses such as diabetes and diabetes-related comorbidities (DRC) [3, 4] or any comorbidity associated with Type 2 diabetes [5].

Association between cART and DRCs has been long suspected [4, 6, 7]. Studies have provided evidence showing that this association may depend not only on the duration of exposure to cART [8] but also the type of cART regimen taken [4]. In other words, while cART is considered a providential life-saving treatment for persons living with HIV (PLWH), it also has the potential of taking away lives by exposing recipients to an array of unwanted illnesses like diabetes mellitus and related complications [3]. This is a disturbing situation that makes the task of controlling and preventing morbidity and mortality among HIV patients more difficult, especially in settings with a high prevalence of HIV like in Botswana [9].

The initiative to provide free cART in Botswana started in 2002 [9]. Recipients of cART have access to standard treatment as defined by the Botswana national HIV & AIDS treatment guidelines [10] and the handbook of the Botswana integrated HIV clinical care guidelines [11] in use between 2002 and 2015. Components of cART regimens in Botswana were as presented in Appendix and described by Rankgoane-Pono and colleagues’ [12]. All these regimens have been to some extent, in different settings, associated with the development of DRCs after some time of exposure [6, 7, 13, 14]. Some experts have reported these associations only with regimens containing protease inhibitors [13], while others have described DRCs association with all types cART regimens [6, 7, 14]. The average time duration from cART initiation to the development of DRCs also significantly differed from one study to another [8, 14]. However, these results are perhaps due to differences in patient characteristics. For example, the patients’ underlying biomedical and demographic risk factors of DRC, such as being overweight or having CD4 cell count below 200 cells/mm³ have the potential to alter outcomes among patients if not well controlled [14]. Attributes associated with the risk of developing DRCs need to be identified earlier before patients’ enrolment in cART programs. This will allow for suitable care which minimizes cART adverse effects and delays or completely prevents DRC development. The availability of models highlighting DRCs’ significant prognostic factors is important to guide clinical care, policy and management of recipients of cARTs in a setting.

The study aimed to: (i) examine cART relationship with DRCs among patients attending HIV clinics in Gaborone, Botswana (which cART regimens are associated with shorter/longer time to the event), (ii) characterize patients’ underlying biomedical and demographic risk factors of DRC and identify the most important (those that significantly affect time to the event), (iii) investigate survival of patients on different cART regimens in the presence of these risk factors.

Study design

The study was a 12-year open retrospective cohort of cART recipients. The outcome event was DRC. This was recorded as diagnosed by the patient’s physician and listed in the 10th revision of the International Statistical Classification of Diseases and Related Health Problems (ICD-10) [5]. Exposures of interest were cART regimens and patients’ underlying biomedical and demographic characteristics. Survival was calculated as the time elapsed from the date of cART initiation until the date of DRC development. The study was conducted in Gaborone, Botswana’s capital city [15]. HIV period prevalence (2008–2013) was estimated at 19% among persons aged 18 months and above [16].Two clinics, Princess Marina Hospital (PMH) HIV clinic and Bontleng HIV clinic, were chosen as study sites out of 14 sites. The two clinics were selected because of their high capacity to cater for HIV patients and were the best at keeping patient records.

Data collection

Data were collected from patient medical records. The records included patient files, referral and discharge booklets. Only data from patients who met the inclusion criteria were considered. A cohort that included: (i) patients who were receiving the first-line cART, (ii) patients who were receiving the second-line cART and (iii) those on the third-line cART was identified and followed between 2002 and 2015. The follow up endpoint was when DRC occurred or the end of the study. We excluded patients who had DRCs upon entry into the treatment program, pregnant women, recipients who were initiated on cART after 2012 (allowing for at least 3 years of follow-up), recipients aged less than 18 years and recipients with discrepant data from different records within the same clinic to assure data accuracies.

From each patient, we collected data on the date of enrolment into the program, the date of cART initiation, weight in kilograms (weight-1) and height in centimetres when entering the program, weight at the time of data collection (weight-2), CD4 cell count at cART initiation (CD4–1), CD4 cell count at the time of data collection (CD4–2), whether or not the treatment regimen was the same for a patient since the enrolment in the program, cART treatment regimen received before if there was a switch and reason for the switch. Only switches between different cART regimens were considered. Data on whether or not adherence to the treatment had been continued or discontinued and information on DRCs as well as the date of the diagnosis were also collected.

The sample size was estimated as by Hennekens and Buring [17], using a sampling error of 0.05 and a beta level of 0.20. The proportion of baseline DRC among HIV patients recipients of cART was 17.6% [18] and the expected magnitude of association between cART and DRC was set at 1.9 odds ratio. This led to an estimated sample size of 483. However, we included 540 participants in the study.

Data analysis

Data were entered into a computer using Microsoft Excel (Redmond, WA) before being imported into IBM SPSS version 21 (Chicago, IL). Proportions of participants on the first-line cART, the second-line, the third-line and on second-, third-line cART regimens were estimated. Potential associations between the different cART regimens and DRCs were investigated by computing unadjusted hazard ratios (UHR) of each regimen and their 95% confidence intervals (CIs) using univariate Cox regression analysis.

Similarly, proportions of participants with underlying biomedical and demographic characteristics known as potential DRC risk factors were assessed, namely CD4 cell count as a continuous variable, CD4 cell count categories, gender, age as a continuous time dependent variable, age groups (< 35 years old/≥35 years old), body mass index (BMI) categories, whether patients adhered to treatment or not and attended PMH clinic or Bontleng clinic, and whether they switched cART regimens or not. Using DRCs as status /event/dependent variable, time in years, potential associations with these factors were assessed by computing UHR and their 95% CIs in univariate Cox regression analysis. Continuous variables namely CD4 cell count at cART initiation, CD4 cell count after cART initiation and age were not included in this analysis unless they were statistically significant in the two independent samples t-test analysis. Covariates that achieved p ≤ 0.9 in the univariate analysis were re-examined in a multivariate Cox regression model to identify those that significantly (p < 0.05) affected the time to event. Time spent in the study was expressed in years. Survival was estimated as the time elapsed from the date of initiation to cART until the date of the first DRC development or the last attendance. Recipients lost to follow-up, those referred to other facilities were left censored at the last date they were seen. Those who got to the end of the study without developing any type of DRC were right censored at the date of data collection. Adjusted hazard ratios (AHRs) and their 95% CIs were computed. The survival functions were plotted and the probability of survival in years from the date of initiation on: first-line cART (cART line1), second-line cART (cART line2), third-line cART (cART line3) or second-, third-line cART (cART line2/3) to the development of DRCs were compared using Log rank X². The Cox regression presented herein is the best model where only covariates with P < 0.05 are kept in. The overall statistical significance of the model or how well the model fits the data was tested by computing the likelihood chi-square statistic.

Ethical approval and consent to participate

Ethical approval to collect data from HIV clinics was sought and obtained from the University of Botswana Review Board and the ethics committee of the Ministry of Health and Wellness, Botswana. Permission to consult clinic record books and systems was also sought and obtained from the clinic management. As this was a record based study, no consent to participate was required.

Nine of the 540 records that were reviewed did not meet the inclusion criteria and were excluded. Thus, 531 patients’ records were included in the analysis. Of these, 368 (69.3%) were females. The mean age [± standard error of the mean (SEM)] of participants was 41.4 ± 8.9 years. The youngest was 19 years old while the oldest aged 82.0 years. The mean (± SEM) weight of patients at cART initiation was 60.6 ± 11.8 kg with the lowest record of 17.5 kg and 101.0 kg as the highest. After cART initiation the mean (± SEM) weight documented was 67.9 ± 14.50 kg, with the minimum of 32.5 kg and 117.7 kg as the maximum. The mean CD4 count (± SEM) at cART initiation was 139.5 ± 5.11 cells/mm³, the minimum count was 0.00 cells/mm³ and the maximum was 889.0 cells/mm³. After cART initiation, the mean CD4 count was 536.0 ± 10.16 cells/mm³ with the minimum count of 25.0 cells/mm³ and the maximum of 1441.0 cells/mm³. At cART initiation, 408 (76.6%) participants had a CD4 count of ≤200 cells /mm³compared to 34 (6.4%) after the initiation of cART. Four hundred and forty two patients or 83.2% did not develop any DRC, whereas 89 (16.8%) of them developed DRCs.

Results presented in Table 1 show proportions of participants by cART regimen and UHR as measures of crude association of each regimen with DRCs.

Patients on first-line cART were 2.03 times more at risk of developing DRCs earlier than those who were not on the first-line cART [UHR = 2.03; (95% CI: 1.26–3.27)]. Patients on second-line cART had a 52% less risk of developing DRCs earlier than patients who were not on the second-line cART [UHR = 0.48; (95% CI: 0.30–0.79)]. Those on third-line cART had 5% less risk of developing a DRC earlier than patients not on the third-line cART but the relationship was not statistically significant [UHR = 0.95; (95% CI: 0.13–6.83) and p = 0.9]. When patients on second-line cART were aggregated with those on third-line (cART line2/3), recipients showed 51% less risk of developing DRCs earlier than patients who were not on cART line2/3 [UHR 0.49, (95% CI 0.31–0.79).

Proportions of patients by biomedical and demographic characteristics and UHR as measures by their crude association with DRCs are shown in Table 2.

UHR for CD4–1 and CD4–2 count as continuous variables are not shown here as both variables did not achieve significance in the exploratory two independent samples t-test analysis. Age as a continuous variable indicated that the hazard of developing DRCs increased by 3% [UHR = 1.03; (95% CI: 1.01–1.05)] for every one year increase in recipients age. Of all the characteristics investigated, seven showed significant crude associations with DRCs (p < 0.05) while others achieved p ≤ 0.9 and were considered for multivariate analysis, details on these results are given in Table 2.

Associations between cART regimens and DRCs after adjusting for biomedical and demographic factors of recipients are presented in Table 3.

HIV patients on second- and third-line cART (cART line2/3) were 47% less likely to develop DRCs earlier than those on first-line cART [AHR = 0.53; (95% CI: 0.33–0.87)]. Recipients of cART who had a CD4–1 count equal to or less than 200 cells/mm³ at cART initiation were 68% more likely to develop DRCs earlier than those who had a CD4–1 count greater than 200 cells/mm³ at initiation of cART [AHR = 1.68; (95% CI: 1.03–2.73)]. Participants who reported adherence to cART had a 3.1 times higher risk of developing DRCs earlier than those who did not adhere [AHR = 3.1; (95% CI: 0.92–10.0), p = 0.06]. Overweight patients at initiation of cART, had a 63% higher risk of developing DRCs earlier than recipients who had normal BMI at cART initiation [AHR = 1.63; (95% CI: 1.07–2.64)]. Males had a 51% lower hazard of the outcome than females [AHR = 0.49; (95% CI: 0.32–0.77)]. Age as a continuous variable did not enter the multivariate model as it did not achieve a p < 0.05.

The risk of new onset DRCs increased each year with cumulative exposure to each cART regimen under study. However, the most significant risk was documented among recipients of the first-line cART regimen while the least risk of new DRC cases per year of cumulative exposure to cART was recorded among recipients of the second-line cART regimen (Fig. 1).

Hazard function for recipients of cART first-line, second-line and third-line: a 12 years retrospective cohort study in two HIV Clinics in Botswana. cART line1 = cART first-line regimen, cART line2 = cART second-line regimen, cART line3 = cART thrid-line regimen

Full size image

The survival function of patients on second- and third-line regimens almost overlapped, with both functions disassociated with the first-line function. Data presented in Fig. 2 show survival functions of aggregated recipients of second- and third-line cART regimens in comparison with those of the first-line cART regimen.

Survival function for recipients of cART first-line, and second-, third-line regimens: a 12 years retrospective cohort study in two HIV Clinics in Botswana. cART 2/3 = cART second-, and third-line regimen, cART line 1 = cART first-line regimen, log rankX² = log rank chi square

Full size image

The log rank (Mantel Cox) test confirmed that recipients of second-, third-line cART had significantly longer survival to DRC development than recipients of first-line cART (X² = 8.97, p = 0.003).

Findings from this study bring new insights to the fight against HIV and related comorbidities. The risk of new DRC onset per year of cumulative exposure to cART was higher among recipients of first-line cART than in recipients of cART second-line and cART third-line regimens. The longest survival/time duration to the development of DRC per year of cumulative exposure to cART was documented among recipients of second- and third-line cART regimens. In other words recipients of the second-, and third-line cART were less likely to develop DRCs earlier than recipients of the first-line cART. This lends support to the claim that the survival of cART recipients to DRC development is not solely a function of the duration that patients are subjected to these drugs but also a function of the specific regimen they are exposed to [19, 20] and the recipients underlying biomedical and demographic characteristics [21]. While the overlapping of the second-, third-line regimen DRC hazard/survival functions may suggest comparable effects between the two cART regimens, they also indicate that recipients of both the second and third-line regimen had a longer survival probability compared to their counterparts on first-line regimen. This implies that had some patients been directly initiated on second- or third-line cART regimen, the outcome might have been much better for all the recipients.

In short, important findings in this study were that recipients of the second-, third-line cART regimen were less likely to develop DRCs earlier than first-line cART regimen recipients. This is good news that may indicate that different cART regimens are differently associated with DRCs and different groups of recipients. Thus, depending on circumstances, the current second- and third-line regimens may be considered for use as first-line regimen while taking into account factors such as recipients’ pre-existing DRC risk factors, drug availability, cost-effectiveness, etc. For instance, patients with characteristics such as overweight, female gender or CD4–1 equal or less than 200 cells/mm³ were identified in this study as at more risk of developing DRCs earlier than those without and should have been put directly on second- or third-line regimen rather than starting with the first-line treatment. Recipients aged 35 years or older also need particular attention compared to those aged less than 35 years before being assigned to a specific cART regimen.

However, the challenge lies in the cost difference between the first-line, the second- and third-line regimens and possible adverse effects of the third-line regimen. While WHO states that treatment programs using Tenofovir or AZT are promising, they may not be the ultimate solution to alleviate the situation since both drugs are known to be expensive [22] and yet to some extent associated with DRCs [6, 7, 23,24,25]. The fact that despite their high cost, these drugs are still associated to a certain degree with DRCs [23,24,25] calls for a meticulous drug selection that takes into account biomedical and demographic characteristics of recipients so as to improve their survival probability to the event.

As with any other study, this study had some limitations. It was a retrospective review of data and may be missing some important DRC risk factors like physical activity, cigarette smoking, alcohol consumption and family history of DRCs. Prospective studies are needed to address this particular issue. Even so, the study has provided useful evidence to inform policy and decision making on alternative approaches for the management and care of HIV patients. The study has also indicated that patients with underlying DRC risk factors such as overweight, female gender, aged 35 years or older or those with CD4 ≤ 200 cells/mm³ should be selectively assigned to treatment regimens for possible improvement of their survival.

The study has shown that the risk of new DRC onset among cART recipients is not only a function of the type of cART regimen or the duration of exposure to these drugs but also a function of patients’ underlying biomedical and demographic DRC risk factors. The study has also presented a survival model highlighting DRCs’ significant prognostic factors to guide clinical care, policy and management of recipients of cARTs. Further studies in the same direction will likely improve the survival to the development of DRC of every cART recipient in this community.

AHR:

Adjusted hazard ratio

AIDS:

Acquired Immune deficiency Syndrome

AZT:

Azidothymidine, also known as zidovudine

BMI:

Body mass index

Bontleng HIV clinic:

one of the study site

cART:

Combination Antiretroviral Therapy

CD4:

Cluster of differentiation 4, glycoprotein found on the surface of T helper cells

CD4–1:

CD4 Count at cART initiation

CD4–2:

CD4 Count at the time of data collection

CI:

Confidence Interval

DRCs:

Diabetes Mellitus-Related Comorbidities

HIV:

Human Immunodeficiency Virus

ICD-10:

The 10th revision of the International Statistical Classification of Diseases and Related Health Problems, a medical classification list by the World Health Organization

Log-rank X² :

Log-rank Chi-square

PLHIV:

People Living with HIV

PMH HIV Clinic:

Princess Marina Hospital HIV Clinic or another study site

SE:

Standard Error

SEM:

Standard Error of the Mean

Survival:

Survival was calculated as the time elapsed from the date of cART initiation to the date of the development of the first DRC at the time of data collection

Time to the event:

Time duration from initiation on cART to the occurrence of DRCs

UHR:

Unadjusted hazard ratio

Weight-1:

Weight in kilograms at cART initiation;

Weight-2:

Weight in kilograms at the time of data collection

WHO:

World Health Organization.

We thank PMH and Bontleng HIV clinic staff for allowing us to collect data at their facilities. We would also like to thank Mr. Kadimo, a librarian at the University of Botswana, for assisting us with literature search.

Funding

The study was financially supported by the Health Resources and Services Administration (HRSA) of the U.S. Department of Health and Human Services (HHS) under grant T84HA22125 (Medical Education Partnership Initiative, $9,400,000). This information or content and conclusions of this study are those of the authors and should not be construed as the official position or policy of, nor should any endorsements be inferred by HRSA, HHS or the US Government.

Availability of data and materials

Data underlying the findings in this study are not publicly available so as to maintain patient confidentiality. The data include potentially identifying demographic and clinical care information. However, the data can be requested from the corresponding author who must first get permission from the management of the HIV clinics where the study was conducted before sharing.

Affiliations

1. Department of Family Medicine and Public Health, Faculty of Medicine, University of Botswana, Gaborone, Botswana

   • Jose Gaby Tshikuka
   • , Goabaone Rankgoane-Pono
   • , Mgaywa Gilbert Mjungu Damas Magafu
   • , Tiny Masupe
   • , Mooketsi Molefi
   • , John Thato Tlhakanelo
   • , Shimeles Genna Hamda
   •  & Vincent Setlhare

2. Department of Health Sciences, National Pedagogic University, Kinshasa I, Democratic Republic of the Congo

   • Jose Gaby Tshikuka
   •  & Maurice Nsikungu-Kalukul

3. Department of Global Health, University of Washington, Seattle, USA

   • Mgaywa Gilbert Mjungu Damas Magafu

Contributions

This study was designed by JGT and GRP. JGT and GRP participated in all stages of this study from proposal preparation to manuscript writing. MGMDM participated in reviewing the study proposal, analysing data and writing the manuscript. TM, MM, SGH, VS, MNK, JTT participated in data analysis and manuscript preparation. All the authors contributed to the writing of the manuscript. All the authors read and approved the final manuscript.

Corresponding author

Correspondence to Jose Gaby Tshikuka.

Ethics approval and consent to participate

Ethical approval to collect data from HIV clinics was sought and obtained from the University of Botswana Review Board and the ethics committee of the Ministry of Health and Wellness, Botswana. Permission to consult clinic record books and systems was also sought and obtained from the clinic management. As this was a record based study, no consent to participate was required.

Consent for publication

Not Applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Reprints and Permissions