Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 11 February 2026

Sec. Precision Medicine

Volume 13 - 2026 | https://doi.org/10.3389/fmed.2026.1764733

This article is part of the Research TopicApplications and Advances of Artificial Intelligence in Medical Image Analysis: PET, SPECT/ CT, MRI, and Pathology ImagingView all 17 articles

Enhanced preoperative prediction for microvascular invasion in hepatocellular carcinoma through an optimized MR Radiomics combination strategy and machine learning predictor

Mengting FengMengting Feng1Yingjian YangYingjian Yang2Zongbo DaiZongbo Dai3Ziran Chen,Ziran Chen1,4Longyu Li,Longyu Li1,4Zewei Wu,Zewei Wu1,4Xuejian LiXuejian Li3Tingwei GuoTingwei Guo3Yiman MengYiman Meng3Qiang LiQiang Li5Zihao Zhao,Zihao Zhao1,4Tao Li,Tao Li1,4Jialin Zhang
Jialin Zhang3*Yan Kang,,,
Yan Kang1,4,6,7*
  • 1College of Medicine and Biological Information Engineering, Northeastern University, Shenyang, China
  • 2Department of Radiological Research and Development, Shenzhen Lanmage Medical Technology Co., Ltd., Shenzhen, Guangdong, China
  • 3Department of Hepatobiliary Surgery, The First Hospital of China Medical University, Shenyang, China
  • 4College of Health Science and Environmental Engineering, Shenzhen Technology University, Shenzhen, China
  • 5School of Data and Computer Science, Shandong Women’s University, Jinan, Shandong, China
  • 6College of Applied Sciences, Shenzhen University, Shenzhen, Guangdong, China
  • 7Faculty of Data Science, City University of Macau, Macao, Macao SAR, China

Background: Preoperative prediction of microvascular invasion (MVI) in hepatocellular carcinoma (HCC) is a crucial step toward personalized treatment, improved treatment outcomes, and enhanced patient survival. However, the disadvantage of existing prediction models of MVI in HCC based on enhanced magnetic resonance imaging (MRI) is that they require combining non-imaging information from enhanced MRI, or determining the perioperative region is highly subjective. These disadvantages are not conducive to the clinical application of predictive models, which hinders clinical decision-making and management for these vulnerable populations.

Methods: To address the problem of combining non-imaging information from enhanced MRI with the highly subjective determination of the perioperative region, we propose an enhanced preoperative prediction model for MVI in HCC using an optimized MR Radiomics combination strategy and a machine learning predictor. First, the HCC was manually segmented from 125 × 512 × 512 × N abdominal enhanced T1-weighted magnetic resonance imaging (T1WI) images during the arterial phase, generating 125 × 512 × 512 × N HCC mask images. Second, 125 × 1,692 MR Radiomics features of HCC are extracted from abdominal enhanced T1WI images based on the HCC mask images. Third, the 125 × N selected and 125 × 10 fused MR Radiomics features are determined using the proposed optimized MR Radiomics combination strategy with 5-fold cross-validation. Finally, the best preoperative prediction model is constructed using a random forest (RF) predictor with 125 × N selected and 125 × 10 fused MR Radiomics features.

Results: The proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) achieves a mean accuracy of 0.7520 ± 0.0867, a mean precision of 0.7354 ± 0.1863, a mean recall of 0.6955 ± 0.2203, a mean F1-score of 0.6943 ± 0.1437, and a mean AUC of 0.7962 ± 0.1700.

Discussion: The proposed best preoperative prediction model can effectively predict MVI in HCC, potentially serving as a strong decision-making tool for these vulnerable populations.

1 Introduction

Hepatocellular carcinoma (HCC) has become the most common type of primary liver cancer, accounting for approximately 85–90% of all cases, resulting in the fourth leading cause of cancer-related deaths globally, with China accounting for 42% of new cases annually in 2022 and a 5-year survival rate of only 18% for advanced patients (13). Microvascular invasion (MVI) is a significant risk factor for postoperative recurrence and metastasis of HCC (4).

Preoperative prediction of MVI in HCC is a crucial step toward personalized treatment, improved treatment outcomes, and enhanced patient survival (47). It also plays an irreplaceable role in optimizing clinical decision-making and improving patient prognosis (8). First, for MVI-positive patients, preoperative prediction can prompt doctors to expand the surgical margin and use a wider resection range (e.g., wide-margin resection) to reduce the risk of tumor cell residue and recurrence (9, 10). Conversely, for MVI-negative patients, the surgical margin can be appropriately reduced to minimize surgical trauma. Second, MVI status can help determine whether more precise surgical methods, such as anatomical liver resection, are needed (10, 11). For MVI-positive patients, anatomical liver resection may be more beneficial for thoroughly clearing the tumor and possible micro-metastases in the surrounding area. Third, patients with a high risk of MVI predicted before surgery may be more suitable for neoadjuvant therapy, such as chemotherapy and targeted therapy, to reduce tumor volume, decrease the degree of MVI, improve the surgical resection rate, and achieve a curative effect (12). Fourth, postoperative adjuvant treatment plans can be formulated in advance based on the preoperative MVI prediction results (13). For MVI-positive or high-risk patients, more aggressive adjuvant therapy, such as hepatic artery infusion chemotherapy and targeted therapy, may be needed after surgery to further reduce the risk of recurrence. In addition, for MVI-negative patients, the intensity and duration of adjuvant therapy can be adjusted appropriately. Finally, MVI is an essential predictor of postoperative recurrence and metastasis in hepatocellular carcinoma (4). Preoperative prediction of MVI status can help doctors more accurately evaluate patient prognosis and provide more reasonable treatment recommendations and follow-up plans. For example, patients with positive MVI have a higher risk of postoperative recurrence and require closer follow-up and monitoring.

Compared with abdominal ultrasound imaging and computed tomography (CT), enhanced T1-weighted magnetic resonance imaging (T1WI) is considered the preferred imaging modality for preoperative prediction of MVI in HCC. Compared with abdominal ultrasound, CT and magnetic resonance imaging (MRI) offer higher spatial resolution and stronger lesion guidance for deep tumors, such as HCC in the liver (14). Compared with abdominal CT imaging, MRI also offers high soft-tissue resolution and is superior to CT and ultrasound in detecting small lesions in solid organs, such as the liver (e.g., liver cancer ≤1 cm) (15). It can also identify tumor boundaries and adjacent nerves and blood vessels, helping detect early liver cancer and determine the relationship between tumors and surrounding tissues. Meanwhile, MRI avoids X-ray exposure and is suitable for children, pregnant women, and patients who require multiple examinations, especially for liver cancer patients who require long-term follow-up, which is safer (16, 17). Compared with enhanced CT, enhanced T1WI can effectively display the anatomical structures of liver organs, and this comparison helps doctors observe their morphology and position (18). Meanwhile, liver tumor lesions, such as HCC, also exhibit different signal characteristics on enhanced T1WI images due to changes in tissue characteristics, indicating that the tumor may be invasive and invade surrounding liver tissue. Therefore, T1WI is an essential sequence in MRI of HCC, providing key information for diagnosis, staging, treatment plan selection, and efficacy evaluation. Finally, selecting abdominal T1WI images acquired during the late arterial phase is crucial for preoperative prediction of MVI in HCC. Specifically, during the arterial phase, mainly the late arterial phase, of dynamic enhanced MRI, liver tumors show uniform or uneven significant enhancement. In contrast, their enhancement in the portal vein phase and/or delayed phase is lower than that of the liver parenchyma (1921). Therefore, the arterial phase of liver tumors shows significant enhancement, while the portal vein phase shows less enhancement than the liver parenchyma.

