1. Introduction
Malignant hepatic lesions are common and are associated with high mortality rates. The role of radiation therapy (RT) has historically been limited to palliation due to the risk of radiation-induced liver disease (RILD), which is of significant concern in a population that often has precarious hepatic function prior to treatment. With advances in delivery techniques and radiobiological understanding, dose-escalated RT can safely be delivered to more individuals, but the population that will benefit from this treatment has not yet been established.
The gold standard treatment of malignant liver lesions is surgery; however, less than 20% of patients are eligible for resection due to medical comorbidities, poor performance status, unresectable tumors, or patient refusal [
1]. The less invasive local therapies of cryoablation, radiofrequency ablation (RFA), percutaneous ethanol injection, and transarterial chemoembolization (TACE) may offer a mild survival advantage, but are not always technically feasible [
2,
3]. RT is a non-invasive treatment option that has few clinical limitations and can be administered to most tumors and patients. Therefore, particularly as many referred patients have exhausted other options, there may be an inordinate pressure to treat as the only option. However, many patients may not survive long enough to benefit from radiation due to early death or progression of disease outside of the liver.
RT to the liver using modern techniques is generally well tolerated [
4]. From the patient perspective, there can be a temporary worsening of appetite and fatigue, but quality of life measures generally return to baseline by 3 months [
5]. Still, there remains a wide range in survival of patients in the literature, and it would be prudent to avoid irradiating patients who are unlikely to survive long enough to benefit from the treatment, thus saving them the potential short-term toxicities. The Child-Pugh score [
6], Barcelona Clinic Liver Cancer score (BCLC) [
1], Model for End-Stage Liver Disease (MELD) [
7], Albumin-Bilirubin model (ALBI) [
8,
9], and Cancer of the Liver Italian Program (CLIP) [
7] are measures that have been developed to predict patient outcomes specific to liver disease; however, none of these were developed for patients being assessed for RT and none focused on identifying patients who should not receive radiation.
The aim of this study is to develop a Clinical Prognostic Index (PILiR score) using easily and commonly obtained pre-treatment patient and tumor factors to improve patient selection for dose escalated radiation treatment of liver malignancies. In particular, the focus of the model development is to identify patients for whom high dose radiation should not be offered. This aim was achieved.
4. Discussion
Malignant liver lesions are common and a significant source of morbidity and mortality. While surgical resection is the gold standard, many patients are not candidates and alternative treatment options need to be considered. Radiotherapy is a valid alternative and is not limited by many factors that preclude other treatments such as percutaneous accessibility, adequate perfusion of the tumor or procedural risks; however, it is not always clear who is most likely to benefit from this treatment since many of these patients have a guarded prognosis. Indeed, radiation may be the last option for many patients and this may result in caregivers overtreating patients. An index using pre-treatment patient and tumor factors was created, allowing for improved patient selection by identifying patients who are more likely to experience early death and thereby not benefit from dose-escalated radiotherapy.
In this study, the presence of extrahepatic disease and poor liver function, as indicated by a CPS of B or C and low serum albumin, were significant predictors of patients who did not benefit from treatment due to their poor prognosis when considering all patients with malignant liver lesions (combined cohort). In patients with HCC, CPS, and tumor volume were important predictive factors. In patients with metastatic liver disease, serum albumin and primary disease site were predictive. By using these factors, the PiLIR score was created (
Figure 2) which has high specificity, allowing the clinician to ensure that patients with a short life expectancy do not receive unnecessary treatment. The tool is easy to use with routinely obtained clinical parameters that can provide immediate information at the time of consultation.
Validated tools for prognostication in patients with HCC exist, including the ALBI nomogram which uses liver function alone to predict survival [
8] and the BCLC staging classification which takes into account both patient factors (performance status, CPS) and disease factors (tumor size, number of tumors, vascular invasion, nodal spread and presence of extrahepatic disease) and suggests that patients with early stage disease benefit from curative intent with transplant, resection, or ablation; whereas those with intermediate or advanced disease benefit from chemoembolization or treatment on clinical trial; and those with terminal disease benefit from supportive care [
1]. Of these, only the ALBI has been validated for use in patients undergoing radiotherapy, but does not specifically identify patients in which radiation should likely be avoided [
8].
General nomograms have previously been created for patients undergoing liver SBRT. Kress et al. (2012) established a prognostic scoring system for use in patients with liver lesions undergoing SBRT where the presence of 2 or more lesions, active systemic disease and poor performance status were found to be predictive of decreased overall survival [
11]. This scoring system, however, is based on a small dataset of only 52 patients with heterogeneous disease histology. Furthermore, measures of hepatic function, which have previously been shown to correlate with patient outcomes [
1,
8], were not included in their analysis. Our nomogram established that factors such as CPS, serum albumin, total bilirubin, presence of ascites, GTV, presence of extrahepatic disease, previous liver resection, and primary colorectal carcinoma were predictive of OS.
Nomograms have also been created for patients with HCC undergoing radiation. Using the Surveillance, Epidemiology and End Results Program (SEER) database, Zhan et al, (2023) created a model for patients with HCC treated with external beam radiotherapy (EBRT) where tumor size, TNM (tumor, nodal, metastasis) staging, AFP, previous surgery, and chemotherapy were predictive of OS [
12]. Similar models have been developed [
13,
14,
15]. These studies found multiple predictive factors including multiple lesions [
14,
15], tumor size [
13,
14,
16,
17,
18,
19,
20,
21,
22,
23], T(tumor) stage [
14,
17,
18,
21], N(nodal) stage [
18], M(metastasis) stage [
17,
18,
19,
21,
24], grade [
17,
18,
23], tumor volume [
16,
23] American Joint Cancer Committee (AJCC) stages T3–4 [
19], CPS [
14,
21], ECOG PS [
20], ALBI grade [
20], Portal Vein Tumor Thrombus(PVTT) classification [
20], BCLC stage [
20], macrovascular invasion [
13,
19], fibrosis score [
17,
19], histological type [
18], presence of cirrhosis [
14], biologic equivalent dose [
13], AFP [
14,
15,
16,
18,
19,
20,
23,
26], C-reactive protein (CRP) [
25,
26] hemoglobin [
14,
26], platelets [
15], alkaline phosphate [
15,
26], certain treatments [
15,
17,
18,
20,
21,
22,
25], patient age [
14,
18,
19], sex [
18,
21], race [
18], year of diagnosis [
18], symptoms related to lymph node metastasis [
24], and uncontrolled intrahepatic disease [
24].
The SEER database was also retrospectively used to create a nomogram for colorectal liver metastasis, using 9800 patients for training: 4897 patients for internal validation, and 60 for external validation [
27]. This nomogram had an AUC of 0.811 and 0.727 for the internal and external cohorts, respectively. Many factors were found to be predictive of OS in this nomogram and others with liver metastases for pancreatic, esophageal, and non-small cell lung cancer in addition to colorectal [
25,
26,
27,
28,
29,
30,
31,
32].These studies also found primary site [
27,
28], T stage [
27,
32], N stage [
27,
32], metastasis of bone [
27,
28,
30,
32,
33], brain [
27,
32,
33], lung [
27,
28,
30,
32,
33], predictive of OS as well as multiple liver metastases [
29], grade [
28,
30,
32,
33], AJCC N status [
28], tumor size [
28,
30,
32], histological type [
28,
30,
33], chemotherapy [
27,
28,
31,
32,
33], carcinoembryonic antigen (CEA) [
28], poor differentiation [
29,
30], adenocarcinoma [
30], and bilobar liver distribution [
28]. Despite our current model incorporating a more heterogenous population than others mentioned above, the predictive AUC values are similar to theirs. The novelty of our nomogram is that it specifically aims to identify patients that would likely not benefit from radiation and other management options should be sought.
There are important considerations and limitations to this study. First, in the current study, median overall survival was 11.9 months in the HCC group and 9.4 months in the metastatic group. This is shorter than other studies looking at survival in patients with malignant liver lesions, where median overall survival ranges from 14.5 months to 34 months [
9,
37,
38,
39,
40]. However, patients included in our study were heavily pre-treated, with 126 out of 195 (64.6%) patients having undergone at least one previous treatment and 17 patients (8.7%) having undergone three or more previous treatment modalities. Therefore, the participants in this study likely exhibited a higher degree of treatment resistance and were in more advanced stages of disease. When compared against studies that included predominantly advanced HCC patients, such as in a trial of sorafenib in unresectable HCC where median overall survival was 9.3 months [
8] and 1 year survival was 29% in untreated patients [
1], our results were more in keeping with their survival data. Thus, clinicians who wish to use this nomogram should determine if the demographics of our patients are generalizable to theirs.
Second, this study was designed to create an index based on pre-treatment factors alone. We have published a separate analysis assessing the role of RT modality and dose distribution on prognosis, along with clinical nomograms for both HCC and liver metastasis patients for predicting overall survival [
40].
Third, the current study produced a model based on analysis of one of the larger published databases of North American patients with malignant liver lesions treated with radiotherapy [
41]. However, the dataset is quite heterogeneous with wide acceptance in parameters such as size of lesion, presence of thrombosis, and presence of extrahepatic disease. In this study, all metastatic patients with a primary disease site other than colorectal were grouped together, resulting in the pooling of diseases with varied prognoses. Other prognostic indices have stratified patients based on the primary disease site [
42] which allows for more homogenous data and, therefore, better predictive ability; however, in doing so the sample size decreases dramatically. A larger database will allow further stratification by primary disease site, allowing for improved predictive ability.
The PILiR scoring system is designed to identify patients that are unlikely to benefit from radiation and the goal was to have a high specificity to avoid incorrectly categorizing patients who would then likely receive only supportive care and be inappropriately excluded from the opportunity to receive a beneficial treatment. However, to avoid this type of error, this nomogram has low sensitivity. Compared to the total population, the number of patients living less than 4 months was relatively low, which resulted in a small number of patients in the higher point groups. This was especially true for the subgroup indices where there were only 3 patients in the HCC subgroup scoring 4 points and only 5 patients in the metastatic subgroup scoring 4 points. By analyzing a larger, more homogenous population, an index with improved specificity and sensitivity can be created. A larger population will also allow us to provide more focused information such as assessing value of treatment for oligometastases, oligoprogression and dominant lesion control. Currently, these patients are grouped into one category with varying goals of treatment.
Use of the PiLiR scoring system can allow for identification of patients who could be considered for aggressive radiotherapy for longer term local control versus patients who may instead benefit from a focus on short-term symptom control alone. By screening at consultation with this tool, patients who would best be served by palliation due to poor life expectancy, as indicated by a high PILiR score, can be identified. In such cases, radiotherapy with palliative dosing such as 800cGy delivered in a single fraction has been shown to improve patient-reported symptoms and quality of life at 1-month post-treatment [
34,
35]. On the other hand, those with low PILiR scores should be considered for aggressive radiation treatment with the goal of controlling tumor progression and future burden of disease.
A large, shared database across several cancer centers is being developed which will allow for further analysis and validation with higher power and reduced bias. External validation will allow for further confidence that the scoring system can provide clinicians with the tools to select patients appropriately and reliably for radiation.
Author Contributions
Conceptualization, L.C. and M.L.; methodology, L.C. M.L.; software, L.C.; validation, L.C., S.G., M.L.; formal analysis, L.C.; investigation, L.C.; resources, L.C., M.L.; data curation, L.C.; writing—original draft preparation, L.C., H.R., N.G., M.L.; writing—review and editing, L.C., H.R., N.G., S.G., E.W., D.T., J.V., J.P., and M.L.; visualization, L.C., H.R., N.G., J.P.; supervision, M.L.; project administration, L.C. All authors have read and agreed to the published version of the manuscript.