- 1Primary Care Pharmacy, Aljarafe-Sevilla Norte District. Andalusian Health Service, Seville, Spain
- 2Coordination of Management and Evaluation, Andalusian Health Service., Seville, Spain
- 3Andalusian Public Foundation for the Management of Health Research in Seville (FISEVI), Seville, Spain
- 4Primary Care Pharmacy, Seville District. Andalusian Health Service, Seville, Spain
- 5Andalusian School of Public Health, Granada, Spain
- 6Andalusian Health Technology Assessment Area, Progress and Health Foundation, Seville, Spain
- 7Center for the Transfusion of Tissues and Cells, Seville, Spain
- 8Research Unit, Andalusian Health Service, Seville, Spain
Objective: The aim of this study was to evaluate the effectiveness and safety of direct oral anticoagulants (DOACs: dabigatran, rivaroxaban, apixaban and edoxaban) and warfarin versus acenocoumarol in patients with atrial fibrillation under real-world clinical practice conditions.
Methods: This was a retrospective, real-world data-based study. The data source was the Andalusian Population Health Database. The study covered the period from January 2012 to December 2020. Effectiveness outcomes were defined as the identification of a first occurrence of ischaemic or bleeding events, or all-cause mortality. The statistical analysis included crude incidence analysis, survival models: Kaplan-Meier curves, propensity score matched pairs analysis, Fine-Gray model, and Cox regression analysis adjusted for possible confounding.
Results: A total of 150,949 patients were included. The mean age of the cohort was 74 years (48.2% female). The mean follow-up time was 3.3 years. The combined effectiveness endpoint of ischaemic events (transient ischaemic attack, systemic embolism, pulmonary embolism, or ischaemic stroke) showed the following results compared to acenocoumarol: warfarin (RR:1.06; 95%CI 0.93–1.22); dabigatran (RR:1.17; 95%CI 1.02–1.33); rivaroxaban (RR:1.15; 95%CI 1.05–1.26); apixaban (RR: 0.96; 95%CI 0.87–1.07) and edoxaban (RR: 1.10; 95%CI 0.79–1.51). Compared to acenocoumarol, the risk of all-cause mortality was lower for dabigatran, rivaroxaban and apixaban (RR:0.77; 95%CI 0.72–0.82; RR:0.79; 95%CI 0.76–0.83; RR:0.85; 95%CI 0.81–0.89, respectively) and higher for warfarin (RR:1.12; 95%CI 1.05–1.20). An increased risk of gastrointestinal bleeding was observed with dabigatran (RR:1.36; 95%CI 1.09–1.70) and a lower risk with rivaroxaban (RR:0.84; 95%CI 0.72–0.98). All 4 DOACs showed a lower risk of intracranial bleeding compared to acenocoumarol. Warfarin carried a higher risk of both gastrointestinal bleeding (RR:1.64; 95%CI 1.31–2.06) and intracranial bleeding (RR:1.61; 95%CI 1.22–2.13) compared to acenocoumarol. An unadjusted analysis of matched groups in a multivariate Cox regression analysis yielded similar results for combined effectiveness and safety outcomes compared to acenocoumarol.
Conclusion: Although DOACs were clearly associated with a lower risk of intracranial bleeding compared to acenocoumarol, our data did not reveal a significant reduction in thromboembolic events. Warfarin was found to be both less effective and less safe than acenocoumarol.
Introduction
Preventing ischaemic or thromboembolic events is one of the primary goals of anticoagulation therapy for patients with atrial fibrillation (AF). Before direct-acting oral anticoagulants (DOACs) became available, vitamin K antagonists (VKAs) such as warfarin or acenocoumarol were the preferred first-line agents for preventing ischaemic events in patients with AF (January et al., 2019; Hindricks et al., 2021).
Clinical trials investigating the use of DOACs in the treatment of AF were conducted with warfarin as the comparator. These trials showed that DOACs had similar efficacy but posed a lower risk of intracranial haemorrhage when compared with warfarin (Connolly et al., 2009; Granger et al., 2011; Patel et al., 2011; Giugliano et al., 2013). In light of the results, and given the fact that DOACs offer simpler dosing, fewer drug and food interactions and do not require anticoagulation monitoring (of the International Normalized Ratio, INR), some scientific societies and clinical practice guidelines now recommend DOACs as the first-line agents in patients with AF, replacing the use of vitamin K antagonists (January et al., 2019; Hindricks et al., 2021).
We decided to conduct this study for several reasons. Firstly, no clinical trials have been identified that compare different DOACs head-to-head. Meanwhile, observational studies have found differences in the outcomes of each DOAC compared to warfarin, as well as in the outcomes of clinical trials (Nielsen et al., 2017; Tepper et al., 2018; Chan, Lee, Chao, et al., 2019; Chan, Lee, See, et al., 2019; Ramagopalan et al., 2019; Xue and Zhang, 2019; Douros et al., 2019; Xue et al., 2020; Domek et al., 2020; Lip et al., 2021; Durand et al., 2021; Ingason et al., 2021; Kim et al., 2021; Li et al., 2021; Mamas et al., 2022; Jaksa et al., 2022; Lau et al., 2022).
Furthermore, while warfarin is the reference VKA drug in the majority of studies conducted in AF patients, acenocoumarol is the drug of choice in certain populations, such as Spain and the Netherlands (Rosa et al., 2018; Rodríguez-Bernal et al., 2021). One possible reason why the results obtained with warfarin in clinical trials and observational studies are not generalisable to acenocoumarol is the difference in their half-lives. When extrapolating data to real-world practice, compliance should also be considered. Research findings from one observational study have indicated that patients with AF who were prescribed DOACs were three times more likely to be non-adherent than those prescribed VKAs (Rodriguez-Bernal et al., 2018). Adherence issues, as well as other contributing factors, such as comorbidities, concomitant medications, age or sex, may affect the effectiveness of oral anticoagulants used in actual clinical practice. Finally, no clinical trials or observational studies have compared the efficacy, effectiveness or safety of edoxaban, which is the most recently marketed DOAC, with acenocoumarol.
The main aim of this study was to evaluate the effectiveness (ischaemic/transient stroke, systemic/pulmonary embolism and mortality) and safety (gastrointestinal and intracranial bleeding) of DOACs (dabigatran, rivaroxaban, apixaban and edoxaban) and warfarin versus acenocoumarol in a large population cohort under real-world clinical practice conditions with 8 years of follow-up.
Materials and methods
Study design
This retrospective, real-world, data-based study of patients with AF aimed to analyse the effectiveness and safety of oral anticoagulant (OAC) therapy (warfarin, dabigatran, rivaroxaban, apixaban and edoxaban) compared to acenocoumarol (Montero-Balosa et al., 2023). In an attempt to simulate a randomised experiment (quasi-experimental study), all outcomes were analysed using intention-to-treat (ITT) and propensity score matching. The study period ran from January 2012 to December 2020.
Population and setting
The study was carried out in the Autonomous Community of Andalusia (Spain). The Andalusian Health Service (AHS) provides approximately 8.5 million people in this region with free universal healthcare at the point of use.
Patients in the study were seen in consultations at any level of healthcare (primary care or hospital) within the AHS and met the following inclusion criteria: they were over 40 years of age, had received a diagnosis of atrial fibrillation (AF) based on ICD-9 or ICD-10 (International Classification of Diseases) codes and were initiating treatment with oral anticoagulants (OACs) in the context of routine clinical practice (naïve population). Only new OAC users were included. This meant those who had previously taken OACs were included as long as provided that they had stopped taking them in the 12 months prior to their inclusion in the study.
Patients were excluded if they had severe mitral stenosis, a diagnosis of valvular heart disease, or were undergoing aortic and/or mitral valve procedures, as defined by the ICD-9 or ICD-10 codes.
A more detailed description of the ICD-9 and ICD-10 codes can be found in the published SIESTA-A study protocol (Montero-Balosa et al., 2023).
Data sources
The data source for this study was Andalusia’s Population Health Database (BPS, Base Poblacional de Salud, Spanish acronym). This population-based health information system collects clinical data of every individual receiving healthcare in the AHS, as well as information on their use of health resources (Muñoyerro-Muñiz et al., 2020). This information is stored in different databases, which are linked by each patient’s unique identification number in the AHS. The information contained in this population health database includes patient affiliations, socio-demographics, diagnosed health problems, cardiovascular risk factors and lifestyles, utilisation of resources in outpatient consultations, primary care, hospital care and emergency care, hospital admissions, mortality, and prescription and pharmacy dispensing data in primary and hospital care (Montero-Balosa et al., 2023).
Outcomes
The effectiveness outcomes (dependent variables) were defined as the incidence of a new first event from the following list: transient ischaemic attack (TIA), systemic embolism, pulmonary embolism, ischaemic stroke, or all-cause mortality. These were considered a combined effectiveness endpoint (Montero-Balosa et al., 2023).
The safety outcomes were defined as the incidence of a new first event of major bleeding leading to hospital admission (gastrointestinal or intracranial bleeding). These were considered a combined safety endpoint (Montero-Balosa et al., 2023). Post-traumatic bleeding was excluded.
All effectiveness and safety events were associated with the date of onset and identified by the relevant ICD-9 or ICD-10 code when the patient was treated in hospital. Diagnoses of TIA with a recorded date were collected from both hospitals and primary care centres.
The index date (i.e., the date of the first treatment) was defined as the date on which a prescribed anticoagulant treatment was dispensed from the pharmacy, provided that no other anticoagulant had been prescribed in the previous 12 months.
Follow-up time was defined as the number of days between the index date of treatment initiation (i.e., pharmacy dispensing) and the date of the first new diagnosis, death or the end of the study period (31 December 2020).
Independent variables
The independent variables included in the multivariate survival analysis were the type of OAC therapy, sociodemographic data (age, sex, place of residence, approximate income based on pharmaceutical co-payment status), comorbidities, medication use, CHADS-VASC (thromboembolic risk assessment) and HAS-BLED (bleeding risk assessment) scores. All of these were assessed at study entry. Frequency of contact with the AHS (primary care, hospital, home visits, and emergency department) was also considered (Montero-Balosa et al., 2023).
Comorbidities were defined as an active diagnosis (ICD-9 or ICD-10) recorded in one of the population health databases (emergency, hospital or primary care) in the 12 months prior to the date of treatment initiation. Medications that increase or reduce bleeding risk were included in the HAS-BLED calculation (Montero-Balosa et al., 2023). All predictor variables were analysed for their association with events, based on the information available in the 12 months prior to the index date.
Statistical analysis
Descriptive statistics were calculated and all outcomes were analysed using the intention-to-treat (ITT) method. Patients were categorised according to the initial OAC treatment started during the study period, provided that they had not taken any other OAC for at least 12 months prior to the study (“washout period”) and regardless of possible changes in OAC therapy during the follow-up period.
When measuring the effectiveness and safety outcomes, only the first recorded instance of each type of event included in the study was considered. If a primary outcome variable was identified in two or more study databases, the first recorded date was used.
Quantitative variables were summarised using means and standard deviation. Non-quantitative variables were summarised using frequency tables and percentages. Point estimates and 95% confidence intervals were obtained for the various statistical analyses. Incidence rates were calculated by dividing the number of events by the corresponding person-time of follow-up.
Sensitivity analysis
In addition to the main analysis, sensitivity analyses were performed to consider the effect of sex and the exclusion of mortality from the combined effectiveness variable. The TIAs were also studied, selecting only those with a confirmed hospital diagnosis. The Fine-Gray model was used to analyse the competing risk of death (Austin and Fine, 2017).
Phases of the analysis
The analysis comprised three phases. The first phase consisted of a descriptive analysis of the total number of individuals included in the study, and of the patients in each of the 6 cohorts, which were grouped according to the type of initial anticoagulant taken. Next, the incidence of effectiveness and safety events over the 8 years of follow-up was analysed.
In the second phase, Propensity Score Matching (PSM) was performed. The nearest neighbor matching approach was used to compare patients taking the 4 DOACs or warfarin to patients taking acenocoumarol (the reference drug).
Propensity Score (PS) values (between 0 and 1) indicated the a priori probability that a patient would receive treatment with acenocoumarol compared to the OAC under comparison. PS scores were calculated using logistic regression models to predict the binary dependent variable. Initially, the following independent variables were included: age, sex, province, pharmaceutical co-payment status, use of healthcare resources by the patient, comorbidities, use of concomitant medications at study entry associated with an increased risk of bleeding (acetylsalicylic acid, other antiplatelet agents, heparins, glucocorticoids, non-steroidal anti-inflammatory drugs, selective serotonin reuptake inhibitors, macrolides) as well as those associated with a reduced risk of bleeding (proton pump inhibitors, H2-receptor antagonists, antihaemorrhagic drugs) and the CHADS-VASC and HAS-BLED scores (Montero-Balosa et al., 2023).
Mean differences for quantitative variables and differences of distribution for categorical variables after PSM were then compared with the differences observed in the original 6 cohorts. The “matchit” function in the Matchit package (R version 4.1.0) was used for the matching analysis to perform pairing, subset selection and subclassification to ensure covariate balance between treatment and control groups.
After matching, the incidence rate of events (person-time in follow-up) was analysed in relation to the effectiveness and safety of the treatment over the 8-year follow-up period. Only the cases selected for each pair of oral anticoagulants were compared (acenocoumarol was used as the reference). Bivariate survival analysis (for time-censored data) was performed using Kaplan-Meier curves to plot the frequency of different events or diseases occurring over time. Survival probabilities were calculated by generating statistics and plotting survival functions for each study group. The log-rank test was then used to test for equality of survival time distributions between the different groups.
Multivariate survival analysis was performed on the combined effectiveness and safety outcomes. Cox regression was used to create time-to-event models using both categorical and continuous predictor variables (covariates) as hypothesised predictors. In addition to Cox regression, a Fine-Gray analysis was conducted to consider the competing risk of death. First, we tested whether the OAC versus acenocoumarol survival curves crossed in the matched groups and whether the proportional hazards assumption was met visually over time in the Kaplan-Meier curves. The risk of outcomes (effectiveness and safety) in DOAC or warfarin users versus acenocoumarol users (reference) was determined by calculating specific instantaneous risk ratios (hazard ratio, HR) for each comparison using multivariate Cox regression models.
A matched-pair analysis was performed for the period 2013–2020. For the acenocoumarol-apixaban and acenocoumarol-edoxaban pairs, which were launched after 2013, the analysis periods were 2014–2020 and 2017–2020, respectively.
The information was extracted from the population health databases and transferred into a single file. This file was then analysed using R (version 4.1.0), and Stata (version 14 for Windows) statistical software. Patients were anonymised for evaluation by the study investigators.
The study was reported in accordance with the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement for cohort studies (Gharaibeh et al., 2014).
Registration details
The study protocol was submitted to the Spanish Agency of Medicines and Medical Devices (AEMPS, Spanish acronym) and classified as an observational study (AEMPS reference number: 0004–2022-OBS; 12 January 2022).
Results
After applying the inclusion and exclusion criteria, a total of 150,949 patients were included in the six study cohorts (Figure 1).
The most common comorbidities were hypertension (77.7%), diabetes (35.8%), cardiac arrhythmias (25.4%), heart failure (21.6%) and chronic obstructive pulmonary disease (COPD, 23.7%). The mean age of the population was 74 years (SD:10.2), and 48.2% of patients were female (n = 72,817) (Table 1).
The mean follow-up time was 3.3 years: 3.8 years for acenocoumarol, 3.6 years for warfarin, 3.7 years for dabigatran, 3.3 years for rivaroxaban, 2.5 years for apixaban and 1.6 years for edoxaban.
At study entry, patients treated with dabigatran were younger (with a higher proportion in the 40–59 age group). They also had the lowest mean CHADS-VASC and HAS-BLED scores, as well as a lower prevalence of hypertension, renal insufficiency, diabetes and heart failure. Patients on apixaban had a higher prevalence of previous stroke, dementia, previous gastrointestinal bleeding and previous intracranial bleeding. Those on edoxaban showed a lower prevalence of myocardial infarction, arrhythmia and previous stroke. Patients on warfarin had a higher prevalence of diabetes, angina, myocardial infarction, arrhythmia, heart failure and peripheral arterial disease. Finally, patients treated with acenocoumarol had a higher prevalence of venous thrombosis, COPD and renal failure compared to those treated with the other five anticoagulants (Table 1).
At the beginning of the study, patients taking acenocoumarol or warfarin were more likely to be taking heparins and antithrombotics (18.6% and 17.7%, respectively) than those taking DOACs (range: 4.8%–8.8%) and were also more likely to be taking antiplatelet agents (27.1% and 31.6%, respectively) than those taking DOACs (range: 21.9%–26.6%). Patients receiving warfarin treatment had a higher consumption of H2 blockers and gastroprotective drugs (9.3%) than those on DOACs (range: 5.0%–6.6%) (Table 2).

