Circulation Reports
Online ISSN : 2434-0790
Imaging
Machine Learning for Evaluating Vulnerable Plaque on Coronary Computed Tomography Using Spectral Imaging
Junji MochizukiYoshiki HataTakeshi Nakaura Katsushi HashimotoHiroyuki UetaniYasunori NagayamaMasafumi KidohYoshinori FunamaToshinori Hirai
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2024 Volume 6 Issue 12 Pages 564-572

Details
Abstract

Background: This study aimed to determine whether spectral imaging with dual-energy computed tomography (CT) can improve diagnostic performance for coronary plaque characterization.

Methods and Results: We conducted a retrospective analysis of 30 patients with coronary plaques, using coronary CT angiography (dual-layer CT) and intravascular ultrasound (IVUS) studies. Based on IVUS findings, patients were diagnosed with either vulnerable or stable plaques. We computed 7 histogram parameters for plaque CT numbers in 120 kVp images and virtual monochromatic images ranging from 40 to 140 keV at 5-keV intervals. A predictive model was developed using histogram data of optimal energy, plaque volume or stenosis, and a combination of both. The model’s performance was evaluated by calculating the area under the receiver operating characteristic curve (AUC) using 5-fold cross-validation. Peak diagnostic performances for each histogram parameter were observed at various energy levels (40–110 keV) in the univariate logistic regression model. The histogram model demonstrated optimal diagnostic performance at 65 keV, with an AUC of 0.81. The combined model, incorporating histogram data and plaque volume, achieved an AUC of 0.85, which was similar to the performance of qualitative CT characteristics (AUC=0.88; P=0.70).

Conclusions: Spectral imaging with dual-energy CT can enhance the diagnostic performance of machine learning using CT histograms for coronary plaque characterization.

Identifying vulnerable plaque in coronary arteries, a contributory factor for myocardial infarction, is crucial. Conventionally, coronary plaques can only be evaluated using invasive methodologies such as intravascular ultrasound (IVUS). Cardiac computed tomography (CT) angiography has emerged as a promising non-invasive technique for detecting vulnerable plaques prior to a subsequent coronary event.17 Characteristics of vulnerable plaque include low attenuation areas, positive remodeling, spotty calcification, ring-like sign, and plaque ulceration with intraplaque contrast penetration.811 The mean CT score is useful for characterizing coronary plaques. Motoyama et al. classified non-calcified plaques into lipid-rich plaque (cut-off <30 HU) and fibrous plaque (30–150 HU).8 However, overlapping lipid and fibrous plaque attenuation posed a limitation.10,1214

The efficacy of CT analysis in plaque diagnosis has been reported to be superior with dual-energy techniques (DECT).1520 Ohta et al. demonstrated that the virtual monochromatic image (VMI) obtained from DECT optimizes image quality at distinct energy levels for various components of coronary plaque.21 It has been suggested that comparing CT numbers at low and high energy settings may facilitate the differential diagnosis of lipid-rich plaque and other plaque types.19 Furthermore, a histogram analysis using machine learning has been reported to surpass conventional cut-off methods in coronary plaque diagnosis using plaque CT numbers in cardiac CT images.22,23 Although these 2 techniques may hold significant diagnostic potential, no current reports exist on their combined usage, and the optimal X-ray energy for plaque evaluation remains undetermined.

Our hypothesis posits that the precision of plaque characterization can be enhanced by machine learning 3-dimensional plaque histogram evaluation using VMIs.

The aim of this study was to assess the diagnostic performance of machine learning on variable plaques based on 3-dimensional histograms of VMI obtained from spectral CT and to investigate the optimal energy of VMI.

Methods

This retrospective study received institutional review board approval; written informed consent was waived. The study included patients who underwent cardiac CT scans in our hospital from January 2017 to March 2019 and had clinical records of severe stenosis. All procedures performed in this study were in accordance with the ethical standards of the institutional research committee and with the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards.

Cardiac CT Image Acquisition

All patients were scanned with a dual-layer spectral detector CT scanner (iQon Spectral CT; Philips Healthcare). Patients received sublingual nitroglycerin. Standard coronary CT angiography was performed using a 13-s intravenous infusion of 18.0 mgI/mL/s (Iopamiron 370; Bayer HealthCare). Contrast medium delivery was followed by flushing with 30 mL of physiological saline at the same injection rate. The acquisition parameters for cardiac CT imaging were as follows (Table 1): detector collimation 64×0.625 mm; tube rotation time 270 ms; tube voltage 120 kVp; tube current 436.1±103 mA (range 290–611 mA); and volume CT dose index 57.9±18.0 mGy (range 33.0–93.8 mGy).

