ORIGINAL RESEARCH article

Front. Oncol., 05 November 2020

Sec. Cancer Imaging and Image-directed Interventions

Volume 10 - 2020 | https://doi.org/10.3389/fonc.2020.574228

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

  • 1. The First Department of Hepatobiliary Surgery, Zhujiang Hospital, Southern Medical University, Guangdong Provincial Clinical and Engineering Center of Digital Medicine, Guangzhou, China

  • 2. Research Laboratory for Medical Imaging and Digital Surgery, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China

Abstract

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 18F-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 (911). 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

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

CharacteristicTraining cohort (N=101)Validation cohort (N=44)P value
Age, median (IQR), years50.00 (41.00–58.00)47.50(43.00–59.50)0.957
Sex, N (%)
 Female19 (18.8)7 (15.9)0.854
 Male82 (81.2)37 (84.1)
Etiology of viral hepatitis, N (%)
None22 (21.8)9 (20.5)0.291
 Hepatitis B73 (72.3)31 (70.5)
 Hepatitis C3 (3.0)2 (4.5)
 Hepatitis E2 (2.0)0 (0.0)
 Hepatitis B + C0 (0.0)2 (4.5)
 Hepatitis B + E1 (1.0)0 (0.0)
Drinking history, N (%)
 No87 (86.1)38 (86.4)1.000
 Yes14 (13.9)6 (13.6)
Alpha fetoprotein content, median (IQR)3.16 (1.46–6.36)2.15 (0.82–4.94)0.067
Number of tumors, median (IQR)1.00 (1.00–1.00)1.00 (1.00–1.00)0.469
Neutrophil counts, median (IQR)3.84 (3.19–4.90)3.93 (2.96–4.95)0.877
Lymphocyte counts, median (IQR)1.74 (1.30–2.20)1.69 (1.33–2.15)0.877
Neutrophil-lymphocyte ratio, median (IQR)2.26 (1.65–2.97)2.41 (1.62–3.33)0.752
Preoperative alanine transaminase, median (IQR)35.00 (20.00–51.00)31.50(20.50–47.25)0.380
Preoperative aspartate aminotransferase, median (IQR)35.00 (23.00–45.00)31.00(23.50–44.00)0.523
Preoperative total bilirubin, median (IQR)11.50 (8.00–16.40)10.30 (8.75–14.60)0.940
Preoperative hemoglobin, median (IQR)135.00(120.00–149.00)134.50(120.75–151.00)0.908
Preoperative albumin, median (IQR)39.00 (36.00–41.60)40.70(35.95–42.73)0.541
Preoperative blood platelet counts, median (IQR)172.00(143.00-246.00)195.00(141.25–256.00)0.563
Preoperative prothrombin time, median (IQR)13.70 (13.00–14.20)13.45(12.67–14.20)0.242
Histopathological differentiation degree, N (%)
 Poor64 (63.4)27 (61.4)0.438
 Moderate22 (21.8)7 (15.9)
 Well15 (14.9)10 (22.7)
Tumor with or without microvascular invasion (%)
 With71 (70.3)28 (63.6)0.550
 Without30 (29.7)16 (36.4)

Demographic and clinicopathologic characteristics of the study participants.

Table 2

CharacteristicPilot (N=18)
Age, median (IQR), years52.00 (48.25–57.25)
Sex, N (%)
 Female3 (16.7)
 Male15 (83.3)
Etiology of viral hepatitis, N (%)
 Hepatitis B10 (55.6)
 None8 (44.4)
Drinking history, N (%)
 No18 (100.0)