Table 2. Utilisation of medication and healthcare resources during follow-up by province, pharmaceutical co-payment status and anticoagulant.
In terms of socioeconomic status, treatment with a DOAC was more common among patients who contributed 40%–50% of the drug’s cost (i.e., younger, employed patients). In contrast, retired patients who contributed less than 10% of the cost were more likely to be treated with a VKA (Table 2).
During the follow-up period, utilisation of health services (outpatient appointments, home visits and primary care emergencies) was higher among patients taking acenocoumarol or warfarin than among those treated with DOACs (Table 2).
After PSM, the differences in means and categorical distributions of the various variables were examined and compared with those observed in the original six cohorts. The distances between variables were reduced in the matched groups, indicating greater similarity between them, thereby allowing for a more accurate comparison of the treatment effect and increasing the validity of the results (Supplementary Table S1). The selection of probability ranges and overlaps for matching patients undergoing treatment with each pair of anticoagulants is shown in Supplementary Figure S1.
Effectiveness results
The first analysis of incidence rates, using number of events per 1,000 person-years of follow-up within matched groups (PSM), showed that dabigatran (RR:0.84; 95%CI 0.79–0.89), rivaroxaban (RR:0.85; 95%CI 0.82–0.88), and apixaban (RR:0.88; 95%CI 0.84–0.91) were more effective on the combined effectiveness endpoint compared to acenocoumarol, while warfarin was less effective (RR:1.11; 95%CI 1.05–1.18) (Table 3).