Currently, the two major mainstream approaches to prediction tasks are machine learning and deep learning (2239). However, deep learning-based prediction tasks are suitable for scenarios with large amount of data, complex problems, and high performance requirements (2224). It often requires more training data during the training process of the prediction network so that the prediction model can achieve a specific performance. Although some scholars have developed transfer learning strategies to overcome the aforementioned technical issues, this remains an insurmountable problem in medical imaging, where training data are scarce. Unlike deep learning, machine learning algorithms have advantages in handling small data sets and in computational resource requirements, making them suitable for scenarios with limited data or computational resources (3739).

Radiomics, combined with machine learning algorithms, a technique for extracting and analyzing high-throughput quantitative features from medical images, has been widely used for auxiliary diagnosis, treatment, and prognosis (2436) and has also opened new avenues for non-invasive preoperative prediction of MVI in HCC based on enhanced MRI (3739). However, the disadvantage of existing prediction models of MVI in HCC based on enhanced MRI is that they require combining non-imaging information from enhanced MRI, such as aminotransferase-to-platelet ratio and gamma-glutamyl transferase-to-platelet ratio (22, 23, 38), or that determining the perioperative region is highly subjective (39). These disadvantages are not conducive to the clinical application of predictive models. Radiomics derived from abdominal enhanced T1WI images during the arterial phase should provide more information to enhance preoperative prediction of MVI in HCC. Therefore, it is necessary to propose a preoperative prediction model for MVI in HCC to address the challenge of combining non-imaging information from enhanced MRI with highly subjective perioperative region delineation. Our contributions in this study are briefly described as follows:

1. We propose a preoperative prediction model for MVI in HCC to eliminate non-imaging information of enhanced MRI and the highly subjective determination of perioperative region, which is conducive to the clinical application of the predictive model.

2. The linear and non-linear features extracted from the MR Radiomics features of HCC are fully exploited by the feature selection and feature fusion algorithms to adapt to different machine learning predictors.

3. The best preoperative prediction model is determined by the proposed optimized MR Radiomics combination strategy to seek the optimal combination of features.

4. The proposed best preoperative prediction model can effectively predict the MVI in HCC, which may become a strong decision-making tool for these vulnerable populations.

2 Materials and methods

This study proposes an enhanced preoperative prediction model for MVI in HCC using an optimized MR Radiomics combination strategy and a machine learning predictor. Based on the above, materials and methods are described in Sections 2.1 and 2.2, respectively.

2.1 Materials

Figure 1 shows the flow diagram for patient selection in this study. Table 1 reports the characteristics of the dataset of the 125 patients. Specifically, all 125 patients underwent abdominal enhanced MRI scanning using the same 3.0 T MRI scanner (GE MEDICAL SYSTEMS, SIGNA Pioneer, Scanning Sequence: GR, Sequence Variant: SS/SK) before surgery. The contrast agent for enhanced MRI scanning is gadopentetamide (Gd-DTPA, Manufacturer: Beijing Beilu Pharmaceutical), which is injected intravenously at a dose of 0.2 mL/kg at a flow rate of 1.0 mL/s and then flushed with 20 mL of physiological saline. Abdominal enhanced MR images at T1-weighted imaging were collected at 30 s (arterial phase), 60 s (portal phase), and 120 s (delayed phase) after injection. After the surgery, two pathologists with ≥10 years of experience double-blindly reviewed the patients’ pathological sections to evaluate the presence of MVI.

Figure 1
Flowchart detailing patient selection for a study at the First Hospital of China Medical University with 436 hepatocellular carcinoma patients. Inclusion and exclusion criteria are listed, resulting in 125 eligible patients for analysis of preoperative MVI prediction in HCC.

Figure 1. Flow diagram for patient selection in this study.

Table 1
www.frontiersin.org

Table 1. Characteristics of the dataset of the 125 patients.

Based on the above findings, 125 preoperative cases with enhanced abdominal T1WI images at the late arterial phase were retrospectively collected, and each case had an explicit postoperative MVI diagnosis (positive or negative MVI). The Medical Ethics Committee of the First Hospital of China Medical University approved this study.

2.2 Methods

Figure 2 shows the construction process of the enhanced preoperative prediction model for MVI in HCC using an optimized MR Radiomics combination strategy and the random forest (RF) predictor. First, the HCC was manually segmented from 125 × 512 × 512 × N abdominal enhanced T1WI images during the arterial phase, generating 125 × 512 × 512 × N HCC mask images. Second, 125 × 1,692 MR Radiomics features of HCC are extracted from abdominal enhanced T1WI images based on the HCC mask images. Third, the 125 × N selected and 125 × 10 fused MR Radiomics features are determined using the proposed optimized MR Radiomics combination strategy with 5-fold cross-validation. Finally, the best preoperative prediction model is constructed using the RF predictor with 125 × N selected and 125 × 10 fused MR Radiomics features.

Figure 2
Flowchart illustrating a four-step pipeline for predicting microvascular invasion (MVI) in hepatocellular carcinoma from abdominal MR images. The steps shown are tumor segmentation, MR radiomics feature extraction, radiomics combination and optimization using LASSO and SPECTRAL, and final MVI prediction using a random forest predictor, with outcomes labeled as positive or negative MVI.

Figure 2. Construction process of the enhanced preoperative prediction model for microvascular invasion (MVI) in hepatocellular carcinoma (HCC) using an optimized MR Radiomics combination strategy and the random forest (RF) predictor.

2.2.1 Step 1: HCC segmentation

