Radiomics analysis for predicting malignant cerebral edema in patients undergoing endovascular treatment for acute ischemic stroke
PDF
Cite
Share
Request
Neuroradiology - Original Article
P: 402-409
March 2023

Radiomics analysis for predicting malignant cerebral edema in patients undergoing endovascular treatment for acute ischemic stroke

Diagn Interv Radiol 2023;29(2):402-409
1. Department of Radiology, Center for Rehabilitation Medicine, Zhejiang Provincial People’s Hospital, Affiliated People’s Hospital, Hangzhou Medical College, Zhejiang, China
2. Department of Radiology, The First People’s Hospital of Daishan, Zhejiang, China
No information available.
No information available
Received Date: 10.01.2022
Accepted Date: 02.02.2023
Publish Date: 29.03.2023
PDF
Cite
Share
Request

ABSTRACT

PURPOSE

Radiomics analysis is a promising image analysis technique. This study aims to extract a radiomics signature from baseline computed tomography (CT) to predict malignant cerebral edema (MCE) in patients with acute anterior circulation infarction after endovascular treatment (EVT).

METHODS

In this retrospective study, 111 patients underwent EVT for acute ischemic stroke caused by middle cerebral artery (MCA) and/or internal carotid artery occlusion. The participants were randomly divided into two datasets: the training set (n = 77) and the test set (n = 34). The clinico-radiological profiles of all patients were collected, including cranial non-contrast-enhanced CT, CT angiography, and CT perfusion. The MCA territory on non-contrast-enhanced CT images was segmented, and the radiomics features associated with MCE were analyzed. The clinico-radiological parameters related to MCE were also identified. In addition, a routine visual radiological model based on radiological factors and a combined model comprising radiomics features and clinico-radiological factors were constructed to predict MCE.

RESULTS

The areas under the curve (AUCs) of the radiomics signature for predicting MCE were 0.870 (P < 0.001) and 0.837 (P = 0.002) in the training and test sets, respectively. The AUCs of the routine visual radiological model were 0.808 (P < 0.001) and 0.813 (P = 0.005) in the training and test sets, respectively. The AUCs of the model combining the radiomics signature and clinico-radiological factors were 0.924 (P < 0.001) and 0.879 (P = 0.001) in the training and test sets, respectively.

CONCLUSION

A CT image-based radiomics signature is a promising tool for predicting MCE in patients with acute anterior circulation infarction after EVT. For clinicians, it may assist in diagnostic decision-making.

Main points

• The development of malignant cerebral edema (MCE) may reduce the benefit–risk ratio of endovascular treatment (EVT).

• Computed tomography-based radiomics analysis may be used to predict MCE.

• A radiomics signature was built to screen patients at high risk of MCE after EVT.

The safety and efficacy of endovascular treatment (EVT) for anterior circulation large vessel occlusive stroke (LVOS) have been verified in randomized controlled trials. However, nearly 45% of the cases with LVOS cannot recover after EVT.1,2 Recent evidence has revealed that malignant cerebral edema (MCE) commonly occurs after EVT, and the development of MCE may reduce the benefit–risk ratio of EVT.3,4,5 MCE is a severe complication of stroke that occurs in approximately 10% of all patients with stroke. The prognosis for MCE is poor, with a mortality rate of 80% when treated conservatively.6 Early decompressive surgery effectively lowers mortality and improves the clinical outcomes of patients with MCE.7 The addition of decompressive craniectomy to the best medical therapy can reduce the mortality rate of MCE by 50%–75%.8 Hence, the early identification of patients at a high risk of MCE is critical for therapeutic decision-making. The imaging predictors for MCE that have been studied so far primarily involve the characterization of the infarction size, neurovascular condition, and brain perfusion of the patient.7

Diffusion-weighted imaging (DWI) is a valuable technique for identifying the ischemic core. Rapid sequence magnetic resonance imaging (MRI) using DWI and fluid-attenuated inversion recovery sequences has gradually emerged as a new screening modality in people with suspected acute ischemic stroke (AIS), especially where there might be delays in obtaining full sequence brain imaging. However, some limitations still exist, such as limited access at some emergency departments and the contraindications of MRI examination for some patients. Computed tomography (CT) is a commonly used imaging technique in clinical practice for diagnosing and treating AIS. However, some reported prognostic variables (i.e., MCE) are imprecise when using conventional CT methods (i.e., cranial CT scan).7 Therefore, new approaches that can promptly identify MCE are needed.

Conventional CT images contain information that cannot be evaluated visually but may be assessed using image analysis tools. Image textures correspond to the brightness value and the position of image pixels. Image texture analysis is a quantitative image analysis method that can improve the accuracy of the diagnostic or classification information extracted from different imaging methods, such as ultrasound, CT, and MRI.9 The correlation between the radiomics signature based on image textures and the progression of MCE after EVT is not fully known. This study hypothesized that the radiomics signature might be an early predictor of MCE in patients with acute anterior circulation infarction after EVT. Furthermore, this study evaluated the predictive potential of a single radiomics signature and a model comprising a radiomics signature as well as clinical factors and routine radiological parameters for MCE.

Methods

This retrospective study was approved by the Zhejiang Provincial People’s Hospital’s Ethics Committee (approval number: 2020QT038), and informed consent was waived due to its retrospective nature. The clinical data and radiological findings of all subjects were obtained from the medical records and the picture archiving and communication system of our hospital.

Study population

The clinical and radiological data of patients diagnosed with AIS caused by middle cerebral artery (MCA) and/or internal carotid artery (ICA) occlusion and who had undergone EVT were retrospectively reviewed from the institutional databases between January 2018 and August 2021. Finally, a total of 111 patients were included (Figure 1).

Figure 1

The inclusion criteria were as follows: (1) patients were diagnosed with AIS and underwent EVT; (2) baseline cranial non-contrast-enhanced CT (NCCT), CT angiography (CTA), CT perfusion (CTP) images, and the follow-up CT/MRI images in case of neurological deterioration or within three days after EVT were available; and (3) patients had a causative occlusion of the MCA and/or ICA. The exclusion criteria were as follows: (1) the presence of AIS in the areas supplied by the posterior circulation or bilateral acute ischemic lesions demonstrated by neuroimaging; (2) a modified Rankin Scale score of over two before admission; (3) the presence of symptomatic intracranial hemorrhage; and (4) difficulty in image interpretation due to incomplete images or artifacts.

Endpoint

The primary endpoint was the occurrence of MCE, defined as a midline shift of over 5 mm at the pineal gland or septum pellucidum on the follow-up imaging, with the need for decompressive hemicraniectomy.3,10,11,12

Imaging protocol

All participants underwent baseline CTP and NCCT scans and had follow-up CT or MRI. CTP and NCCT were performed using a 640-slice CT (Aquilion ONE TSX-301A, Toshiba).

The whole-brain NCCT was performed in one rotation (detector width: 16 cm). Moreover, CTP was acquired after the administration of a contrast agent (50 mL; Omnipaque, GE Healthcare) through intravenous (IV) injection at a rate of 6 mL/s, followed by 50 mL of saline (acquisition parameters: scanning width, 5 mm; scanning coverage, 240 mm; 128 mAs; and 120 kV). At 7 s after the injection of the contrast agent, a pulsed full rotation scan was performed with 18 time points acquired over 1 min and a total image acquisition time of 9.5 s. Multiphase CTAs, including venous peak phase, venous late phase, and arterial peak phase images, were generated from CTP images using MIStar (Apollo Medical Imaging Technology).

Clinical and routine radiological assessment

The baseline neurological status of patients at admission was evaluated using the National Institutes of Health Stroke Scale (NIHSS). The recanalization status after EVT was assessed using the modified Thrombolysis in Cerebral Infarction grading system.

Radiological images, including baseline NCCT, CTP, and CTA, as well as follow-up images, were reviewed and evaluated by two neuroradiologists (H and X, with 6 and 10 years of clinical experience, respectively). Any discordant interpretations between the two raters were resolved by consensus through discussion.

The Alberta Stroke Program Early CT Score (ASPECTS) was calculated by subtracting the score of the low-density regions from a total score of 10 on NCCT images.13

On CTA images, the thrombosis burden was quantified to assess the clot in the anterior circulation using the clot burden score (CBS), which starts at 10 and is deducted based on the loss of contrast opacification on CTA.14 The leptomeningeal collateral status was evaluated on a 0–3 scale: 0, absent collateral filling; 1, ≤50% collateral filling; 2, >50% but <100% collateral filling; and 3, ≥100% collateral filling when compared with the non-occluded side.15

On CTP images, the ischemic core was defined by a decrease in the relative cerebral blood flow to less than 30% of normal. The ischemic lesion was defined as brain tissues with a delay time (DT) of >3 s based on MIStar.

Image preprocessing and radiomics feature selection

Radiomics analysis was conducted based on baseline cranial NCCT (Figure 2).

Figure 2

Image data preprocessing, including gray-level discretization, intensity normalization, and image interpolation, was carried out using the AK software (version 3.0.0. R, Artificial Intelligence Kit, GE Healthcare).

The region of interest (ROI) of the whole MCA territory was manually segmented on a 3D volume of interest on NCCT images using the ITK-SNAP software (www.itksnap.org). The intra- and inter-rater agreement on the ROI segmentation was evaluated using intraclass correlation coefficients. Rater H segmented the ROIs twice within a 2 month interval. Intra-rater ICC was determined by comparing the two measurements of rater H. Inter-rater ICC was calculated by comparing the ROIs in 30 randomly selected patients (measured by raters X and H).

The extracted features included haralick, histogram, formfactor, gray level co-occurrence matrix, run-length matrix, and gray level size-zone matrix. For each patient, 396 features were extracted, which were then standardized to remove the unit limits of the data. Dimension reduction was carried out using the correlation test, Mann–Whitney U test, and analysis of variance. The features were further selected using the least absolute shrinkage and selection operator (LASSO).

Development of the radiomics signature

Multivariate logistic regression analysis was used to construct the radiomics signature, which was then used to predict MCE based on selected features using the LASSO. Subsequently, the radiomics score (rad-score) of each patient was calculated. The formula was derived from the training set data and was then used to determine the rad-score of the patient in the test set. The predictive efficiency of the radiomics signature was evaluated using the receiver operating characteristic (ROC) curves in both the training and test sets.

Development of the prediction model for malignant cerebral edema

All parameters (i.e., sex, age, history of hypertension, atrial fibrillation, hyperlipidemia, diabetes mellitus, alcohol abuse, smoking, baseline NIHSS, hyperdense vessel sign of the MCA, ASPECTS, collateral score, occlusion site, CBS, ischemic lesion volume, ischemic core volume, IV-type plasminogen activator (tPA), time from onset to groin puncture, reperfusion grade, and the radiomics signature) were screened to identify variables associated with MCE. Furthermore, multivariate regression analysis was performed to construct a prediction model for MCE. The model calibration was assessed using the Hosmer–Lemeshow test. The ROC curve was used to evaluate the predictive ability of the model.

Statistical analysis

The MedCalc v15.2.2, SPSS v21.0, and Microsoft R Open v3.3.1 software were used for statistical analyses. Descriptive statistics were presented as median (first quartile and third quartile) for non-normally distributed variables, and the mean ± standard deviation was used for the normally distributed variables. Categorical variables were reported as frequencies (percentages). The LASSO analysis was conducted based on the minimum criterion by 10 fold cross-validation. The Student’s t-test, Mann–Whitney U test, Pearson’s chi-square test, and Fisher’s exact test were used to identify MCE-associated variables. The Youden index in the ROC analysis was used to determine the best cut-off point of CTP parameters and NCCT ASPECTS for identifying patients with MCE. Multivariate regression analysis was performed to establish the radiomics signature, the routine visual radiological model, and the combined model for predicting MCE. The ROC curves were compared using the method developed by DeLong et al.16 The significance level was established as α = 0.05.

Results

The intra- and inter-rater reliability

The intra-rater agreement on the ROI segmentation between the two measurements from the same rater was 0.896 (P < 0.001). The inter-rater agreement between the two raters was 0.842 (P < 0.001).

Baseline characteristics

A total of 111 patients were randomly split into two datasets at a ratio of 7:317 as follows: the training set (n = 77) and the test set (n = 34). There were no significant differences in the clinical or routine radiological characteristics between the two sets (Table 1). The rad-score, ASPECTS, ischemic core volume, DT >3 s lesion volume, and IV-tPA were significantly different between the MCE (+) and MCE (-) groups in the training set (Table 2).

Table 1
Table 2

The radiomics signature

After identifying 187 features using the Mann–Whitney U test and analysis of variance, 22 features were retained after the Spearman correlation analysis. Moreover, the LASSO analysis was performed, and 15 features remained (Supplementary Table 1). The rad-score was obtained using the LASSO model, with a linear combination of the 15 features (Supplementary Materials). The values of the 15 features of each subject were placed into the calculation formula, and the rad-score was obtained to indicate the predictive efficiency for MCE.

Table 1

The sensitivity, specificity, and AUC of the radiomics signature for predicting MCE in the training set were 86.36%, 76.36%, and 0.870 [95% confidence interval (CI): 0.773–0.967], respectively (Figure 3). The radiomics signature showed satisfying goodness-of-fit (P = 0.361). The sensitivity, specificity, and AUC in the test set were 90.00%, 83.33%, and 0.837 (95% CI: 0.655-1), respectively (Figure 3).

Figure 3

The routine visual radiological model

The routine radiological parameters for the prediction of MCE were as follows: 1) ASPECTS with an AUC of 0.800 (95% CI, 0.683–0.917; P < 0.001); 2) ischemic core volume with an AUC of 0.793 (95% CI, 0.675–0.910; P < 0.001); and 3) DT >3 s lesion volume with an AUC of 0.725 (95% CI, 0.580–0.870; P = 0.002). The best cut-off values for the ASPECTS, ischemic core volume, and DT >3 s lesion volume were 7, 37, and 235, respectively (Table 3).

Table 3

The sensitivity, specificity, and AUC of the routine visual radiological model based on the above radiological parameters were 77.27%, 78.18%, and 0.808 (95% CI: 0.687–0.929), respectively, in the training set (Figure 3). The routine visual radiological model showed satisfying goodness-of-fit (P = 0.200). The sensitivity, specificity, and AUC in the test set were 70.00%, 87.50%, and 0.813 (95% CI: 0.649–0.976), respectively (Figure 3).

Figure 3

The combined model

The combined model for predicting MCE included the rad-score, ASPECTS, ischemic core volume, DT >3 s lesion volume, and IV-tPA.

In the training set, the sensitivity, specificity, and AUC of the combined model were 77.27%, 96.36%, and 0.924 (95% CI: 0.850–0.998), respectively (Figure 3). The combined model showed satisfying goodness-of-fit (P = 0.662). In the test set, the sensitivity, specificity, and AUC were 80.00%, 91.67%, and 0.879 (95% CI: 0.712–1), respectively (Figure 3).

Figure 3

The participants were divided into high- and low-risk groups according to the optimal diagnostic cut-off value of the model (0.390). The number of patients with MCE (+) was significantly different between the high- and low-risk groups in the training set (P < 0.001). Additionally, the results showed that the patients categorized into the high-risk group in the test set were all patients with MCE (+) (Figure 4).

Figure 4

The AUCs of the radiomics signature and the combined model were significantly different (P = 0.049). The AUCs of the combined model and the routine visual radiological model were also significantly different (P = 0.020). However, there was no significant difference in the AUC between the radiomics signature and the routine visual radiological model (P = 0.354) (Table 4).

Table 4

Discussion

MCE is a life-threatening condition with high mortality rates. It is usually secondary to the acute occlusion of the proximal MCA trunk or the distal ICA.18 Pathophysiologically, early damage to cerebrovascular autoregulation in the tissues surrounding the infarction, increased vascular permeability, and the loss of the integrity of the endothelial basal layer altogether play important roles in the progression of MCE.19,20 The early detection of MCE and the use of decompressive surgery can improve the clinical outcomes of patients with MCE.21,22 However, the current methods for detecting MCE are mainly reactive: waiting for signs of clinical deterioration and mass effect to appear before surgical intervention. Therefore, the identification of early predictors for MCE is of high importance.

Texture feature analysis is widely used in radiomics, which allows the quantification of internal heterogeneity that may not be perceived by the naked eye, thus promoting the identification and classification of different tissues.23 Until recently, few studies have performed texture analysis of cerebral infarction. CT radiomics based on texture features showed good performance in identifying hyperacute or AIS lesions, assessing the extent of ischemic lesions, and determining the time from symptoms onset in basal ganglia infarction.24,25,26 Kuang et al.27 developed an automated method to compute the ASPECTS on NCCT images from patients with AIS. Additionally, texture analysis based on MRI effectively identified the presence of ischemic stroke lesions and detected hemorrhagic transformation in patients with AIS.28,29 Furthermore, in patients with stroke, the texture features of brain MRI may provide early markers for poststroke cognitive impairment.30 The potential of texture features for predicting MCE in patients with AIS caused by MCA and/or ICA occlusion after EVT was evaluated. By constructing a radiomics model, it was found that the radiomics signature based on texture features was a valuable tool for identifying patients at a high risk of MCE after EVT, with AUCs of 0.870 and 0.837 in the training and test sets, respectively.

By analyzing the correlation between clinico-radiological factors and MCE, it was found that the ASPECTS on NCCT, ischemic core volume, DT >3 s lesion volume, and IV-tPA were correlated with MCE. The ASPECTS, ischemic core volume, and DT >3 s lesion volume can indicate the extent of ischemia or infarct lesion. Compared with patients with MCE (-), patients with MCE (+) showed lower ASPECTS, larger ischemic core volume, and larger DT >3 s lesion volume. This finding suggested that ischemia or infarction core during the early stage was more severe in patients with MCE (+) than in patients with MCE (-), which was consistent with the findings of Tracol et al.31 A midline shift of over 3.9 mm and an infarct volume of more than 220 mL are predictors of severe brain edema and herniation.32 This study’s results further showed that the proportion of the use of IV-tPA in patients with MCE (+) was lower than that in patients with MCE (-). It suggested that IV-tPA could reduce the incidence of MCE, which was in line with the findings of Fuhrer et al.33 Treatment with IV-tPA lowered the odds of experiencing MCE, with a reduced odds ratio of 0.88 (95% CI: 0.83–0.94).

With the development of analytical approaches and high-throughput technology, the multiparameter approach has recently become a useful tool for improving diagnostic performance.34 Therefore, this study built a model combining the radiomics signature, clinical factors, and routine visual radiological factors to predict MCE. Parameters related to MCE, such as the radiomics signature, ASPECTS, ischemic core volume, DT >3 s lesion volume, and use of IV-tPA, were included in this model. The combined model exhibited better differentiation capabilities than the radiomics signature or the routine visual radiological model, with AUCs of 0.924 and 0.879 in the training and test sets, respectively.

The limitations of this study should be noted. First, information and selection biases may exist due to the retrospective nature of this study. Therefore, the validation of the results in a prospective study is required. Second, the sample size was considerably small, and the proposed models were not externally validated. Third, the difference between the combined model and radiomics signature or the routine visual radiological model was small, with only borderline statistical significance between the combined model and radiomics signature. Therefore, a large-scale, multicenter study is required to validate the proposed multimodal diagnostic approach. Despite the abovementioned limitations, a prediction model for MCE after EVT was established in this study by combining the radiomics signature, clinical factors, and routine radiological factors, which may facilitate the accurate early prediction of MCE in patients with acute anterior circulation infarction after EVT.

In conclusion, this study constructed a model for predicting MCE in patients with acute anterior circulation infarction after EVT. These findings may contribute to the accurate and early prediction of MCE and assist clinical decision-making.