Abstract
Background:
Vitamin K antagonists (VKA) are used as prophylaxis for thromboembolic events in patients with cardiovascular diseases. The most common VKA are warfarin and acenocoumarol. These drugs have a narrow therapeutic margin and high inter-individual response variability due to clinical and pharmacogenetic variables.
Objective:
The authors aim to develop an algorithm comprised of clinical and genetic factors to explain the variability in the therapeutic dose of acenocoumarol among Chilean patients
Methodology:
DNA was obtained from 304 patients as a discovery cohort with an international normalized ratio (INR) range of 2.0ā3.0. The non-genetic (demographic and clinical) variables were also recorded. Genotype analyses were performed using real-time PCR for VKORC1 (rs9923231), VKORC1 (rs7294), GGCx (rs11676382), CYP4F2 (rs2108622), ABCB1 (rs1045642), CYP2C9*2 (rs1799853), ApoE (rs429358), and CYP2C9*3 (rs1057910).
Results:
The clinical variables that significantly influenced the weekly therapeutic dose of VKA were age, sex, body mass index (BMI), and initial INR, collectively accounting for 19% of the variability, and the genetic variables with a significant impact were VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910), explaining for another 37% of the variability.
Conclusion:
We developed an algorithm that explains 49.99% of the variability in therapeutic VKA dosage in the Chilean population studied. Factors that significantly affected the dosage included VKORC1, CYP2C9*2, and CYP2C9*3 polymorphisms, as well as age, sex, BMI, and initial INR.
Introduction
Cardiovascular diseases (CVD) have been highlighted as a health priority by several institutions worldwide, including the World Health Organization (WHO) through its Global Action Plan for the Prevention and Control of Non-communicable Diseases 2013ā2020. Oral anticoagulants are indicated for patients who survive a cardiovascular disease in order to prevent new thromboembolic conditions (WHO, 2013).
Coumarin anticoagulants, also called vitamin K antagonists (VKA), include drugs such as warfarin, acenocoumarol, and phenprocoumon. VKA are highly effective antithrombotic agents used to prevent complications associated with atrial fibrillation, artificial heart valves, and thromboembolic diseases (such as deep venous thrombosis and pulmonary embolism). The most commonly used VKA are warfarin and acenocoumarol (Pengo et al., 2006; Guyatt et al., 2012; Verhoef et al., 2013; QuiƱones et al., 2015; Nieto et al., 2019).
To calibrate weekly VKA doses, physicians use prothrombin time [expressed as international normalized ratio (INR) value] empirically to make adjustments at each visit until the patient reaches therapeutic range. Time in therapeutic range (TTR) is then used to assess the quality of anticoagulation. TTR should be above 65% to protect the patient against thrombotic or hemorrhagic risk, according to the National Institute for Health and Care Excellence (NICE) guidelines, or above 70%, according to the European Consensus (NICE, 2014; Esteve-Pastor et al., 2018). Major bleeding is the most concerning adverse effect of anticoagulant therapy, but the risk of bleeding associated with VKA is difficult to estimate. Risk estimates vary according to study design, with an approximate annual incidence of 0.6% for fatal hemorrhage, 3% for major hemorrhage, 9.6% for minor hemorrhage, and 0.2ā0.4% for intracranial bleeding depending on the series (Ageno et al., 2012; Fitzmaurice et al., 2016; Haas et al., 2016; Bosch et al., 2017).
Few studies have characterized Chilean patients treated with coumarin derivatives; however, the GARFIELD-AF study, which assessed Chilean patients with atrial fibrillation (AF), reported in 2017 that the median TTR was 40% for 971 patients treated in several public hospitals and private clinics. In that sample, 36 patients (3.6%) had a cerebrovascular accident as an adverse event (CorbalƔn et al., 2017). According to the same study, the average number of days to reach the desired therapeutic range was 301.6 days, with a median of 204 days, among the patients treated at the Western Metropolitan Health Service (WMHS) in Santiago, Chile. The physicians at the facilities studied relied exclusively on INR values for dose adjustment. Furthermore, the median TTR was only 50% in these AF patients; this low value is worrisome as these patients are at an elevated risk of a new thrombotic pattern or a hemorrhagic complication while out of the therapeutic range (Nieto et al., 2019).
In recent decades, pharmacogenetic research has addressed the relationship between the genetic factors and the required doses of VKA. The most-studied polymorphisms include CYP2C9, VKORC1, CYP4F2, GGCx, and ABCB1 (Leschziner et al., 2007; Caldwell et al., 2008; De Oliveira Almeida et al., 2011; Johnson et al., 2011; Sun et al., 2015; Shendre et al., 2016). Several algorithms with genetic and non-genetic variables have been developed to calculate VKA dosage, improving the efficacy and safety of the treatment according to TTR results (Wu et al., 2008; Ramos et al., 2012; Borobia et al., 2012; Pirmohamed et al., 2013; Wypasek et al., 2014; Santos et al., 2015; Johnson et al., 2017; Galvez et al., 2018). These algorithms seem to be more precise when developed and applied in specific populations as the frequency of the polymorphisms described depends on ethnicity. Furthermore, consumption of green vegetables, which are rich in vitamin K, also varies by geographical location (Visser et al., 2005; Kocael et al., 2019).
Several algorithms have been published for acenocoumarol in diverse populations. Verde et al. (2010) constructed an āacenocoumarol-dose genotype scoreā based on the number of alleles associated with a higher required acenocoumarol dosage carried by each patient for each polymorphism. In addition, two algorithms that include demographic, clinical, and genetic factors have been published for Indian populations, with coefficients of determination of 41 and 61.5% (Rathore et al., 2012; Krishna-Kumar et al., 2015).
Two algorithms have been developed for European populations. The first, designed for a mixed European population, includes CYP2C9 and VKORC1 polymorphisms and clinical variables (age, sex, weight, height, and amiodarone use). The genetic components in the algorithm explained 52.6% of the dosage variance, and the non-genetic variables explained 23.7% (van Schie et al., 2011). The second algorithm was developed in a cohort of 973 patients undergoing anticoagulation therapy and includes clinical factors [age and body mass index (BMI)] and genetic variants (VKORC1, CYP2C9, and CYP4F2 polymorphisms). The genetic and clinical variables explained 50 and 16% of the variance in acenocoumarol dosage, respectively (Cerezo-Manchado et al., 2013).
The aim of this study was to generate a preliminary algorithm with clinical and genetic factors that explains the variability in the therapeutic dose of VKA in Chilean patients. In order to achieve that, we have investigated the association of relevant single nucleotide polymorphisms (SNPs) (Table 1) with acenocoumarol dosage.
Table 1
| Enzyme | Gene | SNP | Change | Effect on VKA dose |
|---|---|---|---|---|
| MDR1 | ABCB1 | rs1045642 | c.3435C > T, exon 26 p.Ile1145Ile silent | Decrease |
| CYP4F2 | CYP4F2 | rs2108622 | c.1297 C > T, exon 11 p.Val433Met missense | Increase |
| CYP2C9 | CYP2C9*2 | rs1799853 | c.3608C > T, exon 3 p.Arg144Cys missense | Decrease |
| CYP2C9*3 | rs1057910 | c.42614 A > C, exon 7 p.Ile359Leu missense | Decrease | |
| GGCx | GGCX | rs11676382 | c.2084+45 C > G Intron 14 | Decrease |
| VKORC1 | VKORC1 | rs9923231 | -1639 G > A promotor | Decrease |
| VKORC1 | rs7294 | 3730 G > A 3ā²UTR | Increase | |
| APOE | ApoE | rs429358 | T > C, exon 4 p.Arg176Cys missense | Decrease |
Genetic variants and their potential effect on vitamin K antagonists (VKA) dosage (modified from Visser et al., 2005; Kocael et al., 2019).
Materials and Methods
Study Design
A retrospective cohort study was carried out between March and December of 2018 among patients treated with acenocoumarol (Coarol, Andrómaco, Santiago, Chile) as an antithrombotic therapy at WMHS in the Santiago and Melipilla provinces of Chile. INR measurements were performed in a capillary sample using CoaguChek pro II® equipment (Roche Mannheim, Germany). The sample size was determined according to the frequency of the carriers with the variant allele carriers in the population under study using PS Power and Sample Size Calculations Version 3.0, January 2009, considering 80% power, α = 0.05, OR = 2.0, and the less frequent minor allele frequencies (MAF) SNP CYP2C9*3 (rs1057910), according to literature (EMBL-EBI, 2019). The minimum patient number obtained was of 284 patients.
Initial Dosage and Dose Adjustment and Frequency of INR Monitoring
The initial dose was one 4-mg tablet of acenocoumarol on day 1. On day 2, the dose was decreased to 50% (half a tablet). The INR [(PTtest/PTnormal)ISI] (Riley et al., 2000) was controlled on day 3; thus, if the INR was higher than 1.8, the dose was again reduced by 50%, and the patients were checked in 2 days for medical control to adjust the dose according to the INR results. The weekly therapeutic dose of acenocoumarol was modified at each control according to the INR value of the patient. The dosage for patients with INR ā¤1.5 was increased by 20%, those of patients with INR >1.5ā < 2.0 was increased by 5%, those of patients with INR >3.0ā3.5 was decreased by 5%, those of patients with INR >3.5ā < 6.0 had their dose discontinued and decreased by 15%, while patients with INR ā„6 had their dose suspended and controlled in 3 days to start again. For INR values within the therapeutic range, the patients were seen in 4 weeks for control (Ansell et al., 2008; Marcatto et al., 2018). Patients having three consecutive INR values in the therapeutic range (2.0ā3.0) at the same dose of acenocoumarol were included in this study.
Ethics Statement
The research was authorized by the Ethics Committees of the University of Chile Faculty of Medicine, Project 222-2015, and the WMHS, Protocol No. 027/2016.
Patient Data
Data were obtained from clinical centers and managed with the statistical module of the anticoagulant treatment dosing software TAONet (Roche Diagnostics, Mannheim, Germany).
Genotypic Analysis
Genomic DNA was isolated from the peripheral blood samples of the subjects using the High Pure PCR Template Preparation Kit (catalog number 11796828001; Roche Diagnostics GmbH, Mannheim, Germany). VKORC1 (rs9923231), VKORC1 (rs7294), GGCx (rs11676382), CYP4F2 (rs2108622), ApoE (rs429358), ABCB1 (rs1045642), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910) were analyzed using the TaqManĀ® SNP Genotyping Assay (catalog number 4362691, Thermo Fisher Scientific, Waltham, MA, United States) in a Stratagene Mx3000p real-time PCR system. For quality assurance purposes, we randomly choose 20% of the samples for (a) repetition of the analysis and (b) PCR-RFLP analysis for coincidence. When the analyses were not coincident, we excluded the samples.
Statistical Analyses
Data analysis was performed with STATA 15.0Ā® software. The ShapiroāWilk test was used to determine whether the sample had a normal distribution for the continuous variable weekly therapeutic dose (WTD; mg/week), that is, the acenocoumarol dosage at which patients were in the therapeutic range. The ladder command from STATA 15.0Ā® was used to choose the best normal distribution expression. Finally, a linear regression analysis was performed among genetic (SNPs) and non-genetic variables with the logarithm of the WTD in the therapeutic range (2.0ā3.0), incorporating adjustment variables (p-value > 0.05). The performance of the algorithm was evaluated by calculating the adjusted coefficient of determination (R2) that represents the variability explained by the model. HardyāWeinberg equilibrium (HWE) was tested for all SNPs using chi2 test.
Results
Characteristics of the Study Population
We enrolled 377 patients on oral anticoagulant treatment with acenocoumarol. A total of 72 patients were excluded for not having three consecutive INR values in the therapeutic range (2.0ā3.0) at the same dose of acenocoumarol, and one was excluded for concomitant treatment with amiodarone (Figure 1). The final sample included 304 patients. As shown in Table 2, 47.4% of the patients were women, and 52.6% were men. The average age was 65.01 ± 13.99 years. The CaucasianāAmerindian admixture was 9.8% Amerindian and 90.2% Caucasian for this population (AcuƱa et al., 2000). No patient had bleeding nor myocardial infarction/stroke during the study.
Figure 1

