Immunosuppressive landscape in hepatocellular carcinoma revealed by single-cell sequencing

Natural killer cell-based cancer immunotherapy: from basics to clinical trials

Cellular immunotherapy exploits the capacity of the human immune system in self-protection and surveillance to achieve the anti-tumor effects. Natural killer (NK) cells are lymphocytes of innate immune system and they display a unique inherent ability to identify and eliminate tumor cells. In this review, we first introduce the basic characteristics of NK cells in the physiological and pathological milieus, followed by a discussion of their effector function and immunosuppression in the tumor microenvironment. Clinical strategies and reports regarding NK cellular therapy are analyzed in the context of tumor treatment, especially against solid tumors. Given the widely studied T-cell therapy in the recent years, particularly the chimeric antigen receptor (CAR) T-cell therapy, we compare the technical features of NK- and T-cell based tumor therapies at the clinical front. Finally, the technical challenges and potential solutions for both T and NK cell-based immunotherapies in treating tumor malignancies are delineated. By overviewing its clinical applications, we envision the NK-cell based immunotherapy as an up-and-comer in cancer therapeutics.

Significance of Immune and Non-Immune Cell Stroma as a Microenvironment of Hepatocellular Carcinoma-From Inflammation to Hepatocellular Carcinoma Progression

Hepatocellular carcinoma (HCC) is the most common liver cancer as well as the most prevalent cause of death in the adult patient population with cirrhosis. The occurrence of HCC is primarily caused by chronic liver inflammation that might occur because of a viral infection, non-alcoholic fatty liver disease (NAFLD), or various lifestyle-associated factors. The objective of this review was to summarize the current knowledge regarding the microenvironment of HCC, indicating how immune- and non-immune-cell stroma might affect the onset and progression of HCC. Therefore, in the following narrative review, we described the role of tumor-infiltrating neutrophils, bone-marrow-derived cells, tumor-associated mast cells, cancer-associated fibroblasts, tumor-associated macrophages, liver-sinusoidal endothelial cells, lymphocytes, and certain cytokines in liver inflammation and the further progression to HCC. A better understanding of the HCC microenvironment might be crucial to introducing novel treatment strategies or combined therapies that could lead to more effective clinical outcomes.

Enhancing cancer immunotherapy using cordycepin and Cordyceps militaris extract to sensitize cancer cells and modulate immune responses

Integrating immunotherapy with natural compounds holds promise in enhancing the immune system's ability to eliminate cancer cells. Cordyceps militaris, a traditional Chinese medicine, emerges as a promising candidate in this regard. This study investigates the effects of cordycepin and C. militaris ethanolic extract (Cm-EE) on sensitizing cancer cells and regulating immune responses against breast cancer (BC) and hepatocellular carcinoma (HCC) cells. Cordycepin, pentostatin and adenosine were identified in Cm-EE. Cordycepin treatment decreased HLA-ABC-positive cells in pre-treated cancer cells, while Cm-EE increased NKG2D ligand and death receptor expression. Additionally, cordycepin enhanced NKG2D receptor and death ligand expression on CD3-negative effector immune cells, particularly on natural killer (NK) cells, while Cm-EE pre-treatment stimulated IL-2, IL-6, and IL-10 production. Co-culturing cancer cells with effector immune cells during cordycepin or Cm-EE incubation resulted in elevated cancer cell death. These findings highlight the potential of cordycepin and Cm-EE in improving the efficacy of cancer immunotherapy for BC and HCC.

Current research insights into the role of CTLA-4 in hepatitis B virus (HBV) infection

