Skip to main content

SYSTEMATIC REVIEW article

Front. Cardiovasc. Med., 04 August 2022
Sec. Cardiovascular Genetics and Systems Medicine

Association between CES1 rs2244613 and the pharmacokinetics and safety of dabigatran: Meta-analysis and quantitative trait loci analysis

  • 1Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Beijing, China
  • 2National Center for Respiratory Medicine, Beijing, China
  • 3Institute of Respiratory Medicine, Chinese Academy of Medical Sciences, Beijing, China
  • 4National Clinical Research Center for Respiratory Diseases, Beijing, China
  • 5Graduate School of Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China
  • 6Shenzhen WeGene Clinical Laboratory, Shenzhen, China
  • 7WeGene, Shenzhen Zaozhidao Technology Co., Ltd., Shenzhen, China
  • 8Hunan Provincial Key Lab on Bioinformatics, School of Computer Science and Engineering, Central South University, Changsha, China
  • 9Graduate School of Capital Medical University, Beijing, China
  • 10Clinical Pharmacology Department, Hospital Universitario de La Princesa, Instituto de Investigación Sanitaria, Universidad Autónoma de Madrid, Madrid, Spain
  • 11Department of Clinical Research and Data Management, Center of Respiratory Medicine, China-Japan Friendship Hospital, Beijing, China
  • 12Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences, Peking Union Medical College, Beijing, China
  • 13VTT Technical Research Centre of Finland Ltd., Espoo, Finland
  • 14Department of Dermatology, China-Japan Friendship Hospital, Beijing, China
  • 15Department of Pharmacy, China-Japan Friendship Hospital, Beijing, China

Objective: To date, the influence of the carboxylesterase 1 (CES1) rs2244613 genotype on the pharmacokinetics (PKs) and safety of dabigatran remains controversial. Hence, a systematic review was performed to study the association between CES1 rs2244613 genotype and the PKs and safety of dabigatran and CES1 relative expression.

Methods: In addition to the three English databases (Web of Science, PubMed, and Embase), two Chinese databases (CNKI and Wanfang) were thoroughly revised. The mean differences (MD) and corresponding 95% confidence intervals (CI) were applied to evaluate the differences in PKs between the CES1 rs2244613 genotype. Odds ratio (OR) was used to study the risk for bleeding events between the CES1 rs2244613 genotypes. Subsequent expression quantitative trait loci (eQTL) analyses were performed to evaluate genotype-specific expressions in human tissues.

Results: Ten studies (n = 2,777) were included. CES1 rs2244613 G allele carriers exhibited significantly lower dabigatran trough concentrations compared to T allele carriers (MD: −8.00 ng/mL; 95% CI: −15.08 to −0.92; p = 0.03). The risk for bleeding events was significantly lower in carriers of the G allele compared to T allele carriers (OR: 0.65; 95% CI: 0.44–0.96; p = 0.03). Subsequent eQTL analysis showed significant genome-wide expressions in two human tissues, whole blood (p = 5.1 × 10–10) and liver (p = 6.2 × 10–43).

Conclusion: Our meta-analysis indicated a definite relation between the CES1 rs2244613 genotype and tolerability variations or pharmacokinetic fluctuations. The carriers of T allele showed higher dabigatran concentrations; therefore, they would benefit from a dose reduction.

