Automated Analysis vs. Expert Reading in Nuclear Cardiology: Correlations with the Angiographic Score
Single photon emission computed tomography (SPECT) myocardial perfusion imaging (MPI) represents one of the most widely used imaging modalities for the non-invasive investigation of coronary artery disease (CAD) . In general, myocardial SPECT results are based on the visual evaluation of the reconstructed images. The reader should also take into account several additional factors, such as the pre-test likelihood of disease, image quality and the potential presence of artefacts. An interesting factor of SPECT MPI is the development of standardised methods for automated quantitation. Currently, automated analysis of three-dimensional myocardial SPECT data is a constant component of practice in nuclear cardiology. To assist clinical decision making, commercially available software packages can provide myocardial perfusion maps and estimate global and segmental measures of stress/rest perfusion. Regional perfusion scores (17-segment model) can be derived using the average defect severity in a given segment . Software packages assign severity scores to segments, according to a 5-point scale . Segmental scores can be summed either per region, or for the whole myocardium [summed stress score (SSS), summed rest score (SRS), summed difference score (SDS)]. Software packages have been developed using normal databases [ , , ]. However, evidence regarding the association between automated quantitation and expert reading (ER) is scarce [ , , ]. Therefore, the actual incremental clinical diagnostic performance of automated SPECT MPI analysis remains rather unclear. In the present study, we aimed to investigate the correlation between SSS, SRS and SDS values, derived by three widely available software packages [Emory Cardiac Toolbox (ECTb), Myovation (MYO), Quantitative Perfusion SPECT (QPS)], with the reader scoring of these parameters. Subsequently, we assessed the associations between the SSS, SRS, and SDS values, as recorded by automated analyses and reader scoring, with angiographic score, using coronary angiography as the gold standard. 

