2.1 Patients
This study adhered to the Declaration of Helsinki, was approved by Lanzhou University Second Hospital's Ethics Committee, exempted from subjects' informed consent, and approval number: 2024A-1267.
We retrospectively collected data from consecutive patients who underwent curative LR for HCC at our hospital between January 2018 and January 2023. The inclusion criteria were as follows: (1) patients aged ≥18 years who underwent initial LR; (2) histopathologically confirmed solitary HCC with R0 resection; (3) preoperative liver spectral CT performed within 2 weeks before surgery; and (4) complete clinical and pathological data. Exclusion criteria included: (1) prior treatment for HCC, such as microwave ablation, radiofrequency ablation, or transarterial chemoembolization (TACE); (2) confirmed distant metastasis before surgery; (3) coexistence of other malignancies; (4) poor CT image quality; and (5) early postoperative death due to severe complications or loss to follow-up. A total of 85 patients were ultimately included in the study. The detailed patient enrollment flowchart is shown in Fig. 1.
2.2 Collection of data
We retrospectively collected baseline admission data, including demographic characteristics, laboratory results, perioperative variables, and histopathological findings. The surgical resection margin was defined as the shortest pathological distance from the tumor edge to the liver transection line. Based on pathological reports, patients were categorized into three groups according to resection margin width: WM (≥1 cm), NM (≥0.5 cm to <1 cm), and ENM (<0.5 cm) [23, 24]. Detailed admission data are provided in Supplementary Material 1.
2.3 Image acquisition and analysis
All patients were scanned using Discovery CT 750 HD (GE Healthcare, Waukesha) and Revolution CT (GE Medical Healthcare). The detailed spectral CT scanning parameters are provided in Supplementary Material 1.
Two abdominal radiologists with 7 and 10 years of diagnostic experience, respectively, independently and blindly evaluated the patients’ CT images. The assessment focused on the presence of liver cirrhosis, splenomegaly, gastroesophageal varices (GEVs), spontaneous portosystemic shunt (SPSS), ascites, peritumoral arterial phase enhancement, tumor margin, tumor capsule, intratumoral necrosis, and radiologic vascular invasion (RVI). In cases of disagreement, a consensus was reached through discussion. The CT diagnostic criteria for clinically significant portal hypertension (CSPH) included the presence of splenomegaly along with at least one of the following: GEVs, SPSS, or ascites [25]. Detailed definitions of CT features are provided in Supplementary Material 1, Table S1.
2.4 Histogram analysis
ID images from the portal venous phase of spectral CT were stored in DICOM format and imported into FireVoxel software (FireVoxel, version 462; https://www.firevoxel.org). Two radiologists, each with over five years of experience in hepatic imaging, independently performed whole-tumor histogram analysis under blinded conditions. Any discrepancies were resolved through consensus. The entire HCC lesion was manually segmented slice by slice along its boundary, with each ROI encompassing as much of the tumor as possible, including necrotic, cystic, and hemorrhagic areas. To minimize partial volume effects, the ROI was drawn slightly smaller than the actual lesion boundary. After ROI delineation, the software automatically generated histogram parameters based on a three-dimensional volume of interest (VOI), including minimum (Min), maximum (Max), mean (Mean), standard deviation (SD), variance, skewness, kurtosis, entropy, and percentiles (1st–99th) (Fig. 2). To ensure the consistency and reliability of the extracted histogram parameters, the intraclass correlation coefficient (ICC) was used to evaluate interobserver agreement.
2.5 Follow-up and Endpoints
All patients were followed up at 1, 3, and 6 months postoperatively, and then every 6 months thereafter. Follow-up assessments included measurements of serum alpha-fetoprotein (AFP) levels, liver function tests, and imaging studies (abdominal ultrasound, contrast-enhanced CT, or MRI). Tumor recurrence was defined as intrahepatic recurrence or extrahepatic metastasis, primarily diagnosed based on imaging findings or confirmed by histopathology through liver biopsy.
The primary endpoint of this study was ER, defined as the occurrence of intrahepatic or extrahepatic tumor recurrence within 2 years after curative resection. The time of confirmed recurrence and the characteristics of the recurrent lesions were recorded. Recurrence-free survival (RFS) was defined as the interval between the date of curative surgery and the date of tumor recurrence or the last follow-up. The follow-up deadline was set as February 1, 2025.
2.6 Statistical analysis
Data processing and analysis were conducted using IBM SPSS Statistics (Version 26.0; IBM, New York, USA), R (Version 4.3.2; https://www.r-project.org/) and Zstats v1.0 (www.zstats.net).
The inter-observer agreement between the two radiologists was assessed using Cohen's Kappa coefficient for categorical variables. The normality of continuous variables was tested using the Shapiro-Wilk test. Continuous variables with a normal distribution are presented as Mean ± SD, while non-normally distributed variables are expressed as Median (Q1, Q3). The differences between continuous variables across groups were compared using the independent samples t-test for normally distributed data and the Mann-Whitney U test for non-normally distributed data. Categorical variables are presented as frequencies (%) and analyzed using the χ² test or Fisher's exact test.
Multicollinearity was evaluated by calculating the variance inflation factor (VIF). Univariate and multivariate Cox regression analyses were performed, and variables with P<0.05 in the multivariate analysis were used to construct the ER prediction model. Continuous variables were dichotomized based on the median, and Kaplan-Meier survival curves were plotted. Differences in survival curves were analyzed using the Log-rank test. The predictive performance of the model was assessed by the area under the time-dependent receiver operating characteristic (td-ROC) curve (AUC). Model calibration was evaluated using calibration curves, and the overall net benefit of the model was assessed using decision curve analysis (DCA). A P-value of <0.05 was considered statistically significant.