- 1Interventional Therapy Center for Oncology, Beijing You’an Hospital, Capital Medical University, Beijing, China
- 2National Center for Infectious Diseases, Beijing Ditan Hospital, Capital Medical University, Beijing, China
- 3Changping Laboratory, Beijing, China
Introduction: This study aimed to develop a prognostic nomogram for predicting the recurrence-free survival (RFS) of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) patients with low preoperative platelet-albumin-bilirubin (PALBI) scores after transarterial chemoembolization (TACE) combined with local ablation treatment.
Methods: We gathered clinical data from 632 HBV-related HCC patients who received the combination treatment at Beijing You’an Hospital, affiliated with Capital Medical University, from January 2014 to January 2020. The patients were divided into two groups based on their PALBI scores: low PALBI group (n=247) and high PALBI group (n=385). The low PALBI group was then divided into two cohorts: training cohort (n=172) and validation cohort (n=75). We utilized eXtreme Gradient Boosting (XGBoost), random survival forest (RSF), and multivariate Cox analysis to pinpoint the risk factors for RFS. Then, we developed a nomogram based on the screened factors and assessed its risk stratification capabilities and predictive performance.
Results: The study finally identified age, aspartate aminotransferase (AST), and prothrombin time activity (PTA) as key predictors. The three variables were included to develop the nomogram for predicting the 1-, 3-, and 5-year RFS of HCC patients. We confirmed the nomogram’s ability to effectively discern high and low risk patients, as evidenced by Kaplan-Meier curves. We further corroborated the excellent discrimination, consistency, and clinical utility of the nomogram through assessments using the C-index, area under the curve (AUC), calibration curve, and decision curve analysis (DCA).
Conclusion: Our study successfully constructed a robust nomogram, effectively predicting 1-, 3-, and 5-year RFS for HBV-related HCC patients with low preoperative PALBI scores after TACE combined with local ablation therapy.
Introduction
Liver cancer, with hepatocellular carcinoma (HCC) as its main subtype, ranks as the third leading cause of cancer-related deaths globally (1). In developing regions such as China and Africa, chronic hepatitis B virus (HBV) infection serves as the predominant factor contributing to HCC (2, 3). Despite the considerable progress in treatment modalities witnessed in recent years, the prognosis of HCC patients remains unsatisfactory, primarily due to its high recurrence rates (4, 5). The prognosis of HCC is often closely related to liver function, so it is crucial to conduct a thorough liver function assessment before treatment (6–8). The traditional method for clinically assessing liver function is the Child-Pugh class, which includes five indicators: albumin, bilirubin, prothrombin time, ascites, and hepatic encephalopathy (9). However, the grading values are solely defined by clinical experience, and ascites and hepatic encephalopathy are susceptible to subjective factors, making it challenging to objectively and accurately evaluate liver function. In 2015, the albumin-bilirubin (ALBI) score, solely based on laboratory parameters, was introduced, and numerous clinical studies have validated its accuracy in predicting prognosis (10–12). However, as HCC often occurs in patients with cirrhosis, the degree of portal hypertension is a crucial prognostic factor. The ALBI score fails to assess the extent of portal hypertension. The platelet-albumin-bilirubin (PALBI) score, incorporating platelet count as a marker for portal hypertension, has been confirmed as a more ideal method for assessing liver function (13–15).
PALBI score categorizes HCC patients into three grades, with higher grades indicating poorer liver function and worse prognosis. Our team has also confirmed that there was a significant statistical difference in prognosis between HBV-related HCC patients with low and high PALBI grades. Specifically, patients with high PALBI scores had a more unfavorable prognosis compared to patients with low PALBI scores. However, although there were differences between the two, the prognosis of low PALBI patients was not very optimistic and their recurrence-free survival (RFS) was not much better than that of high PALBI patients. This finding suggests that for patients with low PALBI scores, caution should still be exercised in managing their prognosis, to fully understand and address potential recurrence risks. Unfortunately, there is currently no specialized prognostic prediction model for low PALBI patients.
Therefore, in this study, we not only investigated the impact of high and low PALBI scores on HBV-related HCC prognosis, but also committed to developing an RFS prediction nomogram particularly for low PALBI patients. In the process of establishing the nomogram, we integrated clinical characteristics, biochemical indicators, and tumor biology features of patients with low PALBI scores, employing advanced machine learning algorithms and statistical methods to systematically analyze and select factors influencing tumor recurrence. The goal is to construct an accurate, reliable, and clinically practical predictive model to assist clinicians in better assessing the recurrence risk of patients, providing essential reference for further optimization of treatment plans.
Materials and methods
Patients enrolled
In this retrospective study, data from 632 patients diagnosed with HBV-related HCC and treated with a combination of transarterial chemoembolization (TACE) and ablation treatment at Beijing You’an Hospital, affiliated with Capital Medical University, were reviewed, from January 2014 to January 2020. All participants enrolled presented with a confirmed diagnosis of HCC according to the criteria stipulated by the American Association for the Study of Liver Disease (AASLD). The combination treatment was undertaken by physicians possessing a minimum of 5 years of expertise in the field. TACE was initially employed to target the cancerous lesions, delivering chemotherapy drugs and embolic agents directly to the tumor-feeding blood vessels. Following the TACE procedure, local ablation therapy was performed within a time frame of 1 to 2 weeks.
Here are the inclusion and exclusion criteria for the study. Inclusion criteria: (1) Patients diagnosed with hepatitis B virus (HBV)-related primary HCC; (2) Received the combination of transarterial chemoembolization (TACE) and ablation treatment; (3) Complete clinical and follow-up data. Exclusion criteria: (1) Non-primary HCC; (2) Not HBV-related HCC; (3) Received other anti-tumor therapies before the combination treatment; (4) Suffering from other malignancies or systemic diseases; (5) Incomplete clinical or follow-up data. Figure 1 illustrates the flowchart outlining the process of selecting qualified patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g001.jpg)
Figure 1 Patient enrollment flowchart. HBV, hepatitis B virus; HCC, hepatocellular carcinoma; TACE, transcatheter arterial chemoembolization; PALBI, platelet-albumin-bilirubin; XGBoost, eXtreme Gradient Boosting; C-index, Harrell’s concordance index; AUC, area under the curve.
The study has obtained approval from the ethics committee of Beijing You’an Hospital, affiliated with Capital Medical University. The committee conducted a thorough review of the research protocol to ensure that the study adheres to ethical guidelines and regulatory requirements. Given the retrospective nature of the study, where data was collected from past medical records, and the stringent measures implemented to protect patient’s privacy, the need for obtaining informed consent from the subjects was deemed unnecessary.
Data collection
The data collected in this study comprises general information, laboratory examination data, and pathological data. General information: age, gender, smoking history, drinking history, family history, antiviral history, hypertension, diabetes and so on. Laboratory examination data: white blood cell (WBC), red blood cell (RBC), lymphocyte (Lym), neutrophil (Neu), monocyte (Mon), eosinophils, basophils, platelet (PLT), hemoglobin (Hb), alanine aminotransferase (ALT), aspartate aminotransferase (AST), AST/ALT, gamma glutamyl transpeptidase (GGT), alkaline phosphatase (ALP), total protein (TP), albumin, globulin, prealbumin, total bilirubin (TBIL), direct bilirubin (DBIL), bile acid, uric acid, BUN, cholesterol, glucose, potassium, natrium, chlorine, creatinine, fibrinogen (Fib), thrombin time (TT), international normalized ratio (INR), activated partial thromboplastin time (APTT), activated partial thromboplastin time ratio (APTTR), prothrombin time (PT), prothrombin time activity (PTA), prothrombin time ratio (PTR) and so on. Pathological data: tumor size, tumor number, Barcelona Clinic Liver Cancer (BCLC) stage and Child-Pugh class.
Patients follow-up
Within one month after the combination treatment, all patients were required to undergo a follow-up check at the outpatient department of our hospital. They were followed every 3 months starting from the 2nd month to 1 year and every six months thereafter. Tumor recurrence refers to the discovery of new lesions in the liver through imaging examinations such as contrast-enhanced computed tomography (CT) and magnetic resonance imaging (MRI) or confirmed by pathological evidence based on tumor biopsy. Recurrence-free survival (RFS) was calculated as the time interval from the date of undergoing treatment to the date of tumor recurrence or the last follow-up date, and was presented in months. The last follow-up date for this study was December 1, 2023 and the median follow-up time was 23.2 months.
Statistical analysis
Categorical variables were expressed as numbers and proportions, and their comparison was executed through the Chi-square test. Continuous variables were represented as mean values accompanied by standard deviations, and their comparison was carried out using Student’s t-test or Mann-Whitney U test. The significance level was P=0.05, where P<0.05 indicated statistical difference. Platelet-albumin-bilirubin (PALBI) score was calculated as: 2.02 * log10 bilirubin - 0.37 * (log10 bilirubin)2 - 0.04 * albumin - 3.48 * log10 platelets + 1.01 * (log10 platelets)2, bilirubin was expressed in μmol/L, albumin in g/L and blood platelet count in 1000/μL. PALBI was divided into three grades, with higher grades indicating poorer liver function. PALBI grade 1: PALBI score≤ -2.53; PALBI grade 2: -2.53 < PALBI score ≤ -2.09; PALBI grade 3: PALBI score > -2.09. In this study, PALBI grade 1 was classified as low PALBI, while grade 2 and 3 were included as high PALBI. Survival analysis was conducted using Kaplan-Meier curves, and the log-rank test was utilized for data comparison. Subsequently, patients categorized under the low PALBI group were randomly split into training and validation cohorts in a 7:3 ratio. Then, we utilized a combination of advanced analytical techniques, including eXtreme Gradient Boosting (XGBoost), random survival forest (RSF), and multivariate Cox analysis, to delve into the factors influencing RFS following treatment for HBV-related HCC with low PALBI score. Based on the screened risk factors, we developed a nomogram to predict the 1-year, 3-year and 5-year RFS of low PALBI HBV-related HCC patients. Furthermore, using the scores generated by the nomogram, patients were effectively categorized into two distinct risk groups: low-risk and high-risk groups. The survival curves of the two groups were also compared by the log-rank test. Moreover, the discrimination ability, calibration, and clinical utility of the nomogram were then comprehensively evaluated to ensure its robustness and reliability.
Results
Comparison of RFS between high PALBI group and low PALBI group
A total of 632 patients diagnosed with HBV-related HCC who underwent TACE combined with ablation therapy were included in this retrospective study. The patients were stratified into two groups according to their preoperative PALBI scores: the low PALBI group comprised 247 individuals, while the high PALBI group consisted of 385 patients. It was confirmed that the high PALBI group exhibited higher recurrence rates compared to low PALBI group. However, delving deeper into the data, we unearthed a rather intriguing revelation: while there existed a disparity between these two groups, the RFS of patients with low PALBI scores did not demonstrate a substantial advantage over that of their high PALBI counterparts, as shown in Figure 2. It implies that despite falling into the low PALBI category, patients were not necessarily immune to recurrence risks and warrant meticulous attention in managing their prognosis. This revelation underscores the imperative of developing specialized prognostic prediction models tailored to this subgroup, enabling clinicians to accurately assess recurrence risks and implement targeted interventions that optimize patient outcomes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g002.jpg)
Figure 2 Kaplan-Meier curves of high PALBI group and low PALBI group. PALBI, platelet-albumin-bilirubin; RFS, recurrence-free survival.
Baseline characteristics of the low PALBI patients
The group of 247 patients with low PALBI scores was divided into two cohorts, the training cohort (n=172) and the validation cohort (n=75), through random assignment, maintaining a ratio of 7 to 3. To ensure the validity of our analysis, a meticulous comparison of baseline characteristics between the training and validation cohorts was conducted. This comparative analysis aimed to ascertain whether there were any differences in the distribution of variables between the two cohorts. The results, as summarized in Table 1, unequivocally demonstrated that there existed no statistically significant difference (P > 0.05) in the baseline characteristics between the training and validation cohorts, indicating their suitability for further studies.
Screening risk factors associated with RFS through XGBoost, RSF, and multivariate Cox analysis
In this study aimed at screening the factors influencing RFS, a multifaceted approach combining machine learning techniques and traditional statistical analysis was employed. At first, we utilized XGBoost, a powerful machine learning algorithm known for its effectiveness in handling complex datasets, to screen for relevant variables affecting RFS (16, 17). XGBoost is particularly adept at capturing nonlinear relationships and interactions within the data, making it suitable for identifying intricate patterns that might influence RFS outcomes. Figure 3A shows the results of XGBoost and presents the top 15 important variables, including age, GGT, cholesterol, PTA, bile acid, prealbumin, DBIL, albumin, AST, PLT, Fib, globulin, creatinine, Lym and potassium.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g003.jpg)
Figure 3 Variable selection process. (A) Variable selection for RFS using eXtreme Gradient Boosting. (B) Variable selection for RFS using random survival forest. (C) Venn diagram showing the intersection of eXtreme Gradient Boosting and random survival forest. RFS, recurrence-free survival; GGT, gamma glutamyl transpeptidase; PTA, prothrombin time activity; DBIL, direct bilirubin; AST, aspartate aminotransferase; PLT, platelet; Fib, fibrinogen; Lym, lymphocyte; PTR, prothrombin time ratio; INR, international normalized ratio; PT, prothrombin time; PTA, prothrombin time activity.
We then employed another machine learning technique - random survival forest (RSF), to scrutinize variables affecting RFS. Figure 3B depicts the outcomes of the RSF analysis, showcasing the top 15 significant variables, which encompass tumor number, age, GGT, PTR, INR, tumor size, PT, PTA, gender, uric acid, albumin, AST, creatinine, Child-Pugh class and glucose.
Subsequently, we took the intersection of the top 15 important variables for each of these two methods, as shown in Figure 3C. It can be seen that age, AST, albumin, GGT, creatinine and PTA are their shared variables. This strategic approach enabled us to pinpoint variables that demonstrated consistent importance across different analytical frameworks. By focusing on the shared variables identified through this process, we ensured a more robust and reliable selection of predictors for further analysis.
Following that, we engaged in a rigorous screening process using multivariate Cox regression analysis. This statistical technique allowed us to assess the significance of each variable while controlling for potential confounding factors. Variables with a significance level below 0.05 were deemed statistically significant and identified as independent risk factors influencing RFS outcomes. Table 2 shows the results of multivariate Cox analysis. We finally identified age, AST, and PTA as independent risk factors. This meticulous approach ensured that only the most relevant and influential variables were considered in our analysis.
Nomogram developed for RFS prediction using the screened factors
Through the above rigorous statistical analysis, factors (age, AST and PTA) significantly associated with RFS were identified and incorporated into the nomogram (Figure 4). This nomogram serves as a practical aid in clinical settings, providing a visual representation of the identified factors’ contributions to the overall risk assessment for RFS. The basic principle of the nomogram, simply put, is to assign scores to each value of each influencing factor based on the contribution (i.e., the values of regression coefficients) of these factors to the RFS in the multivariate Cox regression model. These scores are then summed up to obtain a total score, which is used to predict the 1-year, 3-year and 5-year RFS probabilities of patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g004.jpg)
Figure 4 Nomogram for 1-, 3-, and 5-year RFS prediction in HBV-related HCC patients with low preoperative PALBI score after TACE combined with ablation therapy. Abbreviations: RFS, recurrence-free survival; HBV, hepatitis B virus; HCC, hepatocellular carcinoma; PALBI, platelet-albumin-bilirubin; TACE, transcatheter arterial chemoembolization; AST, aspartate aminotransferase; PTA, prothrombin time activity.
The Kaplan-Meier curves of low and high-risk groups stratified by nomogram scores
In both the training and validation cohorts, patients were categorized into two risk groups—low-risk and high-risk—based on the scores generated by the nomogram. Subsequently, Kaplan-Meier curves were plotted for both the low and high-risk groups to assess their RFS outcomes. In the training cohort, the Kaplan-Meier curves revealed a clear demarcation between the low-risk and high-risk groups (Figure 5). Patients classified as low-risk exhibited a more favorable prognosis, with a higher likelihood of survival over the follow-up period. Conversely, those categorized as high-risk faced a greater risk of recurrence, as evidenced by a steeper decline in survival probability. Similarly, the validation cohort exhibited comparable trends, further validating the prognostic utility of the nomogram-based risk stratification (Supplementary Figure S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g005.jpg)
Figure 5 Kaplan-Meier curves of different risk groups stratified by nomogram-derived points in the training cohort. RFS, recurrence-free survival.
Evaluating the performance of the nomogram
The evaluation of the nomogram’s performance is crucial for assessing its accuracy and reliability in predicting outcomes. Several metrics can be employed to gauge its effectiveness. Harrell’s concordance index (C-index) and Area Under the Curve (AUC) were used to evaluate its discrimination ability. Calibration curve was used to assess the agreement between predicted probabilities generated by the nomogram and observed outcomes. Decision curve analysis (DCA) was employed to evaluate its clinical utility. In the training cohort, we found that the nomogram exhibited strong discrimination ability. Its C-index was 0.689 (95% CI: 0.64–0.738) and its AUC values for 1-year, 3-year and 5-year RFS were 0.684, 0.747 and 0.813, respectively (Figure 6). The 1-year, 3-year and 5-year calibration curves closely approximated the ideal diagonal line, indicating consistent predictions with actual observations (Figure 7). The results of 1-year, 3-year and 5-year DCA curves demonstrated significant utility and clinical relevance (Figure 8). The analysis indicated that the nomogram’s predictions offer substantial net benefit across a range of thresholds, highlighting its potential for guiding clinical decision-making effectively. In the validation cohort, the nomogram also demonstrated excellent predictive efficacy, receiving favorable evaluations across these metrics. The C-index of the nomogram in the validation cohort was 0.642 (95% CI: 0.569–0.714). Its 1-year, 3-year and 5-year AUC values were 0.689, 0.716 and 0.747, respectively (Supplementary Figure S2). The calibration (Supplementary Figure S3) and DCA (Supplementary Figure S4) curves both showcased outstanding performance, affirming the model’s accuracy and significant clinical utility.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g006.jpg)
Figure 6 Receiver operating characteristic (ROC) curves in the training cohort. AUC, area under the curve.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g007.jpg)
Figure 7 Calibration curves in the training cohort. (A) Calibration curve of 1-year RFS. (B) Calibration curve of 3-year RFS. (C) Calibration curve of 5-year RFS. RFS, recurrence-free survival.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1409443/fimmu-15-1409443-HTML/image_m/fimmu-15-1409443-g008.jpg)
Figure 8 Decision curve analysis (DCA) in the training cohort. (A) DCA curve of 1-year RFS. (B) DCA curve of 3-year RFS. (C) DCA curve of 5-year RFS. RFS, recurrence-free survival.
Discussion
The absence of a specialized prognostic prediction model for low PALBI HBV-related HCC patients after receiving TACE combined with ablation therapy poses a considerable challenge in clinical management. To address this gap, our team embarked on comprehensive research aimed at developing a tailored prognostic nomogram specifically designed for this patient subgroup. Compared with past HCC prognostic models, this study has several significant innovations. Our study distinctively focuses on low PALBI score patients, a subgroup typically overlooked in previous research. Instead of treating all patients uniformly, we meticulously explore the unique characteristics and prognostic factors of these individuals, providing more targeted guidance for personalized treatment strategies. Furthermore, we employed cutting-edge machine learning algorithms to develop our prognostic prediction model. Unlike conventional statistical approaches, these algorithms excel in uncovering intricate patterns within data, thereby enhancing the accuracy and stability of our predictions. Through utilizing these advanced technologies, our nomogram offers a more precise recurrence evaluation for low PALBI HBV-related HCC patients.
The nomogram of this study included age, AST, and PTA as predictive indicators. We found that with increasing age, patients with HCC are at a higher risk of recurrence. As age increases, the functionality of the human immune system gradually declines, potentially leading to weakened immune surveillance against tumors and consequently HCC development (18). Furthermore, age-related comorbidities may also affect the prognosis of HCC patients. However, not all studies uniformly support the association between age and HCC development. Some studies have found that, under the same treatment regimens, there is no significant difference in the prognosis of HCC between elderly patients and younger patients (19–22). AST also plays a crucial role in the recurrence of HCC in this study. The elevated serum AST in patients is positively correlated with the degree of liver cell damage. Some scholars have found that persistently elevated AST levels may be associated with poorer survival rates and disease progression (23–26). Therefore, regular monitoring of AST levels may help predict the prognosis of HCC patients. PTA is an indicator used to evaluate coagulation function (27). PTA (%) = (normal person’s prothrombin time - 8.7)/(patient’s prothrombin time - 8.7) * 100. Prothrombin is one of the vitamin K-dependent coagulation factors synthesized by the liver, so it is essential to measure its activity in liver function tests. It is a sensitive indicator for assessing the severity of hepatocyte necrosis and prognosis. When the liver is damaged or impaired, as in conditions like hepatitis, cirrhosis and HCC, the synthesis of prothrombin is compromised, leading to prolonged patient’s prothrombin time and decreased PTA levels. Many studies have pointed out that prolonged prothrombin time is associated with higher morbidity and mortality in patients (28–30), while few studies have focused on the more precise expression of prothrombin time, PTA. This suggests that by monitoring patients’ preoperative PTA levels, it is possible to assess treatment risks in advance and take corresponding preventive measures to improve the success rate of treatment and patient survival.
The future of HCC prognosis models is poised at an exciting convergence of biotechnology, data science, and patient-centered care, heralding a new era of precision oncology. Advances in genomic sequencing and bioinformatics are enabling more nuanced understandings of HCC’s molecular and genetic underpinnings, driving the development of prognosis models that incorporate multi-omics data—genomic, transcriptomic, proteomic, and metabolomic—to offer predictions tailored to the individual genetic and molecular profile of the patient (28, 31–34). This multi-omics approach promises to refine risk stratification, guide therapy selection, and predict therapeutic responses with unprecedented accuracy. Parallelly, artificial intelligence (AI), deep learning (DL) and machine learning (ML) technologies are transforming the landscape of prognostic modeling (35–38). By harnessing vast datasets of patient records, imaging studies, and omics data, AI algorithms can identify patterns and prognostic markers that are imperceptible to the human eye. These AI-driven models are expected to improve continuously through iterative learning, becoming more precise and reliable as they are fed more data. The integration of AI with electronic health records and digital health technologies will further personalize prognosis and treatment plans, incorporating real-time health data from wearable devices and patient-reported outcomes to dynamically adjust prognostic predictions and therapeutic recommendations. Furthermore, the integration of liquid biopsy technologies into HCC prognosis models is a trend with significant potential. By enabling the non-invasive detection of circulating tumor DNA (ctDNA) and circulating tumor cells, liquid biopsies offer a means to monitor tumor burden, detect minimal residual disease, and identify emerging resistance mutations in real time (39–43). This technology could revolutionize HCC management, allowing for earlier intervention, the adaptation of treatment strategies in response to tumor evolution, and improved prognostic accuracy.
Therefore, future models are likely to incorporate patient preferences, quality of life measures, and socio-economic factors alongside traditional clinical and molecular data, facilitating shared decision-making and tailored treatment plans that align with patient values and life goals. In addition, in future research, we will seek to cooperate with other medical institutions to conduct multi-center studies to further verify the effectiveness of the nomogram we developed.
Conclusion
Our study successfully constructed a robust nomogram, effectively predicting 1-, 3-, and 5-year RFS for HBV-related HCC patients with low preoperative PALBI scores after TACE combined with local ablation therapy.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by Ethics Committee of Beijing You’an Hospital, affiliated with Capital Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because Given the retrospective nature of the study, where data was collected from past medical records, and the stringent measures implemented to protect patient’s privacy, the need for obtaining informed consent from the subjects was deemed unnecessary.
Author contributions
QW: Writing – review & editing, Methodology, Data curation, Formal Analysis, Writing – original draft. SS: Formal Analysis, Methodology, Writing – original draft, Writing – review & editing. YX: Writing – review & editing, Data curation, Validation. MH: Writing – review & editing, Conceptualization, Methodology, Project administration. RJ: Conceptualization, Methodology, Project administration, Writing – review & editing. CH: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by Beijing Bethune Charitable Foundation (QZHX-21-ZQN-014).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1409443/full#supplementary-material
Supplementary Figure 1 | Kaplan-Meier curves of different risk groups stratified by nomogram-derived points in the validation cohort. RFS, recurrence-free survival.
Supplementary Figure 2 | Receiver operating characteristic (ROC) curves in the validation cohort. AUC, area under the curve.
Supplementary Figure 3 | Calibration curves in the validation cohort. (A) Calibration curve of 1-year RFS. (B) Calibration curve of 3-year RFS. (C) Calibration curve of 5-year RFS. RFS, recurrence-free survival.
Supplementary Figure 4 | Decision curve analysis (DCA) in the validation cohort. (A) DCA curve of 1-year RFS. (B) DCA curve of 3-year RFS. (C) DCA curve of 5-year RFS. RFS, recurrence-free survival.
Abbreviations
HBV, hepatitis B virus; HCC, hepatocellular carcinoma; RFS, recurrence-free survival; PALBI, platelet-albumin-bilirubin; TACE, transarterial chemoembolization; XGBoost, eXtreme Gradient Boosting; RSF, random survival forest; AST, aspartate aminotransferase; PTA, prothrombin time activity; AUC, area under the curve; DCA, decision curve analysis; ALBI, albumin-bilirubin; AASLD, American Association for the Study of Liver Disease; WBC, white blood cell; RBC, red blood cell; Lym, lymphocyte; Neu, neutrophil; Mon, monocyte; PLT, platelet; Hb, hemoglobin; ALT, alanine aminotransferase; AST, aspartate aminotransferase; GGT, gamma glutamyl transpeptidase; ALP, alkaline phosphatase; TP, total protein; TBIL, total bilirubin; DBIL, direct bilirubin; Fib, fibrinogen; TT, thrombin time; INR, international normalized ratio; APTT, activated partial thromboplastin time; APTTR, activated partial thromboplastin time ratio; PT, prothrombin time; PTA, prothrombin time activity; PTR, prothrombin time ratio; BCLC stage, Barcelona Clinic Liver Cancer stage; CT, computed tomography; MRI, magnetic resonance imaging.
References
1. Rumgay H, Arnold M, Ferlay J, Lesi O, Cabasag CJ, Vignat J, et al. Global burden of primary liver cancer in 2020 and predictions to 2040. J Hepatol. (2022) 77:1598–606. doi: 10.1016/j.jhep.2022.08.021
2. Lin CL, Kao JH. Development of hepatocellular carcinoma in treated and untreated patients with chronic hepatitis B virus infection. Clin Mol Hepatol. (2023) 29:605–22. doi: 10.3350/cmh.2022.0342
3. Llovet JM, Kelley RK, Villanueva A, Singal AG, Pikarsky E, Roayaie S, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. (2021) 7:6. doi: 10.1038/s41572-020-00240-3
4. Tabrizian P, Jibara G, Shrager B, Schwartz M, Roayaie S. Recurrence of hepatocellular cancer after resection: patterns, treatments, and prognosis. Ann Surg. (2015) 261:947–55. doi: 10.1097/SLA.0000000000000710
5. Kim J, Kang W, Sinn DH, Gwak GY, Paik YH, Choi MS, et al. Substantial risk of recurrence even after 5 recurrence-free years in early-stage hepatocellular carcinoma patients. Clin Mol Hepatol. (2020) 26:516–28. doi: 10.3350/cmh.2020.0016
6. Hoekstra LT, de Graaf W, Nibourg GA, Heger M, Bennink RJ, Stieger B, et al. Physiological and biochemical basis of clinical liver function tests: a review. Ann Surg. (2013) 257:27–36. doi: 10.1097/SLA.0b013e31825d5d47
7. Ho SY, Liu PH, Hsu CY, Chiou YY, Su CW, Lee YH, et al. Prognostic performance of ten liver function models in patients with hepatocellular carcinoma undergoing radiofrequency ablation. Sci Rep. (2018) 8:843. doi: 10.1038/s41598-018-19251-y
8. Ho SY, Liu PH, Hsu CY, Hsia CY, Su CW, Lee YH, et al. Comparison of twelve liver functional reserve models for outcome prediction in patients with hepatocellular carcinoma undergoing surgical resection. Sci Rep. (2018) 8:4773. doi: 10.1038/s41598-018-22923-4
9. Kumada T, Toyoda H, Tada T, Yasuda S, Tanaka J. Changes in background liver function in patients with hepatocellular carcinoma over 30 years: comparison of child-pugh classification and albumin bilirubin grade. Liver Cancer. (2020) 9:518–28. doi: 10.1159/000507933
10. Johnson PJ, Berhane S, Kagebayashi C, Satomura S, Teng M, Reeves HL, et al. Assessment of liver function in patients with hepatocellular carcinoma: a new evidence-based approach-the ALBI grade. J Clin Oncol. (2015) 33:550–8. doi: 10.1200/JCO.2014.57.9151
11. Demirtas CO, D’Alessio A, Rimassa L, Sharma R, Pinato DJ. ALBI grade: Evidence for an improved model for liver functional estimation in patients with hepatocellular carcinoma. JHEP Rep. (2021) 3:100347. doi: 10.1016/j.jhepr.2021.100347
12. Hiraoka A, Kumada T, Kudo M, Hirooka M, Tsuji K, Itobayashi E, et al. Albumin-bilirubin (ALBI) grade as part of the evidence-based clinical practice guideline for HCC of the Japan society of hepatology: A comparison with the liver damage and child-pugh classifications. Liver Cancer. (2017) 6:204–15. doi: 10.1159/000452846
13. Liu PH, Hsu CY, Hsia CY, Lee YH, Chiou YY, Huang YH, et al. ALBI and PALBI grade predict survival for HCC across treatment modalities and BCLC stages in the MELD Era. J Gastroenterol Hepatol. (2017) 32:879–86. doi: 10.1111/jgh.13608
14. Ni JY, Fang ZT, An C, Sun HL, Huang ZM, Zhang TQ, et al. Comparison of albumin-bilirubin grade, platelet-albumin-bilirubin grade and Child-Turcotte-Pugh class for prediction of survival in patients with large hepatocellular carcinoma after transarterial chemoembolization combined with microwave ablation. Int J Hyperthermia. (2019) 36:841–53. doi: 10.1080/02656736.2019.1646927
15. Liu R, Li R, Zhang M, Liu W, Li H, Li D. Prognostic value of platelet-albumin-bilirubin grade in child-pugh A and B patients with hepatocellular carcinoma: A meta-analysis. Front Oncol. (2022) 12:914997. doi: 10.3389/fonc.2022.914997
16. Liu W, Zhang L, Xin Z, Zhang H, You L, Bai L, et al. A promising preoperative prediction model for microvascular invasion in hepatocellular carcinoma based on an extreme gradient boosting algorithm. Front Oncol. (2022) 12:852736. doi: 10.3389/fonc.2022.852736
17. Liu HF, Lu Y, Wang Q, Lu YJ, Xing W. Machine learning-based CEMRI radiomics integrating LI-RADS features achieves optimal evaluation of hepatocellular carcinoma differentiation. J Hepatocell Carcinoma. (2023) 10:2103–15. doi: 10.2147/JHC.S434895
18. Wang Y, Dong C, Han Y, Gu Z, Sun C. Immunosenescence, aging and successful aging. Front Immunol. (2022) 13:942796. doi: 10.3389/fimmu.2022.942796
19. Lee CR, Lim JH, Kim SH, Ahn SH, Park YN, Choi GH, et al. A comparative analysis of hepatocellular carcinoma after hepatic resection in young versus elderly patients. J Gastrointest Surg. (2012) 16:1736–43. doi: 10.1007/s11605-012-1966-7
20. Wang CM, Chen ZX, Ma PC, Chen JM, Jiang D, Hu XY, et al. Oncological prognosis and morbidity of hepatectomy in elderly patients with hepatocellular carcinoma: a propensity score matching and multicenter study. BMC Surg. (2023) 23:323. doi: 10.1186/s12893-023-02230-0
21. Chu KKW, Chok KSH. Is the treatment outcome of hepatocellular carcinoma inferior in elderly patients? World J Gastroenterol. (2019) 25:3563–71. doi: 10.3748/wjg.v25.i27.3563
22. Chen CH, Chang TT, Cheng KS, Su WW, Yang SS, Lin HH, et al. Do young hepatocellular carcinoma patients have worse prognosis? The paradox of age as a prognostic factor in the survival of hepatocellular carcinoma patients. Liver Int. (2006) 26:766–73. doi: 10.1111/j.1478-3231.2006.01309.x
23. von Felden J, Wege H, Schulze K. Elevated aspartate aminotransferase to alanine aminotransferase ratio predicts poor outcome in hepatocellular carcinoma. Hepatol Commun. (2020) 4:1382–3. doi: 10.1002/hep4.1570
24. Wan G, Gao F, Chen J, Li Y, Geng M, Sun L, et al. Nomogram prediction of individual prognosis of patients with hepatocellular carcinoma. BMC Cancer. (2017) 17:91. doi: 10.1186/s12885-017-3062-6
25. Xu L, Peng ZW, Chen MS, Shi M, Zhang YJ, Guo RP, et al. Prognostic nomogram for patients with unresectable hepatocellular carcinoma after transcatheter arterial chemoembolization. J Hepatol. (2015) 63:122–30. doi: 10.1016/j.jhep.2015.02.034
26. Ma L, Deng K, Zhang C, Li H, Luo Y, Yang Y, et al. Nomograms for predicting hepatocellular carcinoma recurrence and overall postoperative patient survival. Front Oncol. (2022) 12:843589. doi: 10.3389/fonc.2022.843589
27. Qin G, Bian ZL, Shen Y, Zhang L, Zhu XH, Liu YM, et al. Logistic regression model can reduce unnecessary artificial liver support in hepatitis B virus-associated acute-on-chronic liver failure: decision curve analysis. BMC Med Inform Decis Mak. (2016) 16:59. doi: 10.1186/s12911-016-0302-7
28. Li J, Xu Y, Liu H, Guo B, Guo X, Li Y, et al. Construction and validation of NSMC-ASIL, a predictive model for risk assessment of Malignant hepatic nodules. Am J Cancer Res. (2022) 12:5315–24.
29. Zhan J, Wang J, Zhang Z, Xue R, Jiang S, Liu J, et al. Noninvasive diagnosis of significant liver inflammation in patients with chronic hepatitis B in the indeterminate phase. Virulence. (2023) 14:2268497. doi: 10.1080/21505594.2023.2268497
30. Wang YY, Xiang BD, Ma L, Zhong JH, Ye JZ, Wang K, et al. Development and validation of a nomogram to preoperatively estimate post-hepatectomy liver dysfunction risk and long-term survival in patients with hepatocellular carcinoma. Ann Surg. (2021) 274:e1209–e17. doi: 10.1097/SLA.0000000000003803
31. Li M, Huang J, Zhan G, Li Y, Fang C, Xiang B. The novel-natural-killer-cell-related gene signature predicts the prognosis and immune status of patients with hepatocellular carcinoma. Int J Mol Sci. (2023) 24. doi: 10.3390/ijms24119587
32. Wang Y, Zhang Y, Wang L, Zhang N, Xu W, Zhou J, et al. Development and experimental verification of a prognosis model for cuproptosis-related subtypes in HCC. Hepatol Int. (2022) 16:1435–47. doi: 10.1007/s12072-022-10381-0
33. Zhang Q, Lou Y, Yang J, Wang J, Feng J, Zhao Y, et al. Integrated multiomic analysis reveals comprehensive tumor heterogeneity and novel immunophenotypic classification in hepatocellular carcinomas. Gut. (2019) 68:2019–31. doi: 10.1136/gutjnl-2019-318912
34. Xu LX, He MH, Dai ZH, Yu J, Wang JG, Li XC, et al. Genomic and transcriptional heterogeneity of multifocal hepatocellular carcinoma. Ann Oncol. (2019) 30:990–7. doi: 10.1093/annonc/mdz103
35. Chaudhary K, Poirion OB, Lu L, Garmire LX. Deep learning-based multi-omics integration robustly predicts survival in liver cancer. Clin Cancer Res. (2018) 24:1248–59. doi: 10.1158/1078-0432.CCR-17-0853
36. Saillard C, Schmauch B, Laifa O, Moarii M, Toldo S, Zaslavskiy M, et al. Predicting survival after hepatocellular carcinoma resection using deep learning on histological slides. Hepatology. (2020) 72:2000–13. doi: 10.1002/hep.31207
37. Shi JY, Wang X, Ding GY, Dong Z, Han J, Guan Z, et al. Exploring prognostic indicators in the pathological images of hepatocellular carcinoma based on deep learning. Gut. (2021) 70:951–61. doi: 10.1136/gutjnl-2020-320930
38. Liu Z, Liu Y, Zhang W, Hong Y, Meng J, Wang J, et al. Deep learning for prediction of hepatocellular carcinoma recurrence after resection or liver transplantation: a discovery and validation study. Hepatol Int. (2022) 16:577–89. doi: 10.1007/s12072-022-10321-y
39. Xu RH, Wei W, Krawczyk M, Wang W, Luo H, Flagg K, et al. Circulating tumor DNA methylation markers for diagnosis and prognosis of hepatocellular carcinoma. Nat Mater. (2017) 16:1155–61. doi: 10.1038/nmat4997
40. Ye Q, Ling S, Zheng S, Xu X. Liquid biopsy in hepatocellular carcinoma: circulating tumor cells and circulating tumor DNA. Mol Cancer. (2019) 18:114. doi: 10.1186/s12943-019-1043-x
41. Li J, Han X, Yu X, Xu Z, Yang G, Liu B, et al. Clinical applications of liquid biopsy as prognostic and predictive biomarkers in hepatocellular carcinoma: circulating tumor cells and circulating tumor DNA. J Exp Clin Cancer Res. (2018) 37:213. doi: 10.1186/s13046-018-0893-1
42. Ahn JC, Teng PC, Chen PJ, Posadas E, Tseng HR, Lu SC, et al. Detection of circulating tumor cells and their implications as a biomarker for diagnosis, prognostication, and therapeutic monitoring in hepatocellular carcinoma. Hepatology. (2021) 73:422–36. doi: 10.1002/hep.31165
Keywords: hepatitis B virus (HBV), hepatocellular carcinoma (HCC), recurrence-free survival (RFS), extreme gradient boosting (Xgboost), random survival forest (RSF), nomogram
Citation: Wang Q, Sheng S, Xiong Y, Han M, Jin R and Hu C (2024) Machine learning-based model for predicting tumor recurrence after interventional therapy in HBV-related hepatocellular carcinoma patients with low preoperative platelet-albumin-bilirubin score. Front. Immunol. 15:1409443. doi: 10.3389/fimmu.2024.1409443
Received: 30 March 2024; Accepted: 14 May 2024;
Published: 28 May 2024.
Edited by:
Fanping Meng, Fifth Medical Center of the PLA General Hospital, ChinaCopyright © 2024 Wang, Sheng, Xiong, Han, Jin and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ming Han, hm299792458@126.com; Ronghua Jin, ronghuajin@ccmu.edu.cn; Caixia Hu, hucaixia1217@126.com
†These authors have contributed equally to this work