You're viewing our updated article page. If you need more time to adjust, you can return to the old layout.

ORIGINAL RESEARCH article

Front. Neurol., 27 February 2025

Sec. Applied Neuroimaging

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1471274

Machine learning, clinical-radiomics approach with HIM for hemorrhagic transformation prediction after thrombectomy and treatment

  • 1. Department of Radiology, The Fourth Affiliated Hospital of School of Medicine, and International School of Medicine, International Institutes of Medicine, Zhejiang University, Yiwu, Zhejiang, China

  • 2. Department of Neurology, The Fourth Affiliated Hospital of School of Medicine, and International School of Medicine, International Institutes of Medicine, Zhejiang University, Yiwu, Zhejiang, China

  • 3. Department of Neurosurgery, The Fourth Affiliated Hospital of School of Medicine, and International School of Medicine, International Institutes of Medicine, Zhejiang University, Yiwu, Zhejiang, China

Article metrics

View details

3

Citations

1,7k

Views

579

Downloads

Abstract

Background:

This study aimed to develop a clinical-radiomics model using hyperattenuated imaging markers (HIM), characterized by hyperattenuation on head non-contrast computed tomography immediately after thrombectomy, to predict the risk of hemorrhagic transformation (HT) in patients undergoing endovascular mechanical thrombectomy (MT).

Methods:

A total of 159 consecutive patients with HIM were screened immediately after MT for inclusion. The datasets were randomly divided into training and test cohorts at a ratio of 8:2. An optimal machine learning (ML) algorithm was used for model development. Subsequently, models for clinical, radiomics, and clinical-radiomics were developed. The performance of the models was measured using receiver operating characteristic (ROC) and decision curve analyses (DCA). The interpretability and predictor importance of the model were analyzed using Shapley additive explanations.

Results:

Of the 159 patients, 100 (62.9%) exhibited HT. The support vector machine (SVM) was the optimal ML algorithm for constructing the models. In predicting HT, the areas under the curve (AUCs) of the clinical model were 0.918 (95% confidence interval [CI] = 0.869–0.966) in the training cohort and 0.854 (95% CI = 0.724–0.984) in the test cohort. The AUCs of the radiomics model were 0.869 (95% CI = 0.802–0.936) and 0.829 (95% CI = 0.668–0.990), while those of the clinical-radiomics model were 0.944 (95% CI = 0.905–0.984) and 0.925 (95% CI = 0.832–1.000).

Conclusion:

The suggested clinical-radiomics model based on HIM is a reliable method that can provide a risk evaluation of HT in individuals undergoing MT.

1 Introduction

Endovascular mechanical thrombectomy (MT) is the recommended treatment for acute ischemic stroke (AIS) because it improves reperfusion rates and clinical outcomes (1). Despite advancements in treatment strategies, hemorrhagic transformation (HT) remains a significant determinant of patient prognosis, with over 50% of patients experiencing unfavorable outcomes (2). Symptomatic intracerebral hemorrhage can lead to death in some patients (3). Furthermore, asymptomatic intracerebral hemorrhage and subarachnoid hemorrhage (SAH) contribute to adverse outcomes (4, 5). Appropriate postoperative treatment can improve the prognosis. Therefore, early prediction of HT is crucial for adjusting treatment after intervention. It is particularly important to determine whether blood pressure control and drug therapy should be initiated to prevent early re-occlusion of blood vessels caused by endothelial injury.

The hyperattenuated imaging marker (HIM) detected immediately on non-contrast computed tomography (NCCT) after MT is the most accessible and earliest imaging marker for predicting postoperative complications (6, 7). HIM represents areas of hyperattenuation on post-thrombectomy cranial NCCT scans, often including iodine contrast extravasation (ICE) and occasionally containing hemorrhage. Postoperatively, these extravasated contrast agents are gradually absorbed, while hemorrhage may either start or increase. Numerous studies have reported the postoperative presence of HIM in patients, ranging from 32.9 to 87.5% (8, 9). The presence of HIM indicates disruption of the blood–brain barrier and blood-cerebrospinal fluid barrier, which are closely related to bleeding, making it the most commonly used imaging marker for hemorrhage prediction after surgery (4, 7, 10, 11). Several studies have demonstrated that the presence of metallic hyperdensity signs and cortical involvement in HIM are more sensitive and specific for predicting HT (12, 13). However, there is currently no standardized predictive method with acceptable sensitivity and specificity that can accurately predict HT in all patients with HIM.

Radiomics, a technique that involves the extraction of numerous characteristics from medical images, is becoming increasingly popular for improving the diagnosis and management of ischemic stroke. However, previous studies using radiomic features based on HIM have primarily focused on detecting hemorrhage in the brain parenchyma, overlooking the inclusion of SAH (14, 15). In reality, both intracerebral hemorrhage and SAH can affect the subsequent use of anticoagulants. Accordingly, it is important to consider both types of hemorrhage when using radiomics for the prediction and treatment planning of patients with ischemic stroke.

Machine learning (ML) can be used to integrate radiomic features with high-risk clinical factors to develop a hemorrhage prediction model for stroke (16, 17). Consequently, this study aimed to use radiomic attributes derived from HIM in NCCT and clinical characteristics to develop a predictive model for both intracerebral hemorrhage and SAH in patients undergoing MT using ML algorithms.

2 Materials and methods

2.1 Ethical approval of the study protocol

Approval Number KY2023112 was granted by the Ethics Committee of the Fourth Affiliated Hospital, Zhejiang University College of Medicine. The committee exempted informed consent. All clinical studies followed the principles of the Declaration of Helsinki.

2.2 Patients and study design

A retrospective review was conducted on patients who experienced AIS due to blockage in the large blood vessels inside the skull and received endovascular MT between June 2016 and November 2021 at the Fourth Affiliated Hospital, Zhejiang University School of Medicine. The latest guidelines were followed at the time, and they included the indications and contraindications of MT and thrombolysis. The general clinical features, laboratory examinations, clinical presentations, and imaging data of the patients were collected within the past 90 days.

