Evaluation of antiretroviral effect on mitochondrial DNA depletion among HIV-infected patients in Bali

Abstract

Background

Nucleoside reverse transcriptase inhibitors (NRTIs) are the cornerstone of highly active antiretroviral therapy combination regimens for HIV infection. Unfortunately, NRTIs have been noticeably associated with many adverse effects related to mitochondrial toxicity leading to mitochondrial deoxyribonucleic acid (mtDNA) depletion. However, similar mitochondrial dysfunction has recently been found even in antiretroviral therapy-naïve patients, suggesting HIV itself could contribute to this abnormality. In this study, we determine whether mtDNA depletion was present in either antiretroviral therapy-naïve or NRTI-treated patients at Sanjiwani Hospital, Bali, Indonesia.

Patients and methods

A cross-sectional study was conducted from the peripheral blood mononuclear cells of HIV patients. Specifically, the relative content of mtDNA (mtRNR1 gene) to nuclear DNA (ASPOLG gene) was determined by real-time polymerase chain reaction. Data were analyzed with SPSS 16.0 software and GraphPad Prism 7.02.

Results

A total of 84 samples (67 on NRTIs and 17 HIV-naïve) were suitable for analysis. We identified 21.4% of the samples (18/84) with mtDNA:nDNA ratio <1. Although it was not significant (P=0.121), the median mtDNA:nDNA ratio of HIV-naïve group was slightly higher (median 1.8; interquartile range [IQR]: 1.1–2.1) than NRTI-treated patients (median 1.5; IQR: 1.3–2.85). Tenofovir-based NRTI was more frequently used (73.13%) than zidovudine-based NRTI (26.86%). The period for which NRTI was used probably contributed to the ratio of mtDNA:nDNA. The median ratio of mtDNA:nDNA zidovudine-treated patients was slightly lower (median 1.2; IQR: 1.08–1.98) when compared to tenofovir-based NRTI (median 1.6; IQR: 1.05–2.10), with the median period of former treatment being significantly longer (P<0.001). Although these data overall indicate that NRTI treatment had no effect on mtDNA:nDNA ratios, patients who undergo more than 12 months of NRTIs treatment show a decrease in the ratio; however, further study is required.

Conclusion

Almost one-fourth of the samples showed a lower mtDNA:nDNA ratio. The decreasing of the ratio mtDNA:nDNA was most likely present after 12 months of NRTI treatment.

Keyword:

• mitochondrial
• nuclear DNA ratio
• HIV
• ART

Introduction

Antiretroviral therapy (ART) has changed HIV infection from a lethal disease to a chronic, manageable illness.¹ ART has improved the long-term prognosis and markedly increased the life expectancy of individuals living with HIV. Unfortunately, the increase in life expectancy was found to be associated with several metabolic complications such as diabetes, hypertension, and dyslipidemia,^2,3 or low bone mineral density.^4,5 ART, which consists of a combination of 2 or more reverse transcriptase inhibitors (RTIs) and protease inhibitors, should be used over long periods of time. The first-line ART regimens recommended by the Ministry of Health of Indonesia include the combination of zidovudine (ZDV or AZT), lamivudine (3TC), nevirapine (NVP), or efavirenz (EFV) or a fixed-dose combination consisting of tenofovir (TDF), 3TC or emtricitabine, and EFV. The ZDV-based combination consists of ZDV, 3TC, and NVP/EFV. Meanwhile, the essential components of fixed-dose combination are TDF, 3TC/emtricitabine, and EFV. Second-line ART consists of the same drugs with AZT replaced by one of the protease inhibitors, mostly ritonavir boosted.⁶ Long-term use of ART combination that uses nucleoside RTIs (NRTIs) as the main drug may result in adverse events, such as lipodystrophy syndromes, increased risk of cardiovascular diseases,^7,8 and the life-threatening lactic acidosis.^9–13 NRTIs work by being incorporated into growing viral DNA chain during reverse transcription, which results in chain termination.¹² The NRTIs are also γ polymerase inhibitors, which inhibit mitochondrial DNA (mtDNA) synthesis, and therefore many of those adverse effects are associated with mitochondrial toxicity, specifically mtDNA depletion.^12,14,15 In contrast with other study results, second-line ART, particularly atazanavir/ritonavir-based ART, has been reported to have no favorable effect on plasma lipid profile or adipose tissue distribution.¹⁶

