Objective Post-myocardial infarction (MI) care is crucial to preventing recurrent major adverse cardiovascular events (MACE), but can be complicated to personalise. A tool is needed that effectively stratifies risk of cardiovascular (CV) events 1–3 years after MI but is also clinically usable.
Methods Patients surviving ≥1 year after an index MI with ≥1 risk factor for recurrent MI (ie, age ≥65 years, prior MI, multivessel coronary disease, diabetes, glomerular filtration rate <60 mL/min/1.73 m2) were studied. Cox regression derived sex-specific Intermountain Major Adverse Cardiovascular Events (IMACE) risk scores for the composite of 1-year to 3-year MACE (CV death, MI or stroke). Derivation was performed in 70% of subjects (n=1342 women; 3047 men), with validation in the other 30% (n=576 women; 1290 men). Secondary validations were also performed.
Results In women, predictors of CV events were glucose, creatinine, haemoglobin, platelet count, red cell distribution width (RDW), age and B-type natriuretic peptide (BNP); among men, they were potassium, glucose, blood urea nitrogen, haematocrit, white blood cell count, RDW, mean platelet volume, age and BNP. In the primary validation, in women, IMACE ranged from 0 to 11 (maximum possible: 12) and had HR=1.44 per +1 score (95% CI 1.29 to 1.61; P<0.001); men had IMACE range 0–14 (maximum: 16) and HR=1.29 per +1 score (95% CI 1.20 to 1.38; P<0.001). IMACE ≥5 in women (≥6 in men) showed strikingly higher MACE risk.
Conclusions Sex-specific risk scores strongly stratified 1-year to 3-year post-MI MACE risk. IMACE is an inexpensive, dynamic, electronically delivered tool for evaluating and better managing post-MI patient care.
- clinical decision tool
- precision medicine
- learning healthcare system
- implementation science
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 and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
What is already known about this subject?
Long-term use of dual antiplatelet therapy and other post-myocardial infarction (MI) care is effective at reducing major adverse cardiovascular events (MACE), but limiting risks and costs are increasingly a concern.
What does this study add?
Developed on the methodologies of successfully implemented dynamic clinical risk scores, the sex-specific Intermountain Major Adverse Cardiovascular Events (IMACE) models have a greater ability than other clinical decision tools in predicting 3-year MACE and are more practical to deliver and efficiently use in clinical practice.
How might this impact on clinical practice?
Clinicians may receive IMACE and other post-MI risk prediction scores within the clinical workflow via the electronic health record without changing their practice environment or processes, and use IMACE to direct their use of time and other resources towards patients in need of more precise evaluation and care.
Diagnostic and prognostic risk prediction tools are often created with the goal to improve medical care quality.1–8 In patients with a history of myocardial infarction (MI), the quantification and stratification of risk for future major adverse cardiovascular events (MACE) are challenging, in part because most risk tools use predictors or approaches that add time to a clinician’s workload. Various such tools apply to prognostication after acute coronary syndrome.3 6–11 With the recent approval of dual antiplatelet therapy (DAPT) for treatment periods longer than 12 months following acute MI,12–16 the ability to risk-stratify patients early in the care process to better personalise care plans would be beneficial.17
Use of risk tools to guide clinical care following an acute MI may involve higher precision in tailoring the prescription of antiplatelet agents and other medications to a patient’s needs, as well as gauging the required intensity of clinical follow-up, encouraging adherence to therapy, and other approaches. Unfortunately, most clinical decision tools are focused primarily on their predictive ability at the expense of usability. Those developed to be usable, though, such as the CHADS2score,8 sacrifice some predictive ability18–20 for ease of use.
The Intermountain Mortality Risk Score (IMRS) is a sex-specific decision tool created to solve problems that inhibit risk score utilisation in clinical practice while maintaining a substantial ability to predict outcomes (ie, making them both usable and highly predictive).21 IMRS addressed these challenges via data that are routinely collected or are inexpensive to obtain, and are electronically available from standardised, quantitative measurement methods. Furthermore, methods were applied to categorise variables into data-derived groupings (not overfit as continuous variables or ineffectively applied as clinically based thresholds).21 IMRS and similar 30-day readmission scores apply to heart failure (HF) and MI and eliminate the temporal cost to busy physicians by calculating the scores within an electronic health record (EHR) and providing them electronically, empowering clinicians to practise score-driven medicine without the burden of finding data or entering them into a risk calculator.22
This study sought to use the principles underlying IMRS to create usable but highly predictive Intermountain Major Adverse Cardiovascular Events (IMACE) risk scores for patients requiring long-term post-MI care, with the following purposes: (1) to predict MACE occurring 1–3 years after MI and (2) to be feasible for use in routine clinical care.
The goal of this study was to develop a risk prediction tool to stratify long-term (>1 year) post-MI risk. This risk tool would empower the determination of which patients have a favorable comparison of benefits vs. risks to whom enhanced MACE-ameliorating evaluation and care (including up to 36 months of DAPT) could be applied. Post-MI patients (aged ≥18 years) admitted to Intermountain Medical Center, LDS Hospital or Cottonwood Hospital (Salt Lake City, Utah, USA) for an acute MI (defined as an MI by the attending cardiologist based on symptoms, ECG, and cardiac enzyme and troponin levels) during the era of P2Y12 inhibitors (January 1997 through April 2014) were included as subjects. Subjects were followed until 30 April 2015. In keeping with the criteria of the Prevention of Cardiovascular Events in Patients with Prior Heart Attack Using Ticagrelor Compared to Placebo on a Background of Aspirin-Thrombolysis in Myocardial Infarction 54 (PEGASUS-TIMI 54) trial, inclusion was limited to patients with acute MI who survived ≥1 year and had ≥1 risk factor for recurrent MI (ie, age ≥65 years, prior MI before baseline, ≥2 vessel coronary artery disease, diabetes diagnosis or an estimated glomerular filtration rate <60 mL/min/1.73 m2). Receipt of a P2Y12 inhibitor prescription was not required.
IMACE models were derived in 70% of subjects (women, n=1342; men, n=3047). For validation, IMACE equations were applied to the 30% who were held aside (women, n=576; men, n=1290). Assignment to derivation and validation used a pseudo-random number generator seeded by a long-period Mersenne Twister. A second validation population was composed of subjects free of MI at baseline and seen at the same hospitals who received percutaneous coronary intervention (PCI) with stent placement (women, n=1779; men, n=4977) during the same time frame.
Risk score derivation
IMACE was created based on the principles used previously for derivation of IMRS,21 which has been validated repeatedly.21 22 Its derivation used variables that are almost always available for most patients, making it possible to compute, adding little to no incremental financial expense, and empowering the delivery of IMACE to clinicians via the EHR such that clinical workflow is not interrupted. IMACE was formed from sex-specific multivariable Cox regression models using age, each factor in the complete blood count (CBC) and the comprehensive metabolic profile, troponin I, B-type natriuretic peptide (BNP), haemoglobin A1c (HbA1c), and the standard lipid panel. For further information about why only these variables were used, please see online supplementary methods regarding enhanced IMACE (e-IMACE) risk scores.
Supplementary file 1
Relative weightings of each factor were determined from regression β-coefficients of multivariable Cox models that were evaluated for the composite primary MACE endpoint of cardiovascular (CV) mortality, MI or stroke. Inclusion of variables in Cox modelling used stepwise selection and forced entry based on results for non-missing data (among derivation patients, about 3% were missing CBC data, 2% basic metabolic profile parameters and 38% the additional six comprehensive metabolic profile factors, 20% lipid panel data, 68% BNP, 72% HbA1c, and 35% troponin values). Mortality was determined from hospital records, Utah death certificates (including cause of death) and the national Social Security Death Master File, providing essentially complete follow-up for mortality events.23 Incident MI and stroke were determined based on International Classification of Diseases coding from Intermountain Healthcare electronic records of 22 hospitals throughout Utah and southern Idaho (Intermountain is the healthcare provider for approximately two-thirds of the population in these areas, and the vast majority of patients return to Intermountain facilities for subsequent care) that capture more than 90% of non-fatal events (lost to follow-up for non-fatal CV events is <10% historically).23 In Cox regression, analyses were landmarked at 1 year, and censored at 3 years after the index MI for non-event subjects or at the time of mortality for deaths due to non-CV causes. Primary analyses used survival time until the first event for subjects with more than one MACE event.
Scalar variable weightings were applied to each independent variable to calculate IMACE using previously described methodology.21 Laboratory variables were categorised into quintiles plus a category for subjects with missing data. Age was categorised into decades. Some laboratory factors were excluded due to collinearity of variables: only one of the red blood cell count, haemoglobin and haematocrit triad was used, and one each from the following pairs: mean corpuscular volume/mean corpuscular haemoglobin, sodium/chloride, creatinine/blood urea nitrogen (BUN), and total cholesterol/low-density lipoprotein cholesterol. The Hosmer-Lemeshow statistic was used to evaluate model calibration.
Sex-specific thresholds for IMACE were established as follows: low risk was defined as IMACE values less than the threshold where sensitivity was ≥90% in the derivation population, high risk as scores equal to or above the threshold where specificity was ≥90%, and mild and moderate risk as scores between those low-risk and high-risk thresholds, but further divided at the location where derivation results indicated a substantial increase in MACE risk from mild to moderate.
Risk score validation
Scalar scores for variables contributing to IMACE were applied to post-MI validation population to test whether IMACE applies outside of the derivation population. Cox regression was used with IMACE as the independent and MACE as the dependent variables. Both continuous and categorised IMACE variables were evaluated. HR, 95% CI and P values are reported, along with percentages of events in IMACE risk categories.
Kaplan-Meier survival curves graphically examined the event-free survival by risk category and provided log-rank statistics. A P value of ≤0.05 was designated as significant for sex-specific associations of IMACE with MACE in the validation population (ie, the derivation results were considered hypothesis-generating). Secondary endpoints were the individual events of MI, stroke, CV mortality and all-cause mortality, and required Bonferroni correction for five tests of hypothesis (these four plus the primary hypothesis) at P≤0.01 for significance, with P >0.01 but ≤0.05 requiring further validation in future studies. Adjustment of IMACE associations was performed in Cox regression to evaluate the effect of adjusting for clinical variables not included as risk score components; further, e-IMACE was also created by adding other independent predictors such as HF diagnosis, atrial fibrillation history and medications (online supplementary methods).
Further statistical considerations
To compare the predictive ability of IMACE, the 6-month Global Registry of Acute Coronary Events (GRACE) score,9–111-year IMRS21 and DAPT Study (DAPT-S) risk score17 were calculated in the derivation and validation populations using previously published methodologies. Each risk score was compared with IMACE using measures of discrimination and improvement in risk stratification (ie, the c-statistics and net reclassification improvement (NRI) index) to determine to what extent IMACE provides additional value.24–27 Statistical significance for the NRI was assessed using an asymptotic test of a null hypothesis that NRI=0.25–27 For IMACE, the low-risk and mild-risk groups were combined in assessments of the NRI, while IMRS21 and GRACE28 categorisation used three standard groupings, and tertiles were used for DAPT-S scores. Both the DAPT-S full score that computes risk of events corrected for bleeding risk (by reducing the score for greater age) and the DAPT-S events-only score (excluding age) were calculated.17 The sensitivity, specificity, positive predictive value and negative predictive value (NPV) of each score were calculated for the high-risk group compared with all other groups.
Baseline characteristics are reported as percentages, means with SD, or medians with IQR, as appropriate. Comparisons of baseline variables between IMACE categories used the χ2 test, Student’s t-test or analysis of variance, as appropriate. Continuous variables whose distribution violated the normality assumption were transformed using the natural log prior to analysis.
Variables evaluated for inclusion in IMACE are presented in online supplementary table S1 stratified by IMACE risk categories of low, mild, moderate and high risk: scores of ≤2, 3–4, 5–6 and ≥7 for women, and ≤2, 3–5, 6–7 and ≥8 for men, respectively. The maximum possible IMACE value was 12 for women but actually ranged from 0 to 11, and was 16 for men with a range of 0–14. Online supplementary table S2 presents the results for the variables used to create e-IMACE, including comorbidities and treatments. In derivation and validation groups, PCI was the mode of treatment for 58.7% and 55.4% of women and 66.2% and 67.3% of men, respectively. The average follow-up among non-event subjects was 2.9±0.4 years in women and in men for both derivation and validation populations.
Among women, significant predictors of MACE were age, glucose, creatinine, haemoglobin, platelet count, red cell distribution width (RDW) and BNP. In men, MACE predictors were age, potassium, glucose, BUN, haematocrit, white blood cell count, RDW, mean platelet volume and BNP. Model calibration was good, with Hosmer-Lemeshow goodness of fit having P=0.79 for women and P=0.55 for men. Online supplementary table 3 displays the variable thresholds for each factor used to compose IMACE and the sample sizes of non-missing values, and reveals the substantial differences between predictor variables for women and men.
Survival curves for the association of IMACE with MACE demonstrated significant stratification in the validation groups for both women (figure 1A) and men (figure 1B). In women, comparisons between risk categories had log-rank P<0.001 for IMACE moderate-risk versus low-risk and high-risk versus low-risk categories; men had P<0.001 for both moderate-risk and high-risk versus low-risk categories (HRs and CIs are provided in online supplementary table 4). Among women, high-risk versus moderate-risk categories had log-rank P=0.008, with a trend towards significance for that comparison in men (P=0.054 for high-risk vs moderate-risk categories). The proportion of validation patients in each IMACE risk category for the composite outcome and individually for MI, stroke and CV mortality plus all-cause mortality is shown in table 1. When modelled as a continuous variable, IMACE had HR=1.44 per +1 score in women (95% CI 1.29 to 1.61; P<0.001) and HR=1.29 per +1 score in men (95% CI 1.20 to 1.38; P<0.001). HRs with 95% CI are shown in figure 2 for the validation population, both for IMACE alone and for models adjusting for factors from online supplementary table 2). Findings for the second validation population are shown in the online supplementary results. Association of IMACE with the composite outcome in ST elevation MI versus non-ST elevation MI did not differ, with formal test of interaction showing P interaction=0.19 in validation women and P interaction=0.45 in validation men.
IMACE and other risk score c-statistics in the derivation and validation populations are shown in table 2 (see online supplementary figures S1–S3 for the IMRS, GRACE and DAPT-S Kaplan-Meier survival curves). For women, IMACE had 90% specificity and 96% NPV, while for men it had 90% specificity and 95% NPV. Other measures of predictive value and those of the IMRS, GRACE and DAPT-S risk scores are found in online supplementary table S5. For IMACE, accuracy in the derivation and primary validation groups was 82.1% and 84.2% for women and 84.1% and 85.0% for men, respectively. IMACE significantly reclassified risk stratification of MACE for women and men compared with IMRS, GRACE and DAPT-S risk scores (table 2).
Sex-specific risk scores were created in stable post-MI patients (≥1 year free from CV events) with ≥1 risk factor for recurrent MI. IMACE equations were validated to predict MACE in an independent set of post-MI patients and a second validation of non-MI stent recipients. These dynamic EHR-delivered scores performed well compared with other risk scores. IMACE differentiates itself also in being an easily implemented EHR-based tool that summarises risk using common tests that post-MI patients receive routinely and serially. Developed in part based on previously implemented clinical risk stratification methods,22 IMACE provides additional information to clinicians to facilitate the best use of their expertise in clinical decision making without adding burden to their workflow.
Various studies have evaluated extended DAPT beyond 1 year of therapy. DAPT-S showed that, compared with aspirin alone, clopidogrel (75 mg/day) or prasugrel (10 mg/day) with aspirin reduced stent thrombosis risk by 71% and a composite of all-cause mortality, MI or stroke by 29% in subjects receiving drug-eluting stents.13 Risk of bleeding was increased by DAPT.13 The Optimal Dual Antiplatelet Therapy trial found a non-significant (P=0.17) but similar magnitude 25% reduction in risk of events for long-term DAPT, with no difference in bleeding risk.14
The PEGASUS-TIMI 54 trial provided the evidence base for expansion of the Food and Drug Administration-approved ticagrelor indication to include a twice-daily 60 mg dose (with 75–150 mg of aspirin daily) beyond 1 year post-MI.12 This trial included patients with an MI 1–3 years prior to enrolment and found significant reduction in MACE for long-term DAPT, at the cost of higher rates of major bleeding.12 Given this and other evidence of a modest but important benefit from long-term DAPT above the mild increase in bleeding risk,15 revised guidelines for >1 year of DAPT were released in 2016.16 In particular, they reiterate the need to balance DAPT’s long-term CV benefit with the risk of major bleeding.
To quantify the benefits and costs of long-term DAPT, the DAPT-S risk score was created to identify patients whose potential DAPT benefit is greater than its bleeding risk.17 In contrast to the IMACE derivation population, the DAPT-S score was created among subjects who received PCI with drug-eluting stents and explicitly models both benefits and risks of DAPT. In a validation population, the DAPT-S score discrimination of risk for MI or stent thrombosis was modest (c=0.64).17 In the present study’s post-MI validation, it was poorly predictive of MACE (c<0.60; table 2), but this was improved by removing the DAPT-S score’s age-based correction for bleeding risk (men: c=0.60, women: c=0.65). In the same post-MI validation, IMACE had c=0.734 for women (+0.085 to +0.160 greater than the DAPT-S scores) and c=0.672 for men (+0.044 to +0.121 greater than the DAPT-S score), with IMACE NRIs well above the de facto standard of 0.10.
IMACE, like IMRS, was created for broad application to all types of enhanced evaluation and care, including by harnessing the risk information from emerging clinical risk factors like the RDW.21 22 29 30 IMACE was derived in patients who were not on long-term DAPT; consequently, it does not explicitly model bleeding risk. IMACE can be applied to weighing the risks and benefits of long-term DAPT given its ability to stratify risk of MACE, but also to the costs and benefits of other clinical decisions that may be considered in a patient’s long-term care plan, such as greater frequency of follow-up visits, application of remote health monitoring and electronic nudging to increase compliance with medications and health interventions, and more intense long-term dosing of other post-MI medications.
Altering a patient care plan for reduction of CV events, however, only begins with a precision risk score in a learning healthcare system. The optimal approach combines the risk score with implementation science, including a method of delivery that does not require a physician or staff member’s involvement to obtain the risk information and a care process model (CPM) that is attuned to a patient’s risk stratum. Previously, a successful marriage was performed of usable EHR-based risk scores (ie, IMRS and an Intermountain 30-day readmission risk score for patients with HF) to a CPM derived for patients with higher risk HF.22 Inpatients with HF with a high risk of 30-day readmission or of 30-day mortality were managed according to the high-intensity CPM (lower risk patients received standard care). Pilot results showed improved 30-day mortality (7% after implementation vs 19% preimplementation), with no increase in readmission (12.8% after vs 14.1% preimplementation (with denominators corrected for those alive and at risk of readmission)).22 Development of the CPM was performed by clinicians using published evidence and best care practices.22
Using a similar methodology to the HF programme,22 IMACE can be used in conjunction with IMRS or other scores that predict risk during the first post-MI year and are either fully or partially integrated on the back end of an EHR, where the data processing is outside of the clinical environment but can be then delivered seamlessly to clinicians. IMACE would be used among patients likely to survive that first year to identify who will remain or transition to higher risk in years 1–3 (ie, IMACE score in the moderate or high group) and for whom long-term (>1 year) therapies may provide greater absolute benefit than harm (patients with low or mild risk IMACE scores would likely derive less absolute MACE-reducing benefit and should receive standard care). That is, a higher risk IMACE score could trigger additional evaluations or services at index hospitalisation and during the first year after MI since this long-term post-MI score reveals higher risk that should be ameliorated as soon as possible. Serial IMACE testing (eg, at 1 year post-MI) may also aid in patient evaluation and planning for post-MI years 1–3 because the common laboratory tests underlying IMACE are dynamic, unlike comorbidities or histories of adverse events or procedures, and thus IMACE when measured serially may provide important information about increasing or decreasing risk over time.29
For example, a baseline evaluation of and planning for impacting a patient’s social risk factors (eg, schedule routine patient contact or provide methods for the patient to access an advanced practice clinician rather than just going to the emergency department) would be valuable given the potential continuing adverse impact of these over the long term. A baseline education session by a pharmacist to better instruct the patient about their reconciled list of medications or by an advanced practice clinician regarding other health education needs prior to discharge could be considered and more frequent (than standard) clinic visits during the first year after MI could be valuable. Deployment of higher cost remote health monitoring and electronic nudging to increase compliance with medications and health interventions may be justified for those who will survive for >1 year post-MI but will have potentially more morbidities and health system encounters as indicated by higher IMACE score. Also, a higher IMACE could motivate clinicians to aggressively ensure that the patient attends their 1-year post-MI clinic visit to reassess risk scores (previously we showed that yearly recalculation of IMRS provided added prognostic information29) and biomarkers, and to re-evaluate patient health education, medication needs and compliance, and success of risk factor reduction.
The DAPT-S risk score,17 GRACE model9–11 and IMRS21 were compared with IMACE herein (see online supplementary dicussion), but were created to predict different endpoints or for distinct follow-up periods. Study comparisons may consequently be imperfect and interpretation of results should be performed with caution. At least one implication that is clear, though, is that risk scores created for prognostication during the first year after MI are not as applicable to post-MI years 1–3 as scores created specifically for that long-term period.
Inclusion of subjects for derivation and validation of IMACE was limited to patients who were event-free for 1 year after index MI. Consequently, the setting for application of IMACE is among post-MI patients, although application to non-MI PCI patients herein showed substantial value for applying IMACE. Usefully, IMACE was composed using data available at baseline; thus, modification of a patient care plan based on IMACE can be performed with a perspective on long-term patient health. Further, while the Intermountain data registry used to derive and validate IMACE compares favourably with other international CV populations such as the Reduction of Atherothrombosis for Continued Health and the Study of Heart and Renal Protection registries,31 32 some factors such as smoking and race are substantially different. Subjects included here were from one geographical region, and additional replications of IMACE in geographically and racially distinct populations are required; however, because IMRS was created in a population from the same region and replicated in populations from various locations across the world,21 30 33 it is anticipated that IMACE will replicate in other populations as well.
While ST-elevation MI (STEMI) data were included in e-IMACE, other ECG data were unavailable in this study but would be valuable to consider in future evaluations of IMACE. Due to the passive surveillance methods for capturing non-fatal health events, findings herein may be biased by misspecification of MACE status, although this should only bias towards the null hypothesis. Finally, IMACE does not explicitly account for the bleeding risk from DAPT; therefore, any CPM created to guide care should explicitly account for this.
IMACE sex-specific risk scores effectively stratified 1-year to 3-year post-MI risk of CV events with modest (men) to excellent (women) risk discrimination. IMACE is an inexpensive, efficient, dynamic EHR-delivered tool that empowers the evaluation of which post-MI patients may be the best candidates for more aggressive long-term therapeutic management.
Contributors BDH, JBM and DLL designed the study. BDH, JBM, DB, GLH, NDK and TLB planned the execution of the study. BDH and TLB collected and prepared the data. BDH conducted statistical analyses. BDH, JBM, DB, GLH, NDK and DLL evaluated and interpreted the results. BDH and DLL drafted the manuscript, and JBM, DB, GLH, NDK, TLB and DLL revised the paper for important intellectual content.
Funding This work was supported by an unrestricted educational grant from AstraZeneca to BDH, JBM and DLL, and these Intermountain investigators had final authority over manuscript content.
Competing interests BDH is an inventor of IMRS and two readmission risk scores that are licensed to CareCentra for integration into their products. BDH is PI of grants funded by Intermountain Healthcare’s Foundry Innovation Program, the Intermountain Research and Medical Foundation, CareCentra, AstraZeneca and GlaxoSmithKline for the development and/or implementation of clinical risk scores. DB, GLH and NDK are employees of AstraZeneca.
Ethics approval This study was approved by the Intermountain Healthcare Institutional Review Board and was determined to be a minimal risk, general data-only study that posed minimal risk to subjects and for which a waiver of consent was provided.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement While the data used in this study are HIPAA-protected private information, with the execution of an appropriate institutional contractual agreement, the authors are pleased to collaborate with other scientists on further evaluations of the study data or the larger source data sets available through Intermountain Healthcare.