Methods
Study population
We conducted a retrospective analysis of patients from the Duke Cardiovascular Disease and Nuclear Cardiology Databanks. All patients are followed at 6 months, 1 year and annually thereafter, with recording of major clinical events.11 We identified 278 consecutive patients with angiographically documented CAD and LVSD (ejection fraction (EF) ≤40%), who underwent two serial MPS between 1993 and 2009 (figure 1). All patients had CAD documented by coronary angiography within 180 days of initial MPS.
Study population (CAD, coronary artery disease; LVSD, left ventricular systolic dysfunction; MPS, myocardial perfusion scan; MT, medical therapy; Revasc, revascularisation).
Stress testing and imaging
Patients capable of exercise underwent treadmill stress testing with the Bruce protocol, unless the physician requested an alternative protocol. Patients unable to exercise underwent pharmacological stress testing. MPS was performed according to described protocols.12 ,13 In each patient, the same stress modality (exercise or pharmacological) was performed on the initial and follow-up MPS. Tc-99m was the radiotracer used in all studies and no patients received nitroglycerin during the MPS. The studies were independently interpreted without attenuation correction by three nuclear cardiologists blinded to treatment group.
The MPS were evaluated semiquantitatively for severity and extent of abnormalities with relative perfusion recorded in each myocardial segment (0=no defect, 1=mild defect, 2=moderate defect and 3=severe defect) at rest and stress. The summed difference score (SDS), which is the sum of the differences between the stress and rest perfusion scores (reversible defects), was determined. At the time these data were collected, we used a 12-segment model. We used a previously reported algorithm for conversion of 12-segment perfusion scores to 17-segment scores, which is highly correlated with expert reading of the same studies by the 17-segment model.14 The percentage LV ischaemia for each scan was calculated from the SDS.14
Follow-up and outcomes
Treatment group was assigned on an intention-to-treat basis. Patients were included in the revascularisation group if percutaneous coronary intervention (PCI) or coronary artery bypass grafting (CABG) occurred within 60 days of the first scan. If no revascularisation occurred within this period, they were assigned to the MT group. Follow-up time was initiated at the second MPS (time=0). For all survival analyses, we censored data at the time of participant loss to follow-up, or the end of study event surveillance follow-up. The primary end point was a composite of death/MI. An independent, blinded clinical events committee reviewed major clinical events.
Statistical analysis
Differences in baseline characteristics were compared with the use of the t test or χ2 statistics. When appropriate, paired t tests were used to compare paired data. Continuous variables that were not distributed normally were compared with a Wilcoxon rank-sum test for independent groups and a Wilcoxon signed-rank test for paired data. A threshold change in ischaemia of ≥5% was used, as it represents a change that exceeds test repeatability.15
Kaplan-Meier methods were used to evaluate time to death/MI. Cox proportional hazards regression modelling was used to identify factors that were associated with death/MI. After examining the results of a flexible Cox model-fitting approach involving cubic polynomial spline functions,16 the linearity of the unadjusted relationship between each continuous and ordinal variable and death/MI was assessed. The proportionality assumption was verified using Schoenfeld residuals. For the multivariable model, covariates were chosen based on known clinical risk factors as well as by stepwise selection (and backwards elimination) at p<0.05 from the list of baseline characteristics. Candidate variables for adjustment included age, non-cardiac Charlson index, LVEF, diabetes, smoking, number of diseased vessels, New York Heart Association (NYHA) class, race, history of MI, peripheral vascular disease, cerebrovascular disease and hypertension. To explore potential modifying effects of baseline ischaemia, the interaction between ≥5% ischaemia improvement or worsening and ischaemia on the first scan was tested in the adjusted Cox model. An inverse probability weighting model adjustment to account for propensity for treatment was performed to investigate the association between treatment and outcome. The association between ischaemia change and outcome was also evaluated in this manner. All tests were two-tailed, and statistical significance was declared if the two-sided p value was <0.05. Most variables had very low rates of missingness (ie, <5%). For variables with <5% missingness, we imputed continuous variables to the overall median value, dichotomous variables to ‘no’, and multichotomous variables to the most frequent categorical value. For variables with >5% missingness, we treated the missing values as a separate category. Statistical analyses were performed using SASV.9.2 (SAS Institute, Cary, North Carolina, USA). The study protocol was reviewed and approved by the Duke Institutional Review Board and all participants provided informed consent.