Association between SLCO1B1 genetic polymorphisms and bleeding risk in patients treated with edoxaban

Since SLCO1B1 encodes the uptake transporter OATP1B1, which can influence the pharmacokinetic and pharmacodynamic profiles of edoxaban, polymorphisms in SLCO1B1 may affect the edoxaban response. This study aimed to investigate the association between SLCO1B1 gene polymorphisms and the bleeding risk in patients receiving edoxaban. We genotyped 10 single-nucleotide polymorphisms (SNPs) from the SLCO1B1 gene in patients receiving edoxaban. We also analyzed rs3842 of ABCB1 as a confounder. The odds ratio (OR) and adjusted OR (AOR) were calculated from univariate and multivariable analysis, respectively. The area under the receiver operating characteristic curve (AUROC) was constructed for the discrimination of the model. A total of 159 patients receiving edoxaban were analyzed. Overdose and rs4149056 showed significant association with bleeding complications by around 11- and 5.5-fold, respectively. Additionally, patients with the rs4149057 variant allele (C) had a 3.9-fold increased bleeding risk compared with wild-type homozygote carriers (TT), whereas rs2306283 variant homozygote (GG) carriers had a 0.27-fold reduced bleeding risk compared with wild-type allele (A) carriers. Patients with the variant-type homozygote (CC) of ABCB1 rs3842 had a higher bleeding risk than T allele carriers (AOR = 5.3 and 5.9). The final models for multivariable analyses were acceptable based on the AUROC values (> 0.70). These findings may help predict bleeding risk in patients taking edoxaban and help personalize treatment.

Direct oral anticoagulants (DOACs) are gradually increasing in use due to their rapid anticoagulant effect and similar efficacy to vitamin K antagonists 1 .Despite the advantage of not requiring a routine international normalized ratio test, bleeding including gastrointestinal bleeding, hemorrhage, hematochezia, hematuria, and epistaxis is a major complication of DOAC use 2 .Especially, Asian individuals with DOACs showed a higher risk of bleeding, including intracranial hemorrhage than non-Asians and may require individualized DOAC treatment 3 .
Several studies in pharmacogenomics have focused on individualizing drug treatment with DOACs, primarily investigating genes associated with drug metabolism and transport.For example, in the case of dabigatran, which is the only DOAC that is not metabolized by CYP P450, it has been reported that the CES1 gene polymorphisms were associated with a lower dabigatran concentration due to the significant role of CES1 activity in drug metabolism 4,5 .For rivaroxaban, the impact of ABCB1 gene polymorphism on its clinical significance remains inconclusive.While rs2032582 homozygous mutated genotype was associated with rivaroxabaninduced hemorrhage 6 , rs2032582 was not associated with rivaroxaban pharmacokinetic parameters in healthy volunteers 7 .In the case of apixaban, rs2231142 of ABCG2 gene 8 , and rs776746 of CYP3A5 were associated with increased plasma concentrations 9 .Conversely, there is limited research on pharmacogenomics studies

Study patients and data collection
This study was a retrospective analysis of prospectively collected samples from June 2018 to December 2021.We recruited patients who had been using edoxaban and retrospectively collected the patients' previous records.We collected samples with patient consent on the day of the patient's first outpatient visit after the start of the study.It was conducted at Ewha Womans University Mokdong Hospital and Ewha Womans University Seoul Hospital.
The study subjects were individuals aged ≥ 20 years old who received edoxaban.Patients were excluded if they met the following criteria: (1) had thromboembolic or infarction-related events during the follow-up period, (2) experienced bleeding that was minor or unverified by health professionals while on treatment, (3) experienced bleeding after one year of edoxaban therapy, or (4) treated with edoxaban for less than three months (in the control group).We followed up for bleeding for one year after edoxaban initiation.The primary endpoint was any one-year major bleeding event and clinically relevant non-major bleeding (CRNMB) according to the International Society on Thrombosis and Haemostasis (ISTH) criteria 24,25 .
Among the study population, patients who experienced any one-year major bleeding event or CRNMB were classified as a case group and other patients were classified as a control group.
We obtained the data from electronic medical records.We collected patient demographic data, including sex, age, body mass index, creatinine clearance, prescription dose, concurrent medication, any history of myocardial infarction, stroke, transient ischemic attack, thromboembolism, or bleeding, comorbidities, smoking status, and alcohol status.The CHA 2 DS 2 -VASc stroke assessment score (congestive heart failure, hypertension, age ≥ 75 years, diabetes mellitus, stroke, vascular disease, age 65-74 years, and sex category; range 0-9) was calculated from its component variables 26 .The modified HAS-BLED (hypertension, abnormal renal or liver function, stroke, bleeding history or predisposition, elderly (age ≥ 65 years), concomitant drug and alcohol use; range 0-8) score, which is a specific risk score designed for bleeding risk assessment, were calculated from its component variables 27 .
The studies involving human participants were reviewed and approved by the Institutional Review Board of Ewha Womans University Mokdong Hospital and Ewha Womans University Seoul Hospital in accordance with the 1975 Helsinki Declaration and its later amendments (IRB numbers 2018-04-006 and 2019-05-038, respectively).Written informed consents from all participants were obtained before enrollment.

