Research Article

The Impact of Single Nucleotide Polymorphisms on the Pharmacokinetics of Tacrolimus in Kidney Allograft Recipients of Northern-West, Iran

Introduction

Calcineurin inhibitors (CNI) are administered in 96% of allograft recipients and provided significant benefits in short-term graft survival.¹ However, it is shown that the long-term survival of the graft is limited by the nephrotoxicity of these drugs.² In a 10-year histopathologic follow-up of the transplanted kidney grafts, the use of CNIs was associated with decreased early subclinical rejection while progressive arteriolar hyalinosis, glomerulosclerosis, and tubulointerstitial damage were observed.³ Although tacrolimus, in contrast to cyclosporine, causes less pathological changes in the transplanted allografts within the first year, its pathological insults were similar to those of cyclosporine’s after a long time.⁴ Since tacrolimus has therapeutic benefit with detrimental effects on graft survival, understanding its pharmacokinetics and interactions is important.⁵

Tacrolimus is metabolized by cytochrome P450 (CYP) 3A4 and 3A5 enzymes and variations in their coding genes were shown to be associated with altered drug clearance.⁶ Patients with high clearance of tacrolimus at first 90 days after transplantation are at an increased risk of acute graft rejection.⁷ CYP 3A4*22 T carriers require fewer doses of tacrolimus.⁸ The CYP 3A5*1 allele, which encodes the functional enzyme, causes a two-fold decrease in dose-normalized plasma levels of tacrolimus.⁹ This allele is more frequent among the patients with African-Americans ethnicity and associated with a 50% increase in total dosage requirements of tacrolimus to reach the first therapeutic trough level.¹⁰ In a study, among 337 renal allograft recipients, the presence of the CYP 3A4*22 allele and expression of CYP 3A5 were predictors of 20% and 160% clearance fluctuation, respectively.¹¹ Thus, a personalized medicinal approach is necessary to determine the tacrolimus dosing to maintain therapeutic range in renal allograft recipients.

The A6986G (rs776746) single nucleotide polymorphism (SNP) in the CYP 3A5 gene has been a focus of researchers during recent years. This SNP is located on the 3^rd intron of the gene and impacts the CYP 3A5 expression.¹² However, there are conflicting findings among the studies. Additionally, two other SNPs including the C3435T of multi-drug resistance-1 (MDR-1) and C25385T pregnane X receptor (PXR, also known as nuclear receptor 112) genes were also investigated for their impact on the pharmacokinetics of tacrolimus. The obtained results are widely different across populations with different genetic backgrounds.

In the clinics, kidney recipients differently respond to a given pharmacological therapy, some of whom cannot reach the designated concentration(s) with recommended starting doses. The over-dosing of the drug increases the risks of nephrotoxicity, infection, and other severe drug-specific side effects, while its under-dosing may lead to the acute rejection. Hence, the CNIs management is a challenging issue in transplant patients. Considering that tacrolimus is widely used in transplant immunosuppression and there are inter-individual differences in respond to the drug, understanding the patients’ pharmacogenetics and personalized administration of the drug could be helpful in different populations. Therefore, in this study, the influence of the CYP3A5, MDR-1, and PXR genes SNPs on the pharmacokinetics of tacrolimus and the ratio of concentration per dose of tacrolimus (C/D ratio) were evaluated in a group of kidney recipients in Northwest of Iran.

Methods

Patients

In this cross-sectional study, renal allograft recipients receiving tacrolimus were included. During the study period, renal transplant recipients with an age range of 20-60 years under triple therapy with tacrolimus (Prograf^® manufactured by Astellas Company, in 2 divided doses), mycophenolate mofetil (MMF), and prednisolone were included. Patients with ongoing acute allograft rejection, BK virus nephropathy, hepatitis, and cytomegalovirus (CMV) infections were excluded from the study.

Laboratory examinations