Alpha fetoprotein content, median (IQR)7.70 (4.42–23.05)
Number of tumors, median (IQR)1.00 (1.00–1.00)
Neutrophil counts, median (IQR)3.54 (3.07–4.67)
Lymphocyte counts, median (IQR)2.02 (1.51–2.68)
Neutrophil-lymphocyte ratio, median (IQR)1.89 (1.15–2.38)
Preoperative alanine transaminase, median (IQR)32.00 (21.00–57.75)
Preoperative aspartate aminotransferase, median (IQR)28.00 (24.25–39.25)
Preoperative total bilirubin, median (IQR)14.25 (10.92–15.20)
Preoperative hemoglobin, median (IQR)143.50 (124.50–148.25)
Preoperative albumin, median (IQR)36.80 (34.55–40.28)
Preoperative blood platelet counts, median (IQR)195.00 (172.50–250.75)
Preoperative prothrombin time, median (IQR)13.35 (12.70–13.85)
Histopathological differentiation degree, N (%)
 Poor13 (72.2)
 Moderate1 (5.6)
 Well4 (22.2)
Tumor with or without microvascular invasion, N (%)
 With11 (61.1)
 Without7 (38.9)

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

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

Table 3

CharacteristicUnivariate logistic regressionMultivariate logistic regression
OR (95% CI)P valueOR (95% CI)P value
Score1.093(1.017-1.176)0.016*1.158(1.041-1.29)0.007*
Age0.991(0.958-1.025)0.596NANA
Sex0.897(0.305-2.635)0.843NANA
Size of tumors0.974(0.849-1.117)0.705NANA
Etiology of viral hepatitis1.077(0.371-3.128)0.892NANA
Drinking history1.969(0.618-6.269)0.252NANA
Alpha fetoprotein content1.447(1.216-1.721)<0.001*1.533(1.22-1.927)<0.001*
Number of tumors0.623(0.111-3.507)0.592NANA
Neutrophil counts1.379(1.09-1.745)0.007*1.567(1.138-2.158)0.006*
Lymphocyte counts0.849(0.441-1.634)0.624NANA
Neutrophil-lymphocyte ratio1.212(0.956-1.536)0.113NANA
Preoperative alanine transaminase1.004(0.995-1.013)0.416NANA
Preoperative aspartate aminotransferase1.007(0.999-1.015)0.081NANA
Preoperative total bilirubin0.995(0.967-1.024)0.747NANA
Preoperative hemoglobin0.966(0.944-0.989)0.004*0.966(0.938-0.995)0.022*
Preoperative albumin0.961(0.883-1.046)0.355NANA
Preoperative blood platelet counts1.007(1.003-1.012)0.003*NANA
Preoperative prothrombin time1.035(0.75-1.427)0.835NANA
Histopathological differentiation degree2.765(1.014-7.538)0.047*NANA

Risk factors for the diagnosis of tumors with microvascular invasion.

*P < 0.05.

Figure 4

Figure 5

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 (2729). 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.

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

