Background: The conversion therapy for advanced hepatocellular carcinoma (HCC) shows promise with a triple therapy approach that combines interventional therapy, immune checkpoint inhibitors, and molecular targeted therapy (primarily small-molecule TKIs and the large-molecule bevacizumab). This combination has achieved the highest objective response rates (ORR) along with acceptable safety profiles. The aim of this study was to compare the clinical efficacy of lenvatinib versus bevacizumab, when combined with immune checkpoint inhibitors and interventional triple therapy, as first-line treatments for Chinese patients with unresectable HCC (uHCC). Method: This retrospective multicenter study involved 371 consecutive patients from 21 centers in China, observed between April 2017 and December 2023. The study focused on patients with uHCC at Chinese liver cancer stages IIb to IIIb (Barcelona Clinic Liver Cancer stage B or C) who received lenvatinib or bevacizumab combined with anti-PD-1/L1 and interventional therapy (including TACE and/or HAIC) as first-line treatment. Of the 371 patients, 258 received lenvatinib-based triple therapy, while 113 received bevacizumab-based triple therapy. The primary endpoints were overall survival (OS) and progression-free survival (PFS). To balance baseline clinical characteristics, propensity score matching (PSM) and inverse probability of treatment weighting (IPTW) were applied. Subgroup analysis was also performed based on different clinicopathological characteristics of the enrolled uHCC patients. Results: The median OS in the lenvatinib group was significantly longer than in the bevacizumab group, both before (36.0 vs. 27.9 months; hazard ratio [HR]: 0.536; 95% confidence interval [CI]: 0.344–0.835; p = 0.0016) and after PSM (HR: 0.524; 95% CI: 0.305–0.900; p = 0.01), as well as after IPTW (HR: 0.549; 95% CI: 0.331–0.908; p = 0.01). Before adjustment, PFS in the lenvatinib group was also significantly longer than in the bevacizumab group (20.0 vs. 12.1 months; HR: 0.649; 95% CI: 0.457–0.922; p = 0.0078). However, after PSM (HR: 0.808; 95% CI: 0.535–1.222; p = 0.33) and IPTW, there was no significant difference in PFS between the two groups. Multivariate analysis showed that lenvatinib-based triple therapy was independently associated with improved OS compared to bevacizumab-based triple therapy. Subgroup analysis indicated that patients with age ≤65 years, no history of hepatitis B virus infection, Barcelona Clinic Liver Cancer stage C (BCLC-C), ALT levels ≤40 U/L, platelets ≥100 × 109/L, or log 10 AFP ≥1.40 benefited more from lenvatinib-based triple therapy. Conclusion: Lenvatinib-based triple therapy tends to prolong OS compared to bevacizumab, although the PFS was similar between the two groups. Patients aged ≤65 years, without a history of hepatitis B virus infection, with BCLC-C stage, ALT ≤40 U/L, platelets ≥100 × 109/L, or log 10 AFP ≥1.40 are likely to benefit more from lenvatinib-based triple therapy.

This study explores the effectiveness of two different treatments for advanced liver cancer, known as hepatocellular carcinoma (HCC), in patients who cannot have surgery. The treatments tested combined immune therapy with two drugs: lenvatinib and bevacizumab. These drugs work by targeting different aspects of cancer growth. Lenvatinib is a small molecule that blocks blood vessel growth, while bevacizumab is a larger molecule that also prevents the tumor from getting the nutrients it needs to grow. The research involved 371 patients from 21 hospitals in China, who received treatment between 2017 and 2020. The study compared how well lenvatinib and bevacizumab worked when combined with other therapies, like immune checkpoint inhibitors and interventional treatments (procedures used to treat cancer without surgery). Results showed that both treatments were effective, but the study aimed to determine which one was better in terms of survival and other health outcomes for patients with advanced HCC. The study’s findings could help doctors choose the most effective first-line treatment for these patients in the future.

Hepatocellular carcinoma (HCC) is globally recognized as the sixth most prevalent form of cancer and ranks as the third leading cause of cancer-related mortality [1], while it ranks fourth in incidence and second in mortality in China [2]. Most cases of HCC are diagnosed at a relatively advanced stage, making them ineligible for resection. Conversion therapy has been gaining attention in recent years. For unresectable HCC, systemic therapy options like single-agent lenvatinib and sorafenib are considered first-line treatments [3‒5]. However, over time, novel therapeutic alternatives have emerged, such as immune checkpoint inhibitors (ICIs) and ICIs or tyrosine kinase inhibitors (TKIs), such as atezolizumab plus bevacizumab (T+A) and tremelimumab plus durvalumab and nivolumab plus ipilimumab providing patients with advanced HCC a broader array of treatment options [6‒10].

In recent years, systemic therapy for liver cancer has achieved favorable results and has become an important means of conversion therapy. Among these, the efficacy of lenvatinib monotherapy has been proven to be no inferior to sorafenib, with the highest single-agent objective response rate (ORR) being 18.8% [11]. The efficacy rate of the T+A dual-drug regimen is approximately 30% [12] and has been widely used in clinical practice [13]. Although the combination of lenvatinib and pembrolizumab (the “Combo” regimen) showed negative results in the LEAP-002 study [14], there are also phase II studies (such as Keynote 524, ORR = 36%) [15] that have made this regimen a choice for many physicians in clinical practice. Nevertheless, combining local treatments such as radiotherapy, transarterial chemoembolization (TACE), and hepatic arterial infusion chemotherapy (HAIC) with systemic therapy can achieve an even higher conversion rate. For example, triple therapy using the small-molecule TKI lenvatinib and the large-molecule bevacizumab has received support from extensive phase II clinical research data, with an ORR of 56.7%–74.4% [3, 4, 16, 17]. Additionally, phase III clinical trials such as EMERALD-1 [18] and LEAP-012 [19] have already reported their results, and further analyses are ongoing.

Although lenvatinib and bevacizumab have achieved promising results in the treatment of unresectable HCC, the optimal combination regimen remains uncertain, particularly between lenvatinib-based and bevacizumab-based triple combinations. Direct comparisons between these treatments are lacking. Therefore, this study aimed to compare the prognostic impact of lenvatinib-based and bevacizumab-based triple therapy in patients with unresectable HCC. We believe that this comparison is critical in determining the best treatment options for patients with advanced HCC. To achieve this goal, we planned to conduct a retrospective cohort study using medical record data from patients with unresectable HCC who received triple therapy with either lenvatinib or bevacizumab. The primary outcome will be overall survival (OS) and the secondary outcome will be progression-free survival (PFS). Based on data from multiple centers across the country from the CLEAP research group [20‒25], we designed this study. We hope this study can provide valuable information to guide clinical practice and improve the prognosis of patients with unresectable HCC.

The CLEAP database contains data on 735 uHCC patients treated at 21 clinical centers in China between April 2017 and December 2023. Patients with the following features were included: (1) HCC diagnosed by clinical or pathological features according to the American Association for Study of Liver Diseases (AASLD) practice guidelines [26]; (2) belong to Chinese liver cancer (CNLC) stage IIb, IIIa, and IIIb [27], equivalent to Barcelona clinical liver cancer stage B or C; (3) received lenvatinib-based triple therapy as first-line treatment, or bevacizumab-based triple therapy as first-line treatment; and (4) follow-up data were complete and available. Exclusion criteria were as follows: diagnosis of other malignant tumors; use of interventional therapy other than TACE or HAIC; and use of systemic treatment regimens other than PD-1 or PD-L1 antibodies. Eligible patients were divided into the lenvatinib group and the bevacizumab group according to the treatment regimen adopted for analysis. Treatments were administered according to the product package insert and previous studies, based on current efficacy and safety data, prior treatment, and drug availability [28].

Use of both treatments followed the same inclusion criteria as registered trials. The current study was approved by the Ethics Committee of each center and adhered to clinical practice guidelines, the Declaration of Helsinki, and local laws. All patients provided informed consent to receive treatment for HCC and to have their medical records used for research purposes.

Treatment Overview

Patients were treated with lenvatinib/bevacizumab in combination with PD-1/PD-L1 antibodies and interventional therapies, including TACE and HAIC, as first-line treatment. The PD-1 antibodies used in this study included pembrolizumab, nivolumab, sintilimab, tislelizumab, toripalimab, and camrelizumab, while the PD-L1 antibodies included durvalumab and atezolizumab. Notably, 32.1% of patients initially underwent interventional therapy, while the remaining patients received a combination of the triple therapy regimen concurrently.

TACE is a minimally invasive procedure in which chemotherapy agents are directly delivered to the tumor through the hepatic artery. This is followed by embolization to block the tumor’s blood supply, thereby enhancing the local effectiveness of the drugs. The standard TACE regimen typically includes epirubicin and lipid-based platinum compounds. For cisplatin, the recommended dose is 50 mg/m2 of body surface area, administered at a rate of 4–10 mL/min. TACE treatments are generally performed every 3–6 weeks, with the frequency adjusted based on the patient’s clinical condition and response.