2.1. Study Population The present study was conducted at the Nuclear Medicine Laboratory of the University Hospital of Larissa, Larissa, Greece. Over the 1-year study period (January 2019–January 2020), 1063 patients with known or suspected (medium or high risk stratification) CAD were referred to our laboratory for stress/rest myocardial perfusion SPECT. Out of the patients who underwent coronary angiography prior to or after SPECT MPI (within a 3-month period), 306 consecutive patients were enrolled in the study, as they did not meet any of the exclusion criteria. Before testing, a brief structured interview took place, and data regarding clinical features, medications, previous cardiac events, CAD risk factors, and cardiac or noncardiac comorbidities were collected. Hypertension was defined as a systolic blood pressure of 140 mmHg or greater at rest and/or a diastolic blood pressure of 90 mmHg or greater at rest, or treatment with antihypertensive agents. Furthermore, patients with diabetes mellitus or lipid disorders were defined according to the interview data, including the use of the corresponding medications. Obesity was considered as a condition with body mass index (BMI) of 30.0 or greater (BMI calculated as weight in kilograms divided by height in metres squared). Due to their potential effects on performance during stress testing, MPI, and the associated parameters, cardio-active medications (i.e., b-blockers, calcium channel antagonists, and nitrates) were temporarily withdrawn for approximately five half-lives . Patients without proper withdrawal of cardio-active medications were not included into the study. Moreover, we excluded patients with severe congenital or valvular heart disorder, as well as those with non-ischemic cardiomyopathy. Other exclusion criteria were previous cardiac invasive procedure [percutaneous coronary intervention (PCI) or coronary artery bypass grafting (CABG)], or a history or other evidence of myocardial infarction. Patients with a previous myocardial infarction comprise an heterogeneous group whose myocardial perfusion study is influenced not only by myocardial ischemia, but also by necrosis linked to both episode severity and applied therapy. Finally, patients with qualitatively suboptimal scintigrams, due to artefacts, were excluded from the study. There was a contraindication or inability to achieve a satisfactory exercise level in 183 patients. In these patients, either adenosine or regadenoson pharmacologic testing, combined with low-level exercise, was performed. Moreover, pharmacologic testing without any type of exercise was performed in nine patients with left bundle branch block (LBBB) or an implantable pacemaker. All participants gave informed consent for their complete participation in the study, according to the Hospital Ethics Committee guidelines and in compliance with the ethical guidelines of the Declaration of Helsinki. Moreover, all participants were given written information concerning the appropriate radiation protection measures. 2.2. Stress Testing Patients underwent symptom-limited treadmill tests according to the Bruce Protocol, after cardio-active medication withdrawal, 6 h- to 12 h-fasting and avoiding smoking or engaging in intense physical activity for at least 3 h before the examination, as previously described . Data on symptoms related to the performance of exercise testing, and estimated workload in metabolic equivalents (METs, using standard tables) were recorded. Adenosine or regadenoson pharmacologic testing, combined with low-level exercise or not, was performed based on the European Association of Nuclear Medicine (EANM) procedural guidelines . 2.3. Coronary Angiography and Angiographic Score Cardiac catheterization studies were requested by cardiologists based on patients’ data. All coronary angiographies were blindly interpreted by one experienced observer. Each stenosis of the vessel lumen greater than 50% was considered hemodynamically significant, while the presence of a stenosis in the left mainstem was considered equivalent to a two-vessel disease. Therefore, each angiographic study was scored according to a 4-point scoring system (angiographic score; O: normal study, 1: one-vessel disease, 2: two-vessel disease, 3: three-vessel disease). 2.4. SPECT Myocardial Perfusion Imaging and Semi-Quantification Acquisition and processing protocols were in accordance with the EANM/European Society of Cardiology (EANM/ESC) procedural guidelines . The range of injected activities was 250–400 MBq for stress acquisitions and 625–1000 MBq for rest acquisitions. Studies were performed using technetium 99m (99mTc) tetrofosmin (Myoview, GE Healthcare, Chicago, IL, USA). All acquisitions were carried out in the supine position, without attenuation-scatter correction, using a dual-headed SPECT camera. Polar and three-dimensional mapping were performed (GE Xeleris Software, Chicago, IL, USA), and filtered back projection with the Butterworth Filter was used for tomographic reconstruction. Polar maps, three-dimensional images, and the reconstructed images of both stress and rest studies were blindly evaluated by two independent experienced observers. Left ventricular (LV) myocardium was divided into 17 segments and radiotracer uptake was scored in each of these segments according to a 5-point scoring system (0: normal uptake; 1: mildly decreased uptake; 2: moderately decreased uptake; 3: severely decreased uptake and 4: no uptake) . If counts were reduced in a region and this was attributed to attenuation artefact, the score was 0 . The view of a third observer was requested in eight studies in which discordance between the two observers was detected, and the disagreement was resolved by consensus . SSS and SRS were calculated by adding the scores of each segment in stress and rest studies, and SDS was obtained by subtracting SRS from SSS . The level of inter-rater agreement for SSS, SRS and SDS was significant, with intraclass correlation coefficients (ICCs) ranging from 0.89 to 0.93, ( p < 0.001). Moreover, for each participant, SSS, SRS and SDS values were recorded, as derived by ECTb and QPS software packages. MYO does not provide standardised segmental perfusion scores. Therefore, we converted the average segmental count values (relative to maximum pixel values in the relevant polar plot) to categorical scores according to >70%, 50–69%, 30–49%, 10–29% and <10% thresholds, as previously described . Finally, automated myocardial perfusion measurements were based on the results of the commercially available packages, without any institutional adjustments. 

