Article Text
Abstract
Objective Antithrombotic therapy is essential for patients with atrial fibrillation (AF) and stable coronary artery disease (CAD) because of the high risk of thrombosis, whereas a combination of antiplatelets and anticoagulants is associated with a high risk of bleeding. We sought to develop and validate a machine-learning-based model to predict future adverse events.
Methods Data from 2215 patients with AF and stable CAD enrolled in the Atrial Fibrillation and Ischaemic Events With Rivaroxaban in Patients With Stable Coronary Artery Disease trial were randomly assigned to the development and validation cohorts. Using the random survival forest (RSF) and Cox regression models, risk scores were developed for net adverse clinical events (NACE) defined as all-cause death, myocardial infarction, stroke or major bleeding.
Results Using variables selected by the Boruta algorithm, RSF and Cox models demonstrated acceptable discrimination and calibration in the validation cohort. Using the variables weighted by HR (age, sex, body mass index, systolic blood pressure, alcohol consumption, creatinine clearance, heart failure, diabetes, antiplatelet use and AF type), an integer-based risk score for NACE was developed and classified patients into three risk groups: low (0–4 points), intermediate (5–8) and high (≥9). In both cohorts, the integer-based risk score performed well, with acceptable discrimination (area under the curve 0.70 and 0.66, respectively) and calibration (p>0.40 for both). Decision curve analysis showed the superior net benefits of the risk score.
Conclusions This risk score can predict the risk of NACE in patients with AF and stable CAD.
Trial registration numbers UMIN000016612, NCT02642419.
- atrial fibrillation
- coronary artery disease
- risk factors
Data availability statement
Data are available on reasonable request. The data underlying this article will be shared on reasonable request to the corresponding author.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
A combination of antiplatelet and anticoagulant agents has been an essential antithrombotic therapy for patients suffering from atrial fibrillation (AF) and coronary artery disease (CAD); however, the persistent issue of bleeding complications associated with antithrombotic therapy remains a concern.
WHAT THIS STUDY ADDS
This subanalysis of the Atrial Fibrillation and Ischaemic Events With Rivaroxaban trial developed a machine-learning-based risk prediction model for net adverse clinical events in patients with AF and CAD.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
Managing the risk of both thromboembolic and bleeding events is crucial in the antithrombotic treatment of those patients. Using this risk score can prove beneficial in the decision-making process in clinical settings.
Introduction
Antiplatelet therapy is required for the secondary prevention of cardiovascular events in patients with coronary artery disease (CAD),1–3 and anticoagulant therapy is essential for the prevention of thromboembolic events in patients with atrial fibrillation (AF).4–6 A combination of antiplatelet and anticoagulant agents has been used in patients with AF and CAD; however, the high risk of antithrombotic therapy-related bleeding complications remains a problem. The AFIRE trial (Atrial Fibrillation and Ischaemic Events with Rivaroxaban in Patients with Stable Coronary Artery Disease) showed that rivaroxaban monotherapy was non-inferior in at reducing the risk of cardiovascular events or all-cause death and superior in at reducing the risk of major bleeding compared with a combination of rivaroxaban with a single antiplatelet drug in AF patients with stable CAD at ≥1 year after revascularisation or those with angiographically confirmed CAD not requiring revascularisation.7
Simple conventional risk scores, such as the CHADS2,8 CHA2DS2-VASc9 or HAS-BLED10 scores, have been widely used to estimate the individual risk of thromboembolic or bleeding events in patients with AF. In addition, the CHADS2 score is a useful prediction tool for cardiovascular or cerebrovascular events, even in CAD patients without AF.11 12 However, subanalysis of the AFIRE trial demonstrated that patients with major bleeding had a high risk of subsequent major adverse cardiac and cerebrovascular events, especially within 30 days after major bleeding,13 suggesting that a balanced risk assessment should be established to estimate the integrated risk of both thromboembolic and bleeding events in patients with AF and CAD.
To address this challenging issue, we conducted a post hoc analysis of the AFIRE trial to develop and validate a machine-learning-based risk prediction model for future net adverse clinical events (NACE) among patients with AF and stable CAD.
Methods
Study design and study participants
This study was a post hoc analysis of the AFIRE trial, a multicentre, randomised, open-label, parallel-group trial. The detailed study design, protocol, and results of the primary analysis of the AFIRE trial have been reported previously.7 14 Briefly, the AFIRE trial was conducted in 294 hospitals across Japan between 23 February 2015 and 30 September 2017 and included patients aged ≥20 years, diagnosed with AF with a CHADS2 score ≥1 and stable CAD at ≥1 year after revascularisation or those with angiographically confirmed CAD not requiring revascularisation. Patients with a history of stent thrombosis, coexisting active cancer or poorly controlled hypertension were excluded in the trial. Patients were randomised in a 1:1 ratio to receive either rivaroxaban (10 mg once daily for patients with a creatinine clearance rate (CCR) of 15–49 mL/min or 15 mg once daily for patients with a CCR≥50 mL/min) alone or rivaroxaban plus an antiplatelet drug (either aspirin or P2Y12 inhibitor). Patient follow-up was performed at baseline, 6 months and at the end of the trial, with additional follow-up for routine clinical care if needed. The study follow-up period was at least 24 months and up to 45 months.
The data were reviewed by an independent data and safety monitoring committee. In this post hoc analysis, 2215 patients in the modified intention-to-treat population from the AFIRE trial were analysed.
Outcome
The primary outcome of this study was NACE, consisting of all-cause death, myocardial infarction, stroke or major bleeding according to the criteria of the International Society on Thrombosis and Haemostasis.15
Candidate variables
After excluding of variables with >10% missing values and a correlation coefficient >0.7, 78 candidate variables at baseline were used to develop the risk prediction model as potential predictors of NACE (online supplemental table 1). Missing values were imputed by missForest,16 which is a non-parametric imputation method using a random forest model that can learn non-linearity, easily handle mixed-type data and calculate out-of-bag errors.
Supplemental material
Statistical analysis
Based on the outcome, study participants were randomly assigned to a development (50%) or validation (50%) dataset. The development dataset was further divided into two (a training or tuning set, each 50%) to tune the hyperparameters. Categorical variables are presented as frequencies and percentages, and continuous variables are presented as medians and IQR. Using only data from the development dataset, feature selection and relationship modelling for risk prediction model development were conducted (figure 1).
Feature selection
From all candidate variables, we further identified 13 predictors (age, sex, treatment randomisation, body mass index (BMI), systolic blood pressure, diastolic blood pressure, CCR, type of AF, alcohol consumption, diabetes mellitus, heart failure, prior stent implantation and prior coronary artery bypass grafting) using a feature selection method with the Boruta algorithm, which is a wrapper around a random forest classification algorithm that can provide a numerical estimate of the feature importance without tuning the parameters and removing the variables that are less relevant than random probes by a statistical test.17
Model development and evaluation
The relationships between the outcomes and variables selected by the Boruta algorithm were assessed using Cox proportional hazard (PH) models and random survival forest (RSF) models,18 which is an ensemble tree-based classification method for the analysis of right-censored survival data. Each model was assessed based on Harrell’s concordance index (C-index) for discrimination performance and the Brier score for both discrimination and calibration performance in the development and validation datasets.
Risk score development and validation
For clinical use, a simple integer-based risk score to predict the outcome was developed based on each variable’s unadjusted HR, with statistical significance calculated from the Cox PH model in the development dataset. Continuous variables were categorised as follows: age (≥80 years, <80 years), BMI (<18.5, 18.5–25, ≥25), systolic blood pressure (<90, 90–140, ≥140), diastolic blood pressure (<60, 60–90, ≥90) and CCR (<30, 30–50, ≥50). The discrimination performance of the integer-based risk score was assessed using the area under the receiver operating characteristic curve (AUC-ROC) and its 95% CIs in the development and validation datasets. The calibration performance of the risk score was assessed using a calibration plot and Hosmer-Lemeshow test. Performance comparison of the new risk score with that of conventional risk scores, such as the CHADS2 (Congestive heart failure, Hypertension, Age ≥75 y, Diabetes mellitus, Stroke or TIA), CHA2DS2-VASc (Congestive heart failure, Hypertension, Age ≥75 y, Diabetes mellitus, Stroke or TIA, Vascular disease, Age 65-74 y, Sex category [female]) and HAS-BLED (Hypertension, Abnormal renal/liver function, Stroke, Bleeding, Labile INR, Elderly, Drugs) scores, was performed using ROC analysis, category and category-free net reclassification improvement, integrated discrimination improvement and decision curve analysis to indicate net benefit and clinical utility.19–21
Statistical analyses were performed using SPSS (V.23.0; IBM), R software, V.4.0.5 (The R Project for Statistical Computing) and Python (V.3.7.11, Python Software Foundation). Statistical significance was defined as a two-sided p<0.05.
Results
Study participants
Among the 2215 patients in the AFIRE trial, 1107 were included in the development cohort (median (IQR) age, 75 (69–80) years; 870 (78.6%) males) and 1108 in the validation cohort (median (IQR) age, 75 (69, 80) years; 881 (79.5%) males). During a median (IQR) follow-up of 24.1 (17.3–31.5) months, 215 patients (9.7%) suffered NACE (107 and 108 in the development and validation cohorts, respectively). Patient characteristics for the development and validation cohorts are presented in table 1. There were no statistical differences in demographic and physiological findings, antithrombotic regimen, type of AF, comorbidities, medical history, or prior revascularisation between the cohorts, except for the location of the culprit lesion, interventions other than revascularisation and prior bleeding complications.
Performance of development models and integer-based risk score
The 13 Boruta algorithm-selected variables had high C-index of 0.706 and 0.667 in the development and validation cohorts, respectively, when used with the RSF model (online supplemental table 2). Using the same variables as the Cox PH model, the discrimination performance was acceptable (C-index of 0.680 and 0.650 in the development and validation cohorts, respectively). The RSF model had low Brier scores of 0.081 and 0.080 in the development and validation cohorts, respectively (online supplemental table 2), indicating good calibration. Similarly, the Cox PH model had low Brier scores of 0.079 and 0.081 in the development and validation cohorts, respectively.
The β-coefficients, HRs and 95% CIs for each of the selected variables calculated using the Cox PH model are presented in table 2. Using the variables with statistical significance, an integer-based risk score for NACE was created (figure 2); according to each variable’s HR, a score for the risk of NACE was assigned (table 2). The new risk score model has a theoretical range of 0–23 and is divided into three risk groups: low (0–4 points), intermediate (5–8 points) and high (≥9 points). The new risk score model demonstrated good discrimination with a C-index of 0.70 (95% CI 0.65 to 0.75; figure 3A) and acceptable calibration (Hosmer-Lemeshow test; χ2 of 7.789, p=0.454; online supplemental figure 1a). The cumulative 3-year incidence of NACE increased in a graded fashion across the risk groups in the development cohort (figure 4A).
Validation and decision curve analysis of the integer-based risk score
In the validation cohort, the AUC of the new risk score was 0.66 (95% CI 0.60, 0.72; figure 3B) and the Hosmer-Lemeshow goodness-of-fit χ2 was 3.417 (p=0.906; online supplemental figure 1), indicating moderate discrimination and acceptable calibration, respectively. The cumulative 3-year incidence of NACE also increased in a graded fashion across the three risk groups in the validation cohort (figure 4B). The AUC of conventional risk scores such as the CHADS2, CHA2DS2-VASc and HAS-BLED scores were 0.57 (95% CI 0.51 to 0.62), 0.57 (95% CI 0.51 to 0.62) and 0.51 (95% CI 0.46 to 0.57), respectively (online supplemental figure 2). Compared with these conventional risk scores, the discrimination performance of the new risk score was significantly higher (online supplemental figure 2 and online supplemental table 3). The decision curves for the new and conventional risk scores for predicting NACE in the validation cohort are shown in figure 5A. The new risk score had a greater net benefit, with threshold probabilities of 5%–25% than the conventional risk scores. A 10-fold cross-validated decision curve of the new risk score showed a similar standardised net benefit at each threshold probability (figure 5B), indicating the robustness of the results.
Discussion
This post hoc analysis of the AFIRE trial demonstrated the development and validation of a machine-learning-based risk score that integrated easily available clinical and laboratory data to predict future adverse clinical events in patients with AF and stable CAD. Patients in the high-risk group had a 3-year NACE risk of approximately 40%. This risk score substantially improved the discriminatory ability and clinical usefulness of adverse clinical events compared with conventional risk scores. Antithrombotic therapy is needed for patients with AF and stable CAD because of the high risk of ischaemic and thrombotic events, whereas a combination of antiplatelet therapy with anticoagulants is associated with a high risk of bleeding events. Given this clinically unresolvable dilemma, this machine learning-based risk score is comprehensive, easily available and clinically useful for predicting future adverse events. Based on this risk score, clinicians should reconsider patient management, and close follow-up and cardiovascular healthcare should be provided to high-risk patients.
Machine learning has advantages in prediction model development compared with traditional statistical methods that focus on inference and do not require a prior assumption of causality in variable selection and modelling. Many machine-learning-based risk prediction models have been reported for the diagnosis and prognosis of patients with CAD22–25 or AF,26–28 whereas there has only been one previous report of a machine learning-based prediction model for CAD patients complicated with AF.29 However, only one report showed a machine-learning prediction model for all-cause death among these patients. No established machine-learning-based risk score has been reported for the prediction of NACE combining thromboembolic events with bleeding events; therefore, this was the first study to develop a machine-learning-based risk prediction score for NACE with internal validation. In this study, we used the Boruta algorithm for variable selection and an RSF model in the development process of the prediction model, as well as the traditional Cox PH model for modelling and scoring. The Boruta algorithm is a novel feature selection algorithm to identify all relevant variables for outcomes and a wrapper built around a random forest classification algorithm, which can be performed quickly without tuning the parameters and provides a numerical estimate of feature importance.17 RSF is a random forest method for the analysis of right-censored survival data18 and can mathematically build binary recursive trees for all samples and obtain the maximal survival difference across daughter nodes with the application of bootstrap methods and the log-rank splitting rule.30 31 Although traditional statistical methods for survival data, such as the Cox PH model, rely on restrictive assumptions such as PHs, machine-learning-based methods can manage large multidimensional datasets of right-censored survival data without the need for assumptions of parametric distributions, interaction between variables, linear relationships with outcome and overfitting of models. Therefore, to avoid these mathematical issues, machine-learning approaches may be able to predict complex clinical outcomes such as NACE more accurately. Moreover, the integer-based risk score created in this study would be easy to use in clinical practice because it does not require nesting on a web-based platform or electronic medical record.
In patients with AF, the CHADS28 and CHA2DS2-VASc scores9 have been established as risk scores for thromboembolic events and the HAS-BLED score10 as the risk score for major bleeding events. However, this study showed that these conventional risk scores did not have a good predictive performance for NACE in patients with AF and stable CAD. Here, we developed and validated a machine learning-based risk score for NACE in these patients, and the risk score had modest discrimination and good calibration with good reclassification improvement compared with conventional risk scores. One of the responsible factors may be that the variable selection of the integer-based risk score was conducted by combining machine learning (Boruta algorithm) and classical statistical (Cox PH) methods. The variables selected for the risk score in this study, such as older age, low BMI, female sex, CCR, heart failure, diabetes and combination with antiplatelets, were included in other traditional risk scores for the assessment of thrombotic and bleeding events.10 32 33 Although most variables selected for the risk score can be interpreted and may be useful in understanding the underlying mechanism of clinical adverse events in patients with AF and stable CAD, the association between no alcohol consumption and high risk of NACE may be difficult to interpret. However, this can be explained by the following mechanism: among patients who do not have a drinking habit, some patients cannot consume alcohol because of decreased aldehyde dehydrogenase 2 (ALDH2) activity due to the ALDH2-deficient variant, which presents as alcohol flushing syndrome and is more common in East Asians than in other ethnic populations.34–37 The ALDH2-deficient variant is a risk factor for ischaemic heart diseases, such as coronary spasms38 and acute myocardial infarction.39 40 Based on these findings, alcohol consumption may have been selected as a candidate risk prediction score for NACE in this study. This machine-learning-based risk score is easy to use because it is based on information that is easily available in clinical settings and does not require testing results or information that is difficult to obtain.
This study had several limitations. First, the study population was Japanese and received a rivaroxaban dose of 10 or 15 mg once daily approved by Japan, rather than the once-daily dose of 20 mg approved globally, which may have caused selection bias. Additionally, because the study population included only patients who met the eligibility criteria of the AFIRE trial, we did not verify that the results of this study are applicable to patients in a real-world setting. To confirm generalisability, the risk score should be validated in other ethnic populations or settings. Second, because relevant information, such as haemoglobin,33 platelet count,33 complex percutaneous coronary intervention41 and control of chronic disease (diabetes mellitus and dyslipidaemia), was not captured in the AFIRE trial, the information was not included in the variable selection. Third, although discontinuation of antithrombotic therapy is a risk factor for thrombotic events, data for antithrombotic therapy adherence were not collected in the AFIRE trial and lack of data may have affected the results of this study.
In conclusion, this post hoc analysis of the AFIRE trial demonstrated the development and validation of a machine learning-based risk score that can predict future adverse clinical events in patients with AF and stable CAD. It is important to balance the risk of both thromboembolic and bleeding events in the antithrombotic management of these patients, and the application of this risk score can be useful for decision-making in clinical settings.
Data availability statement
Data are available on reasonable request. The data underlying this article will be shared on reasonable request to the corresponding author.
Ethics statements
Patient consent for publication
Ethics approval
This study involves human participants and the trial was conducted in accordance with the Declaration of Helsinki and approved by the institutional review boards of the Kumamoto University Hospital (Rinri No.2547) and all participating institutions. The data were reviewed by an independent data and safety monitoring committee. In this post hoc analysis, 2215 patients in the modified intention-to-treat population from the AFIRE trial were analysed. Participants gave informed consent to participate in the study before taking part.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Twitter @MasanobuIshii
Contributors All authors have contributed to this manuscript. MI, EN and SN performed the statistical analysis. HO, KoK and SY handled funding and supervision. KoK, SY, MA, JA, TM, MN, KaM, NH, KaK, AH, KuM and HO acquired the data. KoK, SY, MA, JA, TM, MN, KaM, NH, KaK, AH, KuM and HO conceived and designed the study. MI drafted the manuscript. KoK and KT made critical revisions of the manuscript for key intellectual content. KoK accpets full responsibility for the work and the conduct of the study, had access to the data, and controlled the decision to publish.
Funding This work was supported by the Japan Cardiovascular Research Foundation based on a contract with Bayer Yakuhin.
Competing interests KoK reports remuneration for lectures from Bayer Yakuhin, Daiichi-Sankyo, Novartis Pharma AG and Otsuka Pharmaceutical; has received trust research/joint research funds from Bayer Yakuhin, and Daiichi-Sankyo; and has received scholarship funds from Abbott Medical. SY reports grants from Takeda Pharmaceutical, Abbott and Boston Scientific, and personal fees from Daiichi-Sankyo and Bristol Myers Squibb. MA reports grants from the Japan Agency for Medical Research and Development; personal fees from Bristol Myers Squibb and Nippon Boehringer Ingelheim; and grants and personal fees from Bayer Yakuhin and Daiichi-Sankyo. JA reports personal fees from Bayer Yakuhin and Sanofi and grants and personal fees from Daiichi-Sankyo. TM reports grants from the Japan Cardiovascular Research Foundation, and personal fees from Nippon Boehringer Ingelheim, Daiichi-Sankyo, AstraZeneca and Bayer Yakuhin. MN reports grants and personal fees from Bayer Yakuhin, Daiichi-Sankyo and Sanofi and personal fees from Bristol Myers Squibb and Nippon Boehringer Ingelheim. KaM reports personal fees from Amgen Astellas BioPharma, Astellas Pharma, MSD, Bayer Yakuhin, Sanofi, Takeda Pharmaceutical, Daiichi-Sankyo, Nippon Boehringer Ingelheim and Bristol Myers Squibb. NH reports grants and personal fees from Bayer Yakuhin, Nippon Boehringer Ingelheim and personal fees from Bristol Myers Squibb. KaK reports grants from the Japan Cardiovascular Research Foundation grants and personal fees from Bayer Yakuhin, Daiichi-Sankyo, Sanofi, MSD and AstraZeneca, and personal fees from Bristol Myers Squibb and Nippon Boehringer Ingelheim. AH reports grants and personal fees from Boston Scientific Japan, Otsuka Pharmaceutical, Sanofi, Astellas Pharma, Bristol Myers Squibb, Daiichi-Sankyo, Bayer Yakuhin, Fukuda Denshi, Abbott Japan, Japan Lifeline, Takeda Pharmaceutical and Sumitomo Dainippon Pharma, and personal fees from Toa Eiyo, Nippon Boehringer Ingelheim, Amgen Astellas BioPharma, and AstraZeneca. HO reports personal fees from Towa Pharmaceutical, Bristol Meyers Squibb, Pfizer, Toa Eiyo, Bayer Yakuhin, and Novartis Pharma. KT received significant research grant from AMI, Bayer Yakuhin, Bristol-Myers K.K., EA Pharma, Mochida Pharmaceutical, and scholarship fund from AMI, Bayer Yakuhin, Boehringer Ingelheim Japan, Chugai Pharmaceutical, Daiichi-Sankyo, Edwards Lifesciences Corporation, Johnson & Johnson K.K., Ono Pharmaceutical, Otsuka Pharmaceutical, Takeda Pharmaceutical, and honoraria from Amgen K.K., Bayer Yakuhin, Daiichi-Sankyo, Kowa Pharmaceutical, Novartis Pharma K.K., Otsuka Pharmaceutical, Pfizer Japan, and belongs to the endowed departments donated by Abbott Japan, Boston Scientific Japan K.K., Fides-one, GM Medical, ITI, Kaneka Medix, Nipro, Terumo, Abbott Medical, Cardinal Heaith Japan, Fukuda Denshi, Japan Lifeline, Medical Appliance, Medtoronic Japan.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.