The HCC was manually segmented from 125 × 512 × 512 × N abdominal enhanced T1WI images during the arterial phase, generating 125 × 512 × 512 × N HCC mask images. Specifically, two radiologists with ≥5 years of experience used ITK-SNAP software (v3.8.0, http://www.itksnap.org) on 125 × 512 × 512 × N abdominal enhanced T1WI images during the arterial phase to independently manually delineate the region of interest of the HCC and reach consensus through negotiation (intra-group correlation coefficient ICC >0.85). If there is no consensus through negotiation, then an experienced chief radiologist reviews and modifies the 125 × 512 × 512 × N HCC mask images.

2.2.2 Step 2: MR Radiomics feature extraction

125 × 1,692 MR Radiomics features of HCC are extracted from abdominal enhanced T1WI images based on the HCC mask images. Specifically, PyRadiomics (2436), a radiomics feature extraction model, is used to extract 125 × 1,692 MR Radiomics features of HCC from abdominal enhanced T1WI images using HCC mask images. First, six kinds of derived HCC images are generated using the Laplacian of Gaussian (LoG), wavelet, square, square root, logarithm, and exponential filters, respectively, based on the original HCC images. Then, 125 × 1,692 MR Radiomics features are extracted from the original and six derived HCC images based on the classes of radiomic features. The sigma values for LoG are 1.0, 2.0, 3.0, 4.0, and 5.0, and the wavelet filter includes the nine component combinations: LLH, LHL, LHH, HLL, HLH, HHL, HHH, and LLL. The classes of radiomic features include shape-based (2D and 3D) (shape), first-order statistics (first order), gray-level cooccurrence matrix (GLCM), gray-level run length matrix (GLRLM), gray-level size zone matrix (GLSZM), neighboring gray tone difference matrix (NGTDM), and gray-level dependence matrix (GLDM). These 125 × 1,692 MR Radiomics features of HCC are available in Supplementary Table S1.

2.2.3 Step 3: Optimized MR Radiomics combination strategy

The 125 × N selected and 125 × 10 fused MR Radiomics features are determined using the proposed optimized MR Radiomics combination strategy and 5-fold cross-validation. Specifically, the proposed optimized MR Radiomics combination strategy consists of three key steps based on the 5-fold cross-validation. First, the least absolute shrinkage and selection operator (LASSO) (2436, 38), a powerful method for variable selection and regularization in high-dimensional datasets with multicollinearity or a much larger number of variables than the sample size, is performed at each fold cross-validation (25 × 1,692 MR Radiomics features) to determine the 125 × N selected MR Radiomics features. Second, spectral embedding (SPECTRAL) (40, 41), a feature dimensionality reduction method based on graph theory, is performed at each fold cross-validation (the same 25 × 1,692 MR Radiomics features need to be selected by the LASSO) to determine the 125 × 10 fused MR Radiomics features. SPECTRAL constructs an adjacency graph or k-nearest neighbor graph based on the similarity between data points of MR Radiomics features in each fold of cross-validation, calculates the Laplacian matrix of the graph, and finally solves the eigenvectors of the matrix to obtain low-dimensional embeddings. Therefore, the core of SPECTRAL is to preserve the connectivity of data points in MR Radiomics features within the graph, which is very effective for discovering non-convex, streamlined cluster structures. Finally, the 25 × N selected and 25 × 10 fused MR Radiomics features at each fold cross-validation are combined to obtain the optimized MR Radiomics combination vector.

2.2.4 Step 4: Microvascular invasion prediction

RF (2636, 38), an ensemble learning method that improves a model’s generalization and robustness by combining predictions from multiple decision trees, is selected as the predictor for MVI in HCC. Specifically, sub-datasets from the optimized MR Radiomics combination vector are randomly selected with replacement, and different training samples are generated for each tree. Then, during node splitting in each decision tree, a portion of the MR Radiomics features in the sub-dataset is randomly selected to find the optimal splitting point and reduce feature correlation. Secondly, each tree grows independently on the sub-dataset until the stopping condition is met. Finally, the MVI status is predicted using majority voting.

3 Experiments and results

This section comprehensively implements the comparative experiment. Then, the preoperative prediction models based on different predictors and MR Radiomics features are presented.

3.1 Experiments

Figure 3 shows the comparative experimental design for constructing and evaluating the MVI preoperative prediction models. Specifically, nine classic machine learning predictors (2636, 38) are selected to perform the comparative experiment with different MR Radiomics features, including RF, multilayer perceptron (MLP), k-nearest neighbor (KNN), support vector machine (SVM), logistic regression (LR), decision tree (DT), gradient boosting decision tree (GBDT), linear discriminant analysis (LDA), and AdaBoost (Ada). Meanwhile, LASSO is selected as the feature selection algorithm (2636). In addition, principal component analysis (PCA) (34), independent component analysis (ICA) (42), isometric mapping (ISOMAP) (43), uniform manifold approximation and projection (UMAP) (44), and SPECTRAL are selected as the feature fusion algorithms.

Figure 3
Flowchart illustrating three experiments for radiomics feature analysis. Experiment one uses 125 by 1692 radiomics features with five-fold validation and Random Forest for classification. Experiment two applies five-fold validation with LASSO, ICA, PCA, ISOMAP, UMAP, and SPECTRAL for feature selection or dimensionality reduction, followed by MLP, KNN, SVM, LR, and DT predictors. Experiment three combines LASSO with each dimensionality reduction technique, uses five-fold validation, and applies GBDT, LDA, and AdaBoost. Color-coded circles represent positive (purple) and negative (green) MVI outcomes.

Figure 3. Experimental design for constructing and evaluating the microvascular invasion (MVI) preoperative prediction models.

Experiment 1 is designed to evaluate preoperative prediction models based on the nine classic machine learning predictors and 125 × 1,692 Radiomics features using 5-fold cross-validation (training set: validation set = 4:1). Then, Experiment 1 is designed to evaluate the preoperative prediction models based on nine classic machine learning predictors and 125 × N selected radiomics features generated by the LASSO/125 × M (M = 10,20,30,40,50,60), fused radiomics features generated by PCA/ICA/UMAP/SPECTRAL with the 5-fold cross-validation. Finally, Experiment 3 is designed to evaluate the preoperative prediction models based on the nine classic machine learning predictors and 125 × (N + M) selected and fused radiomics features with the 5-fold cross-validation. A standard z-score normalization is applied to these radiomics features before training the nine machine learning predictors. To avoid data leakage, the LASSO algorithm is applied to the training set separately to select the radiomics features. Then, on the validation set, radiomics features consistent with the selected radiomics features from the training set are selected by indexing. Table 2 reports the definition of the feature selection and fused algorithm in Python 3.8. Meanwhile, Table 3 reports the definition of the nine machine learning predictors in Python 3.8.

Table 2
www.frontiersin.org

Table 2. Definition of the feature selection and fused algorithm in Python 3.8.

Table 3
www.frontiersin.org

Table 3. Definition of the nine machine learning predictors in Python 3.8.

3.1.1 Evaluation metrics

To assess performance differences between the preoperative prediction models, five standard evaluation metrics used in this study include the accuracy, precision, recall, F1-score, and area under the curve (AUC) (2436). Specifically, the AUC can be calculated by the receiver operating characteristic (ROC) curve. In addition, the accuracy, precision, recall, and F1-scores are defined by Equations 1–4:

Accuracy = TP + TN TP + TN + FP + FN ,     (1)
Precision = TP TP + FP ,     (2)
Recall = TP TP + FN ,     (3)
F 1 score = 2 × Precision × Recall Precision + Recall ,     (4)

where TP and TN represent the number of positive MVIs predicted to be positive and the number of negative MVIs predicted to be negative. In contrast, FP and FN represent the number of positive MVIs predicted to be negative and the number of negative MVIs predicted to be positive.

Based on Equations 1–4 and the AUC, Equation 5 is determined.

Mean evaluation metric = ( mean accuracy + mean precision + mean recall + mean F 1 score + mean AUC ) / 5 ,     (5)

where mean accuracy, mean precision, mean recall, mean F1-score, and mean AUC separately represent the average of the accuracy, precision, recall, F1-score, and AUC.

3.2 Results

This section reports the results of the comparative experiment to highlight the performance of the proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10). Overall, the proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) achieves a mean accuracy of 0.7520 ± 0.0867, a mean precision of 0.7354 ± 0.1863, a mean recall of 0.6955 ± 0.2203, a mean F1-score of 0.6943 ± 0.1437, and a mean AUC of 0.7962 ± 0.1700.

3.2.1 Comparative results of the MVI preoperative prediction models in Experiment 1

Figure 4 shows the visualized results of the MVI preoperative prediction models constructed by different machine learning predictors and 1,692 MR Radiomics features. Specifically, the performance of these preoperative prediction models in Experiment 1 is reported in Table A1. Figures 4ae show a mean accuracy, a mean precision, a mean recall, a mean F1-score, and a mean AUC of the MVI preoperative prediction models constructed by different machine learning predictors and 1,692 MR Radiomics features. Figure 4f shows the ROC curves of the MVI preoperative prediction models constructed by different machine learning predictors and 125 × 1,692 MR Radiomics features.

