- 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 (1–3). 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 (4–7). 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 (19–21). 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 (22–39). However, deep learning-based prediction tasks are suitable for scenarios with large amount of data, complex problems, and high performance requirements (22–24). 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 (37–39).
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 (24–36) and has also opened new avenues for non-invasive preoperative prediction of MVI in HCC based on enhanced MRI (37–39). 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.
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. 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 (24–36), 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) (24–36, 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 (26–36, 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 (26–36, 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 (26–36). 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. 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.
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) (24–36). 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:
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.
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 4a–e 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. 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. 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 6a–d 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. 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%.
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,
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
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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.
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,
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,
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,
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,
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,
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,
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
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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,
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
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
42. Tharwat, A. Independent component analysis: an introduction. Appl Comput Inform. (2021) 17:222–49. doi: 10.1016/j.aci.2018.08.006
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
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
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,
46. Karabayir, I, and Akbilgic, O. Generalizability of electrocardiographic artificial intelligence. npj Cardiovasc Health. (2025) 2:38. doi: 10.1038/s44325-025-00078-2
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,
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
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
Edited by:
Cheng Wei, University of Dundee, United KingdomCopyright © 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==
Mengting Feng1