Baseline characteristics of TACE + SORA and TACE groups
A total of 354 patients with advanced liver cancer and recurrent primary liver cancer after radical surgery were included in the study.There were 101 cases in TACE + SORA group and 253 cases in TACE group. The median follow-up time was 17.0(4.0-41.0)months.The mean follow-up time was 18.4±6.9 months. There were statistical differences in age (P=0.003), the level of α-L-Fucosidase(AFU) (P=0.016), extrahepatic metastasis whether or not(P<0.1), hemangioma thrombus whether or not(P=0.001) BCLC staging (P=0.016) and ECOG score (P=0.002) between TACE + SORA group and TACE group, as shown in Table 1. After PSM according to the 1:2 (81cases were in TACE + SORA group, 162 cases were in TACE group), the baseline datas were well balanced. After PSM, P values were > 0.05, as shown in Table 1.
Table 1
Comparison of baseline characteristics of patients between TACE+SORA group and TACE group (pre-PSM & post-PSM)
|
Parameters
|
pre-PSM
|
post-PSM
|
TACE+SORA(n=101)
|
TACE
(n=253)
|
P
|
TACE+SORA
(n=81)
|
TACE
(n=162)
|
P
|
Age
|
(years)
|
55.7±11.1
|
59.6±11.2
|
0.003
|
57.43±10.74
|
57.80±10.60
|
0.799
|
Sex
|
Male
|
87(86.1%)
|
213(84.2%)
|
0.645
|
68(84.0%)
|
136(84.0%)
|
1.000
|
Female
|
14(13.9%)
|
40(15.8%)
|
13(16.0%)
|
26(16.0%)
|
Etiologies
|
HBV
|
82(81.2%)
|
210(83.0%)
|
0.476
|
65(80.2%)
|
135(83.3%)
|
0.720
|
HCV
|
3(3.0%)
|
14(5.5%)
|
3(3.7%)
|
7(4.3%)
|
Others
|
16(15.8%)
|
29(11.5%)
|
13(16.0%)
|
20(12.3%)
|
ALT(U/L)
|
50.45±45.97
|
49.94±49.45
|
0.929
|
51.18±50.32
|
53.53±51.92
|
0.737
|
AFU(U/L)
|
34.53±15.55
|
30.24±13.63
|
0.016
|
32.77±14.44
|
31.07±14.22
|
0.385
|
Albmiun(g/L)
|
39.90±5.63
|
39.91±5.09
|
0.998
|
39.67±5.81
|
39.28±5.80
|
0.621
|
AFP levels
|
≤400ng/mL
|
66(65.3%)
|
143(56.5%)
|
0.127
|
50(61.7%)
|
105(64.8%)
|
0.637
|
>400ng/mL
|
35(34.7%)
|
110(43.5%)
|
31(38.3%)
|
57(35.2%)
|
ALBI grades
|
1
|
57(56.4%)
|
137(54.2%)
|
0.472
|
45(55.6%)
|
88(54.3%)
|
0.983
|
2
|
42(41.6%)
|
114(45.1%)
|
35(43.2)
|
72(44.4%)
|
3
|
2(2.0%)
|
2(0.8%)
|
1(1.2%)
|
2(1.2%)
|
Child-Pugh grades
|
5points
|
64(63.4%)
|
166(65.6%)
|
0.727
|
51(63.0%)
|
103(63.6%)
|
0.960
|
6points
|
22(21.8%)
|
56(22.1%)
|
18(22.2%)
|
37(22.8%)
|
B
|
17(16.8%)
|
31(12.3%)
|
12(14.8%)
|
22(13.6%)
|
Intrahepatic metastasis
|
No
|
35(34.7%)
|
89(35.2%)
|
0.926
|
30(37.0%)
|
74(45.7%)
|
0.199
|
Yes
|
66(65.3%)
|
164(64.8%)
|
51(63.0%)
|
88(54.3%)
|
Extrahepatic metastasis
|
No
|
83(82.2%)
|
187(73.9%)
|
0.099
|
66(81.5%)
|
120(82.1%)
|
0.906
|
Yes
|
18(17.8%)
|
66(26.1%)
|
15(18.5%)
|
29(17.9%)
|
Tumor thrombus
|
No
|
45(44.6%)
|
162(64.0%)
|
0.001
|
42(51.9%)
|
85(51.2%)
|
0.928
|
Yes
|
56(55.4%)
|
91(36.0%)
|
39(48.1%)
|
79(48.8%)
|
Tumor burdens
|
≤6
|
19(18.8%)
|
50(19.8%)
|
0.608
|
15(18.5%)
|
34(21.0%)
|
0.890
|
6<and≤12
|
52(51.5%)
|
116(45.8%)
|
42(51.9%)
|
81(0.0%)
|
>12
|
30(29.7%)
|
87(34.4%)
|
24(29.8%)
|
47(29.0%)
|
BCLC staging
|
B
|
36(35.6%)
|
126(49.8%)
|
0.016
|
31(38.3%)
|
71(43.8%)
|
0.408
|
C
|
65(64.4%)
|
127(50.2%)
|
50(61.7%)
|
91(56.2%)
|
ECOG score
|
0
|
88(87.1%)
|
181(71.5%)
|
0.002
|
66(81.5%)
|
122(75.3%)
|
0.278
|
1
|
13(12.9%)
|
72(28.5.0%)
|
15(18.5%)
|
40(24.7%)
|
Distribution of liver cancer patients with different child Pugh grades in different ALBI grades
The distribution of liver cancer patients with different child Pugh grades in different ALBI grades was shown in Table 2. In 243 patients, patients of ALBI grade 1,2,3 contained 133(54.7%) cases、107(44.0%) cases and 3(1.2%)cases. In patients with ALBI grade 1, patients with liver function of Child-Pugh 5points, 6 points and grade B were 111(83.5%) cases, 18(13.5%) and 4(3.0%) cases. In patients with ALBI grade 1, patients with liver function of Child-Pugh 5points, 6 points and grade B were 43(40.2%), 37(34.6%) and 27(25.2%) cases. All of the 3 patients with ALBI grade 3 were Child-Pugh grade B.
In 154 patients with Child-Pugh grade A_5 points, the distribution in ALBI grade 1 was more than that in ALBI grade 2(72.1% vs 27.9%). In 55 patients with Child-Pugh grade A_6 points, patients with ALBI grade 1 were less than those with ALBI grade 2 (32.7% vs 67.3%). In 31patients with Child-Pugh grade B, patients with ALBI grade 2were more than those with ALBI grade 1(12.9% vs 87.1%). There was no significant difference in the distribution of liver cancer patients with different liver function status between the two groups (P=0.851). The number of cases in individual groups was less than 5 (TACE+SORA group: ALBI grade 1_6 points, ALBI grade 1_grade B,ALBI grade 3_grade B; TACE group: ALBI grade 1_grade B, ALBI grade 3_grade B).
Table 2
Distribution of liver cancer patients with different child Pugh grades in different ALBI grades
|
Liver function status
|
TACE+SORA
|
TACE
|
Total
|
P
|
ALBI G1
|
5points
|
39
|
72
|
111
|
0.851
|
6points
|
4
|
14
|
18
|
Grade B
|
2
|
2
|
4
|
ALBI G2
|
5points
|
12
|
31
|
43
|
6points
|
14
|
23
|
37
|
Grade B
|
9
|
18
|
27
|
ALBI G3
|
Grade B
|
1
|
2
|
3
|
Total
|
81
|
162
|
243
|
|
Survival comparison between TACE + SORA group and TACE group and the influence of liver function and tumor load on survival
1. Survival comparison between TACE + SORA group and TACE group
As shown in Figure 2, for PFS and OS, the survival benefit of TACE + SORA group was significantly better than that of TACE group (median PFS 9.0 months vs 6.0 months, median OS 24.0 months vs 17.0 months; P=0.000).
2. Effects of ALBI grades on prognosis of patients receiving TACE + SORA or TACE
The prognosis of 240 patients with liver cancer of grade 1 and 2 was analyzed(Figure 3). In HCC patients with liver function of ALBI grade 1, there were significant differences in median PFS(11.2 months vs 7.0 months, P=0.000) and OS (25.5 months vs 21.0 months, P=0.006) between TACE + SORA group and TACE group. In HCC patients with liver function of ALBI grade 2, there were significant differences in median PFS(6.4 months vs 5.0months, P=0.026) and OS (18.3months vs 13.4 months,P=0.006) between TACE + SORA group and TACE group.
3. Effect of liver function on survival of patients with hepatocellular carcinoma receiving TACE + SORA or TACE
Except for the groups with less than 5 patients with different liver function status in the study population (Table 3), the survival of 232 patients (TACE+SORA group: ALBI grade 1_5 points,ALBI grade 2_5 points, ALBI grade 2_6 points, ALBI grade 2_grade B; TACE group: ALBI grade 1_5 points, ALBI grade 1_6 points, ALBI grade 2_5 points,ALBI grade 2_6 points, ALBI grade 2_grade B) with liver cancer was analyzed (Figure 4).
Patients with liver function of ALBI grade 1_5 points would had a better prognosis of TACE + SORA group than TACE group(median PFS: 11.8 months vs 7.0 months, P=0.000; median OS: 27.6 months vs 23.0 months, P=0.019). In patients with liver function of ALBI grade 2_5 points, the survival benefit of patients receiving TACE + SORA was also higher than that of patients receiving TACE only (median PFS:6.8 months vs 5.0 months, P=0.046; median OS: 22.5 months vs 15.0 months, P=0.005).
In patients with liver function of ALBI grade 2_6 points, there was no significant difference in short-term prognosis between TACE + SORA and TACE (median PFS: 6.3 months vs 5.9 months, P=0.494); however there were significant differences in long-term prognosis(median OS: 16.8 months vs 13.5 months,P=0.001). In patients with liver function of ALBI grade 2_ grade B, there was no significant difference in the prognosis of TACE + SORA group or TACE group (median PFS:6.0 months vs 4.4 months, P=0.370; median OS:14.8 months vs 11.0 months,P=0.131). Also, the median PFS and OS of patients with liver function of ALBI grade 1_6 points in TACE group were 6.5months and 15.8months.
Figure 4 showed that, compared with TACE alone, patients with good liver function (especially ALBI grade 1_5 points) could benefit from TACE + SORA. Compared with TACE alone, combined treatment strategy does not benefit patients with poor liver function (such as ALBI grade 2_grade B) significantly.
4. Effect of tumor burden layers on the prognosis of patients receiving TACE + SORA or TACE
The study analyzed the prognosis of patients with different tumor loads receiving different treatment methods(Figure 5). In patients with HCC of layer 1, there was difference in median PFS (11.9 months vs 7.8months,P=0.000) between TACE + SORA group and TACE group. But there was no significant difference in long-term survival (27.7 months vs 23.9 months, P=0.402). In patients with HCC of layer 2, TACE combined with SORA can significantly improve the prognosis of patients(median PFS: 9.8 months vs 6.0 months,median OS: 24.5 months vs 18.5 months, P<0.05). In patients with HCC of layer 3, combined treatment could also significantly improve the prognosis of patients (median PFS: 6.3 months vs 4.7 months, P=0.118; median OS: 15.0 months vs 12.2 months, P=0.001).
The results showed that the median OS of TACE + SORA or TACE was not significantly different in patients with lower tumor load (layer 1) (P=0.402); patients with high tumor load in TACE + SORA group had a long-term benefit than TACE group (P value<0.05).
5. Survival rates of patients with different liver functions and tumor loads
Table 3 showed the proportions of PFS survival in HCC patients with different liver function and tumor load treated with TACE + SORA or TACE. The 3, 6, 12 months survival rate of patients treated with combined therapy was higher than that of TACE alone in the whole population (98.8% vs 86.3%, 81.5% vs 46.8%, 23.5% vs 7.0%). Based on liver function status and tumor load layers, the PFS rate of TACE + SORA group was higher than that of TACE group.
Table 4 showed the proportions of OS survival in HCC patients with different liver function and tumor load treated with TACE + SORA or TACE. The 12, 24, 36 months survival rate of patients treated with combined therapy was higher than that of TACE alone in the whole population (92.6% vs 79.9%, 49.9% vs 20.9%, 8.9% vs 5.6%). Similarly, based on liver function status and tumor load layers, the OS rate of TACE + SORA group was higher than that of TACE group.
Table 3
The proportion of patients with liver cancer with different liver function and tumor load receiving TACE +SORA or TACE (PFS)
|
|
TACE+SORA
|
TACE
|
3-months PFS(%)
|
6-months PFS(%)
|
12-months PFS(%)
|
3-months PFS(%)
|
6-months PFS(%)
|
12-months PFS(%)
|
Overall patients
|
98.8
|
81.5
|
23.5
|
86.3
|
46.8
|
7.0
|
ALBI grades
|
1
|
100.0
|
95.6
|
37.8
|
93.1
|
55.0
|
8.2
|
|
2
|
100.0
|
65.7
|
5.7
|
83.1
|
36.6
|
2.8
|
Child-Pugh grades
|
5
|
100.0
|
90.2
|
35.3
|
88.2
|
48.4
|
4.0
|
6
|
94.4
|
77.8
|
5.6
|
86.5
|
48.6
|
8.1
|
B
|
91.7
|
50.0
|
0.0
|
77.3
|
36.4
|
9.1
|
ALBI grade1_5points
5分
|
100.0
|
94.9
|
41.0
|
93.0
|
53.3
|
5.8
|
ALBI grade1_6points
|
—
|
—
|
—
|
92.9
|
57.1
|
7.1
|
ALBI grade1_grade B
|
—
|
—
|
—
|
—
|
—
|
—
|
ALBI grade2_5points
5分
|
100.0
|
91.7
|
16.7
|
83.5
|
36.8
|
3.3
|
ALBI grade2_6points
|
100.0
|
71.4
|
0.0
|
91.3
|
43.5
|
8.7
|
ALBI grade2_grade B
|
100.0
|
55.6
|
0.0
|
72.2
|
27.8
|
0.0
|
Tumor burden layer
|
1
|
100.0
|
100.0
|
46.7
|
94.1
|
67.6
|
6.1
|
2
|
100.0
|
85.7
|
31.0
|
86.5
|
45.7
|
6.2
|
3
|
95.8
|
62.5
|
12.5
|
80.0
|
35.6
|
4.4
|
Table 4
The proportion of patients with liver cancer with different liver function and tumor load receiving TACE +SORA or TACE (OS).
Proportional hazard regression analysis of prognosis of patients with hepatocellular carcinoma and establishment of prognosis model
1. Univariate and multivariate Cox regression analysis of survival in patients with hepatocellular carcinoma
Univariate Cox regression analysis showed that treatment regimen, albumin,ALBI grade, Child-Pugh grade, tumor load layer, intrahepatic metastasis, extrahepatic metastasis and tumor thrombus were risk factors for prognosis (P<0.05; Table 5).
The factors with P<0.05 in univariate analysis were included in multivariate Cox regression analysis. The multivariate Cox regression analysis showed that treatment regimen (P=0.002), albumin level (P=0.021),ALBI grade (P=0.005), Child-Pugh grade (P=0.004), tumor load (P=0.027), intrahepatic metastasis (P=0.041) and extrahepatic metastasis (P=0.019) were independent risk factors for the prognosis of HCC patients.
Table 5
Multivariate Cox regression analysis of survival time of HCC
|
|
Univariate analysis
|
Multivariate analysis
|
|
Sig.
|
HR
|
HR(95%CI)
|
Sig.
|
HR
|
HR(95%CI)
|
|
|
|
Lower
|
Upper
|
|
|
Lower
|
Upper
|
Cure(T+S vs T)
|
0.000
|
0.527
|
0.377
|
0.737
|
0.002
|
0.577
|
0.404
|
0.822
|
Sex(Male vs Femal)
|
0.817
|
1.050
|
0.694
|
1.588
|
|
|
|
|
ALT(U/L)
|
0.231
|
1.002
|
0.999
|
1.004
|
|
|
|
|
Alb(g/L)
|
0.000
|
0.841
|
0.815
|
0.869
|
0.021
|
0.929
|
0.872
|
0.989
|
AFP(≤400 vs>400ng/ml)
|
0.127
|
0.788
|
0.580
|
1.070
|
|
|
|
|
ALBI grades(1/2/3)
|
0.001
|
0.113
|
0.033
|
0.387
|
0.005
|
0.065
|
0.009
|
0.442
|
Child-Pugh grades(5/6/B)
|
0.000
|
2.753
|
2.271
|
3.338
|
0.004
|
0.006
|
0.000
|
0.188
|
Tumor burden layers(1/2/3)
|
0.008
|
3.427
|
2.594
|
4.526
|
0.027
|
0.520
|
0.292
|
0.928
|
In-metastasis (No vs Yes)
|
0.002
|
0.529
|
0.388
|
0.720
|
0.041
|
0.696
|
0.491
|
0.985
|
Ex-metastasis (No vs Yes)
|
0.000
|
6.125
|
4.136
|
9.070
|
0.019
|
0.471
|
0.251
|
0.884
|
Tumor thrombus(No vs Yes)
|
0.034
|
2.427
|
1.781
|
3.307
|
0.313
|
0.825
|
0.568
|
1.198
|
2. Establishment of nomogram prediction model for prognosis and risk stratification based on the model
2.1 Establishment of nomogram prediction model
In order to further predict the prognosis of patients in TACE + SORA group and TACE group according to liver function and tumor load, we tried to establish linear prediction model and risk stratification model. Through univariate and multivariate Cox regression analysis, seven independent prognostic risk factors were identified, including treatment regimen, albumin level, ALBI grade, Child-Pugh grade, tumor load layer, intrahepatic metastasis and extrahepatic metastasis. Based on Cox regression analysis, seven variables were included in nomogram function at the same time, and nomogram diagram was constructed (Figure 6). The prognosis of patients was predicted after scoring by constructing nomograms. We could more intuitively reflect the predictive effect of different variables on the prognosis of patients in Cox regression analysis.The higher the total score was, the better the patients' expected survival was.
According to the nomogram, the clinical indicators of patients were scored, and the total score was 180. The specific scores were as follows (unit: points) : TACE+SORA 25 points, TACE 0 points, the score of patients with and without intrahepatic metastasis were 0 and 13 respectively, the score of patients with and without extrahepatic metastasiswere 0 and 16 respectively, patients with ALBI grade 1,2 and grade 3were 0, 4 and 8 points, patients with Child Pugh 5 points,6 points and grade B were 11.3points,0 and 0.5point, tumor load layer1, 2, 3were 32.5 points, 13.5 points and 0 point, The albumin score was as follows: (Alb-20)*(100/35).
2.2 Validation of nomogram prediction model
2.2.1 C-index and AUC of nomogram prediction model
Further, the prediction performance of the prediction model was evaluated and verified. The C-index and AUC of the nomogram model based on ALBI grade and Child-Pugh grade were 0.846 (95%CI: 0.762-0.913) and 0.703(0.631-0.776), respectively. The C-index and AUC of the nomogram model based on ALBI grade or Child-Pugh grade separately were 0.844(0.741-0.915) and 0.839(0.768-0.898). It also performed well in predicting 1-year OS, 2-year OS and 3-year OS,The AUC values of ROC curve were as follows: 0.928(0.891-0.965) , 0.792(0.735-0.849) and 0.707(0.635-0.779)(Figure 7).
The results of Figure 7 showed that the nomogram model combined with ALBI grade and Child-Pugh grade had a good predictive ability for the prognosis of patients with unresectable or recurrent liver cancer after TACE combined with sorafenib or TACE alone, which was better than the nomogram model only included in ALBI grade or Child-Pugh grade; At the same time, the model based on ALBI grade was better than that based on Child-Pugh grade.
ROC curves at different time points show that, the performance of ALBI grade in predicting prognosis was better than Child-Pugh grade(AUC of 1 year: .787 vs 0.774; AUC of 2 years: 0.683 vs 0.668; AUC of 3 years : 0.631 vs 0.582; Figure 9). The ability to predict the survival of patients within 2 years was similar; with the extension of time, ALBI was more effective in predicting the survival (ΔAUC of 1, 2,3-years OS were 0.013,0.015 and 0.049, Figure 8).
2.2.2 Construction of calibration curve of nomogram prediction model
The calibration curves of linear prediction model were further constructed by using RStudio software (Figure 9). According to the calibration curves, the nomogram model could well fit the predicted values when predicting OS of 12 months, 18 months, 24 months and 36 months. The actual predicted value fluctuated around the predicted value in a certain range. This showed the good prediction ability of the model again. The calibration curves showed that within the study population, the prognostic nomogram model could accurately predict OS at 12, 18, 24 and 36 months.
2.3 Nomogram prediction model could be used for risk stratification of prognosis of liver cancer patients
According to the nomogram, seven clinical variables were scored and summed. The total score of the linear prediction model was 180 points, and the score range of 243 patients was 34.1-174.7 points. Using X-tile software to determine the optimal cut-off value (62.3 points and 104.2points, Figure 10B),divide patients into three risk groups. According to the score, the risk stratification of patients was < 62 (high risk), ≥ 62 and < 104 (medium risk), ≥ 104 and < 180 (low risk). The corresponding survival curves of different risk groups were shown in Figure 10B. When the patient's score was ≥ 104 points and < 180 points, the survival curve was like pink curve in Figure 10; When the patient's score was <62 points, the survival curve like blue curve in Figure 10.
The survival curve based on risk stratification (Figure 10) shows that, the higher the total score was, the longer the expected survival was (low-risk group); the higher the total score was, the shorter the expected survival was (high-risk group) .