Next Article in Journal
Cardiac Myxomas Resembling Malignant Neoplasia: Incidentally Diagnosed vs. Cerebral Embolized Myxomas
Next Article in Special Issue
Advanced Magnetic Resonance Imaging (MRI) Techniques: Technical Principles and Applications in Nanomedicine
Previous Article in Journal
RAB20 Promotes Proliferation via G2/M Phase through the Chk1/cdc25c/cdc2-cyclinB1 Pathway in Penile Squamous Cell Carcinoma
Previous Article in Special Issue
Longitudinal CT Imaging to Explore the Predictive Power of 3D Radiomic Tumour Heterogeneity in Precise Imaging of Mantle Cell Lymphoma (MCL)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contrast MR-Based Radiomics and Machine Learning Analysis to Assess Clinical Outcomes following Liver Resection in Colorectal Liver Metastases: A Preliminary Study

1
Division of Radiology, Istituto Nazionale Tumori IRCCS Fondazione Pascale–IRCCS di Napoli, 80131 Naples, Italy
2
Medical Oncology Division, Igea SpA, 80013 Napoli, Italy
3
Department of Medicine and Health Sciences V. Tiberio, University of Molise, 86100 Campobasso, Italy
4
Department of Medicine, Surgery and Dentistry, University of Salerno, 84084 Salerno, Italy
5
Division of Radiology, Università Degli Studi Della Campania Luigi Vanvitelli, 80138 Naples, Italy
6
Division of Abdominal Oncology, Istituto Nazionale Tumori IRCCS Fondazione Pascale–IRCCS di Napoli, 80131 Naples, Italy
7
Division of Radiology, Azienda Ospedaliera Universitaria Careggi, 50134 Florence, Italy
8
Italian Society of Medical and Interventional Radiology (SIRM), SIRM Foundation, Via della Signora 2, 20122 Milan, Italy
9
Division of Epatobiliary Surgical Oncology, Istituto Nazionale Tumori IRCCS Fondazione Pascale–IRCCS di Napoli, 80131 Naples, Italy
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(5), 1110; https://doi.org/10.3390/cancers14051110
Submission received: 7 January 2022 / Revised: 17 February 2022 / Accepted: 21 February 2022 / Published: 22 February 2022
(This article belongs to the Special Issue Radiology and Imaging of Cancer)

Abstract

:

Simple Summary

The objective of the study was to evaluate the radiomics features obtained by contrast MRI studies as prognostic biomarkers in colorectal liver metastases patients to predict clinical outcomes following liver resection. We demonstrated a good performance considering the single textural significant metric in the identification of front of tumor growth (expansive versus infiltrative) and tumor budding (high grade versus low grade or absent), in the recognition of mucinous type and in the detection of recurrences. Moreover, considering linear regression models or neural network classifiers in a multivariate approach was possible to increase the performance in terms of accuracy, sensitivity, and specificity.

Abstract

Purpose: To assess radiomics features efficacy obtained by arterial and portal MRI phase in the prediction of clinical outcomes in the colorectal liver metastases patients, evaluating recurrence, mutational status, pathological characteristic (mucinous and tumor budding) and surgical resection margin. Methods: This retrospective analysis was approved by the local Ethical Committee board, and radiological databases were used to select patients with colorectal liver metastases with pathological proof and MRI study in a pre-surgical setting after neoadjuvant chemotherapy. The cohort of patients included a training set (51 patients with 61 years of median age and 121 liver metastases) and an external validation set (30 patients with single lesion with 60 years of median age). For each segmented volume of interest on MRI by two expert radiologists, 851 radiomics features were extracted as median values using the PyRadiomics tool. Non-parametric Kruskal-Wallis test, intraclass correlation, receiver operating characteristic (ROC) analysis, linear regression modelling and pattern recognition methods (support vector machine (SVM), k-nearest neighbors (KNN), artificial neural network (NNET), and decision tree (DT)) were considered. Results: The best predictor to discriminate expansive versus infiltrative tumor growth front was wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis extracted on portal phase with accuracy of 82%, sensitivity of 84%, and specificity of 77%. The best predictor to discriminate tumor budding was wavelet_LLH_firstorder_10Percentile extracted on portal phase with accuracy of 92%, a sensitivity of 96%, and a specificity of 81%. The best predictor to differentiate the mucinous type of tumor was the wavelet_LLL_glcm_ClusterTendency extracted on portal phase with accuracy of 88%, a sensitivity of 38%, and a specificity of 100%. The best predictor to identify the recurrence was the wavelet_HLH_ngtdm_Complexity extracted on arterial phase with accuracy of 90%, a sensitivity of 71%, and a specificity of 95%. The best linear regression model was obtained in the identification of mucinous type considering the 13 textural significant metrics extracted by arterial phase (accuracy of 94%, sensitivity of 77% and a specificity of 99%). The best results were obtained in the identification of tumor budding with the eleven textural significant features extracted by arterial phase using a KNN (accuracy of 95%, sensitivity of 84%, and a specificity of 99%). Conclusions: Our results confirmed the capacity of radiomics to identify as biomarkers and several prognostic features that could affect the treatment choice in patients with liver metastases in order to obtain a more personalized approach.

1. Introduction

Radiomics is a promising area that investigates the capability of quantitative features extracted by medical images as biomarkers to assess the biology of pathological processes at microscopic levels. These data can be converted into image-based marks to spread diagnostic, prognostic and predictive accuracy in oncological setting [1,2,3,4,5,6,7,8]. Radiomics could theoretically support tumor detection, evaluation of prognosis, estimate treatment response [9,10,11,12,13,14]. Radiomics is designed to be used in decision support of precision medicine, using standard of care images that are routinely acquired in clinical practice. It presents a cost-effective and highly feasible addition for clinical decision support. Moreover, this analysis non-invasively characterize the overall tumor accounting for heterogeneity, interrogating the entire tumor allows the expression of microscopic genomic and proteomics patterns in terms of macroscopic image-based features [15,16,17,18]. Moreover, this analysis gives prognostic and/or predictive biomarker allowing for a fast, low-cost, and repeatable tool for longitudinal monitoring [19,20].
The association of radiomics and molecular features, so named radiogenomics, shows clear effects for management of cancer patients. Although several studies have assessed the rule of radiogenomics in hepatocellular carcinoma, only a few have tested the radiomics rule in colorectal cancer metastatic lesions in the liver [1,2,3]. Today, radiologists play an important role in the multidisciplinary team of colorectal patients with liver metastases. During the staging and surveillance phase, it is critical to identify all liver lesions, since this is related to proper patient management. Additionally, after conversion therapy, all lesions assessed at first exam should be re-evaluated to identify responders and non-responders as soon as possible [21,22,23,24,25]. Although computed tomography (CT) is habitually the diagnostic tool employed for staging and surveillance, magnetic resonance imaging (MRI) is a valuable diagnostic technique in oncologic settings, since it is the only technique that allows evaluating of morphological and functional features of tumor status, providing quantitative parameters that improve the characterization of a lesion and the assessment after therapy [21,22,23,24]. Moreover, several liver-specific contrast agents have been inserted to improve the hepatic lesions detection and characterization. Gadobenate dimeglumine (Gd-BOPTA) and gadolinium ethoxybenzyl diethylenetria-mine pentaacetic acid (Gd-EOB-DTPA) allow obtaining information about the vascularization of lesions in the different phases of contrast circulation and functional parameters in the delayed hepatobiliary phase (EOB-phase).
In this scenario, the possibility to compare radiomic data extracted by MRI in the identification of recurrence, mutational status, pathological characteristics (mucinous and tumor budding), and surgical resection margin could provide significant benefits respect to qualitative evaluation. In fact, radiomics predictors could permit an effective treatment selection in the perspective of personalized medicine, treatment response prediction, in the differentiation of favorable subsets of patients from those with poor prognosis, and selecting patients that may benefit from surgical treatment. In the present study, we assessed the radiomics features efficacy obtained by contrast (arterial and portal phase) MRI to predict clinical outcomes following liver resection in colorectal liver metastases patients.