The inclusion criteria were as follows: (1) Patients scheduled for MT. (2) Patients who underwent head NCCT or magnetic resonance imaging (MRI) after MT; (3) initial postoperative head NCCT was conducted within 30 min after MT; (4) HIM, which is defined as an area of hyperattenuation visible in the brain parenchyma and subarachnoid space, was detected on the initial head NCCT following MT. The exclusion criteria included the following: (1) Patients for whom MT was aborted due to unfavorable vascular anatomy or who were switched to medical therapy; (2) NCCT follow-up time after MT of less than 19 h (18); (3) the judgment of HIM or hemorrhage was influenced by artifacts like metal or motion artifacts; (4) the use of iodinated contrast before preoperative CT affected the determination of HT. Figure 1 displays a flowchart illustrating the study participants.

Figure 1

Figure 1

The flowchart of the patient recruitment pathway. MT, mechanical thrombectomy; NCCT, non-contrast CT; HIM, Hyperattenuated Imaging Marker.

2.3 Imaging

Within 30 min after thrombectomy, initial postoperative head NCCT was performed using a 64-row spiral CT scanner (Somatom Definition AS, Siemens Healthineers, Forchheim, Germany) and a 62-row spiral scanner (Optima CT620, GE Medical Systems, Milwaukee, WI, US). The scanning parameters included an axial mode with a tube voltage of 120 kVp and a tube current of 250–300 mAs. The scanning range was extended from the skull base to the cranial roof with a section thickness of 5 mm. Reconstruction was performed using the standard algorithm. Two neuroradiologists with over a decade of professional experience evaluated the imaging data randomly and were unaware of the clinical circumstances. Two evaluators discussed conflicting images to reach a consensus.

2.4 Reference standard

The HT was classified according to the Heidelberg Classification, adhering to the following criteria: 1. Hyperattenuation remained evident on the 24-h CT scan, and 2. Follow-up CT or MRI conducted within 90 days after surgery revealed hemorrhage within the range of infarction. The maximum Hounsfield unit (HUmax) was defined as the highest HU value within a specified Region of Interest (ROI), measuring 3 × 3 pixels. Two neuroradiologists (S.H. and Z.W.) with >5 years of experience independently evaluated the training and test groups. They were blinded to patient outcomes and assessed the groups separately. Any discrepancy was resolved by consensus.

2.5 Data preprocessing

This study included 159 patients. Patients were randomly divided into a training cohort (127 patients) and a test cohort (32 patients). Data on several baseline and preprocedural factors were collected, including age, gender, diabetes mellitus, coronary artery disease, atrial fibrillation, hypertension, smoking and drinking, thrombolysis, and baseline National Institute of Health Stroke Scale (NIHSS) score. These data were obtained by reviewing the patient’s medical records, procedure notes, and progress and follow-up notes. Two experienced neuroradiologists, who were unaware of the patient’s clinical information, determined the baseline Alberta Stroke Program Early CT Score (ASPECTS) in the anterior and posterior circulation, subarachnoid HIM (sHIM), and HUmax through discussion. Furthermore, data regarding the surgical procedure, including the type of stent, number of stent passes, stent placement, and modified Thrombolysis in Cerebral Infarction scores, were collected.

2.6 Radiomics features extraction

On each initial postoperative CT image, the ROI of the HIM was manually segmented along the HIM contour, incorporating both parenchyma and sulcus. This process was conducted using the ITK-SNAP software (version 3.8.01). Both radiologists were blinded to the clinical information and final outcomes.

To eliminate any potential variations in CT images obtained using different CT scanners, NCCT images were reconstructed using a voxel size of 1 × 1 × 1 mm3 and gray-scale discretization. Two radiologists independently segmented the images of every HIM and measured them using a method that ensured double anonymization. An intraclass correlation coefficient (ICC) ≥ 0.75 was considered robust.

A Pyradiomics in-house feature analysis program2 was used to extract radiomic features from HIM in NCCT. From the NCCT images, 1834 radiomics features were extracted. Various techniques were used to extract texture features, including first-order statistics, shape-based analysis, gray-level co-occurrence matrix (GLCM), gray-level run length matrix (GLRLM), neighboring gray-tone difference matrix (NGTDM), gray-level size zone matrix (GLSZM), and gray-level dependence matrix (GLDM).

2.7 Finding the best ML algorithm to build models

All clinical data and robust radiomics features were combined to develop the combined models. A Mann–Whitney U test and feature screening were performed for all features. The features with a p-value <0.05 were included. Spearman’s rank correlation coefficient was calculated for the features exhibiting high repeatability to evaluate the association between features. If the correlation coefficient between any two characteristics >0.9, one of the characteristics would be preserved. To ensure a comprehensive depiction of the features, a greedy recursive deletion strategy was used for feature filtering. This process includes removing the characteristics with the highest duplication in the existing collection during each cycle. For signature construction, the discovery dataset was subjected to the least absolute shrinkage and selection operator (LASSO) regression model. By applying a regularization weight λ, LASSO reduces the regression coefficients to zero and assigns a value of zero to the coefficients of the irrelevant features. To determine the best λ, a 10-fold cross-validation was applied, and the λ value that resulted in the lowest cross-validation error was selected using the minimum criteria. The included features with coefficients that were not zero were used to fit the regression model and were incorporated into the model. LASSO regression modeling was performed using the Python scikit-learn library.

Following the LASSO feature screening, the selected features were incorporated into the ML models. To develop the HT prediction model, the effects of seven ML algorithms [Support Vector Machine (SVM), Light Gradient Boosting Machine, k-nearest neighbor, random forest, AdaBoost, eXtreme Gradient Boosting, and Gradient Boosting] were compared before modeling. The optimal algorithm was identified through 5-fold cross-validation in the training cohort. The SVM model for ML was identified as having the greatest mean area under the curve (AUC) of the receiver operating characteristic (ROC) curve in the test group (Table 1). Hence, the chosen characteristics were fed into the SVM algorithms to develop the risk model. To obtain the optimized subset of features, a cross-verification was conducted using a 5-fold method.

Table 1

Model Cohort Accuracy AUC 95% CI Sensitivity Specificity PPV NPV Precision Recall F1
SVM Train 0.89 0.944 0.9046–0.9837 0.938 0.809 0.893 0.884 0.893 0.938 0.915
SVM Test 0.875 0.925 0.8325–1.0000 0.9 0.833 0.9 0.833 0.9 0.9 0.9
KNN Train 0.858 0.912 0.8626–0.9619 0.863 0.851 0.908 0.784 0.908 0.863 0.885
KNN Test 0.844 0.842 0.6900–0.9934 0.95 0.727 0.826 0.889 0.826 0.95 0.884
Random-Forest Train 1 1 1.0000–1.0000 1 1 1 1 1 1 1
Random-Forest Test 0.719 0.815 0.6669–0.9623 0.55 1 1 0.571 1 0.55 0.71
XGBoost Train 1 1 1.0000–1.0000 1 1 1 1 1 1 1
XGBoost Test 0.75 0.717 0.5157–0.9176 0.8 0.667 0.8 0.667 0.8 0.8 0.8
LightGBM Train 0.843 0.921 0.8765–0.9647 0.875 0.787 0.875 0.787 0.875 0.875 0.875
LightGBM Test 0.844 0.846 0.7015–0.9902 0.9 0.75 0.857 0.818 0.857 0.9 0.878
Gradient-Boosting Train 0.953 0.994 0.9873–1.0000 0.925 1 1 0.887 1 0.925 0.961
Gradient-Boosting Test 0.75 0.708 0.5033–0.9133 0.85 0.778 0.773 0.7 0.773 0.85 0.81
AdaBoost Train 0.85 0.947 0.9132–0.9799 0.775 0.979 0.984 0.719 0.984 0.775 0.867
AdaBoost Test 0.75 0.815 0.6683–0.9609 0.75 0.75 0.833 0.643 0.833 0.75 0.789

Performance of ML combined models.

2.8 Building models

A radiomics model was developed for the training group using the most important radiomics characteristics that were ultimately chosen following the same procedure as that used for the combined model. Additionally, the radiomics model was validated in the test cohort. In the training cohort, the collected clinical characteristics were integrated into the LASSO regression model to identify the most significant characteristics, retaining only those with coefficients that were not zero. The SVM algorithm was then used to develop a clinical model by employing these characteristics. The efficiency of the model was evaluated using the test cohort. A flowchart of the radiomics model construction is shown in Figure 2.

Figure 2

Figure 2

The process of building the radiomics model. NCCT, non-contrast CT; ROI, region of interest; ICC, intraclass correlation coefficient; LASSO, least absolute shrinkage and selection operator; MSE, mean square error; ROC, receiver operating characteristic; DCA, decision curve analysis; SHAP, Shapley Additive explanations.

The combined model was developed in the training cohort and was subsequently verified in the test cohort. The diagnostic efficacy of the combined model was evaluated using ROC curves. To assess the clinical applicability of the predictive models, a decision curve analysis (DCA) was conducted. Furthermore, the performance of these models was evaluated using Shapley Additive explanations (SHAP) analysis to determine the significance of each feature in predicting HT.

2.9 Statistics

The statistical package for the social sciences software (version 26.0, IBM) was used to perform statistical analyses. For continuous variables, the disparity was determined using either the Student’s t-test or the Mann–Whitney U test, depending on suitability. Continuous data are presented as mean ± standard deviation or median with interquartile range, depending on the distribution of the data. Differences in categorical variables were determined using either the Chi-square or Fisher’s exact test, depending on the circumstances. The significance threshold for all statistical tests was set at p < 0.05. Python software (version 3.11.13) was used to perform feature extraction and screening and to construct the models. The performance of the three models was evaluated using the AUC, 95% CI, accuracy, sensitivity, and specificity.

3 Results

3.1 Comparison of patient clinical characteristics

This study included 159 patients diagnosed with HIM, with 127 and 32 patients assigned to training and test cohorts, respectively. The initial attributes of patients in both the training and test groups are listed in Table 2. HT occurred in 100 of the 159 patients. The training and test cohorts differed significantly in thrombolysis (p = 0.027).

Table 2

Feature name Train Test P-value
Age, mean ± SD 67.22 ± 13.84 66.44 ± 15.59 0.781
Men, n(%) 81(63.78) 20(62.50) 0.894
Atrial fibrillation, n(%) 48(37.80) 9(28.13) 0.311
Hypertension, n(%) 72(56.69) 20(62.50) 0.555
Diabetes mellitus, n(%) 17(13.39) 4(12.50) 0.896
Coronary artery disease, n(%) 20(15.75) 8(25.00) 0.060
Drinking, n(%) 27(21.26) 8(25.00) 0.651
Smoking, n(%) 29(22.83) 9(28.13) 0.534
Thrombolysis, n(%) 48(37.80) 19(59.38) 0.027*
Baseline NIHSS, median (Q1, Q3) 16(12–20) 15(13–18) 0.893
Anterior circulation, n(%) 119(93.70) 30(93.75) 0.992
ASPECTS, median (Q1, Q3) 9(8–10) 9(8–10) 0.504
Humax ≥ 90, n(%) 35(25.55) 14(43.75) 0.077
sHIM, n(%) 46(36.22) 13(40.62) 0.647
mTICI 2b or 3, n(%) 113(88.98) 26(81.25) 0.242
Stent type, n(%) 0.486
SOLITAR 71(55.91) 19(59.38)
TREVO 25(19.69) 5(15,63)
SOLITAR + TREVO 16(12.60) 4(12.50)
others 15(11.81) 4(12.50)
Pass number, median (Q1, Q3) 2(1–3) 2(1–3) 0.958
Stent implantation, n(%) 25(19.69) 8(21.88) 0.784