Table 1.

Scanning Parameters and Contrast Medium Infusion Protocol

Scanning parameter  
Tube voltage (kVp) 120
Tube current (effective mAs) Auto mA
Helical pitch 0.16
Detector collimation (mm) 64×0.625
Rotation time (s) 0.27
Detector collimation (mm) 64×0.625
Rotation time (s) 0.27
Iodine dose (mgI/kg) 240
Bolus tracking trigger (HU) 110
Scan delay (s) 6
Injection duration (s) 13

CT Image Reconstruction

Spectral-based image data were post-processed on a workstation (Spectral Diagnostic Suite; Philips Healthcare) to generate VMIs at 21 distinct energy levels (40–140 keV), with a spectral level of 3 (manufacturer’s recommendation). Conventional CT images reconstructed with IR (iDose level 3; Philips Healthcare) were used as controls. The spectral level and iDose level represent successive approximation reconstructions, and level 3 was chosen based on the manufacturer’s recommendation. The slice thickness/increment of all CT images was 0.8 mm/0.4 mm.

IVUS

The subjects of this study will undergo coronary angiography and FFR to evaluate ischemia. If the FFR is ≤0.80, the diagnosis will be lesion-specific ischemia, and IVUS will then be performed. Target lesions were imaged using IVUS following intracoronary injection of nitroglycerin (1.0 mg). IVUS images were acquired using a 40-MHz IVUS imaging catheter (OptiCross, Boston Scientific Corporation) and an IVUS system (iLAB; Boston Scientific Corporation). Digital copies of the obtained images were preserved for subsequent offline analysis. IVUS analysis was conducted every 1.0 mm, independent of the cardiac cycle, using established software (QIvus 2.1; Medis Medical Imaging Systems, Leiden). The criteria derived from a clinical expert consensus document on guidelines for measurement and evaluation of IVUS from the Japanese Association of Cardiovascular Intervention and Therapeutics were used for grayscale IVUS analysis.24 Specifically, vulnerable plaques were characterized by having at least 2 of the following known attributes of rupture-prone plaques as observed by IVUS: attenuated plaques, positive remodeling, and spotty calcification and a minimal lumen area of <4.0 mm2 by IVUS; or a plaque burden of >70% by IVUS.25,26 The analysis was performed by 2 cardiologists with over 20 years of experience.

Analysis of Coronary CTA Data and Plaque Quantification

Image analysis was conducted by a radiology technologist with 17 years of experience using dedicated cardiac analysis software (Cardiac Viewer and Comprehensive Cardiac Analysis, intelliSpace Portal, Philips Healthcare) and a plaque analysis application for the assessment of non-calcified coronary plaque volume and various characteristics (Figure 1). Semi-automatic lesion identification, characterization, and quantification were performed using a plaque volume analysis program integrated with the dedicated cardiac analysis software.27 The remodeling index (RI) is defined as the ratio of the maximum vessel area to a normal reference vessel area, and plaques are classified as having significant positive remodeling when the RI is >1.1.28 The presence of the napkin-ring sign, characterized by a low-attenuating plaque core surrounded by an area of higher attenuation, was evaluated.29 Spotty calcifications were visually examined as calcifications covering <90° of the lesion circumference and measuring <3 mm in length.30 Low attenuation plaques were defined as CT number >30 HU.31 Societal guidelines using the coronary artery disease reporting and data system (CAD-RADSTM) were used to determine the degree of stenosis: (1) none (0%), or minimal (1–24% stenosis); (2) mild (25–49% stenosis); (3) moderate (50–69% stenosis); (4) severe (70–99% stenosis); and (5) total occlusion (100%).32 Several sequential image processing steps were undertaken, including automatic model-based whole-heart segmentation, automatic extraction of the coronary tree and centerline, automatic luminal and vessel wall contouring, and semi-automatic plaque examination. Plaque volumes were measured and 3-dimensional histograms of plaque extracts at each VMI and 120 kVp were generated. We also calculated 7 histogram parameters (minimum and mean value, standard deviation, maximum value, skewness, kurtosis, and entropy) of the plaque CT number of 120 kVp images and VMI using Python (version 3.8.5; https://www.python.org).

Figure 1.

Plaque analysis. Three-dimensional histogram analysis created using specialized software. Plaque identification and voxel quantification were conducted semi-automatically, and the resultant histogram data was used for analysis.

Training of Machine Learning Models

Due to the limited number of cases, partitioning the training and test groups was infeasible. Consequently, the cross-validation function of scikit-learn was used, which automatically divided the population into 5-fold cross-validation for training and learning purposes. To assess the discriminative ability and optimal energy of each histogram parameter individually for distinguishing between vulnerable and stable plaques, univariate logistic regression was used to train all histogram parameters at all VMIs and 120 kVp, followed by an evaluation. Subsequently, to evaluate the discriminative ability and optimal energy of ML using all histogram parameters of VMI, Random Forest was trained and assessed for all VMIs and 120 kVp. The optimal parameters of Random Forest were determined through the grid search technique.

Furthermore, after determining the optimal energy, we developed a combined model incorporating the histogram parameters at the optimal energy level along with other significant factors. These factors included coronary stenosis and plaque volume measurements that showed statistically significant differences between vulnerable and stable plaques. This combined model aimed to enhance the overall predictive performance by integrating multiple relevant features.

To evaluate the performance of these models, we established a reference standard based on qualitative CT characteristics. This standard was derived from the sum of scores for the presence or absence of previously mentioned features: low attenuation areas, positive remodeling, spotty calcification, and the ring-like sign. This model served as a benchmark to compare the performance of our developed models against the traditional qualitative approach, allowing us to assess potential improvements in plaque vulnerability assessment.

Statistical Analysis

We compared patient’s characteristics and plaque characteristics of vulnerable and stable plaques using the Mann-Whitney U test. The area under the curve (AUC) of our receiver operating characteristic (ROC) analysis of machine learning and conventional CT number-based classifiers was calculated. We compared the diagnostic performance of the machine learning and conventional CT number-based AUCs using the DeLong test. A P value <0.05 was deemed a statistically significant difference. We also calculated the importance of all features. All statistical analyses were performed using Python (version 3.8.5).

Results

Study Population

During the study period, cardiac CT was conducted in 5,679 patients with clinically suspected or confirmed coronary artery disease, following the American College of Cardiology guidelines.33 Of these, 4,451 had no IVUS, 1,129 had no dual-energy data, and 8 were excluded due to motion artifacts and poor image quality. Of the 91 remaining patients, 61 with severe calcified lesions were excluded, and the remaining 30 patients were enrolled in this study. The study group included 25 men and 5 women with a mean age of 68 years (range 48–86 years), a mean body weight of 66.1 kg (range 48.2–99.6 kg), and a mean eGFR of 60.1 mg/dL (range 51.5–69.6; Figure 2). All subjects were diagnosed with chronic coronary artery disease. The disease site was the right coronary artery in 6 patients, left anterior descending artery in 22 patients, and left circumflex artery in 2 patients. Of the total 30 lesions in the study, IVUS analysis classified 17 lesions as vulnerable plaques and 13 lesions as stable plaques. In the comparison between vulnerable and stable plaques, there were no significant differences in the prevalence of hypertension, hyperlipidemia, diabetes, or smoking status (Table 2).

Figure 2.

Flowchart of patient inclusion criteria. CT, computed tomography; IVUS, intravascular ultrasound.

Table 2.

Patient Demographics

  All
(n=30)
Vulnerable plaques
(n=17)
Stable plaques
(n=13)
P value
Age (years) 68.0±10.6 68.8±11.4 66.7±9.2 0.59
Male 24 (80) 13 (76) 11 (85) 0.58
Hypertension 11 (37) 5 (29) 11 (37) 0.35
Hyperlipidaemia 23 (77) 14 (82) 3 (23) 0.40
Diabetes 7 (23) 4 (24) 7 (23) 0.98
Smoking 11 (37) 7 (41) 4 (31) 0.56

Data are presented as n (%), or mean±SD.

Qualitative Assessment of the CT Characteristics

All plaques in both groups exhibited stenosis ≥70%. While plaque volume tended to be larger in vulnerable plaques compared with stable plaques (32.1±17.2 mm3 vs. 21.2±12.0 mm3; P=0.08), this difference did not reach statistical significance. Positive remodeling (RI ≥1.1) was significantly more prevalent in vulnerable plaques compared with stable plaques (65% vs. 15%; P<0.01). The napkin-ring sign was observed more frequently in vulnerable plaques, but this difference was not statistically significant (29% vs. 8%; P=0.13). Spotty calcification showed a trend towards higher prevalence in vulnerable plaques (59% vs. 23%; P=0.05). Similarly, low attenuation plaque (≤30 HU) was more common in vulnerable plaques, although this difference did not reach statistical significance (94% vs. 69%; P=0.07; Table 3).

Table 3.

Comparison of Computed Tomography Characteristics in Vulnerable and Stable Coronary Plaques

  All
(n=30)
Vulnerable plaques
(n=17)
Sable plaques
(n=13)
P value
Stenosis ≥70% 30 (100) 17 (100) 13 (100) 1.00
Plaque volume (mm3) 26.5±10.0 32.1±17.2 21.2±12.0 0.08
Positive remodeling (RI ≥1.1) 13 (43) 11 (65) 2 (15) <0.01
Napkin-ring sign 6 (20) 5 (29) 1 (8) 0.13
Spotty calcification 13 (43) 10 (59) 3 (23) 0.05
Low attenuation plaque (≤30 HU) 25 (83) 16 (94) 9 (69) 0.07

Data are presented as n (%), or mean±SD. RI, remodeling index.

Model Performance Comparison

Using the univariate logistic regression model, the peak diagnostic performance of each histogram parameter was observed at different energy levels (40–140 keV; Figure 3). The AUCs were as follows: mean (0.62; 125–140 keV); minimum (0.73; 120 kVp); maximum (0.71; 40 keV, 45 keV, and 50 keV); standard deviation (0.71; 120 kVp); skewness (0.74; 95 keV); kurtosis (0.61; 95 keV, and 110 keV); and entropy (0.72; 75 keV; Figure 3).

Figure 3.

Area under the curves of each histogram parameter. Area under the curve for each energy level (40–140 keV) in the univariate logistic regression model; however, the optimal virtual monochromatic image energy varied for each histogram parameter.

Using the Random Forest model combining all features, peak diagnostic performance was observed at 65 keV (AUC=0.81; Figure 4A). The AUC of the 65 keV Random Forest model (AUC 0.81) was significantly higher than that of the 120 kVp Random Forest model (AUC 0.72) and the mean CT number (AUC 0.50) using 120 kVp images (P<0.05). Since there was no significant difference in stenosis, a prediction model was created using only plaque volume, resulting in an AUC of 0.73. The combined model, incorporating histogram data and plaque volume, achieved an AUC of 0.85, which was similar to the performance of qualitative CT characteristics (AUC=0.88; P=0.70; Figure 4B). Representative cases are shown in Figure 5.

Figure 4.

Area under the curves (AUC) of various prediction models. AUC of Random Forest at different energy levels (A), and the receiver operating characteristic (ROC) curves (B) of the histogram model at 65 keV (AUC=0.81), plaque volume (AUC=0.73), plaque characteristics (AUC=0.89), and the combined model (AUC=0.85).

Figure 5.

Representative case. Cardiac computed tomography (CT) curved multiplanar image (A) and short-axis image (B) showed stenosis in the left anterior descending artery #6. Plaque histograms (C) at 120 kVp and virtual monochromatic image at different energies showed that the histogram shape changed with energy. Angiography (D) and intravascular ultrasound (E) showed a vulnerable plaque. In this case, the spectral imaging analysis indicated that the plaques were unstable. As a result, distal protection was used to prevent embolization; after percutaneous coronary intervention, there was a transient reduction in blood flow that quickly recovered.

Discussion

The present study suggests that the diagnostic performance of machine learning for coronary vulnerable plaques, utilizing VMI and 3-dimensional histogram parameters, is significantly superior with VMI 65 keV compared with 120 kVp images.

Previous reports have suggested that DECT is a valuable tool for diagnosing coronary plaques. An ex vivo study of 20 human arteries showed that CT number-based characterization of non-calcified plaques using DECT demonstrated selective improvement over conventional CT.19 Furthermore, spectral coronary CTA with low-energy (40–70 keV) post-processing has been reported to enhance the CNR of coronary plaque components.20 At diagnostic energies, Compton scatter and photoelectric absorption, both dependent on the material’s atomic number, are responsible for X-ray attenuation. These interactions vary with the energy provided. Using various X-ray energies, Barreto et al. observed a substantial variation in the attenuation of fibrous plaque but not the lipid-rich necrotic core.15 The precise cause of this phenomenon is challenging to explain, but previous reports have shown that vulnerable plaques are lipid-rich and influenced by blood flow.34 Consequently, it is possible that lower energy did not reveal a difference from stable plaques due to both a reduction in fat absorption and an increase in the contrast effect. If the lipid-rich plaque was minuscule, additional tissue types may have been present in the CT region sample.

In the univariate logistic regression analysis, the optimal VMI energy differed for each histogram parameter; no CT number difference was observed, but the diagnostic performance was superior at lower energies. The minimum CT number demonstrated reasonable results at moderate energies (65–75 keV), while the maximum CT number performed better at lower energies (40–50 keV). This phenomenon may be attributable to the fact that stable plaques have a high fibrous component, causing CT numbers to increase with lower energy and potentially creating the observed difference. The skewness and entropy AUCs exhibited a favorable trend in the intermediate energy range. This outcome may be influenced by the image quality of the VMI, with extremely low energy VMIs increasing noise and possibly not displaying an accurate CT number.

The study findings indicate that the diagnostic efficacy of the histogram analysis approach with machine learning is sensitive to the energy parameters used in dual-energy CT. Histogram analysis can scrutinize the distribution of pixel intensities within an image to appraise fine calcifications and lipids in vulnerable plaques. Given these lesions are minute and characterized by scarce differences in CT values, a highly accurate approach is required to detect them. The 65 keV images, identified as the most suitable energy VMI in this investigation, exhibited superiority over traditional 120 kVp images. In addition to superior spatial and contrast resolution balance, the utilization of VMI benefits, such as beam hardening correction, were believed to have contributed to the positive outcomes obtained.

Another significant finding of this study is that the diagnostic performance of the combined model, which incorporates histogram data at optimal energy and plaque volume (AUC=0.85), is comparable with that of qualitative plaque analysis (AUC=0.88). This observation has important clinical implications. While the model’s value may be limited in centers with extensive experience in cardiac CT evaluation, it offers substantial benefits in other settings. Notably, it enables individuals without specialized knowledge to achieve diagnostic accuracy like that of expert centers by simply measuring coronary plaque volume and analyzing histograms. This democratization of diagnostic capability could significantly impact clinical practice, particularly in areas with limited access to specialized cardiac imaging expertise. Furthermore, this approach lays the groundwork for potential future developments, such as automated systems that could alert referring physicians or radiologists to the presence of vulnerable plaques. Such advancements could enhance early detection and risk stratification in coronary artery disease, potentially improving patient outcomes through timely interventions.

Study Limitations

The limitations of this study should be acknowledged. First, this investigation was conducted as a small-scale, single-center study following a single protocol. The limited sample size prevented us from performing multiple comparisons across the various energy levels examined in this study. Our methodological approach relied on cross-validation techniques without establishing a separate test group, which underscores the need for larger-scale studies with dedicated test sets to fully assess the model’s generalizability. Second, although the study used IVUS to evaluate the results, there was a selection bias in the participating centers. Last, the small sample size precluded further analyses, such as subgroup comparisons based on patient data or multiple comparisons across different parameters. These limitations highlight the need for more comprehensive, multi-center studies with larger sample sizes to validate and extend our findings.

Conclusions

Spectral imaging with dual-energy CT can enhance the diagnostic performance of machine learning using CT histograms for coronary plaque characterization.

Acknowledgment

We thank the team at the Minamino Cardiovascular Hospital for their help with performing the study.

Sources of Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Disclosures

T.N. has received research support from Nemoto Kyorindo Co., Ltd. T.H. has received research support from Canon Medical Systems. The department of diagnostic imaging analysis, to which M.K. belongs, is an endowed chair supported by Philips Healthcare. The Nemoto Kyorindo Co., Ltd, Philips Healthcare and Canon Medical Systems had no control over the interpretation, writing, or publication of this work.

IRB Information

This retrospective study received institutional review board approval (MJ-023, ‘Coronary plaque assessment by Spectral CT’). Written informed consent was waived.

Data Availability

The datasets generated or analyzed during the study are available from the corresponding author on reasonable request.

References
 
© 2024, THE JAPANESE CIRCULATION SOCIETY

This article is licensed under a Creative Commons [Attribution-NonCommercial-NoDerivatives 4.0 International] license.
https://creativecommons.org/licenses/by-nc-nd/4.0/
feedback
Top