Systematic review registration: [https://inplasy.com/inplasy-2022-6-0027/], identifier [NPLASY202260027].

GRAPHICAL ABSTRACT
www.frontiersin.org

GRAPHICAL ABSTRACT.

Introduction

Direct oral anticoagulants (DOACs) are the first alternative to vitamin K antagonists (VKAs) (1). They specifically target a single coagulation protein, including thrombin or coagulation factor Xa. Compared with traditional anticoagulants, the convenience and safety of DOACs is well documented (2). Dabigatran is a representative drug of DOACs widely used to treat atrial fibrillation and pulmonary embolism (3). It is administered as a prodrug–dabigatran etexilate–which is rapidly hydrolyzed into dabigatran, the active moiety, by means of esterases, such as carboxylesterase 1 (CES1) and CES2. Hepatic CES1 mainly catalyzes the conversion of the prodrug dabigatran etexilate to dabigatran, while the intestinal CES2 enzyme plays a compensatory role when CES1 is inhibited (4). This is the reason why we chose CES1 as the subject of this study.

CES1 is a crucial liver enzyme that conduces to the metabolism of drugs containing ester moieties, including dabigatran etexilate or the M1 metabolite (5, 6). As to treatment for atrial fibrillation, CES1 polymorphism may also affect clopidogrel pharmacological metabolism in the body. Up to 85% of the clopidogrel prodrug entering the body is rapidly hydrolyzed into inactive metabolites under the catalysis of CES1, and only 15% of the clopidogrel can exert drug effects. What’s more, CES1 is related to the development of many other thrombotic diseases like venous thromboembolism through regulating the pharmacokinetics of multiple anticoagulants (7, 8).

Single nucleotide polymorphisms (SNPs) in the CES1 gene may lead to interindividual differences in dabigatran pharmacokinetics (PKs), which may affect the metabolism and bioavailability of this drug. In addition, although the tolerability of dabigatran is better than that of VKAs, some serious adverse clinical events such as bleeding or thrombosis may occur.

Due to interindividual variability in PKs, bleeding or thrombotic events may occur in patients taking dabigatran. However, the conclusions of the existing studies on the association between the CES1 SNPs and drug concentration and bleeding risk are controversial due to their small sample sizes (4, 911). For instance, CES1 rs2244613 G allele was related to a reduction in the trough concentration of dabigatran in patients compared to the T allele, and with a reduced risk of bleeding (12, 13). However, Shi et al. (14) observed that this gene locus was unrelated to dabigatran concentration and clinical outcome.

Thence, a systematic review and meta-analysis were conducted with existing studies on the application of dabigatran in atrial fibrillation, cardioembolic stroke, knee arthroplasty, and other diseases. This study explores the relationship between the CES1 rs2244613 variant and patient’s plasma concentration and bleeding risk and determines its clinical relevance to guide individualized dabigatran prescription further.

Materials and methods

We performed this study in the light of the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines (Supplementary Table 1) (15). We have registered our detailed protocol for this systematic review on INPLASY (registration number: INPLASY202260027), and it is available in full on inplasy.com1.

Literature search

A structured search of three English databases (Web of Science, PubMed, and Embase) and two Chinese databases (CNKI and Wanfang) was performed on 16 April 2022. The search terms we applied are as follows: (‘novel oral anticoagulant’ or ‘new oral anticoagulant’ or ‘direct oral anticoagulant’ or ‘target-specific oral anticoagulant’ or NOAC or DOAC or TSOAC or dabigatran) and (CES1 or ‘carboxylesterase 1’ or carboxylesterase-1) and (‘dabigatran concentration’ or bleeding) and (polymorph* or variant* or mutation* or genotyp* or phenotyp* or haplotyp* or SNP or rs2244613).

Data selection and collection

With duplicate studies removed, two researchers (Li and Qiu) excluded irrelevant studies independently, according to the titles and abstracts and assessed the full-text articles for further inclusion. When inconsistencies occur, a team meeting was held with extra researchers, and a consensus would be finally reached.

In the step of data extraction, a predesigned form to obtain information from the included studies was used, which mainly comprised of basic data (including title, author, date, and sample size) and outcome variables (including means and standard bias for dabigatran plasma levels and the number of bleeding events). Then the means and standard deviations was estimated according to Wan’s method and presented the continuous outcomes in the form of medians and interquartile ranges (16).

Quality assessment

The Newcastle–Ottawa scale (NOS) tool, which is based on three domains including the selection of exposed and unexposed subjects (0–4 points), comparability of study groups (0–2 points), and outcome assessment (0–3 points), was used to evaluate the quality of the research (17).

Statistical analysis

The Review Manager software (version 5.3) and STATA software (version 12.0) were used. The MD, OR and 95% CI were used to evaluate the strength of the association. A total of five genetic models were implemented to make an assessment on the association between CES1 rs2244613 and dabigatran PKs and safety, including: homozygote model (GG vs. TT), heterozygote model (GT vs. TT), dominant model (GG + GT vs. TT), recessive model (GG vs. GT + TT), and allele comparison (G vs. T). The Q and I2 statistics were used to evaluate the heterogeneity degree (18). The selection of fixed-effects or random-effects model was based on the degree of heterogeneity (19). I2 < 50% was considered to low heterogeneity, 50 ≤ I2 < 75% was considered to moderate heterogeneity and I2 ≥ 75% was considered to significant heterogeneity. If I2 < 50% and p value > 0.1, the fixed-effects model would be used. If I2 ≥ 50% or P ≤ 0.1, the random-effects model would be used. Multiple populations were enrolled in the present meta-analysis. Therefore, we performed subgroup analysis and evaluated the impact of CES1 rs2244613 on the dabigatran pharmacokinetics and safety based on diverse ethnicities. To validate the credibility of outcomes in this meta-analysis, a sensitivity analysis was performed to identify potentially influential studies. Furthermore, funnel plot and Egger’s test were applied to detect publication bias (20). The funnel plot depends on whether the points on both sides are symmetric, which indicates a possible publication bias. And Egger’s test depends on the Student’s t-test (p < 0.05 suggests a publication bias).

Genotype quantitative trait loci analysis for rs2244613 in human tissues

We assessed the genotype-specific expression of CES1 in 49 human tissues by cis-expression quantitative trait loci (cis-eQTL) and splicing quantitative trait loci (sQTL) analysis through the Genotype-Tissue Expression (GTEx) portal2 (21). Violin plots of the genotype-specific expression were constructed to visualize normalized gene expressions between three variant genotypes (GG, GT, and TT).

Results

Search results and patient characteristics

Fifty four studies were included after the preliminary search, 35 of which remained after removing duplicates. Of 25 removed after full text revision, three were reviews, seven were case reports, six for evaluating other clinical outcomes, and nine for not providing extractable data (Supplementary Figure 1). Finally, ten studies (8, 12, 13, 2228) involving 2,777 subjects were included: Table 1 summarizes the characteristics of them. The earliest year of included literature is 2013, and the latest year is 2021.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of studies included in the systematic review and meta-analysis.

Seven of the included works analyzed the trough plasma concentration of dabigatran in patients with different genotypes, and nine analyzed the bleeding risk. Six of them were conducted with a Caucasian population and four with Asian populations. All publications were evaluated by NOS and scored above seven points.

Association between CES1 rs2244613 and the trough plasma concentration of dabigatran

Meta-analysis showed a statistically significant difference between trough plasma concentrations of dabigatran and rs2244613 genotype. In summary, the CES1 rs2244613 G allele was related to a lower trough plasma concentration of dabigatran when compared with T allele. The following MDs were observed for each model: GG vs. TT, MD = −58.29 ng/mL, 95% CI: −98.64 to −17.94, P = 0.005, I2 = 98%; GT vs. TT: MD = −10.14 ng/mL, 95% CI: −13.21 to −7.07, P < 0.00001, I2 = 0%; GG + GT vs. TT: MD = −12.56 ng/mL, 95% CI: −15.59 to −9.52, P < 0.00001, I2 = 0%; GG vs. GT + TT: MD = −44.86 ng/mL, 95% CI: −79.84 to −9.87, P = 0.01, I2 = 98%; G vs. T: MD = −8.00 ng/mL, 95% CI: −15.08 to −0.92, P = 0.03, I2 = 68% (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. Forest plots comparing dabigatran trough concentrations across rs2244613 genotype polymorphism of the carboxylesterase 1 (CES1) gene: (A) homozygote model, (B) heterozygote model, (C) dominant model, (D) recessive model, and (E) allelic model.

Significant heterogeneity was found for the homozygote model (I2 = 98%, Figure 1), for the recessive model (I2 = 98%, Figure 1), and for the allele contrast model (I2 = 68%, Figure 1). The heterogeneity was lower in Asian population in the homozygote model (I2 = 58%, Figure 2), recessive model (I2 = 67%, Figure 2), and allele contrast model (I2 = 53%, Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Subgroup analyses for the association between carboxylesterase 1 (CES1) rs2244613 polymorphism and the trough plasma concentration of dabigatran: (A) homozygote model, (B) recessive model, and (C) allelic model.

No single study could not influence the overall results qualitatively, indicating robustness and reliability of our results (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Sensitive analysis of carboxylesterase 1 (CES1) rs2244613 polymorphism and the trough plasma concentration of dabigatran: (A) homozygote model, (B) recessive model, and (C) allelic model.

No publication bias was observed, as funnel plots (Figure 4) were relatively symmetrical.

FIGURE 4
www.frontiersin.org

Figure 4. Funnel plots to assess the publication bias related to the dabigatran trough concentrations: (A) homozygote model, (B) heterozygote model, (C) dominant model, (D) recessive model, and (E) allelic model.

Association between CES1 rs2244613 and the risk of bleeding

Meta-analysis showed a statistically significant difference between the risk of developing bleeding and rs2244613 genotype. In summary, the CES1 rs2244613 G allele was related to a lower risk of developing any bleeding when compared with T allele. The following ORs were observed for each model: GG vs. TT, OR = 0.84, 95% CI: 0.40–1.77, P = 0.65, I2 = 40%; GT vs. TT: OR = 0.70, 95% CI: 0.40–1.24, P = 0.22, I2 = 0%; GG + GT vs. TT: OR = 0.64, 95% CI: 0.52–0.78, P < 0.0001, I2 = 0%; GG vs. GT + TT: OR = 0.53, 95% CI: 0.31–0.92, P = 0.02, I2 = 0%; G vs. T: OR = 0.65, 95% CI: 0.44–0.96, P = 0.03, I2 = 0% (Figure 5).

FIGURE 5
www.frontiersin.org

Figure 5. Forest plots comparing the risk of developing any bleeding event across rs2244613 genotype polymorphism of the carboxylesterase 1 (CES1) gene: (A) homozygote model, (B) heterozygote model, (C) dominant model, (D) recessive model, and (E) allelic model.

No publication bias was observed, as funnel plots (Figure 6) were relatively symmetrical.

FIGURE 6
www.frontiersin.org

Figure 6. Funnel plots to assess the publication bias related to the risk of bleeding: (A) homozygote model, (B) heterozygote model, (C) dominant model, (D) recessive model, and (E) allelic model.

Quantitative trait loci analysis of rs2244613 in human tissues

Out of the total 49 genotypic cis-eQTL results for rs2244613, only one cis-eQTLs reached a genome-wide significance threshold in Figure 7A (p = 5.1 × 10–10 in whole blood tissue). Genome-wide cis-eQTLs were upregulated in whole blood tissues in Figure 7B (slope = 0.30). Compared to TT allele patients, the expression of CES1 was significantly lower in GG. sQTLs showed genome-wide significance in seventeen tissues (p < 5 × 10–8) in Figure 7C and Supplementary Figure 2. Particularly, finding the cis-eQTL and sQTLs genotypes implicated the rs2244613 variant as a transcriptional regulatory factor.

FIGURE 7
www.frontiersin.org

Figure 7. (A) Genotype cis-expression quantitative trait loci analysis for rs2244613 in 49 human tissues was obtained from the GTE database. (B) Violin plots of allele-specific cis-eQTLs according to rs2244613 genotypes in whole blood tissue in the Genotype-Tissue Expression (GTEx) dataset. (C) Violin plots of allele-specific splicing quantitative trait loci (sQTL) according to rs2244613 genotypes in liver tissue in the GTEx dataset.

Discussion

Our study comprehensively explored the application of dabigatran in atrial fibrillation, cardioembolic stroke, and knee replacement, and other diseases to explore the relationship between CES1 rs2244613 and dabigatran PKs and bleeding risk. 2,777 patients in 10 articles were included. We found that the bleeding risk of patients taking dabigatran with GG and GT genotypes was significantly lower than that of patients with TT genotype; the bleeding risk of patients with GG genotype was remarkably lower than that of patients with GT + TT genotypes. Moreover, the bleeding risk is lower in patients carrying the G allele compare to T allele carriers. Additionally, we consistently observed that the trough concentrations of dabigatran were notably lower in the G compared to the T allele. Therefore, we conclude that CES1 rs2244613 affects dabigatran plasma concentration and ADR incidence. Moreover, the effect of CES1 rs2244613 on the trough concentrations of dabigatran varied among ethnicities, which is consistent with previous works (29).

Mammalian CES belong to the α, β-hydrolase-fold protein superfamily, which can be divided into five categories in accordance with the homology of the amino acid sequences (CES1CES5). Both CES1 and CES2 are mainly involved in the metabolism of human drugs, and CES1 is mostly found in the human liver (27, 28, 30, 31). Once dabigatran etexilate enters the body, it must be hydrolyzed at two separate sites to form an active thrombin inhibitor. First, in the intestine, the carbamate group is hydrolyzed by CES2, while CES1 hydrolyses the ethyl ester part. After that, it can be converted into dabigatran, which has metabolic activity (5, 14). Then it binds to the specific site of thrombin, inhibiting thrombin activity and preventing fibrin formation, thereby exerting an anticoagulant effect (14).

In fact, apart from CES1 and CES2, there are some other genes encoding enzymes [e.g., UDP-glucuronosyltransferase gene (UGT) and cytochrome P450 gene (CYP)] and genes encoding transporters [e.g., ATP binding cassette subfamily gene (ABC) and solute carriers’ family gene (SLC)]. After oral administration, dabigatran binds to plasma proteins and is catalyzed by three UGTs (UGT1A9, UGT2B7, and UGT2B15) to form acyl glucuronic acid isomers, of which UGT2B15 contains the strongest effect. Particularly, dabigatran 1-O-acylglucuronide, a metabolite of dabigatran, exhibited anticoagulant activity comparable to the parent drug (32). In addition, cytochrome P450 (CYP2D6 and CYP3A5) may metabolize dabigatran after CES esterase’s converting dabigatran to the active moiety. Dabigatran is mainly excreted unchanged in urine (85%) and remains in feces (9). Genes encoding transporters are also reported. P-glycoprotein (P-gp) is a classical transporter encoded by the ABCB1 gene, and dabigatran is one of its substrates. The gene polymorphism of ABCB1 is considered being related to the pharmacokinetics and drug safety of dabigatran, and has been widely confirmed (33). In addition, SLC family transporters are also involved in the metabolism of dabigatran. For example, studies have shown that the SLC22A1 mutant haplotype has higher tmax and t1/2 with dabigatran than heterozygous and wild types, resulting in differences in the pharmacokinetics and safety of dabigatran among users of different genotypes (9).

High interindividual variability in plasma levels of dabigatran was reported, and the coefficient of variation of up to 30% for systemic exposure (34). Genetic variations in drug-metabolizing enzymes, receptors, and transporters have been identified as a major cause of interindividual variability in drug response, potentially leading to differences in responsiveness and adverse reactions to dabigatran therapy among individuals with different genotypes (35). Presently, thousands of SNPs are described in the CES1 gene, such as rs8192935, rs71647871, and rs2244613 (36). The allele frequency of CES1 rs2244613 was previously reported to be different in Chinese vs. Caucasian populations, with a G allele prevalence of 61.1% and 15.3–28.3%, respectively. Furthermore, CES1 rs2244613 G allele was previously associated with reduced trough concentrations and a decreased bleeding risk rather than peak drug concentrations (4, 13, 25, 37). Another study of patients with atrial fibrillation who received oral dabigatran also concluded that the CES1 SNP rs2244613 was remarkably in association with dabigatran trough concentrations (38). In summary, most conclusions in post researches are consistent with ours, except Xu et al. As a meta-analysis, our study has a large sample size and employs data on dabigatran in a variety of disease populations, only for the drug dabigatran rather than a specific disease, so it has a comparatively high reliability. The reason for the large discrepancy between Xu’s research conclusions and ours may be the limitation of their sample size.

This study still has the following limitations. First, the results of our study indicate that SNPs may directly affect the bleeding risk of dabigatran through an internal mechanism and may indirectly influence the occurrence of adverse events by changing the concentration. The specific mechanisms acquire further basic research. Secondly, this study did not control other factors except genotypes, and the heterogeneity cannot be ignored. Thirdly, the blood concentration of dabigatran used in this study is from a single test rather than the average concentration of multiple tests, which may exist to some extent by chance. Fourthly, we have not analyzed other variants within CES1 and CES2, meta-analysis of other variants will be done in the follow-up.

Conclusion

In summary, patients carrying at least one CES1 rs2244613 G allele are associated with decreased dabigatran trough concentrations and lower bleeding risk compared to non-carriers (i.e., with the T/T genotype). This work is of great relevance as it will help eventually in the guidance and individualization of dabigatran prescription.

Data availability statement

The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

XZ, ZZhai, and CW had full access to all the data in the study and took responsibility for the content of the manuscript. HW and ZZhang conceived and designed the study. HL and YQ integrated data, analyzed the data, and wrote the manuscript. GF provided methodological support. YZ, PZ, PY, and A-LV participated in editing of the manuscript. All authors were involved in the revision of the manuscript for important intellectual content and approved the final version.

Funding

This research was funded by Beijing Nova Program, grant number: Z211100002121057 and Elite Medical Professionals project of China-Japan Friendship Hospital, grant number: ZRJY2021-QM12. PZ was financed by a Margarita Salas-UAM postdoctoral contract.

Acknowledgments

We thank Yu Jiang from Peking Union Medical College; Xiaoye Li from Zhongshan Hospital, Fudan University; Sherzod Abdullaev from Research Institute of Molecular and Personalized Medicine; Heze Han from Beijing Tiantan Hospital, Capital Medical University, and Zhu Zhu and Jie Pan from the Second Affiliated Hospital of Soochow University for their inspiration to the study. We also thank the editor and reviewers for their comments and instruction.

Conflict of interest

HW was employed by the Shenzhen Zaozhidao Technology Co., Ltd., and A-LV was employed by the VTT Technical Research Centre of Finland Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.959916/full#supplementary-material

Supplementary Figure 1 | PRISMA flow diagram.

Supplementary Figure 2 | Violin plots of allele-specific sQTLs according to rs2244613 genotypes in 17 human tissues in the GTEx dataset.

Supplementary Table 1 | PRISMA checklist.

Footnotes

  1. ^ https://inplasy.com/inplasy-2022-6-0027/
  2. ^ https://gtexportal.org/home/

References

1. Chan N, Sobieraj-Teague M, Eikelboom JW. Direct oral anticoagulants: evidence and unresolved issues. Lancet. (2020) 396:1767–76. doi: 10.1016/S0140-6736(20)32439-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ruff CT, Giugliano RP, Braunwald E, Hoffman EB, Deenadayalu N, Ezekowitz MD, et al. Comparison of the efficacy and safety of new oral anticoagulants with warfarin in patients with atrial fibrillation: a meta-analysis of randomised trials. Lancet. (2014) 383:955–62. doi: 10.1016/S0140-6736(13)62343-0

CrossRef Full Text | Google Scholar

3. Antonijevic NM, Zivkovic ID, Jovanovic LM, Matic DM, Kocica MJ, Mrdovic IB, et al. Dabigatran-metabolism, pharmacologic properties and drug interactions. Curr Drug Metab. (2017) 18:622–35. doi: 10.2174/1389200218666170427113504

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Liu Y, Yang C, Qi W, Pei Z, Xue W, Zhu H, et al. The impact of ABCB1 and CES1 polymorphisms on dabigatran pharmacokinetics in healthy chinese subjects. Pharmgenomics Pers Med. (2021) 14:477–85. doi: 10.2147/PGPM.S291723

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Laizure SC, Parker RB, Herring VL, Hu ZY. Identification of carboxylesterase-dependent dabigatran etexilate hydrolysis. Drug Metab Dispos. (2014) 42:201–6. doi: 10.1124/dmd.113.054353

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Shi J, Xiao J, Wang X, Jung SM, Bleske BE, Markowitz JS, et al. Plasma carboxylesterase 1 predicts methylphenidate exposure: a proof-of-concept study using plasma protein biomarker for hepatic drug metabolism. Clin Pharmacol Ther. (2022) 111:878–85. doi: 10.1002/cpt.2486

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Baturina O, Andreev D, Fedina L, Mirzaev K, Ivashchenko D, Ryzhikova K, et al. Influence of clinically significant genes on antiplatelet effect of clopidogrel and clinical outcomes in patients with acute coronary syndrome and atrial fibrillation. Pharmacology. (2022) 107:216–26. doi: 10.1159/000521531

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sychev DA, Levanov AN, Shelekhova TV, Bochkov PO, Denisenko NP, Ryzhikova KA, et al. The impact of ABCB1 (rs1045642 and rs4148738) and CES1 (rs2244613) gene polymorphisms on dabigatran equilibrium peak concentration in patients after total knee arthroplasty. Pharmacogenomics Pers Med. (2018) 11:127–37. doi: 10.2147/PGPM.S169277

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zubiaur P, Saiz-Rodríguez M, Ochoa D, Navares-Gómez M, Mejía G, Román M, et al. Effect of sex, use of pantoprazole and polymorphisms in SLC22A1, ABCB1, CES1, CYP3A5 and CYP2D6 on the pharmacokinetics and safety of dabigatran. Adv Ther. (2020) 37:3537–50. doi: 10.1007/s12325-020-01414-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chin PK, Wright DF, Zhang M, Wallace MC, Roberts RL, Patterson DM, et al. Correlation between trough plasma dabigatran concentrations and estimates of glomerular filtration rate based on creatinine and cystatin C. Drugs R D. (2014) 14:113–23. doi: 10.1007/s40268-014-0045-9

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Roşian A-N, Roşian ŞH, Kiss B, Ştefan MG, Trifa AP, Ober CD, et al. Genotype-phenotype correlation for dabigatran in patients with non-valvular atrial fibrillation (a single centre research). Hum Vet Med. (2020) 12: 123–9.

Google Scholar

12. Sychev D, Skripka A, Ryzhikova K, Bochkov P, Shevchenko R, Krupenin P, et al. Effect of CES1 and ABCB1 genotypes on the pharmacokinetics and clinical outcomes of dabigatran etexilate in patients with atrial fibrillation and chronic kidney disease. Drug Metab Pers Ther. (2020) 35:20190029. doi: 10.1515/dmpt-2019-0029

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Paré G, Eriksson N, Lehr T, Connolly S, Eikelboom J, Ezekowitz MD, et al. Genetic determinants of dabigatran plasma levels and their relation to bleeding. Circulation. (2013) 127:1404–12. doi: 10.1161/CIRCULATIONAHA.112.001233

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Shi J, Wang X, Nguyen JH, Bleske BE, Liang Y, Liu L, et al. Dabigatran etexilate activation is affected by the CES1 genetic polymorphism G143E (rs71647871) and gender. Biochem Pharmacol. (2016) 119:76–84. doi: 10.1016/j.bcp.2016.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Page MJ, McKenzie JE, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. The PRISMA 2020 statement: an updated guideline for reporting systematic reviews. BMJ. (2021) 372:n71. doi: 10.1136/bmj.n71

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Wan X, Wang W, Liu J, Tong T. Estimating the sample mean and standard deviation from the sample size, median, range and/or interquartile range. BMC Med Res Methodol. (2014) 14:135. doi: 10.1186/1471-2288-14-135

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wells GA, Shea B, O’Connell D, Peterson J, Welch V, Losos M, et al. The Newcastle-Ottawa Scale (NOS) for Assessing the Quality of Nonrandomised Studies in Meta-Analyses. (2002). Available online at: http://www.ohri.ca/programs/clinical_epidemiology/oxford.asp (accessed on 4 April, 2022).

Google Scholar

18. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. (2003) 327:557–60. doi: 10.1136/bmj.327.7414.557

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. (1959) 22:719–48.

Google Scholar

20. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. (1997) 315:629–34. doi: 10.1136/bmj.315.7109.629

PubMed Abstract | CrossRef Full Text | Google Scholar

21. GTEx Consortium. Human genomics. The genotype-tissue expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. (2015) 348:648–60.

Google Scholar

22. Meshcheryakov YV, Chertovskikh YV, Sychev DA. The pharmacogenetics of the new oral anticoagulant dabigatran - the role of rs2244613 CES1 polymorphism in the adverse reactions development. Pharmacogenet Pharmacogenomics. (2017) 2:18–9.

Google Scholar

23. Xu L. Study on TDM Monitoring, Hemorrhage Risk and Thrombosis Risk of Dabigatran Etexilate in NVAF Patients for Atrial Fibrillation Ablation in China (Chinese). Master’s thesis. Shanghai Jiao Tong University, Shanghai (2018).

Google Scholar

24. Tomek A, Jansky P, Olserova A, Cederquist L. The correlation of through plasmatic concentration of dabigatran and CES1 genotype with major bleeding in stroke patients. Stroke. (2018) 49:ATM110. doi: 10.1161/str.49.suppl_1.TMP110

CrossRef Full Text | Google Scholar

25. Ji Q, Zhang C, Xu Q, Wang Z, Li X, Lv Q. The impact of ABCB1 and CES1 polymorphisms on dabigatran pharmacokinetics and pharmacodynamics in patients with atrial fibrillation. Br J Clin Pharmacol. (2021) 87:2247–55. doi: 10.1111/bcp.14646

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lähteenmäki J, Vuorinen AL, Pajula J, Harno K, Lehto M, Niemi M, et al. Pharmacogenetics of bleeding and thromboembolic events in direct oral anticoagulant users. Clin Pharmacol Ther. (2021) 110:768–76. doi: 10.1002/cpt.2316

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zheng X, Zhang L, Zhang J. Effects of ABCB1 and CES1 genetic polymorphisms on the anticoagulant efficacy of dabigatran etexilate in patients with non-valvular atrial fibrillation (Chinese). Chin J Clin Lab Sci. (2021) 39:903–8.

Google Scholar

28. Xiang J. Correlation Study of ABCB1 and CES1 Gene Polymorphisms with Rivaroxaban/Dabigatran Plasma Concentration and Drug Safety in Patients with Atrial Fibrillation (Chinese). Master’s thesis, ChongQing Medical University, Chongqing (2021).

Google Scholar

29. Sychev DA, Abdullaev SP, Mirzaev KB, Ryzhikova KA, Shuyev GN, Sozaeva ZA, et al. Genetic determinants of dabigatran safety (CES1 gene rs2244613 polymorphism) in the Russian population: multi-ethnic analysis. Mol Biol Rep. (2019) 46:2761–9. doi: 10.1007/s11033-019-04722-w

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Hosokawa M, Furihata T, Yaginuma Y, Yamamoto N, Koyano N, Fujii A, et al. Genomic structure and transcriptional regulation of the rat, mouse, and human carboxylesterase genes. Drug Metab Rev. (2007) 39:1–15. doi: 10.1080/03602530600952164

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Laizure SC, Herring V, Hu Z, Witbrodt K, Parker RB. The role of human carboxylesterases in drug metabolism: have we overlooked their importance? Pharmacotherapy. (2013) 33:210–22. doi: 10.1002/phar.1194

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Ebner T, Wagner K, Wienen W. Dabigatran acylglucuronide, the major human metabolite of dabigatran: in vitro formation, stability, and pharmacological activity. Drug Metab Dispos. (2010) 38:1567–75. doi: 10.1124/dmd.110.033696

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Yamazaki S, Evers R, De Zwart L. Physiologically-based pharmacokinetic modeling to evaluate in vitro-to-in vivo extrapolation for intestinal P-glycoprotein inhibition. CPT Pharmacometrics Syst Pharmacol. (2022) 11:55–67. doi: 10.1002/psp4.12733

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Stangier J, Rathgen K, Stähle H, Gansser D, Roth W. The pharmacokinetics, pharmacodynamics and tolerability of dabigatran etexilate, a new oral direct thrombin inhibitor, in healthy male subjects. Br J Clin Pharmacol. (2007) 64:292–303. doi: 10.1111/j.1365-2125.2007.02899.x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Chen F, Zhang B, Parker RB, Laizure SC. Clinical implications of genetic variation in carboxylesterase drug metabolism. Expert Opin Drug Metab Toxicol. (2018) 14:131–42. doi: 10.1080/17425255.2018.1420164

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Raymond J, Imbert L, Cousin T, Duflot T, Varin R, Wils J, et al. Pharmacogenetics of direct oral anticoagulants: a systematic review. J Pers Med. (2021) 11:37. doi: 10.3390/jpm11010037

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Hamzic S, Kummer D, Milesi S, Mueller D, Joerger M, Aebi S, et al. Novel genetic variants in carboxylesterase 1 predict severe early-onset capecitabine-related toxicity. Clin Pharmacol Ther. (2017) 102:796–804. doi: 10.1002/cpt.641

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Dimatteo C, D’Andrea G, Vecchione G, Paoletti O, Cappucci F, Tiscia GL, et al. Pharmacogenetics of dabigatran etexilate interindividual variability. Thromb Res. (2016) 144:1–5. doi: 10.1016/j.thromres.2016.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: CES1, rs2244613, polymorphism, dabigatran, pharmacokinetics, safety, QTL

Citation: Li H, Zhang Z, Weng H, Qiu Y, Zubiaur P, Zhang Y, Fan G, Yang P, Vuorinen A-L, Zuo X, Zhai Z and Wang C (2022) Association between CES1 rs2244613 and the pharmacokinetics and safety of dabigatran: Meta-analysis and quantitative trait loci analysis. Front. Cardiovasc. Med. 9:959916. doi: 10.3389/fcvm.2022.959916

Received: 02 June 2022; Accepted: 29 June 2022;
Published: 04 August 2022.

Edited by:

Aamir Ahmad, University of Alabama at Birmingham, United States

Reviewed by:

Bihuan Cheng, Wenzhou Medical University, China
Sultan Kadasah, University of Bisha, Saudi Arabia
Jinhui Liu, Nanjing Medical University, China

Copyright © 2022 Li, Zhang, Weng, Qiu, Zubiaur, Zhang, Fan, Yang, Vuorinen, Zuo, Zhai and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xianbo Zuo, zuoxianbo@qq.com; Zhenguo Zhai, zhaizhenguo2011@126.com

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.