The clinical symptoms of mtDNA depletion vary, from mild peripheral neuropathy to life-threatening severe lactic acidosis.¹⁷ The depletion of mtDNA may present prior to the clinical symptoms, and therefore it has been used as a noninvasive tool to evaluate mitochondrial toxicity in HIV-infected patients on ART.¹⁸ The gold standard to assess mitochondrial disorders during HIV infection is a biopsy of muscle or adipose tissue of the heart or liver, or general adipose tissues.^19–21 A less invasive method was developed to assess mitochondrial toxicity by using peripheral blood to monitor mtDNA depletion in HIV-infected individuals on ART.^12,22,23 Additionally, quantification of mtDNA in the peripheral blood mononuclear cells (PBMCs), as measured by quantitative real-time polymerase chain reaction (qRT-PCR), is a relatively inexpensive method to monitor mtDNA depletion in HIV-infected patients.²⁴

In this study, we presented mtDNA:nDNA ratio of HIV-naïve and NRTI-treated patients who visited Sanjiwani Hospital. We also evaluated whether the stage of HIV and period of treatment affect mtDNA:nDNA ratio.

Patients and methods

Sample collection

Ninety-eight patients infected with HIV were recruited from Sinta HIV Clinic of Sanjiwani Hospital at Gianyar Regency and Puskesmas Denpasar Selatan II Sanur from the Denpasar municipality in Bali, Indonesia. Eight milliliters of whole blood samples were collected from the study participants during the period of January to March 2016. PBMCs were isolated by density gradient centrifugation for 10 min at 1,800 rpm using the BD Vacutainer® CPT™ Cell Preparation Tube (Becton Dickinson and Company, Franklin Lakes, NJ, USA). Total DNA was extracted from PBMCs using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). Data on demographic and clinical features as well as disease severity, according to the World Health Organization guidelines, of the study participants, were retrieved from the medical record.

RT-PCR assay

The quantification of mtDNA and nDNA was performed using the CFX96 Touch Real-Time Detection System (Bio-Rad Laboratories, Hercules, CA, USA). The DNA concentration and its purity were determined by measuring the absorbance at 260 nm (A260) and 280 nm (A280) in a spectrophotometer (Nano drop; Thermo Scientific, Waltham, MA, USA) prior to the RT-PCR assay. The RT-PCR assay was validated by assessing the optimal range of primer annealing temperatures, reaction efficiency, and specificity using a standard set of samples.²⁵ The reaction was set up using EVA Green SMX 200R kit (Bio-Rad Laboratories). To evaluate mtDNA depletion, the relative content of mtDNA fragment (mtRNR1 gene) was compared with the nuclear DNA (ASPOLG gene). The reaction contained 5 µL of master mix (2×), 2 µL DNA template (2 ng/µL), 1 µL forward primer (10 µM), 1 µL reverse primer (10 µM), and PCR-grade water to a final volume of 20 mL. The primers of mtRNR1 were as follows: sense 5′-TAGCCCTAAACCTCAACAGT-3′ and mtRNR1 antisense 5′-TGCGCTTACTTTGTAGCCTTCAT-3′;²⁶ the primers for ASPOLG were as follows: sense 5′-GAGCTGTTGACGGAAAGGAG-3′ and ASPOLG antisense 5′-CAGAAGAGAATCCCGGCTAAG-3′. The RT-PCR conditions for mtDNA and nDNA amplification were as follows: one cycle of 3 min at 95°C for denaturation; 40 cycles of 10 s at 95°C for denaturation, 30 s at 64.5°C for annealing, and 1 min at 72°C for extension. Melt curve analysis consisted of 64.5°C–95°C for 5 s with increment 0.5°C. GADPH gene was utilized as a housekeeping gene for normalizing RNA expression in qRT-PCR.

Data and statistical analyses

For descriptive analyses, number and percentage were used for categorical variables and median with interquartile range (IQR) were used for continuous variables. Nonparametric statistical analysis was performed by using Mann–Whitney rank-sum test and Kruskal–Wallis H test. A P-value <0.05 was considered statistically significant. All analyses and graph creations were performed using SPSS 16.0 software (SPSS Inc., Chicago, IL, USA) and GraphPad Prism 7.02 (GraphPad Software, La Jolla, CA, USA).

This study was conducted with approval from the Medical Research Ethics Committees of Faculty of Medicine, Udayana University, Bali, Indonesia (document number 115/UN.14.2/Litbang/2015). All participants enrolled for this study provided written informed consent prior to their participation.

