A sample model established by S-index predicting overall survival after curative resection of primary hepatocellular carcinoma

Abstract

Purpose

Prognostic prediction after curative resection of primary hepatocellular carcinoma (PHCC) remains an arduous task. The S-index calculated from γ-glutamyl transpeptidase, albumin, and platelets is reported to predict the severity of liver fibrosis. We constructed a nomogram for predicting the survival probability of PHCC based on a new indicator, the S-index, combined with other routine clinical parameters.

Patients and methods

We selected 490 patients with PHCC postradical surgery at the First Affiliated Hospital of Wenzhou Medical University between January 2007 and January 2014. The subjects were randomly allocated into the training cohort and the validation cohort in the ratio 7:3 by the digital method. Important variables screened by univariate analysis were included in multivariate analysis to obtain independent risk factors for predicting the prognosis of PHCC. The construction of the nomogram was based on Cox proportional hazard regression models. The concordance index (C-index) was used in the nomogram for evaluating the model performance for prognosis. We drew time-dependent receiver operating characteristic curves to compare our model with other staging systems.

Results

The nomogram based on six independent risk factors after multivariate analyses had good predictive power after radical surgery of PHCC. In the training and validation groups, the C-index of the nomogram was highly consistent for evaluating survival from PHCC. Compared with the traditional scoring system, the areas under time-dependent receiver operating characteristic curves were 0.7382, 0.7293, and 0.7520 for 1-, 3-, and 5-year overall survival, respectively. In summary, the nomogram showed excellent results in terms of prognosis of PHCC.

Conclusion

Based on the S-index and the other clinical indicators, we developed a precise nomogram that predicts the survival probability of patients with PHCC after radical surgery. This tool can provide effective information for surgeons and patients.

Keywords:

• PHCC
• radical resection
• nomogram
• prognosis

Introduction

Primary hepatocellular carcinoma (PHCC) is the most common hepatic cancer, which is the second largest cause of cancer-related death worldwide.¹ Early diagnosis of hepatocellular carcinoma is important to achieve long-term survival. Unfortunately, <30% of patients with hepatocellular carcinoma achieve early diagnosis and treatment.² Partial hepatectomy remains the most effective intervention for patients with hepatocellular carcinoma. The recurrence rates and the 5-year survival probability remain poor, even after curative resection. At present, the internationally accepted criteria for evaluating the outcomes of liver cancer are various staging paradigms, the most common of which is the seventh edition of the TNM (TNM seventh) system and the Barcelona Clinic Liver Cancer (BCLC) staging system. The Okuda staging system and the Cancer of the Liver Italian Program (CLIP) staging system also play important roles in the assessment of the prognosis of liver cancer.³ Unfortunately, these standards are complicated and have some limitations. We attempted to generate a simpler and more clinically relevant way to evaluate the prognosis of PHCC.

It is well known that cirrhosis and hepatic fibrosis are closely related to primary liver cancer. In recent years, many investigators had found that the models for evaluating and diagnosing hepatic fibrosis had a significant impact on the prognosis of liver cancer. Furthermore, Pang et al found that gamma-glutamyl transpeptidase and platelets (PLTs) ratio was an important factor influencing hepatocellular carcinoma associated with hepatitis B virus infections.⁴ Choi et al found that the Forns index was an important influencing factor for determining recurrence and death after curative resection of liver cancer.⁵ Zhou et al found that the S-index obtained by laboratory indicators accurately evaluated the degree of liver cirrhosis and fibrosis for patients with hepatitis B virus infections.⁶ Nevertheless, the relationship between S-index and prognosis after liver cancer resection remains unknown. In recent years, nomograms have shown excellent performance in predicting various tumor prognoses, including hepatocellular carcinoma,⁷ colon cancer,⁸ gastrointestinal stromal tumors, and others.⁹ Here, we establish a nomogram to study the correlation of survival between the S-index and prognosis after curative resection of PHCC.

Patients and methods

Research design and patients

We obtained clinical and laboratory data from 615 patients with liver cancer who were treated in the First Affiliated Hospital of Wenzhou Medical University between January 2007 and January 2014. We developed the following strict exclusion criteria: 1) patients who were diagnosed with nonprimary hepatocellular carcinoma according to clinical data and postoperative pathological findings; 2) patients without complete clinical data; 3) patients who received other treatments before radical surgery, including interventional therapy, radiofrequency ablation, chemical drugs, and radiation therapy; 4) patients with liver cancer secondary to other malignant tumors; 5) patients with other terminal cancers; 6) patients who cannot undergo curative resection, such as those with distant metastasis; 7) patients with basic diseases, such as those with severe vital organ damage; and 8) patients who died due to other noncancer-related factors such as accidents. Based on the exclusion criteria, we finally selected 490 patients with PHCC. All subjects were randomly allocated into the training set and validation set in a 7:3 ratio by the digital method. The training group consisted of 343 patients for establishing a nomogram, and the remaining 147 patients were placed in the validation group for confirming the predicting power of the nomogram. Our study was supported by the Ethical Committee and Institutional Review Board of the First Affiliated Hospital of Wenzhou Medical University. All patients signed informed consent forms, and all study procedures were conducted in accordance with the latest revision of the World Medical Association Declaration of Helsinki.

Data collection and follow-up

All patients diagnosed with PHCC had complete clinical and laboratory variables collected. The clinical variables were as follows: gender, age, history of alcohol abuse, body mass index, tumor number, tumor size, tumor capsule, cirrhosis, ascites, portal vein tumor thrombus (PVTT), microvascular invasion, lymph node metastasis, intrahepatic metastases, invasion of adjacent tissues, and tumor grade. The laboratory data were as follows: white blood cell (WBC), PLT, total bilirubin (TB), albumin (ALB), fibrinogen, aspartate transaminase (AST), alanine aminotransferase (ALT), alpha-fetoprotein (AFP), prothrombin time (PT), and γ-glutamyl transpeptidase (γ-GT). AFP was determined by UniCel DxI-800 analyzer. TB and hepatic function index (ALT, AST, and γ-GT) were tested by Beckman AU5800 analyzer. Blood clotting function (fibrinogen, PT) was measured on a Diagnostica Stago analyzer. Complete blood count, including WBC and PLT, was tested using a Sysmex XE-2100 analyzer.

Based on these variables, we evaluated the staging system for hepatocellular carcinoma, including Child–Pugh classification, BCLC staging system, TNM staging system, the CLIP, and the Okuda staging system.¹⁰ The S-index was calculated from laboratory variables (γ-GT, ALB, and PLTs) with the formula: 1,000× γ-GT/(PLT × ALB).²⁴

All patients were followed up by telephone or outpatient visit every 3 months for the first 2 years, every 6 months after 2 years, and once in every year after 5 years.

Statistical analyses

Continuous variables that conformed to a normal distribution were expressed as mean $\pm \text{SD}(\overline{\mathrm{x}}\pm \text{SD})$; otherwise, they were identified by quartile (P25, P75). Numbers (%) were used to identify categorical variables. Continuous variables were tested by the Student’s t-test or Whitney U test for variables with an abnormal distribution. The chi-squared test or Fisher’s exact test was used to evaluate categorical variables. Overall survival (OS) was calculated by the Kaplan–Meier method and log-rank method, and was defined as the length of time from the date of diagnosis to the date of death or the end of follow-up. Statistically significant variables in univariate analysis (P<0.001) were included in the multivariate analysis. The nomogram that was used to assess the 1-, 3-, and 5-year survival probability was built on variables with P<0.05 in the multivariate analysis. The concordance index (C-index) was used in the nomogram for evaluating the model performance for prognosis in the training and validation sets. The time-dependent receiver operating characteristics (tdROC) curve was used to compare the performance of the nomogram, S-index, BCLC, Okuda, and CLIP for the prognosis of PHCC. Larger C-index and area under the time-dependent receiver operating characteristic (tdAUC) curve indicated higher reliability of prediction and diagnosis. All statistical analyses were carried out with SPSS 22 and RStudio software (version 1.1.447- 2009–2018; RStudio, Inc.)

Results

Baseline characteristics

All patients were diagnosed with PHCC, of which 285 (83.1%) were male in the training group and 128 (87.1%) were male in the validation group. The median survival times were 58.23 and 44.53 months, respectively. Analysis of the baseline characteristics of the 490 patients (343 from the training analysis and 147 from the validation) showed no statistical differences among variables (P>0.05) between the training and validation groups. In the training cohort, the 1-, 3-, and 5-year OS rates were 82.8%, 63.8%, and 47.8%, respectively, corresponding to 77.6%, 58.5%, and 38.8%, respectively, in the validation cohort (Table 1).

Table 1 Baseline demographics and clinical characteristics of patients in the training cohort and the validation cohort

Variables	Training cohort (n=343)	Validation cohort (n=147)	P-value
Age (years)	56.2±11.4	55.9±11.0	0.751
Height (cm)	168.0 (162.0, 170.0)	169.0 (164.0, 172.0)	0.234
Weight (kg)	62.0 (57.0, 69.5)	64.0 (58.0, 70.0)	0.137
BMI (kg/m^2)	22.4 (20.6, 24.5)	23.0 (20.9, 24.5)	0.249
Tumor size (cm)	4.0 (3.0, 6.5)	4.5 (3.0, 7.0)	0.397
AFP (ng/mL)	48.0 (5.7, 348.0)	30.4 (5.5, 448.6)	0.533
PT (seconds)	13.9 (13.3, 14.7)	13.8 (13.2, 14.8)	0.203
Fibrinogen (g/L)	2.8 (2.4, 3.5)	2.9 (2.4, 3.6)	0.782
PLT (×10^9/L)	143.0 (92.0, 175.0)	150.0 (98.0, 192.0)	0.200
ALB (g/L)	40.2 (36.5, 43.6)	40.1 (36.7, 43.7)	0.986
TBil (μmol/L)	12.0 (9.0, 17.0)	11.0 (8.0, 15.0)	0.236
TC (mmol/L)	4.2 (3.5, 4.7)	4.3 (3.6, 4.9)	0.301
ALT (U/L)	37.0 (25.0, 55.0)	36.0 (25.0, 57.0)	0.893
AST(U/L)	44.0 (30.0, 89.0)	41.0 (29.0, 72.0)	0.334
γ-GT(U/L)	52.0 (32.0, 106.0)	65.0 (40.0, 115.0)	0.092
Sex			0.267
Male	285 (83.1)	128 (87.1)
Female	58 (16.9)	19 (12.9)
History of alcohol abuse			0.511
No	192 (56.0)	87 (59.2)
Yes	151 (44.0)	60 (40.8)
Cirrhosis			0.218
No	102 (29.7)	52 (35.4)
Yes	241 (70.3)	95 (64.6)
Ascites			0.474
No	289 (84.3)	120 (81.6)
Yes	54 (15.7)	27 (18.4)
Tumor grade			0.081
3/4	240 (70.0)	91 (61.9)
1/2	103 (30.0)	56 (38.1)
Tumor number			0.419
Single	294 (85.7)	130 (88.4)
Multiple	49 (14.3)	17 (11.6)
Tumor capsule			0.355
No	285 (83.1)	117 (79.6)
Yes	58 (16.9)	30 (20.4)
PVTT			0.628
No	332 (96.8)	141 (95.9)
Yes	11 (3.2)	6 (4.1)
Invasion of adjacent tissues			0.328
No	329 (95.9)	138 (93.9)
Yes	14 (4.1)	9 (6.1)
MVI			0.536
No	301 (87.8)	126 (85.7)
Yes	42 (12.2)	21 (14.3)
Intrahepatic transfer			0.963
No	263 (76.7)	113 (76.9)
Yes	80 (23.3)	34 (23.1)
HBsAg			0.764
Negative	67 (19.5)	27 (18.4)
Positive	276 (80.5)	120 (81.6)
Antiviral therapy			0.140
No	190 (55.4)	92 (62.6)
Yes	153 (44.6)	55 (37.4)
Anti-HCV (IgG)			0.594
Negative	334 (97.4)	145 (98.6)
Positive	9 (2.6)	2 (1.4)
AFP (ng/mL)			0.448
<400	263 (76.7)	108 (73.5)
≥400	80 (23.3)	39 (26.5)
TNM stage			0.188
I/II	298 (86.9)	121 (82.3)
III/IV	45 (13.1)	26 (17.7)
Child–Pugh class			0.784
A	240 (70.0)	102 (69.4)
B	97 (28.3)	41 (27.9)
C	6 (1.7)	4 (2.7)
CLIP			0.266
0	120 (35.0)	39 (26.5)
1	71 (20.7)	40 (27.2)
2	77 (22.4)	31 (21.1)
3	44 (12.8)	26 (17.7)
4	27 (7.9)	9 (6.1)
5	4 (1.2)	2 (1.4)
Okuda			0.556
0	178 (51.9)	66 (44.9)
1	122 (35.6)	59 (40.1)
2	34 (9.9)	17 (11.6)
3	9 (2.6)	5 (3.4)
BCLC			0.494
A	245 (71.4)	106 (72.1)
B	47 (13.7)	14 (9.5)
C	45 (13.1)	23 (15.6)
D	6 (1.7)	4 (2.7)
S-indexa			0.797
<0.28	163 (47.5)	68 (46.3)
≥0.28	180 (52.5)	79 (53.7)
1-year survival probability (n, %)	282 (82.8)	114 (77.6)	0.229
3-year survival probability (n, %)	219 (63.8)	86 (58.5)	0.263
5-year survival probability (n, %)	164 (47.8)	57 (38.8)	0.065

Note:

a The S-index calculated from γ-GT, ALB, and PLTs is a mathematical model for predicting the severity of liver fibrosis. Continuous variables with an abnormal distribution were identified by quartile (P25, P75).

Abbreviations: γ-GT, γ-glutamyl transpeptidase; AFP, alpha-fetoprotein; ALB, albumin; ALT, alanine aminotransferase; AST, aspartate aminotransferase; BCLC, Barcelona Clinic Liver Cancer; BMI, body mass index; CLIP, Cancer of the Liver Italian Program; HCV, hepatitis C virus; MVI, microvascular invasion; PLT, platelet; PT, prothrombin time; PVTT, portal vein tumor thrombus; TBil, total bilirubin; TC, total cholesterol.

Prognostic factors of the training cohort

The optimal cutoff value of S-index, corresponding to the maximum sensitivity and specificity, was 0.28 according to the Youden’s index (Table S1; Figure S1). Patients were stratified into two groups by the optimal cutoff values of S-index (<0.28 and ≥0.28). The results of the univariate analysis are shown in Table 2. Variables of P-values <0.001 in the univariate analysis were included in multivariate analysis. As shown in Table 3, only six variables were independent risk factors for influencing the survival probability of PHCC: ALB, TB, PVTT, tumor grade, tumor size, and S-index (P<0.05).

Table 2 Univariate analysis of the training cohort

Variables	β	SE	P-value
Statistically nonsignificant factors
Age (years)	−0.00214	0.00682	0.7537
Sex	−0.01127	0.20379	0.9559
Height (cm)	0.02084	0.01149	0.0699
Weight (kg)	0.00349	0.00757	0.6449
BMI (kg/m^2)	−0.00873	0.02468	0.7235
PLT (×10^9/L)	−0.00220	0.00123	0.0735
AST (U/L)	0.00034	0.00040	0.3950
HBsAg	0.22708	0.20379	0.2651
Tumor number	0.35162	0.20395	0.0847
Tumor capsule	0.36457	0.18973	0.0547
Cirrhosis	−0.11200	0.17002	0.5101
History of alcohol abuse	−0.00628	0.15346	0.9673
Antiviral therapy	0.16645	0.15298	0.2766
Anti-HCV (IgG)	0.00542	0.45442	0.9905
Statistically significant factors
Tumor size (cm)	0.12838	0.01863	<0.0001
AFP (ng/mL)	0.47461	0.16918	0.0050
PT (seconds)	0.21997	0.05446	<0.0001
Fibrinogen (g/L)	0.19352	0.05718	0.0007
ALB (g/L)	−0.09261	0.01288	<0.0001
TBil (μmol/L)	0.03364	0.00705	<0.0001
TC (mmol/l)	−0.16744	0.08316	0.0441
ALT (U/L)	0.00272	0.00101	0.0071
γ-GT (U/L)	0.00057	0.00027	0.0331
S-index	−1.30610	0.17261	<0.0001
Ascites	0.60656	0.19005	0.0014
PVTT	1.47896	0.33065	<0.0001
Invasion of adjacent	1.01903	1.01903	0.0007
tissues
MVI	0.85045	0.19929	<0.0001
Intrahepatic transfer	0.72404	0.16483	<0.0001
Tumor grade	0.61275	0.15390	<0.0001

Abbreviations: γ-GT, γ-glutamyl transpeptidase; AFP, alpha-fetoprotein; ALB, albumin; ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; HCV, hepatitis C virus; MVI, microvascular invasion; PLT, platelet; PT, prothrombin time; PVTT, portal vein tumor thrombus; SE, standard error; TBil, total bilirubin; TC, total cholesterol.

Table 3 Multivariate analysis of the training cohort

Variables	β	SE	P-value	HR	95% CI
ALB (g/L)	–0.04102	0.01606	0.0106	0.960	0.930	0.990
TBil (μmol/L)	0.02461	0.00785	0.0017	1.025	1.009	1.041
Prothrombin time (seconds)	0.01411	0.07404	0.8488	1.014	0.877	1.173
Fibrinogen (g/L)	0.08868	0.06617	0.1802	1.093	0.960	1.244
Tumor grade	0.47338	0.15924	0.0030	1.605	1.175	2.193
Tumor size (cm)	0.05021	0.02304	0.0293	1.051	1.005	1.100
Ascites	−0.23984	0.22952	0.2960	0.787	0.502	1.234
PVTT	0.86243	0.41234	0.0365	2.369	1.056	5.315
Invasion of adjacent tissues	0.20557	0.34647	0.5530	1.228	0.623	2.422
MVI	0.56170	0.28907	0.0520	1.754	0.995	3.090
Intrahepatic transfer	0.11095	0.24518	0.6509	1.117	0.691	1.807
S-index	0.98779	0.18963	<0.0001	2.685	1.852	3.894

Abbreviations: ALB, albumin; MVI, microvascular invasion; PVTT, portal vein tumor thrombus; SE, standard error; TBil, total bilirubin.

Prognostic nomogram for survival

Based on the independent risk factors that were statistically significant in multivariate analysis, nomograms were constructed to predict the OS of PHCC (Figure 1). In order to construct this nomogram, each subgroup variable was assigned a corresponding score. According to the specific situation of each hepatocellular carcinoma, the scores of each subgroup variable were summarized to obtain a total score, and we then obtained the corresponding 1-, 3-, and 5-year survival probabilities by drawing a straight line vertically down from the total scale to the survival probability scale. In order to compensate for the inability to draw lines on the nomogram to accurately obtain the survival rates, the nomogram scoring system was placed in Table 4, accurately predicting the survival rates. Calibration curves were used to verify the performance of the model in predicting OS. The C-indexes of the model for 1-, 3-, and 5-year OS rates were 0.768 (95% CI=0.695–0.841), 0.757 (95% CI=0.706–0.809), and 0.737 (95% CI=0.696–0.779), respectively, in the training group and 0.679 (95% CI =0.579–0.779), 0.683 (95% CI=0.609–0.757), and 0.679 (95% CI=0.616–0.741), respectively, in the validation group, still performing well in the latter (Figure 2). To further validate the performance of the model for evaluating the prognosis of PHCC, we plotted the tdROC curve for the nomogram, S-index, BCLC, Okuda, and CLIP. The tdAUC curve was used to evaluate the difference between predicted 1-, 3-, and 5-year OS rates of each of the indicators mentioned above (Figure 3). Furthermore, we stratified the training group into three subgroups based on the total points of nomogram (low risk: <200; intermediate risk: 200–250; and high risk: >250). In Figure 4, the Kaplan–Meier curves show excellent prediction results in the nomogram predicting survival.

Figure 1 Nomogram for predicting overall survival after curative resection of PHCC.

Notes: To obtain the corresponding survival probability, sum up each variable score and draw a vertical line from the total points scale to the survival rate scale.

Abbreviations: ALB, albumin; PHCC, primary hepatocellular carcinoma; PVTT, portal vein tumor thrombus; TBil, total bilirubin.

Figure 2 Calibration curves for predicting overall survival probability by the nomogram in the training and validation sets.

Notes: Calibration curves of the prognostic nomogram for 1-year overall survival (A), 3-year overall survival (B) and 5-year overall survival (C) in the training set, correspondingly, calibration curves for 1-year overall survival (D), 3-year overall survival (E), and 5-year overall survival (F) in the validation set.

Figure 3 The tdROC curve of the prognostic nomogram, S-index, CLIP, BCLC and Okuda score in the training cohort.

Note: The tdROC curve for 1-year survival (A), 3-year survival (B), and 5-year survival (C).

Abbreviations: BCLC, Barcelona Clinic Liver Cancer; CLIP, Cancer of the Liver Italian Program; tdROC, time-dependent receiver operating characteristic.

Figure 4 Survival curves stratified by the score calculated by the nomogram in the training cohort (low risk: <200; intermediate risk: 200–250; and high risk: >250).

Table 4 Nomogram scoring system

S-index	Points	PVTT	Points	Tumor grade	Points	TBil	Points	ALB	Points	Tumor size	Points
<28	39	No	28	Grade 1/2	20	0	62	55	73	0	100
≥28	0	Yes	0	Grade 3/4	0	35	30	35	35	15	20
						70	0	15	0	30	0
Total points		1-year survival probability		Total points		3-year survival probability		Total points		5-year survival probability
105		0.1		148		0.1		165		0.1
128		0.2		170		0.2		190		0.2
146		0.3		188		0.3		205		0.3
158		0.4		200		0.4		220		0.4
172		0.5		215		0.5		232		0.5
188		0.6		228		0.6		245		0.6
200		0.7		240		0.7		262		0.7
218		0.8		260		0.8		275		0.8
240		0.9		280		0.9		300		0.9

Abbreviations: ALB, albumin; PVTT, portal vein tumor thrombus; TBil, total bilirubin.

Discussion

A more accurate survival prediction for PHCC after radical surgery is undoubtedly an excellent goal for patients and surgeons. Nevertheless, in terms of predicting survival rates, the conventional staging systems are not satisfactory, including the TNM system and the BCLC staging system.^11,12 If a tool that predicts the prognosis of PHCC can be developed, it will provide better clinical treatment advice regarding prognosis. Based on multivariate analysis, we obtained six independent risk factors for evaluating the prognosis of PHCC. Therefore, we developed a nomogram including S-index, PVTT, tumor grade, tumor size, ALB, and TB. The C-indexes of the nomogram in the training and validation groups performed well in predicting the survival of PHCC. Compared with other traditional staging scores, our model was highly accurate for predicting the prognosis of PHCC. Unfortunately, for predicting the 5-year survival probability of PHCC, the tdAUC of the nomogram was similar to but smaller than the S-index. It is well known that the occurrence and development of tumors are related to many factors, while the S-index combines only laboratory indicators of blood. Therefore, we firmly believe that the nomogram based on multiple factors has greater potential. This also requires much clinical research to validate our model.

According to risk stratification of the nomogram, we can easily divide individuals into different risk subgroups, which may make sense for clinicians to make further treatment plans. For example, patients at high-risk stratification may require further treatment or shorter follow-up periods to improve survival.

As novel tools, nomograms have been regarded by scholars as useful and efficient methods for predicting tumor prognosis in recent years. Neoplastic diseases that have been tracked with nomograms to predict prognosis include gastrointestinal stromal tumors,⁹ colorectal cancer,⁸ pancreatic cancer,¹³ and hepatocellular carcinoma.^3,7

Our nomogram contained six variables in which tumor staging, PVTT, and tumor size were consistent with those of previous studies.^3,7,14 As a carrier for evaluating protein reserves, serum ALB has been associated with many tumors.^15–17 Preoperative low serum ALB was a powerful predictor of prolonged hospital stay and mortality.¹⁸ Many studies have shown that serum ALB represents preoperative nutritional status in cancer patients and can also predict postoperative outcomes.¹⁹ TB is another laboratory indicator included in our model. In recent years, many studies have shown that bilirubin is associated with lymph node metastasis of rectal cancer,²⁰ liver fibrosis of hepatitis B virus,²¹ and invasiveness of hepatocellular carcinoma.²² In other words, cancer patients with elevated bilirubin values have poor prognoses.

To the best of our knowledge, the S-index contained in the model we established is a new indicator for predicting the prognosis of PHCC. The S-index obtained by laboratory indicators accurately evaluates the degree of liver cirrhosis and fibrosis for patients with hepatitis B virus infections.⁶ However, the relationship between the S-index and predicted survival rates after curative resection of PHCC remains unknown. It is well known that cirrhosis and fibrosis are closely related to the occurrence of PHCC. Therefore, we believe the S-index and liver cancer may be closely related and our study confirmed this. The S-index is calculated as 1,000×γ-GT/(PLT × ALB),² containing three important laboratory indicators: γ-GT, PLTs, and ALB. The relationship between ALB and cancer prognosis was described above. Now, we discuss the relationship between the other two variables and PHCC. High γ-GT is closely related to the occurrence of many diseases, including diabetes, liver disease, cardiovascular disease, and all-cause mortality.²³ This may be related to the redox regulation and drug resistance of γ-GT in tumor development.²⁴ Carr et al found that PLTs were closely related to the migration and invasion of liver cancer.²⁵ The survival rates of hepatitis B virus-associated hepatocellular carcinoma with antiplatelet therapy were greatly improved.²⁶

For specific individuals, our nomogram is a very intuitive and accurate predictive tool for survival after curative resection of PHCC. The laboratory indicators included in our nomogram are easily obtained and calculated, including TB, ALB, and the S-index. Tumor size, tumor grade, and PVTT can be obtained by abdominal ultrasound, computed tomography, MRI, and pathological findings. Moreover, the indicators included in our nomogram are objective, avoiding the adverse effects of subjectivity of the surgeon. Nevertheless, our study has certain limitations. First, this was a single-center study, that may limit the scope of its application; second, our sample sizes were relatively small, especially in the validation group; third, in the training and validation groups, the C-index of our model was not highly consistent; finally, our model is only suitable for predicting the survival rates after curative resection of PHCC and is not applicable for patients with incurable advanced liver cancer or patients receiving other treatments.

In summary, we developed an accurate and intuitive tool for predicting the 1-, 3-, and 5-year survival rates after curative resection of PHCC. This nomogram can provide corresponding treatment assistance for patients with PHCC and may provide a basis for clinicians to make treatment decisions. Nevertheless, whether our model is applicable to other inpatients, patients in other regions or even other countries has yet to be verified by a large number of clinical studies by scholars.

Author contributions

Lei Chen and Bin-Bin Cai are co-first authors. All authors contributed to data analysis, drafting and revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Supplementary materials

Figure S1 The tdROC curve for S-index.

Table S1 The optimal cutoff value of S-index

AUC	AUC 95% CI	Cutoff	Sensitivity (%)	Specificity (%)	Youden’s index
0.74	0.68–0.79	0.28	73.26	67.84	0.41

Abbreviation: AUC, area under the curve.

Disclosure

The authors report no conflicts of interest in this work.