Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 06 April 2020
Sec. Pharmacogenetics and Pharmacogenomics
This article is part of the Research Topic Pharmacogenetics Research and Clinical Applications: An International Landscape of the Accomplishments, Challenges, and Opportunities View all 14 articles

A Pharmacogenetically Guided Acenocoumarol Dosing Algorithm for Chilean Patients: A Discovery Cohort Study

Angela Roco,,Angela Roco1,2,3Elena NietoElena Nieto4Marcelo SurezMarcelo Suárez1Mario Rojo,Mario Rojo1,5Maria Paz BertogliaMaria Paz Bertoglia6Gabriel VernGabriel Verón1Francisca TamayoFrancisca Tamayo1Annabella ArredondoAnnabella Arredondo7Daniela CruzDaniela Cruz4Jessica MuozJessica Muñoz4Gabriela BravoGabriela Bravo8Patricio SalasPatricio Salas9Fanny MejíasFanny Mejías10Gerald GodoyGerald Godoy10Paulo VlizPaulo Véliz10Luis Abel Quiones,*Luis Abel Quiñones1,5*
  • 1Laboratory of Chemical Carcinogenesis and Pharmacogenetics, Department of Basic and Clinical Oncology, Faculty of Medicine, University of Chile, Santiago, Chile
  • 2Escuela de Bioquímica Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
  • 3Western Metropolitan Health Service, Santiago, Chile
  • 4San Juan de Dios Hospital, Santiago, Chile
  • 5Latin American Network for Implementation and Validation of Clinical Pharmacogenomics Guidelines (RELIVAF-CYTED), Madrid, Spain
  • 6Institute of Population Health, University of Chile, Santiago, Chile
  • 7Instituto de Salud Pública, Universidad Andrés Bello, Santiago, Chile
  • 8Curacaví Hospital, Curacaví, Chile
  • 9Dr. Salvador Allende G. Reference Health Center, Santiago, Chile
  • 10San José de Melipilla Hospital, Melipilla, Chile

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
www.frontiersin.org

Table 1 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
www.frontiersin.org

Figure 1 Patient recruitment and selection flowchart.

TABLE 2
www.frontiersin.org

Table 2 Characteristics of the study patients.

Genotype Distribution in the Study Population

The analysis of the HWE showed that only CYP2C9*3 (rs1057910) is 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 not 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
www.frontiersin.org

Table 3 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
www.frontiersin.org

Table 4 Linear regression including genetic factors and clinical variables in a single model.