Characteristics of the training and test sets.

*Represents P < 0.05; SD, standard deviation; HU, Hounsfield; NIHSS, National Institutes of Health Stroke Scale; Q1, first quartile; Q3, third quartile; ASPECTS, Alberta Stroke Program Early CT Score; mTICI, modified thrombolysis in cerebral infarction.

3.2 Establishment and performance of the radiomics model

A total of 1834 radiomic characteristics were obtained from the NCCT images after MT, with 1,681 characteristics (ICC ≥ 0.75) demonstrating acceptable consistency between observers. To develop the model, five features per patient were selected by excluding feature pairs with strong correlations.

In the training group, the radiomics model exhibited an AUC of 0.869 (95% confidence interval [CI] = 0.802–0.936) with sensitivity and specificity rates of 0.925 and 0.681, respectively. The test group exhibited an AUC of 0.829 (95% CI = 0.668–0.990) and well-balanced sensitivity (0.800) and specificity (0.833) (Table 3 and Figure 3).

Table 3

Model Cohort Accuracy AUC 95% CI Sensitivity Specificity PPV NPV Precision Recall F1
Combined Train 0.890 0.944 0.9046–0.9837 0.938 0.809 0.893 0.884 0.893 0.938 0.915
Combined Test 0.875 0.925 0.8325–1.0000 0.9 0.833 0.9 0.833 0.9 0.9 0.9
Radiomics Train 0.835 0.869 0.8018–0.9365 0.925 0.681 0.831 0.842 0.831 0.925 0.876
Radiomics Test 0.812 0.829 0.6679–0.9905 0.8 0.833 0.889 0.714 0.889 0.8 0.842
Clinical Train 0.819 0.918 0.8686–0.9665 0.725 0.979 0.983 0.676 0.983 0.725 0.835
Clinical Test 0.781 0.854 0.7242–0.9841 0.65 1 1 0.632 1 0.65 0.788

Performance of the three models.

Figure 3

Figure 3

Receiver operating characteristic (ROC) and DCA curves of the three models. (A) ROC curves for predicting HT in the training cohort. (B) The corresponding ROC curves for the test cohort. (C) DCA curves for the three models. The y-axis represents the net benefit, and the x-axis represents the threshold probability. The blue line, green line, and orange line represent the expected net benefit of predicting HT using the combined model, clinical model, and radiomics model, respectively. ROC, receiver operating characteristic; DCA, decision curve analysis; HT, hemorrhagic transformation.

3.3 Performances of the clinical models

Nine clinical factors, including HUmax ≥90, anterior circulation, sHIM, ASPECTS, male sex, drinking, coronary artery disease, baseline NIHSS score, and pass number, were selected after implementing LASSO feature screening. These factors were used to develop the clinical model.

In the training cohort, the clinical model demonstrated an AUC of 0.918 (95% CI = 0.869–0.966), a sensitivity of 0.725, and a specificity of 0.979. Similarly, the clinical model exhibited an AUC of 0.854 (95% CI = 0.724–0.984) in the test cohort with a sensitivity of 0.650 and specificity of 1.0. These results are presented in Table 3 and Figure 3.

3.4 Performances of the combined model

The models were developed by combining all the clinical data and robust radiomics features. After excluding the feature pairs with strong correlations, 17 features per patient were selected for further analysis. These results are presented in Table 3 and Figure 3.

In the training cohort, the combined model exhibited an AUC of 0.944 (95% CI = 0.905–0.984), with sensitivity and specificity rates of 0.938 and 0.809, respectively; whereas the AUC was 0.925 (95% CI = 0.832–1.000) in the test cohort, and the sensitivity and specificity were 0.900, and 0.833, respectively (Table 3 and Figure 3). The Delong test showed no significant differences between models: combined vs. radiomics (p = 0.172), combined vs. clinical (p = 0.316), and radiomics vs. clinical (p = 0.800). The DCA values of the three models for both the training and test cohorts are illustrated in Figure 3C. The DCA demonstrated that the combined model provided a significant clinical prediction advantage for the majority of threshold probabilities compared to the clinical and radiomics models.

3.5 Explanation and visualization of the combined model

After performing the prediction modeling of HT, the feature importance matrix plot was used to rank the most important variables. The plot revealed the influence of each factor on HT prediction. To demonstrate the influence of each characteristic on the model’s prediction, a summary chart of the SHAP values is depicted in Figure 4, illustrating the correlation between high or low SHAP values and the prediction model. The blue dots, which indicate lower values of logarithm_glrlm_ShortRunLowGrayLevelEmphasis, tend to be more concentrated in areas with a higher likelihood of HT. These results indicate that the SHAP values of this measurement exhibited an inverse relationship with the risk of HT occurrence. Contrarily, the risk of HT was high on the side, with a higher frequency of red dots indicating high Anterior_Circulation, Ibp_3D_m2_glszm_SmallAreaEmphasis, sHIM, and men, implying a positive correlation between the SHAP value of this index and the risk of HT.

Figure 4

Figure 4

The rank of importance for each feature in the combined model for HT prediction. The SHAP summary plot illustrates the distribution of patient features, arranged by importance, with the most critical feature at the top and the least important at the bottom. The color legend on the right side represents the influence of feature values on the outcome, with red indicating high values and blue indicating low values. Positive SHAp values (right on the x-axis) suggest a higher probability of HT, while negative values (left on the x-axis) indicate the opposite. SHAP, Shapley Additive exPlanations; sHIM, subarachnoid hyperattenuated imaging markers; ASPECTS, Alberta Stroke Program Early CT Score.