Figure 4
Composite figure presenting performance comparisons of nine machine learning predictors using bar charts and a line graph. Panels a to e show mean accuracy, mean precision, mean recall, mean F1 score, and mean AUC, respectively, with values labeled atop bars for RF, MLP, KNN, SVM, LR, DT, GBDT, LDA, and Ada predictors. Panel f displays a multi-line ROC curve plot comparing predictors, with AUC values in the legend. Panel g compares the mean evaluation metric for each predictor, highlighting LR as top-ranked. Panel h details LR’s mean scores across five metrics using color-coded bars.

Figure 4. Visualized results of the microvascular invasion (MVI) preoperative prediction models constructed by different machine learning predictors and 125 × 1,692 MR Radiomics features: (a) mean accuracy; (b) mean precision; (c) mean recall; (d) mean F1-score; (e) mean AUC; (f) ROC curves of the MVI preoperative prediction models constructed by different machine learning predictors and 125 × 1,692 MR Radiomics features; (g) model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × 1,692 MR Radiomics features); and (h) mean evaluation metrics of the best performance model (LR + 125 × 1,692 MR Radiomics features).

Overall, Figures 4g,h show that the MVI preoperative prediction model based on the LR and 125 × 1,692 MR Radiomics features performs the best, achieving a mean accuracy of 0.6720 ± 0.0867, a mean precision of 0.6122 ± 0.1615, a mean recall of 0.6840 ± 0.1088, a mean F1-score of 0.6387 ± 0.1151, and a mean AUC of 0.7481 ± 0.1087. In addition, compared with the mean evaluation metric of other MVI preoperative prediction models constructed by RF, MLP, KNN, SVM, DT, GBDT, LDA, and Ada and 125 × 1,692 MR Radiomics features, the best MVI preoperative prediction model improves by 3.35, 4.63, 14.67, 9.86, 16.64, 4.71, 9.02, and 2.76%, respectively.

3.2.2 Comparative results of the MVI preoperative prediction models in Experiment 2

Figure 5 shows the visualized results of the MVI preoperative prediction models constructed by different machine learning predictors with 125 × N selected MR Radiomics features or 125 × M fused MR Radiomics features. Specifically, the performance of these preoperative prediction models in Experiment 2 is reported in Table A2. Figures 5a,b,e show the model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × N selected MR Radiomics features), the mean evaluation metrics of the best performance model (RF + LASSO), and their ROC curves. Meanwhile, Figures 5c,d,f show the model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × M fused MR Radiomics features), the mean evaluation metrics of the best performance model (Ada + UMAP-50), and their ROC curves.

Figure 5
Multi-panel scientific figure presents machine learning predictors and dimensionality reduction performance for a classification task. Panels a and c display bar charts ranking methods by mean evaluation metrics, with RF+LASSO and Ada+UMAP-50 leading their respective groups. Panels b and d show detailed bar charts with five evaluation metrics for RF+LASSO and Ada+UMAP-50, including accuracy, precision, recall, F1 score, and AUC. Panels e and f depict receiver operating characteristic curves comparing predictors’ true positive and false positive rates with reported AUC values.

Figure 5. Visualized results of the microvascular invasion (MVI) preoperative prediction models constructed by different machine learning predictors with 125 × N selected MR Radiomics features, or 125 × M fused MR Radiomics features. (a) Model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × N selected MR Radiomics features); (b) mean evaluation metrics of the best performance model (RF + LASSO); (c) top five model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × M fused MR Radiomics features); (d) mean evaluation metrics of the best performance model (Ada + UMAP-50); (e) ROC curves of the MVI preoperative prediction models constructed by different machine learning predictors and 125 × N selected MR Radiomics features; and (f) ROC curves of the top five MVI preoperative prediction models constructed by different machine learning predictors and 125 × M fused MR Radiomics features.

First, the MVI preoperative prediction model based on the RF and 125 × N selected MR Radiomics features using the LASSO algorithm (RF + LASSO) performs the best among all machine learning predictors using the LASSO algorithm, achieving a mean accuracy of 0.7280 ± 0.0769, a mean precision of 0.6965 ± 0.1481, a mean recall of 0.6788 ± 0.2183, a mean F1-score of 0.6697 ± 0.1324, and a mean AUC of 0.7851 ± 0.1363. In addition, compared with the mean evaluation metric of other MVI preoperative prediction models constructed by MLP, KNN, SVM, LR, DT, GBDT, LDA, and Ada and 125 × N selected radiomics features, the best MVI preoperative prediction model improves by 5.64, 8.65, 2.88, 4.28, 10.76, 4.48, 8.21, and 4.59%, respectively. Compared to the mean evaluation metric of Experiment 1, the mean evaluation metric of the MVI preoperative prediction models constructed by RF, MLP, KNN, SVM, LR, DT, GBDT, LDA, and Ada and 125 × N selected radiomics features in Experiment 2 improves by 7.41, 3.05, 10.09, 11.04, −0.22%, 9.94, 4.29, 4.87, and 2.23%, respectively.

Second, the MVI preoperative prediction model based on the Ada and 125 × 50 fused MR Radiomics features using the UMAP algorithm (Ada + UMAP-50) performs the best among all machine learning predictors using the ICA, PCA, ISOMAP, UMAP, and SPECTRAL algorithms, achieving a mean accuracy of 0.7200 ± 0.0894, a mean precision of 0.6930 ± 0.1795, a mean recall of 0.6516 ± 0.1992, a mean F1-score of 0.6546 ± 0.1434, and a mean AUC of 0.6735 ± 0.1145. In addition, compared to the mean evaluation metric of other best MVI preoperative prediction models among all machine learning predictors using each feature fusion algorithm, LDA + PCA-20, LDA + ICA-20, MLP + ISOMAP-10, and LR + SPECTRAL-10, mean evaluation metric of Ada + UMAP-50 model improves by 1.03, 1.31, 5.51, and 5.38%, respectively. Compared to the mean evaluation metric of the MVI preoperative prediction model based on Ada and 125 × 1,692 MR Radiomics features in Experiment 1, the mean evaluation metric of Ada + UMAP-50 model in Experiment 2 improves by 3.15%. Compared with the mean evaluation metric of the MVI preoperative prediction model based on LDA and 125 × 1,692 MR Radiomics features in Experiment 1, the mean evaluation metric of the LDA + PCA-20 and LDA + ICA-20 models in Experiment 2 improves by 8.74 and 8.46%, respectively. Compared to the mean evaluation metric of the MVI preoperative prediction model based on MLP and 125 × 1,692 MR Radiomics features in Experiment 1, the mean evaluation metric of MLP + ISOMAP-10 model in Experiment 2 improves by 0.23%.

3.2.3 Comparative results of the MVI preoperative prediction models in Experiment 3

Figure 6 shows the visualized results of the MVI preoperative prediction models constructed by different machine learning predictors with 125 × N selected MR Radiomics features and 125 × M fused MR Radiomics features. Specifically, the performance of these preoperative prediction models in Experiment 3 is reported in Table A3. Figures 6ad show the model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × N selected MR Radiomics features + 125 × M fused MR Radiomics features), the mean evaluation metrics of the best performance model (KNN + LASSO + UMAP-50), their ROC curves, and confusion matrix of the best performance model.

Figure 6
Four-panel figure compares machine learning model results. Panel a is a bar chart showing mean evaluation metrics for five model combinations, with RF+LASSO+SPECTRAL-10 performing best. Panel b is a bar chart with five performance metrics for RF+LASSO+SPECTRAL-10, including accuracy, precision, recall, F1 score, and AUC, with mean AUC highest. Panel c is a line chart displaying ROC curves for the five models, with RF+LASSO+ICA-10 having the highest AUC. Panel d is a confusion matrix for RF+LASSO+SPECTRAL-10, showing true and predicted labels with values 0.803, 0.197, 0.315, and 0.685.

Figure 6. Visualized results of the microvascular invasion (MVI) preoperative prediction models constructed by different machine learning predictors with 125 × N selected MR Radiomics features and 125 × M fused MR Radiomics features. (a) Top five model performance ranking based on the mean evaluation metrics (machine learning predictor + 125 × N selected MR Radiomics features + 125 × M fused MR Radiomics features); (b) mean evaluation metrics of the best performance model (RF + LASSO + SPECTRAL-10); (c) ROC curves of the top five MVI preoperative prediction models constructed by different machine learning predictors with 125 × N selected MR Radiomics features and 125 × M fused MR Radiomics features; and (d) confusion matrix of the best performance model (RF + LASSO + SPECTRAL-10).

The MVI preoperative prediction model based on the RF with 125 × N selected MR Radiomics features and 125 × 10 fused MR Radiomics features (RF + LASSO + ICA-10) performs the best among all machine learning predictors using LASSO and ICA, PCA, ISOMAP, UMAP, and SPECTRAL algorithms, achieving a mean accuracy of 0.7520 ± 0.0867, a mean precision of 0.7354 ± 0.1863, a mean recall of 0.6955 ± 0.2203, a mean F1-score of 0.6943 ± 0.1437, and a mean AUC of 0.7962 ± 0.1700. In addition, compared to the mean evaluation metric of other best MVI preoperative prediction models among all machine learning predictors using the LASSO and each feature fusion algorithm, KNN + LASSO + UMAP-50, RF + LASSO + ICA-10, RF + LASSO + PCA-20, and LDA + LASSO + ISOMAP-40, mean evaluation metric of RF + LASSO + ICA-10 model improves by 0.94, 1.43, 1.57, and 2.33%, respectively. Compared to the mean evaluation metric of RF + LASSO, the mean evaluation metric of RF + LASSO + ICA-10, RF + LASSO + ICA-10, and RF + LASSO + PCA-20 models improves by 2.31, 0.87, and 0.73%, respectively. Compared with the mean evaluation metric for KNN + LASSO, the mean evaluation metric of KNN + LASSO + UMAP-50 model improves by 10.02%. Compared to the mean evaluation metric of LDA + LASSO, the mean evaluation metric of LDA + LASSO + ISOMAP-40 model improves by 8.19%.

3.2.4 DeLong’s test

To compare whether there are significant differences in the areas under the two ROC curves of the top five models in Experiment 3, DeLong’s test (45, 46), a non-parametric statistical method, is performed using Statistical Product and Service Solutions (SPSS, launched by IBM, version 27).

Table 4 reports the paired-sample area difference under the ROC curves based on DeLong’s test. Specifically, compared to the significance among the MVI preoperative prediction models constructed by RF vs. RF + LASSO, the significance among the MVI preoperative prediction models constructed by RF + LASSO vs. RF + LASSO + PECTRAL-10/RF vs. RF + LASSO + ICA-10 improves by 0.30%/0.50%. Compared to the significance among the MVI preoperative prediction models constructed by KNN vs. KNN + LASSO, the significance among the MVI preoperative prediction models constructed by KNN vs. KNN + LASSO + UMAP-50 improves by 5.7%. Compared to the significance among the MVI preoperative prediction models constructed by LDA vs. LDA + LASSO, the significance among the MVI preoperative prediction models constructed by LDA vs. LDA + LASSO + ISOMAP-40 improves by 2.2%.

Table 4
www.frontiersin.org

Table 4. Paired-sample area difference under the ROC curves based on DeLong’s test.

4 Discussion

This section conducts the following discussions based on the experimental results. In addition, this section outlines the limitations of this study and its future direction.

4.1 RF predictor

The core advantage of random forest lies in its high robustness and generalization ability, enabled by dual randomness, making it the mainstream algorithm for processing small- to medium-sized data (47), such as the 125 cases of MR Radiomics features. First, by using bootstrap sampling and random selection of MR Radiomics features, inter-tree correlation can be reduced, thereby helping the MVI preoperative prediction model avoid excessive dependence on noise. Second, independent generation across trees supports distributed computing and yields significantly faster training than serial algorithms. Finally, by using self-sampling and majority voting, the problem of category bias can be alleviated, which is particularly important for constructing the MVI preoperative prediction model.

4.2 SPECTRAL

The advantage of SPECTRAL lies in capturing the intrinsic topological structure of data via the spectral decomposition of the Turalaplus matrix, which is particularly suitable for low-dimensional embedding of non-linear, high-dimensional data (40, 41), such as the 25 × 1,692 MR Radiomics features of HCC in each cross-validation fold. First, a similarity matrix (e.g., a Gaussian kernel) is constructed using graph theory, and the 25 × 1,692 MR Radiomics features of HCC are mapped to a low-dimensional space via eigenvalue decomposition, overcoming the limitations of linear dimensionality reduction. In addition, it takes into account both the local neighborhood relationships, as captured by the similarity matrix, and the global topology of the Laplacian matrix, as encoded in the feature vector, to capture the overall connectivity.

4.3 Optimized MR Radiomics combination strategy

The proposed optimized MR Radiomics combination strategy is similar to neural architecture search (48, 49), which determines the optimal network architecture for classification and segmentation tasks. The proposed optimized MR Radiomics combination strategy combines linearly selected MR Radiomics features and non-linearly fused MR Radiomics features and concatenates them to generate an optimized MR Radiomics combination vector that matches the predictor. Since the optimized MR Radiomics combination vector includes both linearly selected MR Radiomics features and non-linearly fused MR Radiomics features, it facilitates RF predictor sampling and feature random selection, thereby reducing inter-tree correlation and avoiding excessive dependence on noise in the MVI preoperative prediction model. In addition, non-linear MR Radiomics features are fused from all MR Radiomics of each case to avoid the problem that LASSO’s selection may miss important non-linear MR Radiomics features.

4.4 Proposed MVI preoperative prediction model

The proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) only needs to extract MR features from abdominal enhanced T1WI images during the artistic phase, without the need for T2-weighted, diffusion-weighted, or other types of images, and without the need to collect non-imaging information of enhanced MRI (such as aminotransferase-to-platelet ratio and gamma-glutamyl transferase-to-platelet ratio) (38). The proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) also avoids the highly subjective determination of the perioperative region. Therefore, the proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) simplifies its clinical application. In addition, compared with the previous contrast agent, Gd-EOB-DTPA (39), the price of the Gd-DTPA contrast agents used in this study is lower, thereby reducing the financial burden of the disease. Meanwhile, the arterial phase images of Gd-DTPA contrast agents are relatively less susceptible to respiratory motion artifacts and maintain stable quality.