2. Materials and Methods

2.1. Dataset Characteristics

This study was approved by the local ethical committee board that renounces the patient informed consent due to the retrospective nature of the study. The study was performed in accordance with relevant guidelines and regulations.
Patient selection was made considering internal radiological databases in the period from January 2018 to May 2021 using the following criteria: (1) liver pathological proven metastases; (2) contrast MRI study in pre-surgical setting after neoadjuvant chemotherapy; (3) MR images of high quality and (4) a follow-up CT scan of at least six months after surgery. The exclusion criteria were: (1) discordance among the imaging diagnosis and the pathological ones, (2) no contrast MRI studies and (c) no high-quality MR images.
The analyzed patients included a training set and an external validation set. The internal training set included 51 patients (18 women and 33 men) with 61 years of median age (range 35–82 years) and 121 liver metastases. The validation cohort, provided by “Careggi Hospital”, Florence, Italy, consisted of a total of 30 patients with single lesion (10 women and 20 men) with 60 years of median age (range 40–78 years). The patient characteristics are summarized in Table 1.

2.2. MR Imaging Protocol

A Magnetom Symphony 1.5 T scanner (Siemens, Erlangen, Germany) and a Magnetom Aera (Siemens) 1.5 T scanner equipped with an 8-element body and phased array coils were used to acquire an MRI study that includes sequences before and after intravenous (IV) contrast agent (CA) injection.
In this study, radiomic features extraction was made on volumetric interpolated breath-hold examination (VIBE) T1-weighted SPAIR with controlled respiration used to acquire images after IV CA injection with a liver-specific CA (0.1 mL/kg of Gd-EOB-BPTA, Primovist, Bayer Schering Pharma, Berlin, Germany) as descripted in [26,27].
The VIBE T1-W sequence was acquired with two different flip angles (10 and 30 degrees). A power injector (Spectris Solaris® EP MR, MEDRAD, Inc., Indianola, IA, USA) was used to administer the CA at an infusion rate of 2 mL/s. VIBE T1-w images were acquired in four different phases: arterial phase (35 s delay), portal venous phase (90 s), late/transitional phase (120 s), and hepatobiliary excretion phase (20 min). MRI protocol details are reported in Table 2.

2.3. Follow-Up CT Scan

CT studies were performed using a scanner with 64 detectors (Optima 660, GE Healthcare, Chicago, IL, USA). The scan data was 120 kVp, 100–470 mA (NI 16.36), slice thickness was 2.5 mm, and table speed/rotation was 0.984/1 mm. The liver protocol included a quadruple phase protocol, counting unenhanced, arterial, portal, and equilibrium phases. A non-ionic contrast agent (120 mL of iomeprol, Iomeron 400, Bracco, Milan, Italy) was injected at a rate of 3 mL/s using an automatic power injector (Empower CTA, EZ-EM Inc., New York, NY, USA). The arterial phase was started 19 s after the descending aorta attenuation reached 100 HU, measured by the bolus localization method.

2.4. Image Processing

Regions of interest (ROIs) were manually drawn slice-by-slice by two expert radiologists with 22 and 15 years of abdominal imaging experience, respectively, first separately and then together and in accordance with each other. Region of interest segmentation was performed avoiding encircling any distortion artefacts. For each volume of interest, radiomics features were extracted as median values, reducing the possible influence by artefacts.
The segmentation was performed on arterial phase and portal phase of VIBE T1-W_FA10. For these reasons, we obtained the results both on arterial phase volume and on portal phase volume.
Manual definition of the ROIs was made using the segmentation tool of 3DSlicer 4.11 (Figure 1) [https://www.slicer.org/, accessed on 20 December 2021].

2.5. MRI Post-Processing with Pyradiomic Tool

Eight hundred fifty-one radiomic features were extracted using PyRadiomics v3.0.1 [28] and included first-order statistics, shape-based (3D) metrics, shape-based (2D) features, gray level co-occurence matrix features, gray level run length matrix features, gray level size zone matrix features, neighboring gray tone difference matrix features and gray level dependence matrix parameters. The extracted features are in compliance with feature definitions as described by the Imaging Biomarker Standardization Initiative (IBSI) [29] and reported in [https://readthedocs.org/projects/pyradiomics/downloads/, accessed on 20 December 2021]. Radiomics analysis was performed blinded to the clinical and pathological data.

2.6. Statistical Analysis

Statistical analysis included univariate and multivariate approaches.

2.6.1. Univariate Analysis

The observer variability assessment was performed by calculating the intraclass correlation coefficient.
A non-parametric Kruskal-Wallis test was performed to identify statistically significant differences among clinical parameters and radiomic metrics of two groups (front of tumor growth: expansive versus infiltrative; tumor budding: high-grade versus low-grade or absent; mucinous type; and presence of recurrence).
Receiver operating characteristic (ROC) analysis was performed, and the Youden index was used to individuate the optimal cutoff value for each feature in order to calculate area under the ROC curve (AUC), sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV) and accuracy.
The McNemar test was used to calculate statistically significant differences among dichotomy data of the performance results.

2.6.2. Multivariate Analysis

A multivariate analysis was performed in order to identify the combinations of variables which best predict the outcomes: (1) front of tumor growth: expansive versus infiltrative; (2) tumor budding: high-grade versus low-grade or absent; (3) mucinous type; and (4) presence of recurrence.
Given the high number of textural features, a first selection of variables was made based on the results obtained from the univariate analysis (Table 3). Therefore, there was no waiting for overfitting in our study because adequate feature selection was made according to sample size.
Linear regression modelling was used to assess the best linear combination of features considered as predictors for each outcome (Table 3). The linear regression model was used to assess the accuracy of linear combination, and ROC analysis with Youden index was used to identify the optimal cut-off value, sensitivity, specificity, PPV, and NPV.
Moreover, pattern recognition methods were used in the context of a multivariate artificial intelligence approach. The tested classifiers with a 10-k fold cross-validation were support vector machine (SVM), k-nearest neighbors (KNN), artificial neural network (NNET), and decision tree (DT)). A description of classifiers can be found in [30]. The best classifier was chosen considering the highest area under ROC curve and highest accuracy. An external validation cohort was used to validate the findings of the best classifier found in the training step.
The statistical analyses were performed using the Statistics Toolbox and Machine Toolbox of MATLAB R2021b (MathWorks, Natick, MA, USA) considering a p value ≤ 0.05 as significant.

3. Results

3.1. Univariate Analysis Findings

