Radiomic Feature-Based Predictive Model for Microvascular Invasion in Patients With Hepatocellular Carcinoma

Objective: This study aimed to build and evaluate a radiomics feature-based model for the preoperative prediction of microvascular invasion (MVI) in patients with hepatocellular carcinoma.

Methods: A total of 145 patients were retrospectively included in the study pool, and the patients were divided randomly into two independent cohorts with a ratio of 7:3 (training cohort: n = 101, validation cohort: n = 44). For a pilot study of this predictive model another 18 patients were recruited into this study. A total of 1,231 computed tomography (CT) image features of the liver parenchyma without tumors were extracted from portal-phase CT images. A least absolute shrinkage and selection operator (LASSO) logistic regression was applied to build a radiomics score (Rad-score) model. Afterwards, a nomogram, including Rad-score as well as other clinicopathological risk factors, was established with a multivariate logistic regression model. The discrimination efficacy, calibration efficacy, and clinical utility value of the nomogram were evaluated.

Results: The Rad-score scoring model could predict MVI with the area under the curve (AUC) of 0.637 (95% CI, 0.516–0.758) in the training cohort as well as of 0.583 (95% CI, 0.395–0.770) in the validation cohort; however, the aforementioned discriminative approach could not completely outperform those existing predictors (alpha fetoprotein, neutrophilic granulocyte, and preoperative hemoglobin). The individual predictive nomogram which included the Rad-score, alpha fetoprotein, neutrophilic granulocyte, and preoperative hemoglobin showed a better discrimination efficacy with AUC of 0.865 (95% CI, 0.786–0.944), which was higher than the conventional methods’ AUCs (nomogram vs Rad-score, alpha fetoprotein, neutrophilic granulocyte, and preoperative hemoglobin at P < 0.001, P = 0.025, P < 0.001, and P = 0.001, respectively). When applied to the validation cohort, the nomogram discrimination efficacy was still outbalanced those above mentioned three remaining methods (AUC: 0.705; 95% CI, 0.537–0.874). The calibration curves of this proposed method showed a satisfying consistency in both cohorts. A prospective pilot analysis showed that the nomogram could predict MVI with an AUC of 0.844 (95% CI, 0.628–1.000).

Conclusions: The radiomics feature-based predictive model improved the preoperative prediction of MVI in HCC patients significantly. It could be a potentially valuable clinical utility.

Introduction

Hepatocellular carcinoma (HCC) is considered as the sixth most high incidence cancer and HCC is also as the world’s third leading cause of cancer deaths. Among the primary cancer of the liver, HCC is the predominant pathological type and has became a significant public health-care concern (1, 2). Among the current treatment strategies for HCC, the optimal and most efficient treatment is radical surgical resection. In the past several decades, along with the development of modern medical imaging technology and improved surgical skills, the success ratio of hepatectomy for HCC has increased continuously (3). In spite of this, the postoperative five-year recurrence ratio of HCC is approximately 70%, which is unsatisfactory and indicates a significant influence on overall survival for HCC patients (4). Microvascular invasion (MVI) is a crucial predictive factor for recurrence among HCC patients who underwent hepatoectomy for tumor resection or received liver transplant surgery, and most patients with early recurrence were pathologically verified as MVI positive (5, 6). Therefore, prediction for MVI before surgery has significant clinical value for decision making, postoperative adjuvant therapy and comprehensive prognostic evaluation of patients with HCC.