Statements

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.

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

    FornerAReigMBruixJ. Hepatocellular carcinoma. Lancet (2018) 391(10127):1301–14. doi: 10.1016/S0140-6736(18)30010-2

  • 2

    MillerKDNogueiraLMariottoABRowlandJHYabroffKRAlfanoCMet al. Cancer treatment and survivorship statistics, 2019. CA Cancer J Clin (2019) 69(5):363–85. doi: 10.3322/caac.21565

  • 3

    FangCZhangPQiX. Digital and intelligent liver surgery in the new era: Prospects and dilemmas. Ebiomedicine (2019) 41:693701. doi: 10.1016/j.ebiom.2019.02.017

  • 4

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

  • 5

    LimKCChowPKAllenJCChiaGSLimMCheowPCet 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

  • 6

    IguchiTShirabeKAishimaSWangHFujitaNNinomiyaMet 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

  • 7

    CostentinCEFerroneCRArellanoRSGanguliSHongTSZhuAX. Hepatocellular carcinoma with macrovascular invasion: Defining the optimal treatment strategy. Liver Cancer (2017) 6(4):360–74. doi: 10.1159/000481315

  • 8

    SumieSNakashimaOOkudaKKuromatsuRKawaguchiANakanoMet 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

  • 9

    RenzulliMBrocchiSCucchettiAMazzottiFGolfieriR. Can current preoperative imaging be used to detect microvascular invasion of hepatocellular carcinoma? Radiology (2015) 279(2):150998. doi: 10.1148/radiol.2015150998

  • 10

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

  • 11

    KornbergAFreesmeyerMRthelEBJandtKKatenkampKSteenbeckJet al. 18F-FDG-Uptake of hepatocellular carcinoma on PET predicts microvascular tumor invasion in liver transplant patients. Am J Transplant (2009) 9(3):592600. doi: 10.1111/j.1600-6143.2008.02516.x

  • 12

    CaiWHeBHuMZhangWXiaoDYuHet al. A radiomics-based nomogram for the preoperative prediction of posthepatectomy liver failure in patients with hepatocellular carcinoma. Surg Oncol (2019) 28:7885. doi: 10.1016/j.suronc.2018.11.013

  • 13

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

  • 14

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

  • 15

    XuXZhangHLLiuQPSunSWZhangJZhuFPet 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

  • 16

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

  • 17

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

  • 18

    LafaroKGrandhiMSHermanJMPawlikTM. The importance of surgical margins in primary malignancies of the liver. J Surg Oncol (2016) 113(3):296303. doi: 10.1002/jso.24123

  • 19

    DuMChenLZhaoJTianFZengHTanYet 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

  • 20

    LeeSKangTWSongKDLeeMWRhimHLimHKet 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

  • 21

    Wen-MingCHBJGuidelineC. 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

  • 22

    HeBHuangCSharpGZhouSHuQFangCet 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

  • 23

    GriethuysenJJMFedorovAParmarCHosnyAAucoinNNarayanVet al. Computational radiomics system to decode the radiographic phenotype. Cancer Res (2017) 77(21):e104–7. doi: 10.1158/0008-5472.CAN-17-0339

  • 24

    MotzerRJBukowskiRMFiglinRAHutsonTEDror MichaelsonMKimSTet al. Prognostic nomogram for Sunitinib in patients with metastatic renal cell carcinoma. Cancer (2008) 113(7):1552–8. doi: 10.1002/cncr.23776

  • 25

    LewisSHectorsSTaouliB. Radiomics of hepatocellular carcinoma. Abdom Radiol (2020). doi: 10.1007/s00261-019-02378-5

  • 26

    LambinPLeijenaarRTHDeistTMPeerlingsJde JongEECvan TimmerenJet 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

  • 27

    FengSTJiaYLiaoBHuangBZhouQLiXet 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

  • 28

    ChenLLiWXianMZhengXLinYLiuBet 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

  • 29

    JiGZhuFXuQWangKWuMTangWet 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

  • 30

    RoayaieSBlumeINThungSNGuidoMFielMIHiotisSet 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

  • 31

    ZhaoHHuaYDaiTHeJTangMFuXet al. Development and validation of a novel predictive scoring model for microvascular invasion in patients with hepatocellular carcinoma. Eur J Radiol (2017) 88(Complete):3240. doi: 10.1016/j.ejrad.2016.12.030

  • 32

    LeiZLiJDongWYongXWangQSiAet 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

Summary

Keywords

hepatocellular carcinoma, microvascular invasion, radiomics, nomogram, computed tomography

Citation

He M, Zhang P, Ma X, He B, Fang C and Jia F (2020) Radiomic Feature-Based Predictive Model for Microvascular Invasion in Patients With Hepatocellular Carcinoma. Front. Oncol. 10:574228. doi: 10.3389/fonc.2020.574228

Received

19 June 2020

Accepted

14 October 2020

Published

05 November 2020

Volume

10 - 2020

Edited by

Dong-Hyun Kim, Northwestern University, United States

Reviewed by

Huijie Jiang, Second Affiliated Hospital of Harbin Medical University, China; Ren Xuequn, Huaihe Hospital of Henan University, China

Updates

Copyright

*Correspondence: Chihua Fang, ; Fucang Jia,

This article was submitted to Cancer Imaging and Image-directed Interventions, a section of the journal Frontiers in Oncology

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics