Abstract
The present study aimed to evaluate the role of genetic polymorphisms in the emergence of lipoatrophy or lipodystrophy in HIV-infected patients with antiretroviral therapy (ART) in Thailand. Position 455 upstream of the Apolipoprotein C3 gene (ApoC3 T–455C, rs2854116), codon 64 of the Beta3 adrenergic receptor gene (ARβ3 Tcod64C, rs4994), and position 670 upstream of the Fas gene (Fas A–670G, rs1800682) were genotyped in 829 HIV-infected Thai patients who had started ART. Crude and adjusted incidence rate ratios (IRR) were calculated using Poisson regression. The serum levels of cholesterol, triglycerides, high-density lipoprotein (HDL), and low-density lipoprotein (LDL) were also analyzed. Multivariate analysis revealed an association between the Fas −670AA genotype, but not the ApoC3 −455 or ARβ3 cod64 genotypes, with the incidence of lipoatrophy after adjusting for gender and stavudine (d4T)-containing regimens (IRR=1.72, 95% CI=1.20–2.45, p=0.003). However, ApoC3 −455C homozygous patients showed elevated serum levels of triglycerides, while this genotype did not affect serum total cholesterol, HDL, or LDL levels in patients with lipoatrophy or lipodystrophy. In contrast, the ARβ3 cod64 genotype did not show any significant association with the serum levels of cholesterol, triglycerides, HDL, or LDL. In conclusion, Fas −670AA affected the incidence of lipoatrophy in HIV-1-infected Thai patients, while the ApoC3 −455C allele affected the serum levels of triglycerides. These results confirmed the role of genetics in the development of ART-related metabolic disorders.
Introduction
M
Morphological alterations are collectively named lipodystrophy (LD); however, in recent years it has been recognized that fat loss and fat gain represent distinct entities. Clinical features of fat loss include loss of subcutaneous fat in the face, arms, legs, and buttocks, while fat accumulation is characterized by excess fat deposition in visceral adipose tissue in the abdomen, breasts, dorsocervical region, and trunk. 4
The underlying mechanism of antiretroviral therapy (ART)-related lipoatrophy (LA) has not been fully elucidated. Many mechanisms have been suggested, including mitochondrial toxicity by nucleoside reverse transcriptase inhibitors (NRTIs), inhibition of adipocyte differentiation by protease inhibitors (PIs), increased levels of inflammatory cytokines, and HIV itself. 5 –7 Mitochondrial toxicity has been clearly shown to play a pivotal role; in particular, thymidine analogs are strong inhibitors of DNA polymerase-γ and lead to mitochondrial depletion and dysfunction 8 ; mainly stavudine (d4T), but also zidovudine (AZT) and didanosine (ddI) are strongly related to LA. 9,10 In previous in vitro and clinical studies, newer nucleoside and nucleotide agents, such as lamivudine (3TC), emtricitabine, abacavir, and tenofovir, appear to be much weaker inhibitors of mitochondrial DNA polymerase-γ or other mitochondrial functions, and appear to be associated with a lower risk of events thought to be related to mitochondrial toxicity. 11 However, a certain degree of mitochondrial toxicity is still present and will probably still affect HIV patients under ART, although to a lesser extent. 12 –14 In fact, switch studies from d4T to other NRTIs and complete switch off of thymidine analogs showed modest even if consistent increases in limb fat. 15 –17
The pathogenesis of lipohypertrophy appears to be multifactorial, with age, sex, HIV itself, immune depression, and duration and type of ART related to its appearance. 18 –20 PIs were suggested to be involved in lipohypertrophy; however, visceral fat accumulation also occurs in the absence of PIs. 21,22 A compartment-specific effect of mitochondrial toxicity within the adipose tissue may be related to lipohypertrophy 23 as well as dysregulation of fatty acid metabolism and altered expression of adipokines. 24 –26
However, LA and lipohypertrophy do not occur in all treated patients, and there is a very large degree of interindividual variability in the timing of emergence and severity of symptoms. A study of identical twins showed that genetic factors are involved in the accumulation of visceral adipose tissue. 27 These data suggest that host genetic factors may play a role, and inherited predispositions may have a significant influence on the appearance of LD and metabolic alterations as well as on the viroimmunological response to the drugs.
Therefore, host genetic polymorphisms may have a significant effect on the response to ART in terms of metabolism and immune response. 28 A study of Italian HIV patients suggested that single nucleotide polymorphisms (SNPs) in genes involved in apoptosis and lipid metabolism mediate the development of LA. 29 In particular, the TT and CT genotypes at position −455 of the Apolipoprotein C3 (ApoC3) gene (ApoC3 T–455C, rs2854116), nonsynonymous (W to R) CC and CT genotypes at codon 64 of the Beta3 adrenergic receptor (ARβ3) gene (ARβ3 Tcod64C, rs4994), and the AA genotype at position −670 of the Fas gene (Fas A–670G, rs1800682) were associated with increased incidence rates of LA. 29 ApoC3 protein, localized mainly in very low-density lipoprotein and high-density lipoprotein (HDL), is involved in fat metabolism and may delay the catabolism of triglyceride-rich particles. ApoC3 inhibits lipoprotein lipase and hepatic lipase and decreases the chylomicrons in hepatic cells. 30,31 ARβ3 protein, expressed mainly in visceral adipose tissue, 32,33 contributes to both lipolysis and the delivery of free fatty acids. Fas is the main gene that controls cell death by inducing apoptosis. A previous study indicated that Fas and Fas ligand influenced the immune reconstitution induced by ART. 34
In Thailand, the National Access to Antiretroviral Care program started in 2002 using standard regimens with GPO-VIR (a combination of d4T, 3TC, and nevirapine) as the first line of ART. The present study was performed to determine whether these genetic polymorphisms found in the Italian cohort also mediate the emergence of LA or LD in HIV-infected patients starting ART in Thailand in consideration of the relevance of these side effects in the currently used regimens and of the differences in ethnicity, alimentation, and way of life of the Thai population.
Materials and Methods
Study design
HIV-1-infected individuals receiving ART at Bamrasnaradura Infectious Diseases Institute, Nonthaburi, Thailand were enrolled in this study. LA was defined as the presence of fat loss at any site of the body (face, arms, legs, or buttocks). LD was defined as LA combined with fat accumulation diagnosed at the same time. Fat accumulation was defined as the presence of at least one alteration, such as buffalo hump, or fat accumulation in the breast, abdomen, neck, or lipomas. To qualify, the alteration had to be recognized by both the patient and the physician. 35 Most patients who developed LA or LD changed regimen for those without d4T. A 200 μl sample of whole blood was collected from each patient and stored at −20°C until DNA extraction. The genotypes of ApoC3 −455, ARβ3 cod64, and Fas −670 were compared between patients who developed LA or LD and those who did not. Patients' age, gender, transmission route, CD4 counts at pre-ART and at 1 year post-ART (cells/μl), drug regimens, and duration of therapy were obtained from their medical records. Lipid profile data were obtained from patients' records of serum biochemical tests. Blood for this test was collected when patients were in the fasting state; however, not all patients had serum biochemical test data as the test incurs extra costs for patients. All participants signed an informed consent form. The study was approved by the institutional ethical committees at the Bamrasnaradura Infectious Diseases Institute and Department of Disease Control, Ministry of Public Health, Thailand.
DNA extraction
Genomic DNA was extracted from 200 μl of each whole blood sample using a QIAamp DNA Blood Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions. 36
Genotyping
Genotyping of Fas gene promoter polymorphism, A–670G (rs1800682)
Polymerase chain reaction (PCR) and restriction fragment length polymorphism analysis for Fas A–670G, modified from a previously published technique, 37 was conducted with the primers 5′-CTACCTAAGAGCTATCTACCGTTC-3′ and 5′-GGCTGTCCATGTTGTGGCTGC-3′. The PCR product was digested with the restriction enzyme MvaI (Roche Biochemicals, Mannheim, Germany) at 37°C overnight. Two polymorphic alleles could be distinguished on the basis of restriction length polymorphisms (188 bp for G and 232 bp for A).
Genotyping of ApoC3 gene promoter polymorphism, T–455C (rs2854116)
TaqMan SNP assay was performed using the StepOne real-time PCR system (Applied Biosystems, Foster City, CA) with primers (5′-GAGCTCAGCCCCTGTAACCAG-3′ and 5′-ACACAGCCTGGAGTAGAGGG-3′) and TaqMan MGB probes (5′-VIC-CTCCAAACA
Genotyping of ARβ3 gene coding polymorphism, Tcod64C (rs4994)
TaqMan SNP assay was performed with a TaqMan SNP assay kit (C___2215549_20; Applied Biosystems).
Statistical analysis
Person-year analysis was conducted to investigate a single endpoint: the incidence of LA or LD. If the patient was diagnosed as both LA and LD, the earlier date of the first diagnosis was used as the endpoint. We focused on the time to the occurrence of the first event. Person-years at risk were calculated from the start date of ART until the last available follow-up or development of the event, whichever occurred first. The Poisson regression univariate model was used to determine whether the genetic polymorphisms as well as patients' age, gender, transmission route, CD4 counts at pre-ART and at 1 year post-ART, and starting drug regimes were predictors of LA or LD. Associations were expressed as incidence rate ratio (IRR) and corresponding 95% confidence interval (CI). Factors showing significant association with LA or LD were then included as covariates in the Poisson regression multivariate model to adjust for their possible confounding effects. Starting drug regimens were grouped into either d4T-containing regimens or other backbone regimens. Continuous variables were compared using the Wilcoxon signed rank test for paired samples and Kruskal–Wallis test for unpaired samples. A variable was considered significant at p<0.05.
Results
Incidences of LA and/or LD in HIV-1-infected Thai patients undergoing ART
A total of 829 HIV-1-infected individuals who had started ART were enrolled in this study at the Bamrasnaradura Infectious Diseases Institute, Nonthaburi, Thailand, from April 2008 to August 2009. Patient demographic data are shown in Table 1. The study population consisted of 478 males (57.7%) and 351 females (42.3%). Most of the patients contracted infection through sexual contact. Almost half of the patients started ART immediately after being diagnosed with HIV infection. Median CD4 count per microliter pre-ART was 45 and that post-ART (after 1 year) was 210.5. The mean follow-up period was 4.85 years (SD=2.42).
SD, standard deviation; IQR, interquartile range.
We first combined patients with LA and those with LD for analysis, as fat loss was commonly observed in both LA and LD cases. Among the 829 patients, 270 (32.6%) were diagnosed with LA or LD within the observation period. We kept on monitoring the other 559 (67.4%) until the end of the observation period and found that they developed neither LA nor LD (non-LA/LD cases). The median duration of therapy until diagnosis was 2.63 years for the LA or LD cases and 5.53 years for the non-LA/LD cases. In a follow-up of 4019.8 person-years, the incidence rate of LA or LD was 6.70 per 100 person-years. LA or LD cases were found more frequently in females (40.7%) than in males (26.6%), and the incidence rate per 100 person-years for females (9.11) was significantly higher than that for males (5.18) (IRR=1.76, 95% CI=1.37–2.25, p<0.001) (Table 2). Patients using d4T/3TC backbone regimens as starting drugs showed a significantly higher incidence rate per 100 person-years (7.71) than those with other non-d4T-containing (mostly AZT/3TC) regimens (3.49) (IRR=2.21, 95% CI=1.53–3.28, p<0.001) (Table 2). These findings were consistent with those of previous studies in other ethnic groups. 9,20,29,38,39 In contrast, age, transmission route, and CD4 counts pre-ART and at 1 year post-ART did not significantly affect the incidence of LA or LD (Table 2).
Diagnosis of lipoatrophy or lipodystrophy.
Gender, NRTIs, and Fas−670 were included as covariates.
Including eight intravenous drug users.
No-d4T group mainly used AZT. One case had received ddI instead of AZT as an initial regimen.
CD4 cells/μl.
IRR, incidence rate ratios; NRTIs, nucleoside reverse transcriptase inhibitors; d4T, stavudine.
Fas −670AA, but not ApoC3 −455 or ARβ3 cod64 genotype, affects incidence of LA or LD in HIV-1-infected Thai patients
To elucidate the possible effects of genetic polymorphisms on the development of LA or LD, we determined the genotypes of Fas −670, ApoC3 −455, and ARβ3 cod64 of 829 HIV-1-infected Thai patients. LA or LD cases were found significantly more frequently among patients with the Fas −670AA genotype than in those with other Fas −670 genotypes (IRR=1.41, 95% CI=1.06–1.86, p=0.015, Table 2). The Fas −670AA genotype was found in 25.9% of the LA or LD cases and in 20.2% of the non-LA/LD cases during the observation period. Increased frequency of the Fas −670AA genotype in LA or LD cases was observed in both males and females (data not shown). In cases of ApoC3 −455 and ARβ3 cod64 genotypes, however, there were no statistically significant differences in the frequency of LA or LD between groups (Table 2). Adjustment for gender and starting drug regimens in multivariate Poisson regression analysis resulted in a slightly higher IRR for Fas −670AA (1.47, 95% CI=1.12-194, p=0.005, Table 2). This difference was statistically significant after correction for multiple testing of three SNPs (p=0.015).
Fas–670AA genotype affects incidence of LA but not that of LD in HIV-1-infected Thai patients
Among 270 LA and/or LD cases, 148 were diagnosed as LA singly and 105 were diagnosed as LD singly. There were 17 cases diagnosed as LA first who then subsequently developed fat accumulation. These 17 cases were grouped as LD cases. The median durations of therapy until diagnosis were 2.58 years for LA and 2.83 years for LD cases. When we analyzed LA and LD separately, females showed a higher IRR in both LA (1.60, 95% CI=1.14–2.24, p=0.005) and LD (2.31, 95% CI=1.59–3.38, p<0.001) (Tables 3 and 4)., Similarly, d4T-containing regimens showed significantly higher IRR in both LA (2.46, 95% CI=1.48–4.36, p<0.001) and LD (2.10, 95%CI=1.23–3.80, p=0.003) (Tables 3 and 4). In the case of the Fas −670AA genotype, however, significantly higher IRR was observed only in LA (1.68, 95% CI=1.15–2.41, p=0.006) but not in LD (1.19, 95%CI=0.74–1.83, p=0.433) (Tables 3 and 4). The statistical significance of higher IRR of LA for Fas −670AA did not change after adjustment for gender and starting drug regimens in multivariate Poisson regression analysis (1.72, 95% CI=1.20–2.45, p=0.003, Table 3). The difference was statistically significant even after correction for multiple testing of three SNPs in two groups (LA or LD) (p=0.018). These results suggested that the Fas −670 genotype affects LA without fat accumulation.
Diagnosis of lipoatrophy.
Gender, NRTIs, and Fas−670 were included as covariates.
CD4 cells/μl.
IRR, incidence rate ratios; NRTIs, nucleoside reverse transcriptase inhibitors; d4T, stavudine.
Diagnosis of lipoatrophy.
Gender, NRTIs, and Fas−670 were included as covariates.
CD4 cells/μl.
Lipodystrophy cases included 17 cases diagnosed with both lipoatrophy and lipodystrophy.
IRR, incidence rate ratios; NRTIs, nucleoside reverse transcriptase inhibitors; d4T, stavudine.
Lipid profiles of HIV-1-infected Thai patients
We analyzed the serum lipid profiles of the patients. From the records of 829 patients, serum cholesterol levels were obtained in 610 and 628 cases at the first and last measurement of the observation period, respectively. Serum triglyceride (TG) levels were obtained in 646 and 716 cases at the first and last measurement, respectively. Serum HDL levels were obtained in 393 and 367 cases at the first and last measurement, respectively. Serum low-density lipoprotein (LDL) levels were obtained in 382 and 438 cases at the first and last measurement, respectively. Serum cholesterol levels showed a normal distribution, while serum TG, HDL, and LDL levels did not. Therefore, we used a nonparametric method to evaluate the statistical significance of differences in patient lipid profiles. Slight increases in serum cholesterol levels were observed in the LA (from median 194.5 to 206 mg/dl, p=0.016), LD (from median 194 to 205 mg/dl, p=0.001), and non-LA/LD cases (from median 195 to 200 mg/dl, p=0.033) during ART (Table 5). Similarly, significant increases in TG were observed in the LA (from median 130 to 183.5 mg/dl, p=0.019) and LD (from median 122 to 165.5 mg/dl, p<0.001), but not in non-LA/LD cases (from median 138 to 148.5 mg/dl, p=0.359; Table 5). Accordingly, the LA and LD cases showed significantly higher TG levels than the non-LA/LD cases at the last measurement (p=0.007; Table 5). These results indicated that lipid metabolism in LA and LD cases was more severely damaged by ART than that in non-LA/LD cases. It should also be noted here that the LA and LD cases showed weak but clear trends toward lower TG levels than non-LA/LD cases at the first measurement (p=0.091; Table 5). In fact, when combined, the LA and LD cases showed significantly lower TG levels (median 122 mg/dl) than the non-LA/LD cases (median 138 mg/dl) at the first measurement (p=0.049). On the other hand, HDL and LDL levels remained almost unchanged during the treatment period in the LA and LD, but non-LA/LD cases showed a slight decrease in LDL (from median 137 to 132.5 mg/dl, p=0.026; Table 5).
Lipoatrophy cases.
Lipodystrophy cases.
Patients without lipoatrophy or lipodystrophy.
Kruskal–Wallis test.
The first measurement after patients started antiretroviral therapy.
The last measurement of the observation period. For LA and LD cases, this corresponded to the diagnosis of LA and LD, respectively.
Wilcoxon signed rank test.
LA, lipoatrophy; LD, lipodystrophy; IQR, interquartile range; TG, triglyceride; HDL, high-density lipoprotein; LDL, low-density lipoprotein.
As the LA and LD cases showed similar lipid profiles, we combined the LA and LD cases again and analyzed the effects of patients' genotypes on the lipid profile together. The results indicated that the ApoC3 −455CC genotype was associated with decreased serum cholesterol levels in non-LA/LD cases at the last measurement (p=0.035; Table 6). However, this association was not observed in LA and LD cases (p=0.393; Table 6) or at the first measurement in both LA and LD cases and non-LA/LD cases (p=0.254 and p=0.583, respectively). The ApoC3 −455CC genotype was also found to be associated with elevated serum TG levels in the LA and LD cases at the time of LA or LD diagnosis (p=0.021; Table 6). In the non-LA/LD cases, the ApoC3 −455CC genotype had virtually no effect on serum TG levels. These results suggested that the ApoC3 −455CC genotype promotes dysregulation of lipid metabolism in LA or LD cases, even though this allele had no deleterious effect on the onset of LA or LD (Table 2). On the other hand, the ApoC3 −455CC genotype affected neither HDL nor LDL serum levels (data not shown). The ARβ3 cod64 genotype had no effect on the serum lipid data of any patient group (Table 6). Similarly, Fas −670AA, which was found to be significantly associated with LA, did not show any effects on serum lipid levels (Table 6).
Lipoatrophy and lipodystrophy cases.
Patients without lipoatrophy or lipodystrophy.
Kruskal–Wallis test.
LA, lipoatrophy; LD, lipodystrophy; IQR, interquartile range.
Discussion
In the present study, we evaluated the influence of SNPs involved in apoptosis and lipid metabolism on developing LA or LD in Thai HIV-1 patients. Specifically, we analyzed the SNPs of the Fas A–670G, ApoC3 T–455C, and ARβ3 Tcod64C in patients undergoing ART.
Of the SNPs investigated, our findings indicated that the AA genotype of Fas −670 was significantly associated with the risk of LA but not that of LD. This result is consistent with that of a previous study in Italian HIV patients by Zanone Poma et al. 29 It was suggested that adipocyte apoptosis induced by mitochondrial toxicity of NRTIs was the primary mechanism involved in LA. 15 Consistent with this suggestion, increased Fas gene expression levels were observed in adipose tissue of HIV-1 patients with NRTI-associated LA. 40 The Fas A–670G polymorphism is located in the γ-interferon activation site, to which transcription factors such as signal transducer and activator of transcription (STAT) bind. 41 Therefore, this SNP may affect the level of Fas gene transcription and may induce apoptosis, consequently increasing the risk of LA. The Fas A–670G polymorphism has been reported to be associated with several autoimmune diseases and inflammatory disorders, 41 –43 which may involve a mechanism similar to LA. The frequency of the AA genotype of Fas −670 in 829 HIV-1-infected Thais was 22.1%, which is similar to the percentages in Europeans and Japanese but is markedly higher than that in Africans. As Fas −670AA caused a 70% increase in risk of LA, this genotype is a candidate for future pretherapy genetic screening, although it is still necessary to find other SNPs that increase the risk of LA for accurate screening. It is also necessary to confirm our findings on Fas −670AA in a larger patient cohort.
On the other hand, the frequency distribution of the ARβ3 cod64 and ApoC3 −455 genotypes was not significantly different between LA or LD cases and non-LA/LD cases. The ARβ3 Tcod64C polymorphism failed to show any effect on lipid metabolism among Thai patients in this study. The results for ARβ3 cod64 were consistent with a previous report that showed no association of this polymorphism with lipid metabolism during weight loss in obese subjects. 44 However, in the Italian cohort, this genotype was associated with the risk of developing LA. 29 Several studies have reported that people with this SNP showed an increased risk of developing obesity and glucose intolerance, 32,33 whereas other studies failed to find any correlation between this mutation and body mass index. 45,46
The results for ApoC3 −455 were also not consistent with those of the Italian cohort, which showed an association of the ApoC3 −455T allele with LA. However, it was also reported that the ApoC3 −455C but not the ApoC3 −455T allele was associated with metabolic syndrome. 47 In fact, a recent study in French HIV patients showed that the ApoC3 −455C but not the ApoC3 −455T allele is associated with the severity of ART-induced dyslipidemia and occurrence of LD. 48,49 Therefore, the precise role of the ApoC3 T–455C polymorphism in LA or LD during ART should be carefully reevaluated in a larger cohort of HIV patients.
Nevertheless, our study also showed that the ApoC3 −455CC genotype was associated with elevated serum levels of TG, which was consistent with the findings of a previous report showing that the ApoC3 −455C allele was associated with metabolic syndrome. 47 These data suggest that the mechanism underlying the promotion of LA and LD in HIV-infected individuals receiving ART may be different from that promoting metabolic syndrome. However, the differences in results may be due to ethnic differences in the cohorts. It was reported that genetic polymorphisms in hepatic lipase showed different effects on plasma lipid levels in different ethnic groups. 50 Foulkes et al. found race- and ethnic-specific differences in plasma lipid levels under ART, as well as differences in the influence of ApoC3 T–455C and the other SNP in ApoC3 gene on the development of protease inhibitor-related hypertriglyceridemia. 51 Further studies including reevaluation of these SNPs in larger patient cohorts are therefore necessary to elucidate the precise mechanisms underlying LA and LD in HIV-infected individuals receiving ART.
In conclusion, the Fas −670AA genotype, but not the ApoC3 −455 or ARβ3 cod64 genotypes, affected the incidence of LA in HIV-1-infected Thai patients. None of these alleles affected the incidence of LD. On the other hand, the ApoC3 −455CC, but not the ARβ3 cod64 genotype, affected the serum levels of TG.
Footnotes
Acknowledgments
We thank Ms. Taksin Pimpak for her help in statistical analysis, Professor Massimo Galli at the University of Milan for his support, and all the participants in this study. This work was supported by grants from the Health Science Foundation, the Ministry of Education, Culture, Sports, Science, and Technology, and the Ministry of Health, Labour, and Welfare, Japan.
Author Disclosure Statement
No competing financial interests exist.