Tumor macrovascular invasion of hepatic vein and/or portal vein in HCC patients can be identified through contrast-enhanced computed tomography (CECT) or magnetic resonance imaging (MRI) examination before surgery (7). Conversely, microvascular invasion (MVI) cannot be detected directly by means of preoperative radiologic examination. In addition, the effectiveness of preoperative pathological biopsy for MVI is unsatisfactory to some extent. Hence, the gold standard to verify the occurrence of MVI is to evaluate the surgical specimen slices through preoperative pathological examination (i.e., H&E staining, etc.) which limits its application as a predictive factor for comprehensive prognosis assessment (8). Imaging markers are emerging as preoperative predictive factors for tumor heterogeneity that have the potential to predict MVI in a high-throughput, high-precision, and noninvasive manner. Current published studies have indicated that the correlation between the MVI and imaging features could be quantitively analyzed for further research on the occurrence and progression of MVI. Imaging features could be extracted from the following medical imaging modalities: contrast-enhanced CT, gadolinium-ethoxybenzyl-diethylenetriamine penta-acetic acid (Gd-EOB-DTPA; EOB Primovist^®)-enhanced magnetic resonance imaging (MRI) and ¹⁸F-fluorodeoxyglucose positron emission tomography (PET). In addition, those imaging features incorporated laborator testing, quantitative or qualitative results to establish predictive models to classify whether patients with HCC are MVI positive or negative, which is of great value for surgeons to evaluate liver transplant recipients or HCC therapy scheme decision making (9–11). In addition, radiomics is an emerging methodology for medical imaging analysis that applies machine learning algorithms and statistical analysis software to obtain high-throughput radiological features for prognosis prediction or tumor-related overall survival analysis. In the past, our research group’s study results indicated that radiomics-based predictive models have the potential to predict the occurrence of posthepatectomy liver failure and the occurrence of postpancreaticoduodenectomy pancreatic fistula among a study cohort with satisfactory predictive performance (12, 13). Recent studies showing predictive radiomics models based on contrast-enhanced CT (CECT) as an easy-to-use, efficient, and noninvasive prediction method for MVI have provided reliable instruction for the definite diagnosis of whether HCC patients had MVI preoperatively (14, 15). Precise prediction of MVI before surgery has significant clinical implications for issues such as efficient stratification of treatment strategies, optimization of treatment schemes, and establishment of comprehensive surgery planning for HCC patients with MVI positivity or negativity. Due to the specific biological characteristics of MVI, the contrast agent in CECT likely enters the vessels of tumors and then diffuses into the microvessels over the range of one centimeter around the primary lesion. An enlarged surgical margin (usually over one centimeter) was verified to decrease the postoperative tumor recurrence ratio, indicating the presence of MVI, among patients with HCC (16, 17). For patients who have a tumor diameter less than three centimeters and are highly suspected to be MVI positive, anatomical hepatectomy is necessary. For non-anatomical hepatectomy, the surgical margin should be over one centimeter from the primary lesion (18, 19). Compared with surgical resection, radiofrequency ablation is not a rational option for small (diameter under three centimeters) HCC patients combined with MVI because of high postoperative recurrence ratio (20).

The radiomic predictive model incorporates clinical predicative factors and imaging characteristic features into a reliable and efficient MVI risk-stratification scoring system, which helps clinical practitioners provide optimal treatment schemes for HCC patients.

In the present study, we aimed at build as well as evaluate a radiomics nomogram based on CT imaging textural features incorporated with clinically related risk factors for MVI occurrence prediction among HCC patients before surgery.

Materials and Methods

Patients

The data originated from 146 patients who received liver resection at Zhujiang Hospital, Southern Medical University, between January 2012 and December 2018, were retrospectively collected and analyzed. The patients with HCC were carefully selected, and the inclusion standards as indicated below: (1) Be at least 18 years old, both male and female were included; (2) contrast-enhanced CT scans were performed within one week before surgery; and (3) all surgical specimens were confirmed to be HCC by MVI evaluation and histopathological examination. The exclusion standards were as follows: (1) previous hepatectomy; (2) HCC lesion boundaries not clearly visible on CT images; (3) preoperative treatment, including transhepatic arterial chemotherapy and embolization, portal vein embolization, radiofrequency ablation therapy, targeted cancer therapy, etc.; and (4) incomplete clinical data or pathology data. Patients were randomly allocated into two independent cohorts as training cohort and validation cohort according to a ratio of 7:3. An additional 18 patients between January and December 2019 were enrolled in a prospective pilot analysis, and all patients signed the informed consent. The Ethics Review Board of Zhujiang Hospital, Southern Medical University, approved this study and supervised the procedures.

Histopathological Analysis

The final diagnostic results, including MVI grade, were based on pathologic reports of surgical specimens. On the basis of the Clinical Practice Guidelines of Chinese Society of Pathology (21), the definite diagnosis of MVI was based on the pathological standard as follows: MVI was graded according to the number of cells seen in the endothelial vascular lumen by microscopic examination, and then MVI was categorized as three additional subgrades, M0 indicating no MVI; M1 for fewer than five MVI occurrence in adjoining liver parenchyma and the distance from tumors was under one centimeter (which is also regarded as the low-risk group); and M2 for more than five MVI in adjoining liver parenchyma and the distance from tumors was greater than one centimeter (regarded as the high-risk group). All specimen slices were reviewed and analyzed by two pathologists who had at least 10 to 15 years of hepatic pathology experience. Aforementioned two specialists were blinded to the related clinicopathological as well medical imaging examination information of the patients who were included in the study cohort. When there was a disagreement, a consistency review was performed.

CT Data Collection