In contrast, HAIC involves the direct infusion of chemotherapy agents into the hepatic artery, achieving high local drug concentrations while minimizing systemic exposure. In this study, the HAIC regimen primarily consisted of the FOLFOX protocol, which includes oxaliplatin, leucovorin, and fluorouracil. Recommended doses for HAIC therapy are 85 mg/m2 for oxaliplatin, 400 mg/m2 for cisplatin, and 2,400 mg/m2 for fluorouracil, administered every 3–4 weeks. As with TACE, the frequency and dosage of HAIC treatments were tailored to the patient’s clinical response and condition.

Treatment with modality continued until disease progression, drug toxicity intolerance, or eligibility for therapeutic liver resection. In clinical practice in China, TACE is predominantly used for patients with CNLC stage IIa–IIb (equivalent to BCLC stage B) HCC, particularly those with multinodular tumors without vascular invasion. Conversely, HAIC is more commonly employed for patients with CNLC stage IIIa–IIIb (BCLC stage C) HCC, especially those with large tumors or portal vein invasion. Additionally, all patients with HBV infection received oral antiviral drugs, such as entecavir or tenofovir, as part of their cancer treatment.

Clinical Assessments

Tumor evaluation was conducted every 6–8 weeks using CT or MRI for intrahepatic tumors and upper abdominal metastases, chest CT for lung metastases, and blood tests to assess liver function and tumor markers. OS was measured from the initiation of treatment until death. PFS was measured from the start of treatment until disease progression or death from any cause. Safety was evaluated based on the occurrence of severe (grade ≥3) treatment-related adverse events (TRAEs).

Statistical Analysis

Categorical variables are summarized as counts (percentages), and continuous variables are summarized as means ± standard deviation. Differences in baseline characteristics between groups were estimated using Fisher’s exact test or Kruskal-Wallis test. The Kaplan-Meier method was used for survival analysis, and the log-rank test was used to compare OS and PFS curves. Survival outcomes were also analyzed by Cox regression, using univariable and multivariable models that included relevant variables. Subgroup analysis was conducted based on age (≤65 or >65 years), gender, HBV infection status, BCLC stage, CNLC stage, ALT level (<40 U/L or ≥40 U/L), platelet level (<100 × 109 or ≥100 × 109/L), and AFP level, with logarithms taken based on commonly used clinical standard values in China (log value <1.40 or ≥1.40).

In order to eliminate differences in potential confounding variables, logistic regression was used to calculate the propensity score (PS), and 1:3 matching was performed on the nearest individuals based on the PS, with the caliper value set to 0.05. Propensity score matching (PSM) was conducted using the “Matchit” package in R software, with covariates including HBV infection status, ECOG PS, Child-Pugh liver function class, and other variables. Sensitivity analysis using inverse probability of treatment weighting was performed on the same variables using the “RISCA” package in R software. All statistical analyses were performed using R software version 4.3.1 (http://www.rproject.org/). p value <0.05 was considered statistically significant.

Baseline Characteristics

This analysis included 371 consecutive patients who received lenvatinib or bevacizumab-based triple therapy between April 1, 2017, and December 31, 2023. The patients’ disposition is shown in Figure 1. In the overall population (Table 1), the mean age was 55.3 years (standard deviation, 10.7), and the majority were male (85%). Most patients had hepatitis B virus (HBV) infection (HBsAg, HBcAb, or HBeAb seropositivity) as the cause of uHCC (85%), and most had ECOG PS 0 (69%) and Child-Pugh class A liver function (87%). Based on CNLC criteria, 24%, 37%, and 39% of patients were diagnosed as stage IIb, stage IIIa, and stage IIIb, respectively.

Fig. 1.

Flowchart of patients selection. CLEAP, China Liver Study Group Young Investigators; CNLC, China liver cancer; ECOG, Eastern, Eastern Cooperative Oncology Group; uHCC, unresectable hepatocellular carcinoma; PSM, propensity score matching; IPTW, inverse probability of treatment weighting.

Fig. 1.

Flowchart of patients selection. CLEAP, China Liver Study Group Young Investigators; CNLC, China liver cancer; ECOG, Eastern, Eastern Cooperative Oncology Group; uHCC, unresectable hepatocellular carcinoma; PSM, propensity score matching; IPTW, inverse probability of treatment weighting.

Close modal
Table. 1.

Baseline demographics and disease characteristics and outcomes of patients before and after PSM

VariableBefore PSMAfter PSM
Overall, N = 3711Beva, N = 1131Lenva, N = 2581SMDp value2Beva, N = 821Lenva, N = 1521SMDp value2
Age, years 55.34 (10.65) 57.11 (9.74) 54.56 (10.95) 0.25 0.051 57.23 (9.96) 54.34 (11.16) 0.27 0.080 
Gender    0.10 0.4   0.03 0.8 
 Female 56/371 (15%) 20/113 (18%) 36/258 (14%)   12/82 (15%) 24/152 (16%)   
 Male 315/371 (85%) 93/113 (82%) 222/258 (86%)   70/82 (85%) 128/152 (84%)   
Intervention    0.10 0.4   0.09 0.5 
 ≤6 times 342/371 (92%) 102/113 (90%) 240/258 (93%)   74/82 (90%) 141/152 (93%)   
 >6 times 29/371 (7.8%) 11/113 (9.7%) 18/258 (7.0%)   8/82 (9.8%) 11/152 (7.2%)   
HBV    0.25 0.021   0.05 0.7 
 Infected 316/371 (85%) 89/113 (79%) 227/258 (88%)   68/82 (83%) 129/152 (85%)   
 Uninfected 55/371 (15%) 24/113 (21%) 31/258 (12%)   14/82 (17%) 23/152 (15%)   
ECOG PS    −0.87 <0.001   −0.28 0.13 
 0 257/371 (69%) 107/113 (95%) 150/258 (58%)   76/82 (93%) 126/152 (83%)   
 1 107/371 (29%) 5/113 (4.4%) 102/258 (40%)   5/82 (6.1%) 23/152 (15%)   
 2 6/371 (1.6%) 1/113 (0.9%) 5/258 (1.9%)   1/82 (1.2%) 2/152 (1.3%)   
 3 1/371 (0.3%) 0/113 (0%) 1/258 (0.4%)   0/82 (0%) 1/152 (0.7%)   
ALBI grade    −0.11 0.6   −0.09 0.7 
 Grade 1 144/371 (39%) 48/113 (42%) 96/258 (37%)   39/82 (48%) 68/152 (45%)   
 Grade 2 216/371 (58%) 62/113 (55%) 154/258 (60%)   42/82 (51%) 79/152 (52%)   
 Grade 3 11/371 (3.0%) 3/113 (2.7%) 8/258 (3.1%)   1/82 (1.2%) 5/152 (3.3%)   
Albumin 38.87 (11.64) 38.64 (5.95) 38.97 (13.40) −0.03 0.6 40.02 (5.44) 40.06 (16.94) 0.00 0.3 
Bilirubin 19.52 (13.76) 20.86 (16.29) 18.93 (12.49) 0.13 0.4 17.91 (9.61) 17.98 (11.30) −0.01 0.6 
ALT 52.14 (58.52) 57.24 (72.85) 49.90 (51.00) 0.12 >0.9 51.21 (47.64) 51.13 (59.17) 0.00 0.9 
AST 76.06 (86.56) 80.88 (116.40) 73.95 (69.78) 0.07 >0.9 65.55 (53.80) 73.89 (76.38) −0.13 0.7 
NLR 4.04 (5.00) 3.88 (3.19) 4.12 (5.61) −0.05 0.7 3.45 (1.80) 4.32 (6.49) −0.18 0.3 
PLT 187.19 (94.63) 176.82 (81.92) 191.74 (99.50) −0.16 0.3 177.77 (81.23) 204.17 (104.18) −0.28 0.13 
Log 10 AFP 2.61 (1.58) 2.55 (1.64) 2.63 (1.55) −0.05 0.5 2.53 (1.65) 2.65 (1.56) −0.08 0.5 
BCLC stage    0.02 0.9   0.06 0.7 
 B 90/371 (24%) 28/113 (25%) 62/258 (24%)   19/82 (23%) 39/152 (26%)   
 C 281/371 (76%) 85/113 (75%) 196/258 (76%)   63/82 (77%) 113/152 (74%)   
CNLC stage    0.28 0.053   0.07 0.9 
 IIb 88/371 (24%) 27/113 (24%) 61/258 (24%)   18/82 (22%) 38/152 (25%)   
 IIIa 137/371 (37%) 51/113 (45%) 86/258 (33%)   33/82 (40%) 59/152 (39%)   
 IIIb 146/371 (39%) 35/113 (31%) 111/258 (43%)   31/82 (38%) 55/152 (36%)   
Child-Pugh    0.34 0.001   0.21 0.14 
 A 321/371 (87%) 88/113 (78%) 233/258 (90%)   78/82 (95%) 136/152 (89%)   
 B 50/371 (13%) 25/113 (22%) 25/258 (9.7%)   4/82 (4.9%) 16/152 (11%)   
VariableBefore PSMAfter PSM
Overall, N = 3711Beva, N = 1131Lenva, N = 2581SMDp value2Beva, N = 821Lenva, N = 1521SMDp value2
Age, years 55.34 (10.65) 57.11 (9.74) 54.56 (10.95) 0.25 0.051 57.23 (9.96) 54.34 (11.16) 0.27 0.080 
Gender    0.10 0.4   0.03 0.8 
 Female 56/371 (15%) 20/113 (18%) 36/258 (14%)   12/82 (15%) 24/152 (16%)   
 Male 315/371 (85%) 93/113 (82%) 222/258 (86%)   70/82 (85%) 128/152 (84%)   
Intervention    0.10 0.4   0.09 0.5 
 ≤6 times 342/371 (92%) 102/113 (90%) 240/258 (93%)   74/82 (90%) 141/152 (93%)   
 >6 times 29/371 (7.8%) 11/113 (9.7%) 18/258 (7.0%)   8/82 (9.8%) 11/152 (7.2%)   
HBV    0.25 0.021   0.05 0.7 
 Infected 316/371 (85%) 89/113 (79%) 227/258 (88%)   68/82 (83%) 129/152 (85%)   
 Uninfected 55/371 (15%) 24/113 (21%) 31/258 (12%)   14/82 (17%) 23/152 (15%)   
ECOG PS    −0.87 <0.001   −0.28 0.13 
 0 257/371 (69%) 107/113 (95%) 150/258 (58%)   76/82 (93%) 126/152 (83%)   
 1 107/371 (29%) 5/113 (4.4%) 102/258 (40%)   5/82 (6.1%) 23/152 (15%)   
 2 6/371 (1.6%) 1/113 (0.9%) 5/258 (1.9%)   1/82 (1.2%) 2/152 (1.3%)   
 3 1/371 (0.3%) 0/113 (0%) 1/258 (0.4%)   0/82 (0%) 1/152 (0.7%)   
ALBI grade    −0.11 0.6   −0.09 0.7 
 Grade 1 144/371 (39%) 48/113 (42%) 96/258 (37%)   39/82 (48%) 68/152 (45%)   
 Grade 2 216/371 (58%) 62/113 (55%) 154/258 (60%)   42/82 (51%) 79/152 (52%)   
 Grade 3 11/371 (3.0%) 3/113 (2.7%) 8/258 (3.1%)   1/82 (1.2%) 5/152 (3.3%)   
Albumin 38.87 (11.64) 38.64 (5.95) 38.97 (13.40) −0.03 0.6 40.02 (5.44) 40.06 (16.94) 0.00 0.3 
Bilirubin 19.52 (13.76) 20.86 (16.29) 18.93 (12.49) 0.13 0.4 17.91 (9.61) 17.98 (11.30) −0.01 0.6 
ALT 52.14 (58.52) 57.24 (72.85) 49.90 (51.00) 0.12 >0.9 51.21 (47.64) 51.13 (59.17) 0.00 0.9 
AST 76.06 (86.56) 80.88 (116.40) 73.95 (69.78) 0.07 >0.9 65.55 (53.80) 73.89 (76.38) −0.13 0.7 
NLR 4.04 (5.00) 3.88 (3.19) 4.12 (5.61) −0.05 0.7 3.45 (1.80) 4.32 (6.49) −0.18 0.3 
PLT 187.19 (94.63) 176.82 (81.92) 191.74 (99.50) −0.16 0.3 177.77 (81.23) 204.17 (104.18) −0.28 0.13 
Log 10 AFP 2.61 (1.58) 2.55 (1.64) 2.63 (1.55) −0.05 0.5 2.53 (1.65) 2.65 (1.56) −0.08 0.5 
BCLC stage    0.02 0.9   0.06 0.7 
 B 90/371 (24%) 28/113 (25%) 62/258 (24%)   19/82 (23%) 39/152 (26%)   
 C 281/371 (76%) 85/113 (75%) 196/258 (76%)   63/82 (77%) 113/152 (74%)   
CNLC stage    0.28 0.053   0.07 0.9 
 IIb 88/371 (24%) 27/113 (24%) 61/258 (24%)   18/82 (22%) 38/152 (25%)   
 IIIa 137/371 (37%) 51/113 (45%) 86/258 (33%)   33/82 (40%) 59/152 (39%)   
 IIIb 146/371 (39%) 35/113 (31%) 111/258 (43%)   31/82 (38%) 55/152 (36%)   
Child-Pugh    0.34 0.001   0.21 0.14 
 A 321/371 (87%) 88/113 (78%) 233/258 (90%)   78/82 (95%) 136/152 (89%)   
 B 50/371 (13%) 25/113 (22%) 25/258 (9.7%)   4/82 (4.9%) 16/152 (11%)   

PSM, propensity score matching; SMD, standardized mean difference; IQR, interquartile range; SD, standard deviation; HCC, hepatocellular carcinoma; HBV, hepatitis B virus; ECOG, Eastern Cooperative Oncology Group; ALBI, albumin-bilirubin grade, BCLC, Barcelona Clinic Liver Cancer; AFP, alpha-fetoprotein.

1Mean (SD); n/N (%).

2Wilcoxon rank sum test, Pearson’s Chi-squared test, and Fisher’s exact test.

Although most baseline variables were comparable between the two groups, significantly more patients were infected with HBV (88% vs. 79%; p = 0.021) and had ECOG PS 1 (40% vs. 4.4%; p < 0.001) in the lenvatinib group compared with the bevacizumab group. There were also significantly more patients with Child-Pugh class A liver function (90% vs. 78%; p = 0.001) in the lenvatinib group.

Primary Outcomes

As of December 2023, a total of 102 (27.5%) patients had died in the overall population. The median OS for all patients was 32.3 months (95% CI, 27.9–36.0 months), and the median PFS was 17.7 months (95% CI, 13.9–22.1 months).

During the follow-up period, 62 patients (24.0%) died in the lenvatinib group and 40 patients (35.4%) died in the bevacizumab group. Compared with the bevacizumab group, the median OS in the lenvatinib group was significantly longer (36.0 months vs. 27.9 months, HR, 0.536; 95% CI, 0.344–0.835; p = 0.0016; Fig. 2a). The OS rates in the lenvatinib group at 12, 24, and 36 months were 85.8%, 62.5%, and 46.5%, respectively. Additionally, the PFS in the lenvatinib group was also significantly longer than that in the bevacizumab group (20.0 months vs. 12.1 months, HR, 0.649; 95% CI, 0.457–0.922; p = 0.0078; Fig. 2b). The PFS rates at 12, 24, and 36 months in the lenvatinib group were 64.2%, 44.0%, and 31.1%, respectively.

Fig. 2.

Kaplan-Meier curves of OS and PFS for patients in the different groups. OS (a) and PFS (b) of patients in the lenvatinib group (N = 258) or the bevacizumab group (N = 113) before PSM. OS (c) and PFS (d) for the lenvatinib group (N = 152) and the bevacizumab group (N = 82) after PSM. OS (e) and PFS (f) for patients in the lenvatinib group (N = 234) and the bevacizumab group (N = 237) weighted by IPTW. PSM, propensity score matching; IPTW, inverse probability of treatment weighting.

Fig. 2.

Kaplan-Meier curves of OS and PFS for patients in the different groups. OS (a) and PFS (b) of patients in the lenvatinib group (N = 258) or the bevacizumab group (N = 113) before PSM. OS (c) and PFS (d) for the lenvatinib group (N = 152) and the bevacizumab group (N = 82) after PSM. OS (e) and PFS (f) for patients in the lenvatinib group (N = 234) and the bevacizumab group (N = 237) weighted by IPTW. PSM, propensity score matching; IPTW, inverse probability of treatment weighting.

Close modal

Multivariate analysis showed that the lenvatinib group had improved OS (HR, 0.51; 95% CI, 0.31–0.84; p = 0.009). In addition, aspartate aminotransferase (AST) ≥40 U/L indicated a poor prognosis (HR, 2.03; 95% CI, 1.13–3.63; p = 0.017) (Table 2).

Table. 2.

The Cox regression analysis is for OS

CharacteristicOS
Univariate analysisMultivariate analysis
p valueHR95% CIp value
Treatment (Lenva vs. Beva) 0.010 0.51 0.31, 0.84 0.009 
Age (>65 vs. ≤65) 0.38    
Gender (male vs. female) 0.54    
HBV (uninfected vs. infected) 0.76    
ALT (≥40 U/L vs. <40 U/L) 0.94    
AST (≥40 U/L vs. <40 U/L) 0.003 2.03 1.13, 3.63 0.017 
PT (≥13 s vs. <13 s) 0.38    
NLR (≥4.00 vs. <4.00) 0.40    
PLT (≥100 × 109 vs. <100 × 1090.14    
Log 10 AFP (≥1.40 vs. <1.40) 0.015 1.95 1.00, 3.80 0.050 
BCLC stage (C vs. B) 0.10    
CNLC stage (IIIa vs. IIb) 0.22    
CNLC stage (IIIb vs. IIb) 0.092    
Child-Pugh class (B vs. A) 0.35    
CharacteristicOS
Univariate analysisMultivariate analysis
p valueHR95% CIp value
Treatment (Lenva vs. Beva) 0.010 0.51 0.31, 0.84 0.009 
Age (>65 vs. ≤65) 0.38    
Gender (male vs. female) 0.54    
HBV (uninfected vs. infected) 0.76    
ALT (≥40 U/L vs. <40 U/L) 0.94    
AST (≥40 U/L vs. <40 U/L) 0.003 2.03 1.13, 3.63 0.017 
PT (≥13 s vs. <13 s) 0.38    
NLR (≥4.00 vs. <4.00) 0.40    
PLT (≥100 × 109 vs. <100 × 1090.14    
Log 10 AFP (≥1.40 vs. <1.40) 0.015 1.95 1.00, 3.80 0.050 
BCLC stage (C vs. B) 0.10    
CNLC stage (IIIa vs. IIb) 0.22    
CNLC stage (IIIb vs. IIb) 0.092    
Child-Pugh class (B vs. A) 0.35    

HR, hazard ratio; CI, confidence interval.

Adjustment

PSM processing was performed on variables such as HBV infection, ECOG PS, and Child-Pugh liver function classification. To retain as many enrolled cases as possible, the matching ratio was set to 1:3. After matching, all baseline variables between the two treatment groups were comparable (all p > 0.05) (Table 1). Notably, after PSM adjustment, the OS in the lenvatinib group was significantly prolonged, with OS rates at 12, 24, and 36 months of 86.6%, 64.5%, and 55.3%, respectively, compared to 65.5%, 55.4%, and 49.9% in the bevacizumab group (HR, 0.524; 95% CI, 0.305–0.900; p = 0.01; Fig. 2c). However, after PSM adjustment, there was no significant difference in PFS between the two groups. The 12-, 24-, and 36-month PFS rates in the lenvatinib group were 59.4%, 41.4%, and 29.5%, respectively, while those in the bevacizumab group were 50.5%, 35.3%, and 23.5%, respectively (HR, 0.808; 95% CI, 0.535–1.222; p = 0.33; Fig. 2d).

After inverse probability of treatment weighting, the OS in the lenvatinib group remained significantly longer, with OS rates at 12, 24, and 36 months of 86.5%, 63.3%, and 54.5%, respectively, compared to 68.1%, 57.2%, and 51.9% in the bevacizumab group (HR, 0.549; 95% CI, 0.331–0.908; p = 0.01; Fig. 2e). The 12-, 24-, and 36-month PFS rates in the lenvatinib group were 59.2%, 40.4%, and 34.3%, respectively, while the PFS rates in the bevacizumab group were 52.3%, 38.3%, and 29.7% (Fig. 2f).

Analysis of Background Factors

The analysis of OS and PFS among different subgroups revealed several noteworthy trends:

Conversion surgery (Fig. 3a): OS was significantly longer in the lenvatinib group compared to the bevacizumab group for patients undergoing conversion surgery (HR, 0.469; 95% CI, 0.256–0.859; p = 0.012). This highlights the potential survival benefit of conversion surgery in the lenvatinib cohort.

Fig. 3.

a–o Kaplan-Meier curves of OS and PFS for patients in the different background groups.

Fig. 3.

a–o Kaplan-Meier curves of OS and PFS for patients in the different background groups.

Close modal

Body mass index (BMI) (Fig. 3b–e): In the lenvatinib group, obese patients demonstrated the best OS outcomes, while underweight patients had the poorest prognosis (p = 0.02, Fig. 3b). In contrast, no significant OS differences were observed among BMI subgroups in the bevacizumab group (p = 0.26, Fig. 3c). Regarding PFS, a trend toward improved outcomes was observed for obese patients in the lenvatinib group, though not statistically significant (p = 0.082, Fig. 3d). No PFS differences were observed among BMI subgroups in the bevacizumab group (p = 0.32, Fig. 3e).

Hypertension (Fig. 3f, g): For OS, patients without hypertension in the bevacizumab group exhibited a rapid mid-term decline, stabilizing later, with no significant differences observed overall (p = 0.09, Fig. 3f). For PFS, no significant differences were detected between hypertensive and nonhypertensive subgroups in either treatment group (p = 0.77, Fig. 3g).

Diabetes mellitus (Fig. 3h, i): In both the lenvatinib and bevacizumab groups, diabetic patients showed a trend toward better OS compared to nondiabetic patients, though this was not statistically significant (p = 0.062, Fig. 3h). Similarly, a trend toward improved PFS was noted in diabetic patients in both groups (p = 0.19, Fig. 3i).

Hyperlipidemia (Fig. 3j, k): No significant differences were observed in OS (p = 0.093, Fig. 3j) or PFS (p = 0.64, Fig. 3k) between patients with and without hyperlipidemia in either treatment group.

Alcohol consumption (Fig. 3l, m): In the lenvatinib group, patients without a history of alcohol consumption exhibited significantly better OS compared to other subgroups (p = 0.0066, Fig. 3l). However, no significant differences in PFS were observed based on alcohol consumption in either group (p = 0.36, Fig. 3m).

Number of interventional treatments (Fig. 3n, o): In the lenvatinib group, patients undergoing 4–6 interventional treatments tended to achieve higher and more stable OS rates, though not statistically significant (p = 0.064, Fig. 3n). In the bevacizumab group, no significant OS differences were observed across subgroups with varying numbers of interventional treatments (p = 0.79, Fig. 3o).

These findings provide insights into how patient background factors and treatment modalities influence survival outcomes, underscoring the importance of personalized approaches in managing HCC.

Primary Outcomes in Subgroups

In subgroup analysis, patients aged ≤65 years significantly benefited from lenvatinib-based triple therapy compared with bevacizumab-based triple therapy (OS HR, 0.48; 95% CI, 0.27–0.85; p = 0.012), whereas no difference was observed in patients older than 65 years (OS HR, 1.02; 95% CI, 0.36–2.93; p = 0.968; Fig. 4). Male patients significantly benefited from lenvatinib combination therapy (HR, 0.48; 95% CI, 0.28–0.82; p = 0.008), whereas no difference was observed in female patients (HR, 1.70; 95% CI, 0.31–9.34). Patients without HBV infection also significantly benefited from lenvatinib combination therapy (HR, 0.18; 95% CI, 0.04–0.80; p = 0.025), whereas no difference was observed in patients with previous HBV infection (HR, 0.63; 95% CI, 0.36–1.11; p = 0.108). Interestingly, patients with extrahepatic metastasis (BCLC stage C, CNLC stage IIIb) (HR, 0.51; 95% CI, 0.30–0.89; p = 0.017) or vascular tumor thrombus (BCLC stage C, CNLC stage IIIa) (HR, 0.39; 95% CI, 0.17–0.87; p = 0.021) showed a significant survival benefit with lenvatinib combination therapy. In addition, patients with ALT <40 U/L (HR, 0.45; 95% CI, 0.23–0.88; p = 0.019), platelets ≥100 × 109/L (HR, 0.51; 95% CI, 0.29–0.91; p = 0.022), and log 10 AFP ≥1.40 ng/mL (HR, 0.55; 95% CI, 0.31–0.98; p = 0.044) also showed a significant survival benefit with lenvatinib combination therapy (Fig. 4).

Fig. 4.

Forest plot of subgroup analysis for OS. HBV, hepatitis B virus; AFP, alpha-fetoprotein; BCLC, Barcelona Clinic Liver Cancer; CNLC, China liver cancer; HR, hazard ratio.

Fig. 4.

Forest plot of subgroup analysis for OS. HBV, hepatitis B virus; AFP, alpha-fetoprotein; BCLC, Barcelona Clinic Liver Cancer; CNLC, China liver cancer; HR, hazard ratio.

Close modal

Tumor Responses

The swimmer plot for OS and the last dose of treatment are shown in Figure 5. As of the last follow-up, 62 patients (40.8%) in the lenvatinib group and 34 patients (41.5%) in the bevacizumab group were still receiving treatment. Among those who discontinued the triple therapy, in the lenvatinib group, 63.3% of patients switched to regorafenib, while 26.7% switched to the T+A regimen, and the remaining patients were lost to follow-up. In the bevacizumab group, 54.2% of patients switched to regorafenib, 25.0% switched to donafenib, 14.6% switched to lenvatinib, and the remaining patients were lost to follow-up. The reasons for discontinuation included successful conversion to resectable disease, adverse events, drug resistance or tumor progression, lack of response, or loss to follow-up.

Fig. 5.

Swimmer plot of overall survival and treatment duration in lenvatinib and bevacizumab groups.

Fig. 5.

Swimmer plot of overall survival and treatment duration in lenvatinib and bevacizumab groups.

Close modal

Regarding the best overall response, 6 patients (4.8%) in the lenvatinib group achieved a complete response (CR), and 47 patients (37.9%) achieved a partial response. In the bevacizumab group, 1 patient (1.9%) achieved CR, and 18 patients (33.3%) achieved partial response (Figure 6). The ORRs were 42.7% in the lenvatinib group and 35.2% in the bevacizumab group. Additionally, the disease control rates (DCRs) were 88.7% for lenvatinib and 96.3% for bevacizumab.

Fig. 6.

Tumor responses rates based on mRECIST in the lenvatinib group and the bevacizumab group before PSM and after PSM.

Fig. 6.

Tumor responses rates based on mRECIST in the lenvatinib group and the bevacizumab group before PSM and after PSM.

Close modal

Conversion Rate

The conversion rate to surgery was evaluated in both groups during the study. The rates were 31.0% in the lenvatinib group and 35.4% in the bevacizumab group (p = 0.406), with no statistically significant difference observed. These results suggest that both therapies have a similar potential to enable conversion to surgery, despite differences in other efficacy metrics such as OS and PFS.

Safety

The safety profile of lenvatinib and bevacizumab was assessed through the incidence of TRAEs (Table 3). The overall incidence of any grade TRAEs was 86% in the lenvatinib group and 91% in the bevacizumab group, with no statistically significant difference observed between the two (p = 0.22). For grade 3 and 4 TRAEs, the incidence was 34% in the lenvatinib group and 38% in the bevacizumab group, also showing no significant difference (p = 0.67).

Table. 3.

Incidence of adverse events

CharacteristicBevacizumab, n (%) (N = 82)Lenvatinib, n (%) (N = 152)p value1
Any gradeGrades 3 and 4Any gradeGrades 3 and 4Any gradeGrades 3 and 4
Overall incidence 75 (91) 31 (38) 130 (86) 52 (34) 0.22 0.67 
Elevated GGT 35 (43) 8 (9.8) 79 (52) 14 (9.2) 0.40 1.00 
Elevated AST 45 (55) 6 (7.3) 79 (52) 10 (6.6) 0.68 0.79 
Decreased HGB 36 (44) 4 (4.9) 69 (45) 10 (6.6) 0.89 0.78 
Elevated ALP 27 (33) 2 (2.4) 55 (36) 4 (2.6) 0.67 1.00 
Decreased PLT 32 (39) 10 (12) 45 (30) 9 (5.9) 0.15 0.13 
Elevated ALT 24 (29) 7 (8.5) 41 (27) 5 (3.3) 0.76 0.12 
Decreased WBC 13 (16) 1 (1.2) 26 (17) 2 (1.3) 0.86 1.00 
Proteinuria 8 (9.8) 3 (3.7) 24 (16) 9 (5.9) 0.24 0.55 
GI bleeding 19 (23) 12 (15) 21 (14) 12 (7.9) 0.10 0.12 
CharacteristicBevacizumab, n (%) (N = 82)Lenvatinib, n (%) (N = 152)p value1
Any gradeGrades 3 and 4Any gradeGrades 3 and 4Any gradeGrades 3 and 4
Overall incidence 75 (91) 31 (38) 130 (86) 52 (34) 0.22 0.67 
Elevated GGT 35 (43) 8 (9.8) 79 (52) 14 (9.2) 0.40 1.00 
Elevated AST 45 (55) 6 (7.3) 79 (52) 10 (6.6) 0.68 0.79 
Decreased HGB 36 (44) 4 (4.9) 69 (45) 10 (6.6) 0.89 0.78 
Elevated ALP 27 (33) 2 (2.4) 55 (36) 4 (2.6) 0.67 1.00 
Decreased PLT 32 (39) 10 (12) 45 (30) 9 (5.9) 0.15 0.13 
Elevated ALT 24 (29) 7 (8.5) 41 (27) 5 (3.3) 0.76 0.12 
Decreased WBC 13 (16) 1 (1.2) 26 (17) 2 (1.3) 0.86 1.00 
Proteinuria 8 (9.8) 3 (3.7) 24 (16) 9 (5.9) 0.24 0.55 
GI bleeding 19 (23) 12 (15) 21 (14) 12 (7.9) 0.10 0.12 

GGT, gamma-glutamyl transferase.

1Fisher’s exact test.

The most frequent TRAEs included elevated gamma-glutamyl transferase, elevated AST, and decreased hemoglobin (HGB). Elevated gamma-glutamyl transferase was observed in 52% of the lenvatinib group and 43% of the bevacizumab group (p = 0.40), with grade 3 and 4 incidences of 9.2% and 9.8%, respectively (p = 1.00). Elevated AST was reported in 52% of the lenvatinib group and 55% of the bevacizumab group (p = 0.68), with grade 3 and 4 incidences of 6.6% and 7.3%, respectively (p = 0.79). Decreased HGB was noted in 45% of the lenvatinib group and 44% of the bevacizumab group (p = 0.89), with grade 3 and 4 incidences of 6.6% and 4.9%, respectively (p = 0.78).

Proteinuria was observed in 15.8% of the lenvatinib group and 9.8% of the bevacizumab group (p = 0.24), with grade 3 and 4 incidences of 5.9% and 3.7%, respectively (p = 0.55). Gastrointestinal (GI) bleeding was reported in 14% of the lenvatinib group and 23% of the bevacizumab group (p = 0.10), with grade 3 and 4 incidences of 7.9% and 15%, respectively (p = 0.12). Notably, the bevacizumab group demonstrated a higher incidence of decreased platelet counts (PLT) compared to the lenvatinib group (39.0% vs. 29.6% for any grade, and 12.2% vs. 5.9% for grades 3 and 4), which may contribute to the higher incidence of GI bleeding in this cohort. This observation suggests a potential link between thrombocytopenia and GI hemorrhagic events (Fig. 7).

Fig. 7.

Pyramid plot of adverse events.

Fig. 7.

Pyramid plot of adverse events.

Close modal

The treatment of liver cancer is progressing rapidly. Although dual therapy is the mainstream, triple therapy combining systemic therapy with locoregional therapy has gained attention due to its high ORR rates [3, 4]. The EMERALD-1 study [29] found that combining interventional therapy with bevacizumab and durvalumab shows a statistically significant improvement in PFS (median PFS 15.0 vs. 8.2 months; HR = 0.77; CI: 0.61–0.98; p = 0.032). The LEAP-012 study [19] demonstrated a significant improvement in PFS with TACE plus lenvatinib and pembrolizumab compared to TACE plus placebo (median PFS 14.6 vs. 10.0 months; HR = 0.66; 95% CI: 0.51–0.84; p = 0.0002). Additionally, several phase III studies combining TACE/HAIC with targeted immunotherapy have been conducted, including TALEN-TACE (NCT04712643), TACE-3 (NCT04268888), and CheckMate-74W (NCT04340193) [30]. This includes previous retrospective studies by CLEAP [23], which also demonstrated the advantages of triple therapy over systemic therapy without TACE/HAIC. Targeted therapy is mainly divided into two categories: large-molecule bevacizumab and small-molecule TKIs. Bevacizumab has shown the best outcomes when combined with atezolizumab [13, 28] or sintilimab [31]. Small-molecule multitarget drugs include lenvatinib, sorafenib, donafenib, apatinib, and second-line regorafenib and cabozantinib. The large-molecule ramucirumab, which targets VEGFR2, is rarely used. Among small molecule multitarget drugs, lenvatinib has the highest ORR of 18.8% [32]. Given the positive results from the LEAP-012 study [19], which demonstrated superior efficacy of TACE combined with lenvatinib and pembrolizumab compared to TACE alone in patients with HCC, the combination of lenvatinib and immune checkpoint inhibitors (ICI) continues to be favored by many clinical experts. This supports the potential of lenvatinib plus ICB to enhance the therapeutic outcomes of TACE, further validated by encouraging results from several phase II studies [3, 4, 17, 33, 34]. Currently, a phase III RCT combining lenvatinib and toripalimab is also in progress. Therefore, comparing the large-molecule bevacizumab with the small-molecule lenvatinib in the context of triple therapy is both important and interesting.

This study is the first real-world, multicenter comparison of the efficacy and safety of triple therapy between lenvatinib and bevacizumab, both combined with PD-1 antibodies and interventional therapy, as first-line treatment in patients with uHCC. It reveals the advantages of lenvatinib combined with immunotherapy and interventional therapy in improving survival, particularly in patients who are not infected with HBV and have HCC with macrovascular invasion, normal alanine transaminase levels, normal platelet count, and low AFP levels. These findings are consistent with previous clinical observations and theories [6, 35‒43]. Notably, our study found that patients without HBV infection and those with macrovascular invasion tend to benefit more from lenvatinib-based triple therapy, resonating with the findings of Camma et al. [44]. Additionally, our study highlights the survival differences in patients with extrahepatic metastasis and tumor thrombus, which aligns with discussions on treatment strategies for metastatic HCC in previous studies [45]. Our findings suggest that multiple factors need to be considered when selecting appropriate patients for treatment, including but not limited to viral status, tumor characteristics, liver inflammation, and AFP levels [8, 27, 41, 42, 46]. This underscores the importance of personalized treatment in uHCC.

The current study revealed a higher ORR in lenvatinib-based triple therapy compared to bevacizumab-based triple therapy. Several factors may contribute to this difference. First, as reported in the literature, the ORR is 18.8% for lenvatinib monotherapy [11], 21.0% for sintilimab combined with bevacizumab biosimilar IBI305 [31], and 30% for T+A [47], while the phase II Keynote-524 study [15] and others have shown a 36% ORR. The higher proportion of sintilimab plus IBI305 in the bevacizumab group may explain part of this outcome. Second, lenvatinib targets VEGFRs, PDGFRs, Kit, FGFRs, and other pathways, and its multitarget mechanism of action makes it more effective in treatment [7, 26, 35, 42, 48], potentially leading to a higher ORR and faster clinical response. Third, lenvatinib has also been reported to modulate antitumoral immunity [49, 50], which may contribute to its higher ORR.

Our study found that OS in the lenvatinib group was significantly better than in the bevacizumab group, which may be attributed to several factors, including a higher depth of response in the lenvatinib group. Literature has indicated a correlation between depth of response and OS [51]. This may also explain the higher ORR in the lenvatinib group. Second, bevacizumab treatment may be associated with adverse events such as bleeding [13], which could negatively impact survival. Also, the similar PFS but higher percentage of death and shorter OS in the bevacizumab may suggest a higher portion of nonprogressive death. However, this hypothesis requires further evidence. Third, for patients undergoing conversion surgery, the withdrawal period for lenvatinib is relatively short (1–2 weeks), whereas bevacizumab requires 4–6 weeks. The shorter discontinuation time may have reduced the risk of tumor progression and spreading due to discontinuation in the lenvatinib group, thereby positively impacting OS. Additionally, available data provide further clues. In the LEAP-002 study [14], the median OS of the lenvatinib combination was 21.2 months, and the OS of the Chinese subgroup was 32.3 months. The OS data from the T+A regimen, which was 24.0 months in the Asian subgroup [13], also support the long-term survival advantage of lenvatinib in triple therapy. Additionally, our study observed that patients with diabetes mellitus and obesity in the lenvatinib group exhibited better OS. This finding aligns with emerging evidence that metabolic factors may positively influence treatment outcomes in advanced HCC. Specifically, the improved OS in these subgroups could be related to altered drug metabolism and enhanced efficacy of lenvatinib, potentially mediated through its impact on the tumor microenvironment. This hypothesis is supported by recent studies, demonstrating that metabolic alterations, such as glucose restriction or inhibition of glucose transporters, can significantly enhance the anticancer effects of targeted therapies in HCC models [52]. Therefore, our findings underscore the need to further explore the role of metabolic status in optimizing personalized treatment strategies for patients with unresectable HCC.

Despite the ORR and OS advantage in the lenvatinib group, the DCR was relatively lower in the lenvatinib group, which may account for the nonsignificant difference in PFS between the two treatment groups. The median PFS for the LEAP002 study [14] was 8.2 months, compared to 4.6 months in the ORIENT-32 study [28] and 6.8 months in the IMBRAVE-150 study [13]. Since PFS was similar in both groups, these data may explain why lenvatinib and bevacizumab performed similarly in triple therapy. This study also revealed that the PD rate was higher for lenvatinib than for bevacizumab, indicating that bevacizumab may be more stable in controlling disease progression [53].

The difference between PFS and OS can be attributed to several factors. First, we speculate that the duration of response for bevacizumab may be relatively short [47], meaning that lenvatinib may produce a more rapid tumor shrinkage effect once it takes effect [54], thus acceleratting conversion to surgical resection, with a relatively short period of drug discontinuation. Additionally, the incidence of bevacizumab-related adverse reactions, such as GI bleeding [55], is relatively higher, which may lead to a scenario where bevacizumab achieves a better DCR but has a higher proportion of deaths.

However, retrospective studies are at risk of selection bias, and our results, although statistically significant, need to be validated by more rigorous prospective studies and randomized controlled trials [56, 57]. Moreover, our study did not directly compare lenvatinib and bevacizumab triple therapy in terms of drug resistance, toxicity, and quality of life, which are important considerations in practical clinical decisions [58, 59].

In conclusion, our study provides real-world evidence for the OS benefits of lenvatinib-based triple therapy over bevacizumab-based therapy in a subset of uHCC patients. Moreover, individual patient characteristics, including viral status and biomarkers, should be considered when selecting the most appropriate treatment. Further research is essential to refine patient stratification and explore biomarkers for predicting treatment response [60, 61].

This study involving human participants was reviewed and approved by Zhongshan Hospital, Fudan University (Approval No. B2022-195R). The patient/participants provided their written informed consent to participate in this study.

The authors have no conflicts of interest to declare.

This research was supported by funds as follows: (1) Tianjin Medical University Cancer Hospital, Precision Treatment Technology Construction Project for Cancer Surgery, ZLWKJZZL14; (2) Zhongguancun Precision Medicine Foundation, Medical and Health Public Welfare – Cancer Medical Research Special Project, ZGC-YXKY-ZL004 & ZGC-YJWC-ZLYXKY-9. (3) Tianjin Key Medical Discipline(Specialty) Construction Project (TJYXZDXK-011A).

Z.P., J.C., and W.Z. designed experiments and drafted the manuscript. D.L., L.F., X.Z., X.Z., Y.P., J.L., C.H., R.J., S.S., X.Z., H.L., X.Y., K.H., X.S., D.W., Y.Z., J.Z., B.X., S.G., T.L., S.Z., L.Z., H.Z., YZ, TW, MK, X.L., T.P., K.W., L.X., H.L., T.S., and H.S. collected the cases. W.Z. approved the final version. The China Liver Cancer Study Group Yonug Investigators (CLEAP) provided the platform for data maintenance. All authors contributed to the article and approved the submitted version.

Additional Information

Zhaolong Pan, Dongming Liu, Junbo Cao, and Xiaodong Zhu contributed equally to this work.

The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding authors.

1.
Singh
SP
,
Madke
T
,
Chand
P
.
Global epidemiology of hepatocellular carcinoma
.
J Clin Exp Hepatol
.
2025
;
15
(
2
):
102446
.
2.
Zhou
J
,
Sun
H
,
Wang
Z
,
Cong
W
,
Zeng
M
,
Zhou
W
, et al
.
Guidelines for the diagnosis and treatment of primary liver cancer (2022 edition)
.
Liver Cancer
.
2023
;
12
(
5
):
405
44
.
3.
Cai
H
,
Chen
S
,
Wu
Z
,
Wang
F
,
Tang
S
,
Chen
L
, et al
.
Atezolizumab plus bevacizumab combined with transarterial embolization plus hepatic arterial infusion chemotherapy for unresectable hepatocellular carcinoma with a diameter >8 cm: a retrospective study
.
J Hepatocell Carcinoma
.
2024
;
11
:
399
409
.
4.
Wang
L
,
Lin
L
,
Zhou
W
.
Efficacy and safety of transarterial chemoembolization combined with lenvatinib and PD-1 inhibitor in the treatment of advanced hepatocellular carcinoma: a meta-analysis
.
Pharmacol Ther
.
2024
;
257
:
108634
.
5.
Duan
Y
,
Zhang
H
,
Tan
T
,
Ye
W
,
Yin
K
,
Yu
Y
, et al
.
The immune response of hepatocellular carcinoma after locoregional and systemic therapies: the available combination option for immunotherapy
.
Biosci Trends
.
2024
;
17
(
6
):
427
44
.
6.
Yang
C
,
Zhang
H
,
Zhang
L
,
Zhu
AX
,
Bernards
R
,
Qin
W
, et al
.
Evolving therapeutic landscape of advanced hepatocellular carcinoma
.
Nat Rev Gastroenterol Hepatol
.
2023
;
20
(
4
):
203
22
.
7.
Faivre
S
,
Rimassa
L
,
Finn
RS
.
Molecular therapies for HCC: looking outside the box
.
J Hepatol
.
2020
;
72
(
2
):
342
52
.
8.
Sonbol
MB
,
Riaz
IB
,
Naqvi
SAA
,
Almquist
DR
,
Mina
S
,
Almasri
J
, et al
.
Systemic therapy and sequencing options in advanced hepatocellular carcinoma: a systematic review and network meta-analysis
.
JAMA Oncol
.
2020
;
6
(
12
):
e204930
.
9.
Melero
I
,
Yau
T
,
Kang
Y-K
,
Kim
T-Y
,
Santoro
A
,
Sangro
B
, et al
.
Nivolumab plus ipilimumab combination therapy in patients with advanced hepatocellular carcinoma previously treated with sorafenib: 5-year results from CheckMate 040
.
Ann Oncol
.
2024
;
35
(
6
):
537
48
.
10.
Wu
J
,
Wu
J
,
Li
S
,
Luo
M
,
Zeng
Z
,
Li
Y
, et al
.
Effect of transcatheter arterial chemoembolization combined with lenvatinib plus anti-PD-1 antibodies in patients with unresectable hepatocellular carcinoma: a treatment with Chinese characteristics
.
Biosci Trends
.
2024
;
18
(
1
):
42
8
.
11.
Kudo
M
,
Finn
RS
,
Qin
S
,
Han
K-H
,
Ikeda
K
,
Piscaglia
F
, et al
.
Lenvatinib versus sorafenib in first-line treatment of patients with unresectable hepatocellular carcinoma: a randomised phase 3 non-inferiority trial
.
Lancet
.
2018
;
391
(
10126
):
1163
73
.
12.
Lee
MS
,
Ryoo
B-Y
,
Hsu
C-H
,
Numata
K
,
Stein
S
,
Verret
W
, et al
.
Atezolizumab with or without bevacizumab in unresectable hepatocellular carcinoma (GO30140): an open-label, multicentre, phase 1b study
.
Lancet Oncol
.
2020
;
21
(
6
):
808
20
.
13.
Finn
RS
,
Qin
S
,
Ikeda
M
,
Galle
PR
,
Ducreux
M
,
Kim
T-Y
, et al
.
Atezolizumab plus bevacizumab in unresectable hepatocellular carcinoma
.
N Engl J Med
.
2020
;
382
(
20
):
1894
905
.
14.
Llovet
JM
,
Kudo
M
,
Merle
P
,
Meyer
T
,
Qin
S
,
Ikeda
M
, et al
.
Lenvatinib plus pembrolizumab versus lenvatinib plus placebo for advanced hepatocellular carcinoma (LEAP-002): a randomised, double-blind, phase 3 trial
.
Lancet Oncol
.
2023
;
24
(
12
):
1399
410
.
15.
Finn
RS
,
Ikeda
M
,
Zhu
AX
,
Sung
MW
,
Baron
AD
,
Kudo
M
, et al
.
Phase ib study of lenvatinib plus pembrolizumab in patients with unresectable hepatocellular carcinoma
.
J Clin Oncol
.
2020
;
38
(
26
):
2960
70
.
16.
Huang
J-T
,
Zhong
B-Y
,
Jiang
N
,
Li
W-C
,
Zhang
S
,
Yin
Y
, et al
.
Transarterial chemoembolization combined with immune checkpoint inhibitors plus tyrosine kinase inhibitors versus immune checkpoint inhibitors plus tyrosine kinase inhibitors for advanced hepatocellular carcinoma
.
J Hepatocell Carcinoma
.
2022
;
9
:
1217
28
.
17.
Zhu
H-D
,
Li
H-L
,
Huang
M-S
,
Yang
W-Z
,
Yin
G-W
,
Zhong
B-Y
, et al
.
Transarterial chemoembolization with PD-(L)1 inhibitors plus molecular targeted therapies for hepatocellular carcinoma (CHANCE001)
.
Signal Transduct Target Ther
.
2023
;
8
(
1
):
58
.
18.
Sangro
B
,
Kudo
M
,
Erinjeri
JP
,
Qin
S
,
Ren
Z
,
Chan
SL
, et al
.
Durvalumab with or without bevacizumab with transarterial chemoembolisation in hepatocellular carcinoma (EMERALD-1): a multiregional, randomised, double-blind, placebo-controlled, phase 3 study
.
Lancet
.
2025
;
405
(
10474
):
216
32
.
19.
Kudo
M
,
Ren
Z
,
Guo
Y
,
Han
G
,
Lin
H
,
Zheng
J
, et al
.
Transarterial chemoembolisation combined with lenvatinib plus pembrolizumab versus dual placebo for unresectable, non-metastatic hepatocellular carcinoma (LEAP-012): a multicentre, randomised, double-blind, phase 3 study
.
Lancet
.
2025
;
405
(
10474
):
203
15
.
20.
Liu
J
,
Zhu
X
,
Pan
Y
,
Zhong
J
,
Jin
R
,
Zheng
X
, et al
.
Prognoses of patients treated with surgical therapy versus continuation of local-plus-systemic therapy following successful down-staging of intermediate-advanced hepatocellular carcinoma: a multicenter real-world study
.
Oncologist
.
2024
;
29
(
4
):
e487
97
.
21.
Huang
C
,
Zhu
X-D
,
Shen
Y-H
,
Xu
B
,
Wu
D
,
Ji
Y
, et al
.
Radiographic and α-fetoprotein response predict pathologic complete response to immunotherapy plus a TKI in hepatocellular carcinoma: a multicenter study
.
BMC Cancer
.
2023
;
23
(
1
):
416
.
22.
Zhang
X
,
Zhu
X
,
Zhong
J
,
Zhao
Y
,
Zhang
X
,
Zhang
W
, et al
.
Surgical treatment improves overall survival of hepatocellular carcinoma with extrahepatic metastases after conversion therapy: a multicenter retrospective study
.
Sci Rep
.
2024
;
14
(
1
):
9745
.
23.
Pan
Y
,
Zhu
X
,
Liu
J
,
Zhong
J
,
Zhang
W
,
Shen
S
, et al
.
Systemic therapy with or without transcatheter intra-arterial therapies for unresectable hepatocellular carcinoma: a real-world, multi-center study
.
Front Immunol
.
2023
;
14
:
1138355
.
24.
Yang
C
,
Pan
Y
,
Ye
F
,
Zhu
X
,
Xue
J
,
Li
X
, et al
.
The pattern of tumor progression on first-line immune checkpoint inhibitor-based systemic therapy for Chinese advanced hepatocellular carcinoma -CLEAP 004 study
.
Front Immunol
.
2024
;
15
:
1310239
.
25.
Xu
B
,
Dong
S
,
Bai
X
,
Song
T
,
Zhang
B
,
Zhou
L
, et al
.
Tumor radiomic features on pretreatment MRI to predict response to lenvatinib plus an anti-PD-1 antibody in advanced hepatocellular carcinoma: a multicenter study
.
Liver Cancer
.
2023
;
12
(
3
):
262
76
.
26.
Bruix
J
,
Reig
M
,
Sherman
M
.
Evidence-based diagnosis, staging, and treatment of patients with hepatocellular carcinoma
.
Gastroenterology
.
2016
;
150
(
4
):
835
53
.
27.
Zhou
J
,
Sun
H
,
Wang
Z
,
Cong
W
,
Wang
J
,
Zeng
M
, et al
.
Guidelines for the diagnosis and treatment of hepatocellular carcinoma (2019 edition)
.
Liver Cancer
.
2020
;
9
(
6
):
682
720
.
28.
Ren
Z
,
Xu
J
,
Bai
Y
,
Xu
A
,
Cang
S
,
Du
C
, et al
.
Sintilimab plus a bevacizumab biosimilar (IBI305) versus sorafenib in unresectable hepatocellular carcinoma (ORIENT-32): a randomised, open-label, phase 2-3 study
.
Lancet Oncol
.
2021
;
22
(
7
):
977
90
.
29.
Lencioni
R
,
Kudo
M
,
Erinjeri
J
,
Qin
S
,
Ren
Z
,
Chan
S
, et al
.
EMERALD-1: a phase 3, randomized, placebo-controlled study of transarterial chemoembolization combined with durvalumab with or without bevacizumab in participants with unresectable hepatocellular carcinoma eligible for embolization
.
J Clin Oncol
.
2024
;
42
(
3_Suppl
):
LBA432
.
30.
Moriguchi
M
,
Kataoka
S
,
Itoh
Y
.
Evolution of systemic treatment for hepatocellular carcinoma: changing treatment strategies and concepts
.
Cancers
.
2024
;
16
(
13
):
2387
.
31.
Ren
Z
,
Xu
J
,
Bai
Y
,
Xu
A
,
Cang
S
,
Du
C
, et al
.
ORIENT-32: updated characterization of response to sintilimab plus bevacizumab biosimilar (IBI305) vs sorafenib for unresectable hepatocellular carcinoma
.
J Clin Oncol
.
2023
;
41
(
4_Suppl
):
570
.
32.
Dong
Y
,
Wong
JSL
,
Sugimura
R
,
Lam
K-O
,
Li
B
,
Kwok
GGW
, et al
.
Recent advances and future prospects in immune checkpoint (ICI)-Based combination therapy for advanced HCC
.
Cancers
.
2021
;
13
(
8
):
1949
.
33.
Gao
B
,
Yang
F
,
Zheng
D
,
Hu
S
,
Liu
J
,
Liu
H
, et al
.
Transarterial chemoembolization combined with tyrosine kinase inhibitors plus immune checkpoint inhibitors for advanced hepatocellular carcinoma: a propensity score matching analysis
.
J Hepatocell Carcinoma
.
2023
;
10
:
2265
76
.
34.
Zhang
W
,
Tong
S
,
Hu
B
,
Wan
T
,
Tang
H
,
Zhao
F
, et al
.
Lenvatinib plus anti-PD-1 antibodies as conversion therapy for patients with unresectable intermediate-advanced hepatocellular carcinoma: a single-arm, phase II trial
.
J Immunother Cancer
.
2023
;
11
(
9
):
e007366
.
35.
Llovet
JM
,
Ricci
S
,
Mazzaferro
V
,
Hilgard
P
,
Gane
E
,
Blanc
J-F
, et al
.
Sorafenib in advanced hepatocellular carcinoma
.
N Engl J Med
.
2008
;
359
(
4
):
378
90
.
36.
Galle
PR
,
Finn
RS
,
Qin
S
,
Ikeda
M
,
Zhu
AX
,
Kim
T-Y
, et al
.
Patient-reported outcomes with atezolizumab plus bevacizumab versus sorafenib in patients with unresectable hepatocellular carcinoma (IMbrave150): an open-label, randomised, phase 3 trial
.
Lancet Oncol
.
2021
;
22
(
7
):
991
1001
.
37.
Yi
C
,
Chen
L
,
Lin
Z
,
Liu
L
,
Shao
W
,
Zhang
R
, et al
.
Lenvatinib targets FGF receptor 4 to enhance antitumor immune response of anti-programmed cell death-1 in HCC
.
Hepatology
.
2021
;
74
(
5
):
2544
60
.
38.
Yang
J
,
Guo
Z
,
Song
M
,
Pan
Q
,
Zhao
J
,
Huang
Y
, et al
.
Lenvatinib improves anti-PD-1 therapeutic efficacy by promoting vascular normalization via the NRP-1-PDGFRβ complex in hepatocellular carcinoma
.
Front Immunol
.
2023
;
14
:
1212577
.
39.
Yau
T
,
Park
J-W
,
Finn
RS
,
Cheng
A-L
,
Mathurin
P
,
Edeline
J
, et al
.
Nivolumab versus sorafenib in advanced hepatocellular carcinoma (CheckMate 459): a randomised, multicentre, open-label, phase 3 trial
.
Lancet Oncol
.
2022
;
23
(
1
):
77
90
.
40.
Girardi
DM
,
Sousa
LP
,
Miranda
TA
,
Haum
FNC
,
Pereira
GCB
,
Pereira
AAL
.
Systemic therapy for advanced hepatocellular carcinoma: current stand and perspectives
.
Cancers
.
2023
;
15
(
6
):
1680
.
41.
Campani
C
,
Nault
J-C
.
Systemic treatment of hepatocellular carcinoma: the times they are a-changin
.
Hepatobiliary Surg Nutr
.
2021
;
10
(
6
):
893
5
.
42.
Fan
Y
,
Xue
H
,
Zheng
H
.
Systemic therapy for hepatocellular carcinoma: current updates and outlook
.
J Hepatocell Carcinoma
.
2022
;
9
:
233
63
.
43.
Dal Bo
M
,
De Mattia
E
,
Baboci
L
,
Mezzalira
S
,
Cecchin
E
,
Assaraf
YG
, et al
.
New insights into the pharmacological, immunological, and CAR-T-cell approaches in the treatment of hepatocellular carcinoma
.
Drug Resist Updat
.
2020
;
51
:
100702
.
44.
Cerreto
M
,
Cardone
F
,
Cerrito
L
,
Stella
L
,
Santopaolo
F
,
Pallozzi
M
, et al
.
The new era of systemic treatment for hepatocellular carcinoma: from the first line to the optimal sequence
.
Curr Oncol
.
2023
;
30
(
10
):
8774
92
.
45.
Huang
L
,
Kang
D
,
Zhao
C
,
Liu
X
.
Correlation between surrogate endpoints and overall survival in unresectable hepatocellular carcinoma patients treated with immune checkpoint inhibitors: a systematic review and meta-analysis
.
Sci Rep
.
2024
;
14
(
1
):
4327
.
46.
Zhang
X-D
,
Wang
Y
,
Ye
L-H
.
Hepatitis B virus X protein accelerates the development of hepatoma
.
Cancer Biol Med
.
2014
;
11
(
3
):
182
90
.
47.
Cheng
A-L
,
Qin
S
,
Ikeda
M
,
Galle
PR
,
Ducreux
M
,
Kim
T-Y
, et al
.
Updated efficacy and safety data from IMbrave150: atezolizumab plus bevacizumab vs. sorafenib for unresectable hepatocellular carcinoma
.
J Hepatol
.
2022
;
76
(
4
):
862
73
.
48.
Konyn
P
,
Ahmed
A
,
Kim
D
.
Current epidemiology in hepatocellular carcinoma
.
Expert Rev Gastroenterol Hepatol
.
2021
;
15
(
11
):
1295
307
.
49.
Yamauchi
M
,
Ono
A
,
Amioka
K
,
Fujii
Y
,
Nakahara
H
,
Teraoka
Y
, et al
.
Lenvatinib activates anti-tumor immunity by suppressing immunoinhibitory infiltrates in the tumor microenvironment of advanced hepatocellular carcinoma
.
Commun Med
.
2023
;
3
(
1
):
152
.
50.
Buttell
A
,
Qiu
W
.
The action and resistance mechanisms of Lenvatinib in liver cancer
.
Mol Carcinog
.
2023
;
62
(
12
):
1918
34
.
51.
Grünwald
V
,
Powles
T
,
Kopyltsov
E
,
Kozlov
V
,
Alonso-Gordoa
T
,
Eto
M
, et al
.
Survival by depth of response and efficacy by international metastatic renal cell carcinoma database consortium subgroup with lenvatinib plus pembrolizumab versus sunitinib in advanced renal cell carcinoma: analysis of the phase 3 randomized CLEAR study
.
Eur Urol Oncol
.
2023
;
6
(
4
):
437
46
.
52.
Zhou
J
,
Feng
J
,
Wu
Y
,
Dai
H-Q
,
Zhu
G-Z
,
Chen
P-H
, et al
.
Simultaneous treatment with sorafenib and glucose restriction inhibits hepatocellular carcinoma in vitro and in vivo by impairing SIAH1-mediated mitophagy
.
Exp Mol Med
.
2022
;
54
(
11
):
2007
21
.
53.
Liu
J
,
Yang
L
,
Wei
S
,
Li
J
,
Yi
P
.
Efficacy and safety of atezolizumab plus bevacizumab versus lenvatinib for unresectable hepatocellular carcinoma: a systematic review and meta-analysis
.
J Cancer Res Clin Oncol
.
2023
;
149
(
17
):
16191
201
.
54.
Takahashi
A
,
Moriguchi
M
,
Seko
Y
,
Shima
T
,
Mitsumoto
Y
,
Takashima
H
, et al
.
Early tumor shrinkage as a predictive factor for outcomes in hepatocellular carcinoma patients treated with lenvatinib: a multicenter analysis
.
Cancers
.
2020
;
12
(
3
):
754
.
55.
Ha
Y
,
Kim
JH
,
Cheon
J
,
Jeon
GS
,
Kim
C
,
Chon
HJ
.
Risk of variceal bleeding in patients with advanced hepatocellular carcinoma receiving atezolizumab/bevacizumab
.
Clin Gastroenterol Hepatol
.
2023
;
21
(
9
):
2421
3.e2
.
56.
Motzer
RJ
,
Rini
BI
,
McDermott
DF
,
Arén Frontera
O
,
Hammers
HJ
,
Carducci
MA
, et al
.
Nivolumab plus ipilimumab versus sunitinib in first-line treatment for advanced renal cell carcinoma: extended follow-up of efficacy and safety results from a randomised, controlled, phase 3 trial
.
Lancet Oncol
.
2019
;
20
(
10
):
1370
85
.
57.
Cappuyns
S
,
Corbett
V
,
Yarchoan
M
,
Finn
RS
,
Llovet
JM
.
Critical appraisal of guideline recommendations on systemic therapies for advanced hepatocellular carcinoma: a review
.
JAMA Oncol
.
2024
;
10
(
3
):
395
404
.
58.
Qin
S
,
Bi
F
,
Gu
S
,
Bai
Y
,
Chen
Z
,
Wang
Z
, et al
.
Donafenib versus sorafenib in first-line treatment of unresectable or metastatic hepatocellular carcinoma: a randomized, open-label, parallel-controlled phase II-III trial
.
J Clin Oncol
.
2021
;
39
(
27
):
3002
11
.
59.
Qin
S
,
Ren
Z
,
Feng
Y-H
,
Yau
T
,
Wang
B
,
Zhao
H
, et al
.
Atezolizumab plus bevacizumab versus sorafenib in the Chinese subpopulation with unresectable hepatocellular carcinoma: phase 3 randomized, open-label IMbrave150 study
.
Liver Cancer
.
2021
;
10
(
4
):
296
308
.
60.
Dyhl-Polk
A
,
Mikkelsen
MK
,
Ladekarl
M
,
Nielsen
DL
.
Clinical trials of immune checkpoint inhibitors in hepatocellular carcinoma
.
J Clin Med
.
2021
;
10
(
12
):
2662
.
61.
Safri
F
,
Nguyen
R
,
Zerehpooshnesfchi
S
,
George
J
,
Qiao
L
.
Heterogeneity of hepatocellular carcinoma: from mechanisms to clinical implications
.
Cancer Gene Ther
.
2024
;
31
(
8
):
1105
12
.