Next Article in Journal
Automated Quality Control Solution for Radiographic Imaging of Lung Diseases
Previous Article in Journal
Perioperative Use of IgM-Enriched Immunoglobulins in Liver Transplantation Recipients at High Risk for Infections: A Preliminary Study
Previous Article in Special Issue
Comparison of Percutaneous Core Needle Biopsy Results in Patients Who Previously Underwent Open and Robot-Assisted Kidney Transplantation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Vitamin D Metabolism Genes and Their Genomic Background in Shaping Cyclosporine A Dosage Parameters after Kidney Transplantation

1
Clinic of Maxillofacial Surgery, Pomeranian Medical University in Szczecin, 70-111 Szczecin, Poland
2
Department of Immunological Diagnostics, Pomeranian Medical University in Szczecin, 70-111 Szczecin, Poland
3
Institute of Physical Culture Sciences, University of Szczecin, 70-453 Szczecin, Poland
4
Department of Personalized Medicine and Cell Therapy, Regional Blood Center, 60-354 Poznan, Poland
5
Department of General Pathology, Pomeranian Medical University in Szczecin, 70-111 Szczecin, Poland
6
Department of Nephrology, Transplantology and Internal Medicine, Pomeranian Medical University in Szczecin, 70-111 Szczecin, Poland
7
Department of General Pharmacology and Pharmacoeconomics, Pomeranian Medical University in Szczecin, 71-210 Szczecin, Poland
8
Department of General Surgery and Transplantation, Pomeranian Medical University in Szczecin, 70-111 Szczecin, Poland
9
Department of Interdisciplinary Dentistry, Pomeranian Medical University, 70-204 Szczecin, Poland
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2024, 13(16), 4966; https://doi.org/10.3390/jcm13164966
Submission received: 13 June 2024 / Revised: 12 August 2024 / Accepted: 19 August 2024 / Published: 22 August 2024
(This article belongs to the Special Issue Kidney Transplantation: Current Challenges and Future Perspectives)

Abstract

:
Background: Kidney transplantation is followed by immunosuppressive therapy involving calcineurin inhibitors (CNIs) such as cyclosporin A. However, long-term high CNIs doses can lead to vitamin D deficiency, and genetic variations influencing vitamin D levels can indirectly impact the necessary CNIs dosage. This study investigates the impact of genetic variations of vitamin D binding protein (DBP) rs2282679 and CYP2R1 hydroxylase rs10741657 polymorphisms on the cyclosporin A dosage in kidney transplant recipients. Additional polymorphisims of genes that are predicted to influence the pharmacogenetic profile were included. Methods: Gene polymorphisms in 177 kidney transplant recipients were analyzed using data mining techniques, including the Random Forest algorithm and Classification and Regression Trees (C&RT). The relationship between the concentration/dose (C/D) ratio of cyclosporin A and genetic profiles was assessed to determine the predictive value of DBP rs2282679 and CYP2R1 rs10741657 polymorphisms. Results: Polymorphic variants of the DBP (rs2282679) demonstrated a strong predictive value for the cyclosporin A C/D ratio in post-kidney transplantation patients. By contrast, the CYP2R1 polymorphism (rs10741657) did not show predictive significance. Additionally, the immune response genes rs231775 CTLA4 and rs1800896 IL10 were identified as predictors of cyclosporin A response, though these did not result in statistically significant differences. Conclusions: DBP rs2282679 polymorphisms can significantly predict the cyclosporin A C/D ratio, potentially enhancing the accuracy of CNI dosing. This can help identify patient groups at risk of vitamin D deficiency, ultimately improving the management of kidney transplant recipients. Understanding these genetic influences allows for more personalized and effective treatment strategies, contributing to better long-term outcomes for patients.

1. Introduction

Chronic Kidney Disease (CKD) affects 100 million people in Europe and is projected to become the fifth leading cause of death worldwide by 2040. The actual number of patients may be even higher, as many remain unaware of the disease due to the lack of early diagnosis. Furthermore, CKD is one of the most costly diseases for healthcare systems, with an estimated annual expenditure of EUR 140 billion in Europe.
As the disease progresses, renal replacement therapy is frequently required, with kidney transplantation representing the optimal treatment, accompanied by the essential inclusion of immunosuppressive therapy. Calcineurin inhibitors (CNIs) were introduced into clinical practice in the early 1980s and quickly became the cornerstone of immunosuppressive therapy following organ transplantation. The first drug of this class to be introduced was cyclosporin A (CsA). At present, the first-line immunosuppressive drug is tacrolimus, which forms the basis of immunosuppressive regimens. However, tacrolimus and other immunosuppressants may induce vitamin D deficiency [1,2].
The latter presents a serious problem in transplant recipients [3,4]. The main causes of vitamin D deficiency in this group include the impaired uptake of calcidiol by the transplanted kidney, the indirect effects of immunosuppressive drugs, hyperparathyroidism, and disruptions to the calcium and phosphate balance [5]. This issue arises because vitamin D has a pleiotropic effect, ranging from the skeletal system to the regulation of over 900 genes involved in innate and adaptive immunity [6,7,8,9,10,11,12]. In kidney transplant recipients, the complications stemming from vitamin D deficiency pose a significant concern regarding alterations in the skeletal system [4]. At the same time, vitamin D deficiency is associated with dysregulated immune responses, the development of insulin resistance, vascular dysfunction and damage, and cardiomyopathy. Hesketh et al. demonstrated that vitamin D helps to modulate immunosuppression [13]. Other studies have linked vitamin D deficiency to an increased risk of kidney transplant failure [14,15].
In humans, the major source of vitamin D is the cutaneous synthesis of cholecalciferol (vitamin D) following exposure to ultraviolet B (UVB) radiation during sun exposure, which is influenced by body temperature. The next step is 25-hydroxylation, which occurs in the liver primarily by P450 2R1 25-hydroxylase (CYP2R1), resulting in the synthesis of calcidiol. Polymorphisms in the CYP2R1 gene are associated with variable 25-hydroxylase activity, which affects serum 25 (OH)D levels [16]. The serum concentration of calcidiol strongly correlates with the single nucleotide polymorphism rs2282679 A>C in the GC gene. The GC gene encodes for vitamin D binding protein (DBP)—a protein crucial for calcidiol and calcitriol transportation. In the renal tubules, CYP27B1 1α-hydroxylase mediates a secondary hydroxylation that produces the active form of vitamin D—calcitriol. Calcitriol can bind intracellularly to the nuclear vitamin D receptor (VDR) and thus exert its target effect [17]. CYP3A4 polymorphisms influence the pharmacokinetic mechanisms involved in the metabolism of both tacrolimus and cyclosporine A [18,19,20].
Genetic variability due to multiple polymorphisms may alter the bioavailability and modulate the effects of vitamin D metabolites. For instance, vitamin D supplementation can affect serum levels differently depending on the genetic variant of the CYP2R1 gene in a healthy individual [16,21]. Cyclosporin A (CsA) and tacrolimus inhibit hepatic CYP3A4 (25-hydroxylase activity), resulting in reduced vitamin D levels [2].
In exploring the interplay between CNIs and vitamin D activity at different stages, we aimed to investigate the association of two major vitamin D polymorphisms with parameters of cyclosporine A dosage. The main objective of our study was to examine the influence of the polymorphisms rs10741657 (CYP2R1) and rs2282679 (GC) on the cyclosporine A C/D ratio, while considering other selected genomic factors potentially associated with calcineurin inhibitor pharmacodynamics. The ontology of other genes involved in our analysis pertained to the immune response and hepatic drug metabolism. Additionally, we sought to develop a predictive algorithm to optimize cyclosporine A dosing. Previous studies have demonstrated success with similar approaches in various post-transplantation patient populations [22,23,24,25,26].
Given our genetic database comprising a large, homogeneous patient cohort and considering the growing, yet unclear, body of information on the influence of vitamin D on CNI therapy, we found it valuable to employ data mining techniques to explore this field.

2. Materials and Methods

The patients were recruited at the Department of Nephrology and Kidney Transplantation of the Independent Public Provincial Hospital in Szczecin and the Department of General and Transplant Surgery of the Pomeranian Medical University in Szczecin. The research was approved by the Bioethics Committee of the Medical University of Karol Marcinkowski in Poznań—resolution No. 560/17. All patients signed an informed consent form.
Only patients with stable graft function who were undergoing routine follow-up at the outpatient clinic were included in the study. Patients with unstable graft function who required hospitalization were excluded from the study group.