4.5 Limitations and future studies

Although we propose an MVI preoperative prediction model for MVI in HCC to eliminate non-imaging information of enhanced MRI and the highly subjective determination of perioperative region, HCC mask images are manually segmented by radiologists. This may hinder the clinical application of the MVI preoperative prediction model. Therefore, the subsequent research aims to develop an HCC automatic segmentation model to fully automate the prediction process. In addition, we do not have sufficient abdominal enhanced MRI images to further validate the proposed MVI preoperative prediction model’s performance. Therefore, we encourage researchers to collect more abdominal-enhanced MRI images to validate the performance of the proposed MVI preoperative prediction models. Finally, the engineering of the constructed prediction model should also be considered for clinical practice.

5 Conclusion

To address the problem of combining non-imaging information from enhanced MRI with the highly subjective determination of the perioperative region, we propose an enhanced preoperative prediction model for MVI in HCC using an optimized MR Radiomics combination strategy and a machine learning predictor. First, the HCC was manually segmented from 125 × 512 × 512 × N abdominal enhanced T1WI images during the arterial phase, generating 125 × 512 × 512 × N HCC mask images. Second, 125 × 1,692 MR Radiomics features of HCC are extracted from abdominal enhanced T1WI images based on the HCC mask images. Third, the 125 × N selected and 125 × 10 fused MR Radiomics features are determined using the proposed optimized MR Radiomics combination strategy with 5-fold cross-validation. Finally, the best preoperative prediction model is constructed using a random forest (RF) predictor with 125 × N selected and 125 × 10 fused MR Radiomics features. The results show that the proposed MVI preoperative prediction model (RF + LASSO + SPECTRAL-10) achieves a mean accuracy of 0.7520 ± 0.0867, a mean precision of 0.7354 ± 0.1863, a mean recall of 0.6955 ± 0.2203, a mean F1-score of 0.6943 ± 0.1437, and a mean AUC of 0.7962 ± 0.1700. The proposed best preoperative prediction model can effectively predict MVI in HCC, potentially serving as a strong decision-making tool for these vulnerable populations.

Data availability statement

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

Author contributions

MF: Formal analysis, Writing – original draft, Software, Visualization, Methodology. YY: Writing – review & editing, Software, Methodology, Writing – original draft, Funding acquisition, Formal analysis, Validation. ZD: Writing – original draft, Investigation, Validation, Formal analysis. ZC: Formal analysis, Methodology, Software, Writing – original draft, Visualization. LL: Software, Writing – original draft, Formal analysis, Visualization, Methodology. ZW: Visualization, Formal analysis, Methodology, Software, Writing – original draft. XL: Validation, Formal analysis, Writing – original draft, Investigation. TG: Writing – original draft, Formal analysis, Validation, Investigation. YM: Validation, Writing – original draft, Investigation, Formal analysis. QL: Validation, Formal analysis, Writing – review & editing, Investigation, Writing – original draft. ZZ: Writing – original draft, Investigation, Validation. TL: Investigation, Writing – original draft. JZ: Project administration, Supervision, Writing – original draft, Conceptualization, Writing – review & editing, Resources, Data curation. YK: Funding acquisition, Writing – original draft, Formal analysis, Writing – review & editing, Visualization, Project administration, Conceptualization.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research was supported by the National Key Research and Development Program of China, Grant Nos. 2022YFF0710800 and 2022YFF0710802; the Heilongjiang Province Natural Science Foundation Joint Guidance Project of China, Grant No. LH2024F044; the Zhongnanshan Medical Foundation of Guangdong Province of China, Grant No. ZNSXS-20230001; the National Natural Science Foundation of China, Grant Nos. 62071311 and 82572349.

Conflict of interest

YY was employed by Shenzhen Lanmage Medical Technology Co., Ltd.

The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2026.1764733/full#supplementary-material

References

1. Gao, Q, Zhu, H, Dong, L, Shi, W, Chen, R, Song, Z, et al. Integrated proteogenomic characterization of HBV-related hepatocellular carcinoma. Cell. (2019) 179:561–77. doi: 10.1016/j.cell.2019.08.052,

PubMed Abstract | Crossref Full Text | Google Scholar

2. Yong, H, Shengxi, H, and Xiufeng, L. Advances and prospects of systemic therapy for hepatocellular carcinoma. J Clin Hepatol. (2025) 41:1491–6. doi: 10.12449/JCH250803

Crossref Full Text | Google Scholar

3. Chan, SL, Sun, HC, Xu, Y, Zeng, H, El-Serag, HB, Lee, JM, et al. The lancet commission on addressing the global hepatocellular carcinoma burden: comprehensive strategies from prevention to treatment. Lancet. (2025) 406:731–78. doi: 10.1016/s0140-6736(25)01042-6,

PubMed Abstract | Crossref Full Text | Google Scholar

4. Yu, Y, Wang, XH, Hu, WJ, Chen, D-H, Hu, Z-L, and Li, S-Q. Patterns, risk factors, and outcomes of recurrence after hepatectomy for hepatocellular carcinoma with and without microvascular invasion. J Hepatocell Carcinoma. (2024) 11:801–12. doi: 10.2147/jhc.s438850,

PubMed Abstract | Crossref Full Text | Google Scholar

5. Zheng, J, Chakraborty, J, Chapman, WC, Gerst, S, Gonen, M, Pak, LM, et al. Preoperative prediction of microvascular invasion in hepatocellular carcinoma using quantitative image analysis. J Am Coll Surg. (2017) 225:778–88. doi: 10.1016/j.jamcollsurg.2017.09.003,

PubMed Abstract | Crossref Full Text | Google Scholar

6. Yang, L, Gu, D, Wei, J, Yang, C, Rao, S, Wang, W, et al. A radiomics nomogram for preoperative prediction of microvascular invasion in hepatocellular carcinoma. Liver Cancer. (2019) 8:373–86. doi: 10.1159/000494099,

PubMed Abstract | Crossref Full Text | Google Scholar

7. Zhang, Y, Chen, J, Yang, C, Dai, Y, and Zeng, M. Preoperative prediction of microvascular invasion in hepatocellular carcinoma using diffusion-weighted imaging-based habitat imaging. Eur Radiol. (2024) 34:3215–25. doi: 10.1007/s00330-023-10339-2,

PubMed Abstract | Crossref Full Text | Google Scholar

8. Wang, F, Liao, HZ, Chen, XL, Lei, H, Luo, G-H, Chen, G-D, et al. Preoperative prediction of microvascular invasion: new insights into personalized therapy for early-stage hepatocellular carcinoma. Quant Imaging Med Surg. (2024) 14:5205. doi: 10.21037/qims-24-44,

PubMed Abstract | Crossref Full Text | Google Scholar

9. Liu, J, Zhuang, G, Bai, S, Hu, Z, Xia, Y, Lu, C, et al. The comparison of surgical margins and type of hepatic resection for hepatocellular carcinoma with microvascular invasion. Oncologist. (2023) 28:e1043–51. doi: 10.1093/oncolo/oyad124,

PubMed Abstract | Crossref Full Text | Google Scholar

10. Yang, S, Ni, H, Zhang, A, Zhang, J, Zang, H, and Ming, Z. Significance of anatomical resection and wide surgical margin for HCC patients with MVI undergoing laparoscopic hepatectomy: a multicenter study. Eur J Surg Oncol. (2025) 51:109353. doi: 10.1016/j.ejso.2024.109353,