The present study was conducted at the Nuclear Medicine Laboratory of the University Hospital of Larissa, Larissa, Greece. Over the 1-year study period (January 2019–January 2020), 1063 patients with known or suspected (medium or high risk stratification) CAD were referred to our laboratory for stress/rest myocardial perfusion SPECT. Out of the patients who underwent coronary angiography prior to or after SPECT MPI (within a 3-month period), 306 consecutive patients were enrolled in the study, as they did not meet any of the exclusion criteria. Before testing, a brief structured interview took place, and data regarding clinical features, medications, previous cardiac events, CAD risk factors, and cardiac or noncardiac comorbidities were collected. Hypertension was defined as a systolic blood pressure of 140 mmHg or greater at rest and/or a diastolic blood pressure of 90 mmHg or greater at rest, or treatment with antihypertensive agents. Furthermore, patients with diabetes mellitus or lipid disorders were defined according to the interview data, including the use of the corresponding medications. Obesity was considered as a condition with body mass index (BMI) of 30.0 or greater (BMI calculated as weight in kilograms divided by height in metres squared). Due to their potential effects on performance during stress testing, MPI, and the associated parameters, cardio-active medications (i.e., b-blockers, calcium channel antagonists, and nitrates) were temporarily withdrawn for approximately five half-lives . Patients without proper withdrawal of cardio-active medications were not included into the study. Moreover, we excluded patients with severe congenital or valvular heart disorder, as well as those with non-ischemic cardiomyopathy. Other exclusion criteria were previous cardiac invasive procedure [percutaneous coronary intervention (PCI) or coronary artery bypass grafting (CABG)], or a history or other evidence of myocardial infarction. Patients with a previous myocardial infarction comprise an heterogeneous group whose myocardial perfusion study is influenced not only by myocardial ischemia, but also by necrosis linked to both episode severity and applied therapy. Finally, patients with qualitatively suboptimal scintigrams, due to artefacts, were excluded from the study. There was a contraindication or inability to achieve a satisfactory exercise level in 183 patients. In these patients, either adenosine or regadenoson pharmacologic testing, combined with low-level exercise, was performed. Moreover, pharmacologic testing without any type of exercise was performed in nine patients with left bundle branch block (LBBB) or an implantable pacemaker. All participants gave informed consent for their complete participation in the study, according to the Hospital Ethics Committee guidelines and in compliance with the ethical guidelines of the Declaration of Helsinki. Moreover, all participants were given written information concerning the appropriate radiation protection measures. 

Patients underwent symptom-limited treadmill tests according to the Bruce Protocol, after cardio-active medication withdrawal, 6 h- to 12 h-fasting and avoiding smoking or engaging in intense physical activity for at least 3 h before the examination, as previously described . Data on symptoms related to the performance of exercise testing, and estimated workload in metabolic equivalents (METs, using standard tables) were recorded. Adenosine or regadenoson pharmacologic testing, combined with low-level exercise or not, was performed based on the European Association of Nuclear Medicine (EANM) procedural guidelines . 

Cardiac catheterization studies were requested by cardiologists based on patients’ data. All coronary angiographies were blindly interpreted by one experienced observer. Each stenosis of the vessel lumen greater than 50% was considered hemodynamically significant, while the presence of a stenosis in the left mainstem was considered equivalent to a two-vessel disease. Therefore, each angiographic study was scored according to a 4-point scoring system (angiographic score; O: normal study, 1: one-vessel disease, 2: two-vessel disease, 3: three-vessel disease). 

Acquisition and processing protocols were in accordance with the EANM/European Society of Cardiology (EANM/ESC) procedural guidelines . The range of injected activities was 250–400 MBq for stress acquisitions and 625–1000 MBq for rest acquisitions. Studies were performed using technetium 99m (99mTc) tetrofosmin (Myoview, GE Healthcare, Chicago, IL, USA). All acquisitions were carried out in the supine position, without attenuation-scatter correction, using a dual-headed SPECT camera. Polar and three-dimensional mapping were performed (GE Xeleris Software, Chicago, IL, USA), and filtered back projection with the Butterworth Filter was used for tomographic reconstruction. Polar maps, three-dimensional images, and the reconstructed images of both stress and rest studies were blindly evaluated by two independent experienced observers. Left ventricular (LV) myocardium was divided into 17 segments and radiotracer uptake was scored in each of these segments according to a 5-point scoring system (0: normal uptake; 1: mildly decreased uptake; 2: moderately decreased uptake; 3: severely decreased uptake and 4: no uptake) . If counts were reduced in a region and this was attributed to attenuation artefact, the score was 0 . The view of a third observer was requested in eight studies in which discordance between the two observers was detected, and the disagreement was resolved by consensus . SSS and SRS were calculated by adding the scores of each segment in stress and rest studies, and SDS was obtained by subtracting SRS from SSS . The level of inter-rater agreement for SSS, SRS and SDS was significant, with intraclass correlation coefficients (ICCs) ranging from 0.89 to 0.93, ( p < 0.001). Moreover, for each participant, SSS, SRS and SDS values were recorded, as derived by ECTb and QPS software packages. MYO does not provide standardised segmental perfusion scores. Therefore, we converted the average segmental count values (relative to maximum pixel values in the relevant polar plot) to categorical scores according to >70%, 50–69%, 30–49%, 10–29% and <10% thresholds, as previously described . Finally, automated myocardial perfusion measurements were based on the results of the commercially available packages, without any institutional adjustments. 

Quantitative variables were expressed as mean (standard deviation) or as median (interquantile range). Qualitative variables were expressed as absolute and relative frequencies. Spearman correlation coefficients were used to explore the association of two continuous variables. Correlation coefficients between 0.1–0.29 were considered “poor”, between 0.30–0.49 “fair”, between 0.50–0.69 “moderate” and between 0.70–1.0 “strong”. ICC was used to assess inter-rater agreement, as well as the agreement between ER and the software packages, concerning SSS, SRS and SDS parameters. Following the recommendations given by Koo and Li, ICC below 0.5 indicates poor agreement, between 0.5 and 0.75 indicates moderate agreement, between 0.75 and 0.90 indicates good agreement, and above 0.9 indicates excellent agreement . Agreement between the expert and the software packages was further assessed by Bland–Altman 95% confidence intervals (CI) for limits of agreement (LOA). The 95% CI for LOA indicates that 95% of the differences fall between these two limits. All reported p values are two-tailed. Statistical significance was set at p < 0.05 and analyses were conducted using SPSS statistical software (version 22.0). 

The sample consisted of 306 patients (62.7% male) with mean age 63.8 years (SD = 9.9 years). Demographic and clinical characteristics of the patients are presented in . The majority of the patients reported having CAD-related symptoms (65.4%). Obesity was recorded in 45.2% of the sample. Smokers were 41.2% of the patients. Arterial hypertension and diabetes mellitus were recorded in 77.8% and 36.6% of the patients, respectively. Median angiographic score was 1 (IQR: 0–2). Mean SSS and mean SRS according to ER were significantly lower than the corresponding SSS and SRS values derived from ECTb, MYO and QPS ( ). Mean SDS based on ER was significantly lower than SDS derived from ECTb and MYO, but significantly higher in comparison to QPS results. ICCs of SSS, SRS and SDS between ER and the software packages are presented in . All ICCs were moderate to excellent, as well as significant (SSS ICCs ranged from 0.82 to 0.91, SRS ICCs ranged from 0.67 to 0.79 and SDS ICCS ranged from 0.69 to 0.82). Bland-Altman plots for SSS, SRS and SDS are presented in . For SSS, limits of agreement between ER and ECTb ranged from −5.07 to 11.82 ( (A1)), between ER and MYO from −4.41 to 11.24 ( (A2)) and between ER and QPS from −5.02 to 4.90 ( (A3)). For SRS, limits of agreement between ER and ECTb ranged from −2.70 to 7.81 ( (B1)), between ER and MYO from −2.54 to 9.31 ( (B2)) and between ER and QPS from −2.83 to 4.35 ( (B3)). For SDS, limits of agreement between ER and ECTb ranged from −6.02 to 7.64 ( (C1)), between ER and MYO from −6.46 to 6.54 ( (C2)) and between ER and QPS from −5.55 to 3.91 ( (C3)). Limits of agreement between ER and software packages were widely indicating that the average discrepancy between expert scoring and software analyses was large enough. SSS values based on ER were significantly correlated with the corresponding results of all software packages ( ) and all coefficients were “strong”. Similarly, SRS and SDS values according to ER were significantly associated with SRS and SDS results from all software packages and correlations for SRS were almost all moderate. presents Spearman’s correlation coefficients of angiographic score with SSS, SRS and SDS based on ER, as well as with the corresponding results derived from software packages. SSS, SRS and SDS values according to either software packages or expert reading were significantly associated with angiographic score. By examining the coefficients, however, expert estimation of SSS, SRS and SDS was more strongly correlated with the angiographic score compared to all three software packages. Specifically, correlation coefficients of the angiographic score with the software packages were fair and less than 0.5, while correlation coefficients of the angiographic score with ER were more than 0.5 and moderate. 

SPECT MPI is a widely used non-invasive imaging modality for the investigation of patients with known or suspected CAD . In general, MPI results are based on the visual evaluation of the reconstructed images. In an effort to restrict the influence of reader’s experience over image interpretation, a number of related algorithms were developed aiming to provide automated measurements of myocardial perfusion. Nowadays, several software packages, such as ECTb, MYO and QPS, are available, and automated analysis of myocardial perfusion has become part of the routine practice in nuclear cardiology. Interestingly, previous studies have demonstrated that software packages can supplement visual evaluation, as well as their high reproducibility [ , , , , ]. However, automated analysis is currently used only as an adjunct to visual interpretation [ , , , ]. The adjunctive role of automated analyses in nuclear cardiology practice is at least partially associated with a main drawback of software packages, which cannot explicitly differentiate between real perfusion abnormalities and artefacts . In addition, although the performance of the software packages has been demonstrated to be similar, certain differences in the magnitudes of the derived values have been reported . In particular, comparing the diagnostic performance of three software packages (ECTb, QPS, and 4DMSPECT), Wolak et al. reported differences in myocardial perfusion quantification, while Knollmann et al. found differences in measurements (derived from QPS and 4DMSPECT), using either manufacturer’s or institutional normal databases . Similarly, Johansson et al. demonstrated considerable differences in performance between ECTb, QPS and 4DMSPECT, and Knollmann et al. noted the influence of heart-axis tilt on automated SSS calculations . Furthermore, studying a large patient population in our laboratory, we reported a favourable concordance between ECTb, MYO and QPS, whereas differences were demonstrated in pair comparisons . Obviously, for each software package, quantitative measurements are algorithm-specific, and the observed differences in the derived values among software packages could be related to the different principles and assumptions on which each algorithm is based. Although several studies exist investigating the cross-correlation of the outputs of different quantitative SPECT MPI algorithms, there is only few published data regarding the association between automated quantitation and ER. In the present study, we aimed to investigate the correlation between automated measurements of SSS, SRS and SDS (using three software packages, as commercially available) with the ER of these parameters. Previously, Arsanjani et al. and Duvall et al. had compared the automated measurements of total perfusion deficit (TPD) with the visual quantitation of stress and rest images, while Driessen at al. had investigated MPI SPECT data, both visually and automatically [SSS, SDS, stress total perfusion deficit (S-TPD), and ischemic total perfusion deficit (I-TPD)] [ , , ]. We demonstrated that the ICCs of SSS, SRS and SDS between ER and the software packages were moderate to excellent, while the average discrepancy between expert scoring and software analyses was large enough (based on the limits of agreement between ER and software packages). Furthermore, although visually defined SSS, SRS and SDS were significantly correlated with the corresponding results of all software packages, expert estimations of SSS, SRS and SDS were more strongly correlated with the angiographic score. On the other hand, Arsanjani et al. reported that automated analysis of attenuation-corrected (AC) and non-corrected (NC) SPECT MPIs was at least equivalent to visual interpretation, when compared to coronary angiography (≥70% luminal stenosis) . Furthermore, Duvall et al. demonstrated that automated analysis and visual quantitation of stress and rest images had similar diagnostic accuracy, in terms of the angiographic detection of ≥70% stenoses, with the most favourable diagnostic accuracy achieved in combined supine and prone stress imaging . Finally, Driessen et al. found that visual quantification slightly outperformed automated analysis in the detection of fractional flow reserve-defined significant CAD, whereas the diagnostic accuracy of automated analysis equalled ER after optimization with an institutional normal database and thresholds . The observed differences between our results and those reported in previous studies could be partially attributed to the differences in study samples, as well as in the definition of the hemodynamically significant stenosis (greater than 50% vs. 70% luminal stenosis). The visual evaluation of myocardial SPECT data continues to play a crucial role in the interpretation of MPI examinations. To the best of our knowledge, this is the first study investigating the correlations between automated analysis and ER of SSS, SRS, and SDS, with regard to the angiographic score. Based on the existing research, ER has demonstrated similar performance in the quantification of myocardial perfusion, in comparison to that derived from automated analyses using different methodologies. According to our results, expert estimations of SSS, SRS and SDS were significantly correlated with the corresponding results of all software packages, but visually defined SSS, SRS and SDS were more strongly correlated with the angiographic score, indicating a better performance of ER when compared to automated analysis. However, automated analysis could have an adjunctive role in the interpretation of SPECT MPI examinations, particularly in patients with ambiguous scintigraphic findings. In addition, automated measurements may help beginners, or less experienced physicians, in their practice. 

In the present study, comparing the performance of automated analysis to expert reading, we demonstrated that visually defined SSS, SRS and SDS were more strongly correlated to angiographic score than the software-derived corresponding values. Our results support the significance of ER for the interpretation of myocardial SPECT data. Measurements derived from software packages can be useful mainly in cases of ambiguous findings, supplementary to visual interpretation. 