2.1. Selection and Characteristics of the Study Group

The study cohort consisted of 177 patients of both genders, ranging in age from 18 to 80, who had undergone kidney transplantation and were receiving cyclosporine (CsA), a calcineurin inhibitor, as part of their immunosuppressive therapy. To assess the immunosuppressant concentration, blood samples were collected from the bend of the elbow vein before the morning drug dose administration for each patient in a fasting state. A brief overview of the characteristics of the study group is presented in Table 1 and Table 2.
The genetic study of the polymorphisms was carried out on the abovementioned biological material, which was stored at a temperature of −20 °C.

2.2. Isolation of Samples

Genomic DNA isolation from blood samples was performed using the QIAamp Mini Blood kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The purity and concentration of each DNA sample were assessed using spectrophotometry (EPOCH Spectrophotometer, Agilent Technologies, Santa Clara, CA, USA). Targeted polymorphisms were analyzed by real-time PCR using LightCycler® 480 (Roche Diagnostics, Rotkreuz, Switzerland). For genotyping, hybridization probes were used as fluorescent dyes labeled with SimpleProbe (Roche, Basel, Switzerland). The results were analyzed based on the melting curve using the LightCycler® 480 Basic Software (16 February 2023). The LightSNiP rs2282679 and rs10741657 collection of polymorphisms contains the appropriate concentrations of specific primers and probes for the amplified DNA fragments, prepared according to the manufacturer’s instructions. The concentration of CsA in the whole blood was determined using the microparticle chemiluminescent immunoassay (CMIA) (ARCHITECT i2000SR analyzer, Abbott laboratories, Green Oaks, IL, USA). The C/D ratio (concentration/dose ratio) was calculated for each patient by dividing the values of the detected blood concentration to the applied dose.
The study of polymorphisms was carried out in the Laboratory of Experimental Pharmacogenetics of the Department of Clinical Pharmacy and Biopharmacy and the Department of Stem Cells and Regenerative Medicine of the Institute of Natural Fibers and Medicinal Plants of the Medical University of Poznań. We involved the data on 16 polymorphisms identified in 15 genes. A complete list of the genes and polymorphisms analyzed is presented in Table 3.
The genetic database used for this analysis was developed gradually over several years. Polymorphisms were selected based on their possible pharmacogenetic or pathophysiological impact on the metabolism of CNI drugs or related processes involving the immune response.

2.3. Statistical Analysis

The analysis of the relationship between the C/D ratio and the genetic profile was performed using data mining methodology. The Random Forest algorithm was used, followed by Classification and Regression Trees (C&RT) in Statistica software (ver. 13.1, StatSoft Polska, Cracow, Poland). The Random Forest analysis is frequently used in pharmacology research and was successfully employed by others [22,24,27,28,29,30,31]. The analysis was supplemented with ANOVA tests together with Tukey’s post-hoc test for various n, in the case of variables satisfying the condition of normal distribution. The Shapiro–Wilk test was used to assess the normality of the distribution.

3. Results

In the first stage, we performed an analysis using the Random Forest algorithm to determine the importance of individual predictors. Random Forest analysis for CsA revealed the 49A>G RS231775 polymorphism in the CTLA4 gene as the highest rank predictor. The rs2282679GC polymorphisms of the vitamin D binding protein gene and the rs10741657 of CYP2R1 25-hydroxylase gene polymorphisms revealed similar, lower ranks.
Table 3. Importance of individual predictors for CsA. Random Forest analysis of DNA samples isolated from the peripheral blood of patients post-kidney transplantation treated with cyclosporin A as an immunosuppressant revealed that the rs2282679GC polymorphism had a high predictive rank, while the rs10741657 polymorphism of the CYP2R1 25-hydroxylase gene showed a similar but lower predictor importance.
Table 3. Importance of individual predictors for CsA. Random Forest analysis of DNA samples isolated from the peripheral blood of patients post-kidney transplantation treated with cyclosporin A as an immunosuppressant revealed that the rs2282679GC polymorphism had a high predictive rank, while the rs10741657 polymorphism of the CYP2R1 25-hydroxylase gene showed a similar but lower predictor importance.
PolymorphismRank 1
rs231775 CTLA4100
rs602662 FUT294
rs2282679 GC90
rs1800896 IL1089
rs1800795 IL688
rs8175347 UGT1A1*2886
rs4654748 NBPF383
rs2032582 MDR80
rs1045642 MDR79
rs2069762 IL276
rs10741657 CYP2R172
rs1800468 TGFB63
rs1800629 TNF-alfa53
rs1227200453
rs33972313 SLC23A141
rs2740574 CYP3A422
1 Predictor importance for CsA with C/D ratio as a dependent variable.
In the next stage, the Classification and Regression Trees (C&RT) algorithm was used to create binary regression models for cyclosporine A C/D ratio values. The hierarchy of individual ranks obtained from the Random Forest algorithm was considered. Predictors with importance values of less than 0.7 for CsA were excluded. A 10-fold cross-validation was applied. The following validation parameters were used: the cost of resubstitution, the cost of cross-check, and the standard error.
Figure 1 shows the regression model for the C/D CsA values. The following attributes were included in the analysis: rs231775 CTLA4, rs602662 FUT2, rs2282679 GC, rs1800896 IL10, rs1800795 IL6, rs8175347 UGT1A1*28, rs4654748 NBPF3, rs2032582 MDR, rs1045642 MDR2, and rs1045642 MDR. Detailed tree characteristics and validation parameters are presented in Figure 2. The model revealed the rs2282679 GC polymorphism as an overarching predictor and a series of leaves described by the mean C/D value and coexisting polymorphisms. The rs10741657 CYP2R1 polymorphism was not identified as a predictor in the model. The mean C/D values of individual leaves and nodes are shown in boxes. Considering the binary nature of the tree, four final leaves were distinguished qualitatively, characterized by at least one specific polymorphism. The genetic characteristics of all four abovementioned leaves is listed below.
  • Leaf ID 6: rs2282679 CA, IL 10-1082A>G RS1800896 GA or GG, MDR2677G>T RS2032582 GT
  • Leaf ID 5: rs2282679 CA, IL10-1082A>G RS1800896 AA
  • Leaf ID 10: rs2282679 CC or AA, rs602662 FUT2 GG or AG, IL10-1082A>G RS1800896 AA
  • Leaf ID 9: rs2282679 CC or CA, rs602662 FUT2 AA
IL101082A>GRS1800896 AA appeared as a defining single polymorphism for two of the leaves (ID 5 and ID 10). However, MDR2677G>Trs2032582 (leaf ID 6) and rs602662FUT2 (leaf ID 9) were decisive for statistical significance in cyclosporine C/D ratios (p = 0.006, Figure 2 and Table 4).

4. Discussion