4 Discussion

Hemorrhagic transformation is a significant complication of MT in patients with AIS. Identifying patients at risk of developing HT is crucial for optimizing the management of stroke patients. We aimed to address this requirement by developing and evaluating interpretable ML models. The occurrence of HT after MT was predicted using baseline clinical information and imaging features based on HIM in these models. Our combined model, incorporating both clinical factors and radiomics features, demonstrated superior predictive performance. This model yielded an AUC of 0.944 in the training group and 0.925 in the test group. These results highlight the potential of our model to accurately predict the risk of post-MT HT in patients with AIS. The predictive capacity of the model was improved by combining clinical and radiomics characteristics, providing valuable information to optimize patient care and minimize the incidence of this severe complication.

The occurrence of HT is common after MT, with the most common type being intracerebral hemorrhage. Recent studies have demonstrated an increase in the number of reports of SAH as the indications for MT have expanded. In a study of anterior circulation MT, the incidence of SAH was 13.2%, and the incidence of intracerebral hemorrhage was 39.6%. SAH accounted for 1/4 of all hemorrhages (19). Intracerebral hemorrhage and SAH after thrombectomy mutually influence each other. Disruption of the blood–brain barrier or injury to the vascular wall in the brain parenchyma leads to SAH (4, 20), and vascular wall injury can also cause parenchymal damage (21, 22). Patients who develop SAH after thrombectomy have a poor prognosis, which is further exacerbated by the presence of intracerebral hemorrhage (23–25). SAH, like intracerebral hemorrhage, also affects timely postoperative treatment. However, multiple studies have separated them for analysis and have not included SAH as an outcome measure in predictive models (12, 26, 27). Optimal management and treatment strategies for patients with AIS can be achieved through the precise and prompt prediction of hemorrhage after thrombectomy. Patients without hemorrhage should receive appropriate perfusion maintenance and prevention of vessel re-occlusion to salvage more brain tissue. Patients with hemorrhage should have their blood pressure lowered and antiplatelet aggregation drugs discontinued to prevent further bleeding. This study included both intracerebral hemorrhage and SAH as outcome measures and achieved satisfactory predictive results.

HIM-based radiomics has demonstrated the HT prediction potential. Chen et al. (14) conducted a study using radiomic features of HIM within 1 h after MT to determine the presence of hemorrhage, achieving an AUC of 0.826. Ma et al. (15) developed a model that combined radiomics features of intracerebral HIM within 24 h post-surgery with clinical characteristics to predict HT. Their model demonstrated an AUC of 0.926 for intracerebral hemorrhage occurring 48 h after surgery. However, they used a broad definition of HIM and hemorrhage timing, which may have inadvertently excluded certain instances of HIM and hemorrhage that have been resolved.

Contrarily, our study followed stricter inclusion criteria that restricted the examination time of HIM to a maximum of 30 min. But the strict 30-min post-MT imaging window may limit generalizability. Future studies should broaden this timeframe to minimize potential bias. Moreover, unlike their study, we included SAH as an outcome measure, which is more applicable in clinical practice. Our clinical radiomics model demonstrated the best predictive performance, with AUCs of 0.944 and 0.925 in the training and test cohorts, respectively, confirming its superior predictive capability further, even though the Delong test did not reveal statistically significant differences between models.

Support vector machine, a widely utilized ML algorithm, is frequently used for tasks involving recognition and classification. The SVM has become a prominent algorithm for supervised learning and is widely used in the fields of pattern recognition and ML. The algorithm uses a linear decision boundary known as a hyperplane to identify the optimal hyperplane that maximizes the gap distance between two categories (28). Furthermore, SVM can address non-linear classification problems by transforming data into a space with more dimensions through diverse kernel functions. In previous studies, SVM has been successfully applied in neuroimaging to predict SAH and prognosis after treatment (16, 29). This study involved the development of seven different predictive models for HT, with the SVM model exhibiting superior performance. Moreover, we used SHAP, a game-theory-based method, to visually represent the impact of individual features in the combined model for prognostic forecasting. Our SHAP analysis identified several radiomic features with clinical relevance (30). Among the positively correlated features, elevated SmallAreaEmphasis values suggest densely clustered small-sized hyperdense foci within HIM regions. These micro-foci likely represent early microhemorrhages or ICE, and their aggregation indicates exacerbated blood–brain barrier disruption, which may predispose patients to HT following reperfusion therapy. Elongation exhibited a positive correlation with hemorrhage risk. HIM with an elongated shape predominantly localized to cortical regions may reflect the high metabolic demand of cortical tissue, which can amplify ischemia–reperfusion injury (8). Conversely, negatively correlated features revealed protective imaging patterns. Reduced DifferenceVariance signifies a homogenized grayscale distribution within HIM, suggesting potential alignment with infarct core regions characterized by severe blood–brain barrier damage, thereby elevating hemorrhage susceptibility (31). The negative correlations observed for ShortRunLowGrayLevelEmphasis and ShortRunEmphasis indicate diminished low-density areas and a lack of short-run texture structures within HIM. This further aligns with studies suggesting that lesions dominated by homogeneous hyperdensity are more prone to hemorrhagic complications (12).