Log WTD = 3.081 + (0.167 × men) - (age × 0.081) - (initial INR × 0.55) + (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
www.frontiersin.org

Table 5 Comparison between the present study and the vitamin K antagonists (VKA) dosage algorithms published for other populations.

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.

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. The patients/participants provided their written informed consent to participate in this study.

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.

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.

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.

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.

References

Acuña, M., Llop, E., Rothhmmer, F. (2000). Genetic composition of Chilean population: rural communities of Elqui, Limari and Choapa valleys. Rev. Med. Chil. 128, 593–600.

PubMed Abstract | Google Scholar

Ageno, W., Gallus, A. S., Wittkowsky, A., Crowther, M., Hylek, E. M., Palareti, G. (2012). Oral anticoagulant therapy: Antithrombotic Therapy and Prevention of Thrombosis, 9th ed: American College of Chest Physicians Evidence-Based Clinical Practice Guidelines. Chest 141, e44S–e488. doi: 10.1378/chest.11-2292

PubMed Abstract | CrossRef Full Text | Google Scholar

Ansell, J., Hirsh, J., Hylek, E., Jacobson, A., Crowther, M., Palareti, G. (2008). Pharmacology and management of the vitamin K antagonists: american college of chest physicians evidence-based clinical practice guidelines (8th edition). Chest 133 (6 Suppl.), 160S–198S. doi: 10.1378/chest.08-0670

PubMed Abstract | CrossRef Full Text | Google Scholar

Borobia, A. M., Lubomirov, R., Ramírez, E., Lorenzo, A., Campos, A., Muñoz-Romo, R., et al. (2012). An Acenocoumarol Dosing Algorithm Using Clinical and Pharmacogenetic Data in Spanish Patients with Thromboembolic Disease. PloS One 7 (7), e41360. doi: 10.1371/journal.pone.0041360

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosch, J., Eikelboom, J. W., Connolly, S. J., Bruns, N. C., Lanius, V., Yuan, F., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Botton, M. R., Bandinelli, E., Rohde, L. E., Amon, L. C., Hutz, M. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Caldwell, M. D., Awad, T., Johnson, J. A., Gage, B. F., Falkowski, M., Gardina, P., et al. (2008). CYP4F2 genetic variant alters required warfarin dose. Blood 111 (8), 4106–4112. doi: 10.1182/blood-2007-11-122010

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerezo-Manchado, J. J., Rosafalco, M., Antón, A. I., Pérez-Andreu, V., Garcia-Barberá, N., Martinez, A. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Corbalán, R., Conejeros, C., Rey, C., Stockins, B., Eggers, G., Astudillo, C., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

De Oliveira Almeida, V. C., De Souza Ferreira, A. C., Ribeiro, D. D., Gomes Borges, K. B., Salles Moura Fernandes, A. P., Brunialti Godard, A. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Esteve-Pastor, M. A., Rivera-Caravaca, J. M., Roldán-Rabadán, I., Roldán, V., Muñiz, J., Raña-Míguez, P., 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 Europace 20 (9), 1435–1441. doi: 10.1093/europace/eux314

CrossRef Full Text | Google Scholar

Fitzmaurice, D. A., Accetta, G., Haas, S., Kayani, G., Lucas Luciardi, H., Misselwitz, F., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Galvez, J. M., Restrepo, C. M., Contreras, N. C. (2018). Creating and validating a warfarin pharmacogenetic dosing algorithm for Colombian patients. Pharmacogenomics pers. Med. 11, 169–178. doi: 10.2147/PGPM.S170515

CrossRef Full Text | Google Scholar

Guyatt, G. H., Akl, E. A., Crowther, M., Gutterman, D. D., Schuünemann, H. J. (2012). Executive summary: Antithrombotic therapy and prevention of thrombosis, 9th ed: American College of Chest Physicians evidence-based clinical practice guidelines. Chest 141 (2 Suppl.), 7–47. doi: 10.1378/chest.1412S3

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, S., Ten Cate, H., Accetta, G., Angchaisuksiri, P., Bassand, J. P., Camm, A. 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 One 11 (10), e0164076. doi: 10.1371/journal.pone.0164076

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J. A., Gong, L., Whirl-Carrillo, M., Gage, B. F., Scott, S. A., Stein, C. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J., Caudle, K., Gong, L., Whirl-Carrillo, M., Stein, C. M., Scott, S. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Kocael, A., Eronat, A. P., Tüzüner, M. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Krishna-Kumar, D., Shewade, D. G., Loriot, M. A., Beaune, P., Sai Chandran, B. V., Balachander, J., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Leschziner, G. D., Andrew, T., Pirmohamed, M., Johnson, M. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcatto, L. R., Sacilotto, L., Tavares, L. C., Facin, M., Olivetti, N., Strunz, C. 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.

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Nieto, E., Suarez, M., Roco, A., Rubilar, J. C., Tamayo, F., Rojo, M., et al. (2019). Anticoagulation Management With Coumarinic Drugs in Chilean Patients. Clin. Appl. Thrombosis/Hemostasis. 25, 1–7. doi: 10.1177/1076029619834342

CrossRef Full Text | Google Scholar

Pengo, V., Pegoraro, C., Cucchini, U., Iliceto, S. (2006). Worldwide management of oral anticoagulant therapy: The ISAM study. J. Thromb. Thrombolysis 21 (1), 73–77. doi: 10.1007/s11239-006-5580-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirmohamed, M., Burnside, G., Eriksson, N., Jorgensen, A. L., Toh, C. H., Nicholson, T., et al. (2013). A randomized trial of genotype-guided dosing of warfarin. N Engl. J. Med. 369, 2294–2303. doi: 10.1056/NEJMoa1311386

PubMed Abstract | CrossRef Full Text | Google Scholar

Quiñones, L., Roco, A., Miranda, C. (2015). Farmacogenómica: Aplicaciones Cardiovasculares. Rev. Médica Clínica las Condes 196 (2), 198–209. doi: 10.1016/j.rmclc.2015.04.008

CrossRef Full Text | Google Scholar

Ramos, A. S., Seip, R. L., Rivera-Miranda, G., Felici-Giovanini, M. E., Garcia-Berdecia, R., Alejandro-Cowan, Y., et al. (2012). Development of a pharmacogenetic-guided warfarin dosing algorithm for Puerto Rican patients. Pharmacogenomics 13 (16), 1937–1950. doi: 10.2217/pgs.12.171

PubMed Abstract | CrossRef Full Text | Google Scholar

Rathore, S., Agarwal, S. K., Pande, S., Singh, S. K., Mittal, T., Mittal, B. (2012). Therapeutic dosing of acenocoumarol: Proposal of a Population Specific Pharmacogenetic Dosing Algorithm and Its Validation in North Indians. PloS One 7 (5), e37844. doi: 10.1371/journal.pone.0037844

PubMed Abstract | CrossRef Full Text | Google Scholar

Riley, R. S., Rowe, D., Fisher, L. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, P. C., Marcatto, L. R., Duarte, N. E., Gadi Soares, R. A., Cassaro Strunz, C. M., Scanavacca, M., et al. (2015). Development of a pharmacogenetic based warfarin dosing algorithm and its performance in Brazilian patients: highlighting the importance of population-specific calibration. Pharmacogenomics 16 (8), 865–876. doi: 10.2217/pgs.15.48

PubMed Abstract | CrossRef Full Text | Google Scholar

Shendre, A., Brown, T. M., Liu, N. (2016). Race-Specific Influence of CYP4F2 on Dose and Risk of Hemorrhage Among Warfarin Users. Pharmacotherapy 36 (3), 263–272. doi: 10.1002/phar.1717

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Wu, Z., Li, S. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Tong, H. Y., Dávila-Fajardo, C. L., Borobia, A. M., Martínez-González, L. J., Lubomirov, R., Perea-León, L. 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 2016 11 (3), e0150456. doi: 10.1371/journal.pone.0150456. eCollection 2016.

CrossRef Full Text | Google Scholar

van Schie, R. M., Wessels, J. A., le Cessie, S., de Boer, A., Schalekamp, T., van der Meer, F. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Verde, Z., Ruiz, J. R., Santiago, C., Valle, B., Bandrés, F., Calvo, E., 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

CrossRef Full Text | Google Scholar

Verhoef, T. I., Ragia, G., de Boer, A., Barallon, R., Kolovou, G., Kolovou, V., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Visser, L. E., Trienekens, P. H., De Smet, P. A., Vulto, A. G., Hofman, A., van Duijn, C. M., et al. (2005). Patients with an APOE epsilon4 allele require lower doses of coumarin anticoagulants. Pharmacogenet. Genomics 15 (2), 69–74. doi: 10.1097/01213011-200502000-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Wu, A. H., Wang, P., Smith, A., Haller, C., Drake, K., Linder, M., et al. (2008). Dosing algorithm for warfarin using CYP2C9 and VKORC1 genotyping from a multi-ethnic population: comparison with other equations. Pharmacogenomics 9 (2), 169–178. doi: 10.2217/14622416.9.2.169

PubMed Abstract | CrossRef Full Text | Google Scholar

Wypasek, E., Branicka, A., Awsiuk, M., Sadowski, J., Undas, A. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

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

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, bHF1aW5vbmVAbWVkLnVjaGlsZS5jbA==

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.