Patient recruitment and selection flowchart.
Table 2
| Characteristic | N (%) |
|---|---|
| Female | 144 (47.4%) |
| Male | 160 (52.6%) |
| Total | 304 (100%) |
| % mixed Caucasianāaboriginal ethnicity | 9.80% |
| Age ± SD [range] (years) | 65.01 ± 13.99 [22ā95] |
| Body mass index (BMI) (kg/m2) ± SD (median) | 29.2 ± 5.7 (28.4) |
| INR range | 2.0ā3.0 |
| Acenocoumarol | 100% |
| Weekly therapeutic dose of acenocoumarol (mg/week) ± SD | 14.6 ± 2.2 |
| Acenocoumarol dosage, mg/week; range (median) | 3.5ā46 (13) |
| Time to reach therapeutic range (days); average ± SD | 308 ± 343 |
| Time to reach therapeutic range (days); range (median) | 3ā353 (206) |
| Primary diagnosis | N (%) |
| Rhythm disorder | 156 (51.3%) |
| Venous thrombosis with/without pulmonary thromboembolism | 64 (21.1%) |
| Occlusive arterial disease | 24 (7.9%) |
| Stroke | 17 (5.6%) |
| Others | 43 (14.1%) |
| Total | 304 (100%) |
| Secondary diagnosis | N (%) |
| Arterial hypertension | 64 (25.6%) |
| Diabetes mellitus | 27 (10.8%) |
| Cardiomyopathy | 26 (10.4%) |
| Other | 88 (42.9%) |
| Total | 205 (100%) |
Characteristics of the study patients.
SD, standard deviation.
Genotype Distribution in the Study Population
The analysis of the HWE showed that only CYP2C9*3 (rs1057910) is not in HWE (chi2 = 4.67). All other SNPs, VKORC1 (rs9923231) (chi2 = 0.09), VKORC1 (rs7294) (chi2 = 0.62), GGCx (rs11676382) (chi2 = 1.2), CYP4F2 (rs2108622) (chi2 = 1.02), ApoE (rs429358) (chi2 = 0.08), ABCB1 (rs1045642) (chi2 = 0.68), and CYP2C9*2 (rs1799853) (chi2 = 2.33), are in HWE.
The MAF were as follows: 0.467 for VKORC1 (rs9923231), 0.311 for VKORC1 (rs7294), 0.229 for CYP4F2 (rs2108622), 0.081 for CYP2C9*2 (rs1799853), 0.041 for CYP2C9*3 (rs1057910), 0.036 for GGCx (rs11676382), 0.092 for APOE (rs429358), and 0.627 for ABCB1 (rs1045642). In Table 3, it is possible to see that the MAF of VKORC1 (rs9923231) is similar to that of the Colombian population, slightly higher than those of Spain, Puerto Rico, and Europe, higher than the African-American, and lower than East-Asian populations. VKORC1 (rs7294) MAF is similar to Spain, Puerto Rico, Colombia, and European populations, lower than the African-American, and higher than East Asia. CYP2C9*2 (rs1799853) MAF is slightly lower than Spain, Puerto Rico, Colombia, and Europe and higher than African-American and East Asians. CYP2C9*3 (rs1057910) MAF is similar to Puerto Rico and East Asia, lower than Spain, Colombia, and Europe and higher than African-Americans; CYP4F2 (rs2108622) MAF is similar to Puerto Rico, Colombia, Europe, and East Asia but lower than Spain and higher than African-Americans. ABCB1 (rs1045642) MAF is higher than all the populations described but near to the European MAF. GGCx (rs11676382) MAF is similar to Puerto Rico and Colombia, lower than Spain and Europe, and higher than African-American and East Asian. Finally, APOE (rs429358) MAF is similar to Puerto Rico and East Asia and lower than Spain, Colombia, African-American, and European.
Table 3
| Polymorphic variant | Chile (this study) | Spain | Puerto Rico | Colombia | African-American | European | East Asia | Effect on VKA dose |
|---|---|---|---|---|---|---|---|---|
| VKORC1 (rs9923231) | 0.467 | 0.360 | 0.389 | 0.420 | 0.054 | 0.388 | 0.885 | Decrease |
| VKORC1 (rs7294) | 0.311 | 0.355 | 0.341 | 0.356 | 0.454 | 0.366 | 0.112 | Increase |
| CYP2C9*2 (rs1799853) | 0.081 | 0.140 | 0.139 | 0.122 | 0.008 | 0.124 | 0.001 | Decrease |
| CYP2C9*3 (rs1057910) | 0.041 | 0.084 | 0.043 | 0.064 | 0.002 | 0.073 | 0.034 | Decrease |
| CYP4F2 (rs2108622) | 0.229 | 0.355 | 0.288 | 0.282 | 0.082 | 0.290 | 0.214 | Increase |
| ABCB1 (rs1045642) | 0.627 | 0.463 | 0.428 | 0.441 | 0.150 | 0.518 | 0.398 | Decrease |
| GGCx (rs11676382) | 0.036 | 0.093 | 0.034 | 0.037 | 0.002 | 0.094 | 0.000 | Decrease |
| APOE (rs429358) | 0.092 | 0.140 | 0.106 | 0.154 | 0.268 | 0.155 | 0.086 | Decrease |
Comparison of minor allele frequencies obtained in this study and those conducted among in Spain, Puerto Rico, Brazil, Colombia, and African American population (http://www.ensembl.org/Homo_sapiens/Info/Index; modified from Visser et al., 2005; Kocael et al., 2019).
Genotype and VKA Dose Ratio in the Study Population
The continuous variable WTD did not have a normal distribution. The logarithm of the WTD was determined to be the optimal transformation for this analysis. The clinical variables that influenced the logarithm of WTD were sex (p << 0.0001), age (p < 0.0001), BMI (p < 0.0002), and INR at the beginning of treatment (p < 0.0011). The pharmacogenetic variables, that is, the polymorphisms analyzed which influenced the logarithm of the WTD, were CYP2C9*2 (rs1799853) (p < 0.0342), CYP2C9*3 (rs1057910) (p < 0.0001), and VKORC1 (rs9923231) (p < 0.0001).
Algorithm for Acenocoumarol Dosing in the Chilean Population
After the clinical and pharmacogenetic variables that influenced the logarithm of WTD were selected, a linear regression was performed using only the clinical variables, resulting in a model with an R2 of 0.2013, adjusted R2 of 0.19, and model p-value of <0.0001. A linear regression was then performed using only the pharmacogenetic variables VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910), producing a model with a R2 of 0.3790, adjusted R2 of 0.3685, and model p-value <0.0001. Finally, a linear regression was performed using the genetic factors VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910) and the clinical variables age, sex, BMI, and initial INR to produce a single model explaining the variability in the logarithm of acenocoumarol WTD, with R2 of 0.5147, adjusted R2 of 0.4999, and model p-value <0.0001 (Table 4). Therefore, the algorithm equation developed in this study is the following:
Table 4
| N observed | 287 | ||||
|---|---|---|---|---|---|
| Model p-value* | < 0.0001 | ||||
| R2 | 0.5147 | ||||
| Adjusted R2 | 0.4999 | ||||
| Variable | Coefficient | Standard error | p-value* | CI (95%) | |
| Sex (men) | 0.1668786 | 0.0407027 | 0.000 | 0.0867528 | 0.2470045 |
| Age | -0.008101 | 0.001472 | 0.000 | -0.0109987 | -0.0052034 |
| Initial INR | -0.0547186 | 0.0168253 | 0.001 | -0.0878404 | -0.0215969 |
| BMI | 0.0125554 | 0.0035861 | 0.001 | 0.0054959 | 0.0196149 |
| CYP2C9*2 (rs1799853) | |||||
| *1/*2 | -0.1067491 | 0.0538426 | 0.048 | -0.2127418 | -0.0007565 |
| CYP2C9*3 (rs1057910) | |||||
| *1/*3 | -0.3227895 | 0.0806461 | 0.000 | -0.4815465 | -0.1640324 |
| *3/*3 | -0.7465348 | 0.2416193 | 0.002 | -1.222178 | -0.2708915 |
| VKORC1 (rs9923231) | |||||
| G/A | -0.2704925 | 0.0479039 | 0.000 | -0.3647945 | -0.1761906 |
| A/A | -0.7008277 | 0.0583063 | 0.000 | -0.8156074 | -0.586048 |
| Constant | 3.080551 | 0.1622701 | 0.000 | 2.761112 | 3.33999 |
Linear regression including genetic factors and clinical variables in a single model.
*P < 0.05 is considered as statistically significant.
BMI, body mass index; INR, international normalized ratio.
Log WTD = 3.081 + (0.167 Ć men) - (age Ć 0.0081) - (initial INR Ć 0.055) + (BMI Ć 0.013) - (CYP2C9*1/*2 Ć 0.107) - (CYP2C9*1/*3 Ć 0.323) - (CYP2C9*3/*3 Ć 0.746) - (VKORC1 G/A Ć 0.270) - (VKORC1 A/A Ć 0.701).
Discussion
The patients enrolled in this study are a representative sample of the Chilean population, which is predominantly AmerindianāCaucasian admixture (9.8% of this sample) (AcuƱa et al., 2000). The group had an average BMI of 29.2 and median BMI of 28.4, classifying these patients as overweight (BMI 25ā29.9). These data are consistent with the results of the most recent National Health Survey 2016ā2017 (NHS, 2018), which indicated that 39.8% of the population was overweight (43.3% of males and 36.4% of females) and that among Chilean people in the age range of our study patients (65 years or older), 41.2% were overweight and 34.5% were obese. Notably, the time to reach therapeutic range in these patients was 308 ± 343 days on average, with a median of 206 days (Table 2). The patients are at high risk for stroke or hemorrhage while out of the therapeutic range. As relying exclusively on INR for dose adjustment is known to delay the time to reach the therapeutic range, the proposed pharmacogenetic dosage algorithm might be quite useful for the Chilean population.
As noted above, other published algorithms differ in the number of variables included and the weight of these variables, as well as in the type of population and methods used to develop the predictive model. The clinical variables included in these algorithms differ essentially in terms of inclusion of non-genetic variables such as sex, BMI, and use of amiodarone or enzyme-inducer drugs. In terms of the genetic variants, all algorithms published to date have included CYP2C9 and VKORC1 polymorphisms, whereas CYP4F2 and ApoE are used only in some models. In addition, several models have been designed exclusively for patients with deep vein thrombosis and/or pulmonary embolism (Borobia et al., 2012), while others have included patient cohorts with phenprocoumon and acenocoumarol treatment. The genetic variables included CYP2C9 and VKORC1 genes and the clinical variables include weight, height, sex, age, and amiodarone use and explained up to 76% of stable dose (van Schie et al., 2011). Another algorithm for acenocoumarol included clinical factors (age, body mass index, and body surface area) and genetic variants (VKORC1, CYP2C9*, and CYP4F2 polymorphisms) and explained up to 50% of stable dose (Cerezo-Manchado et al., 2013).
The original Clinical Pharmacogenetics Implementation Consortium (CPIC) algorithm published in the United States accounted for 47% of warfarin dose variability and included the clinical variables age, amiodarone use, weight, height, use of CYP2C9 inducers, and race/ethnicity as well as the pharmacogenetic factors VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910). CPIC suggests using CYP2C9*5 (rs28371686) if the patient is African-American and added CYP4F2 as an optional factor in the most recent update in 2017 (Table 4) (Johnson et al., 2017). Three algorithms have been published for Latin American populations. The algorithm for the population of Puerto Rico (Ramos et al., 2012), which explained 51% of the variability in warfarin dosage, was performed only in men, includes non-genetic variables such as age, initial INR, and use of amiodarone and the genetic variables VKORC1 (rs9923231), CYP2C9*2 (rs1799853), CYP2C9*3 (rs1057910), and CYP2C9*5 (rs28371686). The CYP2C9*5 (rs28371686) polymorphism was included due to the presence of a large African-American component in this population (Table 5). The Brazilian algorithm, in turn, accounted for 40% of the variability in warfarin dosage and includes the non-genetic variables age, sex, use of amiodarone or CYP2C9 inducers, and self-declared race, which, according to the Brazilian census criteria, includes white, mixed race, or black. The genetic variables included were VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910) (Botton et al., 2011; Santos et al., 2015), similar to our model. Finally, the Colombian model explained 45.9% of the variability in warfarin dosage. This model included non-genetic variables (age, use of amiodarone, weight, height, use of CYP2C9 inducers, and race/ethnicity), and the genetic variables were the same as those in our model (Table 4) (Galvez et al., 2018).
Table 5
| Algorithm | Chile (this study) | Spain (Borobia et al., 2012) | Germany (van Schie et al, 2011) | CPIC (Johnson et al., 2017) | Puerto Rico (Ramos et al., 2012) | Brazil (Santos et al., 2015) | Colombia (Galvez et al., 2018) | |
|---|---|---|---|---|---|---|---|---|
| Drug | VKA | Acenocoumarol | Acenocoumarol | Acenocoumarol | Warfarin | Warfarin | Warfarin | Warfarin |
| Clinical variables | Age | X | X | X | X | X | X | X |
| Sex | X | X | X | |||||
| Initial INR | X | X | ||||||
| Amiodarone use | X | X | X | X | X | |||
| Weight | X | X | X | X | ||||
| Height | X | X | X | |||||
| Body mass index | X | X | ||||||
| CYP2C9 inducer use | X | X | X | X | X | |||
| Race/ethnicity | X | X | X | |||||
| % contribution to the final model | 19% | 22% | 23.7% | N.D | 19% | N.D | 15.9% | |
| Genetic variables | VKORC1 (rs9923231) | X | X | X | X | X | X | X |
| CYP2C9*2 (rs1799853) | X | X | X | X | X | X | X | |
| CYP2C9*3 (rs1057910) | X | X | X | X | X | X | X | |
| CYP2C9*5 (rs28371686) | African-Americans | X | ||||||
| ApoE (rs429358) | X | |||||||
| CYP4F2 (rs2108622) | X | Optional | ||||||
| % contribution to the final model | 36,9% | 39% | 52.6% | N.D | 32% | N.D | 30% | |
| Percentage of variability in VKA dosage explained | 49.99% | 60.6% | 76.3% | 47% | 51% | 40% | 45.9% | |
Comparison between the present study and the vitamin K antagonists (VKA) dosage algorithms published for other populations.
N.D., no data published.
All of the above models were developed for warfarin. In Chile, however, as in Spain, the Ministry of Health indicates that acenocoumarol should be used in preference to any other coumarin. The Spanish model included the non-genetic variables sex, age, BMI, and initial INR value, and the genetic variables were the VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910) polymorphisms. The model explained 60% of the variability in acenocoumarol dosage, similar to our results (49.99%). Our model showed no association between WTD and CYP4F2 or ApoE polymorphisms, which are included in the Spanish study (Borobia et al., 2012). Therefore, the final algorithm for this Amerindian-Caucasian admixture includes the genetic factors VKORC1 (rs9923231), CYP2C9*2 (rs1799853), and CYP2C9*3 (rs1057910) and non-genetic variables age, sex, BMI, and initial INR, explaining almost 50% of the variability in WTD in the Chilean population studied. The algorithm equation developed for the group of Chilean patients explains a similar percentage of dose variability as the Puerto Rican, Brazilian, and Colombian algorithms (Table 5).
There are a number of limitations in this study. A number of parameters that affect coumarin dosage were not included, such as smoking status and use of other concomitant medications. These are important factors to keep in mind when establishing a stabilized dosage of acenocoumarol. Moreover, as this is a discovery cohort (also called a derivation or retrospective cohort), the next step is to perform a clinical application of this algorithm in a well-designed prospective validation cohort (also called a test cohort) to obtain sensitivity, specificity, and, of course, predictive values (Tong et al., 2016) before the algorithm is used routinely for acenocoumarol dose adjustment in Chilean patients.
Conclusion
Establishing appropriate coumarin dosage is challenging due in part to significant inter-individual variability in the dose required to achieve a stable range of anticoagulation (INR 2.0ā3.0). Various genetic and non-genetic factors have been associated with coumarin dosage requirements, and pharmacogenetic-guided dosing algorithms for warfarin and acenocoumarol have been developed for mixed populations with different predictive values. Here, we have developed the first acenocoumarol dosage algorithm for this Chilean mixed AmerindianāCaucasian population, which explains about 50% of dose variability. After clinical validation, this algorithm could provide a new tool for adjusting VKA dosage, considerably improving TTR, and thereby reducing thrombotic and hemorrhagic risks in Chilean patients.
Funding
The authors declare that this study received funding from Roche Diagnostics Chile (Probes Donation). The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary files, further inquiries can be directed to the corresponding author.
Ethics statement
The research was authorized by the Ethics Committees of the University of Chile Faculty of Medicine, Project 222-2015, and the WMHS, Protocol No. 027/2016.
Author contributions
AR participated in the conception of the research, analysis of data, and writing of the manuscript. EN participated in the conception of the research and writing of the manuscript. MS, MR, GV, and FT conducted the experimental analyses. MB, AA, DC, JM, GB, PS, FM, GG, and PV facilitated the enrolment of patients. LQ participated in the conception of the research, the analysis of data, and the writing of the manuscript.
Acknowledgments
The authors wish to thank the patients from the Western Metropolitan Health Service in the Santiago and Melipilla provinces for their altruistic collaboration in pursuit of the common welfare. We also thank the Latin American Society of Pharmacogenomics and Personalized Medicine (SOLFAGEM) for sponsoring this article.
Conflict of interest
The 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.
References
1
AcuƱaM.LlopE.RothhmmerF. (2000). Genetic composition of Chilean population: rural communities of Elqui, Limari and Choapa valleys. Rev. Med. Chil.128, 593ā600.
2
AgenoW.GallusA. S.WittkowskyA.CrowtherM.HylekE. M.PalaretiG. (2012). Oral anticoagulant therapy: Antithrombotic Therapy and Prevention of Thrombosis, 9th ed: American College of Chest Physicians Evidence-Based Clinical Practice Guidelines. Chest141, e44Sāe488. doi:Ā 10.1378/chest.11-2292
3
AnsellJ.HirshJ.HylekE.JacobsonA.CrowtherM.PalaretiG. (2008). Pharmacology and management of the vitamin K antagonists: american college of chest physicians evidence-based clinical practice guidelines (8th edition). Chest133 (6 Suppl.), 160Sā198S. doi:Ā 10.1378/chest.08-0670
4
BorobiaA. M.LubomirovR.RamĆrezE.LorenzoA.CamposA.MuƱoz-RomoR.et al. (2012). An Acenocoumarol Dosing Algorithm Using Clinical and Pharmacogenetic Data in Spanish Patients with Thromboembolic Disease. PloS One7 (7), e41360. doi:Ā 10.1371/journal.pone.0041360
5
BoschJ.EikelboomJ. W.ConnollyS. J.BrunsN. C.LaniusV.YuanF.et al. (2017). Rationale, Design and Baseline Characteristics of Participants in the Cardiovascular Outcomes for People Using Anticoagulation Strategies (COMPASS) Trial. Can. J. Cardiol.33 (8), 1027ā1035. doi:Ā 10.1016/j.cjca.2017.06.001
6
BottonM. R.BandinelliE.RohdeL. E.AmonL. C.HutzM. H. (2011). Influence of genetic, biological and pharmacological factors on warfarin dose in a Southern Brazilian population of European ancestry. Br. J. Clin. Pharmacol.72 (3), 442ā450. doi:Ā 10.1111/j.1365-2125.2011.03942.x
7
CaldwellM. D.AwadT.JohnsonJ. A.GageB. F.FalkowskiM.GardinaP.et al. (2008). CYP4F2 genetic variant alters required warfarin dose. Blood111 (8), 4106ā4112. doi:Ā 10.1182/blood-2007-11-122010
8
Cerezo-ManchadoJ. J.RosafalcoM.AntónA. I.PĆ©rez-AndreuV.Garcia-BarberĆ”N.MartinezA. B.et al. (2013). Creating a genotype-based dosing algorithm for acenocoumarol steady dose. Thromb. Haemost.109 (1), 146ā153. doi:Ā 10.1160/TH12-08-0631
9
CorbalĆ”nR.ConejerosC.ReyC.StockinsB.EggersG.AstudilloC.et al. (2017). Features, management and prognosis of Chilean patients with non valvular atrial fibrillation: GARFIELD AF registry. Rev. Med. Chil.145 (8), 963ā971. doi:Ā 10.4067/s0034-98872017000800963
10
De Oliveira AlmeidaV. C.De Souza FerreiraA. C.RibeiroD. D.Gomes BorgesK. B.Salles Moura FernandesA. P.Brunialti GodardA. L. (2011). Association of the C3435T polymorphism of the MDR1 gene and therapeutic doses of warfarin in thrombophilic patients. J. Thromb. Haemost.9 (10), 2120ā2122. doi:Ā 10.1111/j.1538-7836.2011.04452.x
11
EMBL-EBI (2019). The European Bioinformatics Institute. rs1057910 SNP. Ensembl. 2016 Available from: http://www.ensembl.org/Homo_sapiens/Variation/Population?db=core;r=10:94980796-94981796;v=rs1057910;vdb=variation;vf=768358.
12
Esteve-PastorM. A.Rivera-CaravacaJ. M.RoldĆ”n-RabadĆ”nI.RoldĆ”nV.MuƱizJ.RaƱa-MĆguezP.et al. (2018). Quality of oral anticoagulation with vitamin K antagonists in āreal-worldā patients with atrial fibrillation: a report from the prospective multicentre FANTASIIA registry. EP Europace20 (9), 1435ā1441. doi:Ā 10.1093/europace/eux314
13
FitzmauriceD. A.AccettaG.HaasS.KayaniG.Lucas LuciardiH.MisselwitzF.et al.(2016). Comparison of international normalized ratio audit parameters inpatients enrolled in GARFIELD-AF and treated with vitamin K antagonists. Br. J. Haematol.174 (4), 610ā623. doi:Ā 10.1111/bjh,14084
14
GalvezJ. M.RestrepoC. M.ContrerasN. C. (2018). Creating and validating a warfarin pharmacogenetic dosing algorithm for Colombian patients. Pharmacogenomics pers. Med.11, 169ā178. doi:Ā 10.2147/PGPM.S170515
15
GuyattG. H.AklE. A.CrowtherM.GuttermanD. D.SchuünemannH. J. (2012). Executive summary: Antithrombotic therapy and prevention of thrombosis, 9th ed: American College of Chest Physicians evidence-based clinical practice guidelines. Chest141 (2 Suppl.), 7ā47. doi:Ā 10.1378/chest.1412S3
16
HaasS.Ten CateH.AccettaG.AngchaisuksiriP.BassandJ. P.CammA. J.et al.(2016). Quality of Vitamin K Antagonist Control and 1-Year Outcomes in Patients with Atrial Fibrillation: A Global Perspective from the GARFIELD-AF Registry. PloS One11 (10), e0164076. doi:Ā 10.1371/journal.pone.0164076
17
JohnsonJ. A.GongL.Whirl-CarrilloM.GageB. F.ScottS. A.SteinC. M.et al. (2011). Clinical Pharmacogenetics Implementation Consortium Guidelines for CYP2C9 and VKORC1 genotypes and warfarin dosing. Clin. Pharmacol. Ther.90 (4), 625ā629. doi: 10.1038/clpt.2011.185
18
JohnsonJ.CaudleK.GongL.Whirl-CarrilloM.SteinC. M.ScottS. A.et al. (2017). Supplement to: Clinical Pharmacogenetics Implementation Consortium (CPIC) Guidelines for Pharmacogenetics-guided Warfarin Dosing: 2016 Update. Clin. Pharmacol. Ther.102 (3), 397ā404. doi:Ā 10.1002/cpt.668
19
KocaelA.EronatA. P.TüzünerM. B. (2019). Interpretation of the effect of CYP2C9, VKORC1 and CYP4F2 variants on warfarin dosing adjustment in Turkey. Mol. Biol. Rep.46 (2), 1825ā1833. doi:Ā 10.1007/s11033-019-04634-9
20
Krishna-KumarD.ShewadeD. G.LoriotM. A.BeauneP.Sai ChandranB. V.BalachanderJ.et al. (2015). An acenocoumarol dosing algorithm exploiting clinical and genetic factors in South Indian (Dravidian) population. Eur. J. Clin. Pharmacol.71 (2), 173ā181. doi:Ā 10.1007/s00228-014-1791-x
21
LeschzinerG. D.AndrewT.PirmohamedM.JohnsonM. R. (2007). ABCB1 genotype and PGP expression, function and therapeutic drug response: a critical review and recommendations for future research. Pharmacogenomics J.7, 154ā179. doi: 10.1038/sj.tpj.6500413
22
MarcattoL. R.SacilottoL.TavaresL. C.FacinM.OlivettiN.StrunzC. M. C.et al. (2018). Pharmaceutical Care Increases Time in Therapeutic Range of Patients With Poor Quality of Anticoagulation With Warfarin. Front. Pharmacol.9, 1052. doi:Ā 10.3389/fphar.2018.01052. eCollection 2018.
23
NHS (2018). National Health Survey 2016-2017, Chile. Web site: http://epi.minsal.cl/wp-content/uploads/2018/03/1%C2%BA-Resultados-ENS_DEPTO.EPIDEMIOLOGIA.MINSAL.pdf.
24
NICE (2014). Atrial fibrillation: the management of atrial fibrillation (CG180), Issued: June 2014, last modified: August 2014. Web site: https://www.nice.org.uk/guidance/cg180.
25
NietoE.SuarezM.RocoA.RubilarJ. C.TamayoF.RojoM.et al. (2019). Anticoagulation Management With Coumarinic Drugs in Chilean Patients. Clin. Appl. Thrombosis/Hemostasis.25, 1ā7. doi:Ā 10.1177/1076029619834342
26
PengoV.PegoraroC.CucchiniU.IlicetoS. (2006). Worldwide management of oral anticoagulant therapy: The ISAM study. J. Thromb. Thrombolysis21 (1), 73ā77. doi:Ā 10.1007/s11239-006-5580-y
27
PirmohamedM.BurnsideG.ErikssonN.JorgensenA. L.TohC. H.NicholsonT.et al. (2013). A randomized trial of genotype-guided dosing of warfarin. N Engl. J. Med.369, 2294ā2303. doi:Ā 10.1056/NEJMoa1311386
28
QuiƱonesL.RocoA.MirandaC. (2015). Farmacogenómica: Aplicaciones Cardiovasculares. Rev. MĆ©dica ClĆnica las Condes196 (2), 198ā209. doi:Ā 10.1016/j.rmclc.2015.04.008
29
RamosA. S.SeipR. L.Rivera-MirandaG.Felici-GiovaniniM. E.Garcia-BerdeciaR.Alejandro-CowanY.et al. (2012). Development of a pharmacogenetic-guided warfarin dosing algorithm for Puerto Rican patients. Pharmacogenomics13 (16), 1937ā1950. doi:Ā 10.2217/pgs.12.171
30
RathoreS.AgarwalS. K.PandeS.SinghS. K.MittalT.MittalB. (2012). Therapeutic dosing of acenocoumarol: Proposal of a Population Specific Pharmacogenetic Dosing Algorithm and Its Validation in North Indians. PloS One7 (5), e37844. doi:Ā 10.1371/journal.pone.0037844
31
RileyR. S.RoweD.FisherL. M. (2000). Clinical utilization of the international normalized ratio (INR). J. Clin. Lab. Anal.14, 101ā114. doi:Ā 10.1002/(SICI)1098-2825(2000)14:3<101::AID-JCLA4>3.0.CO;2-A
32
SantosP. C.MarcattoL. R.DuarteN. E.Gadi SoaresR. A.Cassaro StrunzC. M.ScanavaccaM.et al. (2015). Development of a pharmacogenetic based warfarin dosing algorithm and its performance in Brazilian patients: highlighting the importance of population-specific calibration. Pharmacogenomics16 (8), 865ā876. doi: 10.2217/pgs.15.48
33
ShendreA.BrownT. M.LiuN. (2016). Race-Specific Influence of CYP4F2 on Dose and Risk of Hemorrhage Among Warfarin Users. Pharmacotherapy36 (3), 263ā272. doi:Ā 10.1002/phar.1717
34
SunY.WuZ.LiS. (2015). Impact of gamma-glutamyl carboxylase gene polymorphisms on warfarin dose requirement: a systematic review and meta-analysis. Thromb. Res.135 (4), 739ā747. doi:Ā 10.1016/j.thromres.2015.01.029
35
TongH. Y.DĆ”vila-FajardoC. L.BorobiaA. M.MartĆnez-GonzĆ”lezL. J.LubomirovR.Perea-LeónL. M.et al. (2016). A New Pharmacogenetic Algorithm to Predict the Most Appropriate Dosage of Acenocoumarol for Stable Anticoagulation in a Mixed Spanish Population. PloS One 201611 (3), e0150456. doi:Ā 10.1371/journal.pone.0150456. eCollection 2016.
36
van SchieR. M.WesselsJ. A.le CessieS.de BoerA.SchalekampT.van der MeerF. J.et al. (2011). Loading and maintenance dose algorithms for phenprocoumon and acenocoumarol using patient characteristics and pharmacogenetic data. Eur. Heart J.32 (15), 1909ā1917. doi:Ā 10.1093/eurheartj/ehr116
37
VerdeZ.RuizJ. R.SantiagoC.ValleB.BandrésF.CalvoE.et al. (2010). A novel, single algorithm approach to predict acenocoumarol dose based on CYP2C9 and VKORC1 allele variants. PloS One 2010;5 (6), e11210. doi: 10.1371/journal.pone.0011210
38
VerhoefT. I.RagiaG.de BoerA.BarallonR.KolovouG.KolovouV.et al. (2013). A randomized trial of genotype-guided dosing of acenocoumarol and phenprocoumon. N Engl. J. Med.369, 2304ā2312. doi:Ā 10.1056/NEJMoa1311388
39
VisserL. E.TrienekensP. H.De SmetP. A.VultoA. G.HofmanA.van DuijnC. M.et al. (2005). Patients with an APOE epsilon4 allele require lower doses of coumarin anticoagulants. Pharmacogenet. Genomics15 (2), 69ā74. doi: 10.1097/01213011-200502000-00002
40
WHO (2013). Global Action Plan for the Prevention and Control of NCDs 2013-2020. 05/03/2019, web site: https://apps.who.int/iris/bitstream/10665/94384/1/9789241506236_eng.pdf?ua=1.
41
WuA. H.WangP.SmithA.HallerC.DrakeK.LinderM.et al. (2008). Dosing algorithm for warfarin using CYP2C9 and VKORC1 genotyping from a multi-ethnic population: comparison with other equations. Pharmacogenomics9 (2), 169ā178. doi:Ā 10.2217/14622416.9.2.169
42
WypasekE.BranickaA.AwsiukM.SadowskiJ.UndasA. (2014). Genetic determinants of acenocoumarol and warfarin maintenance dose requirements in Slavic population: a potential role of CYP4F2 and GGCX polymorphisms. Thromb. Res.134 (3), 604ā609. doi:Ā 10.1016/j.thromres.2014.06.022
Summary
Keywords
acenocoumarol, coumarins, algorithm, pharmacogenetics, pharmacogenomics, anticoagulation
Citation
Roco A, Nieto E, SuĆ”rez M, Rojo M, Bertoglia MP, Verón G, Tamayo F, Arredondo A, Cruz D, MuƱoz J, Bravo G, Salas P, MejĆas F, Godoy G, VĆ©liz P and QuiƱones LA (2020) A Pharmacogenetically Guided Acenocoumarol Dosing Algorithm for Chilean Patients: A Discovery Cohort Study. Front. Pharmacol. 11:325. doi: 10.3389/fphar.2020.00325
Received
04 December 2019
Accepted
05 March 2020
Published
06 April 2020
Volume
11 - 2020
Edited by
Nathalie K. Zgheib, American University of Beirut, Lebanon
Reviewed by
Andrew A. Somogyi, University of Adelaide, Australia; Ming Ta Michael Lee, Geisinger Health System, United States
Updates
Copyright
Ā© 2020 Roco, Nieto, SuĆ”rez, Rojo, Bertoglia, Verón, Tamayo, Arredondo, Cruz, MuƱoz, Bravo, Salas, MejĆas, Godoy, VĆ©liz and QuiƱones.
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: Luis Abel QuiƱones, lquinone@uchile.cl
This article was submitted to Pharmacogenetics and Pharmacogenomics, a section of the journal Frontiers in Pharmacology
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.