A Multi-Parametric Radiomics Nomogram for Preoperative Prediction of Microvascular Invasion Status in Intrahepatic Cholangiocarcinoma

Background: Intrahepatic cholangiocarcinoma (ICC) is the second most common primary liver cancer with increasing incidence in the last decades. Microvascular invasion (MVI) is a poor prognostic factor for patients with ICC, which correlates early recurrence and poor prognosis, and it can affect the selection of personalized therapeutic regime.

Purpose: This study aimed to develop and validate a radiomics-based nomogram for predicting MVI in ICC patients preoperatively.

Methods: A total of 163 pathologically confirmed ICC patients (training cohort: n = 130; validation cohort: n = 33) with postoperative Ga-DTPA-enhanced MR examination were enrolled, and a time-independent test cohort (n = 24) was collected for external validation. Univariate and multivariate analyses were used to determine the independent predictors of MVI status, which were then incorporated into the MVI prediction nomogram. Least absolute shrinkage and selection operator logistic regression was performed to select optimal features and construct radiomics models. The prediction performances of models were assessed by receiver operating characteristic (ROC) curve analysis. The performance of the MVI prediction nomogram was evaluated by its calibration, discrimination, and clinical utility.

Results: Larger tumor size (p = 0.003) and intrahepatic duct dilatation (p = 0.002) are independent predictors of MVI. The final radiomics model shows desirable and stable prediction performance in the training cohort (AUC = 0.950), validation cohort (AUC = 0.883), and test cohort (AUC = 0.812). The MVI prediction nomogram incorporates tumor size, intrahepatic duct dilatation, and the final radiomics model and achieves excellent predictive efficacy in training cohort (AUC = 0.953), validation cohort (AUC = 0.861), and test cohort (AUC = 0.819), fitting well in calibration curves (p > 0.05). Decision curve and clinical impact curve further confirm the clinical usefulness of the nomogram.

Conclusion: The nomogram incorporating tumor size, intrahepatic duct dilatation, and the final radiomics model is a potential biomarker for preoperative prediction of the MVI status in ICC patients.

1 Introduction

Intrahepatic cholangiocarcinoma (ICC) is the second most common primary liver cancer and accounts for 10%–15% of all cases, which arises from cholangiocytes of intrahepatic bile ducts or bile ductules (1–3). ICC has three predominant macroscopic growth patterns: mass-forming type, periductal infiltrating type, and intraductal papillary type (4). Several studies have reported increasing incidence of ICC in the last decades (5, 6) and the 5-year survival rate is still lower than 10% (7). At present, hepatectomy is still the most effective treatment for long-term survival of ICC patients (8, 9), and several poor prognostic factors have been reported, including lymph node metastasis, microvascular invasion (MVI), tumor size ≥5 cm, and multiple nodules (10).

MVI is an important histopathological feature and refers to the cancer cell nest in vessels of the surrounding hepatic tissue lined with endothelial cells (11). As a poor prognostic factor, MVI correlates early recurrence and poor outcomes and is an independent factor for overall survival in ICC patients (12, 13). However, the status of MVI is difficult to detect by radiographic images and can only be determined by histological evaluation after hepatectomy (14). In addition, MVI can affect the selection of personalized therapeutic regime, for instance, ICC patients without MVI do not need to receive adjuvant chemotherapy after R0 resections (10). Herein, the preoperative determination of MVI status is of great value in ICC patients, and it holds promises for effective patient management and estimation of outcomes.

Radiomics is a powerful and sophisticated image mining tool, and it can improve diagnostic accuracy and predict prognosis by high-throughput selecting imaging features from medical images (15). Also, several studies have constructed radiomics-based nomogram in distinguishing different pathological types of primary liver cancer (16) and predicting MVI of hepatocellular carcinoma preoperatively (17, 18). Recently, radiomics nomograms have been established for the prediction of lymph node metastasis (19), early recurrence (20), and prognosis after hepatectomy (21) in ICC patients. For MVI prediction of ICC, Zhou et al. showed the promise of seven wavelet features extracted from preoperative dynamic contrast-enhanced (DCE) MR images with an area under curve (AUC) of 0.873 (22).

In the present study, we focused on mass-forming ICCs and aimed to develop and validate a radiomics nomogram integrating clinical, imaging, and radiomics features for preoperative prediction of MVI in ICC. In order to verify the accuracy, the radiomics nomogram will be further validated by a test cohort.

2 Materials and Methods

2.1 Patients