PubMed Abstract | Crossref Full Text | Google Scholar

11. Yang, S, Ni, H, Zhang, A, Zhang, J, Liang, H, Li, X, et al. Grading severity of MVI impacts long-term outcomes after laparoscopic liver resection for early-stage hepatocellular carcinoma: a multicenter study. Am J Surg. (2024) 238:115988. doi: 10.1016/j.amjsurg.2024.115988,

PubMed Abstract | Crossref Full Text | Google Scholar

12. Hu, H, Qi, S, Zeng, S, Zhang, P, He, L, Wen, S, et al. Importance of microvascular invasion risk and tumor size on recurrence and survival of hepatocellular carcinoma after anatomical resection and non-anatomical resection. Front Oncol. (2021) 11:621622. doi: 10.3389/fonc.2021.621622,

PubMed Abstract | Crossref Full Text | Google Scholar

13. Li, J, Yang, F, Li, J, Huang, ZY, Cheng, Q, and Zhang, EL. Postoperative adjuvant therapy for hepatocellular carcinoma with microvascular invasion. World J Gastrointest Surg. (2023) 15:19. doi: 10.4240/wjgs.v15.i1.19,

PubMed Abstract | Crossref Full Text | Google Scholar

14. Candita, G, Rossi, S, Cwiklinska, K, Fanni, SC, Cioni, D, Lencioni, R, et al. Imaging diagnosis of hepatocellular carcinoma: a state-of-the-art review. Diagnostics. (2023) 13:625. doi: 10.3390/diagnostics13040625,

PubMed Abstract | Crossref Full Text | Google Scholar

15. Tan, B, Shen, F, Yang, C, Zeng, M, Wang, J, Zuo, C, et al. Imaging techniques for liver cancer diagnosis In: Diagnosis and treatment of liver cancer in China: past, present and prospects. Singapore: Springer (2025). 49–69.

Google Scholar

16. Albakri, AA, Alzahrani, MM, and Alghamdi, SH. Medical imaging in pregnancy: safety, appropriate utilization, and alternative modalities for imaging pregnant patients. Cureus. (2024) 16:e54346. doi: 10.7759/cureus.54346,

PubMed Abstract | Crossref Full Text | Google Scholar

17. LeJeune, CL, Nougaret, S, Massera, RT, Jha, P, Bonanno, C, Bosmans, H, et al. Abdominal shielding not recommended for diagnostic imaging with ionising radiation during pregnancy. Lancet Oncol. (2025) 26:1413–6. doi: 10.1016/s1470-2045(25)00368-7,

PubMed Abstract | Crossref Full Text | Google Scholar

18. Vernuccio, F, Cannella, R, Bartolotta, TV, Galia, M, Tang, A, and Brancatelli, G. Advances in liver US, CT, and MRI: moving toward the future. Eur Radiol Exp. (2021) 5:52. doi: 10.1186/s41747-021-00250-0,

PubMed Abstract | Crossref Full Text | Google Scholar

19. Marrero, JA, Kulik, LM, Sirlin, CB, Zhu, AX, Finn, RS, Abecassis, MM, et al. Diagnosis, staging, and management of hepatocellular carcinoma: 2018 practice guidelines by the American Association for the Study of Liver Diseases. Hepatology. (2018) 68:723–50. doi: 10.1002/hep.29913,

PubMed Abstract | Crossref Full Text | Google Scholar

20. Vogel, A, Cervantes, A, Chau, I, Daniele, B, Llovet, JM, Meyer, T, et al. Hepatocellular carcinoma: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. (2018) 29:iv238–55. doi: 10.1093/annonc/mdy308,

PubMed Abstract | Crossref Full Text | Google Scholar

21. Omata, M, Cheng, AL, Kokudo, N, Kudo, M, Lee, JM, Jia, J, et al. Asia–Pacific clinical practice guidelines on the management of hepatocellular carcinoma: a 2017 update. Hepatol Int. (2017) 11:317–70. doi: 10.1007/s12072-017-9799-9,

PubMed Abstract | Crossref Full Text | Google Scholar

22. Tang, L, Shao, H, Yang, J, Xu, J, Li, J, Feng, Y, et al. (2024). Multi-modal learning for predicting the progression of transarterial chemoembolization therapy in hepatocellular carcinoma. Chinese Conference on Pattern Recognition and Computer Vision (PRCV). 178–193

Google Scholar

23. Sun, Z, Shi, Z, Xin, Y, Zhao, S, Jiang, H, Wang, D, et al. Artificial intelligent multi-modal point-of-care system for predicting response of transarterial chemoembolization in hepatocellular carcinoma. Front Bioeng Biotechnol. (2021) 9:761548. doi: 10.3389/fbioe.2021.761548,

PubMed Abstract | Crossref Full Text | Google Scholar

24. Yang, Y, Li, W, Guo, Y, Zeng, N, Wang, S, Chen, Z, et al. Lung radiomics features for characterizing and classifying COPD stage based on feature combination strategy and multi-layer perceptron classifier. Math Biosci Eng. (2022) 19:7826–55. doi: 10.3934/mbe.2022366,

PubMed Abstract | Crossref Full Text | Google Scholar

25. Yang, Y, and Guo, Y. Ischemic stroke outcome prediction with diversity features from whole brain tissue using deep learning network. Front Neurol. (2024) 15:1394879. doi: 10.3389/fneur.2024.1394879,

PubMed Abstract | Crossref Full Text | Google Scholar

26. Guo, Y, Yang, Y, Cao, F, Wang, M, Luo, Y, Guo, J, et al. A focus on the role of DSC-PWI dynamic radiomics features in diagnosis and outcome prediction of ischemic stroke. J Clin Med. (2022) 11:5364. doi: 10.3390/jcm11185364,

PubMed Abstract | Crossref Full Text | Google Scholar

27. Guo, Y, Yang, Y, Cao, F, Li, W, Wang, M, Luo, Y, et al. Novel survival features generated by clinical text information and radiomics features may improve the prediction of ischemic stroke outcome. Diagnostics. (2022) 12:1664. doi: 10.3390/diagnostics12071664,

PubMed Abstract | Crossref Full Text | Google Scholar

28. Guo, Y, Yang, Y, Wang, M, Luo, Y, Guo, J, Cao, F, et al. The combination of whole-brain features and local-lesion features in DSC-PWI may improve ischemic stroke outcome prediction. Life. (2022) 12:1847. doi: 10.3390/life12111847,

PubMed Abstract | Crossref Full Text | Google Scholar

29. Guo, Y, Yang, Y, Cao, F, Liu, Y, Li, W, Yang, C, et al. Radiomics features of DSC-PWI in time dimension may provide a new chance to identify ischemic stroke. Front Neurol. (2022) 13:889090. doi: 10.3389/fneur.2022.889090,

PubMed Abstract | Crossref Full Text | Google Scholar

30. Yang, Y, Wang, S, Zeng, N, Duan, W, Chen, Z, Liu, Y, et al. Lung radiomics features selection for COPD stage classification based on auto-metric graph neural network. Diagnostics. (2022) 12:2274. doi: 10.3390/diagnostics12102274,

PubMed Abstract | Crossref Full Text | Google Scholar