Selection of SNPs and genotyping
We selected SNPs through the following process.First, the Haploreg program was used to identify minor allele frequencies and linkage disequilibrium (LD) data of each SNP in the Asian population 28 .The tagger function within the Haploview of v4.2 was used to assign tag SNPs.Tag SNPs of the SLCO1B1 gene were assigned with a condition MAF ≥ 20% and an r 2 threshold of 0.8 in Asian populations.Second, among them, nine SNPs, including one synonymous SNP (rs4149057 29,30 ) and eight intronic SNPs (rs11045879 31,32 , rs12317268 33,34 , rs4149081 31,35 , rs999278 36 , rs2306283 13,37 , rs10841753 14,38 , rs2417957 39,40 , and rs4149042 41 ) of SLCO1B1 were selected based on previous studies and PharmGKB 42 , which is a pharmacogenomics Knowledge Base that offers information about how human genetic variation impacts drug response.Although rs4149056 (missense SNP) had a MAF of 0.13, it was included because it has been investigated for various drug-induced toxicity and studied with edoxaban 10,11,13,14,32,38 .Finally, a total of 10 SNPs were included in this study.We included ABCB1 rs3842, a 3'-untranslated region (UTR) SNP, as a known confounder based on previous studies 23,43,44 .

Statistical analysis
We compared continuous variables between patients who experienced bleeding and those who did not, using unpaired t-tests.The Kolmogorov-Smirnov method was used to test for the normality of the continuous variables.We analyzed categorical variables with the chi-squared test and Fisher's exact test.For genetic association analysis, we included both dominant and recessive models, and we selected the most appropriate model based on effect size and statistical significance.We identified independent risk factors for bleeding after adjusting for variables with P < 0.1 in the univariate analysis in addition to age, sex, and ABCB1 rs3842 via a multivariable logistic regression model.The unadjusted odds ratio (OR) and adjusted OR (AOR) with the 95% confidence interval (CI) were calculated from univariate and multivariable analyses, respectively.For the selection of the best model, the method of backward hierarchical elimination was used.To test the fit of the prediction model, we performed the Hosmer-Lemeshow goodness-of-fit test.We further evaluated the model discrimination by calculating the area under the receiver operating characteristic curve (AUROC).The predictive power of the logistic regression model was calculated.
Time to a bleeding event was analyzed using the Kaplan-Meier survival curves and the log-rank test.The Cox proportional-hazards model was sued for the multivariable analysis.Factors having P < 0.1 from the univariate analysis along with strong confounders of sex and age were included in the multivariable analysis.Hazard ratio (HR) and adjusted HR were calculated from the univariate and multivariable analyses, respectively.
All analyses were based on two-tail statistics and were performed using the Statistical Package for Social Sciences version 20.0 (IBM Corp., Armonk, NY, USA).P < 0.05 was considered statistically significant.