Chronic hepatitis B virus (HBV) infection is a significant global public health concern, and the clearance of HBV is closely linked to the activity of HBV-specific T cells, which is regulated by various co-suppressor molecules. Cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) is among these co-suppressor molecules which induces T cell exhaustion by competitively inhibiting CD28 and dampening the function of HBV-specific T cells. CTLA-4 also plays a role in the regulation of T helper (Th) cell differentiation and influences cytokine release. In addition, CTLA-4 can impact glucose metabolism in hepatocellular carcinoma through its interaction with T regulatory (Treg) cells. This review aims to provide a comprehensive overview of the existing literature related to the role of CTLA-4 in HBV patients across different subsets of T cells. Additionally, we propose a discussion on the possible mechanisms through which CTLA-4 may contribute to HBV infection, as well as the development of HBV-induced cirrhosis and hepatocellular carcinoma.

CD74 facilitates immunotherapy response by shaping the tumor microenvironment of hepatocellular carcinoma

BACKGROUND

CD74 is ectopically expressed in many tumors and can regulate tumor immunity. However, there are many gaps in the study of the prognostic value of CD74 expression and immune infiltration in hepatocellular carcinoma (HCC).

METHODS

An online tumor database was searched to obtain data on gene/protein expression. Immune infiltration analysis was performed using the Tumor Immune Estimation Resource and Comprehensive Analysis on Multi-Omics of Immunotherapy in Pan-cancer databases. Single-cell data were obtained from the Tissue-specific Gene Expression and Regulation, Single-cell Transcriptomes of Tumor Immune Microenvironment and Tumor Immune Single-cell Hub 2 databases.

RESULTS

CD74 was highly expressed in HCC patients. HCC patients with high CD74 expression who consumed alcohol or were negative for hepatitis virus had a better prognosis than patients with low CD74 expression. CD74 was mainly enriched in immune response regulation pathways. Both copy number variations in CD74 and CD74 expression patterns affected the infiltration levels of immune cells. Interestingly, CD74 regulated the differentiation of myeloid cells. CD74 in macrophages and dendritic cells (DCs) forms complex networks with malignant cells and hepatic progenitor cell (HPC)-like cells, respectively. High CD74 expression in HPC-like cells and malignant cells significantly decreased the fraction of C-type lectin domain family 9 A (CLEC9A)-cDC1+ DCs and IL-1B+ macrophages, respectively. Their crosstalk subsequently shaped the tumor microenvironment of HCC, possibly through the CD74-MIF axis. Importantly, patients with high CD74 expression presented higher immune scores and achieved good outcomes after receiving immunotherapy.

CONCLUSION

High CD74 expression is associated with the abundance of a variety of immune cell types, mediating interactions among tumor and immune cells and shaping the malignant behavior of HCC. In summary, CD74 may be a hallmark for determining the prognosis and immune cell infiltration levels of HCC patients.

Immunosuppressive tumor microenvironment in the progression, metastasis, and therapy of hepatocellular carcinoma: from bench to bedside

Hepatocellular carcinoma (HCC) is a highly heterogeneous malignancy with high incidence, recurrence, and metastasis rates. The emergence of immunotherapy has improved the treatment of advanced HCC, but problems such as drug resistance and immune-related adverse events still exist in clinical practice. The immunosuppressive tumor microenvironment (TME) of HCC restricts the efficacy of immunotherapy and is essential for HCC progression and metastasis. Therefore, it is necessary to elucidate the mechanisms behind immunosuppressive TME to develop and apply immunotherapy. This review systematically summarizes the pathogenesis of HCC, the formation of the highly heterogeneous TME, and the mechanisms by which the immunosuppressive TME accelerates HCC progression and metastasis. We also review the status of HCC immunotherapy and further discuss the existing challenges and potential therapeutic strategies targeting immunosuppressive TME. We hope to inspire optimizing and innovating immunotherapeutic strategies by comprehensively understanding the structure and function of immunosuppressive TME in HCC.

CanCellVar: A database for single-cell variants map in human cancer

Numerous variants, including both single-nucleotide variants (SNVs) in DNA and A>G RNA edits in mRNA as essential drivers of cellular proliferation and tumorigenesis, are commonly associated with cancer progression and growth. Thus, mining and summarizing single-cell variants will provide a refined and higher-resolution view of cancer and further contribute to precision medicine. Here, we established a database, CanCellVar, which aims to provide and visualize the comprehensive atlas of single-cell variants in tumor microenvironment. The current CanCellVar identified ∼3 million variants (∼1.4 million SNVs and ∼1.4 million A>G RNA edits) involved in 2,754,531 cells of 5 major cell types across 37 cancer types. CanCellVar provides the basic annotation information as well as cellular and molecular function properties of variants. In addition, the clinical relevance of variants can be obtained including tumor grade, treatment, metastasis, and others. Several flexible tools were also developed to aid retrieval and to analyze cell-cell interactions, gene expression, cell-development trajectories, regulation, and molecular structure affected by variants. Collectively, CanCellVar will serve as a valuable resource for investigating the functions and characteristics of single-cell variations and their roles in human tumor evolution and treatment.

IMbrave152/SKYSCRAPER-14: a Phase III study of atezolizumab, bevacizumab and tiragolumab in advanced hepatocellular carcinoma

Atezolizumab plus bevacizumab is a standard of care, first-line therapy for advanced hepatocellular carcinoma (HCC). Myeloid and T regulatory cells are key immunosuppressive cell types within the hepatic tumor microenvironment associated with clinical resistance to atezolizumab and bevacizumab therapy for HCC and overall poor prognosis. Therapeutic targeting of TIGIT, which is highly expressed in these cells, with tiragolumab may overcome the immunosuppressive environment and improve clinical benefit, a hypothesis supported by positive efficacy signals in the Phase Ib/II MORPHEUS-Liver study. This paper describes the rationale and design of IMbrave152/SKYSCRAPER-14, a randomized, double-blind, placebo-controlled Phase III study comparing atezolizumab and bevacizumab with tiragolumab or placebo in patients with HCC and no prior systemic treatment.Clinical Trial Registration: NCT05904886 (ClinicalTrials.gov).

Single cell analysis unveils the commonality and heterogeneity between nasopharyngeal and oropharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) and oropharyngeal carcinoma (OPC) are subtypes of head and neck cancer with different treatment effects due to the heterogeneity of tumor microenvironments. This study was to investigate the distinctive tumor microenvironments of NPC and OPC. Analyzing single-cell data from 10 cases of each subtype, we reveal significant differences in cellular composition, with NPC microenvironment dominated by T/NK and B cells, and OPC characterized by prevalent epithelial cells and fibroblasts. Dynamic transitions of CD8 T cells are observed in both tumor types, involving shifts from naivety to cytotoxicity, proliferation, and eventual exhaustion/exhausted states. Additionally, Tregs exhibit heightened proliferative abilities in later developmental stages, concomitant with exhaustion. These highly proliferative T cells and Tregs manifest elevated glycolysis and lactate metabolism activities. Furthermore, we explore intercellular communication between glycolytic malignant epithelial cells and these proliferative T cells. These findings offer comprehensive insights into the heterogeneity of tumor microenvironments and provide a solid foundation for future therapeutic strategies and targeted interventions.

Regulatory T cell: a double-edged sword from metabolic-dysfunction-associated steatohepatitis to hepatocellular carcinoma

Metabolic-dysfunction-associated steatotic liver disease (MASLD) is becoming a leading cause of end-stage liver disease globally. Metabolic-dysfunction-associated steatohepatitis (MASH) represents a progressive inflammatory manifestation of MASLD. MASH underlies a versatile and dynamic inflammatory microenvironment, accompanied by aberrant metabolism and ongoing liver regeneration, establishing itself as a significant risk factor for hepatocellular carcinoma (HCC). The mechanisms underlying the escape and survival of malignant cells within the extensive inflammatory microenvironment of MASH remain elusive. Regulatory T cells (Tregs) play a crucial role in maintaining homeostasis and preventing excessive immune responses in the liver. Paradoxically, Tregs have been implicated in inhibiting tumour-promoting inflammation and facilitating the evasion of cancer cells. Recent studies have unveiled distinct behaviours of Tregs at different stages of MASLD, suggesting a dual role in the pathogenesis. In this review, we explore the fate of Tregs from MASLD to HCC, offering recent insights into potential targets for clinical intervention.

Single-Cell RNA Sequencing Revealed That the Enrichment of TPI1+ Malignant Hepatocytes Was Linked to HCC Metastasis and Immunosuppressive Microenvironment

Background

Tumor metastasis is the leading cause of high mortality in hepatocellular carcinoma (HCC). The metastasis-related HCC microenvironment is characterized by high heterogeneity. Single-cell RNA sequencing (scRNA-seq) may aid in determining specific cell clusters involved in regulating the immune microenvironment of HCC.

Methods

The scRNA-seq data of 10 HCC samples were collected from the Gene Expression Omnibus (GEO) database GSE124395. Correlations between key gene expression and clinicopathological data were determined using public databases. HCC tissues and matched tumor-adjacent and normal tissue samples were obtained by surgical resection at Sichuan Cancer Hospital. Immune cell infiltration analysis was performed and verified by immunohistochemistry and immunofluorescent staining.

Results

Nine malignant hepatocyte clusters with different marker genes and biological functions were identified. C3_Hepatocyte-SERF2 and C6_Hepatocyte-IL13RA2 were mainly involved in the regulation of the immune microenvironment, which was also a significant pathway in regulating HCC metastasis. Key genes in malignant hepatocyte clusters that associated with HCC metastasis were further screened by LASSO regression analysis. TPI1, a key gene in C6_Hepatocyte-IL13RA2 and HCC metastasis, could participate in regulating the HCC immune microenvironment in The Cancer Genome Atlas (TCGA) and Tumor Immune Estimation Resource (TIMER) databases. Moreover, immunohistochemistry analysis demonstrated that TPI1 expression was positively correlated with HCC metastasis and poor prognosis, while negatively correlated with CD8+ T cell infiltration. The negative correlation between TPI1 expression and CD8+ T cell infiltration was further confirmed by immunofluorescence staining.

Conclusion

In summary, a cluster of TPI1+ malignant hepatocytes was associated with the suppression of CD8+ T cell infiltration and HCC metastasis, providing novel insights into potential biomarkers for immunotherapy in HCC.

Key oncogenic signaling pathways affecting tumor-infiltrating lymphocytes infiltration in hepatocellular carcinoma: basic principles and recent advances

The incidence of hepatocellular carcinoma (HCC) ranks first among primary liver cancers, and its mortality rate exhibits a consistent annual increase. The treatment of HCC has witnessed a significant surge in recent years, with the emergence of targeted immune therapy as an adjunct to early surgical resection. Adoptive cell therapy (ACT) using tumor-infiltrating lymphocytes (TIL) has shown promising results in other types of solid tumors. This article aims to provide a comprehensive overview of the intricate interactions between different types of TILs and their impact on HCC, elucidate strategies for targeting neoantigens through TILs, and address the challenges encountered in TIL therapies along with potential solutions. Furthermore, this article specifically examines the impact of oncogenic signaling pathways activation within the HCC tumor microenvironment on the infiltration dynamics of TILs. Additionally, a concise overview is provided regarding TIL preparation techniques and an update on clinical trials investigating TIL-based immunotherapy in solid tumors.

Single-cell RNA sequencing reveals the immune microenvironment landscape of osteosarcoma before and after chemotherapy

Chemotherapy, a primary treatment for osteosarcoma (OS), has limited knowledge regarding its impact on tumor immune microenvironment (TIME). Here, tissues from 6 chemotherapy-naive OS patients underwent single-cell RNA sequencing (scRNA-seq) and were analyzed alongside public dataset (GSE152048) containing 7 post-chemotherapy OS tissues. CD45+ (PTPRC+) cells were used for cell clustering and annotation. Changes in immune cell composition pre- and post-chemotherapy were characterized. Totally, 28,636 high-quality CD45+ (PTPRC+) cells were extracted. Following chemotherapy, the proportions of regulatory T cells (Tregs) and activated CD8 T cells decreased, while CD8 effector T cells increased. GO analysis indicated that differentially expressed genes (DEGs) in T cells were associated with cell activation, adaptive immune response, and immune response to tumor cells. Furthermore, the proportions of plasma cells increased, while naive B cells decreased. B cell surface receptors expression was upregulated, and GO analysis revealed DEGs of B cells were mainly enriched in B cell-mediated immunity and B cell activation. Moreover, M2 polarization of macrophages was suppressed post-chemotherapy. Overall, this study elucidates chemotherapy remodels the OS TIME landscape, triggering immune heterogeneity and enhancing anti-tumor properties.

Crosstalk of cuproptosis-related prognostic signature and competing endogenous RNAs regulation in hepatocellular carcinoma

BACKGROUND

Cuproptosis is a new type of programmed cell death involved in the regulation of neuroendocrine tumors, immune microenvironment, and substance metabolism. However, the role of cuproptosis-related genes (CRGs) in Hepatocellular carcinoma (HCC) remains unclear.

METHOD

Through multiple bioinformatics analysis, we constructed a prognostic gene model and competing endogenous RNA (ceRNA) network. The correlation between CRGs and prognosis, immune infiltration, immune checkpoints, microsatellite instability (MSI) and tumor mutational burden (TMB) was analyzed by Kaplan-Meier curve, univariate Cox, multivariate regression, and Spearman's analysis in HCC patients. Besides, the qRT-PCR and immunohistochemistry assays were used to determine prognostic CRGs mRNA and protein expression in HCC.

RESULTS

We established a novel 3-gene signature related to CRGs for evaluating the prognosis of HCC patients. HCC patients with high risk scores had a poor prognosis with an area under the curve of 0.737, 0.646, and 0.634 on 1-year, 3-year, and 5-year receiver operating characteristic curves. Significant correlation was observed between prognostic CRGs and immune infiltration, immune checkpoints, MSI and TMB. We also developed five ceRNA networks to regulate the occurrence and progression of HCC. CDKN2A, DLAT, and PDHA1 protein expression was up-regulated in HCC versus normal tissues. Besides, the mRNA expression levels of CDKN2A, DLAT, GLS, and PDHA1 were elevated in the HCC cell lines compared to the normal liver cell lines.

CONCLUSIONS

This novel prognostic CRGs signature could be accurately predict the prognosis of patients with HCC. The ceRNA regulatory network might be potential prognostic biomarkers and therapeutic targets for HCC patients.

Molecular and immune landscape of hepatocellular carcinoma to guide therapeutic decision-making

Liver cancer, primarily HCC, exhibits highly heterogeneous histological and molecular aberrations across tumors and within individual tumor nodules. Such intertumor and intratumor heterogeneities may lead to diversity in the natural history of disease progression and various clinical disparities across the patients. Recently developed multimodality, single-cell, and spatial omics profiling technologies have enabled interrogation of the intertumor/intratumor heterogeneity in the cancer cells and the tumor immune microenvironment. These features may influence the natural history and efficacy of emerging therapies targeting novel molecular and immune pathways, some of which had been deemed undruggable. Thus, comprehensive characterization of the heterogeneities at various levels may facilitate the discovery of biomarkers that enable personalized and rational treatment decisions, and optimize treatment efficacy while minimizing the risk of adverse effects. Such companion biomarkers will also refine HCC treatment algorithms across disease stages for cost-effective patient management by optimizing the allocation of limited medical resources. Despite this promise, the complexity of the intertumor/intratumor heterogeneity and ever-expanding inventory of therapeutic agents and regimens have made clinical evaluation and translation of biomarkers increasingly challenging. To address this issue, novel clinical trial designs have been proposed and incorporated into recent studies. In this review, we discuss the latest findings in the molecular and immune landscape of HCC for their potential and utility as biomarkers, the framework of evaluation and clinical application of predictive/prognostic biomarkers, and ongoing biomarker-guided therapeutic clinical trials. These new developments may revolutionize patient care and substantially impact the still dismal HCC mortality.

The role of myeloid-derived suppressor cells in liver cancer

MDSCs are immature myeloid immune cells, which accumulate in models of liver cancer to reduce effector immune cell activity, contribute to immune escape and treatment resistance. The accumulation of MDSCs suppresses the role of CTL and the killing effects of NK cells, induces the accumulation of Treg cells, and blocks the antigen presentation of DCs, thus promoting the progression of liver cancer. Recently, immunotherapy has emerged a valuable approach following chemoradiotherapy in the therapy of advanced liver cancer. A considerable increasing of researches had proved that targeting MDSCs has become one of the therapeutic targets to enhance tumor immunity. In preclinical study models, targeting MDSCs have shown encouraging results in both alone and in combination administration. In this paper, we elaborated immune microenvironment of the liver, function and regulatory mechanisms of MDSCs, and therapeutic approaches to target MDSCs. We also expect these strategies to supply new views for future immunotherapy for the treatment of liver cancer.

Decoding Immune Signature to Detect the Risk for Early-Stage HCC Recurrence

Hepatocellular carcinoma (HCC) is often recognized as an inflammation-linked cancer, which possesses an immunosuppressive tumor microenvironment. Curative treatments such as surgical resection, liver transplantation, and percutaneous ablation are mainly applicable in the early stage and demonstrate significant improvement of survival rate in most patients. However, 70-80% of patients report HCC recurrence within 5 years of curative treatment, representing an important clinical issue. However, there is no effective recurrence marker after surgical and locoregional therapies, thus, tumor size, number, and histological features such as cancer cell differentiation are often considered as risk factors for HCC recurrence. Host immunity plays a critical role in regulating carcinogenesis, and the immune microenvironment characterized by its composition, functional status, and density undergoes significant alterations in each stage of cancer progression. Recent studies reported that analysis of immune contexture could yield valuable information regarding the treatment response, prognosis and recurrence. This review emphasizes the prognostic value of tumors associated with immune factors in HCC recurrence after curative treatment. In particular, we review the immune landscape and immunological factors contributing to early-stage HCC recurrence, and discuss the immunotherapeutic interventions to prevent tumor recurrence following curative treatments.

Tumor-Infiltrating B Lymphocytes: Promising Immunotherapeutic Targets for Primary Liver Cancer Treatment

Hepatocellular carcinoma and cholangiocarcinoma are the fourth most lethal primary cancers worldwide. Therefore, there is an urgent need for therapeutic strategies, including immune cell targeting therapies. The heterogeneity of liver cancer is partially explained by the characteristics of the tumor microenvironment (TME), where adaptive and innate immune system cells are the main components. Pioneering studies of primary liver cancers revealed that tumor-infiltrating immune cells and their dynamic interaction with cancer cells significantly impacted carcinogenesis, playing an important role in cancer immune evasion and responses to immunotherapy treatment. In particular, B cells may play a prominent role and have a controversial function in the TME. In this work, we highlight the effect of B lymphocytes as tumor infiltrates in relation to primary liver cancers and their potential prognostic value. We also present the key pathways underlying B-cell interactions within the TME, as well as the way that a comprehensive characterization of B-cell biology can be exploited to develop novel immune-based therapeutic approaches.

TIGIT Expression on Intratumoral Lymphocytes Correlates with Improved Prognosis in Oral Squamous Cell Carcinoma

(1) Background: T-cell immunoglobulin and ITIM domain (TIGIT) is a potential immunotherapeutic target in a variety of malignant entities, and antibody-based treatments are currently under investigation in clinical trials. While promising results were observed in patients with lung cancer, the role of TIGIT in oral squamous cell carcinoma (OSCC) as a biomarker as well as a therapeutic target remains elusive. Therefore, we evaluated the role of TIGIT as a prognostic factor in OSCC. (2) Methods: Here, we describe the results of a retrospective tissue microarray (TMA) OSCC cohort. Using immunohistochemistry, TIGIT expression was correlated with overall and recurrence-free survival (OAS and RFS, respectively). Additionally, in silico analysis was performed based on the TCGA Head and Neck Squamous Cell Carcinoma (HNSCC) cohort in order to correlate patients' survival with TIGIT and CD274 (encoding for PD-L1) gene expression levels. (3) Results: Database analysis revealed a beneficial outcome in OAS for tumor patients with high intraepithelial CD3-TIGIT-expression (n = 327). Hereby, OAS was 53.9 months vs. 30.1 months for patients with lower TIGIT gene expression levels (p = 0.033). In our retrospective OSCC-TMA cohort, elevated TIGIT levels on CD3+ cells correlated significantly with improved OAS (p = 0.025) as well as distant RFS (p = 0.026). (4) Conclusions: This study introduces TIGIT as a novel prognostic factor in OSCC, indicating the improved outcome of OSCC patients relative to their increased TIGIT expression. TIGIT might provide therapeutic implications for future immunotherapy in advanced-stage OSCC patients.

Sorafenib, Lenvatinib, or Lenvatinib Combining PD-1 Inhibitors Plus TACE in Unresectable Hepatocellular Carcinoma: A Retrospective Analysis

Introduction: This retrospective study aimed to compare the efficacy and safety of transarterial chemoembolization plus lenvatinib and programmed death 1 (PD-1) inhibitors versus transarterial chemoembolization plus lenvatinib or sorafenib in patients with unresectable hepatocellular carcinoma. Methods: Consecutive patients with unresectable hepatocellular carcinoma who received transarterial chemoembolization plus lenvatinib and PD-1 inhibitors, lenvatinib, or sorafenib were retrospectively identified in our institution between January 2018 and August 2020. The primary endpoint was overall survival. Results: A total of 84 patients were included in this analysis. The median overall survival was significantly improved in the transarterial chemoembolization plus lenvatinib and PD-1 inhibitor group compared with the transarterial chemoembolization plus sorafenib group (26.7 months [95% confidence interval 25.2-31.6] vs 14.4 months [95% confidence interval 9.5-18.9]; hazard ratio 0.39 [95% confidence interval 0.17-0.72]; P  =  .007) or the transarterial chemoembolization plus lenvatinib group (26.7 months [95% confidence interval 25.2-31.6] vs 17.9 [95% confidence interval 13.4-22.2] months; hazard ratio 0.45 [95% confidence interval 0.17-0.87]; P  =  .031). Transarterial chemoembolization plus lenvatinib and PD-1 inhibitor also significantly prolonged median progression-free survival compared with transarterial chemoembolization plus sorafenib group (8.2 months [95% confidence interval 3.3-13.0] vs 6.0 months [95% confidence interval 4.2-7.8]; hazard ratio 0.47 [95% confidence interval 0.24-0.74]; P  =  .005) or the transarterial chemoembolization plus lenvatinib group (8.2 months [95% confidence interval 3.3-13.0] vs 6.6 [95% confidence interval 4.3-7.9] months; hazard ratio 0.58 [95% confidence interval 0.31-0.96]; P  =  .047). No significant difference was seen between groups in the incidence of an adverse event or grade 3 or higher adverse event. Conclusion: Transarterial chemoembolization plus lenvatinib, and PD-1 inhibitor was associated with better survival benefits and acceptable toxicities, which may provide an additional therapeutic option for unresectable hepatocellular carcinoma.