31. Deng, X, Li, W, Yang, Y, Wang, S, Zeng, N, Xu, J, et al. COPD stage detection: leveraging the auto-metric graph neural network with inspiratory and expiratory chest CT images. Med Biol Eng Comput. (2024) 62:1733–49. doi: 10.1007/s11517-024-03016-z,

PubMed Abstract | Crossref Full Text | Google Scholar

32. Yang, Y, Li, W, Guo, Y, Liu, Y, Li, Q, Yang, K, et al. Early COPD risk decision for adults aged from 40 to 79 years based on lung radiomics features. Front Med. (2022) 9:845286. doi: 10.3389/fmed.2022.845286,

PubMed Abstract | Crossref Full Text | Google Scholar

33. Wang, S, Li, W, Zeng, N, Xu, J, Yang, Y, Deng, X, et al. Acute exacerbation prediction of COPD based on auto-metric graph neural network with inspiratory and expiratory chest CT images. Heliyon. (2024) 10:e28724. doi: 10.1016/j.heliyon.2024.e28724,

PubMed Abstract | Crossref Full Text | Google Scholar

34. Yang, Y, Chen, Z, Li, W, Zeng, N, Guo, Y, Wang, S, et al. Multi-modal data combination strategy based on chest HRCT images and PFT parameters for intelligent dyspnea identification in COPD. Front Med. (2022) 9:980950. doi: 10.3389/fmed.2022.980950,

PubMed Abstract | Crossref Full Text | Google Scholar

35. Yang, H, Guo, Y, Lu, J, Hassan, H, Cao, A, Yang, Y, et al. Optimizing acute ischemic stroke outcome prediction by integrating radiomics features of DSC-PWI and perfusion parameter maps. Front Neurol. (2025) 16:1528812. doi: 10.3389/fneur.2025.1528812,

PubMed Abstract | Crossref Full Text | Google Scholar

36. Yang, Y, Li, W, Kang, Y, Guo, Y, Yang, K, Li, Q, et al. A novel lung radiomics feature for characterizing resting heart rate and COPD stage evolution based on radiomics feature combination strategy. Math Biosci Eng. (2022) 19:4145–65. doi: 10.3934/mbe.2022191,

PubMed Abstract | Crossref Full Text | Google Scholar

37. Li, Z, Xu, L, Zhu, S, Qi, X, Zhang, W, and Tang, Y. Current advances in classification, prediction and management of microvascular invasion in hepatocellular carcinoma. J Cell Mol Med. (2025) 29:e70746. doi: 10.1111/jcmm.70746,

PubMed Abstract | Crossref Full Text | Google Scholar

38. Wang, F, Yan, CY, Qin, Y, Wang, ZM, Liu, D, He, Y, et al. Multiple machine-learning fusion model based on Gd-EOB-DTPA-enhanced MRI and aminotransferase-to-platelet ratio and gamma-glutamyl transferase-to-platelet ratio to predict microvascular invasion in solitary hepatocellular carcinoma: a multicenter study. J Hepatocell Carcinoma. (2024) 11:427–42. doi: 10.2147/JHC.S449737,

PubMed Abstract | Crossref Full Text | Google Scholar

39. Feng, ST, Jia, Y, Liao, B, Huang, B, Zhou, Q, Li, X, et al. Preoperative prediction of microvascular invasion in hepatocellular cancer: a radiomics model using Gd-EOB-DTPA-enhanced MRI. Eur Radiol. (2019) 29:4648–59. doi: 10.1007/s00330-018-5935-8,

PubMed Abstract | Crossref Full Text | Google Scholar

40. Rubin-Delanchy, P, Cape, J, Tang, M, and Priebe, CE. A statistical interpretation of spectral embedding: the generalised random dot product graph. J R Stat Soc B. (2022) 84:1446–73. doi: 10.1111/rssb.12509

Crossref Full Text | Google Scholar

41. Gallagher, I, Jones, A, Bertiger, A, Priebe, CE, and Rubin-Delanchy, P. Spectral embedding of weighted graphs. J Am Stat Assoc. (2024) 119:1923–32. doi: 10.1080/01621459.2023.2225239

Crossref Full Text | Google Scholar

42. Tharwat, A. Independent component analysis: an introduction. Appl Comput Inform. (2021) 17:222–49. doi: 10.1016/j.aci.2018.08.006

Crossref Full Text | Google Scholar

43. Wang, Y, Zhang, Z, and Lin, Y. Multi-cluster feature selection based on isometric mapping. IEEE CAA J Autom Sin. (2021) 9:570–2. doi: 10.1109/JAS.2021.1004398

Crossref Full Text | Google Scholar

44. Healy, J, and McInnes, L. Uniform manifold approximation and projection. Nat Rev Methods Primers. (2024) 4:82. doi: 10.1038/s43586-024-00363-x

Crossref Full Text | Google Scholar

45. Demler, OV, Pencina, MJ, and D’Agostino, RB Sr. Misuse of DeLong test to compare AUCs for nested models. Stat Med. (2012) 31:2577–87. doi: 10.1002/sim.5328,

PubMed Abstract | Crossref Full Text | Google Scholar

46. Karabayir, I, and Akbilgic, O. Generalizability of electrocardiographic artificial intelligence. npj Cardiovasc Health. (2025) 2:38. doi: 10.1038/s44325-025-00078-2

Crossref Full Text | Google Scholar

47. Langsetmo, L, Schousboe, JT, Taylor, BC, Cauley, JA, Fink, HA, Cawthon, PM, et al. Advantages and disadvantages of random forest models for prediction of hip fracture risk versus mortality risk in the oldest old. JBMR Plus. (2023) 7:e10757. doi: 10.1002/jbm4.10757,

PubMed Abstract | Crossref Full Text | Google Scholar

48. Elsken, T, Metzen, JH, and Hutter, F. Neural architecture search: a survey. J Mach Learn Res. (2019) 20:1–21. doi: 10.1007/978-3-030-05318-5_11

Crossref Full Text | Google Scholar

49. Yang, G, Chen, Z, Guo, P, Ma, J, Huang, J, Jin, C, et al. BLoss-DDNet: bending loss and dual-task decoding network for overlap-ping cell nucleus segmentation of cervical clinical LBC images. Front Artif Intell. (2025) 8:1649452. doi: 10.3389/frai.2025.1649452,

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: enhanced T1-weighted magnetic resonance imaging, hepatocellular carcinoma, machine learning algorithms, microvascular invasion, preoperative prediction, radiomics

Citation: Feng M, Yang Y, Dai Z, Chen Z, Li L, Wu Z, Li X, Guo T, Meng Y, Li Q, Zhao Z, Li T, Zhang J and Kang Y (2026) Enhanced preoperative prediction for microvascular invasion in hepatocellular carcinoma through an optimized MR Radiomics combination strategy and machine learning predictor. Front. Med. 13:1764733. doi: 10.3389/fmed.2026.1764733

Received: 10 December 2025; Revised: 15 January 2026; Accepted: 27 January 2026;
Published: 11 February 2026.

Edited by:

Cheng Wei, University of Dundee, United Kingdom

Reviewed by:

Na Zhao, Southeast University, China
Runnan He, Tianjin University, China

Copyright © 2026 Feng, Yang, Dai, Chen, Li, Wu, Li, Guo, Meng, Li, Zhao, Li, Zhang and Kang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jialin Zhang, amx6MjAwMEB5ZWFoLm5ldA==; Yan Kang, a2FuZ3lhbkBzenR1LmVkdS5jbg==

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.