A baseline laboratory examination including plasma creatinine, blood urea nitrogen, lipid profile, total bilirubin, and 25-hydroxy vitamin D levels were measured. The measurement of plasma levels of tacrolimus trough level was performed exactly 12 hours after the last dose of the drug and immediately before taking the next dose. Tacrolimus metabolism rate (C/D ratio) was calculated by dividing the blood concentration of the drug (C, ng/mL) to the daily dose of Tacrolimus (D, mg) at 2, 4 and 8 weeks after transplantation.

The genomic DNA was extracted from blood with the magnetic nanoparticle-based method using ZiAViZ^® DNA extraction kit (Itan). Amplification refractory mutation systempolymerase chain reaction (ARMS-PCR) method was performed using primers (Table 1) designed based on previous studies.^13,14 The commercial master mix kit (Sinaclon-Iran, #CatNo. MM2062) was used for the PCR procedure (Mastercycler Eppendorf, Germany).

Statistical analysis

The Hardy-Weinberg equilibrium was checked for each SNP. IBM SPSS Statistics version 24.0 was used for statistical analysis. Descriptive statistics demonstrates the frequency, percentage, mean, and standard deviation (SD). A logistic regression model with a confidence interval (CI) of 95% was applied for the evaluation of the correlations among the variables. ANOVA or Kruskal-Wallis tests were applied to compare the differences between tacrolimus C/D ratio at different time points (2, 4, and 8 weeks). The chi-square test was applied to check the correlation between the viral infections and graft rejection vs. SNP genotypes. A P value of < 0.05 was considered significant.

Results and Discussion

In total, 65 renal recipients were included in the study with a male/female ratio of 37:28. The average age was 38 ± 14.17 years old. The demographic data of the patients are presented in Table 2. It is reported that in the CYP3A5 nonexpressers, genotyping of NR1I2/ MDR1/ CYP3A4 polymorphisms can be useful for controlling tacrolimus dosing.¹⁵ In the present study, the variant allele frequencies of CYP3A5*3, MDR-1 C3435T, and PXR C25385T were 95.38, 20.77, and 26.92%, respectively.

Allele frequency and genotype distribution of the CYP3A5*3

Most of the kidney recipients (90.8%, n = 59) were homozygotes for the CYP3A5*3 allele (CYP3A5 *3/*3, non-expressers CYP3A5), while 9.2% of the participants (n = 6) were heterozygous for the CYP3A5*3 (CYP3A5 *1/*3 genotype, functional CYP3A5 expressers). Across different populations, the frequency of CYP3A5*3 SNP differs widely. In our kidney recipient population, the frequency of the CYP3A5*3 allele was 95%. The result of a meta‐analysis has stated that the allele frequency of CYP3A5*3 is higher in European (94%), admixed American (80%), East Asian (71%), and South Asian (67%) subjects, while it is low (33%) in African subjects.¹⁶ Another meta-analysis (2015) concluded that the CYP3A5 A6986G SNP could affect the pharmacokinetics of tacrolimus.¹⁷ In the carriers of the wild-type allele (CYP3A5*1, expresser genotype, fast metabolizers), the expression of CYP3A5 leads to an elevated metabolism and clearance of tacrolimus and a lower C/D ratio compared to the other genotypes (CYP3A5*1/*3 and CYP3A5*3/*3 carriers). South African kidney transplant recipients in comparison with global transplant populations have a high rate of CYP3A5 expression that significantly affects the pharmacokinetics of tacrolimus.¹⁸ However, some studies did not find a statistically significant correlation between the A6986G SNP and higher clearance of tacrolimus¹⁹ being in agreement with our study.

Allele frequency and genotype distribution of the MDR-1 C3435T

The frequency of the C and T alleles in MDR-1 C3435T gene was 79.23 and 20.77%, respectively. The genotypes distribution and allele frequency of the studied polymorphism are presented in Table 3. The MDR-1 affects the pharmacodynamics of tacrolimus by active transporting of a wide variety of drugs to the outsides of the cells.²⁰ Similar to the CYP3A5, wild-type carriers of at least one MDR-1 C allele are expected to express more P-gp. The MDR1 C3435T SNP is a synonymous mutation in which the ATC codon changes to ATT (Ile1145Ile), causing alteration in P-gp activity.^21-24 A higher activity of the MDR-1 limits the enteral absorption of tacrolimus. Similar to the results of Loh et al,²⁵ 67.7% of our recipients were homozygote C allele carriers in contrast to other reports.²⁶ In addition, study by Tada et al²⁷ revealed that the MDR1 C3435T polymorphism did not affect tacrolimus pharmacokinetics of renal transplanted patients. Among the Asian population, SNPs in the MDR-1 and CYP3A5 genes were shown to influence the plasma levels of tacrolimus but not cyclosporine. Moreover, it has been proposed that diltiazem, a non-dihydropyridine calcium channel blocker, may help achieve the optimum levels of tacrolimus by blocking the MDR-1 and CYP 3A4 proteins.²⁵

Allele frequency and genotype distribution of the PXR C25385T SNP

The frequency of the C allele in the PXR genes was 73.08%. The genotypes distribution and allele frequency of this SNP are presented in Table 3. The PXR gene, encoding a nuclear receptor, is involved in the regulation of CYP3A enzymes, MDR-1, and other several proteins and therefore, plays a crucial role in the pharmacokinetics of tacrolimus.²⁸ However, there are contradictions about its effect on the pharmacokinetics of tacrolimus.²⁹ This study revealed that the frequency of the PXR C25385T allele was 26.92%.

The correlation between clinical variables and gene polymorphisms

The median (min-max) C/D ratio (ng/mL/mg) for the kidney recipients at 2, 4, and 8 weeks were 1.2 (0.13-5), 1.68 (0.32-10.67), and 2 (0.33-13), respectively. There was a significant difference in C/D ratios at 2 and 8 weeks (P = 0.015). In addition, the C/D ratios of patients with different genotypes in the studied genes are given in Table 4. The mean of C/D ratios in the non-expresser group compared to the expresser group were 1.56 ± 1.1 vs. 0.83 ± 0.16 ng/mL/mg at 2 weeks, 2.01 ± 1.0 vs. 1.75 ± 0.91 ng/mL/mg at 4 weeks, and 2 vs. 1.2 at 8 weeks after renal transplantation, respectively. There was a significant difference when C/D ratios of homozygote CYP3A5 *3/*3 carriers were compared between 2 and 8 weeks (P = 0.012).

Six patients experienced acute graft rejection, none of the SNPs seemed to have a significant impact on rejection rates (P> 0.05). In addition, 10 patients experienced viral infections (CMV, BK virus, or both), however, none of the SNP genotypes had a significant predisposing influence on them (P> 0.05).

Pharmacogenetics and tacrolimus metabolism

Based on the studied genotypes and C/D ratios after 2, 4, and 8 weeks, renal transplant recipients were categorized into three groups (Table 5) based on the rate of tacrolimus metabolism. The groups were included C/D ratio < 1.05 (ng/mL/mg)/ (mg/kg/d) as fast metabolizers, 1.05 < C/D ratio < 1.55 (ng/mL/mg)/ (mg/kg/d) as intermediate metabolizers, and C/D ratio > 1.55 (ng/mL/mg)/ (mg/kg/d) as slow metabolizers, Thölking et al.³⁰

It is reported that the expresser genotype (with a C/D ratio < 1.05) is correlated with a higher risk of chronic nephrotoxicity and acute rejection and a lower estimated glomerular filtration rate (eGFR) compared to slow and intermediate tacrolimus metabolizers.^17,31,32 In the present study, about half of the kidney recipients with CYP3A5*3/*3 genotype were slow-metabolizers, confirming the impact of other genes in the pharmacokinetics of tacrolimus.

It is shown that the 3435 C/C genotype of MDR-1 gene is associated with 40% lower C/D ratios.²⁹ Since drug concentrations accomplished with a standard dose will be lower for carriers of the wild-type C allele, C/D ratios is expected to be lower. Dissimilar to this expectation, our patients with the C allele had higher C/D ratios at week 2; 1.57 ± 1.03 for homozygotes and 1.25 (min-max 0.29-5) for heterozygotes versus 0.8 ± 0.4 T allele. The SNPs in CYP3As and MDR-1 genes are attracting attention in clinical practice.³³ Our results showed that the frequency of the T allele was 20.77% among renal transplanted, 3 of the carriers were fast metabolizers with a 0.91 ± 0.11 C/D ratio and 2 patients were slow metabolizers with a 1.7 ± 0.16 C/D ratio. However, this polymorphism did no significant effect on the metabolism of tacrolimus.

Kurzawsk et al demonstrated that minor T allele carriers (CT or TT genotype) of the PXR C25385T require significantly low tacrolimus dose administration particularly in the first 6 months after renal transplantation compared to CC genotypes.³⁴ In the present study, 2 of fast metabolizers with PXR SNP TT genotype had a 0.86 ± 0.19 C/D ratio and 6 patients were slow metabolizers with a 2.74 ± 1.42 C/D ratio. There were no significant differences between the studied polymorphism and fast/ intermediate/slow metabolizers (Table 5).

The nephrotoxicity of CNIs has led to the assumption that perhaps lowering the dosage of CNIs or utilizing a less toxic agent (e.g. belatacept, sirolimus) would be a safer option.^2,35 On the contrary, the triple regimen of tacrolimus, MMF, and a corticosteroid is ubiquitously used at present.³⁶ It is shown that lower plasma levels of tacrolimus with a higher dosage of MMF are associated with better outcomes.³⁷ Scholten et al developed a flexible method for accurate dosing of tacrolimus regardless of the genetic variables which is an ongoing field of research.³⁸ There are somewhat controversial findings from other studies. We believe that the differences among the methodologies, tacrolimus release (immediate or extended), and perhaps ethnic differences may lead to such results.

Conclusion

In this study, the frequency of the CYP3A5*3 allele was 95%. There was a significant difference in the C/D ratios at 2 and 8 weeks among homozygote CYP3A5 *3/*3 carriers, suggested that choosing the initial dosage according to the CYP3A5 genotype may result in a better outcome. No significant correlations were found between the MDR-1, and PXR gene SNPs and the tacrolimus C/D ratios. These result will be a step toward personalized medicine and may prolong the graft survival in renal allograft recipients. Pre-transplant genetic study of recipient could prevents low drug level rejection or high level nephrotoxicity to occur.

Acknowledgments

The authors greatly appreciate the Kidney Research Center, Tabriz University of Medical Sciences, Tabriz, Iran for scientific and financial support (Grant number: 59763). This paper was extracted from the residential thesis of Elaheh Jabbari-Hagh, MD. Authors also thanks Transplantation Ward of Imam Reza Hospital of Tabriz University of Medical Science, Tabriz, Iran.

Competing Interests

Authors declare no conflicts of interest regarding this manuscript.

Ethical Approval

The Helsinki Declaration of ethics in medical research was honored in this study. The informed written consent was obtained from all patients. The study was approved by the ethics committee of Tabriz University of Medical Sciences, Tabriz, Iran and is registered with the following ethical code: IR.TBZMED.REC.1397.128.

References

1. The impact of calcineurin inhibitors on graft survival. Transplant Rev (Orlando) 2013; 27(3):93-5. doi: 10.1016/j.trre.2013.04.003 [Crossref]
2. Calcineurin inhibitors in kidney transplantation: friend or foe?. CurrOpin Nephrol Hypertens 2011; 20(6):610-5. doi: 10.1097/MNH.0b013e32834b4343 [Crossref]
3. The natural history of chronic allograft nephropathy. N Engl J Med 2003; 349(24):2326-33. doi: 10.1056/NEJMoa020009 [Crossref]
4. Calcineurin inhibitor nephrotoxicity through the lens of longitudinal histology: comparison of cyclosporine and tacrolimus eras. Transplantation 2016; 100(8):1723-31. doi: 10.1097/tp.0000000000001243 [Crossref]
5. Clinical pharmacokinetics of tacrolimus. Clin Pharmacokinet 1995; 29(6):404-30. doi: 10.2165/00003088-199529060-00003 [Crossref]
6. Tacrolimus population pharmacokinetics and multiple CYP3A5 genotypes in black and white renal transplant recipients. J Clin Pharmacol 2018; 58(9):1184-95. doi: 10.1002/jcph.1118 [Crossref]
7. High tacrolimus clearance is a risk factor for acute rejection in the early phase after renal transplantation. Transplantation 2017; 101(8):e273-e9. doi: 10.1097/tp.0000000000001796 [Crossref]
8. The role of pharmacogenetics in the disposition of and response to tacrolimus in solid organ transplantation. Clin Pharmacokinet 2014; 53(2):123-39. doi: 10.1007/s40262-013-0120-3 [Crossref]
9. Tacrolimus pharmacogenetics: the CYP3A5*1 allele predicts low dose-normalized tacrolimus blood concentrations in whites and South Asians. Transplantation 2005; 79(4):499-502. doi: 10.1097/01.tp.0000151766.73249.12 [Crossref]
10. Prevalence of CYP3A5 genomic variances and their impact on tacrolimus dosing requirements among kidney transplant recipients in eastern North Carolina. Pharmacotherapy 2017; 37(9):1081-8. doi: 10.1002/phar.1970 [Crossref]
11. A population pharmacokinetic model to predict the individual starting dose of tacrolimus in adult renal transplant recipients. Br J Clin Pharmacol 2019; 85(3):601-15. doi: 10.1111/bcp.13838 [Crossref]
12. Capability of utilizing CYP3A5 polymorphisms to predict therapeutic dosage of tacrolimus at early stage post-renal transplantation. Int J Mol Sci 2015; 16(1):1840-54. doi: 10.3390/ijms16011840 [Crossref]
13. Influence of CYP3A5 and ABCB1 polymorphism on tacrolimus drug dosing in South Indian renal allograft recipients. Indian J Nephrol 2019; 29(4):261-6. doi: 10.4103/ijn.IJN_97_18 [Crossref]
14. Genetic polymorphisms influencing xenobiotic metabolism and transport in patients with primary biliary cirrhosis. Hepatology 2005; 41(1):55-63. doi: 10.1002/hep.20516 [Crossref]
15. Interactive effects of CYP3A4, CYP3A5, MDR1 and NR1I2 polymorphisms on tracrolimus trough concentrations in early postrenal transplant recipients. Pharmacogenomics 2015; 16(12):1355-65. doi: 10.2217/pgs.15.78 [Crossref]
16. Worldwide distribution of cytochrome P450 alleles: a meta-analysis of population-scale sequencing projects. Clin PharmacolTher 2017; 102(4):688-700. doi: 10.1002/cpt.690 [Crossref]
17. Effect of CYP3A5*3 on kidney transplant recipients treated with tacrolimus: a systematic review and meta-analysis of observational studies. Pharmacogenomics J 2015; 15(1):38-48. doi: 10.1038/tpj.2014.38 [Crossref]
18. CYP3A5 polymorphisms and their effects on tacrolimus exposure in an ethnically diverse South African renal transplant population. S Afr Med J 2020; 110(2):159-66. doi: 10.7196/SAMJ.2020.v110i2.13969 [Crossref]
19. The impact of CYP3A5 and MDR1 polymorphisms on tacrolimus dosage requirements and trough concentrations in pediatric renal transplant recipients. Saudi J Kidney Dis Transpl 2013; 24(6):1125-36. doi: 10.4103/1319-2442.121268 [Crossref]
20. Therapeutic implications of the MDR-1 gene. J Vet PharmacolTher 2004; 27(5):257-64. doi: 10.1111/j.1365-2885.2004.00607.x [Crossref]
21. Functional polymorphisms of the human multidrug-resistance gene: multiple sequence variations and correlation of one allele with P-glycoprotein expression and activity in vivo. Proc Natl Acad Sci U S A 2000; 97(7):3473-8. doi: 10.1073/pnas.97.7.3473 [Crossref]
22. MDR1 gene polymorphisms and disposition of the P-glycoprotein substrate fexofenadine. Br J Clin Pharmacol 2002; 53(5):526-34. doi: 10.1046/j.1365-2125.2002.01591.x [Crossref]
23. C3435T polymorphism in the MDR1 gene affects the enterocyte expression level of CYP3A4 rather than Pgp in recipients of living-donor liver transplantation. Pharmacogenetics 2002; 12(6):451-7. doi: 10.1097/00008571-200208000-00005 [Crossref]
24. P-glycoprotein: from genomics to mechanism. Oncogene 2003; 22(47):7468-85. doi: 10.1038/sj.onc.1206948 [Crossref]
25. Significant impact of gene polymorphisms on tacrolimus but not cyclosporine dosing in Asian renal transplant recipients. Transplant Proc 2008; 40(5):1690-5. doi: 10.1016/j.transproceed.2008.04.010 [Crossref]
26. MDR1 haplotype frequencies in Japanese and Caucasian, and in Japanese patients with colorectal cancer and esophageal cancer. Drug MetabPharmacokinet 2006; 21(2):126-32. doi: 10.2133/dmpk.21.126 [Crossref]
27. Impact of CYP3A5 and MDR1(ABCB1) C3435T polymorphisms on the pharmacokinetics of tacrolimus in renal transplant recipients. Transplant Proc 2005; 37(4):1730-2. doi: 10.1016/j.transproceed.2005.02.073 [Crossref]
28. NR1I2 polymorphisms are related to tacrolimus dose-adjusted exposure and BK viremia in adult kidney transplantation. Transplantation 2012; 94(10):1025-32. doi: 10.1097/TP.0b013e31826c3985 [Crossref]
29. Genetic polymorphisms and individualized tacrolimus dosing. Transplant Proc 2010; 42(8):3031-3. doi: 10.1016/j.transproceed.2010.08.001 [Crossref]
30. The tacrolimus metabolism rate influences renal function after kidney transplantation. PLoS One 2014; 9(10):e111128. doi: 10.1371/journal.pone.0111128 [Crossref]
31. Tacrolimus dose requirements and CYP3A5 genotype and the development of calcineurin inhibitor-associated nephrotoxicity in renal allograft recipients. Ther Drug Monit 2010; 32(4):394-404. doi: 10.1097/FTD.0b013e3181e06818 [Crossref]
32. Influence of tacrolimus metabolism rate on renal function after solid organ transplantation. World J Transplant 2017; 7(1):26-33. doi: 10.5500/wjt.v7.i1.26 [Crossref]
33. Design of allele specific PCR for rapid detection of CYP3A5 (A6986G) and Mdr-1 (C3435T) polymorphisms. Indian J Clin Biochem 2011; 26(1):18-21. doi: 10.1007/s12291-010-0085-z [Crossref]
34. Analysis of common polymorphisms within NR1I2 and NR1I3 genes and tacrolimus dose-adjusted concentration in stable kidney transplant recipients. Pharmacogenet Genomics 2017; 27(10):372-7. doi: 10.1097/fpc.0000000000000301 [Crossref]
35. Leas BF, Uhl S, Sawinski DL, Trofe-Clark J, Tuteja S, Kaczmarek JL, et al. Calcineurin inhibitors for renal transplant. Rockville, MD: Agency for Healthcare Research and Quality (US); 2016.
36. Immunosuppression: evolution in practice and trends, 1994-2004. Am J Transplant 2006; 6(5 Pt 2):1111-31. doi: 10.1111/j.1600-6143.2006.01270.x [Crossref]
37. Relationship of tacrolimus exposure and mycophenolate mofetil dose with renal function after renal transplantation. Transplantation 2011; 92(1):82-7. doi: 10.1097/TP.0b013e31821fad06 [Crossref]
38. AUC-guided dosing of tacrolimus prevents progressive systemic overexposure in renal transplant recipients. Kidney Int 2005; 67(6):2440-7. doi: 10.1111/j.1523-1755.2005.00352.x [Crossref]