Table 3. Matched group comparisons of each anticoagulant versus acenocoumarol (propensity score matching): Relative risk and 95% confidence intervals for combined effectiveness, combined safety, and specific event types.
Following the second analysis, which excluded all-cause mortality from the combined effectiveness outcome, dabigatran and rivaroxaban were found to be associated with an increased risk of ischaemic events (RR:1.17; 95%CI 1.02–1.33; RR:1.15; 95%CI 1.05–1.26, respectively) (Table 4). After sensitivity analysis (using the Fine-Gray model) to account for the competing risk of death, the results were consistent for all the anticoagulants under study (Table 5).

Table 4. Sensitivity analysis of matched groups of anticoagulants versus acenocoumarol (propensity score matching). Relative risk and 95% confidence intervals for all ischaemic events and transient ischaemic attack for the total population and according to sex.

Table 5. Sensitivity analysis of matched group of anticoagulants versus acenocoumarol (propensity score matching): Results of the Fine-Gray model for competing risk of death, with subdistribution hazard ratios and 95% confidence intervals for ischaemic events and combined safety events.
An analysis of specfic event types in this study revealed an increased risk of TIA with rivaroxaban (RR: 1.18; 95%CI 1.04–1.34) and apixaban (RR: 1.17; 95%CI 1.01–1.35), and a lower risk of ischaemic stroke with apixaban (RR: 0.86; 95%CI 0.75–0.98) compared to acenocoumarol. Warfarin was less effective than acenocoumarol in preventing systemic embolism (RR: 2.05; 95%CI 1.16–3.76) (Table 3).
The risk of death was lower for dabigatran, rivaroxaban and apixaban (RR: 0.77; 95%CI 0.72–0.82; RR: 0.79; 95%CI 0.76–0.83; RR: 0.85; 95%CI 0.81–0.89, respectively) and higher for warfarin compared to acenocoumarol (RR: 1.12; 95%CI 1.05–1.20) (Table 3; Figure 2).

Figure 2. Effectiveness outcomes. Relative risk and 95% confidence intervals for combined effectiveness and specific event types. Effectiveness (combined): transient ischaemic attack, systemic embolism, pulmonary embolism, ischaemic stroke or all-cause mortality.
Safety results
The analysis of incidence rates by matched cohorts (PSM) showed a lower incidence of combined safety events with rivaroxaban (RR: 0.73; 95%CI 0.64–0.82) and apixaban (RR: 0.85; 95%CI 0.74–0.97) compared to acenocoumarol. An increased risk of gastrointestinal bleeding was observed with dabigatran (RR: 1.36; 95%CI 1.09–1.70) and a lower risk with rivaroxaban (RR: 0.84; 95%CI 0.72–0.98). The risk of intracranial bleeding was lower for all 4 DOACs compared to acenocoumarol. Warfarin had a higher risk of combined safety events (RR: 1.63; 95%CI 1.37–1.95), gastrointestinal bleeding (RR: 1.64; 95%CI 1.31–2.06) and intracranial bleeding (RR: 1.61; 95%CI 1.22–2.13) compared to acenocoumarol (Table 3; Figure 3).

Figure 3. Safety outcomes. Relative risk and 95% confidence intervals for combined safety and specific event types. Safety (combined): gastrointestinal or intracranial bleeding.
Figures 4–8 show the cumulative incidence (Kaplan-Meier curves) of all-cause mortality, ischaemic stroke and major bleeding (intracranial or gastrointestinal) for each pair of anticoagulants.

Figure 4. Acenocoumarol vs. warfarin. Crude cumulative incidence or failure curves of (A) “all-cause mortality”, (B) “ischaemic stroke”, and (C) “major bleeding” (intracranial or gastrointestinal bleeding) according to initiated treatment.

Figure 5. Acenocoumarol vs. dabigatran. Crude cumulative incidence or failure curves of (A) “all-cause mortality”, (B) “ischaemic stroke”, and (C) “major bleeding” (intracranial or gastrointestinal bleeding) according to initiated treatment.

Figure 6. Acenocoumarol vs. rivaroxaban. Crude cumulative incidence or failure curves of (A) “all-cause mortality”, (B) “ischaemic stroke”, and (C) “major bleeding” (intracranial or gastrointestinal bleeding) according to initiated treatment.

Figure 7. Acenocoumarol vs. apixaban. Crude cumulative incidence or failure curves of (A) “all-cause mortality”, (B) “ischaemic stroke”, and (C) “major bleeding” (intracranial or gastrointestinal bleeding) according to initiated treatment.

Figure 8. Acenocoumarol vs. edoxaban. Crude cumulative incidence or failure curves of (A) “all-cause mortality”, (B) “ischaemic stroke”, and (C) “major bleeding” (intracranial or gastrointestinal bleeding) according to initiated treatment.
A comparison of incidence rates for combined effectiveness and safety by anticoagulant pair is shown in Supplementary Tables S2-S6. Incidence rates for specific event types are shown in Supplementary Tables S7-S11.
Analysis of other confounding factors
The results obtained using the described methodology (calculation of relative risks, prior PSM, Table 3) are robust, as demonstrated by their concordance with those obtained using multivariate modelling (see Supplementary Table S12) and hazard ratios (HR). This was observed for the combined effectiveness events: RR = 1.11 vs. HR = 1.05 for warfarin; RR = 0.84 vs. HR = 0.82 for dabigatran; RR = 0.85 vs. HR = 0.91 for rivaroxaban; and RR = 0.88 vs. HR = 0.96 for apixaban.
Similar results were obtained for the combined safety events: RR = 1.63 vs. HR = 1.64 for warfarin; RR = 0.88 vs. HR = 0.84 for dabigatran; RR = 0.73 vs. HR = 0.71 for rivaroxaban; and RR = 0.85 vs. HR = 0.87 for apixaban.
These data can be explained by controlling for confounding factors in the models, which had little influence on the outcome. Propensity score matching of groups prior to multivariate analysis reduced the distances between variables. Relative risks and hazard ratios are typically very similar when follow-up times are the same for the entire sample (Martínez-González and Alonso 2008) and in our study, follow-up times were largely homogeneous between pairs after matching cases on propensity score.
Sensitivity analysis
Consistency and possible differences in outcomes by sex were assessed by performing new analyses of incidence density rates and RR by propensity score-matched groups. The statistically significant reduction in the combined effectiveness risk of ischaemic events was maintained with dabigatran, rivaroxaban and apixaban in males, and with rivaroxaban and apixaban in females. An increased risk of ischaemic stroke was observed with dabigatran in females (Tables 6, 7).

Table 6. Study in Females: Matched groups of anticoagulants versus acenocoumarol (propensity score matching): Relative risk and 95% confidence intervals for combined effectiveness, combined safety and specific event types.