Zhongshan Hospital, Fudan University, and Xuzhou Central Hospital ethics committees approved this retrospective study, and patient informed consent was waived. Between June 2015 and June 2019, 163 pathologically confirmed ICC patients (118 men and 45 women; 60.48 ± 11.42 years) after hepatectomy with postoperative Ga-DTPA-enhanced MRI examination from Zhongshan Hospital were enrolled by the following inclusion criteria (Figure 1): (a) without previous history of liver cancer treatment (including hepatectomy, transcatheter arterial chemoembolization, radiofrequency ablation, chemotherapy, radiotherapy, and immunosuppressive therapy); (b) single mass-forming type ICC with the longest diameter ≥1.0 cm, and without macrovascular invasion and lymphatic metastasis; (c) complete histopathologic description of ICC; (d) MRI scanned within 30 days before surgery; (e) sufficient MR image quality satisfied the diagnostic criteria. All enrolled patients from June 2015 to June 2019 were divided into training cohort (n = 130, 38 MVI positive and 92 MVI negative) and validation cohort (n = 33, 10 MVI positive and 23 MVI negative) in a ratio of 8:2. Importantly, between July 2019 and October 2021, a time-independent test cohort (n = 24, 12 MVI positive and 12 MVI negative) from Zhongshan Hospital and Xuzhou Central Hospital was collected for external validation (Supplementary Table S1).

FIGURE 1

Figure 1 Study flowchart of the enrolled patients.

2.2 Laboratory Tests and Histopathology

Demographic and preoperative laboratory indexes (Table 1) including serum alpha-fetoprotein (AFP), carcinoembryonic antigen (CEA), carbohydrate antigen 19-9 (CA199), hepatitis B virus (HBV), and HBV-DNA loads were collected from our electronic medical records system. ICC specimens were sampled using 7-point baseline sampling protocol (11). Pathological characteristics including tumor number, Edmondson-Steiner grade, and MVI status were assessed in consensus by two experienced abdominal pathologists. MVI was defined as the presence of tumor cell nest (the number of suspended tumor cell more than 50) in the portal vein, hepatic vein, or large capsular vessel of the surrounding hepatic tissue that was only visible by microscopy (11, 23).

TABLE 1

Table 1 Comparison of MVI status and clinicoradiologic characteristics in ICC patients of training and validation cohorts.

2.3 Gd-DTPA MR Imaging

All patients underwent MR imaging with intravenous injection of 0.2 mmol/kg gadopentetate dimeglumine (Gd-DTPA; Bayer HealthCare, Berlin, Germany) and immediately followed by a 20-ml saline flush at a speed of 2 ml/s. Taking Magnetom Aera 1.5T scanner (Siemens Healthcare, Erlangen, Germany) as an example, imaging sequences included axial T2-weighted imaging with fat suppression (T2WI-FS), diffusion-weighted imaging (DWI), in-phase and opposed-phase T1-weighted imaging (IP-OP T1WI), axial precontrast three-dimensional volumetric-interpolated breath-hold examination T1-weighted imaging (3D-VIBE T1WI) with fat suppression, and postcontrast dynamic-enhanced 3D-VIBE-T1WI at arterial phase (AP, 20–30 s), portal venous phase (PVP, 60–70 s), and delayed phase (DP, 180 s). Detailed parameters of each acquisition sequence are shown in Supplementary Table S2.

2.4 MR Images Analysis

All MR images were reviewed independently on a picture archiving and communication system (PACS; Pathspeed, GE Medical Systems Integrated Imaging Solutions, Chicago, IL, USA) by 2 experienced abdominal radiologists (XM and XL with 10 and 15 years of experience, respectively). Both radiologists were blinded to all demographic, clinical, laboratory, and histopathologic information. If any discrepancies occurred, a consensus was reached after discussion. The following imaging features were assessed by 2 abdominal radiologists: (a) tumor size, defined as the maximum diameter on transverse T1WI image; (b) tumor morphology, including spherical/hemispherical/oval, lobulated and irregular; (c) signal intensity on T1WI, T2WI-FS, and DWI, including hypointense, isointense, and hyperintense; (d) target sign on T2WI-FS and DWI, defined as peripheral hyperintense with central isointense/hypointense (24); (e) rim enhancement on AP, defined as peripheral enhancement of the lesion on AP, including complete and incomplete rim; (f) enhancement pattern, including gradual and filling, arterial and persistent, and wash-in and wash-out enhancement; (g) the liver imaging reporting and data system (LI-RADS), defined as LR category based on LI-RADSv2018 (25); (h) intrahepatic duct dilatation, defined as intrahepatic duct dilatation within or outside of the lesion; (i) hepatic capsular retraction, defined as retraction of hepatic capsular adjacent to the lesion; (j) visible vessel penetration, defined as the presence of penetrating vessels in the lesion, including hepatic artery, portal vein, and hepatic vein (26); and (k) peripherally hepatic enhancement, defined as enhancement of liver parenchyma around the lesion on any phase.

2.5 Radiomics Analysis

2.5.1 Workflow

The workflow of the radiomics analysis included tumor segmentation, feature extraction, feature selection, model construction, model analysis, and evaluation (Figure 2).

FIGURE 2

Figure 2 Study flowchart of the radiomics analysis.

2.5.2 Image Segmentation

The whole tumor segmentation was manually delineated in ITK-SNAP (http://www.itksnap.org/pmwiki/pmwiki.php) by an abdominal radiologist with 6 years of experience (XQ) and validated by a senior abdominal radiologist with 15 years of experience (XL). Volumes of interests (VOIs) were drawn on 6 sequences: DWI with b-values of 500 s/mm², T2WI-FS, 3D-VIBE T1WI, AP, PVP, and DP.

2.5.3 Feature Extraction

Radiomics features were extracted from the VOIs using uAI Research Portal (Version: 20210730), and 2,600 radiomics features were extracted from each sequence (Supplementary Table S3). These radiomics features were classified into First Order, Shape, Gray-Level Co-occurrence Matrix (GLCM), Gray-Level Size Zone Matrix (GLSZM), Gray-Level Run Length Matrix (GLRLM), Neighboring Gray Tone Difference Matrix (NGTDM), and Gray-Level Dependence Matrix (GLDM) features.

2.5.4 Feature Selection

To eliminate index dimension difference, the extracted radiomics features of each sequence were standardized into a normal distribution with z-scores. A test-retest procedure was performed on 30 randomly selected tumors, reproducible radiomics features were considered features with intraclass correlation coefficient greater than 0.75 and included in the following feature selection procedures. The variance threshold, SelectKBest and least absolute shrinkage and selection operator (LASSO) methods were performed to select optimal features of each sequence (Supplementary Table S3 and Supplementary Figure S1). The corresponding radiomics score (Rad-score) of each sequence in the training and validation cohorts was calculated (Supplementary Table S4 and Supplementary Figure S2).

2.5.5 Model Construction

To construct clinical and imaging models, the univariate analysis was used to assess the potential predictors of MVI status, and the multivariate logistic regression analysis was used to determine the independent predictors of MVI status. Radiomics models of each sequence were constructed by the corresponding optimal features. The sequences with Rad-scores showed significant differences between MVI-positive and MVI-negative ICCs in both the training cohort and validation cohort were selected for the final radiomics model construction. The MVI prediction model incorporated imaging model and the final radiomics model. All models were constructed with logistic regression (LR), random forest (RF), and support vector machine (SVM) classifiers respectively for comparison.

2.5.6 Model Analysis and Evaluation

The receiver operating characteristic curves were plotted, and the AUC, accuracy, sensitivity, specificity, F1-score, and precision were calculated to quantify the predictive efficacy of each model in training, validation, and test cohorts. The comparison of predictive performances between multiple models was performed by the Delong test. A radiomics nomogram was built on the MVI prediction model with the LR classifier. Hosmer-Lemeshow test was performed to identify the agreement between nomogram-predicted MVI status and actual MVI status, and calibration curves in the training and validation cohorts were plotted. Decision curve and clinical impact curve were plotted for assessing the clinical usefulness of the nomogram by quantifying the net benefits at different risk thresholds.

2.6 Statistical Analysis

Clinical and imaging features were analyzed for statistical differences in the training, validation, and test cohorts by Student’s t-test, Mann-Whitney U test, Wilcoxon test, Chi-square test, or Fisher’s exact test, as appropriate. The statistical analyses were conducted using the IBM SPSS Statistics (version 20) and R software (version 3.4.1). A two-tailed p-value of less than 0.05 was considered statistically significant.

3 Results

3.1 Clinicoradiologic Characteristics and Performances

Comparison of MVI status and clinicoradiologic characteristics in training and validated ICC patients are shown in Table 1. The final cohort of 163 patients with single ICC was divided into training cohort (n = 130, 92 patients were MVI negative and 38 patients were MVI positive) and validation cohort (n = 33, 23 patients were MVI negative and 10 patients were MVI positive). There is no significant difference of the status of MVI between training and validation cohort (p = 0.904). Univariate analysis of clinicoradiologic characteristics indicates that serum CEA level (p = 0.035; OR = 2.716, 95% CI: 1.065–6.918), serum CA199 level (p = 0.030; OR = 2.346, 95% CI: 1.092–5.139), Edmondson-Steiner grade (p = 0.021; OR = 3.126, 95% CI: 1.254–8.977), tumor size (p < 0.001; OR = 1.033, 95% CI: 1.016–1.052), tumor morphology (p= 0.071; OR = 1.604, 95% CI: 0.964–2.708), and intrahepatic duct dilatation (p < 0.001; OR =5.610, 95% CI: 2.505–13.308) are significantly associated with MVI. At the multivariate analysis, tumor size (p = 0.003; OR = 1.032, 95% CI: 1.011–1.055) and intrahepatic duct dilatation (p = 0.002; OR = 4.552, 95% CI: 1.777–12.259) are independent predictors of MVI (Table 2). The imaging model constructed with two predictors has an AUC of 0.726 in the training cohort and 0.522 in the validation cohort (Table 3). Examples of representative radiological characteristics of MVI-positive and MVI-negative ICCs are shown in Figure 3.

TABLE 2

Table 2 Univariate and multivariate analyses of predictive characteristics related with MVI status in ICC.

TABLE 3

Table 3 The performance of imaging, radiomics of single and multiple MR sequences, and final fusion models for predicting MVI in ICC patients.

FIGURE 3

Figure 3 Two examples of representative MVI-negative and MVI-positive ICCs. (A–D) A 62-year-old man with a well-circumscribed MVI-negative ICC in hepatic segment II (arrows). DWI image showed target sign (A), axial arterial phase image showed rim enhancement (B), and portal vein phase image (C) and delayed phase image (D) showed the typical enhancement type of ICC: gradual and filling enhancement. (E–H) A 62-year-old man with a lobulated MVI-positive ICC in hepatic segment IV (arrows). DWI image showed hyperintensity (E), axial arterial phase image showed marginal moderate enhancement with no internal enhancement (F) and dilated bile ducts next to tumor (arrowheads), and portal vein phase image (G) and delayed phase image (H) showed the typical enhancement type of ICC: gradual and filling enhancement.

3.2 Performance of Radiomics Features From Single MR Sequence

Robust radiomics features of each single MR sequence were selected by reproducibility analysis and LASSO analysis (Supplementary Table S3 and Supplementary Figure S1). The predictive performance of radiomics features on each MR sequence is shown in Table 3. However, all single sequences show overfit predictive performance in the training cohort (AUCs = 1.000) and poor predictive performance in the validation cohort (AUCs: 0.422–0.665). The Rad-score of each MR sequence in the training and validation cohorts are shown in Supplementary Table S4, and Rad-scores of the diffusion-weighted imaging, precontrast T1-weighted imaging, and delayed phase imaging show significant differences between MVI-positive and MVI-negative ICCs in both the training cohort (p < 0.001) and validation cohort (p_DWI = 0.025, p_T1 = 0.003, p_T1D = 0.001) (Supplementary Table S4 and Supplementary Figure S2). Therefore, further analysis about diffusion-weighted imaging, precontrast T1-weighted imaging, and delayed phase imaging was conducted. The detailed information of 22 radiomics features in diffusion-weighted imaging, 17 radiomics features in precontrast T1-weighted imaging, and 5 radiomics features in delayed phase imaging are shown in Supplementary Table S5.

3.3 Performance of Radiomics Features From Multiple MR Sequences

Three pairwise combination models of the diffusion-weighted images, precontrast T1-weighted images, and delayed phase images are constructed, and all show a satisfying performance in both the training cohort (AUC = 0.883–0.941) and validation cohort (AUC = 0.817–0.874) (Table 3). The final radiomics model incorporates these three sequences, and it performs better in predicting MVI of ICC in both the training cohort (AUC = 0.950, accuracy = 0.862, sensitivity = 0.921, and specificity = 0.837) and validation cohort (AUC = 0.883, accuracy = 0.788, sensitivity = 0.900, and specificity = 0.739) than three pairwise combination models (Tables 3, 4). Notably, the final radiomics model also performs desirably and stably in the test cohort with AUC, accuracy, sensitivity, and specificity of 0.812 (95% CI: 0.617–1.000), 0.792, 0.750, and 0.833, respectively (Table 3).

TABLE 4

Table 4 The comparison of models in training, validation, and test cohorts.

3.4 Performance of MVI Prediction Model

The MVI prediction model includes imaging model and final radiomics model, and it achieves excellent predictive efficacy in the training cohort (AUC = 0.953, accuracy = 0.892, sensitivity = 0.974, and specificity = 0.859), validation cohort (AUC = 0.861, accuracy = 0.818, sensitivity = 0.900, and specificity = 0.783), and test cohort (AUC = 0.819, accuracy = 0.875, sensitivity = 0.833, and specificity = 0.917) (Table 3).

Moreover, the MVI prediction model performs better than imaging model in the training cohort (AUCs: 0.953 vs. 0.726, p < 0.001) and validation cohort (AUCs: 0.861 vs. 0.522, p = 0.018). However, there is no statistical difference between the MVI prediction model and final radiomics model in the training cohort (AUCs: 0.953 vs. 0.950, p = 0.629), validation cohort (AUCs: 0.861 vs. 0.883, p = 0.202), and test cohort (AUCs: 0.819 vs. 0.812, p = 0.732) (Table 4, Figure 4 and Supplementary Figure S3).

FIGURE 4

Figure 4 Comparison of receiver operating characteristic (ROC) curves for prediction of MVI in ICC. ROC curves of imaging model constructed with tumor size and intrahepatic duct dilatation, radiomics model constructed with diffusion-weighted image, precontrast T1-weighted image, and delayed phase image, and MVI prediction model constructed imaging model and radiomics model in the (A) training, (B) validation, and (C) test cohorts.

3.5 Development and Verification of the Nomogram

The nomogram of the MVI prediction model is presented in Figure 5A, and the formula is as follows. It achieves satisfying performance with AUCs of 0.953 in the training cohort and 0.861 in the validation cohort. Calibration curves (Figures 5B, C) of the nomogram exhibiting satisfactory predictive performances are aligned with the actual MVI estimates in the training (p = 0.364) and validation (p = 0.543) cohorts. Decision curve (Figure 5D) of the nomogram shows that the net benefit is higher than that assuming all patients have MVI. Clinical impact curve (Figure 5E) shows that the predicted probabilities of the nomogram is close to actual high risk with event when risk threshold is 0.2–0.7 and is consistent when risk threshold is over 0.7.

$$Y=-11.420+0.020\times {\text{Imaging}}_{\text{Tumor size}}+0.723\times {\text{Imaging}}_{\text{Intrahepatic duct dilatation}}+9.130\times {\text{Rad Score}}_{\text{DWI}}+10.835\times {\text{Rad Score}}_{\text{T}1}+4.690\times {\text{Rad Score}}_{\text{T}1\text{D}}$$

FIGURE 5

Figure 5 Nomogram of MVI prediction model, calibration curves of the nomogram in the training and validation cohort, decision curve, and clinical impact curve in the overall ICC patients. (A) A nomogram integrates imaging factors including tumor size and intrahepatic duct dilatation, and radiomics factors includes Rad-scores of diffusion-weighted images, precontrast T1-weighted images, and delayed phase images. (B, C) Calibration curves of the nomogram in the training and validation cohort. x-axis is a nomogram-predicted risk of MVI. y-axis is actual risk of MVI, and the diagonal dashed line indicates the ideal prediction by an ideal model. (D) Decision curve for the nomogram in the overall patients. The gray line is the net benefit of assuming that all patients have MVI; the black line is the net benefit of assuming no patients have MVI; and the red line is the expected net benefit of per patient based on the nomogram. (E) Clinical impact curve for the nomogram in the 1,000 simulated samples. The blue dashed line is the actual number of high risk, and the red line is the number of high risk based on nomogram.

Discussion

In this study, we established a radiomics-based nomogram for preoperatively predicting MVI in patients with ICC. The final MVI prediction model achieves a satisfying prediction performance, and it incorporates Ga-DTPA-enhanced MRI-based radiomics features of the diffusion-weighted images, precontrast T1-weighted images, and delayed phase images and imaging features including tumor size and intrahepatic duct dilatation.

High serum level of CEA and CA199 (>100 U/ml) can predict the diagnosis of ICC (27, 28), but the predictive value in predicting MVI status of ICC is still unknown. Although univariate analysis shows elevated serum CEA and CA199 level are significant clinical features, there are no statistical difference between the MVI-positive and MVI-negative groups by multivariate analysis, which are consistent with the studies of Zhou et al. (22) and Ma et al. (29). As other studies on preoperative MVI prediction in HCC (30–32) and ICC (14, 22, 33), tumor size is also an independent predictor of MVI in ICC patients in our study, but intrahepatic duct dilatation is another important imaging predictor, which is different from the studies of Zhou et al. (22). This may be due to the different MVI predictor screening methods; in his study, intrahepatic duct dilatation was excluded because there was no significant difference between MVI-positive and MVI-negative groups in the validation cohort (p = 0.279), but a p-value with 0.097 in the training cohort indicated this predictor should be further analyzed by multivariate analysis. Although hepatic capsular retraction and progressive centripetal enhancement in the venous phase have been described as classical imaging features of ICC (34), there are no statistical difference between MVI-positive and MVI-negative groups. In general, the imaging model constructed with tumor size and intrahepatic duct dilatation yields a good performance in the training cohort, but an unsatisfying performance in the validation cohort limits its application in preoperatively predicting MVI. Therefore, a combined model based on radiomics is necessary.

Of 44 radiomics features, 6 are considered optimal features with the absolute values of LASSO coefficients being greater than 0.1, including First Order _Uniformity, GLCM _Maximum Probability, GLCM _Inverse Difference, GLCM _ Informational Measure of Correlation 1, GLDM _Dependence Variance, and GLRLM _Long Run Emphasis. Histologically, ICC often shows aggressive trait, and inflammation, necrosis, and fibrosis are common (2); therefore, heterogeneous signal intensity in MR images are found to be more frequent for MVI-positive ICC, which is in concordance with First Order _Uniformity. Additionally, other 5 radiomics features indicate the higher the neighboring intensity value, variance, and gray-level value of VOI, the higher is the probability of MVI. Compared with the 100 radiomics features pool constructed by Zhou et al. (22), 42 of 44 radiomics features in our study are repeatable, and 3 of 7 optimal radiomics features in his study are also detected in our study. Both three pairwise combination models and the final radiomics model based on the diffusion-weighted images, precontrast T1-weighted images, and delayed phase images have solved the defect of overfit predictive performance in the training cohort and poor predictive performance in the validation cohort of single MR sequence models.

As far as we know, this study is the first to establish a nomogram based on clinicoradiologic and radiomics features. The MVI prediction model incorporates the final radiomics model, and imaging model exhibits excellent performance in both the training cohort (AUC = 0.953) and validation cohort (AUC = 0.861) with good calibration, which are better than the previous study (AUC_training = 0.873, AUC_validation = 0.850) (22). Also, the MVI prediction model achieves better efficacy than the imaging model in our study, suggesting radiomics features are indispensable in MVI prediction. More importantly, the nomogram in our study also exhibits a desirable prediction performance in the test cohort (AUC = 0.819) and performs better than clinical factor model constructed by Tang et al. (35) with AUCs of 0.739, 0.717, and 0.709 in training, validation, and test cohorts. Hence, the use of our nomogram preoperatively may be a noninvasive and robust method of assisting personalized treatment making, and patients with ICCs may suffer a higher net benefit from it.

There are some limitations in our study. Firstly, selection bias is inevitable in our retrospective study. Secondly, compared with MVI-negative ICCs, the number of MVI-positive ICCs is relatively small. Thirdly, the association between radiomics features and complex tumor biological features needs to be further explained. Fourthly, the study about whether the preoperative prediction of overall survival (OS) and disease-free survival (DFS) based on our MVI prediction model has the same effect as the postoperative prediction of OS and DFS based on pathological MVI after surgery is needed. Finally, although the nomogram has achieved a desirable prediction performance in the test cohort, larger cohorts from other centers are needed to be collected for the prospective validation of our nomogram.

In conclusion, radiomics features extracted from diffusion-weighted images, precontrast T1-weighted images, and delayed phase images of Ga-DTPA-enhanced MR imaging can assist in predicting MVI status of ICC patients. The MVI prediction nomogram incorporating radiomics features and imaging features including tumor size and intrahepatic duct dilatation is a potential biomarker and clinical tool in MVI stratification of ICC patients preoperatively.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by the Zhongshan Hospital, Fudan University ethics committee and Xuzhou Central Hospital ethics committee. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

Conceptualization: MZ and YS. Methodology: MZ, XQ, and FW. Software: XQ, XL, and FW. Validation: XM, YZ, and CZ. Formal analysis: XQ and XL. Resources: XQ, XL, and MZ. Data curation: XQ and XL. Writing—original draft preparation: XQ and XL. Writing—review and editing: MZ and YS. Supervision: MZ and YS. Project administration: MZ. Funding acquisition: MZ and YS. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

This study was supported by grants from the National Natural Science Foundation of China (No. 82171897), Shanghai Municipal Key Clinical Specialty (No. shslczdzk03202), Clinical Research Plan of SHDC (No. SHDC2020CR1029B), Clinical Research Project of Zhongshan Hospital, Fudan University (No. 2020ZSLC61), “the Development Research Foundation of Affiliated Hospital” of Xuzhou Medical University (No. XYFM2020020), and the Medical Scientific Research Program of Jiangsu Commission of Health (No. M2021014).

Conflict of Interest

Author FW is employed by Shanghai United Imaging Intelligence Co., Ltd. Shanghai, China.

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

Publisher’s Note

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

Acknowledgments

We sincerely thank Chun Yang for the assistance of radiomics analysis.

Supplementary Material

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

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sempoux C, Jibara G, Ward SC, Fan C, Qin L, Roayaie S, et al. Intrahepatic Cholangiocarcinoma: New Insights in Pathology. Semin Liver Dis (2011) 31(1):49–60. doi: 10.1055/s-0031-1272839

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Massarweh NN, El-Serag HB. Epidemiology of Hepatocellular Carcinoma and Intrahepatic Cholangiocarcinoma. Cancer Control (2017) 24(3):1073274817729245. doi: 10.1177/1073274817729245

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Brindley PJ, Bachini M, Ilyas SI, Khan SA, Loukas A, Sirica AE, et al. Cholangiocarcinoma. Nat Rev Dis Primers (2021) 7(1):65. doi: 10.1038/s41572-021-00300-2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Utada M, Ohno Y, Tamaki T, Sobue T, Endo G. Long-Term Trends in Incidence and Mortality of Intrahepatic and Extrahepatic Bile Duct Cancer in Japan. J Epidemiol (2014) 24(3):193–9. doi: 10.2188/jea.je20130122

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bergquist A, von Seth E. Epidemiology of Cholangiocarcinoma. Best Pract Res Clin Gastroenterol (2015) 29(2):221–32. doi: 10.1016/j.bpg.2015.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Everhart JE, Ruhl CE. Burden of Digestive Diseases in the United States Part III: Liver, Biliary Tract, and Pancreas. Gastroenterology (2009) 136(4):1134–44. doi: 10.1053/j.gastro.2009.02.038

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang K, Zhang H, Xia Y, Liu J, Shen F. Surgical Options for Intrahepatic Cholangiocarcinoma. Hepatobiliary Surg Nutr (2017) 6(2):79–90. doi: 10.21037/hbsn.2017.01.06

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Rahnemai-Azar AA, Weisbrod AB, Dillhoff M, Schmidt C, Pawlik TM. Intrahepatic Cholangiocarcinoma: Current Management and Emerging Therapies. Expert Rev Gastroenterol Hepatol (2017) 11(5):439–49. doi: 10.1080/17474124.2017.1309290

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tsukamoto M, Yamashita YI, Imai K, Umezaki N, Yamao T, Okabe H, et al. Predictors of Cure of Intrahepatic Cholangiocarcinoma After Hepatic Resection. Anticancer Res (2017) 37(12):6971–5. doi: 10.21873/anticanres.12164

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Cong WM, Bu H, Chen J, Dong H, Zhu YY, Feng LH, et al. Practice Guidelines for the Pathological Diagnosis of Primary Liver Cancer: 2015 Update. World J Gastroenterol (2016) 22(42):9279–87. doi: 10.3748/wjg.v22.i42.9279

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Surov A, Pech M, Omari J, Fischbach F, Damm R, Fischbach K, et al. Diffusion-Weighted Imaging Reflects Tumor Grading and Microvascular Invasion in Hepatocellular Carcinoma. Liver Cancer (2021) 10(1):10–24. doi: 10.1159/000511384

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Shao C, Chen J, Chen J, Shi J, Huang L, Qiu Y. Histological Classification of Microvascular Invasion to Predict Prognosis in Intrahepatic Cholangiocarcinoma. Int J Clin Exp Pathol (2017) 10(7):7674–81.

PubMed Abstract | Google Scholar

14. Zhou Y, Wang X, Xu C, Zhou G, Liu X, Gao S, et al. Mass-Forming Intrahepatic Cholangiocarcinoma: Can Diffusion-Weighted Imaging Predict Microvascular Invasion? J Magn Reson Imaging (2019) 50(1):315–24. doi: 10.1002/jmri.26566

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: The Bridge Between Medical Imaging and Personalized Medicine. Nat Rev Clin Oncol (2017) 14(12):749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Wang X, Wang S, Yin X, Zheng Y. MRI-Based Radiomics Distinguish Different Pathological Types of Hepatocellular Carcinoma. Comput Biol Med (2021) 141:105058. doi: 10.1016/j.compbiomed.2021.105058

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang Q, Li C, Zhang J, Hu X, Fan Y, Ma K, et al. Radiomics Models for Predicting Microvascular Invasion in Hepatocellular Carcinoma: A Systematic Review and Radiomics Quality Score Assessment. Cancers (Basel) (2021) 13(22):5864. doi: 10.3390/cancers13225864

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yang Y, Fan W, Gu T, Yu L, Chen H, Lv Y, et al. Radiomic Features of Multi-ROI and Multi-Phase MRI for the Prediction of Microvascular Invasion in Solitary Hepatocellular Carcinoma. Front Oncol (2021) 11:756216. doi: 10.3389/fonc.2021.756216

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ji GW, Zhu FP, Zhang YD, Liu XS, Wu FY, Wang K, et al. A Radiomics Approach to Predict Lymph Node Metastasis and Clinical Outcome of Intrahepatic Cholangiocarcinoma. Eur Radiol (2019) 29(7):3725–35. doi: 10.1007/s00330-019-06142-7

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Liang W, Xu L, Yang P, Zhang L, Wan D, Huang Q, et al. Novel Nomogram for Preoperative Prediction of Early Recurrence in Intrahepatic Cholangiocarcinoma. Front Oncol (2018) 8:360. doi: 10.3389/fonc.2018.00360

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Wang Y, Li J, Xia Y, Gong R, Wang K, Yan Z, et al. Prognostic Nomogram for Intrahepatic Cholangiocarcinoma After Partial Hepatectomy. J Clin Oncol (2013) 31(9):1188–95. doi: 10.1200/jco.2012.41.5984

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhou Y, Zhou G, Zhang J, Xu C, Wang X, Xu P. Radiomics Signature on Dynamic Contrast-Enhanced MR Images: A Potential Imaging Biomarker for Prediction of Microvascular Invasion in Mass-Forming Intrahepatic Cholangiocarcinoma. Eur Radiol (2021) 31(9):6846–55. doi: 10.1007/s00330-021-07793-1

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Lei Z, Li J, Wu D, Xia Y, Wang Q, Si A, et al. Nomogram for Preoperative Estimation of Microvascular Invasion Risk in Hepatitis B Virus-Related Hepatocellular Carcinoma Within the Milan Criteria. JAMA Surg (2016) 151(4):356–63. doi: 10.1001/jamasurg.2015.4257

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lewis S, Besa C, Wagner M, Jhaveri K, Kihira S, Zhu H, et al. Prediction of the Histopathologic Findings of Intrahepatic Cholangiocarcinoma: Qualitative and Quantitative Assessment of Diffusion-Weighted Imaging. Eur Radiol (2018) 28(5):2047–57. doi: 10.1007/s00330-017-5156-6

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Marrero JA, Kulik LM, Sirlin CB, Zhu AX, Finn RS, Abecassis MM, et al. Diagnosis, Staging, and Management of Hepatocellular Carcinoma: 2018 Practice Guidance by the American Association for the Study of Liver Diseases. Hepatology (2018) 68(2):723–50. doi: 10.1002/hep.29913

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Haradome H, Unno T, Morisaka H, Toda Y, Kwee TC, Kondo H, et al. Gadoxetic Acid Disodium-Enhanced MR Imaging of Cholangiolocellular Carcinoma of the Liver: Imaging Characteristics and Histopathological Correlations. Eur Radiol (2017) 27(11):4461–71. doi: 10.1007/s00330-017-4811-2

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Squadroni M, Tondulli L, Gatta G, Mosconi S, Beretta G, Labianca R. Cholangiocarcinoma. Crit Rev Oncol Hematol (2017) 116:11–31. doi: 10.1016/j.critrevonc.2016.11.012

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet (2014) 383(9935):2168–79. doi: 10.1016/s0140-6736(13)61903-0

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ma X, Liu L, Fang J, Rao S, Lv L, Zeng M, et al. MRI Features Predict Microvascular Invasion in Intrahepatic Cholangiocarcinoma. Cancer Imaging (2020) 20(1):40. doi: 10.1186/s40644-020-00318-x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Xu X, Zhang HL, Liu QP, Sun SW, Zhang J, Zhu FP, et al. Radiomic Analysis of Contrast-Enhanced CT Predicts Microvascular Invasion and Outcome in Hepatocellular Carcinoma. J Hepatol (2019) 70(6):1133–44. doi: 10.1016/j.jhep.2019.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Chong HH, Yang L, Sheng RF, Yu YL, Wu DJ, Rao SX, et al. Multi-Scale and Multi-Parametric Radiomics of Gadoxetate Disodium-Enhanced MRI Predicts Microvascular Invasion and Outcome in Patients With Solitary Hepatocellular Carcinoma ≤ 5 Cm. Eur Radiol (2021) 31(7):4824–38. doi: 10.1007/s00330-020-07601-2

PubMed Abstract | CrossRef Full Text | Google Scholar

32. 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(5):373–86. doi: 10.1159/000494099

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Spolverato G, Ejaz A, Kim Y, Sotiropoulos GC, Pau A, Alexandrescu S, et al. Tumor Size Predicts Vascular Invasion and Histologic Grade Among Patients Undergoing Resection of Intrahepatic Cholangiocarcinoma. J Gastrointest Surg (2014) 18(7):1284–91. doi: 10.1007/s11605-014-2533-1

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Oliveira IS, Kilcoyne A, Everett JM, Mino-Kenudson M, Harisinghani MG, Ganesan K. Cholangiocarcinoma: Classification, Diagnosis, Staging, Imaging Features, and Management. Abdom Radiol (NY) (2017) 42(6):1637–49. doi: 10.1007/s00261-017-1094-7

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Tang Z, Liu WR, Zhou PY, Ding ZB, Jiang XF, Wang H, et al. Prognostic Value and Predication Model of Microvascular Invasion in Patients With Intrahepatic Cholangiocarcinoma. J Cancer (2019) 10(22):5575–84. doi: 10.7150/jca.32199

PubMed Abstract | CrossRef Full Text | Google Scholar