The top nine predictive indicators in the combined model included clinical factors such as anterior circulation, sHIM, male gender, ASPECTS score, and HUmax ≥90. Morhard et al. (32) conducted a study comparing the incidence of HIM in patients with anterior and posterior circulation stroke using dual-energy CT. The results revealed that the incidence of HIM was significantly higher in patients with anterior circulation stroke than in those with posterior circulation stroke, leading to an increased rate of HT. This difference can be attributed to the relatively small volume of the posterior circulation area and limited resolution of the posterior fossa CT caused by beam-hardening artifacts. Manual ROI delineation may miss subtle HIM regions. Future work should prioritize automated segmentation tools to improve consistency. Kim et al. (33) discovered that the infarcted core was larger in the sHIM group, indicating a poor collateral circulation status, thereby leading to the expansion of the progressive infarcted core and an increased risk of HT (34). Besides, sHIM could result from arterial injury caused by the procedure, indicating the presence of significant intracranial hemorrhage. This finding aligns with those of the present study. The risk of HT after MT was comparable between men and women. However, in Asian populations, there is a stronger association between male gender and HT (35–37). Further evidence is required to confirm these observations. Numerous studies have demonstrated that a low preoperative ASPECTS may indicate severe damage to the blood–brain barrier due to ischemia (7, 38–41). This damage increases the risk of HT after MT. Additionally, the maximum CT value, especially when the threshold is set at 90 HU, indicates impairment of the blood–brain barrier and acts as a reliable indicator of the risk of HT (7, 8, 12, 42).

Our study has several limitations that should be addressed. First, the retrospective nature of the study, its single-center design, and its relatively small sample size may have introduced selection bias and overestimation of our findings. We are actively collaborating with external institutions to conduct prospective multicenter validation. Additionally, we plan to implement local prospective internal validation in ongoing studies. Until such external corroboration is achieved, clinical application of the model requires cautious interpretation, and its generalizability remains to be further verified. Second, the manual segmentation process for lesions in our study was time-consuming and complex, especially for lesions with indistinct boundaries. Consequently, future research endeavors should prioritize the improvement of automatic segmentation technology to achieve a satisfactory level of reliability and reproducibility. Third, the small sample size and the imbalance between the two groups in our study may have affected the statistical power and generalizability of the results. Larger sample sizes and improved balance between the groups should be considered in future studies. Fourth, our approach relies on conventional radiomics characteristics. The prediction ability of the model can be improved by incorporating a deep-learning approach. Finally, interpreting the radiomics features used in our study presents a formidable challenge to understanding their fundamental biological significance. To improve our understanding of their biological significance, future research should focus on investigating the association between the histological attributes of thrombi and radiomics features.

5 Conclusion

This study involved the development of a model that integrates radiomics characteristics obtained from NCCT images of HIM with clinical characteristics to evaluate the risk of HT in patients who have undergone MT. This reliable model can assist frontline doctors in identifying patients with a significantly increased risk of HT and provide support during clinical decision-making.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Fourth Affiliated Hospital, Zhejiang University College of Medicine (Approval number: KY2023112). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin.

Author contributions

SH: Conceptualization, Funding acquisition, Writing – original draft. JL: Formal analysis, Methodology, Validation, Writing – review & editing. JYH: Methodology, Writing – review & editing. YC: Formal analysis, Methodology, Validation, Writing – original draft. ZW: Investigation, Visualization, Writing – review & editing. JBH: Supervision, Writing – review & editing. SG: Resources, Validation, Writing – review & editing. XY: Data curation, Writing – review & editing. JF: Conceptualization, Data curation, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Science and Technology program of Jinhua Science and Technology Bureau (Grant no. 2021–4-184) and the Fundamental Research Funds for the Central Universities (Grant no. 226202400226).

Acknowledgments

Some of our experiments were carried out on the Onekey AI platform. Thank Onekey AI and its developers for help in this scientific research work.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    Yoshimura S Sakai N Yamagami H Uchida K Beppu M Toyoda K et al . Endovascular therapy for acute stroke with a large ischemic region. N Engl J Med. (2022) 386:1303–13. doi: 10.1056/NEJMoa2118191

  • 2.

    Powers WJ Rabinstein AA Ackerson T Adeoye OM Bambakidis NC Becker K et al . Guidelines for the early Management of Patients with Acute Ischemic Stroke: 2019 update to the 2018 guidelines for the early Management of Acute Ischemic Stroke: a guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke. (2019) 50:e344–418. doi: 10.1161/STR.0000000000000211

  • 3.

    van der Steen W van der Ende NAM Luijten SPR Rinkel LA van Kranendonk KR van Voorst H et al . Type of intracranial hemorrhage after endovascular stroke treatment: association with functional outcome. J NeuroIntervent Surg. (2023) 15:971–6. doi: 10.1136/jnis-2022-019474

  • 4.

    Zidan M Ghaei S Bode FJ Weller JM Krueger N Lehnen NC et al . Clinical significance and prevalence of subarachnoid hyperdensities on flat detector CT after mechanical thrombectomy: does it really matter?J NeuroIntervent Surg. (2024). 16:966–973. doi: 10.1136/jnis-2023-020661

  • 5.

    He J Fu F Zhang W Zhan Z Cheng Z . Prognostic significance of the clinical and radiological haemorrhagic transformation subtypes in acute ischaemic stroke: a systematic review and meta-analysis. Eur J Neurol. (2022) 29:3449–59. doi: 10.1111/ene.15482

  • 6.

    Hu S Hong J Liu F Wang Z Li N Wang S et al . An integrated nomogram combining clinical and radiomic features of hyperattenuated imaging markers to predict malignant cerebral edema following endovascular thrombectomy. Quant Imaging Med Surg. (2024) 14:4936–49. doi: 10.21037/qims-24-99

  • 7.

    Han N Zhang G Li Y Ma H Ge H Zhang X et al . Hyperattenuated lesions on immediate non-contrast CT after endovascular therapy predict intracranial hemorrhage in patients with acute ischemic stroke: a retrospective propensity matched study. Front Neurol. (2021) 12:664262. doi: 10.3389/fneur.2021.664262

  • 8.

    Kim J-M Park K-Y Lee WJ Byun JS Kim JK Park M-S et al . The cortical contrast accumulation from brain computed tomography after endovascular treatment predicts symptomatic hemorrhage. Eur J Neurol. (2015) 22:1453–8. doi: 10.1111/ene.12764

  • 9.

    Kau T Hauser M Obmann SM Niedermayer M Weber JR Hausegger KA . Flat detector Angio-CT following intra-arterial therapy of acute ischemic stroke: identification of hemorrhage and distinction from contrast accumulation due to blood-brain barrier disruption. AJNR Am J Neuroradiol. (2014) 35:1759–64. doi: 10.3174/ajnr.A4021

  • 10.

    Maïer B Brauner R Escalard S Gory B Lapergue B Sibon I et al . Association of Contrast Enhancement after Reperfusion with Outcomes According to blood pressure lowering in acute ischemic stroke patients. Neurology. (2022) 99:e2385–94. doi: 10.1212/WNL.0000000000201173

  • 11.

    Pinckaers FM Mentink MM Boogaarts HD van Zwam WH van Oostenbrugge RJ Postma AA . Early post-endovascular treatment contrast extravasation on dual-energy CT is associated with clinical and radiological stroke outcomes: a 10-year single-Centre experience. Eur Stroke J. (2023) 8:508–16. doi: 10.1177/23969873231157901

  • 12.

    Xu C Zhou Y Zhang R Chen Z Zhong W Gong X et al . Metallic hyperdensity sign on non-contrast CT immediately after mechanical Thrombectomy predicts parenchymal hemorrhage in patients with acute large-artery occlusion. AJNR Am J Neuroradiol. (2019) 40:661–7. doi: 10.3174/ajnr.A6008

  • 13.

    Yang S-J Lu Y-H Huang Y-C Chan L Ting W-Y . Immediate CT change after thrombectomy predicting symptomatic hemorrhagic transformation. Chin Med Associat. (2023) 86:854–8. doi: 10.1097/JCMA.0000000000000958

  • 14.

    Chen X Li Y Zhou Y Yang Y Yang J Pang P et al . CT-based radiomics for differentiating intracranial contrast extravasation from intraparenchymal haemorrhage after mechanical thrombectomy. Eur Radiol. (2022) 32:4771–9. doi: 10.1007/s00330-022-08541-9

  • 15.

    Ma Y Wang J Zhang H Li H Wang F Lv P et al . A CT-based radiomics nomogram for classification of intraparenchymal hyperdense areas in patients with acute ischemic stroke following mechanical thrombectomy treatment. Front Neurosci. (2023) 16:1061745. doi: 10.3389/fnins.2022.1061745

  • 16.

    Shan D Wang J Qi P Lu J Wang D . Non-contrasted CT Radiomics for SAH prognosis prediction. Bioengineering. (2023) 10:967. doi: 10.3390/bioengineering10080967

  • 17.

    Ren H Song H Wang J Xiong H Long B Gong M et al . A clinical–radiomics model based on non-contrast computed tomography to predict hemorrhagic transformation after stroke by machine learning: a multicenter study. Insights Imaging. (2023) 14:52. doi: 10.1186/s13244-023-01399-5

  • 18.

    Dekeyzer S Nikoubashman O Lutin B De Groote J Vancaester E De Blauwe S et al . Distinction between contrast staining and hemorrhage after endovascular stroke treatment: one CT is not enough. J NeuroIntervent Surg. (2017) 9:394–8. doi: 10.1136/neurintsurg-2016-012290

  • 19.

    Suzuki K Matsumaru Y Takeuchi M Morimoto M Kanazawa R Takayama Y et al . The impact of SAH finding on CT to the clinical outcome after mechanical thrombectomy for large vessel occlusion. J Neurol Sci. (2023) 453:120797. doi: 10.1016/j.jns.2023.120797

  • 20.

    Bernardo-Castro S Sousa JA Brás A Cecília C Rodrigues B Almendra L et al . Pathophysiology of blood–brain barrier permeability throughout the different stages of ischemic stroke and its implication on hemorrhagic transformation and recovery. Front Neurol. (2020) 11:594672. doi: 10.3389/fneur.2020.594672

  • 21.

    Renú A Laredo C Lopez-Rueda A Llull L Tudela R San-Roman L et al . Vessel Wall enhancement and blood–cerebrospinal fluid barrier disruption after mechanical Thrombectomy in acute ischemic stroke. Stroke. (2017) 48:651–7. doi: 10.1161/STROKEAHA.116.015648

  • 22.

    Chen Z Zhang Y Su Y Sun Y He Y Chen H . Contrast extravasation is predictive of poor clinical outcomes in patients undergoing endovascular therapy for acute ischemic stroke in the anterior circulation. J Stroke Cerebrovasc Dis. (2020) 29:104494. doi: 10.1016/j.jstrokecerebrovasdis.2019.104494

  • 23.

    Lee H Qureshi AM Mueller-Kronast NH Zaidat OO Froehler MT Liebeskind DS et al . Subarachnoid hemorrhage in mechanical Thrombectomy for acute ischemic stroke: analysis of the STRATIS registry, systematic review, and meta-analysis. Front Neurol. (2021) 12:663058. doi: 10.3389/fneur.2021.663058

  • 24.

    Qureshi AI Saleem MA Aytac E . Postprocedure subarachnoid hemorrhage after endovascular treatment for acute ischemic stroke. J Neuroimaging. (2017) 27:493–8. doi: 10.1111/jon.12430

  • 25.

    Kim DY Baik SH Jung C Kim JY Han S-G Kim BJ et al . Predictors and impact of Sulcal SAH after mechanical Thrombectomy in patients with isolated M2 occlusion. AJNR Am J Neuroradiol. (2022) 43:1292–8. doi: 10.3174/ajnr.A7594

  • 26.

    Cai J Zhou Y Zhao Y Xu C Yan S Ding X et al . Comparison of various reconstructions derived from dual-energy CT immediately after endovascular treatment of acute ischemic stroke in predicting hemorrhage. Eur Radiol. (2021) 31:4419–27. doi: 10.1007/s00330-020-07574-2

  • 27.

    Liu K Jiang L Ruan J Xia W Huang H Niu G et al . The role of dual energy CT in evaluating hemorrhagic complications at different stages after thrombectomy. Front Neurol. (2020) 11:583411. doi: 10.3389/fneur.2020.583411

  • 28.

    Sidey-Gibbons JAM Sidey-Gibbons CJ . Machine learning in medicine: a practical introduction. BMC Med Res Methodol. (2019) 19:64. doi: 10.1186/s12874-019-0681-4

  • 29.

    Zhang L Wu J Yu R Xu R Yang J Fan Q et al . Non-contrast CT radiomics and machine learning for outcomes prediction of patients with acute ischemic stroke receiving conventional treatment. Eur J Radiol. (2023) 165:110959. doi: 10.1016/j.ejrad.2023.110959

  • 30.

    Tomaszewski MR Gillies RJ . The biological meaning of Radiomic features. Radiology. (2021) 298:505–16. doi: 10.1148/radiol.2021202553

  • 31.

    Ma C Xu D Hui Q Gao X Peng M . Quantitative intracerebral iodine extravasation in risk stratification for intracranial hemorrhage in patients with acute ischemic stroke. AJNR Am J Neuroradiol. (2022) 43:1589–96. doi: 10.3174/ajnr.A7671

  • 32.

    Morhard D Ertl L Gerdsmeier-Petz W Ertl-Wagner B Schulte-Altedorneburg G . Dual-energy CT immediately after endovascular stroke intervention: prognostic implications. Cardiovasc Intervent Radiol. (2014) 37:1171–8. doi: 10.1007/s00270-013-0804-y

  • 33.

    Kim H Lee S-J Lee T-K Jung K-O . Subarachnoid contrast accumulation and Alberta stroke program early computed tomography score applied to contrast accumulation after Thrombectomy as predictors of symptomatic hemorrhage. World Neurosurg. (2020) 138:e847–58. doi: 10.1016/j.wneu.2020.03.102

  • 34.

    Lee S-J Do T Lee T-K . Postinterventional contrast accumulation early predicts malignant stroke in successfully recanalized patients with emergent large vessel occlusion. Interdiscip Neurosurg. (2022) 30:101621. doi: 10.1016/j.inat.2022.101621

  • 35.

    Wen L Zhang S Wan K Zhang H Zhang X . Risk factors of haemorrhagic transformation for acute ischaemic stroke in Chinese patients receiving intravenous thrombolysis. Medicine. (2020) 99:e18995. doi: 10.1097/MD.0000000000018995

  • 36.

    Iancu A Buleu F Chita DS Tutelca A Tudor R Brad S . Early hemorrhagic transformation after reperfusion therapy in patients with acute ischemic stroke: analysis of risk factors and predictors. Brain Sci. (2023) 13:840. doi: 10.3390/brainsci13050840

  • 37.

    Won Y-D Kim J-M Cheong J-H Ryu J-I Koh S-H Han M-H . Effect of possible osteoporosis on parenchymal-type hemorrhagic transformation in patients with cardioembolic stroke. JCM. (2021) 10:2526. doi: 10.3390/jcm10112526

  • 38.

    Ebashi R Ogata A Nishihara M Inoue K Yoshioka F Takase Y et al . Significance of simulated conventional images on dual energy CT after endovascular treatment for ischemic stroke. J NeuroIntervent Surg. (2019) 11:898–902. doi: 10.1136/neurintsurg-2018-014486

  • 39.

    Ahn S Mummareddy N Roth SG Jo J Bhamidipati A Ko Y et al . The clinical utility of dual-energy CT in post-thrombectomy care: part 1, predictors and outcomes of subarachnoid and intraparenchymal hemorrhage. J Stroke Cerebrovasc Dis. (2023) 32:107217. doi: 10.1016/j.jstrokecerebrovasdis.2023.107217

  • 40.

    Hao Z Yang C Xiang L Wu B Liu M . Risk factors for intracranial hemorrhage after mechanical thrombectomy: a systematic review and meta-analysis. Expert Rev Neurother. (2019) 19:927–35. doi: 10.1080/14737175.2019.1632191

  • 41.

    Sun J Lam C Christie L Blair C Li X Werdiger F et al . Risk factors of hemorrhagic transformation in acute ischaemic stroke: a systematic review and meta-analysis. Front Neurol. (2023) 14:1079205. doi: 10.3389/fneur.2023.1079205

  • 42.

    Yoon W Seo JJ Kim JK Cho KH Park JG Kang HK . Contrast enhancement and contrast extravasation on computed tomography after intra-arterial thrombolysis in patients with acute ischemic stroke. Stroke. (2004) 35:876–81. doi: 10.1161/01.STR.0000120726.69501.74

Summary

Keywords

hemorrhagic transformation, machine learning, thrombectomy, acute ischemic stroke, multi-detector CT

Citation

Hu S, Liu J, Hong J, Chen Y, Wang Z, Hu J, Gai S, Yu X and Fu J (2025) Machine learning, clinical-radiomics approach with HIM for hemorrhagic transformation prediction after thrombectomy and treatment. Front. Neurol. 16:1471274. doi: 10.3389/fneur.2025.1471274

Received

27 July 2024

Accepted

08 February 2025

Published

27 February 2025

Volume

16 - 2025

Edited by

Xindao Yin, Nanjing Medical University, China

Reviewed by

Zhenyu Shu, Zhejiang Provincial People’s Hospital, China

Duchang Zhai, Second Affiliated Hospital of Soochow University, China

Qianqian Gao, Wuxi People’s Hospital Affiliated to Nanjing Medical University, China

Updates

Copyright

*Correspondence: Jingjing Fu,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics