1 Introduction

Global cancer statistics 2020 reported that primary liver cancer was the sixth most common cancer and the third leading cause of cancer-related death worldwide, with hepatocellular carcinoma (HCC) accounting for 75–85% of primary liver cancer cases [1]. The age-standardized incidence rate of HCC was 7.3 cases per 100,000 person-year, including 11.6 cases per 100,000 in males and 3.4 cases per 100,000 in females [2]. The incidence rate of HCC had been increasing rapidly in recent decades, with a concomitant increase in mortality rate due to poor prognosis [3]. The main risk factors for HCC included cirrhosis, hepatitis B virus (HBV) and hepatitis C virus (HCV) infection, smoking, drinking, and consumption of aflatoxin-infected food. At the same time, overweight, diabetes and family history also increased the risk of HCC [4, 5].

Patients with early-stage or intermediate-stage HCC could be treated with local therapies, including liver resection, liver transplantation or trans-arterial chemoembolization (TACE). However, most patients were diagnosed only when they had reached the advanced stage. Before the clinical use of immunotherapies, tyrosine kinase inhibitors (TKIs) such as sorafenib and lenvatinib have been recommended to patients who received unresectable HCC in the last few years. Unfortunately, the prognosis for unresectable HCC patients remained unsatisfactory [6, 7]. Recently, the success of the use of immune checkpoint inhibitors (ICIs) in metastatic melanoma, lung cancer, and triple-negative breast cancer has brought new hope for unresectable HCC patients [8]. ICIs was one of the fastest-growing immunotherapeutic strategies and became a hot spot for anti-tumour drug research. Studies have reported good therapeutic effects of ICIs for solid tumours. Most of them blocked the relevant receptors to help T cells overcome the immune regulation and achieve anti-tumour effects [9].

As for unresectable HCC patients, several clinical trials showed that ICIs were effective and well-tolerated, and suggested that ICIs could be used as a treatment choice or in combination with other drugs [10,11,12]. In particular, the large clinical trial IMBrave150 reported that ICIs-based therapy, compared to sorafenib, resulted in significant benefits in terms of overall survival (OS) and progression-free survival (PFS) [13].

The National Comprehensive Cancer Network (NCCN) guidelines had included atezolizumab plus bevacizumab in the first-line treatment option for unresectable HCC and recommended the clinical use of immunotherapy for eligible patients [14].In addition, the Chinese Society of Clinical Oncology (CSCO) guidelines also recommended atezolizumab plus bevacizumab as well as other immunotherapies, such as lenvatinib plus camrelizumab or nivolumab, camrelizumab plus oxaliplatin-based chemotherapy, and camrelizumab plus apatinib as first-line treatment strategies for unresectable HCC patients [15].

At present, a research gap in comparing different ICIs-based regimens for unresectable HCC patients remains. Most of the available studies on the use of ICIs-based therapies for HCC patients only included some monotherapy or focused on both first-line and non first-line treatments at the same time [16,17,18,19,20]. With the rapid development of ICIs, evaluating them as first-line treatment in terms of efficacy and safety is needed to support evidence-based clinical practice.

Thus, we used systematic review, meta-analysis and network meta-analysis to evaluate the results of all clinical trials using ICIs-based treatment as first-line therapy for unresectable HCC patients. We hope that the results of this research could provide evidence for informing physicians’ choice of treatment and keeping researchers abreast of the latest evidence development in this field.

2 Methods

This systematic review, meta-analysis and network meta-analysis complied with the PRISMA statement [21, 22]. The study protocol was also registered on the PROSPERO website, and the registration number was CRD42021288188.

2.1 Inclusion and exclusion criteria

Inclusion criteria for clinical trials were developed in accordance with PICOS principles [23]. Clinical studies which investigated ICIs-based regimens for unresectable HCC were eligible for inclusion. The involved population must be adults with unresectable HCC who never received systematic treatment before. Patients with metastatic tumors were excluded from this study. Efficacy results included PFS, OS, objective response rate (ORR) and disease control rate (DCR), and the incidence of treatment-rated adverse events (TRAEs). The specific eligibility criteria for inclusion were shown in Table 1.

Table 1 Study inclusion and exclusion criteria

2.2 Literature search

Four databases (PubMed, Web of Science, Embase, and clinicaltrials.gov) were searched for eligible studies. The search time was from the inception of the database to 11st August 2022. The main search terms were “hepatocellular carcinoma”, “immune checkpoint inhibitor”, “PD-1”, “PD-L1”, “CTLA-4”, and their synonyms. Abstracts of recent major oncology conferences and included studies' references are also screened. The specific search strategy is shown in Supplementary materials.

2.3 Literature selection

The searched literature was screened using the literature management software Endnote X9. After removing duplicates, two researchers completed screening independently, and the selection included two rounds. In the first round, we reviewed titles and abstracts, and excluded irrelevant literature or literature that did not meet the inclusion criteria. We read the remaining literature in full text in the second round to complete the selection. In case of disagreement, a third and more senior researcher was to make a final judgement.

2.4 Data extraction

Two researchers extracted data separately, with a third researcher making the final judgement in case of any disagreement. The extracted information included the following: (1) basic information of the clinical trial, including authorship, publication date and clinical trial registration number; (2) study design of the clinical trial, including sample size, location, allocation, intervention model, masking, primary purpose, and follow-up time; (3) basic characteristics of included patients, including gender ratio, median age, and basic liver condition; (4) treatments of intervention group including drug and doses; (5) treatment of control group including drug and doses if available; and (6) outcomes of this study, including PFS, OS, ORR, DCR, TRAEs.

2.5 Quality assessment

Eligible studies were assessed by two authors using Cochrane risk of bias (ROB) tool 2.0 for randomized controlled trials (RCT) [24] and Newcastle–Ottawa scale (NOS) for single-arm clinical trials [25]. ROB tool 2.0 judged clinical trials from five perspectives. As long as one of these perspectives was identified as high risk, the study would be evaluated as high risk. NOS took a score value, with a score of 7 or more considered to be of high quality.

2.6 Statistical Analysis

The 'meta' packages in the R language was used to complete this meta-analysis with the frequency methods [26]. The pooled results of continuous variables were shown in standard deviation and 95% confidence interval (CI), and the results of dichotomous variables were shown in OR value and 95% CI [27]. I2 was used to evaluate heterogeneity. When I2 was smaller than 50% indicating a low level of heterogeneity, a fixed-effect model was used. Otherwise, a random-effect model was used whenever the heterogeneity was considered high. The significance level was set at 0.05, with P less than 0.05 indicating a significant difference. Subgroup analysis and sensitivity analysis were used to explain heterogeneity with sufficient data [28].

The 'netmeta' package in the R language was used to perform the network meta-analysis with the frequency methods. The methods for synthesizing results, determining whether to use a random-effect model or a fixed-effect model and evaluating whether there existed a significant difference, were the same as meta-analysis.

Forest plot was done using the 'forest' function in R to visualize the difference in effectiveness among interventions in both meta-analysis and network meta-analysis [27]. Software Stata 15 was used to draw a network of relationships to represent the comparative relationships that existed among interventions. The size of the vertices indicated the sample size included in each intervention, and the thickness of the lines indicated the number of studies included [29].

Using the 'netrank' function in the 'netmeta' package in R to calculate P-scores, we get established all interventions' rank. The P-score depended on the estimating result and confidence interval of the effect value, from 0 to 1, indicating a progressive increase in the efficacy of the included drugs [30].

2.7 Publication bias

Funnel plots were drawn using the 'funnel' function in 'meta' and 'netmeta' package in R to visualize publication bias and then quantified by Egger's test if there were significant asymmetries in the funnel plots [31].

3 Results

3.1 Included study

Initially, 8579 studies were identified through database searching. Of them, 2900 were excluded because of duplication. The first round of selection of 5679 remaining studies was completed by reviewing titles and abstracts. After the 1st round screening, 4520 studies were excluded because of the following reasons: publication types of 1584 studies were out of scope; study designs of 697 studies were out of scope; 12 studies were not English publications; and 2227 studies did not focus on unresectable HCC. Then, 1159 studies were reviewed by full text in the second round. Of which, 1135 studies were excluded for the following reasons: outcomes of 383 studies were not of interest; 139 studies did not report latest clinical trial data; 89 studies did not focus on unresectable HCC; 225 studies did not take ICIs as first-line treatment; and 299 studies had not completed yet. Finally, 24 studies were included in this review, including 12 full texts and 12 conference abstracts (Fig. 1).

Fig. 1
figure 1

PRISMA flowchart of study selection

3.2 Study characteristics

Study characteristics of all the included studies were listed in Table 2. A total of 10 studies were multi-center, nine were conducted in mainland China, and of the remaining trials, three were in Japan and one each in the USA and Korea.

Table 2 Characteristics of included study

3.2.1 Study design

All clinical trials were open-label and designed for treatment. Among them eight were randomized, five were non-randomized, and the allocation method of the remaining eleven clinical trials was not provided. Six trials were phase I, nine were phase II, and two were phase I/II trials, while the remaining seven clinical trials were phase III.

Thirteen of them were single-arm trials, which took patients’ pre-drug condition as comparator. Three trials compared ICIs-based therapy as first-line and second-line treatment, five selected one kind of TKIs as direct comparators, and three other clinical trials were set up with three different treatment regimens, including a direct comparison among three different ICIs-based therapies, a direct comparison among two ICIs-based therapies and TKIs, and a direct comparison among ICIs-based and two TKIs.

3.2.2 Study population

The sample size of included clinical trials ranged from 16 to 1171. Most of clinical trials were predominantly male, the gender ratio ranged from 73.8 to 91.7%. The median age of included population varied from 49 to 76 years. Eight clinical trials reported the HCC aetiology of patients, most of whom were HBV or HCV infection. The majority of patients were class A on the Child–Pugh, and the performance status score on Eastern Cooperative Oncology Group (ECOG) was of 0 or 1.

3.2.3 Study drugs

ICIs as first-line treatment in clinical trials included in this study were of three types, including anti-programmed cell death protein 1 (PD-1) antibodies, anti-programmed cell death ligand 1 (PD-L1) antibodies, and anti-cytotoxic T-lymphocyte antigen 4 (CTLA-4) antibodies. Nivolumab and atezolizumab was the most investigated drug of all treatment regimens (n = 4), followed by pembrolizumab (n = 3), toripalimab (n = 3), camrelizumab (n = 2), sintilimab (n = 2), durvalumab (n = 2), AK105 (n = 1), avelumab (n = 1), tremelimumab (n = 1), AK104 (n = 1), and ezabenlimab (n = 1).

Three studies with unresectable HCC patients were treated with monotherapy therapy of ICIs [32,33,34]. Two studies focused on the efficacy and safety of nivolumab, and the remaining one focused on pembrolizumab. A total of seven studies used combination therapy of ICIs and angiogenesis inhibitory monoclonal antibodies (AI mAbs) as the treatment of intervention group [13, 35,36,37,38,39,40]. Three studies examined atezolizumab plus bevacizumab, followed by sintilimab plus IBI305 (n = 1), toripalimab plus bevacizumab (n = 1), ezabenlimab plus BI 836,880 (n = 1), and durvalumab plus bevacizumab (n = 1). Eleven eligible studies were about combination therapy of ICIs and TKIs [41,42,43,44,45,46,47,48,49,50,51] involving 9 kinds of ICIs. The remaining three studies chose ICIs with ICIs, chemotherapy, and TKIs plus hepatic arterial infusion chemotherapy (HAIC) as first-line regimens, respectively [52,53,54].

3.3 Quality assessment

Thirteen clinical trials were single-arm studies, of which the quality was evaluated with NOS from three perspectives. As shown in Table 3, the score of these thirteen studies were all the same indicating a common concern on the comparability because the study design was a single group assignment.

Table 3 Quality assessment of single-arm clinical trials

Eleven controlled clinical trials were involved in the network meta-analysis. Firstly, we evaluated them using the Cochrane ROB Tool (2.0) for RCTs. Considering that all of them were open-label, they were marked high risk in randomization process domains. Thus, all these ten controlled trials were of high risk. The specific results are shown showed in Fig. 2.

Fig. 2
figure 2

Quality assessment of controlled clinical trials

3.4 Survival outcomes

3.4.1 Pooled results of PFS

Eighteen studies reported median PFS. Two studies compared efficacy of first-line and second-line treatment, and both of them reported a longer median PFS displayed in first-line group [32, 42].

The pooled result was 5.76 months (95% CI 4.82–6.69) (Fig. 3A), including four monotherapy and 13 combination therapy. The median PFS of monotherapy and combination therapy was 3.68 months (95% CI 3.45–3.92) and 6.45 months (95% CI 5.48–7.42) (Fig. 3B, C) respectively. For subgroups of patients who received ICIs plus TKIs and ICIs plus AI mAbs, the median PFS was 6.62 months (95% CI 5.50–7.74) vs 6.58 months (95% CI 4.95–8.22) (Fig. 3C, D). In the pooled analysis for PFS, all subgroups within the combination therapy group showed a median PFS exceeding the 6 months compared to the monotherapy group.

Fig. 3
figure 3

Pooled results of PFS

As showed in Fig. 4, six cohorts from five studies compared ICIs-based immunotherapy with TKIs, and all of them chose sorafenib as the control group [13, 33, 36, 50, 54]. The pooled results showed that in these four randomized control trials, the median PFS of ICIs-based group was 4.84 months (95% CI 3.65–6.03), and that of TKIs-based therapy was 3.72 months (95% CI 3.11–4.33). The P value between two regimens was 0.10, which meant no significant difference existed. It implied that the efficacy of ICIs-based therapy and TKIs-based was the same in prolonging PFS.

Fig. 4
figure 4

Pooled results of PFS of ICIs vs TKIs

3.4.2 Pooled results of OS

Only six clinical trials reported median OS. Two studies compared the efficacy of first-line and two-line treatment. One reported median OS of first-line group was 9.8 months (95% CI 3.7–14.3), and 7.4 months (95% CI 2.3–12.1) of second-line group indicating that first-line ICIs-based treatment could give a longer OS for HCC patients [32]. However, another study reported that the OS was similar between two groups in patients who had a complete/partial response, stable disease, or progressive disease. In addition, as summarized by the authors, unresectable HCC patients could benefit in OS from the use of ICIs-based treatment regardless first-line or second-line [42].

The pooled result of five studies was 16.35 months (95% CI 15.19–17.51) (Fig. 5). Due to the follow-up time was not long enough, it was immature to include the median OS from the other studies.

Fig. 5
figure 5

Pooled results of median OS

Three studies reported comparing ICIs-based treatment and TKIs [13, 33, 54]. As Fig. 6 showed, the median OS was 17.40 months (95% CI 15.45–19.36) vs 14.06 months (95% CI 13.11–15.01), and the p-value was smaller than 0.01. The results illustrated unresectable HCC patients could benefit in survival by choosing ICIs as their first-line treatment.

Fig. 6
figure 6

Pooled results of OS of ICIs vs TKIs

3.5 Best response rate

As shown in Supplementary Table S1, 25 groups reported ORR of the trials while 19 reported DCR of the patients. ORR and DCR were 25.1% (95% CI 20.8%-29.5%) and 75.2% (95% CI 70.3%-80.2%) as evaluated with response evaluation criteria in solid tumours (RECIST) v1.1. When measured with modified RECIST (mRECIST) v1.1, the pooled outcomes were 40.2% (95% CI 31.7%-48.6%) with 75.2% (95% CI 68.3%-82.1%).

We also calculated the pooled results for subgroups in Supplementary Table S1. When measured with RECIST v1.1, ORR and DCR of monotherapy was 16.1% (95% CI 13.7–18.5%) and 54.1% (95% CI 45.8–62.5%), and for the combination group, ORR and DCR were 27.8% (95% CI 22.8–32.8%) and 78.2% (95% CI 74.3–82.2%) respectively. When measured by mRECIST v1.1, because of the lack of results of monotherapy, we only got pooled results of combination therapy and its subgroups. The pooled result of ORR was 41.9% (95% CI 33.7–50.2%) for patients who received combination therapy, and DCR was 77.5% (95% CI 72.1–83.0%).

The pooled results of response rate revealed that ICIs plus TKIs group was more beneficial to unresectable HCC patients than ICIs plus AI mAbs regimens.

3.6 Pooled result of TRAEs

All 24 studies reported the rate of TRAEs. The incidence of any adverse events was 82.3% (95% CI 73.9%-90.7%) (Fig. 7A). Calculated according to different grades of adverse events, including grade1-2, grade3-4 and grade 5, the pooled results were 36.5% (95% CI 22.5–50.5%), 37.7% (95% CI 26.9–48.5%) and 3.6% (95% CI 0.6–6.7%) respectively (Fig. 7B, C, D). Moreover, we summarized the incidence of each adverse event in Supplementary Table S2. It was easy to discover that the most common adverse events differed in different grades. The most common adverse event of any grade was dysphonia, while the most common adverse in grade 1–2 and grade 3–4 were proteinuria and hypertension, respectively.

Fig. 7
figure 7

Pooled results of TRAE

3.7 Network meta-analysis

3.7.1 Network meta-analysis of ORR

Six clinical trials involved 3,741 patients had reported ORR. Based on these seven treatment options, we set up a network plot displayed as shown in Fig. 8. It clearly indicates that, compared to sorafenib as first-line treatment, all ICIs-based treatment regimens showed a significant improvement in ORR (Fig. 9).

Fig. 8
figure 8

Network plot of ORR

Fig. 9
figure 9

Forest plot for ORR compared to sorafenib

We also did further analysis to indirectly compare the results between every two regimens (Table 4). It showed the same conclusion as meta-analysis that using ICIs-based treatment as first-line therapy could benefit patients more than sorafenib.

Table 4 Results of network meta-analysis

3.7.2 Network meta-analysis of DCR

Different from ORR, only four clinical trials reported DCR. As seen in Fig. 10, four treatment options from three clinical studies were involved in this analysis of DCR. After using random effect model to evaluate the difference of DCR among four therapies, as shown in Fig. 11, combination therapy including atezolizumab plus bevacizumab 3.55 (95% CI 3.11–4.05) and sintilimab plus IBI305 3.09 (95% CI 2.58–3.71) both showed significant differences on DCR compared to sorafenib, while the result atezolizumab monotherapy was 0.91 (95% CI 0.63–1.32). The indirect results revealed that combination therapy would be a better choice than monotherapy.

Fig. 10
figure 10

Network plot of DCR

Fig. 11
figure 11

Forest plot for DCR compared to sorafenib

3.8 Publication bias

Publication bias existed in the analysis of median PFS via funnel plot shown in Supplementary Fig S1 and egger’s test shown in Supplementary Table S3, while no obvious publication bias of median OS was identified. For best responses, as shown in Supplementary Fig S2 and Supplementary Table S4, publication bias was observed in the pooled results of CR and ORR, and when measured by RECIST v1.1, the result of PR also showed a publication bias, while SD, PD and DCR did not show any significant publication bias as measured with funnel plots and Egger’s test.

4 Discussion

This systematic review, meta-analysis and network meta-analysis included 24 clinical trials, including 13 single-arms and 11 controlled trials, to evaluate the efficacy and safety of ICIs-based therapy as first-line treatment option for unresectable HCC patients. The results of this research indicated that the use of ICIs-based therapy could benefit unresectable HCC patients in both efficacy and safety. Our research findings also raise some issues that are worth discussing as below.

The NCCN guideline currently recommended 2 first-line treatment options including ICIs for specific unresectable HCC patients who provided category 1 evidence, while the CSCO guideline recommended 4 options. They were recommended because of the good efficacy and safety shown in clinical trials. Now, in this research, we have found that other ICIs-based therapies would also be promising as first-line treatments for patients. For example, sintilimab plus IBI305 are new drugs developed in recent years, of which the efficacy as first-line treatment exceeds all existing regimens underscoring the potential role of this combination therapy.

Before the clinical use of ICIs-based treatment, patients with unresectable HCC mostly accepted sorafenib as first-line therapy, which demonstrated a median OS of 11.4 months (95% CI 8.46–14.47) and median PFS of 4.17 months (95% CI 3.08–5.25) [55]. As a pervious review concluded, the potentiality of ICIs-based therapy brought big hope for both patients and researchers, and quite a lot clinical trials worked on different regimens based on their biological rationale [56]. We verified the great potential of ICIs-based in this study: compared with TKI therapy, ICIs-based therapy had significantly better efficacy in prolonging patients’ survival. Moreover, the combination therapy including TKIs and AI mAbs all demonstrated good efficacy. The anti-angiogenic effects of TKIs and AI mAbs could act synergistically with ICIs to destruct the vasculature of malignant tumours and thereby inhibit their growth [57, 58]. Further analysis revealed that different kinds of therapies had their own advantages: sintilimab plus IBI305 showed the best efficacy in ORR compared with other ICIs-based therapy, while atezolizumab plus bevacizumab performed better in terms of DCR. In general, the addition of ICIs as first-line treatment had a favourable impact on unresectable HCC patient survival. As another review suggested, positive outcomes for the use of ICIs-based therapies showed great potentials for developing more drugs against unresectable HCC [59].

Our research found that all different cohorts reported the rate of adverse events after taking ICIs-based treatment with dysphonia being most commonly reported across all grades of adverse events. For TRAEs of grade 1–2, the increased incidence of raised ALT and proteinuria suggests that clinicians need to pay attention to the patient's liver and kidney function during treatment. For TRAEs of grade 3–4, commonly reported hypertension and pyrexia warrant continuous monitoring of patients’ physical condition and adjustments in the treatment whenever necessary. Additionally, it is suggested that researchers could report the incidence of immune-related adverse events for a more detailed analysis.

Regarding clinical research design, it is worth noting that most clinical trials were open label, which might pertain high risks of bias leading to low reliability of the results. Also, 18 of 24 were single-arm trials or could not be used to build a meta-analysis, which meant no sufficient evidence to support the results to fully evaluate the performance of different treatment regimens. In addition, the follow-up time in some included studies was not long enough, resulting in data immature for calculating the median OS. As most existing studies were non-comparative, open-label, and small, it could lead to high heterogeneity and publication bias. Therefore, randomized controlled trials with masking and high-quality are urgently needed. What’s more, more comparative studies involving head-to-head comparisons are required to build a strong and reliable network to build a more robust conclusion about the advantages of different regimens. Since that the stage of liver function of most included participants were Child Pugh A, the limitation of participants may lead we missed the good efficacy for more HCC patients. just as other researchers did to verify the efficacy of metronomic capecitabine, we suggested do similar research for Child Pugh B patients [60, 61].

To our knowledge, this is the first study combined systematic review, meta-analysis and network meta-analysis about clinical trials of ICIs-based therapy as first-line treatment for unresectable HCC patients. However, there are some limitations worth attention. Firstly, because of the lack of aetiological data, we could not do further subgroup analysis to evaluate the efficacy and safety of ICIs-based therapy for patients of different conditions. Secondly, the quality of included clinical trials was evaluated as poor or fair because of the study design. While it is reasonable considering the urgent needs at clinical, it may lead to some bias in outcomes, especially in the long-term evaluation of OS. Future studies with the improvement of research design will be required to provide more extensive evidence of ICIs for HCC. Thirdly, this study only included clinical trials. It should be under consideration to include relevant real-world data to conduct a further and more reliable evaluation.

5 Conclusions

This study showed that ICIs-based therapy could benefit unresectable HCC patients in survival and tumour response rates with manageable adverse events. Furthermore, the combination regimens of atezolizumab plus bevacizumab and sintilimab plus IBI305 presented good outcomes in the evaluation, implying that combination therapy could become a new trend in clinical practice for unresectable HCC patients. However, due to the shortage of study design for existing clinical trials, more comparable and high-quality studies of ICIs-based regimens as first-line treatment are needed.