The intraclass correlation coefficients median value for extracted features was 0.94 (range 0.88–0.98). The size of the lesion did not affect the values of the extracted metrics (p-value > 0.05 at the Kruskal-Wallis test considering lesions < 2 cm and ≥ 2 cm).
Among significant features to differentiate the tumor growth front in the arterial phase, 7 textural parameters obtained an accuracy ≥ 75% Among these 7 features, the best performance to discriminate expansive versus infiltrative front of tumor growth was obtained by the wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis with accuracy of 79%, sensitivity of 95%, specificity of 51%, PPV and NPV of 77% and 85%, respectively, and a cut-off value of 0.12 (Table 4).
Among significant features to differentiate the front of tumor growth in portal phase, 9 textural parameters obtained an accuracy ≥ 80%. Among these 9 features, the best performance to discriminate expansive versus infiltrative front of tumor growth was obtained by the wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis (the same feature of previous case) with accuracy of 82%, sensitivity of 84%, specificity of 77%, PPV and NPV of 85% and 74%, respectively, and a cut-off value of 0.12 (Table 4).
Among significant features to differentiate the tumor budding on arterial phase, 11 textural parameters obtained an accuracy ≥ 80%. Among these 11 features, the best performance to discriminate tumor budding was obtained by the wavelet_LHH_firstorder_Minimum with accuracy of 86%, sensitivity of 98%, specificity of 52%, PPV and NPV of 85% and 89%, respectively, and a cut-off value of −41.76 (Table 4).
Among significant features to differentiate the tumor budding in the portal phase, 13 textural parameters obtained an accuracy ≥ 85%. Among these 13 features, the best performance to discriminate tumor budding was obtained by the wavelet_LLH_firstorder_10Percentile with accuracy of 92%, sensitivity of 96%, specificity of 81%, PPV and NPV of 93% and 86%, respectively, and a cut-off value of −37.14 (Table 4).
Among significant features to differentiate the mucinous type of tumor in the arterial phase, 13 textural parameters obtained an accuracy ≥ 80%. Among these 13 features, the best performance to differentiate the mucinous type of tumor was obtained by the wavelet_HLH_glszm_LargeAreaHighGrayLevelEmphasis with accuracy of 85%, sensitivity of 35%, specificity of 99%, PPV and NPV of 90% and 85%, respectively, and a cut-off value of −0.02 (Table 4).
Among significant features to differentiate the mucinous type of tumor in the portal phase, 12 textural parameters obtained an accuracy ≥ 85%. Among these 12 features, the best performance to differentiate the mucinous type of tumor was obtained by the wavelet_LLL_glcm_ClusterTendency with accuracy of 88%, sensitivity of 38%, specificity of 100%, PPV and NPV of 100% and 86%, respectively, and a cut-off value of 408.22 (Table 4).
Among significant features to identify tumor recurrence in the arterial phase, 10 textural parameters obtained an accuracy ≥ 80%. Among these 10 features, the best performance to identify tumor recurrence was obtained by the wavelet_HLH_ngtdm_Complexity with accuracy of 90%, sensitivity of 71%, specificity of 95%, PPV and NPV of 79% and 90%, respectively, and a cut-off value of 3.34 (Table 4).
Among significant features to identify tumor recurrence in the portal phase, 11 textural parameters obtained an accuracy ≥ 85%. Among these 11 features, the best performance to identify tumor recurrence was obtained by the LHL_glcm_Correlation with accuracy of 89%, sensitivity of 71%, specificity of 94%, PPV and NPV of 81% and 90%, respectively, and a cut-off value of 1.54 (Table 4).

3.2. Multivariate Analysis Findings

3.2.1. Linear Regression Analysis Findings

Linear regression models obtained good results in each considered classification problem (1. Front of tumor growth: expansive versus infiltrative; 2. tumor budding: high-grade versus low-grade or absent; 3. mucinous type and 4. presence of recurrence) with accuracy in the range of 83−94% (Table 5 and Table 6, Figure 2 and Figure 3). The best linear regression model was obtained in the identification of mucinous type considering the 13 textural significant metrics extracted by the arterial phase (AUC of 0.93, accuracy of 94%, sensitivity of 77%, and specificity of 99%) and in the identification of tumor budding considering the 11 textural significant metrics extracted by the arterial phase (AUC of 0.92, accuracy of 93%, sensitivity of 94%, and specificity of 90%).
The coefficients of these linear models are reported in the Table 7.

3.2.2. Pattern Recognition Approaches Findings

Considering significant texture metrics tested with pattern recognition approaches, the best performance for each outcome (1. front of tumor growth: expansive versus infiltrative; 2. tumor budding: high-grade versus low-grade or absent; 3. mucinous type and 4. presence of recurrence) was reached by a KNN as classifier given the features extracted by the arterial phase (Figure 4). Instead, considering the features extracted by the portal phase, the best performance was reached by a KNN as classifier in the identification of the font of tumor growth, mucinous type, and for the detection of recurrences, and by a decision tree for the tumor budding identification (Table 5 and Table 6).
The accuracy was always greater than 88% (Table 5 and Table 6, Figure 5) in both the training and validation sets, and the best results were obtained in the identification of tumor budding, with the eleven textural significant features extracted by the arterial phase (AUC of 0.95, accuracy of 95%, sensitivity of 84%, and specificity of 99%).
The best performance in terms of accuracy obtained considering the extracted radiomics features on arterial phase was significantly superior at the best performance obtained in the portal phase (p value < 0.05 at McNemar test).

4. Discussion

The present study confirmed the possibility of radiomics to recognize as biomarkers several features that could influence the treatment choice in patients with liver metastases in order to obtain a more personalized approach. Our data were verified by an external validation dataset. We obtained a good performance considering the single textural significant metric in the identification of front of tumor growth (expansive versus infiltrative) and tumor budding (high-grade versus low-grade or absent), in the recognition of mucinous type, and in the detection of recurrences.
At univariate analysis, the best predictor to discriminate expansive versus infiltrative tumor growth front was wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis by portal phase with accuracy of 82%. The best predictor to discriminate tumor budding was wavelet_LLH_firstorder_10Percentile by portal phase with accuracy of 92%. The best predictor to differentiate the mucinous type was the wavelet_LLL_glcm_ClusterTendency by portal phase with accuracy of 88%. The best predictor to identify the recurrence was wavelet_HLH_ngtdm_Complexity by arterial phase with accuracy of 90%.
The best linear regression model was obtained in the identification of mucinous type considering the 13 textural significant metrics extracted by the arterial phase (AUC of 0.93, accuracy of 94%, sensitivity of 77%, and specificity of 99%). The best results were obtained in the identification of tumor budding with the eleven textural significant features extracted by the arterial phase using a KNN (AUC of 0.95, accuracy of 95%, sensitivity of 84%, and specificity of 99%). Therefore, the best performance was reached considering the radiomics features extracted on arterial phase.
Several studies demonstrated the correlation between radiomics parameters and prognosis [31,32,33,34,35,36,37,38,39,40]. An association between homogeneity and worse overall survival (OS) was demonstrated by Andersen et al. [32]. According to Rahmim et al., radiomic parameters of heterogeneity obtained by FDG PET were predictors of lower OS [37]. Lubner et al. demonstrated that the degree of skewness was inversely correlated to KRAS, while the entropy was related with OS [34]. In addition to the survival advantages, the possibility to predict recurrence in the liver has been demonstrated [37,38,39,40]. According to our results, Ravanelli et al. related high CT uniformity and low OS and PFS in patients with CRC and liver metastasis [39].
Radiomics and radiogenomics are emerging tools with significant limits. The major weakness is the heterogeneity of software employed in different investigations due to the variety of imaging devices in different hospitals. This clearly hinders the interpretation of different data for multicenter studies. In addition, the segmentation may affect results [41].
Many previous studies have shown that the lower the degree of differentiation of the primary tumor (mainly manifested by more aggressive and malignant phenotype), the worse the prognosis; the number of metastases (≤ 4) and the diameter of metastases (≤5 cm) are the main prognostic factors affecting the prognosis of patients with liver metastases [42,43,44]. Ma et al. showed that the lower the degree of tumor differentiation, the greater the number of metastases, and the larger the diameter of the metastases, the heavier the tumor burden throughout the body, and the shorter the survival period. The distance between metastases and great vessels could affect the recurrence-free survival of patients [45]. However, the prognostic and predictive value of radiomics in colorectal cancer metastases to the liver have been well studied, demonstrating its higher utility in predicting clinical outcomes compared to other clinical data [2].
The present study had several limitations: (1) the small size of the population considered, although the analysis was performed on a homogeneous sample and on all individual lesions; furthermore the patients analyzed included a training set and an external validation set for a total of 151 liver lesions analyzed. The external dataset was used to validate the results obtained during the training phase; however, the results of this study were considered preliminary, and the future goal is to broaden the dataset to evaluate the generalization of the results; (2) the retrospective nature of the study; (3) a manual segmentation, that, although several studies support automatic segmentation to avoid inter-observer variability, in our opinion, the manual approach is more realistic. Moreover, we did not assess the impact of chemotherapy. However, we assessed the impact of the different phases of contrast study (arterial and portal), although we have not evaluated these results with respect to transitional and EOB-phase due to morphological sequences, such as T2-weigthed, T1-weigthed, or diffusion-weighted imaging. Data that we plan to evaluate in a future study.

5. Conclusions

Our results confirmed the capacity of radiomics to identify as biomarkers several prognostic features that could affect the treatment choice in patients with liver metastases in order to obtain a more personalized approach. We obtained a good performance considering the single textural significant metric in the identification of front of tumor growth (expansive versus infiltrative) and tumor budding (high-grade versus low-grade or absent), in the recognition of mucinous type and in the detection of recurrences.

Author Contributions

Data curation, V.G.; Formal analysis, R.F. and C.C.; Investigation, V.G., F.D.M., S.V.S., F.d.A., A.O., A.A., G.N., F.G., V.P., V.M., L.B., F.I. and A.P.; Methodology, V.G., R.F., F.D.M., S.V.S., F.d.A., A.O., A.A., G.N., F.G., V.P., V.M., L.B., F.I. and A.P.; Project administration, C.C.; Writing–original draft, V.G.; Writing–review and editing, V.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the Local Ethical Committee board that renounces the patient informed consent due to the retrospective nature of the study (Deliberation N.8 of 3 January 2019 of National Cancer Institute of Naples Pascale Foundation). The study was performed in accordance with relevant guidelines and regulations.

Informed Consent Statement

Patient consent was waived due to the retrospective nature of the study.

Data Availability Statement

Data and material are available at https://zenodo.org/deposit/6198703.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Granata, V.; Fusco, R.; Costa, M.; Picone, C.; Cozzi, D.; Moroni, C.; La Casella, G.; Montanino, A.; Monti, R.; Mazzoni, F.; et al. Preliminary Report on Computed Tomography Radiomics Features as Biomarkers to Immunotherapy Selection in Lung Adenocarcinoma Patients. Cancers 2021, 13, 3992. [Google Scholar] [CrossRef] [PubMed]
  2. Granata, V.; Fusco, R.; Barretta, M.L.; Picone, C.; Avallone, A.; Belli, A.; Patrone, R.; Ferrante, M.; Cozzi, D.; Grassi, R.; et al. Radiomics in hepatic metastasis by colorectal cancer. Infect. Agents Cancer 2021, 16, 1–9. [Google Scholar] [CrossRef]
  3. Fusco, R.; Piccirillo, A.; Sansone, M.; Granata, V.; Rubulotta, M.; Petrosino, T.; Barretta, M.; Vallone, P.; Di Giacomo, R.; Esposito, E.; et al. Radiomics and Artificial Intelligence Analysis with Textural Metrics Extracted by Contrast-Enhanced Mammography in the Breast Lesions Classification. Diagnostics 2021, 11, 815. [Google Scholar] [CrossRef] [PubMed]
  4. Fusco, R.; Granata, V.; Mazzei, M.A.; Di Meglio, N.; Del Roscio, D.; Moroni, C.; Monti, R.; Cappabianca, C.; Picone, C.; Neri, E.; et al. Quantitative imaging decision support (QIDSTM) tool consistency evaluation and radiomic analysis by means of 594 metrics in lung carcinoma on chest CT scan. Cancer Control 2021, 28, 1073274820985786. [Google Scholar] [CrossRef] [PubMed]
  5. Granata, V.; Fusco, R.; Avallone, A.; De Stefano, A.; Ottaiano, A.; Sbordone, C.; Brunese, L.; Izzo, F.; Petrillo, A. Radiomics-Derived Data by Contrast Enhanced Magnetic Resonance in RAS Mutations Detection in Colorectal Liver Metastases. Cancers 2021, 13, 453. [Google Scholar] [CrossRef] [PubMed]
  6. Granata, V.; Fusco, R.; Risi, C.; Ottaiano, A.; Avallone, A.; De Stefano, A.; Grimm, R.; Grassi, R.; Brunese, L.; Izzo, F.; et al. Diffusion-Weighted MRI and Diffusion Kurtosis Imaging to Detect RAS Mutation in Colorectal Liver Metastasis. Cancers 2020, 12, 2420. [Google Scholar] [CrossRef]
  7. Petralia, G.; Summers, P.E.; Agostini, A.; Ambrosini, R.; Cianci, R.; Cristel, G.; Calistri, L.; Colagrande, S. Dynamic contrast-enhanced MRI in oncology: How we do it. Radiol. Med. 2020, 125, 1288–1300. [Google Scholar] [CrossRef]
  8. Ria, F.; Samei, E. Is regulatory compliance enough to ensure excellence in medicine? Radiol. Med. 2020, 125, 904–905. [Google Scholar] [CrossRef]
  9. Zhang, A.; Song, J.; Ma, Z.; Chen, T. Combined dynamic contrast-enhanced magnetic resonance imaging and diffusion-weighted imaging to predict neoadjuvant chemotherapy effect in FIGO stage IB2–IIA2 cervical cancers. Radiol. Med. 2020, 125, 1233–1242. [Google Scholar] [CrossRef]
  10. Crimì, F.; Capelli, G.; Spolverato, G.; Bao, Q.R.; Florio, A.; Rossi, S.M.; Cecchin, D.; Albertoni, L.; Campi, C.; Pucciarelli, S.; et al. MRI T2-weighted sequences-based texture analysis (TA) as a predictor of response to neoadjuvant chemo-radiotherapy (nCRT) in patients with locally advanced rectal cancer (LARC). Radiol. Med. 2020, 125, 1216–1224. [Google Scholar] [CrossRef]
  11. Kirienko, M.; Ninatti, G.; Cozzi, L.; Voulaz, E.; Gennaro, N.; Barajon, I.; Ricci, F.; Carlo-Stella, C.; Zucali, P.; Sollini, M.; et al. Computed tomography (CT)-derived radiomic features differentiate prevascular mediastinum masses as thymic neoplasms versus lymphomas. Radiol. Med. 2020, 125, 951–960. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, L.; Kang, L.; Li, G.; Zhang, X.; Ren, J.; Shi, Z.; Li, J.; Yu, S. Computed tomography-based radiomics model for discriminating the risk stratification of gastrointestinal stromal tumors. Radiol. Med. 2020, 125, 465–473. [Google Scholar] [CrossRef] [PubMed]
  13. Gurgitano, M.; Angileri, S.A.; Rodà, G.M.; Liguori, A.; Pandolfi, M.; Ierardi, A.M.; Wood, B.J.; Carrafiello, G. Interventional Radiology ex-machina: Impact of Artificial Intelligence on practice. Radiol. Med. 2021, 126, 998–1006. [Google Scholar] [CrossRef] [PubMed]
  14. Scapicchio, C.; Gabelloni, M.; Barucci, A.; Cioni, D.; Saba, L.; Neri, E. A deep look into radiomics. Radiol. Med. 2021, 126, 1296–1311. [Google Scholar] [CrossRef]
  15. Wei, J.; Jiang, H.; Gu, D.; Niu, M.; Fu, F.; Han, Y.; Song, B.; Tian, J. Radiomics in liver diseases: Current progress and future opportunities. Liver Int. 2020, 40, 2050–2063. [Google Scholar] [CrossRef]
  16. Saini, A.; Breen, I.; Pershad, Y.; Naidu, S.; Knuttinen, M.G.; Alzubaidi, S.; Sheth, R.; Albadawi, H.; Kuo, M.; Oklu, R. Radiogenomics and Radiomics in Liver Cancers. Diagnostics 2018, 9, 4. [Google Scholar] [CrossRef] [Green Version]
  17. De la Pinta, C.; Castillo, M.E.; Collado, M.; Galindo-Pumariño, C.; Peña, C. Radiogenomics: Hunting Down Liver Metastasis in Colorectal Cancer Patients. Cancers 2021, 13, 5547. [Google Scholar] [CrossRef]
  18. Nardone, V.; Reginelli, A.; Grassi, R.; Boldrini, L.; Vacca, G.; D’Ippolito, E.; Annunziata, S.; Farchione, A.; Belfiore, M.P.; Desideri, I.; et al. Delta radiomics: A systematic review. Radiol. Med. 2021, 126, 1–13. [Google Scholar] [CrossRef]
  19. Brunese, L.; Brunese, M.C.; Carbone, M.; Ciccone, V.; Mercaldo, F.; Santone, A. Automatic PI-RADS assignment by means of formal methods. Radiol. Med. 2021, 127, 83–89. [Google Scholar] [CrossRef]
  20. Van der Lubbe, M.F.J.A.; Vaidyanathan, A.; de Wit, M.; Burg, E.L.V.D.; Postma, A.A.; Bruintjes, T.D.; Bilderbeek-Beckers, M.A.L.; Dammeijer, P.F.M.; Bossche, S.V.; Van Rompaey, V.; et al. A non-invasive, automated diagnosis of Menière’s disease using radiomics and machine learning on conventional magnetic resonance imaging: A multicentric, case-controlled feasibility study. Radiol. Med. 2021, 127, 72–82. [Google Scholar] [CrossRef]
  21. Granata, V.; Fusco, R.; Avallone, A.; Cassata, A.; Palaia, R.; Delrio, P.; Grassi, R.; Tatangelo, F.; Grazzini, G.; Izzo, F.; et al. Abbreviated MRI protocol for colorectal liver metastases: How the radiologist could work in pre surgical setting. PLoS ONE 2020, 15, e0241431. [Google Scholar] [CrossRef] [PubMed]
  22. Granata, V.; Fusco, R.; Setola, S.V.; Raso, M.M.; Avallone, A.; De Stefano, A.; Nasti, G.; Palaia, R.; DelRio, P.; Petrillo, A.; et al. Liver radiologic findings of chemotherapy-induced toxicity in liver colorectal metastases patients. Eur. Rev. Med. Pharmacol. Sci. 2019, 23, 9697–9706. [Google Scholar] [PubMed]
  23. Granata, V.; Fusco, R.; Maio, F.; Avallone, A.; Nasti, G.; Palaia, R.; Albino, V.; Grassi, R.; Izzo, F.; Petrillo, A. Qualitative assessment of EOB-GD-DTPA and Gd-BT-DO3A MR contrast studies in HCC patients and colorectal liver metastases. Infect. Agents Cancer 2019, 14, 1–9. [Google Scholar] [CrossRef] [PubMed]
  24. Granata, V.; Fusco, R.; Castelguidone, E.D.L.D.; Avallone, A.; Palaia, R.; Delrio, P.; Tatangelo, F.; Botti, G.; Grassi, R.; Izzo, F.; et al. Diagnostic performance of gadoxetic acid–enhanced liver MRI versus multidetector CT in the assessment of colorectal liver metastases compared to hepatic resection. BMC Gastroenterol. 2019, 19, 1–11. [Google Scholar] [CrossRef] [Green Version]
  25. Granata, V.; Fusco, R.; Avallone, A.; Catalano, O.; Piccirillo, M.; Palaia, R.; Nasti, G.; Petrillo, A.; Izzo, F. A radiologist’s point of view in the presurgical and intraoperative setting of colorectal liver metastases. Futur. Oncol. 2018, 14, 2189–2206. [Google Scholar] [CrossRef]
  26. Granata, V.; Fusco, R.; Catalano, O.; Avallone, A.; Palaia, R.; Botti, G.; Tatangelo, F.; Granata, F.; Cascella, M.; Izzo, F.; et al. Diagnostic accuracy of magnetic resonance, computed tomography and contrast enhanced ultrasound in radiological multimodality assessment of peribiliary liver metastases. PLoS ONE 2017, 12, e0179951. [Google Scholar] [CrossRef] [Green Version]
  27. Granata, V.; Fusco, R.; Catalano, O.; Filice, S.; Amato, D.M.; Nasti, G.; Avallone, A.; Izzo, F.; Petrillo, A. Early Assessment of Colorectal Cancer Patients with Liver Metastases Treated with Antiangiogenic Drugs: The Role of Intravoxel Incoherent Motion in Dif-fusion-Weighted Imaging. PLoS ONE 2015, 10, e0142876. [Google Scholar] [CrossRef]
  28. Radiomics Community. Radiomic Features. 2016. Available online: https://pyradiomics.readthedocs.io/en/latest/features.html (accessed on 22 December 2021).
  29. Zwanenburg, A.; Vallières, M.; Abdalah, M.A.; Aerts, H.J.W.L.; Andrearczyk, V.; Apte, A.; Ashrafinia, S.; Bakas, S.; Beukinga, R.J.; Boellaard, R.; et al. The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Im-age-based Phenotyping. Radiology 2020, 295, 328–338. [Google Scholar] [CrossRef] [Green Version]
  30. Fusco, R.; Sansone, M.; Filice, S.; Carone, G.; Amato, D.M.; Sansone, C.; Petrillo, A. Pattern Recognition Approaches for Breast Cancer DCE-MRI Classification: A Systematic Review. J. Med. Biol. Eng. 2016, 36, 449–459. [Google Scholar] [CrossRef] [Green Version]
  31. Beckers, R.C.J.; Trebeschi, S.; Maas, M.; Schnerr, R.S.; Sijmons, J.M.L.; Beets, G.L.; Houwers, J.B.; Beets-Tan, R.G.H.; Lambregts, D.M.J. CT texture analysis in colorectal liver metastases and the surrounding liver parenchyma and its potential as an imaging biomarker of disease aggressiveness, response and survival. Eur. J. Radiol. 2018, 102, 15–21. [Google Scholar] [CrossRef]
  32. Andersen, I.R.; Thorup, K.; Andersen, M.B.; Olesen, R.; Mortensen, F.V.; Nielsen, D.T.; Rasmussen, F. Texture in the monitoring of regorafenib therapy in patients with colorectal liver metastases. Acta Radiol. 2019, 60, 1084–1093. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, H.; Li, W.; Hu, F.; Sun, Y.; Hu, T.; Tong, T. MR texture analysis: Potential imaging biomarker for predicting the chemotherapeutic response of patients with colorectal liver metastases. Abdom. Radiol. 2018, 44, 65–71. [Google Scholar] [CrossRef] [PubMed]
  34. Lubner, M.G.; Stabo, N.; Lubner, S.J.; Del Rio, A.M.; Song, C.; Halberg, R.B.; Pickhardt, P.J. CT textural analysis of hepatic metastatic colorectal cancer: Pre-treatment tumor heterogeneity correlates with pathology and clinical outcomes. Gastrointest. Radiol. 2015, 40, 2331–2337. [Google Scholar] [CrossRef] [PubMed]
  35. Simpson, A.L.; Doussot, A.; Creasy, J.M.; Adams, L.B.; Allen, P.J.; DeMatteo, R.P.; Gönen, M.; Kemeny, N.E.; Kingham, T.P.; Shia, J.; et al. Computed Tomography Image Texture: A Noninvasive Prognostic Marker of Hepatic Recurrence After Hepatectomy for Metastatic Colorectal Cancer. Ann. Surg. Oncol. 2017, 24, 2482–2490. [Google Scholar] [CrossRef] [PubMed]
  36. Ganeshan, B.; Miles, K.A.; Young, R.C.; Chatwin, C.R. Hepatic Enhancement in Colorectal Cancer: Texture Analysis Correlates with Hepatic Hemodynamics and Patient Survival. Acad. Radiol. 2007, 14, 1520–1530. [Google Scholar] [CrossRef] [PubMed]
  37. Rahmim, A.; Bak-Fredslund, K.P.; Ashrafinia, S.; Lu, L.; Schmidtlein, C.; Subramaniam, R.M.; Morsing, A.; Keiding, S.; Horsager, J.; Munk, O.L. Prognostic modeling for patients with colorectal liver metastases incorporating FDG PET radiomic features. Eur. J. Radiol. 2019, 113, 101–109. [Google Scholar] [CrossRef]
  38. Dercle, L.; Lu, L.; Schwartz, L.H.; Qian, M.; Tejpar, S.; Eggleton, P.; Zhao, B.; Piessevaux, H. Radiomics Response Signature for Identification of Metastatic Colorectal Cancer Sensitive to Therapies Targeting EGFR Pathway. JNCI J. Natl. Cancer Inst. 2020, 112, 902–912. [Google Scholar] [CrossRef]
  39. Ravanelli, M.; Agazzi, G.M.; Tononcelli, E.; Roca, E.; Cabassa, P.; Baiocchi, G.L.; Berruti, A.; Maroldi, R.; Farina, D. Texture features of colorectal liver metastases on pretreatment contrast-enhanced CT may predict response and prognosis in patients treated with bevacizumab-containing chemotherapy: A pilot study including comparison with standard chemotherapy. Radiol. Med. 2019, 124, 877–886. [Google Scholar] [CrossRef]
  40. Taghavi, M.; Staal, F.C.; Simões, R.; Hong, E.K.; Lambregts, D.M.; van der Heide, U.A.; Beets-Tan, R.G.; Maas, M. CT radiomics models are unable to predict new liver metastasis after successful thermal ablation of colorectal liver metastases. Acta Radiol. 2021, 2841851211060437, Epub ahead of print. [Google Scholar] [CrossRef]
  41. Rizzetto, F.; Calderoni, F.; De Mattia, C.; Defeudis, A.; Giannini, V.; Mazzetti, S.; Vassallo, L.; Ghezzi, S.; Sartore-Bianchi, A.; Marsoni, S.; et al. Impact of inter-reader contouring variability on textural radiomics of colorectal liver metastases. Eur. Radiol. Exp. 2020, 4, 1–12. [Google Scholar] [CrossRef]
  42. Wang, Y.; Liu, Y.; Cheng, Y.; Yi, D.; Li, P.; Song, W.; Fu, D.; Wang, X. Prognosis of colorectal cancer with liver metastasis: Value of a prognostic index. Braz. J. Med. Biol. Res. 2010, 43, 1116–1122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kuo, I.-M.; Huang, S.-F.; Chiang, J.-M.; Yeh, C.-Y.; Chan, K.-M.; Chen, J.-S.; Yu, M.-C. Clinical features prognosis in hepatectomy for colorectal cancer with centrally located liver metastasis. World J. Surg. Oncol. 2015, 13, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Scheele, J.; Stang, R.; Altendorf-Hofmann, A.; Paul, M. Resection of colorectal liver metastases. World J. Surg. 1995, 19, 59–71. [Google Scholar] [CrossRef] [PubMed]
  45. Ma, Z.-H.; Wang, Y.-P.; Zheng, W.-H.; Ma, J.; Bai, X.; Zhang, Y.; Wang, Y.-H.; Chi, D.; Fu, X.-B.; Hua, X.-D. Prognostic factors and therapeutic effects of different treatment modalities for colorectal cancer liver metastases. World J. Gastrointest. Oncol. 2020, 12, 1177–1194. [Google Scholar] [CrossRef]
Figure 1. An example of manual definition of the ROIs made using the segmentation tool of 3DSlicer on VIBE T1-W_FA10.
Figure 1. An example of manual definition of the ROIs made using the segmentation tool of 3DSlicer on VIBE T1-W_FA10.
Cancers 14 01110 g001
Figure 2. ROC curves of linear regression analysis with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by arterial phase.
Figure 2. ROC curves of linear regression analysis with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by arterial phase.
Cancers 14 01110 g002
Figure 3. ROC curves of linear regression analysis with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by portal phase.
Figure 3. ROC curves of linear regression analysis with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by portal phase.
Cancers 14 01110 g003
Figure 4. ROC curves of KNN with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by arterial phase.
Figure 4. ROC curves of KNN with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by arterial phase.
Cancers 14 01110 g004
Figure 5. ROC curves of KNN with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by portal phase.
Figure 5. ROC curves of KNN with respect to the front of tumor growth (A), tumor budding (B), tumor mucinous type (C), and the recurrence presence (D) obtained considering significant features extracted by portal phase.
Cancers 14 01110 g005
Table 1. Characteristics of the study population (81 patients).
Table 1. Characteristics of the study population (81 patients).
Patient DescriptionNumbers (%)/Range
Gender Men 53 (65.4%)
Women 28 (34.6%)
Age61 y; range: 35–82 y
Primary cancer site
Colon52 (64.2%)
Rectum29 (35.8%)
Prior Chemotherapy81 (100%)
Hepatic metastases description
Patients with single nodule52 (64.2%)
Patients with multiple nodules 29 (35.8%)/range: 2–13 metastases
Nodule size (mm)mean size 36.4 mm; range 7–58 mm
Front of tumor growth
Expansive30 (37.0%)
Infiltrative 51 (63.0%)
Tumor Budding
Absent12 (14.8%)
Low grade14 (17.3%)
High grade55 (67.9%)
Mucinous carcinoma25 (30.9%)
Recurrence
(new liver metastases)
19 (23.5%)
RAS mutation42 (51.9%)
Table 2. MR Sequence parameters.
Table 2. MR Sequence parameters.
SequenceOrientationTR/TE/FA
(ms/ms/deg.)
AT
(min)
Acquisition MatrixST/Gap (mm)FS
Trufisp T2-WCoronal4.30/2.15/800.46512 × 5124/0without
HASTE T2-W Axial1500/90/1700.36320 × 3205/0without and with (SPAIR)
HASTE T2wCoronal1500/92/1700.38320 × 3205/0without
In-Out phase T1-WAxial160/2.35/700.33256 × 1925/0without
VIBE
T1-W_FA10
Axial4.80/1.76/100.18320 × 2603/0with (SPAIR)
VIBE
T1-W_FA30
Axial4.80/1.76/300.18320 × 2603/0with (SPAIR)
Note: W = weighted, TR = repetition time, TE = echo time, FA = flip angle, AT = acquisition time, SPAIR = spectral adiabatic inversion recovery, VIBE = volumetric interpolated breath hold examination, HASTE = half-Fourier-acquired single-shot turbo spin echo.
Table 3. (Sub)datasets, variable selection criteria and predictors combinations.
Table 3. (Sub)datasets, variable selection criteria and predictors combinations.
DatasetOutcome VariablePredictorsAccuracy Threshold on Univariate Analysis
Dataset 1Front of tumor growthRadiomic metrics on lesion by VIBE_FA10≥0.75
Dataset 2Tumor budding Radiomic metrics on lesion by VIBE_FA10≥0.80
Dataset 3Mucinous TypeRadiomic metrics on lesion by VIBE_FA10≥0.80
Dataset 4Recurrence presenceRadiomic metrics on lesion by VIBE_FA10≥0.80
Dataset 5Front of tumor growthRadiomic metrics on lesion by VIBE_FA30≥0.80
Dataset 6Tumor budding Radiomic metrics on lesion by VIBE_FA30≥0.85
Dataset 7Mucinous TypeRadiomic metrics on lesion by VIBE_FA30≥0.85
Dataset 8Recurrence presenceRadiomic metrics on lesion by VIBE_FA30≥0.85
Table 4. Findings by univariate analysis with ROC performance results.
Table 4. Findings by univariate analysis with ROC performance results.
Significant Textural Features Extractedby Arterial Phase Respect to the Front of Tumor Growth by Portal Phase Respect to the Front of Tumor Growthby Arterial Phase Respect to the Tumor Budding by Portal Phase Respect to the Tumor Budding by Arterial Phase Respect to the Mucinous Type by Portal Phase respect to the Mucinous Type by Arterial Phase Respect to Recurrence by Portal Phase Respect to Recurrence
wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasiswavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasiswavelet_LHH_firstorder_Minimumwavelet_LLH_firstorder_10Percentilewavelet_HLH_glszm_LargeAreaHighGrayLevelEmphasiswavelet_LLL_glcm_ClusterTendencywavelet_HLH_ngtdm_Complexitywavelet_LLH_glcm_DifferenceEntropy
AUC0.690.800.710.800.590.700.740.74
Sensitivity0.950.840.980.960.350.380.710.71
Specificity0.510.770.520.810.991.000.950.94
PPV0.770.850.850.930.901.000.790.81
NPV0.850.740.890.860.850.860.900.90
Accuracy0.790.820.860.920.850.880.900.89
Cut-off0.120.12−41.76−37.14−0.02408.223.341.54
Table 5. Linear regression and pattern recognition analysis with significant features from the arterial phase.
Table 5. Linear regression and pattern recognition analysis with significant features from the arterial phase.
Linear Regression of Significant Features Extracted by the Arterial PhaseAUCSensitivitySpecificityPPVNPVAccuracyCut-off
respect to the front of tumor growth0.740.890.890.930.830.891.45
respect to the budding0.920.940.900.970.850.931.38
respect to the mucinous type0.930.770.990.950.940.940.37
respect to the recurrence0.810.580.970.860.870.870.43
Pattern Recognition Analysis with Significant FeaturesDatasetAUCAccuracySensitivitySpecificityTraining
Time [sec]
Model Type and Parameters
The best classifier is a KNN considering significant features extracted on arterial phase respect each of outcome (front of tumor growth, budding, mucinous type, recurrence)Training set0.970.910.910.912.34Weighted KNN; number of neighbors:10; distance metric: Euclidean; distance weight: squared inverse
Validation set0.960.890.850.91
Training set0.950.950.840.994.27
Validation set0.950.950.81
Training set0.870.880.970.568.55
Validation set0.910.910.960.73
Training set0.960.920.970.7710.38
Validation set0.930.9210.66
Table 6. Results of linear regression and pattern recognition analysis with significant features from the portal phase.
Table 6. Results of linear regression and pattern recognition analysis with significant features from the portal phase.
Linear Regression of Significant Features Extracted by The Portal PhaseAUCSensitivitySpecificityPPVNPVAccuracyCut-off
respect to the front of tumor growth0.880.800.890.920.730.831.58
respect to the budding0.820.930.670.830.860.831.50
respect to the mucinous type0.880.770.960.830.940.920.36
respect to the recurrence0.920.940.820.640.970.850.28
Pattern recognition analysis resultsDatasetAUCAccuracySensitivitySpecificityTraining
time [sec]
Model Type and parameters
The best classifier is a KNN considering significant features extracted on portal phase respect to the front of tumor growthTraining set0.960.900.910.8913.4Weighted KNN; number of neighbors:10; distance metric: Euclidean; distance weight: squared inverse
Validation set0.970.920.840.979.74
The best classifier is a decision tree considering significant features extracted on portal phase respect to the buddingTraining set0.990.910.810.96 Maximum number of splits: 100
Split criterion: Gini’s diversity index
Surrogate decision splits: Off
Hyperparameter options disabled
Validation set0.970.930.840.973.4
The best classifier is a KNN considering significant features extracted on portal phase respect to the to the mucinous typeTraining set0.890.930.81 Weighted KNN; number of neighbors:10; distance metric: Euclidean; distance weight: squared inverse
Validation set0.920.910.990.6211.8
Training set0.980.9210.62
The best classifier is a KNN considering significant features extracted on portal phase respect to the recurrenceValidation set0.940.930.990.7710.1
Table 7. Linear regression model coefficients and intercept with respective p value.
Table 7. Linear regression model coefficients and intercept with respective p value.
Linear Regression of the Textural Features Extracted by the Arterial Phase with Respect to the Front of Tumor GrowthCoefficientsp Value p Value
Intercept−1.990.31<0.000
wavelet_LHH_gldm_SmallDependenceLowGrayLevelEmphasis33.140.19
wavelet_LHH_firstorder_Minimum0.010.02
wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis−1.320.14
wavelet_LHH_glrlm_ShortRunEmphasis−3.320.14
wavelet_LLH_glszm_SmallAreaLowGrayLevelEmphasis2.110.03
wavelet_HLH_glcm_MaximumProbability19.520.00
wavelet_HHH_gldm_SmallDependenceHighGrayLevelEmphasis5.170.39
wavelet_HHH_glrlm_ShortRunHighGrayLevelEmphasis0.060.70
Linear regression of the textural features extracted by the arterial phase with respect to the tumor buddingCoefficientsp value p value
Intercept−12.520.00<0.000
original_glcm_Idn31.700.00
original_glcm_Idm42.600.00
original_glcm_Id−56.440.00
wavelet_LHH_firstorder_Minimum0.020.00
wavelet_LHH_firstorder_10Percentile−0.060.40
wavelet_LLH_glcm_MaximumProbability1.880.16
wavelet_LLH_glcm_Imc18.920.01
wavelet_LLH_firstorder_10Percentile0.000.74
wavelet_LLH_glrlm_GrayLevelNonUniformityNormalized−4.570.05
wavelet_LLH_glszm_SmallAreaLowGrayLevelEmphasis1.670.11
wavelet_HLH_firstorder_10Percentile0.440.00
Linear regression of the textural features extracted by the arterial phase with respect to the mucinous typeCoefficientsp value p value
Intercept−2.180.01<0.000
original_glszm_ZoneVariance0.000.14
original_glszm_LargeAreaEmphasis0.000.11
original_glszm_LargeAreaLowGrayLevelEmphasis0.000.01
wavelet_HLL_glcm_InverseVariance4.620.01
wavelet_HLL_glrlm_RunLengthNonUniformity0.000.01
wavelet_LHH_glszm_LargeAreaEmphasis0.000.08
wavelet_LHH_glszm_ZonePercentage0.000.01
wavelet_LHH_glszm_LargeAreaLowGrayLevelEmphasis17.350.00
wavelet_LHH_glszm_HighGrayLevelZoneEmphasis0.000.00
wavelet_LLH_glcm_InverseVariance0.000.95
wavelet_HLH_glcm_Imc10.610.64
wavelet_HLH_glszm_LargeAreaHighGrayLevelEmphasis11.350.00
wavelet_HHH_glszm_ZonePercentage0.000.00
Linear regression of the textural features extracted by the arterial phase with respect to the recurrence presenceCoefficientsp value p value
Intercept0.440.110.030
wavelet_LHL_glcm_JointAverage0.00-
wavelet_LHL_glcm_SumAverage−0.200.08
wavelet_LHL_glcm_MCC0.260.65
wavelet_LHL_glszm_SmallAreaHighGrayLevelEmphasis−0.030.42
wavelet_LHL_glszm_HighGrayLevelZoneEmphasis0.070.04
wavelet_LHL_ngtdm_Complexity−0.020.48
wavelet_LLH_firstorder_InterquartileRange0.110.20
wavelet_LLH_firstorder_RobustMeanAbsoluteDeviation−0.250.22
wavelet_LLH_ngtdm_Contrast8.370.04
wavelet_HLH_ngtdm_Complexity0.030.07
Linear regression of the textural features extracted by the portal phase with respect to the front of tumor growthCoefficientsp value p value
Intercept−5.360.09<0.000
wavelet_LHH_gldm_SmallDependenceLowGrayLevelEmphasis−11.710.34
wavelet_LHH_glrlm_ShortRunLowGrayLevelEmphasis1.470.01
wavelet_LHH_glszm_GrayLevelNonUniformityNormalized0.140.78
wavelet_LLH_firstorder_10Percentile0.000.57
wavelet_HLH_glcm_JointEnergy23.110.06
wavelet_HLH_glcm_MCC1.220.11
wavelet_HHH_glcm_MCC16.450.00
wavelet_HHH_glcm_Imc2−9.750.04
wavelet_LLL_firstorder_Uniformity−0.500.47
Linear regression of the textural features extracted by the portal phase with respect to the tumor buddingCoefficientsp value p value
Intercept29.690.06<0.000
original_glrlm_GrayLevelNonUniformityNormalized2.160.52
original_glszm_ZoneVariance0.000.02
original_glszm_SmallAreaLowGrayLevelEmphasis1.380.39
wavelet_LHH_firstorder_10Percentile0.180.00
wavelet_LHH_ngtdm_Busyness0.000.44
wavelet_LLH_firstorder_10Percentile0.020.00
wavelet_LLH_glszm_LargeAreaLowGrayLevelEmphasis0.000.23
wavelet_LLH_glszm_SmallAreaLowGrayLevelEmphasis5.370.06
wavelet_HHH_glcm_JointEnergy−111.340.09
wavelet_HHH_glcm_MCC16.230.00
wavelet_LLL_glrlm_GrayLevelNonUniformityNormalized−8.050.15
wavelet_LLL_glszm_ZoneVariance0.000.33
wavelet_LLL_glszm_LargeAreaEmphasis0.000.38
Linear regression of the textural features extracted by the portal phase with respect to the mucinous typeCoefficientsp value p value
Intercept−0.100.51<0.000
original_gldm_GrayLevelVariance−2.920.05
original_glcm_SumSquares2.400.32
original_glcm_ClusterProminence0.000.16
original_glcm_ClusterTendency0.180.82
original_firstorder_Variance0.000.81
original_glrlm_GrayLevelVariance−0.620.00
wavelet_LLL_gldm_GrayLevelVariance1.730.03
wavelet_LLL_glcm_SumSquares−0.300.35
wavelet_LLL_glcm_ClusterProminence0.000.10
wavelet_LLL_glcm_ClusterTendency−0.020.87
wavelet_LLL_firstorder_Variance0.000.10
wavelet_LLL_glszm_GrayLevelVariance0.020.00
Linear regression of the textural features extracted by the portal phase h respect to the recurrence presenceCoefficientsp value p value
Intercept−0.230.81<0.000
wavelet_LLH_gldm_GrayLevelVariance6.150.00
wavelet_LLH_glcm_JointEntropy−0.250.48
wavelet_LLH_glcm_Contrast−2.960.01
wavelet_LLH_glcm_DifferenceEntropy−4.970.05
wavelet_LLH_glcm_DifferenceVariance4.990.03
wavelet_LLH_glcm_DifferenceAverage9.930.00
wavelet_LLH_firstorder_MeanAbsoluteDeviation0.090.14
wavelet_LLH_firstorder_RootMeanSquared0.050.14
wavelet_LLH_firstorder_Variance−0.010.00
wavelet_LLH_firstorder_Mean0.040.00
wavelet_LLH_glrlm_GrayLevelVariance−1.060.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Granata, V.; Fusco, R.; De Muzio, F.; Cutolo, C.; Setola, S.V.; dell’ Aversana, F.; Ottaiano, A.; Avallone, A.; Nasti, G.; Grassi, F.; et al. Contrast MR-Based Radiomics and Machine Learning Analysis to Assess Clinical Outcomes following Liver Resection in Colorectal Liver Metastases: A Preliminary Study. Cancers 2022, 14, 1110. https://doi.org/10.3390/cancers14051110

AMA Style

Granata V, Fusco R, De Muzio F, Cutolo C, Setola SV, dell’ Aversana F, Ottaiano A, Avallone A, Nasti G, Grassi F, et al. Contrast MR-Based Radiomics and Machine Learning Analysis to Assess Clinical Outcomes following Liver Resection in Colorectal Liver Metastases: A Preliminary Study. Cancers. 2022; 14(5):1110. https://doi.org/10.3390/cancers14051110

Chicago/Turabian Style

Granata, Vincenza, Roberta Fusco, Federica De Muzio, Carmen Cutolo, Sergio Venanzio Setola, Federica dell’ Aversana, Alessandro Ottaiano, Antonio Avallone, Guglielmo Nasti, Francesca Grassi, and et al. 2022. "Contrast MR-Based Radiomics and Machine Learning Analysis to Assess Clinical Outcomes following Liver Resection in Colorectal Liver Metastases: A Preliminary Study" Cancers 14, no. 5: 1110. https://doi.org/10.3390/cancers14051110

APA Style

Granata, V., Fusco, R., De Muzio, F., Cutolo, C., Setola, S. V., dell’ Aversana, F., Ottaiano, A., Avallone, A., Nasti, G., Grassi, F., Pilone, V., Miele, V., Brunese, L., Izzo, F., & Petrillo, A. (2022). Contrast MR-Based Radiomics and Machine Learning Analysis to Assess Clinical Outcomes following Liver Resection in Colorectal Liver Metastases: A Preliminary Study. Cancers, 14(5), 1110. https://doi.org/10.3390/cancers14051110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop