Abstract
Tenofovir disoproxil fumarate (TDF) associates with renal tubular dysfunction (RTD) in some people living with HIV (PLWH). We studied clinical and genetic factors associated with RTD in Thai PLWH receiving TDF. RTD was diagnosed in 13 of 65 (20%) patients. The median (interquartile range) age and CD4 cell counts were 43.8 (40.4-50.9) years and 554 (437-716) cells/mm3, respectively. The median duration of TDF use was 46.9 (31.5-54.1) months. Univariate logistic regression demonstrated body mass index (BMI), concomitant use of protease inhibitor (PI), hyperlipidemia, and homozygous C/C SNP rs1059751 of ABCC4 gene as predisposing factors of RTD. In multivariate model, concomitant use of PI [adjusted odds ratio (aOR) 11.39; 95% confidence interval (CI), 1.59- 81.56; P = 0.015], hyperlipidemia (aOR 8.59; 95% CI, 1.46-50.40; P = 0.017), and BMI (aOR 0.76; 95% CI, 0.59-0.98; P = 0.037) remained associated with RTD in patients receiving TDF. PLWH receiving TDF with the presence of these factors should be closely monitored for RTD.
Introduction
Tenofovir disoproxil fumarate (TDF) is a nucleotide reverse transcriptase inhibitor that is widely used for the treatment of people living with human immunodeficiency virus type 1 (HIV-1). TDF has been recommended as a component of the nucleoside backbone of the preferred HIV treatment regimens in Thailand since 2010. 1 The drug is an orally bioavailable prodrug of tenofovir (TFV). 2 TFV is excreted by glomerular filtration and active tubular secretion. Although TFV is considered safe and tolerable in clinical trials,3,4 multiple cases of proximal renal tubulopathy and renal failure have been reported in practical settings.5-7 A study in the Asia-Pacific reported TDF-nephrotoxicity with an incidence of 1.75 per 100 person-years. 8
Although the subsequent more favorable renal safety profile pro-drug formulation, tenofovir alafenamide (TAF), is available, some concerns on weight and metabolic changes, 9 as well as an economic value may limit the use of TAF. 10
The mechanism of TFV-induced kidney injury is not completely understood. It is suggested that mitochondria are a target of TFV toxicity. Mitochondrial damage in the proximal renal tubular cells was observed in people living with HIV (PLWH) with renal tubular dysfunction (RTD) receiving TDF.11,12 TFV is a weak inhibitor of mitochondrial DNA polymerase-γ, and mitochondrial DNA depletion has been demonstrated in animal models of TFV toxicity. 13
Transporter proteins in the renal proximal tubules play a part in controlling the intracellular concentration of TFV and hence alteration in their expression could associate with TFV-related RTD. TFV enters renal tubule epithelial cells through the basolateral membrane and involves organic anion transporters (OATs), OAT1, and OAT314,15 that are encoded by the genes SLC22A6 and SLC22A8, respectively. OAT4, encoded by the gene SLC22A11, is expressed in the apical membrane of renal proximal tubule cells. A polymorphism in this gene (rs11231809) has been shown to interfere with diuretic torsemide excretion 16 and could affect the renal clearance of other compounds. Another transporter on the luminal membrane; multidrug resistance protein (MRP) 4, encoded by the adenosine triphosphate–binding cassette (ABC) gene ABCC417,18 contributes to TFV excretion into the urine. Associations between polymorphisms at the ABCC2 gene and TFV-associated RTD have been demonstrated,19,20 although the involvement of MRP2 with TFV-associate RTD remains unclear. MRP7, encoded by ABCC10, has also been shown to play a role in TFV transport in the proximal tubule and the SNPs were associated with tubular toxicity. 21 P-glycoprotein, encoded by the ABCB1 gene, is a membrane protein expressed on the renal proximal tubular cells, intestinal cells, and hepatocytes.22-24 P-glycoprotein transports TDF, therefore the alteration of P-glycoprotein expression in enterocytes could affect TFV exposure.
TFV-associated nephrotoxicity is multifactorial. Genetic polymorphisms of renal tubular cell transporter proteins have been shown to affect the risk of developing TFV-related RTD.19,20,25 Clinical factors, including increasing age, low body weight, pre-existing decreased kidney function, coexisting diabetes mellitus, and concomitant use of protease inhibitor (PI) or nephrotoxic drugs are known to be risk factors for nephrotoxicity.26-28
We aimed to study the association between RTD and clinical factors as well as 57 single nucleotide polymorphisms (SNPs) in the ABCC2, ABCC4, ABCC10, SCL22A6, SCL22A11, and ABCB1 genes in Thai PLWH receiving TDF. Characterizing risk factors could guide frequent monitoring of renal tubular function in PLWH with particular risk factors whom TDF is needed or switching to other proper antiretroviral regimens.
Methods
Study Population
A cross-sectional study was conducted to investigate the association between clinical factors and SNPs in genes encoding renal tubular transporter and TFV-associated RTD in Thai PLWH. Thai adults (age >15 years) living with HIV receiving TDF-containing antiretroviral regimens between September 2012 and January 2013 were invited to participate in this study. PLWH invited were the subset 149 PLWH from the previous study on the association between SNPs and nevirapine-induced rash, who were treated with stavudine, lamivudine, and nevirapine regimen, 29 and stavudine was switched to TDF. The exclusion criteria included PLWH with fasting plasma glucose >200 mg/dL or diagnosed with diabetes mellitus and those denied to participate in the study. Demographic data were extracted from the medical records.
Measurement and Definition of Renal Proximal Tubular Dysfunction
Blood and spot urine sample were collected together with body weight and height measurement. RTD was diagnosed by the presence of at least three of the following criteria: β2-microglobulinuria (urine β2-microblobulin >1000 μg/g Cr), 20 nondiabetic glucosuria (urine glucose >25 mg/dL), 30 increased fractional excretion of uric acid (FEUricAcid): (uric acid(urine) /uric acid(serum))/(creatinine(urine)/creatinine(serum)); abnormal >0.1, increased fractional excretion of phosphate FEPhos) : (phosphate(urine) / phosphate(serum)) / (creatinine(urine) / creatinine(serum)); abnormal >0.2, (>0.1 with serum phosphate <2.5 mg/dL), and renal tubular acidosis (serum bicarbonate <21 mmol/L and urinary pH >5.5). 31
Selection of Variations in Genes Related to TFD Pharmacodynamic
SNPs in genes encoding tubular transporters were chosen for the study based on their functional significance from previously published literatures.20,25 Fifty-seven SNPs in the ABCC2, ABCC4, ABCC10, SCL22A6, SCL22A11, and ABCB1 genes data retrieved from the previous genome-wide association study (GWAS) samples 29 were collected. GenABEL package 32 in the R statistical program was used to estimate the significance of candidate SNPs. Bonferroni adjustment was used for the correction of multiple testing. SNPs with complete homozygosity and minimal allele frequency less than 0.05 were excluded for further analysis.
Statistical Analysis
Median (interquartile range, IQR) and frequencies (%) were used to describe characteristics of PLWH. Chi-square test or Fisher's exact test and Mann-Whitney U test were used to compare categorical variables and continuous variables between the two groups, respectively. Associations between clinical factors, alleles, and RTD were tested by univariate and multivariate logistic regression analyses. Variables that presented P <0.1, were considered in a multivariate logistic regression model with forward and backward stepwise logistic regression. Odds ratio (OR) and adjusted odds ratio (aOR) and their 95% confidence interval (CI) were calculated. A P-value of <0.05 was considered statistically significant. The receiver operating characteristic (ROC) curve and the Youden index were used to determine optimal cut-off value of body mass index (BMI) to classify low versus high risk for RTD. The statistical analyses were conducted using Stata statistical software version 17.0 (StataCorp, College Station, TX, USA).
Ethical Approval and Informed Consent
The protocol of this study was reviewed and approved by the Committee on Human Rights Related to Research Involving Human Subjects, Faculty of Medicine Ramathibodi Hospital, Mahidol University (MURA2012/217). All participants provided written informed consent prior to enrollment in the study.
Results
A total of 65 PLWH who were able to provide blood and urine samples and fulfilled the study criteria were enrolled in the study. The baseline characteristics and laboratory data for PLWH with and without RTD are shown in Table 1. RTD was diagnosed in 13 of the 65 (20%) PLWH. Thirty-three (50.5%) of the PLWH were female. Median (IQR) age and CD4 cell counts were 43.8 (40.4-50.9) years and 554 (437-716) cells/mm3, respectively. Three PLWH had HIV viral load above the lower limits of detection of the assays. Of the total PLWH in the study, 54 (83.1%), 10 (15.4%), and 1 (1.5%) PLWH receiving non-nucleoside reverse transcriptor, PI, and integrase inhibitor-based regimens, respectively. The median duration of TDF using was 46.9 (31.5-54.1) months. Although serum creatinine was comparable in the two groups, creatinine clearance calculation based on both Cockroft-Gault and Modification of Diet in Renal Disease (MDRD) was lower in PLWH with RTD. PLWH with RTD had lower BMI compared to those without RTD. There were significant numbers of individuals with hyperlipidemia and receiving PI-containing regimens in PLWH with RTD. Renal tubular markers and the contribution of each parameter to RTD in this study are shown in Table 1. Univariate analysis of clinical factors showed a significant association between RTD and BMI, concomitant use of PI-containing regimens, and diagnosed hyperlipidemia (Table 2).
Characteristics of Patients with and Without Renal Tubular Dysfunction.
Univariate Analysis of Clinical Factors of Renal Tubular DysFunction in PLWH Receiving TDF.
SNPs with significant levels after conservative and their distribution of genotypes and alleles are shown in Table 3. All polymorphisms were in Hardy-Weinberg equilibrium in all samples. The single SNP analysis showed a greater percentage of individuals with RTD among the SNP 4976T→C, rs1059751 of ABCC4. The effect of this SNP was more evident in PLWH with C/C homozygote than heterozygote following the additive genetic model. OR of having RTD in PLWH carrying C/C homozygote was 13.13 times compared to those who harbored T/T homozygote (Table 4). No significant association was found between other SNPs and RTD.
Genotype, Allele Frequencies and Univariate Analysis of SNPs in PLWH Receiving TDF.
*Not applicable because of complete discrimination between the two groups of patients with and without renal tubular dysfunction, no odds ratio could be estimated for the SNP rs2125739.
Multivariate Analysis of Significant Clinical and Genetic Risks Factor Associated with RTD in PLWH Receiving TDF.
Each variable was adjusted for the SNP rs1059751.
Multivariate logistic regression model combining both clinical and genetic factors was performed, only clinical factors remained independently associated with RTD in PLWH receiving TDF; concomitant use of PI (aOR 11.39; 95% CI, 1.59- 81.56; P = 0.015), hyperlipidemia (aOR 8.59; 95% CI, 1.46-50.40; P = 0.017), and BMI (aOR 0.76; 95% CI, 0.59-0.98; P = 0.037). Youden index showed predictive ability of BMI to predict low versus high risk for RTD at a cut-off of 18.6 kg/m2 (a sensitivity, specificity, area under the ROC curve of 58.3%, 80.9%, and 0.70, respectively).
Discussion
We investigated the risk factors of Thai PLWH receiving TDF. RTD in this group of PLWH is not uncommon. We demonstrated low BMI, concomitant use of PI, and hyperlipidemia as associating factors of RTD in Thai PLWH receiving TDF. When clinical and genetic factors were combined in the multivariate regression model, only clinical factors remained statistical significant.
TFV-associated nephrotoxicity is multifactorial. Clinical and genetic polymorphisms of renal tubular cell transporter proteins have been shown to influence the risk of developing TFV-related RTD.3,19,20,26,33 Previous studies also showed low BMI as one of the risk factors for TFV-associated renal dysfunction.33,34 Pharmacokinetic studies demonstrated that small body weight is associated with reduced plasma TFV clearance and thus high plasma TFV concentrations, which could result in RTD.35,36 Concomitant use of PI has been shown to associate with risk factors for tubular dysfunction and/or nephrotoxicity.26,37 The mechanism could be explained by possible inhibition of TDF excretion from the proximal tubule and increasing the intracellular concentration. A previous study demonstrated lower renal clearance of TFV in PLWH on lopinavir-ritonavir combination than in those not taking a protease inhibitor. 38 Our findings were consistent with earlier studies. We also revealed an association between hyperlipidemia and proximal tubular dysfunction. Dyslipidemia is known to accelerate the progression of chronic kidney disease leading to tubule-interstitial injury. In cell culture, low-density lipoprotein (LDL) was demonstrated to induce stress in human proximal tubular cells by multiple mechanisms. 39
Although genetic polymorphism in this study was not associated with RTD by the multivariate model in the study, OR of having C/C genotype at SNP rs1059751 was 13 times compared to T/T genotype. The mechanism by which polymorphism in the ABCC4 gene influences the risk of RTD could be explained by its encoding protein, MRP4, function. MRP4 is the transporter for TFV excretion at the luminal membrane of the kidney proximal tubule cells. The efflux at the luminal membrane is a rate-limiting step for many antiviral drugs and occasionally results in intracellular accumulation. The SNP 4976T→C, rs105975 locates in the 3′UTR region of the gene. Polymorphism in this region may affect translation efficiency or mRNA stability of the protein and hence could affect tubular excretion of TFV. A previous study in Thai population also showed the association of ABCC4 4976T→C variation with β2 microgloburinuria. 40 Other studies demonstrated the association of ABCC4 3463A→G and ABCC4 4131T→G polymorphisms in association with higher intracellular 36 and plasma 41 concentrations of TFV, respectively.
Previous studies have demonstrated genetic polymorphisms in association with RTD. Izzedine et al reported the role of 1249 G-->A SNP and CATC haplotype ( − 24C, 1249A, 3563T, 3972C) of ABCC2 in RTD. 19 This study defined renal toxicity based on metabolic acidosis, urine potassium loss, hypophosphatemia, low uric acid levels, and aminoaciduria within one month after TDF initiation. Rodriguez-Novoa et al demonstrated genotype CC at ABCC2 position −24 as one of the risk factors of RTD. 25 Nishijima et al also demonstrated genotype CC at position − 24 and additional genotype AA at position 1249 of the ABCC2 gene associated with RTD in the Japanese population. 20 Rodriguez-Novoa et al studied the PLWH with a median duration of exposure to TDF of 34 months and the presence of at least 2 criteria of the proximal tubular abnormalities. Nishijima et al collected data in the Japanese population with a median 35-month follow-up time in PLWH with the presence of at least three criteria of the proximal tubular dysfunction. We did not find a significant association of previously identified SNPs with RTD in the present study in Thai PLWH. The previous study on factors associated with TDF discontinuation also showed no significant association of genetic variants in the multivariable analyses. 42
The results of these studies, however, are not exactly comparable owing to inconsistency in the criteria for RTD diagnosis. Our cross-sectional study collected PLWH with a median TDF treatment duration of 47 months. The criteria for the diagnosis of RTD in our study were different from those of previous studies. The discrepancy in diagnosis criteria suggested the need for uniformly established criteria for the diagnosis of TDF-associated RTD.
The strength of our study is the genetic factors that were retrieved from GWAS, in which common polymorphisms in genes of the candidate transporter proteins were investigated. This study, however, has some limitations. First, the study was a small cross-sectional study among PLWH receiving TDF in a subset of the previously studied cohort. Second, PLWH with more significant renal toxicity at earlier time points could have been missed, and baseline renal function was not collected. These factors might lead to selection bias. Third, other transporters that were not considered in this study could be involved in the transportation of the drug. Prospective studies with a larger sample size and full genetic analysis are needed.
In conclusion, concomitant use of PI-containing antiretroviral regimen, hyperlipidemia, and low BMI are significantly associated with RTD in Thai PLWH receiving TDF. Closed monitoring of RTD should be warranted in PLWH who have these factors.
Footnotes
Acknowledgements
We thank all participants in the study for their kind cooperation.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethical Approval and Informed Consent
The protocol of this study was reviewed and approved by the Committee on Human Rights Related to Research Involving Human Subjects, Faculty of Medicine Ramathibodi Hospital, Mahidol University (MURA2012/217). All patients provided written informed consent prior to enrollment in the study.
Funding
This study was supported by research grants from Pharmacogenomics Projects, the collaboration between Ramathibodi Hospital, Mahidol University and Thailand Center of Excellence of Life Sciences (TCELS).
