Nomogram for predicting overall survival in patients with large (>5 cm) hepatocellular carcinoma based on real-world practice
Article information
Abstract
Background/Aim
Patients with large (>5 cm) hepatocellular carcinoma (HCC) have limited treatment options, thus necessitating the identification of prognostic factors and the development of predictive tools. This study aimed to identify prognostic factors and to construct a nomogram to predict survival outcomes in patients with large HCC.
Methods
A cohort of 438 patients, who were diagnosed with large HCC at a tertiary hospital between 2015 and 2018, was analyzed. Cox proportional hazards models were used to identify key prognosticators of overall survival (OS), and an independent set of prognostic factors was used to develop a nomogram. The discrimination and calibration abilities of the nomogram were assessed and internal validation was performed using cross-validation and bootstrapping methods.
Results
During a median follow-up of 9.3 months, the median OS was 9.9 months, and the 1-year OS rate was 43.9%. Multivariable Cox regression analysis revealed that performance status, modified albumin-bilirubin grade, tumor size, extent of portal vein tumor thrombosis, and initial treatment significantly affected OS. The newly developed nomogram incorporating these variables demonstrated favorable accuracy (Harrell’s concordance index, 0.807).
Conclusions
The newly developed nomogram facilitated the estimation of individual survival outcomes in patients with large HCC, providing an acceptable level of accuracy.
INTRODUCTION
Despite the identification of high-risk groups and risk factors for hepatocellular carcinoma (HCC), along with active recommendations for early and regular screening, including tumor markers and ultrasonography, patients are often diagnosed at an advanced stage with large (>5 cm) tumors with or without symptoms.1 Large tumor size can limit the implementation of curative treatment, such as radiofrequency ablation, surgical resection, or even liver transplantation, even if it is diagnosed as Barcelona Clinic Liver Cancer (BCLC) stage A without any other prognostic factors, such as portal vein tumor thrombosis (PVTT).2,3 Despite the potential for long-term survival, surgical resection in these patients is challenging due to the complexity of the required procedure(s) and their association with an elevated risk for morbidity and mortality.4-6 According to data from the Korean Nationwide Cancer Registry, 19.8% of all patients and 34.5% of patients with BCLC stage A underwent surgical resection between 2012 and 2014.7
Although various local treatments, such as transarterial chemoembolization (TACE) or radiotherapy (RT), are considered when surgery (resection or liver transplantation) is not feasible, there are limited data and little consensus regarding treatment strategies for large HCC. Patients with large HCC may present with varying tumor extent, vascular invasion, PVTT, and liver function abnormalities.8,9 As such, identifying prognostic factors is essential for predicting survival outcomes.
Recently, a nomogram, which is a graphical representation of a statistical predictive model, has been proposed as a tool for the individualized prediction and stratification of patients.10 Accordingly, the objective of this study was to identify prognostic factors in patients with large HCC (>5 cm) based on a data from a tertiary cancer center registry and to develop a nomogram incorporating these predictors.
Methods
1. Patients
From the HCC registry, 438 patients with newly diagnosed, previously untreated, large (>5 cm) HCC treated at Samsung Medical Center between 2015 and 2018, were identified.11 This single-institutional retrospective cohort study was approved by the institutional review board of Samsung Medical Center (No. 2023-06-092), and requirements for informed consent were waived because only anonymized, routinely collected data gathered during hospital visits, were used. This study is reported in accordance with the strengthening the reporting of observational studies in epidemiology guidelines (Supplementary Table 1).
2. Data collection
Data extracted included (1) patients, (2) tumors, and (3) treatment-related factors. Patient factors included age (years), sex, Eastern Cooperative Oncology Group performance status (ECOG PS), Child-Pugh class, and albuminbilirubin index (ALBI) score. Tumor-related factors included etiology, levels of tumor markers (alpha-fetoprotein [AFP] and protein induced by vitamin K absence-II [PIVKA-II]), tumor size, vessel/bile duct invasion, extent of PVTT, and tumor stage (as assessed according to the modified Union for International Cancer Control [mUICC] and BCLC).3,12,13 Treatment-related factors included details of the initial treatment.
Baseline liver function was assessed using both the Child-Pugh classification and the ALBI score. The ALBI score was calculated using the following equation.14
ALBI=(-0.085×albumin [g/L]) + (0.66×log[bilirubin] [µ mol/L]).
Patients were categorized into four groups based on modified ALBI (mALBI) grade, as follows: grade 1, ≤-2.6; grade 2a, -2.60 to -2.27; grade 2b, -2.27 to -1.39; and grade 3, >-1.39.15,16 The extent of PVTT was classified into five grades based on the Japanese classification, as follows: (1) Vp0, no tumor thrombus in the portal vein; (2) Vp1, presence of a tumor thrombus distal to (but not in) the second order branches of the portal vein; (3) Vp2, presence of a tumor thrombus in the second-order branches of the portal vein; (4) Vp3, presence of a tumor thrombus in the first-order branches of the portal vein; and (5) Vp4, presence of a tumor thrombus in the main trunk of the portal vein and/or a portal vein branch contralateral to the primarily involved lobe.12 Information regarding the initial treatment performed in the patients is summarized in Supplementary Table 2. Eighty-eight (20.1%) patients underwent curative surgical resection (hepatectomy), while more than half (n=231 [52.7%]) underwent non-surgical local treatments, as follows: TACE (n=104 [23.7%]); TACE followed by scheduled RT (n=78 [17.8%]); RT with sorafenib (n=29 [6.6%]); and RT alone (n=20 [4.6%]).
3. Statistical analysis
The primary endpoint was overall survival (OS), calculated from the date of diagnosis to the date of death or the last follow-up visit, and analyzed using the Kaplan–Meier method. Progression-free survival (PFS) was calculated from the date of diagnosis to the date of disease progression, last follow-up visit, or death. To identify the prognostic factors related to OS and PFS, multivariable analysis was performed using the Cox proportional hazard model for statistically significant variables in the univariate analysis. The variance inflation factor was also calculated to assess multicollinearity among the factors included in the multivariable analysis. All factors exhibited a variance inflation factor <10.
The results of the coefficients introduced into the multivariable Cox proportional hazards model were used to build a new nomogram for predicting OS at 1 and 2 years. Regarding liver function, a quantitative mALBI system was selected instead of the Child-Pugh classification to eliminate multicollinearity. The final nomogram was internally validated using 1,000 bootstrap simulations. Then predictive performance of the nomogram was determined by calculating the discriminatory potential using Harrell’s concordance index (C-index) and Heagerty’s integrated area under the curve (iAUC) and by plotting calibration curves for survival probabilities at 1 and 2 years. Finally, random survival forest analysis was performed to identify the importance of the variables. In all analyses, a two-sided P<0.05 was considered statistically significant. All analyses were performed using the R package version 4.2.3 (R Foundation for Statistical Computing, Vienna, Austria).
Results
1. Patients and tumor characteristics
The baseline characteristics of the patients and tumors are summarized in Table 1. The median patient age was 57 years (interquartile range [IQR], 51-67). Underlying liver diseases included the following: hepatitis B virus infection, 293 (66.9%); non-viral infection, 121 (27.6%); hepatitis C virus infection, 17 (3.9%); and hepatitis B and C viral infections, seven (1.6%). Regarding baseline liver function, 357 (81.5%) and 244 patients (55.7 %) were classified as Child-Pugh class A and mALBI grade 1, respectively. More than one-half of patients (n=233, 53.2%) exhibited multiple lesions. The median tumor size was 15.0 cm (IQR, 9.8-15.0), and 297 patients (67.8%) exhibited tumors ≥10 cm in size. Based on the mUICC staging system, 134 patients (30.6%) had stage IVA disease, 119 (27.2%) had stage III disease, and 101 (23.1%) had stage IVB disease; BCLC stage C was the most common (n=303, 69.2%).
2. Survival outcomes and prognostic factors
The median follow-up was 9.3 months (IQR, 3.5-31.4) for the entire cohort and 71.3 months (IQR, 61.2-83.1) for survivors. A total of 357 patients (81.5%) died during the study period. Median PFS and OS were 3.7 and 9.9 months, respectively. The 1-year PFS and OS rates were 18.9% and 43.9%, respectively (Fig. 1). The 2-year PFS and OS rates were 12.2% and 30.5%, respectively.
Regarding patient-related factors, ECOG PS (hazard ratio [HR], 4.10; 95% confidence interval [CI], 2.09-8.03]) and mALBI grade (grade 1 vs. grade 2a; HR, 1.49; 95% CI, 1.12-1.98 vs. grade 2b/3; HR, 2.41; 95% CI, 1.70-3.45) were associated with OS outcomes. In addition, tumor size (5-10 vs. >10 cm; HR, 1.04; 95% CI, 1.01-1.08), PVTT classification (Vp0 vs. Vp3-4; HR, 1.43; 95% CI, 1.02-2.00), and BCLC stage (A-B vs. C-D; HR, 1.65; 95% CI, 1.03-2.82) were associated with OS outcomes. Finally, multivariable analysis also demonstrated that the type of initial treatment significantly affected OS (all P<0.05) (Table 2).
Multivariable analyses revealed that patient-related factors, including ECOG PS (HR, 3.70; 95% CI, 1.99-6.84) and mALBI grade 3-4 (HR, 1.68; 95% CI, 1.18-2.41) and tumorrelated factors of AFP (HR, 1.00; 95% CI, 1.00-1.00), PIVKA-II (HR, 1.00; 95% CI, 1.00-1.00), Vp3-4 classification (HR, 1.42; 95% CI, 1.01-1.99), and lymph node metastasis (HR, 1.41; 95% CI, 1.05-1.90) significantly contributed to inferior PFS outcomes (Supplementary Table 3).
3. Nomogram
Considering the multiple factors related to OS outcomes, a nomogram for predicting the 1- and 2-year OS outcomes was developed (Fig. 2). In the random survival forest analysis, initial treatment was the most important prognostic factor, followed by BCLC stage and mALBI grade (Table 3). Tumor size and etiology were less influential than other factors. The C-index and the iAUC of the nomogram were 0.807 (95% CI, 0.749-0.855) and 0.880 (95% CI, 0.865-0.911), respectively (Table 4). The internal calibration plot for predicting the 1- and 2-year survival probability to validate the nomogram confirmed that the predicted survival rates correlated well with the actual survival rates at 1 and 2 years (Fig. 3). Harrell’s C-indices for 1- and 2-year OS were 0.808 and 0.798, respectively. The iAUC for 1- and 2-year OS were 0.873 and 0.885, respectively.
![Figure 2.](/upload//thumbnails/jlc-2023-08-10f2.jpg)
Nomogram for predicting 1-, and 2-year survival rates in patients with large hepatocellular carcinoma. The nomogram summed the points identified on the scale for each variable. The total points projected on the bottom scale indicate the probabilities of survival. ECOG PS, Eastern Cooperative Oncology Group performance status; mALBI, modified albumin-bilirubin index; AFP, alpha fetoprotein; BSC, best supportive care; Tx, treatment; BCLC, Barcelona Clinic Liver Cancer.
4. Survival outcomes according to the mUICC and BCLC systems
According to the mUICC system, patients with stages IVA and IVB disease exhibited similar OS outcomes (Fig. 4A), whereas significant differences in OS were observed based on the BCLC system (Fig. 4B). However, 303 patients had BCLC stage C disease. To identify differences in OS based on the new nomogram, the patients were divided into quartiles according to total points. A significant difference in OS among these classes was identified compared to the classes categorized based on mUICC (Fig. 4C). Additionally, a significant difference was observed in the OS of the 303 patients with BCLC stage C when categorizing patients based on the new nomogram (Fig. 4D). A significant difference in PFS was observed when patients were stratified according to the new nomogram (Supplementary Fig. 1).
Discussion
In this study, we developed a nomogram for predicting OS in patients with large HCC (>5 cm), which yielded accurate predictions based on internal validation. Although several studies have identified prognostic factors in this patient population, estimating individual survival outcomes remains difficult. This nomogram can help physicians predict survival probability based on patient- and tumor-specific covariates for individualized treatment decision making.
In the current analysis, we identified that an initial treatment modality favoring surgical resection was an independent prognostic factor and the most relevant factor affecting OS. Although there are no absolute size contraindications for surgical resection, large tumors are generally considered to be unresectable. Recently, several reports have described the benefits of surgical resection in this population with a tumor size exceeding the Milan criteria.17,18 Even in tumors ≥10 cm, liver resection yielded better OS outcomes than TACE, even after adjusting for clinical factors.19-22 Factors involved that may have influenced these retrospective results by overestimating the benefits of surgical resection.
In addition, non-surgical local treatment was associated with better OS outcomes than sorafenib alone in the developed nomogram. Peng et al.23 reported that the combination of TACE and radiofrequency ablation with sorafenib yielded better OS outcomes (median OS, 12 vs. 8 months; P=0.001) than sorafenib alone for recurrent HCC >5 cm. Although a phase III trial failed to improve OS outcomes after sorafenib with TACE compared to sorafenib alone in advanced HCC, combination therapy significantly improved PFS, tumor response rate, and time to progression.24,25 Regarding RT, a recent phase II trial demonstrated favorable median OS outcomes of 24.6 months after concurrent chemoradiotherapy with 5-fluorouracil and leucovorin followed by sorafenib. The median tumor size among the 47 patients was 8.4 cm (IQR, 6.5-12.0).26 Additionally, a recent randomized trial reported that the addition of RT to TACE was associated with better PFS (3 months, 86.7% vs. 34.3%), longer time to progression (31 vs. 12 weeks), and OS (median, 55 vs. 43 weeks) than sorafenib alone in patients with macroscopic vascular invasion.27 Several meta-analyses have shown that concurrent RT with sorafenib has improved OS compared to sorafenib monotherapy.28,29 Although we could not analyze the difference among non-surgical treatments due to the limited number of patients, Li et al.28 reported that RT may be the best choice for combination therapy compared with TACE or hepatic artery infusion chemotherapy with sorafenib. In summary, we suggest that surgical resection is the preferred initial treatment whenever possible. However, for patients who are ineligible for surgical resection, other local treatments (e.g., TACE and RT) should be considered in combination with systemic therapy rather than with systemic therapy alone.
Liver functional reserve is important for determining the initial treatment option and treatment outcomes in patients with HCC. Serum albumin and bilirubin levels are reliable markers of decline in liver function. The ALBI score is one of the most useful markers for estimating liver function.14,15,30 Fang et al.31 found that large HCCs with ALBI grade 1 exhibited comparable outcomes to those with BCLC stage A disease, and large HCCs with ALBI grade 2 or 3 exhibited similar outcomes to those with BCLC stage B disease. Therefore, the ALBI grade appears to help stratify this population. Furthermore, ALBI grade 2-3 was related to inferior OS outcomes in 143 patients with HCC ≥10 cm.22 Considering the relative importance of mALBI score, liver function assessed by mALBI has an acceptable discriminative potential to predict the prognosis of patients, even among those with large HCC.
In random survival forest analysis, PVTT classification was another important factor affecting OS outcomes. Although various treatment options have been considered for PVTT, the presence of PVTT was associated with dismal outcomes.32-34 The presence of PVTT could promote intrahepatic tumor progression, treatment failure, and deterioration of liver function.35-37 Mähringer-Kunz et al.38 reported that the extent of PVTT and OS were significantly associated. The median OS for Vp4 was 4.8 months compared to 14.6 months for Vp1. Moreover, Li et al.39 demonstrated that main trunk involvement in PVTT was related to OS outcomes in patients with main PVTT treated with TACE. Given the suboptimal outcomes for PVTT, it is essential to explore combination treatment strategies to improve survival outcomes.40
The present study had several limitations, the first of which were the inherent limitations of its single-center, retrospective design. For example, selection bias in the treatment modality could significantly affect OS outcomes (Supplementary Table 4). Because treatment modality is mainly decided based on tumor stage in clinical practice, it should not be selected based on this nomogram until further validation is performed. However, using the total points from the current nomogram (Supplementary Fig. 2) enabled discrimination between patients in each treatment modality. However, because most patients in the surgery group were allocated to the first quartile of the developed nomogram, additional validation with an external dataset is needed to predict outcomes in the surgery group. Second, the current nomogram should be validated externally in other centers, where the selection criteria for initial treatment modalities may vary substantially. However, the strength of our nomogram is the relatively large sample size and real-world data. In addition, we observed acceptable results in internal validation with bootstrapping. Finally, the nomogram did not include subsequent treatments after recurrence, which could have affected OS outcomes. However, our nomogram results, based on initial diagnosis and treatment, could be informative and helpful to both patients and physicians. Incorporating patient and treatment factors with tumor factors, such as BCLC stage, PVTT classification, AFP value, and tumor size, demonstrated an improved prediction of OS compared with the traditional mUICC and BCLC systems. Further studies, however, are required to validate the use of this nomogram in clinical practice.
In conclusion, we identified several prognostic factors related to survival outcomes in patients with large HCC. In addition, the nomogram, based on patient and tumor characteristics, demonstrated acceptable accuracy in this patient population. Continued efforts to refine and validate nomogram prediction tools are necessary to assist physicians in appraising the anticipated course of the disease and aid in decision-making for individualized treatment.
Notes
Conflict of Interest
The authors have no conflicts of interest to disclose.
Ethics Statement
This single-institutional retrospective cohort study was approved by the institutional review board of Samsung Medical Center (No. 2023-06-092), and requirements for informed consent were waived because only anonymized, routinely collected data gathered during hospital visits, were used.
Funding Statement
This study was supported by the Korean Liver Cancer Association Research Award (2021).
Data Availability
The datasets generated and analyzed during the current study are not publicly available due to institutional data protection law and confidentiality of patient data but are available from the corresponding author on reasonable request in person.
Author Contribution
Conceptualization: JIY, HCP
Data curation: JIY, HCP, UYH, HYL, MJG, YHP
Formal analysis: NK
Investigation: JIY, HCP, UYH, HYL, MJG, YHP
Methodology: NK, JIY
Project administration: JIY, HCP
Resources: JIY, HCP
Software: NK
Supervision: JIY, HCP
Validation: NK, JIY, HCP
Visualization: NK
Writing-original draft: NK, JIY
Writing-review & editing: all authors
Approval of final manuscript: all authors
Supplementary Material
Supplementary data can be found with this article online https://doi.org/10.17998/jlc.2023.08.10.
Acknowledgements
None.