All CT examinations were finished in our hospital spiral CT (128 -slice, Siemens Medical system, Erlangen, Germany). To ensure the consistency of the image data, all the scanning parameters were set as follows: 120 kV, 250–300 mAs, and a matrix of 512 × 512. The slice thickness and interval were set as one mm and zero mm, respectively, and image reconstruction was set to 1 mm. To obtain contrast-enhanced image data, 90 ml iodine contrast agent (iodine concentration: 350 mg/ml) was administered via a vein by a contrast injector at a speed of 3 milliliter per second, and dynamic contrast-enhanced imaging data acquisition was performed at fixed time points: for the arterial phase, acquisition occurred at approximately 30–33 s after administration; for the portal vein phase, it was 67–70 s, and for the equilibrium phase, 177–180 s. After the entire acquisition procedure was completed, we selected the portal venous phase raw data for texture analysis. All the medical images in this study were collected from the Pictures Archiving and Communication System at our hospital radiology department.

Radiomics Feature Extraction

The holistic radiomics modeling procedure was shown in Figure 1. We extracted a total of 1,231 CT image features from the liver parenchyma without tumors. The extracted image features were categorized into four types as follows: diagnostic intensity for liver parenchyma (N=13), original hepatic textural feature (N=100), tumor size and shape-based feature (N=430), and wavelet (N=688). The mathematical descriptions of each corresponding feature above can be found in the website of the open-source application PyRadiomics (https://pyradiomics.readthedocs.io/en/latest/).

FIGURE 1

Figure 1 Flow chart of the study. (A) The portal vein contrast-enhanced computed tomography (CECT) image was used as input, and the liver was segmented by homemade segmentation software (22) and manually corrected. The tumor ROI was contoured by a radiologist with more than 10 years of experience. (B) The liver parenchyma without tumor was set as the ROI image, and features containing as shape related features, texture related features and wavelet filtered features which were extracted by PyRadiomics software (23). (C) The extracted features were reduced via least absolute shrinkage and selection operator (LASSO) and correlation matrix selection. (D) The radiomics features and other clinicopathological variables were applied to establish a nomogram for predict microvascular invasion, and ROI curve analysis was conducted. (E) The nomogram model was evaluated with calibration curve and decision curve analysis.

Data Analysis

One hundred forty-five patients were included in the study pool, therefore, those patients were randomly allocated into two independent cohorts according to a ratio of 7:3 (training cohort: n = 101, validation cohort: n = 44). An additional 18 patients were also enrolled into this study for a prospective pilot analysis which was aimed to evaluate its underlying clinical practice potency.

A histologic diagnosis of tumor with microvascular invasion or tumor without microvascular invasion was entered as a dependent variable in a least absolute shrinkage and selection operator (LASSO) logistic regression algorithm. As to the LASSO logistic regression, we applied 10-fold cross validation method to identified the optimal regularization parameter lambda, which reduces the mean cross-validation error at maximum extent. those image features who were considered as non-zero coefficients within the LASSO regression model were regarded as most predictive radiomics features and were chosen for calculating the radiomics score. The radiomics score was calculated through a formula which in brief was a semplice linear combination of the selected radiomics based features multiplied by their corresponding coefficients.

The radiomics score as well as other clinicopathological risk variables, including age, sex, tumor sizes, viral hepatitis etiology, drinking history, alpha fetoprotein content, tumor numbers, neutrophil counts, lymphocyte counts, neutrophil-lymphocyte ratio, preoperative alanine transaminase, preoperative aspartate aminotransferase, preoperative total bilirubin, preoperative hemoglobin, preoperative albumin, preoperative blood platelet counts, preoperative prothrombin time, and histopathological differentiation degree, were used in a multivariate logistic regression model as predictive factors for tumors with microvascular invasion or tumors without microvascular invasion among training cohort. The optimal model selection was performed applied a backward stepwise selection procedure through utilizing likelihood ratio testing with Akaike information criterion (AIC). Nomogram was constructed based on those significant features in the multivariate logistic regression model among which radiomics score as well as clinicopathological risk factors were included.

The nomogram’s predictive accuracy was evaluated by three methods. First, discrimination efficacy was assessed through the result of the area under the area under curve (AUC) with values ranging from 0.5 (which means mentioned model possess no discrimination potent) to 1 (which means mentioned model possess a perfect discrimination potent). Second, by plotting the predicted probability of the nomogram relative to the ratio of the observed events on a series of equidistant values within the range of predicted probabilities, the power of calibration can be evaluated visually. In addition, decision curve analysis (DCA) was carried out to tell physicians the interval of threshold probabilities to judge whether the prediction model would be of clinical benefit or not. If the typical threshold probabilities of tumors with microvascular invasion or tumors without microvascular invasion at which surgeons would opt for comprehensive treatment lie within a certain benefit range, the model is of clinical value.

Statistical Analysis

For statistical analysis in this study an open source software R (version 3.6.1) was applied. (https://www.r-project.org). consider the various type of data, we employed different statistical test as below: for continuous variables, the Mann-Whitney U test were used for the comparison. Correspondingly, as for categorical variables, the Chi-square test was used for the comparison. in order to conduct LASSO logistic regression analysis we used “glmnet” package, “pROC” package was applied for calculation and comparisons of AUC and ROC plotting of each model as well as for nomogram training and calibration “rms” package was employed, in addition for decision curve analysis we chosen “rmda” package to conduct it. A two-sided P value < 0.05 was considered as indicative of statistical significance. For the final independent predictors included in this prediction model selection, we applied both univariant logistic regression and multivariant logistic regression and the factors whose P value < 0.05 were considered as indicative of final independent predictors.H

Results

Patient Demographics and Clinicopathological Characteristics

The demographics and clinicopathological characteristics of the study participants are displayed in Table 1. It is clearly that there were no differences between the training cohort and the validation cohort. The median ages of these two cohorts were 50.0 and 47.5 years, respectively. As for patient gender in both cohort male was ponderance proportion (81.2% and 84.1%, respectively). In addition, for positive occurrence ratio there were 30 (29.7%) and 16 (36.4%) patients who experienced tumors with microvascular invasion among training cohort as well validation cohort, respectively. For patients among pilot study, those data are displayed in Table 2.

TABLE 1

Table 1 Demographic and clinicopathologic characteristics of the study participants.

TABLE 2

Table 2 Demographic and clinicopathologic characteristics of the pilot participants.

Development of the Rad-Score

A total of 1,231 features were extracted from the image data for each patient. These 1,231 features were reduced to 2 possible indicators based on the 101 patients among training cohort. The Rad-score formula constructed through a simple linear combination of aforementioned 2 indicators multiplied by LASSO coefficients as follows: Rad-score=55186.42 + 37645.32*log-sigma-3-0-mm-3D_gldm_SmallDependenceLowGrayLevelEmphasis - 46900.48 * log-sigma-5-0-mm-3D_glcm_Idmn (Figure 2). Patients who had tumors without MVI obtained relatively lower Rad-scores generally than patients who suffered from tumors combined with microvascular invasion. This discrepancy was detected between on Rad-scores (median (interquartile range)) among training cohort (−2.567 (-6.731~1.592) vs 0.164 (-5.160~4.068), respectively, P=0.013) also among the validation cohort (-2.477 (-7.465~0.864) vs 0.123 (-6.434~3.899), respectively, P=0.369).

FIGURE 2

Figure 2 The least absolute shrinkage and selection operator (LASSO) logistic regression algorithm was employed to determine the most significant radiomics features in terms of microvascular invasion (MVI) prediction as well calculated for the Rad-score. (A) a ten-fold cross validation method was employed to select the optimal lambda parameter contained in the LASSO algorithm that minimized the mean squared error. The mean-squared errors between upper bound and lower bound standard deviations were plotted according to the lambda changing sequence. The vertical dotted lines in the graph stands for the two selected lambdas. The lambda on the left corner of graph provided the minimum mean cross-validation error, meanwhile, the lambda on the right corner of graph provided the most optimal regularized model which the error within one minimal standard error, therefore, the optimal lambda value 0.0063, whose log(lambda)= -5.0637, was employed to feature selection. (B) A coefficient profile plot of the LASSO model was produced. Every single curve represents the trajectory changing of each independent predictor. The greater lambda is, the more the coefficients are reduced. The vertical dotted line indicates the optimal lambda, which results in two coefficients.

Development, Validation, and Assessment of the Nomogram

Both univariate and multivariate logistic regression analysis were conducted combining the Rad-score as well as clinical stage among training cohort. The backward stepwise selection process selected Rad-score, alpha fetoprotein content, neutrophil counts, and preoperative hemoglobin as the final independent predictors. A prediction model was constructed and is shown as a nomogram in Figure 3. The odds ratios of the prediction model are presented in Table 3. Among the training cohort, the AUC of the nomogram was 0.865 (95% CI 0.786–0.944). among validation cohort, the AUC was 0.705 (95% CI: 0.537–0.874). The AUC value of the ROC curve of the nomogram was significantly higher than that of each independent predictor (nomogram vs Rad-score, alpha fetoprotein, neutrophilic granulocyte, and preoperative hemoglobin at P < 0.001, P = 0.025, P < 0.001, and P = 0.001, respectively). The particular results of the discriminatory efficiency are shown in Figures 3C, D). This indicates that the Rad-score incorporated with another three clinicopathological factors possessed a greater diagnostic efficiency than any single independent predictors. The calibration curve analysis of this purposed predictive model presents good agreement among our research for either the training cohort or validation cohort (Figures 4 and 5).

FIGURE 3

Figure 3 Nomogram for the prediction of the tumor with microvascular invasion occurring or not and its discrimination performance. (A) data among the training cohort were employed to establish the purposed monogram which the Rad-score, alpha fetoprotein, neutrophilic granulocyte and preoperative hemoglobin were incorporated. Instructions for reading the nomogram: Find the Rad-score in the Rad-score axis; a straight line was drew to the point axis and the intersection point stands for actual points of the predicted probability which a patient with microvascular invasion (MVI) occurrence or not according to the own Rad-score; repeat above mentioned procedures for another predictors each straight line was drew to the points axis; calculated the total points obtained from each predictor and located it on the total points axis; draw a line perpendicular to the predict risk axis to determine the patient’s risk (24). (B, C) display the result of ROC curves analysis for the nomogram and each individual predictor in predicting whether tumors with microvascular invasion occurred among the training and validation cohorts, respectively.

TABLE 3

Table 3 Risk factors for the diagnosis of tumors with microvascular invasion.

FIGURE 4

Figure 4 The result of calibration curve analysis as well decision curve analysis for the purposed nomogram. (A, B) show the calibration curves analysis for the nomogram in predicting whether tumors with microvascular invasion occurrence among training cohort as well as validation cohort, respectively. the dashed line stands for perfect prediction, the dotted line represents apparent estimates of predicted vs. observed values, meanwhile the solid line on behalf of bias correcting estimates employing 1,000 bootstrap sampling. (C, D) show the result of decision curve analysis for the nomogram in predicting whether tumors with microvascular invasion occurrence among training cohort as well as validation cohort, respectively. The result of decision curve analysis indicated that assume the threshold probability that is located in an interval of 0.2 to 0.86, microvascular invasion (MVI) prediction using the nomogram can give more net benefit than by treating either no or all patients in all cohorts.

FIGURE 5

Figure 5 The area under the curve (AUC) (A), calibration curve (B), and decision curve analysis (DCA) (C) of the nomogram to prospectively analyze the pilot cohort. This purposed radiomics feature-based nomogram was able to predict tumors combined microvascular invasion (MVI) occurrence with an AUC of 0.844 (95% confidence interval:0.628–1.000) (A) the result calibration curve showed that the radiomics nomogram was with the ideal line (B) and as to result of the DCA represented that at the situation where the threshold probability was within a range from 0.2 to 1.0 (C), this purposed nomogram obtained greater net benefit than either “treat none” or “treat all” scheme, suggesting that this purposed nomogram is still a good clinical utility.

The DCA curve showed that as threshold probability was within a range from 0.1 to 1.0, this nomogram could obtain a greater net benefit than either “treat none” or “treat all” scheme. The nomogram also received higher net benefit than the Rad-score alone across the reasonable threshold probabilities range.

Discussion

In this study, we established as well assessed radiomics model based on CT imaging textural features and clinically related risk factors for accurate prediction of MVI of HCC before surgical resection. The efficacy of this predictive model was average-to-good but did not totally outperform other current predictive models or scoring systems. Therefore, we developed and evaluated a predictive nomogram based on radiological characteristics to predict the incidence of MVI. Our predictive nomogram included factors such as AFP, NLR, routine blood tests, liver function, and pathological examination results and a combination of potential MVI predictive imaging features. Our results were consistent with previous studies; higher serum AFP level, NLR or poorly differential tumor cell type had a positive correlation with the incidence of MVI among patients with HCC. Through the results of decision curve analysis, independent predictive risk factors such as AFP, NLR, or radiomics features had beneficial effects on improving the predictive power of MVI, which was consistent with current studies; however, the diagnostic efficacy of the model did not outperform other existing predictive models or scoring systems. When we combined the above independent predictive risk factors into a comprehensive diagnostic nomogram, the diagnostic efficacy exceeded the predictive power of any single predictive factor or a combination of predictive factors, and its performance showed satisfactory consistency in both the training cohort and validation cohort. The aforementioned nomogram could improve the individual predictive accuracy of MVI for patients with HCC preoperatively, and the results of decision curve analysis showed that this practical nomogram had the potential to be applied in daily clinical practice.

The diagnosis of HCC depends on radiological examination, such as CT imaging or MRI, to a large extent, which makes it preferable to perform quantitative or qualitative imaging analysis based on the radiomics method. The radiomics method can extract a large number of specific imaging features from radiological images, such as CT images, magnetic resonance images or PET-CT images, in a high-throughput and automatic manner and then utilize those radiomics imaging features with clinicopathological data to establish statistical models for judgment of tumor characteristics or optimal treatment choice and prognosis for individual patients (25). Therefore, radiomics can connect medical imaging and personalized medicine. However, due to the lack of a comprehensive standardized assessment of clinical significance and scientific integrity in previously radiomics studies, it is difficult to promote the application of radiomics image analysis methods among different institutions (26). In current clinical practice, for the reason that there are no known single and highly reliable predictive factors for the incidence of MVI, combining multiple MVI status-related factors with radiomics signatures to build multivariate radiomic models is a feasible option. Previous published studies have utilized multiple imaging modal images (e.g., CT, MR, ultrasound) with radiomics algorithms to predict the incidence of MVI preoperatively and then achieve the ultimate goal of assisting in choosing an optimal surgical treatment or predicting the early recurrence of HCC to provide a reasonable treatment scheme for those patients (27–29). With the rapid development and improvement of radiomics technology, we are able to identify high risky HCC patients with MVI with higher accuracy and better efficacy before surgery which have great significance in terms such as select appropriate surgical therapy strategy, give assistance to confirm rational resection scope to achieve radical hepatoectomy for HCC patients rather than barely refer to tumor location, size or two-dimensional spatial relationship between intrahepatic vasculatures as conventional method did consider of the high postoperative recurrence ratio, for patients whose tumor size lower than three millimeters also susceptible to MVI positive anatomical hepatoectomy are recommended as first choice than radiofrequency ablation, for patients susceptible to MVI positive when select to perform non-anatomical hepatoectomy, the distance from resection margin tumor should beyond one millimeter to decrease the risk of postoperative recurrence (20). Accurately prediction of MVI occurrence has contributed to achieve radical tumor resection for patients with HCC in addition prolong the disease-free survival to the maximum extent. MVI as a significant factor which influence long term survival rate, radiomics features based prediction model applied to this crucial issue will provide more effective guidance for subsequent clinical researches and therapy scheme making.

Related studies have shown that approximately 85% of MVI is located within one millimeter (which is called the peritumoral area) of the margin of tumors (30). In this research, we extracted both peritumoral and intratumoral CT radiomics features to establish a nomogram model. Zhao et al. (31) established a scoring system using predictive factors such as intratumoral arteries, nonnodular HCC type also the absence of a tumor vessel based on CECT images to predict MVI preoperatively, take no account of tumor size or diameter. The AUCs of aforementioned scoring system were 0.872 and 0.856 among training cohort as well validation cohort, respectively. Lei et al. (32) built a nomogram scoring model with risk factors such as tumor diameter, number, vessel condition, serum AFP content, platelet level, HBV-DNA loading in addition with representative dynamical magnetic resonance imaging features for predicting MVI occurrence among patients suffered from HBV-related HCC under the Milan criteria, which shows satisfactory predictive efficiency whose AUC was 0.81 among training cohort also 0.80 among validation cohort. Similarly, Renzulli et al. (9) applied integration of three imaging features of predictive significance, gaining an AUC of 0.90. In terms of clinicopathological risk factors among our study, AFP, neutrophilic granulocyte, and preoperative hemoglobin content were of greater importance than the Rad-score in the regression model. When applied in the validation cohort, the discernibility ability of proposed model outperformed three conventional method mentioned in this article. For calibration curve analysis, the result in both cohorts are satisfying, furthermore, the result of decision curve analysis showed this proposed model possessed clinical application potential. The result of pilot study demonstrated that the accuracy of this model is preferably among study data set (whose AUC was 0.844, 95% confidential interval was 0.628–1.000).

Limitations

There are certainly some underlying limitations to this study. First, the analysis results suffer from inherent selection biases without a large enough number of included patients. Second, due to the study cohort selection from a single center, which could not offer sufficient diversity of characteristics to be representative of the general population, the lack of external validation from other centers was a shortcoming of this study. Moreover, heterogeneity of CT images will give rise to bias during the procedures as image features extraction and segmentation for region of interest. In the follow-up study, a larger study cohort and external validation will test and improve the efficacy and feasibility of this model for patient management and clinical decision-making for MVI prediction in HCC. Therefore, unified standards in addition with elaborated and rigorous assessment for quality of study images should be established in order to acquire more reliable study data.

Conclusions

We have built and evaluated a novel radiomics feature-based nomogram to predict the incidence of MVI among patients suffered from HCC. the purposed nomogram has provided one novel method for risk identification of MVI incidence and of radiomics features as a supplementary resource to clinicopathological data and conventional CT images that could lower the medical costs and improve the availability of clinical decision-making and treatment scheme guidance.

Data Availability Statement

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

Ethics Statement

The studies involving human participants were reviewed and approved by The Ethics Review Board of Zhujiang Hospital, Southern Medical University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

CF and FJ conceived of the presented idea. PZ collected the data. MH, XM, BH, and FJ analyzed the data. MH drafted the manuscript. All authors reviewed the manuscript and FJ made corrections to the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Key R&D Program, China (Nos. 2019YFC0118100 and 2016YFC0106500), the Major Instrument Project of National Natural Science Fund of China (No. 81627805), the Shenzhen Key Basic Science Program (JCYJ20170413162213765 and JCYJ20180507182437217), the NSFC-Shenzhen Union Program (U1613221), the Youth Program of Natural Science Foundation of Anhui Province (2008085QH418), the Key-Area Research and Development Program of Guangdong Province (2020B010165004), and the National Key Research and Development Program (2017YFC0110903).

Conflict of Interest

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

References

1. Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet (2018) 391(10127):1301–14. doi: 10.1016/S0140-6736(18)30010-2

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Miller KD, Nogueira L, Mariotto AB, Rowland JH, Yabroff KR, Alfano CM, et al. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin (2019) 69(5):363–85. doi: 10.3322/caac.21565

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fang C, Zhang P, Qi X. Digital and intelligent liver surgery in the new era: Prospects and dilemmas. Ebiomedicine (2019) 41:693–701. doi: 10.1016/j.ebiom.2019.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Vilarinho S, Calvisi DF. New advances in precision medicine for hepatocellular carcinoma recurrence prediction and treatment. Hepatology (2014) 60(6):1812–4. doi: 10.1002/hep.27311

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lim KC, Chow PK, Allen JC, Chia GS, Lim M, Cheow PC, et al. Microvascular invasion is a better predictor of tumor recurrence and overall survival following surgical resection for hepatocellular carcinoma compared to the Milan criteria. Ann Surg (2011) 254(1):108–13. doi: 10.1097/SLA.0b013e31821ad884

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Iguchi T, Shirabe K, Aishima S, Wang H, Fujita N, Ninomiya M, et al. New pathologic stratification of microvascular invasion in hepatocellular carcinoma: Predicting prognosis after living-donor liver transplantation. Transplantation (2015) 99(6):1236–42. doi: 10.1097/TP.0000000000000489

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Costentin CE, Ferrone CR, Arellano RS, Ganguli S, Hong TS, Zhu AX. Hepatocellular carcinoma with macrovascular invasion: Defining the optimal treatment strategy. Liver Cancer (2017) 6(4):360–74. doi: 10.1159/000481315

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sumie S, Nakashima O, Okuda K, Kuromatsu R, Kawaguchi A, Nakano M, et al. The significance of classifying microvascular invasion in patients with hepatocellular carcinoma. Ann Surg Oncol (2014) 21(3):1002–9. doi: 10.1245/s10434-013-3376-9

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Renzulli M, Brocchi S, Cucchetti A, Mazzotti F, Golfieri R. Can current preoperative imaging be used to detect microvascular invasion of hepatocellular carcinoma? Radiology (2015) 279(2):150998. doi: 10.1148/radiol.2015150998

CrossRef Full Text | Google Scholar

10. Lee S, Kim SH, Ji EL, Dong HS, Park CK. Preoperative Gadoxetic acid-enhanced MRI for predicting microvascular invasion in patients with single hepatocellular carcinoma. J Hepatol (2017) 67(3):526–34. doi: 10.1016/j.jhep.2017.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kornberg A, Freesmeyer M, Rthel EB, Jandt K, Katenkamp K, Steenbeck J, et al. 18F-FDG-Uptake of hepatocellular carcinoma on PET predicts microvascular tumor invasion in liver transplant patients. Am J Transplant (2009) 9(3):592–600. doi: 10.1111/j.1600-6143.2008.02516.x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cai W, He B, Hu M, Zhang W, Xiao D, Yu H, et al. A radiomics-based nomogram for the preoperative prediction of posthepatectomy liver failure in patients with hepatocellular carcinoma. Surg Oncol (2019) 28:78–85. doi: 10.1016/j.suronc.2018.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhang W, Cai W, He B, Xiang N, Fang C. Jia F. A radiomics-based formula for the preoperative prediction of postoperative pancreatic fistula in patients with pancreaticoduodenectomy. Cancer Manag Res (2018) 10:6469–78. doi: 10.2147/CMAR.S185865

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Jie P, Jing Z, Qifan Z, Yikai X, Jie Z, Li L. A radiomics nomogram for preoperative prediction of microvascular invasion risk in hepatitis B virus-related hepatocellular carcinoma. Diagn Interventional Radiol (2018), 121–7. doi: 10.5152/dir.2018.17467

CrossRef Full Text | Google Scholar

15. 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

16. Orcutt ST, Anaya DA. Liver resection and surgical strategies for management of primary liver cancer. Cancer Control J Moffitt Cancer Center (2018) 25(1):544014718. doi: 10.1177/1073274817744621

CrossRef Full Text | Google Scholar

17. Nakashima Y, Nakashima O, Tanaka M, Okuda K, Nakashima M, Kojiro M. Portal vein invasion and intrahepatic micrometastasis in small hepatocellular carcinoma by gross type. Hepatol Res (2003) 26(2):142–7. doi: 10.1016/s1386-6346(03)00007-x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lafaro K, Grandhi MS, Herman JM, Pawlik TM. The importance of surgical margins in primary malignancies of the liver. J Surg Oncol (2016) 113(3):296–303. doi: 10.1002/jso.24123

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Du M, Chen L, Zhao J, Tian F, Zeng H, Tan Y, et al. Microvascular invasion (MVI) is a poorer prognostic predictor for small hepatocellular carcinoma. BMC Cancer (2014) 14(1):38. doi: 10.1186/1471-2407-14-38

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Lee S, Kang TW, Song KD, Lee MW, Rhim H, Lim HK, et al. Effect of microvascular invasion risk on early recurrence of hepatocellular carcinoma after surgery and radiofrequency ablation. Ann Surg (2019). doi: 10.1097/SLA.0000000000003268

CrossRef Full Text | Google Scholar

21. Wen-Ming CHBJ, Guideline C. Practice guidelines for the pathological diagnosis of primary liver cancer: 2015 update. World J Gastroentero (2016) 22(42):9279–87. doi: 10.3748/wjg.v22.i42.9279

CrossRef Full Text | Google Scholar

22. He B, Huang C, Sharp G, Zhou S, Hu Q, Fang C, et al. Fast Automatic 3D liver segmentation based on a three-level AdaBoost-guided active shape model. Med Phys (2016) 43(5):2421–34. doi: 10.1118/1.4946817

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational radiomics system to decode the radiographic phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Motzer RJ, Bukowski RM, Figlin RA, Hutson TE, Dror Michaelson M, Kim ST, et al. Prognostic nomogram for Sunitinib in patients with metastatic renal cell carcinoma. Cancer (2008) 113(7):1552–8. doi: 10.1002/cncr.23776

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Lewis S, Hectors S, Taouli B. Radiomics of hepatocellular carcinoma. Abdom Radiol (2020). doi: 10.1007/s00261-019-02378-5

CrossRef Full Text | Google Scholar

26. 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

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

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Chen L, Li W, Xian M, Zheng X, Lin Y, Liu B, et al. Preoperative prediction of tumour deposits in rectal cancer by an artificial neural network–based US radiomics model. Eur Radiol (2020) 30(4):1969–79. doi: 10.1007/s00330-019-06558-1

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ji G, Zhu F, Xu Q, Wang K, Wu M, Tang W, et al. Radiomic features at contrast-enhanced CT predict recurrence in early stage hepatocellular carcinoma: A multi-institutional study. Radiology (2020) 294(3):568–79. doi: 10.1148/radiol.2020191470

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Roayaie S, Blume IN, Thung SN, Guido M, Fiel MI, Hiotis S, et al. A system of classifying microvascular invasion to predict outcome after resection in patients with hepatocellular carcinoma. Gastroenterology (2009) 137(3):850–5. doi: 10.1053/j.gastro.2009.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zhao H, Hua Y, Dai T, He J, Tang M, Fu X, et al. Development and validation of a novel predictive scoring model for microvascular invasion in patients with hepatocellular carcinoma. Eur J Radiol (2017) 88(Complete):32–40. doi: 10.1016/j.ejrad.2016.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Lei Z, Li J, Dong W, Yong X, 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