Results
We selected a total of 212 patients for the study, excluding 12 patients treated with edoxaban for < three months in the control group, 13 patients who experienced thromboembolic or infarction-related events during the follow-up period, 11 patients who reported minor bleeding during edoxaban treatment, 15 patients who had any bleeding at least one year after edoxaban therapy, one patient with a sample insufficient for DNA analysis, and one patient who withdrew informed consent (Fig. 1).Finally, 159 patients were included in the analysis.A total of 18 patients (11.3%) experienced bleeding complications, 7 of which were major and 11 of which were CRNMB events.The time (mean ± standard deviation) to a bleeding event was 111.72 ± 118.44 days.The indication for edoxaban treatment in the study patients was atrial fibrillation or secondary stroke prevention.Of the 18 patients, half the patients visited the emergency room or were hospitalized for bleeding.Two of them received red blood cell transfusions.All hospitalized patients recovered and were discharged.
Table 1 shows the demographic and clinical characteristics of the study population taking edoxaban.The mean age of the included patients was 71 years, and 92 patients (57.9%) were male.Approximately one-third of the patients received an underdose of edoxaban.The most common co-medications were statins, followed by beta-blockers.Approximately 97% of the patients had atrial fibrillation, and approximately 64% had hypertension.There was no significant factor for the incidence of bleeding complications.
In the genotype analysis, rs4149057, rs999278, rs2306283, and rs4149056 of SLCO1B1 and rs3842 of ABCB1 were significantly associated with bleeding risk (Table 2).
We performed a multivariable logistic regression analysis using variables with P < 0.1, age, and sex.We constructed two models for the multivariable analyses of the edoxaban subgroup, as one pair of SNPs (rs4149057 and rs999278) was in LD (r 2 = 0.97) in study population (Table 3).Model I included age, sex, overdose, rs999278, rs2306283, rs4149056, rs2417957, and rs3842.Model II included rs4149057 instead of rs999278 in Model I. Patients overdosed and variant-type allele (C) carriers of rs4149056 had an 11.0-11.4and 5.5-5.7 fold increased bleeding risk, respectively, after adjusting for confounders.Additionally, patients with rs4149057 variant allele (C) had a 3.9-fold increased bleeding risk compared with wild-type homozygote (TT) carriers, whereas rs2306283 variant homozygote (GG) carriers had a 0.27-fold reduced bleeding risk compared with wild-type allele (A) carriers.ABCB1 rs3842 variant homozygote carriers (CC) also showed a 5.3-5.9fold increased bleeding risk.
Vol:.( 1234567890 were also conducted (Table 3).Among them, SLCO1B1*1A and SLCO1B1*15 were significant for bleeding events.In the multivariable analysis, patients with SLCO1B1*15 showed an increased risk of bleeding events by 5.7-5.8folds, whereas SLCO1B1*1A failed to demonstrate statistical significance.RS999278 was a new SNP that increased the risk of bleeding complications by 3.7 times.
After Bonferroni correction of multiple comparisons for ten SNPs analyses, SLCO1B1 rs4149056 and ABCB1 rs3842 were the sole SNPs that reached statistical significance in the Table 2.When the simplest genetic model was constructed using these two SNPs, SLCO1B1 rs4149056 and ABCB1 rs3842 were associated with 4.7-and 5.7-fold increased bleeding risk, respectively.This genetic model still held a high predictive value (89.2%), similar to the predictive value (91.1%) of the prediction model shown initially in our study.
Since overdose can drive associations, further analyses were performed excluding overdose cases.Regardless of the model, two SLCO1B1 SNPs (rs2306283, rs4149056) and one ABCB1 SNP (rs3842) still showed statistical significance.Patients with the rs4149056 variant allele (C) had a 4.5-fold increased bleeding risk compared with wild-type homozygote carriers (TT), whereas rs2306283 variant homozygote (GG) carriers had a 0.25-fold reduced bleeding risk compared with wild-type allele (A) carriers.Patients with the variant-type homozygote (CC) of ABCB1 rs3842 had a higher bleeding risk than T allele carriers (AOR = 5.7).
Time to a bleeding event was analyzed.Among the clinical characteristics, there were no factors affecting the time to bleeding events (Supplementary Table 1).For genetic factors, SLCO1B1 rs4149056, rs4149057, rs999278 and ABCB1 rs3842 were significantly related with time to bleeding events in both univariate and multivariable analyses (Supplementary Tables 2 and 3).
Two models were constructed for the multivariable analyses, as one pair of SNPs (rs4149057 and rs999278) was in LD (r 2 = 0.97) in the study population (Supplementary Table 3).Rs4149056 and rs4149057 were associated with an increased hazard of time to bleeding by around 3.9 times and 3.0 times, respectively.Patients with the rs999278 variant allele (A) had a 2.9-fold increased hazard of time to bleeding compared with wild-type homozygote carriers (CC).Patients with the variant-type homozygote (CC) of ABCB1 rs3842 had a 3.7-3.8-timesthe hazard of time to bleeding compared to T allele carriers.

Discussion
This study revealed the clinical and genetic risk factors associated with bleeding complications during edoxaban therapy.Edoxaban overdose, SLCO1B1 rs4149057, rs2306283, and rs4149056, and ABCB1 rs3842 were significant factors for bleeding risk.www.nature.com/scientificreports/OATP1B1, a drug transporter expressed in the liver, plays an important role in transporting drugs and endogenous substrates from the blood into the hepatocytes 48 .The SLCO1B1 gene polymorphism, which encodes OATP1B1, may affect transporter activity.Previous studies on the association between SLCO1B1 and DOACs have been mainly limited to edoxaban 10,11 .
Among the SLCO1B1 SNP rs4149057, rs2306283, and rs4149056 that were significant in our edoxaban study, rs4149056 is the only SNP whose association with edoxaban has been studied.In the case of rs4149056, the relationship with edoxaban pharmacokinetics was investigated, but the relationship with edoxaban pharmacodynamics including bleeding has not yet been studied.
Rs4149056 (SLCO1B1*5) is a functional polymorphism in exon 5 and a well-studied SNP associated with drug toxicities such as stain-induced myopathy and methotrexate toxicity 13,14 .Rs4149056 was typically associated www.nature.com/scientificreports/with reduced transporter activity, resulting in increased systemic drug exposure and increased toxicity risk 13,38 .
In our study, rs4149056 significantly increased the bleeding risk.A pharmacogenomic analysis that combined 14 phase 1 studies investigated the relationship between SLCO1B1 gene polymorphism (rs4149056) and edoxaban pharmacokinetics 10 .Although the effect of rs4149056 on overall edoxaban exposure was insignificant, carriers with variant allele (C) showed an increased exposure to M4, the most-abundant metabolite of edoxaban.M4 was considered an insignificant metabolite for the overall anticoagulant effect because it accounted for less than 10% of the total anticoagulant exposure 49 .In another edoxaban study on SLCO1B1*15 haplotype, the correlation between prothrombin time and M4 concentration was significant in non-valvular atrial fibrillation patients with repeated administration, and the clinical contribution of M4 was confirmed 11 .Unlike rs4149056, studies on effects of rs4149057 and rs2306283 in patients taking edoxaban are limited.In the case of rs2306283, the relationship between rs2306283 and statin-induced myopathy has been primarily investigated.However, the role of rs2306283 polymorphism in transport remains controversial.The variant allele (G) of rs2306283 was not associated with statin-induced myopathy 50 , but in another study, SLCO1B1*15 carriers (rs4149056 allele C and rs2306283 allele G) showed higher rosuvastatin plasma concentration than SLCO1B1 haplotype (rs4149056 allele T and rs2306283 allele A), which may increase the risk of statin-induced myopathy 51 .Although rs4149057 is a synonymous SNP located in exon 5 of SLCO1B1, it is associated with the clearance change of mitotane and irinotecan, substrates of OATP1B1 29,30 .This study suggests that they play a critical role in edoxaban therapy.Future studies should elucidate the mechanisms of these SNPs.
The ABCB1 gene, which encodes the P-gp efflux pump, has been one of the frequently studied genes 52 .In our previous pharmacogenomics study, we found ABCB1 rs3842 was a significant factor associated with increased bleeding risk in patients with DOACs 23 .ABCB1 rs3842 is a 3'-UTP SNP that might influence protein expression by disrupting or creating microRNA binding sites 22 .Previous research also have shown that carriers of the ABCB1 rs3842 variant allele (C) had high bioavailability or lower disease activity scores 43,44 .Accordingly, we analyzed including ABCB1 rs3842 as a confounder, and rs3842 was also found to be a factor affecting bleeding complications in this study.www.nature.com/scientificreports/In this study, we observed an elevated bleeding risk among patients who experienced an overdose of edoxaban.To determine edoxaban overdose, we considered dose adjustment criteria based on the patient's renal function.DOAC overdose increases the risk of bleeding events 53 .In a systematic review regarding off-label dosages of DOACs, the incidence of adverse reactions was higher in patients administered an overdose than in those administered the recommended dose 54 .DOAC overdosing including edoxaban was associated with not only bleeding events but also stroke/systematic embolism, all-hospitalization, and all-cause mortality.
Several limitations were observed in this study, including relatively small sample size and exclusive focus on Asian individuals residing in Korea, which could potentially restrict the generalizability of the findings.Although we collected samples prospectively, there is a risk of bias associated with a retrospective design because patient data were collected retrospectively.Furthermore, the follow-up period was limited to one year to eliminate the likelihood of non-drug-related bleeding.Since this study was conducted on edoxaban, a substrate of SLCO1B1, the results of this study cannot be applied to other drugs or healthy populations.However, to our knowledge, this is the first study to investigate the association between SLCO1B1 polymorphisms and edoxaban-related bleeding complications.The results must be validated in different populations to generalize and apply them in clinical practice.

Figure 2 .
Figure 2. The area under the receiver operating characteristics (AUROC) curves for bleeding complications for edoxaban users.

Table 1 .
Baseline characteristics of patients who administered edoxaban.The CHA 2 DS 2 -VAS C score is a point-based system used to stratify the risk of stroke in atrial fibrillation patients.It stands for congestive heart failure, hypertension, age, diabetes mellitus, stroke, vascular disease, and sex category.ACEIs angiotensin converting enzyme inhibitors, ARBs angiotensin II receptor blockers, BMI body mass index, CYP cytochrome P450 family, DOACs direct oral anticoagulants, NA not available, TIA transient ischemic attack.a Standard dose was defined according to the FDA-approved labeling.

Table 2 .
Effects of gene polymorphisms on bleeding complications in patients who administered edoxaban.ABCB1 ATP binding cassette subfamily B member 1.