Table 7. Study in Males: Matched groups of anticoagulants versus acenocoumarol (propensity score matching): Relative risk and 95% confidence intervals for combined effectiveness, combined safety and specific event types.
The increased risk of TIA associated with rivaroxaban and apixaban remained higher in males. The results of the sensitivity analysis of TIAs with a confirmed hospital diagnosis were consistent with those of the main analysis. These results indicated a higher risk of TIA with apixaban in males (Table 4).
The risk of mortality was lower in both males and females when using dabigatran, rivaroxaban and apixaban (Tables 6, 7).
However, excluding all-cause mortality from the combined effectiveness variable in the sensitivity analysis revealed an increased risk of ischaemic events with dabigatran and rivaroxaban. These results were confirmed in female patients (Table 4).
The subdistribution hazard ratios from the Fine-Gray model of combined effectiveness (excluding mortality) showed the same trend in ischaemic event results when mortality was integrated as a competing risk in survival analysis (Table 5). The results of the combined safety variable also exhibited the same trend after the Fine-Gray model was applied (Table 5).
The statistically significant decreased risk of combined safety was maintained with rivaroxaban in both sexes. The reduced risk of intracranial bleeding was confirmed for dabigatran and rivaroxaban in both sexes (Tables 6, 7). In contrast, warfarin was less effective than acenocoumarol in males and showed a higher risk of major bleeding (intracranial and gastrointestinal bleeding) in both sexes.
Discussion
Our study was conducted on a large population sample (>150,000 patients) with long-term patient follow-up (real world data).). Using this sample, we compared patient health outcomes for acenocoumarol, (a VKA for which little data from clinical trials or observational studies is available) and warfarin (another VKA). We also estimated the effectiveness and safety of four direct oral anticoagulants (DOACs) compared to acenocoumarol.
Although a clear tendency towards reduced intracranial bleeding was observed with DOACs compared to acenocoumarol, our data did not demonstrate a significant decrease in thromboembolic events when acenocoumarol was used as the reference drug. Warfarin appeared to be less effective and safe than acenocoumarol.
As edoxaban was the most recently marketed DOAC (2017), the mean follow-up time was only 1.6 years and the number of events occurring was small (Supplementary Table S11). There were no statistically significant differences in effectiveness or safety between acenocoumarol and edoxaban. Consequently, firm conclusions cannot be drawn. The SIESTA research team has received funding to extend the follow-up period of the study by 3 years and it is hoped that this will allow us to draw clear conclusions about the effectiveness and safety of edoxaban compared to acenocoumarol.
One factor that lends robustness to our study is that the results for the total population are consistent with the trends observed when the data are stratified by sex. In the latter case, however, some of these results are not statistically significant due to smaller sample sizes when patients are disaggregated.
Different prescribing patterns were observed for oral anticoagulants. Among DOACs, dabigatran was more commonly used among younger patients, whereas apixaban was more frequently prescribed to older patients (over 70 years of age). This pattern of use is consistent with that reported by other investigators (Larsen et al., 2016; Durand et al., 2021; Rodríguez-Bernal et al., 2021). The greater use of dabigatran among younger patients may be related to the recommendation in the package leaflet to reduce the dose (from 150 mg to 110 mg every 12 h) for patients over 80 years of age.
Among patients taking VKAs, those receiving warfarin had a higher incidence of previous heart disease than those treated with acenocoumarol.
Statistical analysis based on PS matching enabled the creation of groups with variables with more similar means or categorical distributions between the different pairs of anticoagulants studied, which allowed us to assess health outcomes in everyday clinical practice using directly comparable groups. This strengthened the conclusions of our study (Ciminata et al., 2022).
DOACs vs. acenocoumarol: effectiveness
No clear reduction in thromboembolic risk was observed with the use of any DOAC compared to acenocoumarol. The results of the combined effectiveness analysis were significantly influenced by the high proportion of deaths among the total number of first effectiveness events recorded. This led to different outcomes depending on whether all-cause mortality was included in the combined outcome. When all-cause mortality was excluded from the combined effectiveness assessment, dabigatran and rivaroxaban were associated with an increased risk of ischaemic events (at least 2% and 5%, respectively, when considering the lower limit of the confidence interval).
We found no immediate explanation for the increased risk of TIA (at least 1%) alongside the reduced risk of ischaemic stroke (at least 2%) with apixaban. Rivaroxaban was also associated with an increased risk of TIA (at least 4%). This slight increase in the risk of TIA was unexpected, given that DOACs were thought to be as effective, if not more so, than VKAs (Rodríguez-Bernal et al., 2021). One possible explanation for these results may be that the TIA diagnoses were extracted not only from hospital databases but also from notes made by primary care physicians. In the absence of diagnostic confirmation, the incidence of TIA may have been overestimated. Many patients presenting with disorientation and focal neurological signs may have been diagnosed as possible TIA in primary care and referred to hospital for diagnostic confirmation. However, focal neurological signs, disorientation, or consciousness deficit may disappear without being confirmed by neurological imaging tests. This would suggest that the transient neurological condition was not thrombotic in origin. The results of our sensitivity analysis, which included only TIAs with confirmed hospital diagnosis, were consistent with those of the main analysis, and showed a small increase in the risk of TIA for rivaroxaban and apixaban (at least 1%).
These data are inconsistent with those of Rodríguez-Bernal et al. (2021), who found no significant overall difference in the incidence of TIA between the three DOACs evaluated (dabigatran, rivaroxaban, and apixaban) compared to acenocoumarol.
In our study, apixaban was the only DOAC that showed a reduced risk of ischaemic stroke compared to acenocoumarol. The other DOACs showed a trend towards an increased risk of ischaemic stroke, though this was not statistically significant. These results are consistent with those of Rodríguez-Bernal et al. (2021).
In a Spanish cohort of 2,178 patients, Anguita et al. (2020) observed a trend towards a lower incidence of stroke, systemic embolism, and all-cause mortality with DOACs compared to acenocoumarol, although the differences were not statistically significant. The study’s substantial data set allowed for statistically significant but small RR values close to 1.
In our study, all-cause death was the most common event recorded, with 29,475 cases corresponding to 19.5% of those under follow-up and 61.8% of total events recorded, which is consistent with the results obtained by Rodríguez-Bernal et al. (2021), in which mortality was the most frequently recorded event (accounting for around 70% of events). Three DOACs (dabigatran, rivaroxaban, and apixaban) showed a decreased risk of all-cause mortality compared to acenocoumarol. These data contradict those obtained by Rodríguez-Bernal et al. (2021) and Anguita et al. (2020) in similar populations, who found no difference in mortality for DOACs compared to acenocoumarol. The differences in our study may be due to the longer follow-up period (8 years) and the higher number of all-cause deaths. Furthermore, the statistical analysis used by Rodríguez-Bernal et al. (2021) was based on (inverse probability weighting (IPW), whereas our study used propensity score matching. Nevertheless, when other investigators (Austin, 2011; Ciminata et al., 2022) used these two analytical models to evaluate DOACs in patients with AF, they found that there were no statistically significant differences.
With regard to mortality (which accounted for almost one-fifth of the patients under follow-up), we emphasise the systematic and thorough collection of death data from the official statistical sources of the Autonomous Community of Andalusia. We do not know which factors had a greater impact on the risk of mortality in patients with AF. Although this study considered many demographic and health-related variables, it would have been impossible to record every factor or describe the complex interplay of causes implicated in such a wide-ranging issue as the probability of death. The use of medications was just one of the variables involved.
The paired analysis by sex revealed a mortality trend similar to that observed in the overall sample. Rodríguez-Bernal et al. (2021) also observed a reduced risk of mortality among males taking dabigatran compared to acenocoumarol, but not among females.
DOACs vs. acenocoumarol: safety
The comparative safety evaluation showed that the overall risk of bleeding events was significantly lower with rivaroxaban and apixaban than with acenocoumarol. All four DOACs evaluated against acenocoumarol showed a significant decrease in the incidence of intracranial bleeding. However, rivaroxaban was the only DOAC to significantly reduce the incidence of gastrointestinal bleeding, while dabigatran was the least safe DOAC, increasing the risk.
These results differ somewhat from those of Rodríguez-Bernal et al. (2021), who found no differences in gastrointestinal bleeding. However, they observed a reduction in intracranial bleeding with dabigatran and rivaroxaban compared to acenocoumarol. Anguita et al. (2020). Observed a trend towards a lower incidence of severe bleeding with DOACs compared to VKAs (primarily acenocoumarol), although it was not statistically significant.
Our results are partly consistent with those obtained by van den Ham et al. (2021), who conducted a pooled analysis of databases from four European countries and six Canadian provinces. They found that apixaban showed a lower risk of gastrointestinal bleeding, whereas dabigatran and rivaroxaban showed an increased risk compared to VKAs. They also confirmed a decreased risk of intracranial bleeding for all DOACs compared to VKAs.
With regard to the sensitivity analysis by sex, our results are consistent with those obtained in the primary analysis. The risk of intracranial bleeding is consistent with that obtained by Rodríguez-Bernal et al. (2021), who found a statistically significant reduction for dabigatran and rivaroxaban in females, and for dabigatran in males, compared to acenocoumarol.
As some investigators have pointed out (Ingason et al., 2021), the effectiveness and safety of DOACs in real-world practice is greatly influenced by differences in patient adherence. Although DOACs have similar half-lives, dabigatran and apixaban require twice-daily dosing. Durand et al. (2021) have suggested that this could explain their observation of an increase in ischaemic and bleeding events with rivaroxaban at the recommended once-daily dose, since this is associated with higher peak plasma concentrations that decline progressively after dosing. However, our results do not support this theory. While it could be argued that lower adherence to apixaban (which requires two daily doses) might explain the trend towards an increase in the risk of TIA observed in our study, we could find no justification for the same increase in TIA with rivaroxaban (with a once-daily dose) nor for the reduced risk of ischaemic stroke with apixaban. As we have not conducted a comparative study between different DOACs or evaluated patient adherence, we are unable to draw conclusions about the effect of the dosing regimen on health outcomes.
The RE-LY study (Connolly et al., 2009) revealed that the proportion of patients who discontinued dabigatran treatment due to gastrointestinal complications was three times higher than with warfarin. These results could also explain those obtained in our study, in which dabigatran was the only DOAC to show an increased risk of major gastrointestinal bleeding.
It is important to note that we did not compare the effectiveness and safety of DOACs with patients receiving treatment with VKA with good anticoagulation control (high time in therapeutic range, TTR), which may have overestimated the effect of DOACs in the study. Wallentin et al. (2010) investigated the primary and secondary outcomes of the RE-LY trial (Connolly et al., 2009) in relation to the mean TTR in each centre in warfarin patients. They concluded that dabigatran offered greater advantages for all vascular events, non-haemorrhagic events, and mortality at sites with poor INR control than at sites with good INR control in warfarin patients. Joosten et al. (2024) conducted a clinical trial in frail patients with AF (aged ≥75 years with a Groningen Frailty Indicator score ≥3) and concluded that switching from INR–guided VKA treatment to a DOAC in frail older patients with AF was associated with more bleeding complications compared with continuing VKA treatment, without reducing thromboembolic complications. Doni et al. (2024) conducted a systematic review of randomized controlled trials on DOACs to assess the safety of long-term intake of DOACs in older adults with AF. They found that elderly patients (aged 85 years and older) who use DOACs may be at an increased risk of major or clinically relevant bleeding events compared to those who use VKAs.
Based on these considerations, our study’s findings would support the recommendation in the 2024 European Society of Cardiology Guidelines (Van Gelder et al., 2024) that, for older patients aged ≥75 years with polypharmacy, who are clinically stable, and on therapeutic VKA, maintaining VKA treatment rather than switching to a DOAC may be considered in order to minimise bleeding risk.
Warfarin vs. acenocoumarol
Warfarin has been the traditional comparator drug used in clinical trials and observational studies against the DOACs (Connolly et al., 2009; Granger et al., 2011; Patel et al., 2011; Giugliano et al., 2013; Larsen et al., 2016; Nielsen et al., 2017; Domek et al., 2020). The results of our study revealed statistically significant and clinically relevant differences in both effectiveness and safety between the two VKA drugs. Compared to acenocoumarol, warfarin was associated with an increased risk of systemic embolism (by at least 16%), an increased risk of all-cause mortality (by 5%), and an increased risk of major bleeding events (by 37%), taking the lower limit of the confidence interval for the relative risk (Table 3).
These data suggest that conclusions drawn from studies using warfarin as the reference drug should not be automatically extrapolated to acenocoumarol.
Strengths
Firstly, the large sample size and long follow-up period of our study allow extrapolation to the general population, providing high external validity.
Secondly, our study provides a comparative evaluation of two VKAs with different half-lives, acenocoumarol and warfarin, in a routine clinical practice setting. This may explain the differing outcomes obtained when comparing them with DOACs and with each other. We also emphasise that the study was conducted by a multidisciplinary team and that there were no conflicts of interest with the pharmaceutical industry.
A third strength of our study is the systematic and thorough collection of mortality data from official records, which contributes to the dependability and reliability of the data. This approach sets our study apart from others, as demonstrated by the numbers and proportions of event types included.
Fourth, our methodology enabled us to select various sources, types of data, as well as essential variables for use in future studies, and to assess the usefulness of each variable in terms of reliability, consistency, and proportion of missing data. In addition, we standardised and normalised disease diagnoses found in various classifications (ICD-9, ICD-10, and the International Classification of Primary Care) and developed integrated scales and clinical evaluations of patients based on different demographic information, comorbidities, use of resources, etc.
Fifthly, a dual approach was adopted using R and Stata applications. This was done firstly to ensure more robust results, and secondly to use the most efficient, user-friendly, and widely-recognised instructions and utilities. We consider the establishment of a methodology and essential information for future drug evaluation studies, including those from other therapeutic groups within our organisation, to be a significant achievement. Consequently, we have preserved and clearly explained all syntax and commands for each step. With our large and diverse multidisciplinary team, we envisage this study laying the foundation of a series of prioritised studies to be developed over time. Any necessary improvements and specific issues related to each study can be incorporated as needed.
Finally, a methodology was chosen that enabled us to transition from a purely observational cohort study of patients to a quasi-experimental study emulating a clinical trial in the comparative evaluation of six oral anticoagulants.
Limitations
The first limitation is the risk of bias inherent in any observational study. We attempted to control for this by performing a propensity score-matched cohort analysis (PSM).
A second limitation of observational studies is the possibility of missing data in the patients’ clinical records. In our study, we extracted each patient’s clinical information from several databases (primary care, hospital, outpatient, emergency), in order to facilitate the collection of essential variables for the analysis.
Thirdly, we did not include all relevant information, as not all patients had it available. This included information such as smoking, alcohol consumption, diet or the patient’s index of dependence on others for performing basic activities of daily living (Barthel test). Therefore, as with any study, we cannot rule out the possibility that these and other unaccounted for variables may have influenced the results.
Fourthly, there is a possible risk of selection bias for oral anticoagulants, since DOACs require prior authorisation and their dosage depends on the patient’s renal function. Depending on the referring hospital, there may also be differences in the choice of VKA. This issue was addressed through the implementation of propensity score-matched cohort analysis (PSM).
Fifthly, we operated on the assumption that the medication dispensed was the same as the one taken by the patient, as determined by other investigators (Rosa et al., 2018; Durand et al., 2021; Ingason et al., 2021; Lip et al., 2021). However, discrepancies between the prescribed medication and the medication actually taken cannot be entirely ruled out in some cases.
Sixthly, for patients who changed medication during the follow-up period (“per-protocol analysis”), we did not conduct a health outcomes analysis based on the most recent OAC the patient was taking when the event was recorded. However, our analysis using the intention-to-treat approach is consistent with standard methods employed in other observational studies (Larsen et al., 2016; Nielsen et al., 2017; Durand et al., 2020; Rodríguez-Bernal et al., 2021; Lau et al., 2022) and in clinical trials (Connolly et al., 2009; Granger et al., 2011; Patel et al., 2011; Giugliano et al., 2013). In the near future, we will perform a complementary pre-treatment analysis, examining outcomes based on the treatment that patients were receiving immediately before an effectiveness or safety event occurred, rather than looking only at their initial treatment assignment.
Seventhly, we did not consider an analysis of low-dose versus standard dose of DOACs.
We also did not assess the relative efficiency or cost implications of the switch in prescription practice from VKA to DOACs. These are open questions for future analysis.
Eighthly, it was not possible to access the INR registry databases in order to evaluate the proper control of patients undergoing VKA treatment. We are aware that poor INR monitoring in acenocoumarol users may result in an artificial inflation of the risks of bleeding or mortality, thereby favouring DOACs. Therefore, our findings should be interpreted with caution in light of this limitation. This aspect will be the subject of future investigations by our research team.
Ninthly, we did not take into account the time that elapsed between the AF diagnosis and the initiation of anticoagulant therapy. This may contribute to potential confounding by disease duration, which could differ systematically across drug groups. Nevertheless, the study exclusively included new users of oral anticoagulants, defined as patients who had not used OACs during the 12 months preceding their enrolment in the study, regardless of any prior use.
Conclusion
Our study has shown that, in general, there are some differences in the effectiveness and safety of DOACs compared to acenocoumarol. Due to a lack of sufficient data, a clear evaluation of edoxaban could not be made. The other three DOACs (dabigatran, rivaroxaban and apixaban) did not demonstrate a significant reduction in ischaemic events.
The three DOACs appear to be safer, with a lower incidence of severe bleeding compared to acenocoumarol, except for dabigatran, which achieved good results in the incidence of intracranial bleeding but had a higher incidence of gastrointestinal bleeding compared to acenocoumarol. These results are supported by the analysis of paired data disaggregated by sex. The values observed for both sexes are comparable to and consistent with the global data. Warfarin was less effective and less safe in comparison with acenocoumarol.
The results provide evidence from real-world clinical practice which needs to be confirmed in well-designed clinical trials involving acenocoumarol. Both types of study, RWD and clinical trials, complement each other and allow us to expand our knowledge of treatments and determine the role of new drugs in the therapeutic arsenal. The focus of our future research will be on analysing patient subgroups with certain comorbidities (e.g., diabetes, renal insufficiency), older patients, with inadequate anticoagulation (time in therapeutic range) for VKA drugs, and the use of low-dose DOACs compared to acenocoumarol. Our results are expected to be of great value when updating clinical practice guidelines in healthcare systems. This information will be important for clinicians in selecting and personalising treatments for patients with atrial fibrillation.
Data availability statement
Data cannot be shared publicly because they are based on the Andalusian Health Population Database. Data are available from the Subdirección Técnica Asesora de Gestión de la Información, Andalusian Health System, Seville, Spain, for researchers who meet the criteria for access to confidential data.
Ethics statement
The studies involving humans were approved by Andalusian Biomedical Research Ethics Coordinating Committee. Andalusia, Spain. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
MM-B: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. JL-M: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing. AL-A: Data curation, Formal Analysis, Methodology, Software, Validation, Writing – review and editing. BG-R: Conceptualization, Investigation, Methodology, Project administration, Supervision, Writing – review and editing. PS-V: Formal Analysis, Methodology, Validation, Writing – review and editing. RI-G: Writing – review and editing. MA-R: Conceptualization, Supervision, Writing – review and editing. LL: Conceptualization, Methodology, Supervision, Writing – review and editing. MM-L: Conceptualization, Methodology, Supervision, Visualization, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Andalusian Public Foundation for Health Research Management in Seville (FISEVI) of the Andalusian Ministry of Health (AP-0222-2019 and AP-0379-2023-C4-F2) and by the Instituto de Salud Carlos III (PI20/01281), Spanish Ministry of Health (co-financed by the European Regional Development Funds).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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/fphar.2025.1548298/full#supplementary-material
References
Anguita, M., Bertiomeu, V., Ruiz, M., Cequier, A., Roldan, I., Muñiz, J., et al. (2020). Direct oral anticoagulants versus vitamin K antagonists in real-world patients with nonvalvular atrial fibrillation. The FANTASIIA study. Rev. Esp. Cardiol. 73 (1), 14–20. doi:10.1016/j.rec.2019.02.021
Austin, P. C. (2011). An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivar. Behav. Res. 46 (3), 399–424. doi:10.1080/00273171.2011.568786
Austin, P. C., and Fine, J. P. (2017). Practical recommendations for reporting fine-gray model analyses for competing risk data. Statistics Med. 36 (27), 4391–4400. doi:10.1002/sim.7501
Chan, Y. H., Lee, H. F., and Chao, T. F. (2019b). “Real-world comparisons of direct oral anticoagulants for stroke prevention in Asian patients with non-valvular atrial fibrillation: a systematic review and meta-analysis,” in Cardiovascular drugs and therapy. Springer, 701–710. doi:10.1007/s10557-019-06910-z
Chan, Y. H., Lee, H. F., See, L. C., Tu, H. T., Chao, T. F., and Yeh, Y. H. (2019a). Effectiveness and safety of four direct oral anticoagulants in Asian patients with nonvalvular atrial fibrillation. Chest 156 (3), 529–543. doi:10.1016/j.chest.2019.04.108
Ciminata, G., Geue, C., Wu, O., Deidda, M., Kreif, N., and Langhorne, P. (2022). Propensity score methods for comparative-effectiveness analysis: a case study of direct oral anticoagulants in the atrial fibrillation population. PloS one 17 (1), e0262293. doi:10.1371/journal.pone.0262293
Connolly, S. J., Ezekowitz, M. D., Yusuf, S., Eikelboom, J., Oldgren, J., Parekh, A., et al. (2009). Dabigatran versus warfarin in patients with atrial fibrillation. N. Engl. J. Med. 361 (12), 1139–1151. doi:10.1056/nejmoa0905561
Domek, M., Gumprecht, J., Ding, W. Y., Lip, G. Y. H., and Lane, D. A. (2020). Practice-derived data on non-vitamin K antagonist oral anticoagulant therapy to complement observations from randomized trials. Eur. Heart J. 22, I1-I12–I12. doi:10.1093/EURHEARTJ/SUAA100
Doni, K., Bühn, S., Weise, A., Mann, N. K., Hess, S., Sönnichsen, A., et al. (2024). Safety outcomes of direct oral anticoagulants in older adults with atrial fibrillation: a systematic review and meta-analysis of (subgroup analyses from) randomized controlled trials. GeroScience 46 (1), 923–944. doi:10.1007/s11357-023-00825-2
Douros, A., Durand, M., Doyle, C. M., Yoon, S., Reynier, P., and Filion, K. B. (2019). “Comparative effectiveness and safety of direct oral anticoagulants in patients with atrial fibrillation: a systematic review and meta-analysis of observational studies,” Drug Saf. 42. 1135–1148. doi:10.1007/s40264-019-00842-1
Durand, M., Schnitzer, M. E., Pang, M., Carney, G., Eltonsy, S., Filion, K. B., et al. (2020). Comparative effectiveness and safety of direct oral anticoagulants versus vitamin K antagonists in nonvalvular atrial fibrillation: a Canadian multicentre observational cohort study. CMAJ open 8 (4), E877–E886. doi:10.9778/cmajo.20200055
Durand, M., Schnitzer, M. E., Pang, M., Carney, G., Eltonsy, S., Filion, K. B., et al. (2021). Effectiveness and safety among direct oral anticoagulants in nonvalvular atrial fibrillation: a multi-database cohort study with meta-analysis. Br. J. Clin. Pharmacol. 87 (6), 2589–2601. doi:10.1111/bcp.14669
Gharaibeh, A., Koppikar, S., and Bonilla-Escobar, J. (2014). Strengthening the reporting of observational studies in epidemiology (STROBE) in the international journal of medical students. Int. J. Med. Students 2 (2), 36–37. doi:10.5195/ijms.2014.76
Giugliano, R. P., Ruff, C. T., Braunwald, E., Murphy, S. A., Wiviott, S. D., Halperin, J. L., et al. (2013). Edoxaban versus warfarin in patients with atrial fibrillation. N. Engl. J. Med. 369 (22), 2093–2104. doi:10.1056/nejmoa1310907
Granger, C. B., Alexander, J. H., McMurray, J. J. V., Lopes, R. D., Hylek, E. M., Hanna, M., et al. (2011). Apixaban versus warfarin in patients with atrial fibrillation. N. Engl. J. Med. 365 (11), 981–992. doi:10.1056/nejmoa1107039
Hindricks, G., Potpara, T., Dagres, N., Arbelo, E., Bax, J. J., Blomström-Lundqvist, C., et al. (2021). 2020 ESC guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the european association for cardio-thoracic surgery (EACTS): the task force for the diagnosis and management of atrial fibrillation of the european society of cardiology (ESC) developed with the special contribution of the european heart rhythm association (EHRA) of the ESC. Eur. Heart J. 42 (5), 373–498. doi:10.1093/eurheartj/ehaa612
Ingason, A. B., Hreinsson, J. P., Ágústsson, A. S., Lund, S. H., Rumba, E., Pálsson, D. A., et al. (2021). Rivaroxaban is associated with higher rates of gastrointestinal bleeding than other direct oral anticoagulants a nationwide propensity score-weighted study. Ann. Intern. Med. 174 (11), 1493–1502. doi:10.7326/M21-1474
Jaksa, A., Gibbs, L., Kent, S., Rowark, S., Duffield, S., Sharma, M., et al. (2022). Using primary care data to assess comparative effectiveness and safety of apixaban and rivaroxaban in patients with nonvalvular atrial fibrillation in the UK: an observational cohort study. BMJ Open 12 (10), e064662. doi:10.1136/bmjopen-2022-064662
January, C. T., Wann, L. S., Calkins, H., Chen, L. Y., Cigarroa, J. E., Cleveland, J. C., et al. (2019). 2019 AHA/ACC/HRS focused update of the 2014 AHA/ACC/HRS guideline for the management of patients with atrial fibrillation: a report of the American college of cardiology/american heart association task force on clinical practice guidelines and the heart rhythm society in collaboration with the society of thoracic surgeons. Circulation 140 (2), e125–e151. doi:10.1161/CIR.0000000000000665
Joosten, L. P. T., van Doorn, S., van de Ven, P. M., Köhlen, B. T. G., Nierman, M. C., Koek, H. L., et al. (2024). Safety of switching from a vitamin K antagonist to a non-vitamin K antagonist oral anticoagulant in frail older patients with atrial fibrillation: results of the FRAIL-AF randomized controlled trial. Circulation 149 (4), 279–289. doi:10.1161/CIRCULATIONAHA.123.066485
Kim, S. M., Jeon, E. T., Jung, J. M., and Lee, J. S. (2021). “Real-world oral anticoagulants for Asian patients with non-valvular atrial fibrillation: a PRISMA-Compliant article,” Med. (United States) 100. E26883. doi:10.1097/MD.0000000000026883
Larsen, T. B., Skjøth, F., Nielsen, P. B., Kjældgaard, J. N., and Lip, G. Y. H. (2016). Comparative effectiveness and safety of non-Vitamin K antagonist oral anticoagulants and warfarin in patients with atrial fibrillation: propensity weighted nationwide cohort study. BMJ (Online) 353, i3189. doi:10.1136/bmj.i3189
Lau, W. C. Y., Torre, C. O., Man, K. K. C., Stewart, H. M., Seager, S., Van Zandt, M., et al. (2022). Comparative effectiveness and safety between apixaban, dabigatran, edoxaban, and rivaroxaban among patients with atrial fibrillation: a multinational population-based cohort study. Ann. Intern. Med. 175 (11), 1515–1524. doi:10.7326/M22-0511
Li, W. J., Archontakis-Barakakis, P., Palaiodimos, L., Kalaitzoglou, D., Tzelves, L., Manolopoulos, A., et al. (2021). Dabigatran, rivaroxaban, and apixaban are superior to warfarin in Asian patients with non-valvular atrial fibrillation: an updated meta-analysis. World J. Cardiol. 13 (4), 82–94. doi:10.4330/wjc.v13.i4.82
Lip, G. Y. H., Keshishian, A. V., Zhang, Y., Kang, A., Dhamane, A. D., Luo, X., et al. (2021). Oral anticoagulants for nonvalvular atrial fibrillation in patients with high risk of gastrointestinal bleeding. JAMA Netw. Open 4, e2120064. doi:10.1001/jamanetworkopen.2021.20064
Mamas, M. A., Batson, S., Pollock, K. G., Grundy, S., Matthew, A., Chapman, C., et al. (2022). Meta-analysis comparing apixaban versus rivaroxaban for management of patients with nonvalvular atrial fibrillation. Am. J. Cardiol. 166, 58–64. doi:10.1016/j.amjcard.2021.11.021
Martínez-González, M. A., and Alonso, A. (2008). Qué es una hazard ratio? Med. Clin. Barc. 131 (2), 65–72. doi:10.1157/13123495
Montero-Balosa, M. C., Limón-Mora, J. A., Leal-Atienza, A., Luque-Romero, L. G., Aguado-Romeo, M. J., Isabel-Gómez, R., et al. (2023). Effectiveness and safety of oral anticoagulant therapy in a real-world cohort with atrial fibrillation: the SIESTA-A study protocol. PloS one 18 (11), e0294822. doi:10.1371/journal.pone.0294822
Muñoyerro-Muñiz, D., Goicoechea-Salazar, J. A., García-León, F. J., Laguna-Téllez, A., Larrocha-Mata, D., and Cardero-Rivas, M. (2020). Health record linkage: andalusian health population database. Gac. Sanit. 34 (2), 105–113. doi:10.1016/j.gaceta.2019.03.003
Nielsen, P. B., Skjøth, F., Søgaard, M., Kjældgaard, J. N., Lip, G. Y. H., and Larsen, T. B. (2017). Effectiveness and safety of reduced dose non-vitamin K antagonist oral anticoagulants and warfarin in patients with atrial fibrillation: propensity weighted nationwide cohort study. BMJ (Online) 356, j510. doi:10.1136/bmj.j510
Patel, M. R., Mahaffey, K. W., Garg, J., Pan, G., Singer, D. E., Hacke, W., et al. (2011). Rivaroxaban versus warfarin in nonvalvular atrial fibrillation. N. Engl. J. Med. 365 (10), 883–891. doi:10.1056/NEJMoa1009638
Ramagopalan, S. V., Sicras-Mainar, A., Polanco-Sanchez, C., Carroll, R., and de Bobadilla, J. F. (2019). Patient characteristics and stroke and bleeding events in nonvalvular atrial fibrillation patients treated with apixaban and vitamin K antagonists: a Spanish real-world study. J. Comp. Eff. Res. 8 (14), 1201–1212. doi:10.2217/cer-2019-0079
Rodriguez-Bernal, C. L., Peiró, S., Hurtado, I., García-Sempere, A., and Sanfélix-Gimeno, G. (2018). Primary nonadherence to oral anticoagulants in patients with atrial fibrillation: real-world data from a population-based cohort. J. Manag. Care and Specialty Pharm. JMCP May 24, 440–448. doi:10.18553/jmcp.2018.24.5.440
Rodríguez-Bernal, C. L., Santa-Ana-Téllez, Y., García-Sempere, A., Hurtado, I., Peiró, S., and Sanfélix-Gimeno, G. (2021). Clinical outcomes of nonvitamin K oral anticoagulants and acenocoumarol for stroke prevention in contemporary practice: a population-based propensity-weighted cohort study. Br. J. Clin. Pharmacol. 87 (2), 632–643. doi:10.1111/bcp.14430
Rosa, M. (2018). Anticoagulation control with acenocoumarol or warfarin in non-valvular atrial fibrillation in primary care (Fantas-TIC study). doi:10.3390/ijerph
Tepper, P. G., Mardekian, J., Masseria, C., Phatak, H., Kamble, S., Abdulsattar, Y., et al. (2018). Real-world comparison of bleeding risks among non-valvular atrial fibrillation patients prescribed apixaban, dabigatran, or rivaroxaban. PLoS ONE 13 (11), e0205989. doi:10.1371/journal.pone.0205989
van den Ham, H. A., Souverein, P. C., Klungel, O. H., Platt, R. W., Ernst, P., Dell'Aniello, S., et al. (2021). Major bleeding in users of direct oral anticoagulants in atrial fibrillation: a pooled analysis of results from multiple population-based cohort studies. Pharmacoepidemiol. Drug Saf. 30 (10), 1339–1352. doi:10.1002/pds.5317
Van Gelder, I. C., Rienstra, M., Bunting, K. V., Casado-Arroyo, R., Caso, V., Crijns, H. J. G. M., et al. (2024). 2024 ESC Guidelines for the management of atrial fibrillation developed in collaboration with the European Association for Cardio-Thoracic Surgery (EACTS). Eur. Heart J. 45 (36), 3314–3414. doi:10.1093/eurheartj/ehae176
Wallentin, L., Yusuf, S., Ezekowitz, M. D., Alings, M., Flather, M., Franzosi, M. G., et al. (2010). Efficacy and safety of dabigatran compared with warfarin at different levels of international normalised ratio control for stroke prevention in atrial fibrillation: an analysis of the RE-LY trial. Lancet 376 (9745), 975–983. doi:10.1016/S0140-6736(10)61194-4
Xue, Z., and Zhang, H. (2019). Non–vitamin K antagonist oral anticoagulants versus warfarin in Asians with atrial fibrillation meta-analysis of randomized trials and real-world studies. Stroke 50 (10), 2819–2828. doi:10.1161/STROKEAHA.119.026054
Keywords: oral anticoagulant agents, oral direct anticoagulants, safety, health outcomes, effectiveness
Citation: Montero-Balosa MC, Limón-Mora JA, Leal-Atienza A, García-Robredo B, Sánchez-Villegas P, Isabel-Gómez R, Aguado-Romeo MJ, Luque Romero LG and Molina-López MT (2025) Comparative clinical outcomes of acenocoumarol versus direct oral anticoagulants (DOACs) and warfarin in patients with atrial fibrillation: real-world-evidence (SIESTA-A study). Front. Pharmacol. 16:1548298. doi: 10.3389/fphar.2025.1548298
Received: 19 December 2024; Accepted: 18 July 2025;
Published: 01 August 2025.
Edited by:
Christos Kontogiorgis, Democritus University of Thrace, GreeceReviewed by:
Ju-Yeun Lee, Seoul National University, Republic of KoreaJulian Librero, NavarraBiomed, Spain
Copyright © 2025 Montero-Balosa, Limón-Mora, Leal-Atienza, García-Robredo, Sánchez-Villegas, Isabel-Gómez, Aguado-Romeo, Luque Romero and Molina-López. 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: Ma Carmen Montero-Balosa, bWMubW9udGVyby5zc3BhQGp1bnRhZGVhbmRhbHVjaWEuZXM=