4.1. Demographics and tumor characteristics
The present study included 347 patients, which 83.6% of patients (290 patients) were male and 16.4% were female (57 patients). Among the 347 patients, 300 patients’ (86.5%) hepatitis B surface antigen (HBsAg) were positive, 222 patients (64%) had liver cirrhosis, 188 patients (54.2%) received postoperative local ablation therapy or TACE, and 33 patients (9.5%) used multiple kinase inhibitors (MKI). In terms of tumor parameters, 304 patients (87.6%) had a tumor capsule; 305 patients (87.9%) had a single tumor; 168 patients’ (48.4%) tumor size was > 5 cm; 172 patients (49.6%) had microvascular invasion (MVI); and 29 (8.4%), 283 (81.6%), and 35 (10.1%) patients had well, moderately and poorly differentiated tumor cells, respectively. To facilitate the analysis, 9 patients with BCLC stage 0 were merged into the BCLC stage A group, and 5 patients with AJCC stage IV were merged into the AJCC stage III group (Table 1).
Table 1
Baseline characteristics of the 347 patients
Variable
|
|
Cases (%)
/Median (P25-P75)a
|
Variable
|
|
Cases (%)
/Median (P25-P75)
|
Sex
|
Male
|
290 (83.6)
|
Child class
|
A
|
335 (96.5)
|
|
Female
|
57 (16.4)
|
|
B
|
12 (3.5)
|
Age (years)
|
≥60
|
95 (27.4)
|
Performance status
|
1
|
144 (41.3)
|
|
< 60
|
252 (72.6)
|
|
0
|
203 (58.2)
|
HBsAg
|
Positive
|
300 (86.5)
|
BCLC staging
|
Ad
|
175 (50.4)
|
|
Negative
|
47 (13.5)
|
|
B
|
17 (4.9)
|
Liver cirrhosis
|
Yes
|
222 (64)
|
|
C
|
155 (44.7)
|
|
No
|
125 (36)
|
AJCC staging
|
Ⅰ
|
166 (47.8)
|
Portal vein invasion
|
Yes
|
21 (6.1)
|
|
Ⅱ
|
136 (39.2)
|
|
No
|
326 (93.9)
|
|
Ⅲe
|
45 (13)
|
Ascites
|
Yes
|
11 (3.2)
|
CNLC staging
|
Ⅰ
|
291 (83.9)
|
|
No
|
336 (96.8)
|
|
Ⅱ
|
33 (9.5)
|
Ablation or TACE
|
Yes
|
188 (54.2)
|
|
Ⅲ
|
23 (6.6)
|
|
No
|
159 (45.8)
|
Neutrophil count
|
|
3.6 (2.56–4.72)
|
AFP (ng/ml)
|
> 400
|
132 (38)
|
Lymphocyte count
|
|
1.45 (1.03–1.84)
|
|
≤400
|
215 (62)
|
Platelet count
|
|
150 (109–205)
|
Tumor capsule
|
Yes
|
304 (87.6)
|
γ-GT
|
|
60 (36–115)
|
|
No
|
43 (12.4)
|
AST
|
|
37 (27–53)
|
Tumor number
|
≥2
|
42 (12.1)
|
Total bilirubin
|
|
15.2 (12-21.2)
|
|
1
|
305 (87.9)
|
Prothrombin time
|
|
11.7 (11-12.4)
|
Tumor size (cm)
|
>5
|
168 (48.4)
|
Albumin
|
|
44 (40.7–46.7)
|
|
≤5
|
179 (51.6)
|
Fibrinogen
|
|
2.61 (2.13–3.27)
|
MVIb
|
Yes
|
172 (49.6)
|
NLR
|
|
2.44 (1.71–3.37)
|
|
No
|
175 (50.4)
|
PLR
|
|
104.46(74.69-150.52)
|
Cell differentiation
|
Poor
|
35 (10.1)
|
GPR
|
|
0.45 (0.22–0.88)
|
|
Moderate
|
283 (81.6)
|
ALR
|
|
26.97 (17.02–44.74)
|
|
Well
|
29 (8.4)
|
FAR
|
|
0.06 (0.05–0.08)
|
MKIc
|
Yes
|
33 (9.5)
|
|
|
|
|
No
|
314 (90.5)
|
|
|
|
a. P25-P75: (lower quartile - upper quartile) b. MVI: microvascular invasion c: MKI: multiple kinase inhibitor |
d. Nine patients with BCLC stage 0 were merged into the BCLC stage A group |
e. Five patients with AJCC stage IV were merged into the AJCC stage III group |
4.2. Optimal cutoff values of the inflammatory markers
The optimal cutoff values of NLR, PLR, GPR, ALR, and FAR were 2.33, 117.09, 0.48, 31, and 0.06, respectively. The areas under the curve (AUC) of those markers were 0.569, 0.553, 0.680, 0.647, and 0.632 (Fig. 2).
4.3. OS and DFS rates
The median follow-up time was 45 months. During the follow-ups, 216 patients (62.2%) experienced recurrence, and 147 patients (42.4%) passed away. The 1-, 3-, and 5-year DFS rates were 69.8%, 41.5%, and 30.8%, and the 1-, 3-, and 5-year OS rates were 84.6%, 59.4%, and 52.2%, respectively. For the DFS rates, the P values of different groups of NLR, PLR, GPR, ALR and FAR values were 0.088, 0.082, < 0.001, < 0.001 and < 0.001, respectively. For the OS rate, the P values of different groups of NLR, PLR, GPR, ALR and FAR were 0.004, 0.008, < 0.001, < 0.001 and < 0.001, respectively (Fig. 3).
4.4. Independent prognostic factors for DFS and OS rates
In the univariate analysis for DFS, variables with P < 0.1 including sex, age, portal vein invasion, ascites, AFP, tumor capsule, tumor size, tumor number, MVI, cell differentiation, MKI, NLR, PLR, GPR, ALR and FAR were selected for the multivariate analysis. According to the results of multivariate analysis, AFP > 400 ng/ml, ALR > 31, GPR > 0.48, MVI, absence of tumor capsule, and tumor size>5cm were independent prognostic factors for DFS (Table 2).
Table 2
Univariate and multivariate analysis of DFS
Variable
|
|
Univariate analysis
|
Multivariate analysis
|
|
HR (95%CI)
|
P value
|
HR (95%CI)
|
P value
|
Sex
|
Male/Female
|
1.401 (0.947–2.072)
|
0.092
|
|
|
Age (years)
|
≥60/<60
|
0.767 (0.560–1.050)
|
0.097
|
|
|
HBsAg
|
Positive/Negative
|
1.317 (0.874–1.985)
|
0.189
|
|
|
Liver cirrhosis
|
Yes/No
|
1.080 (0.819–1.425)
|
0.586
|
|
|
Portal vein invasion
|
Yes/No
|
2.710 (1.686–4.357)
|
< 0.001
|
|
|
Ascites
|
Yes/No
|
2.997 (1.621–5.539)
|
< 0.001
|
|
|
Ablation or TACE
|
Yes/No
|
0.833 (0.637–1.089)
|
0.181
|
|
|
AFP (ng/ml)
|
> 400/≤400
|
1.674 (1.279–2.190)
|
< 0.001
|
1.377 (1.042–1.821)
|
0.025
|
Tumor capsule
|
No/Yes
|
2.146 (1.497–3.078)
|
< 0.001
|
1.864 (1.289–2.697)
|
0.001
|
Tumor number
|
≥2/1
|
1.481 (1.006–2.181)
|
0.047
|
|
|
Tumor size (cm)
|
>5/≤5
|
1.766 (1.348–2.313)
|
< 0.001
|
1.440 (1.090–1.901)
|
0.010
|
MVIa
|
Yes/No
|
1.649 (1.259–2.160)
|
< 0.001
|
1.352 (1.022–1.789)
|
0.035
|
Cell differentiation
|
Moderate/Well
|
2.248 (1.189–4.251)
|
0.013
|
|
|
Poor/Well
|
2.869 (1.364–6.031)
|
0.005
|
|
|
MKIb
|
No/Yes
|
0.596 (0.395-0.900)
|
0.014
|
|
|
Child class
|
B/A
|
1.683 (0.891–3.177)
|
0.109
|
|
|
NLR
|
> 2.33/≤2.33
|
1.264 (0.965–1.656)
|
0.089
|
|
|
PLR
|
> 117.09/≤117.09
|
1.270 (0.969–1.664)
|
0.083
|
|
|
GPR
|
> 0.48/≤0.48
|
2.327 (1.770–3.060)
|
< 0.001
|
1.931 (1.445–2.581)
|
< 0.001
|
ALR
|
> 31/≤31
|
1.903 (1.455–2.488)
|
< 0.001
|
1.438 (1.083–1.910)
|
0.012
|
FAR
|
> 0.06/≤0.06
|
1.617 (1.234–2.117)
|
< 0.001
|
|
|
a: MVI: microvascular invasion b: MKI: multiple kinase inhibitor |
In the univariate analysis for OS, variables with P < 0.1 including HBV, portal vein invasion, ascites, AFP, tumor capsule, tumor size, tumor number, MVI, cell differentiation, NLR, PLR, GPR, ALR and FAR, were then entered into the Cox multivariate analysis. The results demonstrated that PLR > 117.09, ALR > 31, GPR > 0.48, MVI, absence of tumor capsule, and tumor size>5cm were independent prognostic factors for OS (Table 3).
Table 3
Univariate and multivariate analysis of OS
Variable
|
|
Univariate analysis
|
Multivariate analysis
|
|
HR (95%CI)
|
P value
|
HR (95%CI)
|
P value
|
Sex
|
Male/Female
|
1.478 (0.912–2.393)
|
0.113
|
|
|
Age (years)
|
≥60/<60
|
0.670 (0.451–0.996)
|
0.048
|
|
|
HBsAg
|
Positive/Negative
|
1.630 (0.940–2.736)
|
0.083
|
|
|
Liver cirrhosis
|
Yes/No
|
1.036 (0.742–1.448)
|
0.834
|
|
|
Portal vein invasion
|
Yes/No
|
3.757 (2.256–6.256)
|
< 0.001
|
|
|
Ascites
|
Yes/No
|
3.564 (1.804–7.039)
|
< 0.001
|
|
|
Ablation or TACE
|
Yes/No
|
1.047 (0.757–1.449)
|
0.780
|
|
|
AFP (ng/ml)
|
> 400/≤400
|
1.844 (1.334–2.549)
|
< 0.001
|
|
|
Tumor capsule
|
No/Yes
|
2.710 (1.811–4.056)
|
< 0.001
|
2.119 (1.394–3.222)
|
< 0.001
|
Tumor number
|
≥2/1
|
1.893 (1.230–2.915)
|
0.004
|
|
|
Tumor size (cm)
|
>5/≤5
|
2.170 (1.553–3.033)
|
< 0.001
|
1.646 (1.157–2.342)
|
0.006
|
MVIa
|
Yes/No
|
2.084 (1.495–2.906)
|
< 0.001
|
1.678 (1.195–2.355)
|
0.003
|
Cell differentiation
|
Moderate/well
|
2.718 (1.112–6.647)
|
0.028
|
|
|
Poor/well
|
3.703 (1.375–9.976)
|
0.010
|
|
|
MKIb
|
No/Yes
|
0.853 (0.492–1.480)
|
0.572
|
|
|
Child class
|
B/A
|
1.380 (0.609–3.126)
|
0.440
|
|
|
NLR
|
> 2.37/≤2.37
|
1.615 (1.157–2.253)
|
0.005
|
|
|
PLR
|
> 117.09/≤117.09
|
1.542 (1.115–2.133)
|
0.009
|
1.465 (1.024–2.096)
|
0.037
|
GPR
|
> 0.48/≤0.48
|
3.002 (2.131–4.230)
|
< 0.001
|
2.554 (1.757–3.712)
|
< 0.001
|
ALR
|
> 31/≤31
|
2.255 (1.628–3.122)
|
< 0.001
|
1.553 (1.098–2.197)
|
0.013
|
FAR
|
> 0.06/≤0.06
|
1.984 (1.425–2.764)
|
< 0.001
|
|
|
a: MVI: microvascular invasion b: MKI: multiple kinase inhibitor |
4.5. The relationship between independent prognostic inflammatory markers and clinicopathological features
As ALR, GPR were independent prognostic factors for DFS; ALR, GPR, and PLR were independent prognostic factors for OS, the correlations of commonly clinicopathological variables with different groups of ALR, GPR and PLR were computed by the χ2 test, respectively. The results presented that ALR was correlated with portal vein invasion, ascites, AFP, tumor number, tumor size, MVI, BCLC, AJCC and CNLC staging; GPR was correlated with sex, portal vein invasion, ascites, AFP, tumor capsule, tumor number, tumor size, MVI, BCLC, AJCC and CNLC staging; PLR was correlated with sex, postoperative ablation or TACE, tumor capsule, tumor size, and BCLC staging [see Additional file 1].
4.6. Creation and comparison of inflammatory scoring models for DFS and OS
We generated models for the ALR, GPR and PLR score. At the same time, we created an ALR-GPR score model for DFS, and the ALR-GPR, ALR-PLR, GPR-PLR and ALR-GPR-PLR score (A-G-P score) models for OS. For simplicity of calculation, ALR≤31, GPR≤0.48 and PLR≤117.09 were defined as a score of 0, and ALR>31, GPR>0.48 and PLR>117.09 were defined as a score of 1.
ALR, GPR and PLR score models consisted of scores of 0 and 1, models for the ALR-GPR, ALR-PLR and GPR-PLR score consisted of scores of ≤1 and 2; the A-G-P score model consisted of scores of ≤1, 2 and 3 (Table 4).
Table 4
Models of inflammatory markers
Model for OS
|
|
Score
|
Model for DFS
|
|
Score
|
ALR score
|
ALR > 31
|
1
|
ALR score
|
ALR > 31
|
1
|
|
ALR≤31
|
0
|
|
ALR≤31
|
0
|
GPR score
|
GPR > 0.48
|
1
|
GPR score
|
GPR > 0.48
|
1
|
|
GPR≤0.48
|
0
|
|
GPR≤0.48
|
0
|
PLR score
|
PLR > 117.09
|
1
|
ALR-GPR score
|
ALR > 31 and GPR > 0.48
|
2
|
|
PLR≤117.09
|
0
|
|
others
|
≤1
|
ALR-GPR score
|
ALR > 31 and GPR > 0.48
|
2
|
|
|
|
|
others
|
≤1
|
|
|
|
ALR-PLR score
|
ALR > 31 and PLR > 117.09
|
2
|
|
|
|
|
others
|
≤1
|
|
|
|
GPR-PLR score
|
GPR > 0.48 and PLR > 117.09
|
2
|
|
|
|
|
others
|
≤1
|
|
|
|
A-G-P scorea
|
ALR > 31, GPR > 0.48
and PLR > 117.09
|
3
|
|
|
|
|
ALR > 31, GPR > 0.48,
and PLR≤117.09;
or ALR > 31, PLR > 117.09,
and GPR≤0.48;
or PLR > 117.09, GPR > 0.48,
and ALR≤31
|
2
|
|
|
|
|
others
|
≤1
|
|
|
|
a: A-G-P score: ALR-GPR-PLR score |
We further verified whether the above combined scoring models (as categorical variables) w ere independent predictors of prognosis through univariate and multivariate analysis. Obviously, single inflammatory marker model ALR, GPR and PLR score were independent predictors. When verifying one combined scoring model, we no longer put the individual inflammatory markers that make up the model into analyses in order to exclude their interactions. The results demonstrated that all of the ALR-GPR, ALR-PLR, GPR-PLR and A-G-P score models were independent predictors for OS, and the ALR-GPR score was an independent predictive factor for DFS [see Additional file 2–6].
Then, we compared the single inflammatory marker models and the combined scoring models by the C-index, AIC and likelihood ratio. For OS, compared with other models, the A-G-P score model had the smallest AIC value (1569.94), the largest C-index value (0.653, 95% CI: 0.610–0.696) and the largest likelihood ratio (50.48), suggesting that the A-G-P score model has a better prediction accuracy, goodness of fit, and uniformity in predicting the survival of patients who underwent resection. In terms of DFS, the single inflammatory marker model, the GPR score had the smallest AIC value (2264.32), the largest C-index value (0.605, 95%CI: 0.572–0.638) and largest likelihood ratio (37.39), suggesting that the GPR score model has a better prediction (Table 5).
Table 5
Comparison of models for OS/DFS
Model
|
AIC
|
C-index
|
Likelihood ratio
|
OS
|
|
|
|
ALR score
|
1596.53
|
0.61
|
23.89
|
GPR score
|
1578.07
|
0.64
|
42.34
|
PLR score
|
1613.68
|
0.56
|
6.73
|
ALR-GPR score
|
1574.60
|
0.61
|
30.96
|
ALR-PLR score
|
1598.11
|
0.59
|
22.30
|
GPR-PLR score
|
1585.71
|
0.60
|
34.70
|
A-G-P scorea
|
1569.94
|
0.65
|
50.48
|
DFS
|
|
|
|
ALR score
|
2280.10
|
0.59
|
21.61
|
GPR score
|
2264.32
|
0.61
|
37.39
|
ALR-GPR score
|
2272.73
|
0.59
|
28.98
|
a: A-G-P score: ALR-GPR-PLR score |
4.7. The effect of stratification of the A-G-P score model and GPR score model on different stages of HCC
The A-G-P score model has a good discriminating ability for OS in the whole population with a statistically significant difference (score ≤1/score 2, P < 0.001; score ≤1/ score 3, P < 0.001; score 2/score 3, P < 0.001). For the pairwise comparison of the A-G-P score, we used Bonferroni correction, and P < 0.0167 (α = 0.05/3) was considered statistically significant between the different scores. Regarding BCLC staging, score ≤1/score 3 (P = 0.007) in stage A; score ≤1/score 2 (P = 0.004), score ≤1/score 3 (P < 0.001), and score 2/score 3 (P < 0.004) in stage C were significantly different. Regarding AJCC staging, score ≤1/score 3 (P < 0.001) and score 2/score 3 (P = 0.006) in stage I, score ≤1/score 2 (P < 0.001) and score ≤1/score 3 (P < 0.001) in stage II were significantly different. For CNLC staging, score ≤1/score 2 (P < 0.001), score ≤1/score 3 (P < 0.001), and score 2/score 3 (P < 0.001) in stage Ⅰ showed significant differences for the prognosis of OS (Fig. 4).
The GPR score model was a good differentiator for DFS in the whole population, with P < 0.001. For different stages of HCC, in BCLC stage C (P < 0.001), AJCC stage II (P < 0.001), AJCC stage III (P = 0.004), CNLC stage I (P < 0.001), and CNLC stage III (P = 0.027), different scores of GPR score model showed significant differences for DFS (Fig. 5).