Results

Clinical data

The clinical and epidemiological characteristics of the study participants and therapeutic data are summarized in Table 1. Eighty-four of a total of 98 samples were eligible for analysis. Almost 90% of the samples were Balinese, and more than 50% practiced unsafe sexual intercourse with multiple partners. The 84 eligible samples consisted of 17 patients who were NRTI-naïve and 67 who were on ART. Most of the participants had full-blown AIDS (76.1%). In addition, 51 of 84 individuals (60%) had contracted opportunistic infections. Thirty-three percent of the samples were on TDF-based ART. Median CD4⁺ baseline was 53 cells/µL (1–480 cells/µL). The mean duration of TDF-based ART was significantly shorter than ZDV-based ART in the participants (P<0.001). We could not evaluate the clinical symptoms of mitochondrial toxicity, such as varying degrees of fatigue, because a standard tool to quantitatively measure the symptoms was not available. Other clinical symptoms, including rapid onset of weight loss, were also not evaluated. Only a minority of the patients presented with unclear numbness (data not shown).

Table 1 Demographical data and mtDNA:nDNA ratio of the on ART and ART-naïve samples

Characteristics	On ART (N=67)	ART-naïve (N=17)	P-value*
Age (years)	35 (29–42)	30 (27–34)	0.012
Sex
Male (f, %)	44 (66%)	7 (41%)
Female (f, %)	23 (34%)	10 (59%)
Duration of HIV diagnosis (months)	22 (8–36)	1 (1–1)	<0.001
mtDNA:nDNA ratio	1.5 (1.10–2.10)	1.8 (1.30–2.85)	0.121

Notes: Data are reported as median (interquartile range).

* Mann–Whitney rank-sum test.

Abbreviations: ART, antiretroviral therapy; mtDNA, mitochondrial DNA.

The median of treatment duration roughly showed a tendency of decreasing mtDNA:nDNA ratio after 1 year of ART treatment. However, a Kruskal–Wallis H test showed that this was not a statistically significant decrease, χ²(4) =2.38, P=0.66, with a mean rank mtDNA:nDNA ratio of 34.95 after 1 year or less of treatment, 39.23 following 1–2 years treatment, 34.53 with 2–3 years of treatment, 30.45 after 3–4 years of treatment, and 25.33 following treatment >4 years. Interval treatment duration of ART was 1=0–12 month, 2=12.1–24 month, 3=24.1–36 months, 4=36.1–48 months, 5=48.1–50 months.

mtDNA:nDNA ratio between 3 different groups

Overall, most of the participants (78.6%, 66/84) had ratios more than 1 (Figure 1), and most of those samples were from the ART-naïve group. Our data showed that 2.35% (2/84) of the NRTI-naïve samples had mtDNA:nDNA ratio <1. As expected, the median mtDNA:nDNA ratio was higher in NRTI-naïve samples (Table 1) than NRTI-treated samples. Although the median among the groups was not statistically significant, there was a tendency toward a lower ratio in the NRTI-treated groups compared to the NRTI-naïve group.

Figure 1 Distribution mtDNA:nDNA ratio by duration of treatment.

Abbreviation: mtDNA, mitochondrial DNA.

Among 3 different groups – NRTI-naïve, ZDV-based, and TDF-based, the ZDV-based group had the lowest mtDNA:nDNA ratio.

Additional statistical analyses were carried out to determine the profile of ZDV- and TDF-based treatment (Table 2). The number of patients who underwent ZDV-based ART was 18/67 (26.9%). This group had a lower mtDNA:nDNA ratio than the TDF-based group, and this might be as a result of the median treatment duration, which was significantly longer at 41 months (IQR 31–47) versus 16 months (IQR 7–30). Although these data indicate that NRTI treatment had no effect on mtDNA:nDNA ratios, the median of treatment duration roughly showed a tendency of decreasing mtDNA:nDNA ratio after 1 year of ART treatment (Figure 1).

Table 2 Profile comparison of ZDV-based and TDF-based antiretroviral type

Profile	ZDV-based ART (N=18)	TDF-based ART (N=49)	P-value*
Duration of ART (months)	41 (31–47)	16 (7–30)	<0.001
HIV stage
Non-AIDS (f, %)	2 (11%)	14 (29%)
AIDS (f, %)	16 (89%)	35 (71%)
mtDNA:nDNA ratio	1.2 (1.08–1.98)	1.6 (1.05–2.10)	0.220

Notes: Data are reported as median (interquartile range).

* Mann–Whitney rank-sum test.

Abbreviations: ART, antiretroviral therapy; mtDNA, mitochondrial DNA; TDF, tenofovir; ZDV, zidovudine.

Discussion

Some in vivo and in vitro studies indicated the ratio of mtDNA:nDNA was significantly lower in ART-treated HIV participants,²² and this treatment was associated with symptomatic hyperlactatemia. A study of mtDNA:nDNA ratio in adipose tissue also showed significantly lower in ART-treated HIV patients.¹⁵ In addition, an in vitro study in primary human lymphocytes stimulated with NRTIs demonstrated a lower mtDNA level on day 10 of exposure, followed by increasing lactic acid production.²⁴ Mitochondrial AC/TG mutations were found in a prospective cohort study in the blood of a newborn from a mother who had been exposed to ART during pregnancy.²³ The reason for this discrepancy may be related to different study methods, including the technique used to evaluate mitochondrial toxicity in HIV-treated participants. A study to determine the effect of HIV infection and ART on placental mitochondria with qRT-PCR also found the ratio to be significantly reduced in ART and HIV-1 exposed placentas in comparison with uninfected controls.²⁷ This study, which found a significant effect of combined ART containing ZDV to mtDNA depletion on placentas of HIV-infected pregnant women, showed increased oxidative stress level, apoptosis suggestive of secondary mitochondrial failure, and a potential basis of associated adverse perinatal outcome.²⁸

In contrast with the aforementioned studies, our result indicated that there was no significant difference between mtDNA:nDNA ratio of HIV-infected patients who were on ART or in the ART-naïve groups. This finding supports some studies which showed that NRTI treatment had no effect on mtDNA and mitochondrial number.¹⁰ In this study, we also found the ratio of mtDNA:nDNA was <1 in a small number of ART-naïve patients. These findings may indicate that HIV itself influences the mtDNA contents and toxicity, although the underlying mechanism of this remains to be established. This could be because mitochondria may have been exposed chronically to reactive oxygen species that is released in chronic HIV inflammation.²⁹

Furthermore, the HIV-proteins, HIV-Tat protein,³⁰ and Vpr protein³¹ have been shown to influence mitochondrial integrity. Casula et al³² found in a cross-sectional study that mtDNA reduces 1 year after seroconversion of HIV infection, while no significant decrease was observed during the subsequent 4 years with ART. Another study which was conducted by Maagaard et al³³ reported that there was no significant difference between depletion of mtDNA on HIV + ART and HIV + ART-naïve respondents. Based on their result, they concluded that there is no correlation between the mtDNA:nDNA ratio in PBMC and muscle.³³ Our study did not evaluate the clinical symptoms of mitochondrial toxicity because there was no accurate tool available to define the associated symptoms.

Limitations

The major limitation of this study is the sample size, particularly the HIV-ART-naïve group, which was probably too small. Clinical symptoms of mitochondrial toxicity such as varying degrees of fatigue could not be evaluated because a standard tool to quantitatively measure the symptoms was not available. Other clinical symptoms including rapid onset of weight loss were not evaluated.

Conclusion

This study reported that there was a depletion pattern in mtDNA:nDNA ratio after 12 months of ART treatment. Longer duration of ART treatment may further inhibit DNA polymerase-γ; therefore, mitochondrial toxicities were more likely to develop in the patients. Insignificant reduction of mtDNA in PBMCs with ART treatment may not reflect the mtDNA ratio in tissue. Because mitochondrial toxicity is tissue specific, mtDNA:nDNA ratio on PBMCs may not represent the real condition of mitochondrial toxicity in a patient.

Author contributions

The study was conceptualized by SM. The data was curated by DGB. Laboratory support was provided by SQK, NN, and IA. Formal analysis was carried out by HS, SAA, and DM. SM and DM aided with manuscript writing. All authors contributed toward data analysis, drafting and critically revising the paper and agree to be accountable for all aspects of the work.

Acknowledgments

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. This research was supported by the Faculty of Medicine and Health Science, Warmadewa University. We would like to express our gratitude to the patients for their participation in this study. We particularly thank the medical doctors, nurses, and the staff in Sanjiwani Hospital Gianyar, Bali, for their great support.

Disclosure

The authors report no conflicts of interest in this work.