Secondary Logo

Journal Logo

Selection Between Liver Resection Versus Transarterial Chemoembolization in Hepatocellular Carcinoma

A Multicenter Study

Fu, Sirui MD1,2; Wei, Jingwei PhD2,3; Zhang, Jie MD4; Dong, Di PhD2,3; Song, Jiangdian PhD2,3; Li, Yong MD1; Duan, Chongyang PhD5; Zhang, Shuaitong PhD2,3; Li, Xiaoqun MD6; Gu, Dongsheng PhD2,3; Chen, Xudong MD7; Hao, Xiaohan PhD2,3; He, Xiaofeng MD8; Yan, Jianfeng MD9; Liu, Zhenyu PhD2,3; Tian, Jie PhD2,3,10; Lu, Ligong MD1

Clinical and Translational Gastroenterology: August 2019 - Volume 10 - Issue 8 - p e00070
doi: 10.14309/ctg.0000000000000070
ARTICLE
Open
SDC

OBJECTIVES: Models should be developed to assist choice between liver resection (LR) and transarterial chemoembolization (TACE) for hepatocellular carcinoma.

METHODS: After separating 520 cases from 5 hospitals into training (n = 302) and validation (n = 218) data sets, we weighted the cases to control baseline difference and ensured the causal effect between treatments (LR and TACE) and estimated progression-free survival (PFS) difference. A noninvasive PFS model was constructed with clinical factors, radiological characteristics, and radiomic features. We compared our model with other 4 state-of-the-art models. Finally, patients were classified into subgroups with and without significant PFS difference between treatments.

RESULTS: Our model included treatments, age, sex, modified Barcelona Clinic Liver Cancer stage, fusion lesions, hepatocellular carcinoma capsule, and 3 radiomic features, with good discrimination and calibrations (area under the curve for 3-year PFS was 0.80 in the training data set and 0.75 in the validation data set; similar results were achieved in 1- and 2-year PFS). The model had better accuracy than the other 4 models. A nomogram was built, with different scores assigned for LR and TACE. Separated by the threshold of score difference between treatments, for some patients, LR provided longer PFS and might be the better option (training: hazard ratio [HR] = 0.50, P = 0.014; validation: HR = 0.52, P = 0.026); in the others, LR provided similar PFS with TACE (training: HR = 0.84, P = 0.388; validation: HR = 1.14, P = 0.614). TACE may be better because it was less invasive.

DISCUSSION: We propose an individualized model predicting PFS difference between LR and TACE to assist in the optimal treatment choice.

1Zhuhai Interventional Medical Center, Zhuhai Precision Medical Center, Zhuhai People's Hospital, Zhuhai Hospital of Jinan University, Zhuhai, China;

2Key Laboratory of Molecular Imaging, Institute of Automation, Chinese Academy of Sciences, Beijing, China;

3University of Chinese Academy of Sciences, Beijing, China;

4Department of Radiology, Zhuhai People's Hospital, Zhuhai Hospital of Jinan University, Zhuhai, China;

5Department of Biostatistics, School of Public Health, Southern Medical University, Guangzhou, China;

6Department of Interventional Treatment, Zhongshan City People's Hospital, Zhongshan, China;

7Department of Radiology, Shenzhen People's Hospital, Shenzhen, China;

8Interventional Diagnosis and Treatment Department, Nanfang Hospital, Southern Medical University, Guangzhou, China;

9Department of Radiology, Yangjiang People's hospital, Yangjiang, China;

10Beijing Advanced Innovation Center for Big Data-Based Precision Medicine, Beihang University, Beijing.

Correspondence: Jie Tian, PhD. E-mail: tian@ieee.org. Ligong Lu, MD. E-mail: llg0902@sina.com.

SUPPLEMENTARY MATERIAL accompanies this paper at http://links.lww.com/CTG/A84, http://links.lww.com/CTG/A85, and http://links.lww.com/CTG/A86

Received February 25, 2019

Accepted June 29, 2019

Online date: August 1, 2019

This is an open access article distributed under the Creative Commons Attribution License 4.0 (CCBY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Back to Top | Article Outline

INTRODUCTION

Worldwide, liver cancer ranks high in terms of incidence and mortality; 70%–90% of these cases are hepatocellular carcinoma (HCC) (1). HCCs are often accompanied by underlying liver dysfunction; therefore, death is caused not only by tumor burden but also by deterioration of liver function (2), which makes it important to avoid unnecessary trauma. For HCC, liver resection (LR) is curative but is highly traumatic, whereas transarterial chemoembolization (TACE) is minimally invasive but may leave some residual tumor. With the advancement of technology, their adaptation has expanded and overlapped (3–8). The definition of “unresectable HCC” is somewhat ambiguous in the guidelines of the European Association for the Study of the Liver (9) and the American Association for the Study of Liver Diseases (10). The latest National Comprehensive Cancer Network guideline states that the performance of LR (in some previous TACE-treated cohorts) is “controversial” (11). Researchers suggest that we should move from an approach of “what can be done” to a “what is worth doing” process (12,13). Thus, selecting treatments for patients based on predecided parameters may be rational (3,14). The existing methods for treatment selection have some limitations, such as subjective outcomes and a lack of validation. A new model that can help making the treatment choice between LR and TACE is required. Furthermore, the model should use proper method to ensure that the survival difference between LR and TACE is causal effect (15). To this end, individualized models to assistant choice between LR and TACE are needed.

During HCC management, biopsy is not a necessity (9,10). Hence, the choice between LR and TACE would be better assisted by noninvasive methods. To collect information for the noninvasive methods, in addition to traditional clinical factors and radiological characteristics (16), radiomics may also be helpful. Radiomics can quantitatively mine data from medical images, including but not limited to size/shape, heterogeneity/texture, and relationships with surrounding tissues. Radiomics can improve diagnostic and prognostic accuracy (17) and is valuable for decision making in cancer (18,19) and liver disease (20).

Therefore, we conducted this multicenter study on patients with HCC, aiming to establish noninvasive individualized models by combining clinical factors, radiological characteristics, and radiomic features. Assisted by the model, we planned to predict prognostic difference associated with LR and TACE and then classified the population for the optimal treatment accordingly.

Back to Top | Article Outline

METHODS

Patients

HCC cases from Nanfang Hospital, Shenzhen People's Hospital, Yangjiang People's Hospital, Zhongshan City People's Hospital, and Zhuhai People's Hospital in China were reviewed. After screening the electronic medical record system from 2008 to 2016, we identified 520 HCC cases for which computed tomography (CT) records at diagnosis were available. Inclusion criteria were as follows: (i) HCC confirmed pathologically or clinically; (ii) patients initially treated by LR or TACE; and (iii) at least one radiological disease progression (PD) confirmed by CT/MRI or a follow-up of more than 1 year before the end date. Exclusion criteria were as follows: (i) initially treated by ablation or other treatments; (ii) irregular follow-up; (iii) administration of other combination therapies before PD (such as ablation or systemic therapy); (iv) incomplete initial treatment of intrahepatic lesion (such as positive margin in LR and untreated lesion in TACE); and (v) severe medical comorbidities such as extensive cardiovascular disease (see Figure S1, Supplemental Digital Content 3, http://links.lww.com/CTG/A86).

The 2018 American Association for the Study of Liver Diseases guideline recommended the modified Barcelona Clinic Liver Cancer (BCLC) staging system rather than the TNM system (10). Single HCC > 5 cm was not an optimal candidate for LR and appeared to have a worse prognosis than HCC ≤ 5 cm (10), which resulted in the adoption of stage AB in research and discussed in the European Association for the Study of the Liver guideline (3,9). Thus, we tested whether the application of modified BCLC and stage AB would improve the performance of our model.

The study protocols were approved by the Ethics Committee of Zhuhai People's Hospital. The requirement for informed consent to use the patients' data for medical researches was waived because the data were collected retrospectively. All patient records and information were anonymized and deidentified before analysis.

Back to Top | Article Outline

Follow-ups and treatments

For the included patients, the follow-up was scheduled every 4–6 weeks before PD or at least for 1 year (without PD), after that, the interval might be increased from 3 to 6 months. Every follow-up included chest radiograms, abdominal CT or MRI, and other necessary laboratory tests. Additional CT or MRI scans were obtained if extrahepatic metastasis was suspected.

The initial treatment between LR and TACE was decided based on the same criteria across hospitals, according to tumor characteristics referring to guidelines (9–11), liver functional status, and the patients' choice. LR was performed with the aim of removing all intrahepatic macroscopic tumors with a margin confirmed by pathological examination. Conventional TACE was performed using the standard method as selectively as possible (9–11). Before PD, no additional LR, TACE, or ablation procedures were performed. Inspired by a previous study (3), for BCLC stage A, some patients with early-stage HCC preferred TACE over LR because TACE was less invasive. For BCLC stage C, patients with suspicious extrahepatic nodule proven to be metastases during the follow-ups, and patients with macrovascular invasion at diagnosis were not excluded because of (i) the limited efficacy of targeted therapies recommended by guidelines (21); (ii) the potential benefit derived from LR and TACE for these patients (3,5,6,8,11,22–24) and treatments should not be delayed (especially for those with macrovascular invasion in portal vein branches or suspicious regional lymph node metastasis). After informing patients of the potential consequences and obtaining informed consent for LR or TACE procedures, the physicians performed LR/TACE as the initial treatment.

Back to Top | Article Outline

Outcomes

Based on ethical considerations, after PD, to more effectively control the disease, patients could accept other therapies, even those beyond first-line treatments recommended by guidelines (25); however, this might introduce bias in the analysis of overall survival. Therefore, we used progression-free survival (PFS) as the end point (26). PFS was defined as the time from diagnosis to PD or death. PD was defined according to the modified Response Evaluation Criteria in Solid Tumors criteria (27), referring to both target lesion and nontarget lesions (such as macrovascular invasion or extrahepatic metastases).

Back to Top | Article Outline

Candidate clinical factors and radiological characteristics

Besides the abovementioned baseline factors and radiomic features, we also included several other factors during model construction: (i) neutrophil-to-lymphocyte ratio; (ii) HCC location: lobe (classified as left, right, and cross); surface (whether close to the liver capsule or not, classified as negative or positive); and (iii) radiological characteristics: fusion lesions, boundary fusion of >2 lesions, classified as absent (Figure 1a) or present (Figure 1b); invasive shape, protrusions like crab foot, classified as noninvasive (Figure 1c) or invasive (Figure 1d); HCC capsule, classified as absent (Figure 1e), unintegral (Figure 1f), or integral (Figure 1g); HCC capsule breakthrough (Figure 1h), HCC growth beyond a preexisting capsule, classified as absent or present; corona enhancement (16) (Figure 1i), a transient zone or rim of enhancement around HCC, classified as absent or present; corona with low attenuation (Figure 1j), classified as absent or present; mosaic architecture (16) (Figure 1k), a mass of randomly distributed internal nodules or compartments differing in enhancement, classified as absent or present; nodule-in-nodule architecture (16) (Figure 1l), the presence of a nodule within a larger nodule or mass, classified as absent or present; and HCC with enhancement, classified as <25% (Figure 1m), 25%–50% (Figure 1n), 50%–75% (Figure 1o), or >75% (Figure 1p).

Figure 1

Figure 1

Back to Top | Article Outline

Radiomic feature extraction

Baseline CT images were obtained using the scanners in the collaborative hospitals (see Table S1, Supplemental Digital Content 2, http://links.lww.com/CTG/A85). All CT images were retrieved from the picture archiving and communication system.

Because HCC capsule was identified in the portal phase rather than the arterial phase (16), which could increase the accuracy of lesion segmentation, the portal phase was used for feature extraction. The target lesion was selected according to the modified Response Evaluation Criteria in Solid Tumor (27) (see Text S1, Supplemental Digital Content 1, http://links.lww.com/CTG/A84). In total, we extracted 708 radiomic features (see Text S2, Supplemental Digital Content 1, http://links.lww.com/CTG/A84). After randomly sampling 20% cases per hospital, we used morphologic perturbations to test the robustness and redundancy (see Text S3, Supplemental Digital Content 1, http://links.lww.com/CTG/A84) (28). The related program, instruction, and example of patients can be downloaded from http://www.radiomics.net.cn/post/108.

Back to Top | Article Outline

Statistical analysis

Continuous variables were expressed as mean (SD) or median (range) based on whether they were normally distributed or not, and groups were compared using the t-test or Wilcoxon rank-sum test as appropriate. Categorical variables were expressed as percentages, and groups were compared using the Pearson χ2 test or Fisher exact test.

For model construction, after randomly splitting the patients into the training and validation data sets, we used the following steps to develop our individualized model. First, we collected the clinical factors, radiological characteristics, and radiomic features for our model. Second, inverse probability of treatment weighting (IPTW) based on the propensity score was used to weight the patients in the training and validation data sets separately (see Text S4, Supplemental Digital Content 1, http://links.lww.com/CTG/A84) (29,30), because our model aimed to ensure the causal effect (31) in treatment outcome differences between LR and TACE. Third, the Cox proportional hazard model was applied to predict the PFS associated with LR and TACE using the weighted training data set. The candidate factors and their interactions with treatments were selected by the Least Absolute Shrinkage and Selection Operator method (32) and the backward stepwise selection method using the Akaike information criterion. We constructed 3 models to evaluate whether the combination of clinical factors, radiological characteristics, and radiomic features was necessary: ModelCR included clinical factors (C) and radiological characteristics (R); ModelR, radiomic features (R); and ModelCRR, clinical factors (C), radiological characteristics (R), and radiomic features (R). We also compared our developed models with 2 other existing models (ITA.LI.CA (33) and CLIP (34)) reported to be superior. The predictive accuracy of the developed models was assessed by both discriminations measured by the time-dependent receiver operating characteristic curve for 1-, 2-, and 3-year PFS and Harrell concordance index (C-index), and calibration evaluated by the calibration plot. A nomogram of the best developed model was obtained. Fourth, we calculated the score difference between LR and TACE of the best model; the cutoff value of the score difference was used to subdivide the patients for classification. We used Kaplan-Meier plots to compare the difference between the subgroups. The abovementioned steps are illustrated in Figure 2.

Figure 2

Figure 2

All statistical tests performed were two sided, and P values <0.05 were considered statistically significant. We analyzed data by the R statistical package (http://www.r-project.org/).

Back to Top | Article Outline

RESULTS

Characteristics of the study population

Among the 520 patients, the training and validation data sets had 302 and 218 patients, respectively (Table 1). The patients initially treated by LR and TACE were 97 (32%) and 205 (68%) in the training data set and 65 (30%) and 153 (70%) in the validation data set, respectively. During follow-up, 338 patients (training: 193; validation: 145) showed PD, and 182 patients (training: 102; validation: 80) died. Baseline differences between LR and TACE groups before IPTW were shown in Table S2 (see Supplemental Digital Content 2, http://links.lww.com/CTG/A85).

Table 1

Table 1

Back to Top | Article Outline

Radiomic features and IPTW results

For the 708 radiomic features, 607 features with intraclass correlation coefficient > 0.75 were used for the final analysis (see Figure S2, Supplemental Digital Content 3, http://links.lww.com/CTG/A86). IPTW controlled the baseline differences between LR and TACE groups such as tumor burden, liver function, and chronic liver disease (see Text S5, Supplemental Digital Content 1, http://links.lww.com/CTG/A84) in both data sets (see Figure S3a,b, Supplemental Digital Content 3, http://links.lww.com/CTG/A86) and most standardized differences <10% (see Table S2, Supplemental Digital Content 2, http://links.lww.com/CTG/A85). After IPTW, there were 512 patients included for PFS analysis (training: 298 patients; validation: 214 patients).

Back to Top | Article Outline

Model development and validation

Among the 5 models, ModelCRR, containing treatments, age, sex, BCLC stage, fusion lesions, HCC capsule, and 3 radiomic features (FOS_Kurtosis, CO_IV, and POF_entropy), showed the best area under the curve for 1-, 2- and 3-year PFS (training: 0.78, 0.80, and 0.80; validation: 0.73, 0.74, and 0.75, respectively, Figure 3). The differences between ModelCRR and the other 4 models were statistically significant (see Table S3, Supplemental Digital Content 2, http://links.lww.com/CTG/A85), with ModelCRR having the highest Harrell C-index (ModelCR: 0.696; ModelR: 0.626; ModelCRR: 0.707; ITA.LI.CA: 0.650; CLIP: 0.620). ModelCRR showed good calibration in both data sets (Figure 4a,b). In addition, ModelCRR showed an improved performance than the model without modified BCLC and stage AB (see Text S6, Supplemental Digital Content 1, http://links.lww.com/CTG/A84, and see Figure S4, Figure S5, Supplemental Digital Content 3, http://links.lww.com/CTG/A86).

Figure 3

Figure 3

Figure 4

Figure 4

Based on the above results, a nomogram of ModelCRR was constructed (Figure 4c). In the nomogram, 3 radiomic features were included, with their calculation formulas and schematic diagram displayed (Figure 4d–f) and a detailed explanation provided in the discussion. CO_IV was identified to have an interaction with treatments, and therefore, different scores were assigned for CO_IV, facilitating the predicted PFS for LR and TACE each.

Back to Top | Article Outline

Population classification

In the IPTW-weighted patients, the PFS difference between LR and TACE group was significant in the training data set: hazard ratio (HR) = 0.67 (95% confidence interval [CI]: 0.47–0.91), P = 0.012 (Figure 5a; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85); but not in the validation data set: HR = 0.82 (95% CI: 0.57–1.19), P = 0.304 (Figure 5b; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85). The inconsistency between data sets indicated that not all patients got PFS benefit from LR, and patients should be classified into subgroups for more reasonable treatment decision. Thus, we calculated the score difference of ModelCRR (ΔModelCRR) by the formula ΔModelCRR = ModelCRR (LR score)−ModelCRR (TACE score). We used −5.00 as the threshold because we wanted to ensure PFS benefit for patients identified to be suitable for LR.

Figure 5

Figure 5

In patients with ΔModelCRR ≤ −5.00 (LR− and TACE− subgroups), LR provided longer PFS. We obtained an HR = 0.50 (95% CI: 0.29–0.87), with P = 0.014 in the training data set (Figure 5c; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85), and HR = 0.52 (95% CI: 0.29–0.93), with P = 0.026 in the validation data set (Figure 5d; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85). Thus, LR might be the first option for these patients.

In patients with ΔModelCRR > −5.00 (LR+ and TACE+ subgroups), there were no differences in PFS. In the training data set, the HR was 0.84 (95% CI: 0.57–1.25), with P = 0.388 (Figure 5e; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85); in the validation data set, the HR was 1.14 (95% CI: 0.69–1.85), with P = 0.614 (Figure 5f; see Table S4, Supplemental Digital Content 2, http://links.lww.com/CTG/A85). Thus, TACE might be a better choice because it was less invasive and did not cause decreased PFS.

We tested whether the classification was still informative in patients without extrahepatic metastasis as in a previous study (3), and the subgroup demonstrated similar results (see Figure S6, Supplemental Digital Content 3, http://links.lww.com/CTG/A86). Even if we excluded patients in BCLC stages A and C, our model could still classify patients either for LR or TACE (Figure 6). In patients with ΔModelCRR ≤ −5.00, LR provided longer PFS. We obtained an HR = 0.41 (95% CI: 0.20–0.85), with P = 0.016 in the training data set, and HR = 0.32 (95% CI: 0.14–0.70), with P = 0.004 in the validation data set. For patients with ΔModelCRR > −5.00, there were no differences in PFS: training data set, HR = 0.78 (95% CI: 0.43–1.42), with P = 0.422, validation dataset, HR = 1.57 (95% CI: 0.76–3.23), with P = 0.225. Representative patients of the 4 subgroups are shown in Figure 7.

Figure 6

Figure 6

Figure 7

Figure 7

Back to Top | Article Outline

DISCUSSION

In this study, based on IPTW-weighted data, we used a noninvasive model to predict PFS of LR and TACE, which had good discrimination and calibration in both training and validation data sets. Then, subdivided by the threshold of the score difference, for some patients, LR provided better PFS than TACE, which suggested LR to be a potential better choice for increased PFS. However, for the other patients, because LR and TACE had similar PFS, TACE seemed a better option to control unnecessary trauma and risks. These conclusions were especially useful for patients in BCLC stages AB and B, in which there were more controversies between the choice of LR and TACE.

The choice between LR and TACE mainly depends on the difference in survival benefit; however, the patients' preference between controlling trauma/risks and survival difference should also be considered (35). Thus, besides classification, doctors and patients also need to know the predicted PFS associated with each treatment. For example, with a 6-month predicted PFS difference, a young father may tend to accept LR, whereas a 60-year-old patient may want TACE to control sufferings. However, when the PFS difference is 1-year, both may choose LR as the first option. Thus, after knowing the exact PFS difference predicted by the nomogram, the doctors may adjust the classification threshold of −5.00 according to clinical situation and patients' preference between risks/trauma and PFS benefit.

Compared with previous studies on HCC prognosis (36), our study not only performed PFS prediction on both LR and TACE cohorts but also established a procedure to compare the PFS between LR and TACE. Besides traditional clinical factors, such as BCLC stage including maximum diameter and lesion number, we also combined radiological characteristics and radiomic features. Compared with other clinical factors such as alpha-fetoprotein, these factors facilitated data mining from conventional CT images and provided more information of intratumor heterogeneity and were included in the final model. Compared with previous studies on treatments decision (3,14), instead of based on “post-treatment regret of doctors”, our model used more objective PFS end point, and we constructed our model on IPTW-weighted data set, which increased the comparability between LR and TACE groups, ensuring the causal effect between treatments and PFS difference. Because all our factors could be extracted through noninvasive methods, no additional biopsy would be required. Thus, treatment decision could be assisted with minimal disturbance of current HCC management.

In this study, besides conventional radiomic features (17), we also included “peer-off” features indicating HCC heterogeneity from the outside to the inside. Similar to CT attenuation, radiomic features reflect the granular changes in radiological images. As we cannot provide a definite causal relationship between CT attenuation and pathological results, we also cannot explain radiomic features pathologically. However, we may provide an explanation for the identified features based on their extraction methods. FOS_Kurtosis is Kurtosis of first-order statistics (0) and represents the first-order features on an original image. If the kurtosis is >3, the distribution of the intensity is sharper than a normal distribution (Figure 4d). CO_IV is inverse variance of co-occurrence (1,2) and represents the co-occurrence of the textural features of inverse variance on the x-direction high-pass and the y-direction low-pass filtered image. CO_IV emphasizes more edge information in the x-direction, and greater inverse variance indicates a more heterogeneous texture (Figure 4e). In our study, we found that when CO_IV decreases, the PFS between LR and TACE decreased at different rates, so independent scores of CO_IV were assigned for LR and TACE separately. Finally, POF_entropy is a peer-off feature of entropy (9) and represents the peel-off feature entropy in the innermost layer. It measures the texture randomness or irregularity of the layer (Figure 4f).

Our study has some limitations. First, because we perform IPTW in the training and validation separately, giving them a certain degree of independence, the potential bias of the analysis procedures can be controlled. Still, considering the complexity of the issue, an external validation data set is needed to further test our model. Second, for the differences of HCC between Asia and the Western countries, validation in a non-Chinese population is needed before our conclusions are applied in western cohorts. Third, as patients are divided into 8 subgroups for PFS analysis, the sample size of the study population may be relatively limited, and we were unable to perform more detailed analysis according to different BCLC stages. But ΔModelCRR still successfully classifies patients based on PFS difference, which may also be helpful in controlling unnecessary surgical trauma and risks while ensuring better treatment for patients. Hopefully, these shortcomings can be addressed with further international research and cooperation.

In conclusion, our individualized model predicts the PFS difference between LR and TACE and assists the treatment decision: in selected patients, LR can be used to increase PFS benefit; in the others, TACE can be preferable to avoid unnecessary trauma and risks.

Back to Top | Article Outline

CONFLICTS OF INTEREST

Guarantors of the article: Jie Tian, PhD and Ligong Lu, MD.

Specific author contributions: Sirui Fu, MD, Jingwei Wei, PhD, Jie Zhang, MD, and Di Dong, PhD, contributed equally to this work. S.F., J.W., J.Z., and D.D. conceived and designed the project with supervision from J.T. and L.L. J.S., Y.L., X.L., X.C., X.H., and J.Y. acquired the data. J.Z. and J.Y. segmented the data manually. S.Z., D.G., X.H., and Z.L. extracted radiomic features. C.D. performed the statistical analysis. All authors were involved in drafting and reviewing the manuscript and approved the final manuscript for submission.

Financial support: This work was supported by the National Key R&D Program of China (Nos. 2017YFA0205200, 2017YFC1308700, 2017YFC1308701, 2017YFC1309100, and 2016CZYD0001), the National Natural Science Foundation of China (Nos. 81571785, 81771957, 81227901, 61231004, 81771924, 81501616, 81671851, 81527805, and 81501549), the Natural Science Foundation of Guangdong Province, China (Nos. 2016A030311055 and 2016A030313770), the Science and Technology Service Network Initiative of the Chinese Academy of Sciences (KFJ-SW-STS-160), the Instrument Developing Project of the Chinese Academy of Sciences (YZ201502), the Beijing Municipal Science and Technology Commission (Z161100002616022), and the Youth Innovation Promotion Association CAS.

Potential competing interests: None.

Back to Top | Article Outline

Study Highlights

WHAT IS KNOWN

  • ✓ LR and TACE are first-line treatments with overlapped adaptation in HCC.
  • ✓ Choice between them should be assisted by noninvasive models to avoid unnecessary biopsies.
Back to Top | Article Outline

WHAT IS NEW HERE

  • ✓ Our model predicted and compared PFS between LR and TACE.
  • ✓ Subdivided by the score difference, LR provides longer PFS than TACE in some patients.
  • ✓ In the other patients, the 2 treatments provided similar PFS.
Back to Top | Article Outline

TRANSLATIONAL IMPACT

  • ✓ Our model predicted PFS difference between LR and TACE and assisted in choosing the optimal treatment.
Back to Top | Article Outline

REFERENCES

1. Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2018;68:394–424.
2. Dhir M, Melin AA, Douaiher J, et al. A review and update of treatment options and controversies in the management of hepatocellular carcinoma. Ann Surg 2016;263:1112–25.
3. Vitale A, Burra P, Frigo AC, et al. Survival benefit of liver resection for patients with hepatocellular carcinoma across different Barcelona clinic liver cancer stages: A multicentre study. J Hepatol 2015;62:617–24.
4. Yin L, Li H, Li AJ, et al. Partial hepatectomy vs. transcatheter arterial chemoembolization for resectable multiple hepatocellular carcinoma beyond Milan criteria: A RCT. J Hepatol 2014;61:82–8.
5. Kokudo T, Hasegawa K, Matsuyama Y, et al. Liver resection for hepatocellular carcinoma associated with hepatic vein invasion: A Japanese nationwide survey. Hepatology 2017;66:510–7.
6. Zhu K, Chen J, Lai L, et al. Hepatocellular carcinoma with portal vein tumor thrombus: Treatment with transarterial chemoembolization combined with sorafenib—A retrospective controlled study. Radiology 2014;272:284–93.
7. Yang T, Lin C, Zhai J, et al. Surgical resection for advanced hepatocellular carcinoma according to Barcelona Clinic Liver Cancer (BCLC) staging. J Cancer Res Clin Oncol 2012;138:1121–9.
8. Pinter M, Hucke F, Graziadei I, et al. Advanced-stage hepatocellular carcinoma: Transarterial chemoembolization versus sorafenib. Radiology 2012;263:590–9.
9. Galle PR, Forner A, Llovet JM, et al. EASL clinical practice guidelines: Management of hepatocellular carcinoma. J Hepatol 2018;69:182–236.
10. Marrero JA, Kulik LM, Sirlin CB, et al. Diagnosis, staging, and management of hepatocellular carcinoma: 2018 practice guidance by the by the American Association for the study of liver diseases. Hepatology 2018;68:723–50.
11. National Comprehensive Cancer Network. Hepatobiliary Cancers (Version 1.2018). (https://www.nccn.org/professionals/physician_gls/pdf/hepatobiliary.pdf). Accessed February 14, 2018.
12. Dhir M, Are C. Response: we still in an update process of the BCLC system. Ann Surg 2018;267:e24–5.
13. Romagnoli R, Mazzaferro V, Bruix J. Surgical resection for hepatocellular carcinoma: Moving from what can be done to what is worth doing. Hepatology 2015;62:340–2.
14. Cucchetti A, Djulbegovic B, Tsalatsanis A, et al. When to perform hepatic resection for intermediate-stage hepatocellular carcinoma. Hepatology 2015;61:905–14.
15. Ma J, Hobbs BP, Stingo FC. Statistical methods for establishing personalized treatment rules in oncology. Biomed Res Int 2015;2015:670691.
16. Choi JY, Lee JM, Sirlin CB. CT and MR imaging diagnosis and staging of hepatocellular carcinoma: Part II. Extracellular agents, hepatobiliary agents, and ancillary imaging features. Radiology 2014;273:30–50.
17. Lambin P, Leijenaar RTH, Deist TM, et al. Radiomics: The bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol 2017;14:749–62.
18. Sun R, Limkin EJ, Vakalopoulou M, et al. A radiomics approach to assess tumour-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: An imaging biomarker, retrospective multicohort study. Lancet Oncol 2018;19:1180–91.
19. Song J, Shi J, Dong D, et al. A new approach to predict progression-free survival in stage IV EGFR-mutant NSCLC patients with EGFR-TKI therapy. Clin Cancer Res 2018;24:3583–92.
20. Wang K, Lu X, Zhou H, et al. Deep learning radiomics of shear wave elastography significantly improved diagnostic performance for assessing liver fibrosis in chronic hepatitis B: A prospective multicentre study. Gut 2019;68:729–41.
21. Kudo M, Finn RS, Qin S, et al. Lenvatinib versus sorafenib in first-line treatment of patients with unresectable hepatocellular carcinoma: A randomised phase 3 non-inferiority trial. Lancet 2018;391:1163–73.
22. Zhong JH, Ke Y, Gong WF, et al. Hepatic resection associated with good survival for selected patients with intermediate and advanced-stage hepatocellular carcinoma. Ann Surg 2014;260:329–40.
23. Zhong JH, Ke Y, Wang YY, et al. Liver resection for patients with hepatocellular carcinoma and macrovascular invasion, multiple tumours, or portal hypertension. Gut 2015;64:520–1.
24. Mao K, Yan Y, Zhang J, et al. The impact of liver resection on survival outcomes of hepatocellular carcinoma patients with extrahepatic metastases: A propensity score matching study. Cancer Med-Us 2018;7:4475–84.
25. Chen K, Chen G, Wang H, et al. Increased survival in hepatocellular carcinoma with iodine-125 implantation plus radiofrequency ablation: A prospective randomized controlled trial. J Hepatol 2014;61:1304–11.
26. Wilson MK, Karakasis K, Oza AM. Outcomes and endpoints in trials of cancer treatment: The past, present, and future. Lancet Oncol 2015;16:e32–42.
27. Lencioni R, Llovet JM. Modified RECIST (mRECIST) assessment for hepatocellular carcinoma. Semin Liver Dis 2010;30:52–60.
28. Cui Y, Song J, Pollom E, et al. Quantitative analysis of (18)F-fluorodeoxyglucose positron emission tomography identifies novel prognostic imaging biomarkers in locally advanced pancreatic cancer patients treated with stereotactic body radiation therapy. Int J Radiat Oncol Biol Phys 2016;96:102–9.
29. Mohkam K, Dumont PN, Manichon AF, et al. No-touch multibipolar radiofrequency ablation vs. surgical resection for solitary hepatocellular carcinoma ranging from 2 to 5 cm. J Hepatol 2018;68:1172–80.
30. Wahl DR, Stenmark MH, Tao Y, et al. Outcomes after stereotactic body radiotherapy or radiofrequency ablation for hepatocellular carcinoma. J Clin Oncol 2016;34:452–9.
31. Ma H, Jm R. Causal Inference. Chapman & Hall/CRC: Boca Raton, 2018.
32. Callegaro A, Spiessens B, Dizier B, et al. Testing interaction between treatment and high-dimensional covariates in randomized clinical trials. Biom J 2017;59:672–84.
33. Farinati F, Vitale A, Spolverato G, et al. Development and validation of a new prognostic system for patients with hepatocellular carcinoma. PLoS Med 2016;13:e1002006.
34. Liu PH, Hsu CY, Hsia CY, et al. Prognosis of hepatocellular carcinoma: Assessment of eleven staging systems. J Hepatol 2016;64:601–8.
35. Detry MA, Lewis RJ. The intention-to-treat principle: How to assess the true effect of choosing a medical treatment. JAMA 2014;312:85–6.
36. He W, Peng B, Tang Y, et al. Nomogram to predict survival of patients with recurrence of hepatocellular carcinoma after surgery. Clin Gastroenterol Hepatol 2018;16:756–64.

Supplemental Digital Content

Back to Top | Article Outline
© 2019 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of The American College of Gastroenterology