In the presented study, using the Random Forest method, it was shown that the studied polymorphisms, which have a strong influence on vitamin D metabolism, are associated with a high predictive rank for the C/D ratio of cyclosporine. The significance of the rs2282679 GC polymorphism was then confirmed in the C&RT model.
There is currently no reliable data in the literature regarding the variability of effective drug dosage over time, which may be influenced by factors such as epigenetic changes. Considering this, the concentration/dose ratio, which is regarded as the most independent measure of therapy efficacy, was used for the analysis in the presented study.
The GWAS study on a large cohort of subjects confirmed the significance of the polymorphisms rs2282679, rs12785878, and rs10741657 in relation to serum concentrations of 25 (OH) D [32,33,34]. It was suggested that certain variants might serve as serum deficiency predictors [35]. Homozygosity for the C allele in rs2282679 has been associated with lower serum concentrations of vitamin D compared to AC heterozygosity and AA homozygosity [32,35]. Consequently, this has been found to increase the risk of vitamin D deficiency and a suboptimal response to vitamin D supplementation [36], in both adults and children, irrespective of the season [37]. Perna et al. showed that carriers of at least one copy of the less frequent allele exhibit a lower efficiency of skin vitamin D synthesis [38]. Dong et al. reported that rs2282679 alleles negatively affect serum vitamin D levels in pregnant women. Using decision tree methodology, they created a model to assess the risk of deficiency, and based on the results, the clinical procedure was modified [39]. This study effectively integrates the results of genetic determinations with data mining methods, which ultimately become tools of clinical significance.
Similarly to the rs2282679 GC, the rs10741657 CYP2R1 polymorphism is associated with vitamin D serum concentration. It is believed to regulate gene expression and, therefore, may influence the level of 25-hydroxylase activity and expression. Based on genome-wide association studies by Ahn et al. [40], Anderson et al. [41], O’Brien et al. [42], and Jiang et al. [43], more than 20 SNPs of the CYP2R1 gene were found to affect vitamin D concentration in the European population. Carriers of two G alleles in rs10741657 are characterized by significantly lower vitamin D concentrations compared to AG heterozygotes and AA homozygotes. The GG genotype shows decreased CYP2R1 enzyme activity, which translates into a marked reduction in 25 (OH) D concentrations compared to the AA genotype. In several populations, carriers of the G allele exhibited the lowest serum concentrations of 25 (OH) D, demonstrating that the tested SNP correlates with 25 (OH) D levels. By contrast, the AA and AG genotypes of CYP2R1 rs10741657 are associated with a greater increase in 25 (OH) D after UVB exposure compared to GG homozygotes [43]. The relationship between rs10741657 and the increased risk of vitamin D deficiency was significant in the dominant model (GG + AG/AA), but not in the recessive model (GG/AG + AA) [44].
Vitamin D deficiency due to genetic variation can be effectively addressed through supplementation, and genotyping can aid in identifying at-risk groups [45]. Many studies investigating the association between genetic variation and vitamin D concentration focus on polymorphisms that affect the functioning of DBP and vitamin D hydroxylases [46]. Our study also addresses this issue. The analysis using the Random Forest algorithm confirmed that both studied polymorphisms (rs2282679 GC and rs10741657 CYP2R1) had a high predictive rank with respect to the C/D ratio in patients after kidney transplantation. However, the C&RT algorithm analysis did not confirm the predictive significance of the rs10741657 CYP2R1 polymorphism in the tested model; therefore, the discussion will focus on the role of the rs2282679 GC polymorphism.
Kidney transplantation is the gold standard in the treatment of end-stage renal disease (ESRD), providing patients with statistically longer survival compared to other renal-replacement therapies. Any improvements that enhance and prolong organ function and improve the quality and length of life for recipients are of particular interest. One of the methods proposed in this work is the genotyping of polymorphisms involved in the vitamin D biosynthesis pathway as a significant factor modulating recipients’ homeostasis. Our results, supported by studies conducted in various centers, demonstrate that the GC polymorphism (rs2282679) in the gene for vitamin D binding protein impacts therapy with calcineurin inhibitors.
There is a relationship between rs2282679 DBP polymorphism and carcinogenesis. As mentioned, the C allele, in comparison to AA homozygotes, increases the risk of a reduced DBP concentration and vitamin D deficiency. It has also been shown that DBP levels are associated with a shorter interval between relapses in colorectal cancer and an increased risk of several other cancers [47,48,49]. DBP is a precursor molecule for a potent macrophage activating factor, and AA homozygotes exhibit higher serum DBP concentrations, which may exert a more substantial immunomodulatory effect than in the AC/CC genotypes. These data could be relevant for transplantation patients with an increased cancer risk [50].
Moreover, the rs2282679 polymorphism and lower levels of vitamin D are associated with greater disease severity and a higher risk of death in patients with SARS-CoV-2 infection. These results demonstrate that DBP plays an active role in inflammation, including a chemotactic effect on neutrophils and the activation of macrophages. In addition, DBP has a protective effect on thromboembolic complications during COVID-19 [51].
There is also a significant correlation between DBP concentration and bone mineral density (BMD). While low BMD is not a disease marker per se, it reflects overall bone health. A healthy skeletal and muscular system determines the maintenance of mobility—an important factor in preventing chronic diseases. Low serum DBP is associated with disturbed calcium homeostasis, decreased bone mineral density, and an increased incidence of osteoporotic fractures, particularly in postmenopausal women [52]. Postmenopausal recipients require special attention, as reduced concentrations of vitamin D and DBP are often observed in this patient group [53,54]. Genotyping could contribute to the development of individualized and more effective treatment plans, thereby reducing the risk of complications related to calcium and phosphate metabolism disorders [55]. DBP may also influence BMD through another mechanism: by serving as a reservoir for vitamin D metabolites, suggesting that DBP can regulate their bioavailability and the half-life of vitamin D [52]. The association of the two major GC variants with serum DBP and BMD has been confirmed [52]. Therefore, identifying a sensitive and predictive biomarker for detecting osteoporosis would benefit patients.
The work of Ponticelli and Passerini described numerous gastrointestinal complications in patients after kidney transplantation, which may lead to the loss of the transplanted organ [56]. In the context of kidney transplantation, immunosuppressive therapy can disrupt the intestinal immune environment, increasing susceptibility to damage. This disruption compromises the epithelial barrier of the intestinal mucosa, allowing prolonged exposure to antigens. Calcitriol enhances the intestinal barrier by stimulating the expression of tight junction proteins between enterocytes [57,58]. The loss of the intestinal barrier function is also thought to contribute to the development of GvHD. The vitamin D-dependent release of cathelicidin and the protection of epithelial barriers (including those in other parts of the digestive and genitourinary tracts) may affect the transplanted organ’s functionality during infection. Furthermore, vitamin D inhibits inflammatory mediators released by neutrophils, such as TNF-α, IL-6, and IL-8 [59].
Vitamin D and VDR may regulate the synthesis of endothelial nitric oxide (NO) [60]. Nitric oxide has a vasodilating effect, which protects blood vessels. Vitamin D supplementation may reduce endothelial dysfunction, as confirmed by long-term studies [61,62,63,64,65]. Adequate levels of vitamin D have a protective effect under oxidative stress conditions, which can occur during transplantation [66]. In addition, cyclosporin A and tacrolimus inhibit NO synthesis, leading to a reduced endothelial response to NO and a secondary lack of renal vasodilation [67,68]. There is also an increased risk of post-transplant diabetes mellitus (PTDM) associated with low vitamin D levels and impaired glucose metabolism [69]. It has been shown that pancreatic β cells express the enzyme 1α-hydroxylase and vitamin D receptors, which may influence insulin secretion [70].
The immunomodulatory properties of vitamin D affect the doses of immunosuppressants used by organ recipients [71]. Following heart transplantation, patients receiving low-dose calcitriol required significantly lower cumulative doses of cyclosporin A without affecting the incidence of rejection episodes, infection, or death [23]. Vitamin D enhances the expression of the CYP3A4 gene. Lindh et al. showed that this effect varies with the season, as vitamin D concentration is strongly related to insolation. Furthermore, as reported, individual changes in vitamin D levels may influence the metabolism of drugs including cyclosporin A, tacrolimus, and sirolimus. Notably, the authors also used the C/D ratio for their analysis [71]. Vitamin D affects the expression of genes dependent on it through its receptor [72]. The effect of vitamin D on CYP3A4 transcription is related to its interaction with a specific VDR receptor [71,73]. Therefore, it is possible that the CYP3A4 cytochrome might be the key to understanding the interaction between vitamin D and CNIs [74]. Both vitamin D and calcineurin inhibitors are CYP3A4 substrates. Additionally, vitamin D has a cytochrome-inducing effect [75,76,77]. In our study, however, the polymorphism rs2740574 CYP3A4 had a very low predictive rank for the cyclosporin A C/D ratio.
Most immunosuppressants have a narrow therapeutic range. The occurrence of most side effects is dose-dependent and correlates with the drug concentration in the blood serum. The main goal in developing immunomodulatory strategies is to establish drug combinations that produce synergistic effects while minimizing toxicity. There is a synergistic relationship between vitamin D function and immunosuppressive therapy [78,79,80,81,82,83]. Notably, Gim et al. analyzed the effect of the variability of rs776746 (CYP3A5) and rs1137115 (CYP2A6) on tacrolimus dosing using data mining methods. They emphasized that combining decision tree algorithms with Random Forest algorithms is an effective tool for pharmacogenetic analysis [84]. The decision tree algorithm has been shown to outperform other data mining methods [26].
In our study, rs2282679GC demonstrated a superior predictive capability regarding cyclosporine A dosage parameters compared to other polymorphisms. Interestingly, immune response-related genes were also identified as predictors—rs231775 CTLA4 in Random Forest and rs1800896 IL10 in both Random Forest and the C&RT model. However, these predictors were not decisive for statistically significant differences. This may be attributed to the limited number of observations, which represents a limitation of this study.
Given the numerous side effects of CNI therapy, efforts have been made to minimize or exclude their use from therapeutic regimens. In the absence of a better solution, research should focus on personalizing CNI therapy [85]. Determining single nucleotide polymorphisms (SNPs) is a crucial tool in developing a personalized approach to transplantation. Nevertheless, genome-wide multicenter association studies are needed to advance personalized transplant medicine [86].
The initial concept behind pharmacogenetics aimed to create a genetic ID for individual patients, allowing physicians to make optimal therapeutic decisions. Assuming that several genes directly affect drug metabolism in an exclusive way, this could be a milestone in therapy. Subsequent years of research have cast a shadow of uncertainty over this idea. Drug metabolism is influenced by many variable factors, not only genetic ones. Although we managed to demonstrate a strong relationship between the studied polymorphisms and CsA dosing, the study has certain limitations. As mentioned earlier, drug metabolism can be affected by various factors, including environmental conditions, diet, interactions with concomitant drugs, etc. Additionally, factors affecting the metabolism and bioavailability of vitamin D, such as the quality of supplementation and the accuracy of vitamin D concentration measurements during periodic check-ups, also play a role [87]. Creating a model that incorporates all of these variables exceeds the capabilities of our database. Moreover, accessing kidney transplant recipients is challenging due to their relatively small numbers compared to other patient groups, such as those with hypertension. It should be concluded that further, extensive, highly standardized multi-center studies are necessary, including much larger groups of patients. Data mining methods are an ideal analysis tool for this purpose.

5. Conclusions

The analysis using the random forest algorithm revealed that among other polymorphisms, both studied polymorphisms (rs2282679 GC and rs10741657 CYP2R1) had high predictive ranks for the C/D ratio-dependent variable. Among these, rs2282679 GC represented an overreaching predictive value in the C&RT model, whereas rs10741657 CYP2R1 was not found to be predictive. Significant differences in cyclosporine C/D ratios were associated with combinations of rs2282679 GC and two other polymorphisms. Immune response-related polymorphisms were also identified as predictive by both algorithms; however, these did not result in statistically significant differences in cyclosporine C/D ratios. Precise pharmacogenomic profiling holds promise for enhancing the efficiency and safety of immunosuppressive therapy following kidney transplantation. Given the pleiotropic effects of vitamin D, the genetic factors influencing its metabolism are potentially important in this context.

Author Contributions

Conceptualization, K.K., M.K. and J.S.; methodology, M.K., A.B. and B.C.; software, B.W. and K.T.; validation, J.S. and M.K.; formal analysis, K.K. and M.K.; investigation, K.K.; resources, M.K. and J.S.; data curation, M.K., K.K. and B.W.; writing—original draft preparation, K.K.; writing—review and editing, I.S., K.K., S.D. and G.T.; visualization, K.T. and B.W.; supervision, M.K. and J.S.; project administration, B.M. and B.C.; funding acquisition, B.M. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

Genetic tests were carried out on the material collected as part of the scientific grant “The importance of pharmacogenetic factors in the individualization of therapy with the use of immunosuppressive drugs in kidney transplants” funded by National Science Center of Poland, grant number 2011/03/B/NZ7/06550.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by Bioethics Committee of the Medical University of Karol Marcinkowski in Poznań—approval No. 83/17 (5 January 2017) and No. 574/18 (16 May 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Filipov, J.J.; Zlatkov, B.K.; Dimitrov, E.P.; Svinarov, D. Relationship between Vitamin D Status and Immunosuppressive Therapy in Kidney Transplant Recipients. Biotechnol. Biotechnol. Equip. 2015, 29, 331–335. [Google Scholar] [CrossRef] [PubMed]
  2. Eyal, O.; Aharon, M.; Safadi, R.; Elhalel, M.D. Serum vitamin D levels in kidney transplant recipients: The importance of an immunosuppression regimen and sun exposure. Isr. Med. Assoc. J. 2013, 15, 628–633. [Google Scholar] [PubMed]
  3. Sadlier, D.M.; Magee, C.C. Prevalence of 25(OH) Vitamin D (Calcidiol) Deficiency at Time of Renal Transplantation: A Prospective Study. Clin. Transplant. 2007, 21, 683–688. [Google Scholar] [CrossRef]
  4. Querings, K.; Girndt, M.; Geisel, J.; Georg, T.; Tilgen, W.; Reichrath, J. 25-Hydroxyvitamin D Deficiency in Renal Transplant Recipients. J. Clin. Endocrinol. Metab. 2006, 91, 526–529. [Google Scholar] [CrossRef] [PubMed]
  5. Yuan, D.; Fang, Z.; Sun, F.; Chang, J.; Teng, J.; Lin, S.; Liu, X. Effect of Vitamin D and Tacrolimus Combination Therapy on IgA Nephropathy. Med. Sci. Monit. 2017, 23, 3170–3177. [Google Scholar] [CrossRef]
  6. Navas-Nazario, A.; Li, F.Y.; Shabanova, V.; Weiss, P.; Cole, D.E.C.; Carpenter, T.O.; Bazzy-Asaad, A. Effect of Vitamin D–Binding Protein Genotype on the Development of Asthma in Children. Ann. Allergy Asthma Immunol. 2014, 112, 519–524. [Google Scholar] [CrossRef]
  7. Freitas, A.T.; Calhau, C.; Antunes, G.; Araújo, B.; Bandeira, M.; Barreira, S.; Bazenga, F.; Braz, S.; Caldeira, D.; Santos, S.C.R.; et al. Vitamin D-Related Polymorphisms and Vitamin D Levels as Risk Biomarkers of COVID-19 Disease Severity. Sci. Rep. 2021, 11, 20837. [Google Scholar] [CrossRef]
  8. Flamann, C.; Peter, K.; Kreutz, M.; Bruns, H. Regulation of the Immune Balance During Allogeneic Hematopoietic Stem Cell Transplantation by Vitamin D. Front. Immunol. 2019, 10, 2586. [Google Scholar] [CrossRef]
  9. Fabri, M.; Stenger, S.; Shin, D.-M.; Yuk, J.-M.; Liu, P.T.; Realegeno, S.; Lee, H.-M.; Krutzik, S.R.; Schenk, M.; Sieling, P.A.; et al. Vitamin D Is Required for IFN-Gamma-Mediated Antimicrobial Activity of Human Macrophages. Sci. Transl. Med. 2011, 3, 104ra102. [Google Scholar] [CrossRef]
  10. Wang, T.-T.; Nestel, F.P.; Bourdeau, V.; Nagai, Y.; Wang, Q.; Liao, J.; Tavera-Mendoza, L.; Lin, R.; Hanrahan, J.W.; Mader, S.; et al. Cutting Edge: 1,25-Dihydroxyvitamin D3 Is a Direct Inducer of Antimicrobial Peptide Gene Expression. J. Immunol. 2004, 173, 2909–2912. [Google Scholar] [CrossRef] [PubMed]
  11. Gombart, A.F.; Borregaard, N.; Koeffler, H.P. Human Cathelicidin Antimicrobial Peptide (CAMP) Gene Is a Direct Target of the Vitamin D Receptor and Is Strongly up-Regulated in Myeloid Cells by 1,25-Dihydroxyvitamin D3. FASEB J. 2005, 19, 1067–1077. [Google Scholar] [CrossRef]
  12. Martineau, A.R.; Wilkinson, K.A.; Newton, S.M.; Floto, R.A.; Norman, A.W.; Skolimowska, K.; Davidson, R.N.; Sørensen, O.E.; Kampmann, B.; Griffiths, C.J.; et al. IFN-Gamma- and TNF-Independent Vitamin D-Inducible Human Suppression of Mycobacteria: The Role of Cathelicidin LL-37. J. Immunol. 2007, 178, 7190–7198. [Google Scholar] [CrossRef] [PubMed]
  13. Hesketh, C.C.; Knoll, G.A.; Molnar, A.O.; Tsampalieros, A.; Zimmerman, D.L. Vitamin D and Kidney Transplant Outcomes: A Protocol for a Systematic Review and Meta-Analysis. Syst. Rev. 2014, 3, 64. [Google Scholar] [CrossRef]
  14. Zeng, S.; Yang, Y.; Li, S.; Hocher, C.F.; Chu, C.; Wang, Z.; Zheng, Z.; Krämer, B.K.; Hocher, B. 25(OH)D-but not 1,25(OH)2D-Is an independent risk factor predicting graft loss in stable kidney transplant recipients. Front. Med. 2023, 10, 1141646. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  15. Thorsen, I.S.; Bleskestad, I.H.; Åsberg, A.; Hartmann, A.; Skadberg, Ø.; Brede, C.; Ueland, T.; Pasch, A.; Reisaeter, A.V.; Gøransson, L.G. Vitamin D as a risk factor for patient survival after kidney transplantation: A prospective observational cohort study. Clin. Transplant. 2019, 33, e13517. [Google Scholar] [CrossRef] [PubMed]
  16. Thacher, T.D.; Fischer, P.R.; Singh, R.J.; Roizen, J.; Levine, M.A. CYP2R1 Mutations Impair Generation of 25-Hydroxyvitamin D and Cause an Atypical Form of Vitamin D Deficiency. J. Clin. Endocrinol. Metab. 2015, 100, E1005–E1013. [Google Scholar] [CrossRef] [PubMed]
  17. Yu, O.B.; Arnold, L.A. Calcitroic Acid-A Review. ACS Chem. Biol. 2016, 11, 2665–2672. [Google Scholar] [CrossRef] [PubMed]
  18. Hernandez, S.; Aquilante, C.; Deininger, K.; Lindenfeld, J.; Schlendorf, K.; Van Driest, S.; Liu, M. Composite CYP3A (CYP3A4 and CYP3A5) phenotypes and influences on tacrolimus dose adjusted concentration in adult heart transplant recipients. Res. Sq. 2023; preprint. https://doi.org/10.21203/rs.3.rs-2921796/v1; Update in: Pharmacogenomics J. 2024, 24, 4. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  19. Wang, C.E.; Lu, K.P.; Chang, Z.; Guo, M.L.; Qiao, H.L. Association of CYP3A4*1B genotype with Cyclosporin A pharmacokinetics in renal transplant recipients: A meta-analysis. Gene 2018, 664, 44–49. [Google Scholar] [CrossRef] [PubMed]
  20. Żochowska, D.; Wyzgał, J.; Pączek, L. Impact of CYP3A4*1B and CYP3A5*3 polymorphisms on the pharmacokinetics of cyclosporine and sirolimus in renal transplant recipients. Ann. Transplant. 2012, 17, 36–44. [Google Scholar] [CrossRef] [PubMed]
  21. Casella, A.; Long, C.; Zhou, J.; Lai, M.; O’Lear, L.; Caplan, I.; Levine, M.A.; Roizen, J.D. Differential Frequency of CYP2R1 Variants Across Populations Reveals Pathway Selection for Vitamin D Homeostasis. J. Clin. Endocrinol. Metab. 2020, 105, 1302–1315. [Google Scholar] [CrossRef]
  22. Greco, R.; Papalia, T.; Lofaro, D.; Maestripieri, S.; Mancuso, D.; Bonofiglio, R. Decisional Trees in Renal Transplant Follow-Up. Transplant. Proc. 2010, 42, 1134–1136. [Google Scholar] [CrossRef] [PubMed]
  23. Briffa, N.K.; Keogh, A.M.; Sambrook, P.N.; Eisman, J.A. Reduction of Immunosuppressant Therapy Requirement in Heart Transplantation by Calcitriol. Transplantation 2003, 75, 2133–2134. [Google Scholar] [CrossRef] [PubMed]
  24. Burlacu, A.; Iftene, A.; Jugrin, D.; Popa, I.V.; Lupu, P.M.; Vlad, C.; Covic, A. Using Artificial Intelligence Resources in Dialysis and Kidney Transplant Patients: A Literature Review. BioMed Res. Int. 2020, 2020, 9867872. [Google Scholar] [CrossRef] [PubMed]
  25. Yoo, K.D.; Noh, J.; Lee, H.; Kim, D.K.; Lim, C.S.; Kim, Y.H.; Lee, J.P.; Kim, G.; Kim, Y.S. A Machine Learning Approach Using Survival Statistics to Predict Graft Survival in Kidney Transplant Recipients: A Multicenter Cohort Study. Sci. Rep. 2017, 7, 8904. [Google Scholar] [CrossRef]
  26. Tang, J.; Liu, R.; Zhang, Y.-L.; Liu, M.-Z.; Hu, Y.-F.; Shao, M.-J.; Zhu, L.-J.; Xin, H.-W.; Feng, G.-W.; Shang, W.-J.; et al. Corrigendum: Application of Machine-Learning Models to Predict Tacrolimus Stable Dose in Renal Transplant Recipients. Sci. Rep. 2018, 8, 46936. [Google Scholar] [CrossRef]
  27. Chandak, P.; Tatonetti, N.P. Using Machine Learning to Identify Adverse Drug Effects Posing Increased Risk to Women. Patterns 2020, 1, 100108. [Google Scholar] [CrossRef]
  28. Ma, Z.; Wang, P.; Gao, Z.; Wang, R.; Khalighi, K. Ensemble of Machine Learning Algorithms Using the Stacked Generalization Approach to Estimate the Warfarin Dose. PLoS ONE 2018, 13, e0205872. [Google Scholar] [CrossRef]
  29. Li, X.; Liu, R.; Luo, Z.-Y.; Yan, H.; Huang, W.-H.; Yin, J.-Y.; Mao, X.-Y.; Chen, X.-P.; Liu, Z.-Q.; Zhou, H.-H.; et al. Comparison of the Predictive Abilities of Pharmacogenetics-Based Warfarin Dosing Algorithms Using Seven Mathematical Models in Chinese Patients. Pharmacogenomics 2015, 16, 583–590. [Google Scholar] [CrossRef] [PubMed]
  30. Cilluffo, G.; Fasola, S.; Ferrante, G.; Malizia, V.; Montalbano, L.; La Grutta, S. Machine Learning: An Overview and Applications in Pharmacogenetics. Genes 2021, 12, 1511. [Google Scholar] [CrossRef]
  31. Cai, N.; Zhang, X.; Zheng, C.; Zhu, L.; Zhu, M.; Cheng, Z.; Luo, X. A Novel Random Forest Integrative Approach Based on Endogenous CYP3A4 Phenotype for Predicting Tacrolimus Concentrations and Dosages in Chinese Renal Transplant Patients. J. Clin. Pharm. Ther. 2020, 45, 318–323. [Google Scholar] [CrossRef]
  32. Wang, T.J.; Zhang, F.; Richards, J.B.; Kestenbaum, B.; van Meurs, J.B.; Berry, D.; Kiel, D.P.; Streeten, E.A.; Ohlsson, C.; Koller, D.L.; et al. Common Genetic Determinants of Vitamin D Insufficiency: A Genome-Wide Association Study. Lancet 2010, 376, 180–188. [Google Scholar] [CrossRef] [PubMed]
  33. Ahn, J.; Yu, K.; Stolzenberg-Solomon, R.; Simon, K.C.; McCullough, M.L.; Gallicchio, L.; Jacobs, E.J.; Ascherio, A.; Helzlsouer, K.; Jacobs, K.B.; et al. Genome-Wide Association Study of Circulating Vitamin D Levels. Hum. Mol. Genet. 2010, 19, 2739–2745. [Google Scholar] [CrossRef] [PubMed]
  34. Revez, J.A.; Lin, T.; Qiao, Z.; Xue, A.; Holtz, Y.; Zhu, Z.; Zeng, J.; Wang, H.; Sidorenko, J.; Kemper, K.E.; et al. Genome-Wide Association Study Identifies 143 Loci Associated with 25 Hydroxyvitamin D Concentration. Nat. Commun. 2020, 11, 1647. [Google Scholar] [CrossRef] [PubMed]
  35. Cheung, C.-L.; Lau, K.-S.; Sham, P.-C.; Tan, K.C.B.; Kung, A.W.C. Genetic Variant in Vitamin D Binding Protein Is Associated with Serum 25-Hydroxyvitamin D and Vitamin D Insufficiency in Southern Chinese. J. Hum. Genet. 2013, 58, 749–751. [Google Scholar] [CrossRef] [PubMed]
  36. Muindi, J.R.; Adjei, A.A.; Wu, Z.R.; Olson, I.; Huang, H.; Groman, A.; Tian, L.; Singh, P.K.; Sucheston, L.E.; Johnson, C.S.; et al. Serum Vitamin D Metabolites in Colorectal Cancer Patients Receiving Cholecalciferol Supplementation: Correlation with Polymorphisms in the Vitamin D Genes. Horm. Cancer 2013, 4, 242–250. [Google Scholar] [CrossRef]
  37. Suaini, N.H.A.; Koplin, J.J.; Ellis, J.A.; Peters, R.L.; Ponsonby, A.-L.; Dharmage, S.C.; Matheson, M.C.; Wake, M.; Panjari, M.; Tan, H.-T.T.; et al. Environmental and Genetic Determinants of Vitamin D Insufficiency in 12-Month-Old Infants. J. Steroid Biochem. Mol. Biol. 2014, 144 Pt B, 445–454. [Google Scholar] [CrossRef]
  38. Perna, L.; Felix, J.F.; Breitling, L.P.; Haug, U.; Raum, E.; Burwinkel, B.; Schöttker, B.; Brenner, H. Genetic Variations in the Vitamin D Binding Protein and Season-Specific Levels of Vitamin D among Older Adults. Epidemiology 2013, 24, 104–109. [Google Scholar] [CrossRef]
  39. Dong, J.; Zhou, Q.; Wang, J.; Lu, Y.; Li, J.; Wang, L.; Wang, L.; Meng, P.; Li, F.; Zhou, H.; et al. Association between Variants in Vitamin D-Binding Protein Gene and Vitamin D Deficiency among Pregnant Women in China. J. Clin. Lab. Anal. 2020, 34, e23376. [Google Scholar] [CrossRef]
  40. Ahn, J.; Albanes, D.; Berndt, S.I.; Peters, U.; Chatterjee, N.; Freedman, N.D.; Abnet, C.C.; Huang, W.-Y.; Kibel, A.S.; Crawford, E.D.; et al. Vitamin D-Related Genes, Serum Vitamin D Concentrations and Prostate Cancer Risk. Carcinogenesis 2009, 30, 769–776. [Google Scholar] [CrossRef]
  41. Anderson, D.; Holt, B.J.; Pennell, C.E.; Holt, P.G.; Hart, P.H.; Blackwell, J.M. Genome-Wide Association Study of Vitamin D Levels in Children: Replication in the Western Australian Pregnancy Cohort (Raine) Study. Genes. Immun. 2014, 15, 578–583. [Google Scholar] [CrossRef]
  42. O’Brien, K.M.; Sandler, D.P.; Shi, M.; Harmon, Q.E.; Taylor, J.A.; Weinberg, C.R. Genome-Wide Association Study of Serum 25-Hydroxyvitamin D in US Women. Front. Genet. 2018, 9, 67. [Google Scholar] [CrossRef] [PubMed]
  43. Jiang, X.; O’Reilly, P.F.; Aschard, H.; Hsu, Y.-H.; Richards, J.B.; Dupuis, J.; Ingelsson, E.; Karasik, D.; Pilz, S.; Berry, D.; et al. Genome-Wide Association Study in 79,366 European-Ancestry Individuals Informs the Genetic Architecture of 25-Hydroxyvitamin D Levels. Nat. Commun. 2018, 9, 260. [Google Scholar] [CrossRef] [PubMed]
  44. Duan, L.; Xue, Z.; Ji, H.; Zhang, D.; Wang, Y. Effects of CYP2R1 Gene Variants on Vitamin D Levels and Status: A Systematic Review and Meta-Analysis. Gene 2018, 678, 361–369. [Google Scholar] [CrossRef]
  45. Nissen, J.; Vogel, U.; Ravn-Haren, G.; Andersen, E.W.; Madsen, K.H.; Nexø, B.A.; Andersen, R.; Mejborn, H.; Bjerrum, P.J.; Rasmussen, L.B.; et al. Common Variants in CYP2R1 and GC Genes Are Both Determinants of Serum 25-Hydroxyvitamin D Concentrations after UVB Irradiation and after Consumption of Vitamin D₃-Fortified Bread and Milk during Winter in Denmark. Am. J. Clin. Nutr. 2015, 101, 218–227. [Google Scholar] [CrossRef]
  46. Lafi, Z.M.; Irshaid, Y.M.; El-Khateeb, M.; Ajlouni, K.M.; Hyassat, D. Association of Rs7041 and Rs4588 Polymorphisms of the Vitamin D Binding Protein and the Rs10741657 Polymorphism of CYP2R1 with Vitamin D Status Among Jordanian Patients. Genet. Test. Mol. Biomarkers 2015, 19, 629–636. [Google Scholar] [CrossRef] [PubMed]
  47. Zhu, Y.; Wang, P.P.; Zhai, G.; Bapat, B.; Savas, S.; Woodrow, J.R.; Campbell, P.T.; Li, Y.; Yang, N.; Zhou, X.; et al. Association of Rs2282679 A>C Polymorphism in Vitamin D Binding Protein Gene with Colorectal Cancer Risk and Survival: Effect Modification by Dietary Vitamin D Intake. BMC Cancer 2018, 18, 155. [Google Scholar] [CrossRef]
  48. Mondul, A.M.; Weinstein, S.J.; Moy, K.A.; Männistö, S.; Albanes, D. Vitamin D-binding Protein, Circulating Vitamin D and Risk of Renal Cell Carcinoma. Int. J. Cancer 2014, 134, 2699–2706. [Google Scholar] [CrossRef]
  49. Weinstein, S.J.; Mondul, A.M.; Kopp, W.; Rager, H.; Virtamo, J.; Albanes, D. Circulating 25-hydroxyvitamin D, Vitamin D-binding Protein and Risk of Prostate Cancer. Int. J. Cancer 2013, 132, 2940–2947. [Google Scholar] [CrossRef]
  50. Bouquegneau, A.; Salam, S.; Delanaye, P.; Eastell, R.; Khwaja, A. Bone Disease after Kidney Transplantation. Clin. J. Am. Soc. Nephrol. 2016, 11, 1282–1296. [Google Scholar] [CrossRef]
  51. Grant, W.B.; Lahore, H.; McDonnell, S.L.; Baggerly, C.A.; French, C.B.; Aliano, J.L.; Bhattoa, H.P. Evidence That Vitamin D Supplementation Could Reduce Risk of Influenza and COVID-19 Infections and Deaths. Nutrients 2020, 12, 988. [Google Scholar] [CrossRef] [PubMed]
  52. Rivera-Paredez, B.; Hidalgo-Bravo, A.; León-Reyes, G.; Antuna-Puente, B.; Flores, Y.N.; Salmerón, J.; Velázquez-Cruz, R. Association of GC Variants with Bone Mineral Density and Serum VDBP Concentrations in Mexican Population. Genes 2021, 12, 1176. [Google Scholar] [CrossRef] [PubMed]
  53. Goltzman, D. Functions of Vitamin D in Bone. Histochem. Cell Biol. 2018, 149, 305–312. [Google Scholar] [CrossRef]
  54. Ngai, M.; Lin, V.; Wong, H.C.; Vathsala, A.; How, P. Vitamin D Status and Its Association with Mineral and Bone Disorder in a Multi-Ethnic Chronic Kidney Disease Population. Clin. Nephrol. 2014, 82, 231–239. [Google Scholar] [CrossRef]
  55. Ezura, Y.; Nakajima, T.; Kajita, M.; Ishida, R.; Inoue, S.; Yoshida, H.; Suzuki, T.; Shiraki, M.; Hosoi, T.; Orimo, H.; et al. Association of Molecular Variants, Haplotypes, and Linkage Disequilibrium within the Human Vitamin D-Binding Protein (DBP) Gene with Postmenopausal Bone Mineral Density. J. Bone Miner. Res. 2003, 18, 1642–1649. [Google Scholar] [CrossRef]
  56. Ponticelli, C.; Passerini, P. Gastrointestinal Complications in Renal Transplant Recipients. Transpl. Int. 2005, 18, 643–650. [Google Scholar] [CrossRef] [PubMed]
  57. Kong, J.; Zhang, Z.; Musch, M.W.; Ning, G.; Sun, J.; Hart, J.; Bissonnette, M.; Li, Y.C. Novel Role of the Vitamin D Receptor in Maintaining the Integrity of the Intestinal Mucosal Barrier. Am. J. Physiol. Gastrointest. Liver Physiol. 2008, 294, G208–G216. [Google Scholar] [CrossRef] [PubMed]
  58. Zhang, Y.; Wu, S.; Lu, R.; Zhou, D.; Zhou, J.; Carmeliet, G.; Petrof, E.; Claud, E.C.; Sun, J. Tight Junction CLDN2 Gene Is a Direct Target of the Vitamin D Receptor. Sci. Rep. 2015, 5, 10642. [Google Scholar] [CrossRef]
  59. Nalle, S.C.; Turner, J.R. Intestinal Barrier Loss as a Critical Pathogenic Link between Inflammatory Bowel Disease and Graft-versus-Host Disease. Mucosal Immunol. 2015, 8, 720–730. [Google Scholar] [CrossRef]
  60. Legarth, C.; Grimm, D.; Wehland, M.; Bauer, J.; Krüger, M. The Impact of Vitamin D in the Treatment of Essential Hypertension. Int. J. Mol. Sci. 2018, 19, 455. [Google Scholar] [CrossRef]
  61. Victor, V.M.; Rocha, M.; Solá, E.; Bañuls, C.; Garcia-Malpartida, K.; Hernández-Mijares, A. Oxidative Stress, Endothelial Dysfunction and Atherosclerosis. Curr. Pharm. Des. 2009, 15, 2988–3002. [Google Scholar] [CrossRef] [PubMed]
  62. Kim, D.-H.; Meza, C.A.; Clarke, H.; Kim, J.-S.; Hickner, R.C. Vitamin D and Endothelial Function. Nutrients 2020, 12, 575. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, Q.-Y.; Jiang, C.-M.; Sun, C.; Tang, T.-F.; Jin, B.; Cao, D.-W.; He, J.-S.; Zhang, M. Hypovitaminosis D Is Associated with Endothelial Dysfunction in Patients with Non-Dialysis Chronic Kidney Disease. J. Nephrol. 2015, 28, 471–476. [Google Scholar] [CrossRef] [PubMed]
  64. Napoli, C.; de Nigris, F.; Williams-Ignarro, S.; Pignalosa, O.; Sica, V.; Ignarro, L.J. Nitric Oxide and Atherosclerosis: An Update. Nitric Oxide 2006, 15, 265–279. [Google Scholar] [CrossRef]
  65. Siasos, G.; Tousoulis, D.; Oikonomou, E.; Maniatis, K.; Kioufis, S.; Kokkou, E.; Vavuranakis, M.; Zaromitidou, M.; Kassi, E.; Miliou, A.; et al. Vitamin D3, D2 and Arterial Wall Properties in Coronary Artery Disease. Curr. Pharm. Des. 2014, 20, 5914–5918. [Google Scholar] [CrossRef]
  66. Campise, M.; Bamonti, F.; Novembrino, C.; Ippolito, S.; Tarantino, A.; Cornelli, U.; Lonati, S.; Cesana, B.M.; Ponticelli, C. Oxidative Stress in Kidney Transplant Patients. Transplantation 2003, 76, 1474–1478. [Google Scholar] [CrossRef]
  67. Hämäläinen, M.; Lahti, A.; Moilanen, E. Calcineurin Inhibitors, Cyclosporin A and Tacrolimus Inhibit Expression of Inducible Nitric Oxide Synthase in Colon Epithelial and Macrophage Cell Lines. Eur. J. Pharmacol. 2002, 448, 239–244. [Google Scholar] [CrossRef]
  68. Mangray, M.; Vella, J.P. Hypertension After Kidney Transplant. Am. J. Kidney Dis. 2011, 57, 331–341. [Google Scholar] [CrossRef] [PubMed]
  69. Quach, K.; Abdelmasih, M.; Chen, P.X.; Li, Y.; Famure, O.; Nash, M.; Prasad, R.; Perkins, B.A.; Yip, P.M.; Kim, S.J. Vitamin D Levels and the Risk of Posttransplant Diabetes Mellitus After Kidney Transplantation. Prog. Transplant. 2021, 31, 133–141. [Google Scholar] [CrossRef]
  70. Bland, R.; Markovic, D.; Hills, C.E.; Hughes, S.V.; Chan, S.L.F.; Squires, P.E.; Hewison, M. Expression of 25-Hydroxyvitamin D3-1alpha-Hydroxylase in Pancreatic Islets. J. Steroid Biochem. Mol. Biol. 2004, 89–90, 121–125. [Google Scholar] [CrossRef]
  71. Lindh, J.D.; Andersson, M.L.; Eliasson, E.; Björkhem-Bergman, L. Seasonal Variation in Blood Drug Concentrations and a Potential Relationship to Vitamin D. Drug Metab. Dispos. 2011, 39, 933–937. [Google Scholar] [CrossRef] [PubMed]
  72. Haussler, M.R.; Whitfield, G.K.; Haussler, C.A.; Hsieh, J.C.; Thompson, P.D.; Selznick, S.H.; Dominguez, C.E.; Jurutka, P.W. The Nuclear Vitamin D Receptor: Biological and Molecular Regulatory Properties Revealed. J. Bone Miner. Res. 1998, 13, 325–349. [Google Scholar] [CrossRef]
  73. Drocourt, L.; Ourlin, J.-C.; Pascussi, J.-M.; Maurel, P.; Vilarem, M.-J. Expression of CYP3A4, CYP2B6, and CYP2C9 Is Regulated by the Vitamin D Receptor Pathway in Primary Human Hepatocytes. J. Biol. Chem. 2002, 277, 25125–25132. [Google Scholar] [CrossRef]
  74. Prytuła, A.; Cransberg, K.; Raes, A. CYP3A4 Is a Crosslink between Vitamin D and Calcineurin Inhibitors in Solid Organ Transplant Recipients: Implications for Bone Health. Pharmacogenomics J. 2017, 17, 481–487. [Google Scholar] [CrossRef]
  75. Wang, Z.; Schuetz, E.G.; Xu, Y.; Thummel, K.E. Interplay between Vitamin D and the Drug Metabolizing Enzyme CYP3A4. J. Steroid Biochem. Mol. Biol. 2013, 136, 54–58. [Google Scholar] [CrossRef]
  76. Thummel, K.E.; Brimer, C.; Yasuda, K.; Thottassery, J.; Senn, T.; Lin, Y.; Ishizuka, H.; Kharasch, E.; Schuetz, J.; Schuetz, E. Transcriptional Control of Intestinal Cytochrome P-4503A by 1alpha,25-Dihydroxy Vitamin D3. Mol. Pharmacol. 2001, 60, 1399–1406. [Google Scholar] [CrossRef]
  77. Qin, X.; Wang, X. Role of Vitamin D Receptor in the Regulation of CYP3A Gene Expression. Acta Pharm. Sin. B 2019, 9, 1087–1098. [Google Scholar] [CrossRef] [PubMed]
  78. Van Etten, E.; Branisteanu, D.D.; Verstuyf, A.; Waer, M.; Bouillon, R.; Mathieu, C. Analogs of 1,25-Dihydroxyvitamin D3 as Dose-Reducing Agents for Classical Immunosuppressants. Transplantation 2000, 69, 1932–1942. [Google Scholar] [CrossRef] [PubMed]
  79. Rigby, W.F.; Denome, S.; Fanger, M.W. Regulation of Lymphokine Production and Human T Lymphocyte Activation by 1,25-Dihydroxyvitamin D3. Specific Inhibition at the Level of Messenger RNA. J. Clin. Investig. 1987, 79, 1659–1664. [Google Scholar] [CrossRef]
  80. Tokuda, N.; Mizuki, N.; Kasahara, M.; Levy, R.B. 1,25-Dihydroxyvitamin D3 down-Regulation of HLA-DR on Human Peripheral Blood Monocytes. Immunology 1992, 75, 349–354. [Google Scholar]
  81. Branisteanu, D.D.; Mathieu, C.; Bouillon, R. Synergism between Sirolimus and 1,25-Dihydroxyvitamin D3 in Vitro and in Vivo. J. Neuroimmunol. 1997, 79, 138–147. [Google Scholar] [CrossRef]
  82. Mathieu, C.; Waer, M.; Laureys, J.; Rutgeerts, O.; Bouillon, R. Activated Form of Vitamin D [1,25(OH)2D3] and Its Analogs Are Dose-Reducing Agents for Cyclosporine in Vitro and in Vivo. Transplant. Proc. 1994, 26, 3048–3049. [Google Scholar]
  83. Mathieu, C.; Bouillon, R.; Rutgeerts, O.; Vandeputte, M.; Waer, M. Potential Role of 1,25(OH)2 Vitamin D3 as a Dose-Reducing Agent for Cyclosporine and FK 506. Transplant. Proc. 1994, 26, 3130. [Google Scholar]
  84. Gim, J.-A.; Kwon, Y.; Lee, H.A.; Lee, K.-R.; Kim, S.; Choi, Y.; Kim, Y.K.; Lee, H. A Machine Learning-Based Identification of Genes Affecting the Pharmacokinetics of Tacrolimus Using the DMETTM Plus Platform. Int. J. Mol. Sci. 2020, 21, 2517. [Google Scholar] [CrossRef] [PubMed]
  85. Wiebe, C.; Ho, J.; Gibson, I.W.; Rush, D.N.; Nickerson, P.W. Carpe Diem-Time to Transition from Empiric to Precision Medicine in Kidney Transplantation. Am. J. Transplant. 2018, 18, 1615–1625. [Google Scholar] [CrossRef] [PubMed]
  86. Kelava, T.; Turcic, P.; Markotic, A.; Ostojic, A.; Sisl, D.; Mrzljak, A. Importance of Genetic Polymorphisms in Liver Transplantation Outcomes. World J. Gastroenterol. 2020, 26, 1273–1285. [Google Scholar] [CrossRef]
  87. Battaglia, Y.; Cojocaru, E.; Fiorini, F.; Granata, A.; Esposito, P.; Russo, L.; Bortoluzzi, A.; Storari, A.; Russo, D. Vitamin D in kidney transplant recipients. Clin. Nephrol. 2020, 93, 57–64. [Google Scholar] [CrossRef] [PubMed]
Figure 1. C&RT regression decision tree model for C/D ratio of cyclosporin A and the analyzed SNP. Model parameters: node count: 11; knot counts 15; trim measure: variance; standard error: 0.01; resubstitution cost 0.07; SK cost 0.09. Terminal leaves are marked in red, and four of these are characterized with at least one SNP. Leaves’ IDs, mean and variance values are reported inside the leaves and nodes. The genotypes relevant for particular leaves are described below.
Figure 1. C&RT regression decision tree model for C/D ratio of cyclosporin A and the analyzed SNP. Model parameters: node count: 11; knot counts 15; trim measure: variance; standard error: 0.01; resubstitution cost 0.07; SK cost 0.09. Terminal leaves are marked in red, and four of these are characterized with at least one SNP. Leaves’ IDs, mean and variance values are reported inside the leaves and nodes. The genotypes relevant for particular leaves are described below.
Jcm 13 04966 g001
Figure 2. ANOVA results for C/D CsA. ANOVA analysis of the mean value of C/D ratio ±95% confidence interval belonging to the final four leaves of the regression model shown in Figure 1. The results show the statistically significant difference between mean C/D ratio value of recipients with polymorphisms: leaf ID 10: rs2282679 CC or AA, rs602662 FUT2 GG or AG, IL-10-1082A>G RS1800896 AA; and leaf ID 9: rs2282679 CC or CA, rs602662 FUT2 AA, p = 0.006. Statistical analysis was performed using Statistica software (ver. 13.1).
Figure 2. ANOVA results for C/D CsA. ANOVA analysis of the mean value of C/D ratio ±95% confidence interval belonging to the final four leaves of the regression model shown in Figure 1. The results show the statistically significant difference between mean C/D ratio value of recipients with polymorphisms: leaf ID 10: rs2282679 CC or AA, rs602662 FUT2 GG or AG, IL-10-1082A>G RS1800896 AA; and leaf ID 9: rs2282679 CC or CA, rs602662 FUT2 AA, p = 0.006. Statistical analysis was performed using Statistica software (ver. 13.1).
Jcm 13 04966 g002
Table 1. Biochemical characteristics of the study group.
Table 1. Biochemical characteristics of the study group.
VariableMeanMedianMinimumMaximumStd. Dev.Conf. SD
−95%
Conf. SD
+95%
Standard
Error
WBC7.677.321.8517.312.432.202.710.18
RBC4.234.212.865.530.590.530.650.04
HGB12.8312.808.8017.101.701.541.900.13
HCT37.8937.8027.0051.404.814.355.370.36
PLT222.11213.0045.00467.0071.4664.6979.825.39
Na141.16141.00120.00149.003.573.233.990.27
K4.174.173.145.420.440.400.490.03
bilirubin0.610.540.161.690.290.260.320.02
BUN29.9826.458.3071.1013.0311.8014.560.98
creatinine1.611.440.674.310.640.580.710.05
uric acid6.987.003.7010.901.481.341.660.11
AlAT19.3316.003.0073.0011.1510.0912.450.84
AspAT20.6618.0010.0070.008.918.079.950.67
total CH187.81185.0080.00344.0039.0635.3243.712.99
HDL62.9960.6027.60118.2020.1217.8323.091.86
LDL95.6991.0036.00239.0036.4832.3241.903.39
TG142.87127.0048.00404.0066.4658.9076.266.14
lipids total639.37622.00430.001183.00133.75118.53153.4912.37
Table 2. Clinical characteristics of the study group.
Table 2. Clinical characteristics of the study group.
VariableMeanMedianMinimumMaximumStd. Dev.Conf. SD
−95%
Conf. SD
+95%
Standard
Error
Age [y]54.2755.0022.0082.0012.4111.2313.850.93
Weight [kg]78.8978.0040.00139.0017.2515.6219.271.30
Heihght [cm]169.59169.00150.00193.009.438.5410.540.71
BMI [kg/m2]27.3426.3717.3142.905.124.635.710.38
Time from KTx [m]10.3810.001.0027.006.245.656.970.47
Dose [mg]171.30150.0050.00500.0061.2755.4868.414.60
Concentration [ng/mL]116.37113.002.00315.9048.0343.4953.633.61
Concentration/dose0.720.670.011.750.280.250.310.02
Table 4. The differences in C/D ratio for CsA between particular genotypes in a post-hoc Tukey test. Statistical significances have been bolded.
Table 4. The differences in C/D ratio for CsA between particular genotypes in a post-hoc Tukey test. Statistical significances have been bolded.
Genotypers2282679GC CA
rs1800896 IL10 AA
(Leaf ID 5)
rs1800896 IL10 GG GA
rs2032582 MDR2677bGT
rs2282679 GC CC CA
(Leaf ID 6)
rs602662 FUT2 AA
rs2282679 GC
(Leaf ID 9)
rs2282679 GC CA
rs1800896 IL10 AA
p = 0.08p = 0.54
(Leaf ID 5)
rs1800896 IL10 GGGA
rs2032582 MDR2677bGT
p = 0.08 p = 0.006
rs2282679 GC CC or CA
(Leaf ID 6)
rs602662 FUT2 AA
p = 0.57p = 0.006
rs2282679 GC
(Leaf ID 9)
rs602662 FUT2 GG AG
rs1800896 IL10 AA
(Leaf ID 10)
p = 0.62p = 0.73p = 0.08
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kotowska, K.; Wojciuk, B.; Sieńko, J.; Bogacz, A.; Stukan, I.; Drożdżal, S.; Czerny, B.; Tejchman, K.; Trybek, G.; Machaliński, B.; et al. The Role of Vitamin D Metabolism Genes and Their Genomic Background in Shaping Cyclosporine A Dosage Parameters after Kidney Transplantation. J. Clin. Med. 2024, 13, 4966. https://doi.org/10.3390/jcm13164966

AMA Style

Kotowska K, Wojciuk B, Sieńko J, Bogacz A, Stukan I, Drożdżal S, Czerny B, Tejchman K, Trybek G, Machaliński B, et al. The Role of Vitamin D Metabolism Genes and Their Genomic Background in Shaping Cyclosporine A Dosage Parameters after Kidney Transplantation. Journal of Clinical Medicine. 2024; 13(16):4966. https://doi.org/10.3390/jcm13164966

Chicago/Turabian Style

Kotowska, Katarzyna, Bartosz Wojciuk, Jerzy Sieńko, Anna Bogacz, Iga Stukan, Sylwester Drożdżal, Bogusław Czerny, Karol Tejchman, Grzegorz Trybek, Bogusław Machaliński, and et al. 2024. "The Role of Vitamin D Metabolism Genes and Their Genomic Background in Shaping Cyclosporine A Dosage Parameters after Kidney Transplantation" Journal of Clinical Medicine 13, no. 16: 4966. https://doi.org/10.3390/jcm13164966

APA Style

Kotowska, K., Wojciuk, B., Sieńko, J., Bogacz, A., Stukan, I., Drożdżal, S., Czerny, B., Tejchman, K., Trybek, G., Machaliński, B., & Kotowski, M. (2024). The Role of Vitamin D Metabolism Genes and Their Genomic Background in Shaping Cyclosporine A Dosage Parameters after Kidney Transplantation. Journal of Clinical Medicine, 13(16), 4966. https://doi.org/10.3390/jcm13164966

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop