Abstract:Background: Various immune cell types in the tumor microenvironment (TME) of hepatocellular carcinoma (HCC) have been identified as important parameters associated with prognosis and responsiveness to immunotherapy. However, how various factors influence immune cell infiltration remains incompletely understood. Hence, we investigated the single cell multi-omics landscape of immune infiltration in HCC, particularly key gene and cell subsets that influence immune infiltration, thus potentially linking the immunotherapy response and immune cell infiltration. Methods: We grouped patients with HCC according to immune cell infiltration scores calculated by single sample gene set enrichment analysis (ssGSEA). Differential expression analysis, functional enrichment, clinical trait association, gene mutation analysis, tumor immune dysfunction and exclusion (TIDE) and prognostic model construction were used to investigate the immune infiltration landscape through multi-omics. Stepwise regression was further used to identify key genes regulating immune infiltration. Single cell analysis was performed to explore expression patterns of candidate genes and investigate associated cellular populations. Correlation analysis, ROC analysis, Immunotherapy cohorts were used to explore and confirm the role of key gene and cellular population in predicting immune infiltration state and immunotherapy response. Immunohistochemistry and multiplexed fluorescence staining were used to further validated our results. Results: Patients with HCC were clustered into high and low immune infiltration groups. Mutations of CTNNB1 and TTN were significantly associated with immune infiltration and altered enrichment of cell populations in the TME. TIDE analysis demonstrated that T cell dysfunction and the T cell exclusion score were elevated in the high and low infiltration groups, respectively. Six risk genes and five risk immune cell types were identified and used to construct risk scores and a nomogram model. CXCR6 and LTA, identified by stepwise regression, were highly associated with immune infiltration. Single cell analysis revealed that LTA was expressed primarily in tumor infiltrating T lymphocytes and partial B lymphocytes, whereas CXCR6 was enriched predominantly in T and NK cells. Notably, CXCR6 + CD8 T cells were characterized as tumor enriched cells that may be potential predictors of high immune infiltration and the immune-checkpoint blockade response, and may serve as therapeutic targets. Conclusion: We constructed a comprehensive single cell and multi-omics landscape of immune infiltration in HCC, and delineated key genes and cellular populations regulating immune infiltration and immunotherapy response, thus providing insights into the mechanisms of immune infiltration and future therapeutic control.

Tumor-associated endothelial cell prognostic risk model and tumor immune environment modulation in liver cancer based on single-cell and bulk RNA sequencing: Experimental verification

Identification and Validation of a Novel Tumor Microenvironment-Related Prognostic Signature of Patients With Hepatocellular Carcinoma

Identification of a Novel Signature Based on Macrophage-Related Marker Genes to Predict Prognosis and Immunotherapeutic Effects in Hepatocellular Carcinoma.

Identification of Novel Tumor Microenvironment-Related Long Noncoding RNAs to Determine the Prognosis and Response to Immunotherapy of Hepatocellular Carcinoma Patients

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on NK cell marker genes to predict prognosis and immunotherapy response in hepatocellular carcinoma

Identification of a Distinct Tumor Endothelial Cell-Related Gene Expression Signature Associated with Patient Prognosis and Immunotherapy Response in Multiple Cancers.

Prognostic prediction signature and molecular subtype for liver cancer: A CTC/CTM‑related gene prediction model and independent external validation

Prognostic prediction using a gene signature developed based on exhausted T cells for liver cancer patients

A novel prognostic model based on single-cell RNA sequencing data for hepatocellular carcinoma

Integrated single cell and bulk sequencing analysis identifies tumor reactive CXCR6+ CD8 T cells as a predictor of immune infiltration and immunotherapy outcomes in hepatocellular carcinoma

Risk predictive model based on three immune-related gene pairs to assess prognosis and therapeutic sensitivity for hepatocellular carcinoma

Prognostic risk model of LIHC T-cells based on scRNA-seq and RNA-seq and the regulation of the tumor immune microenvironment

T cell-related prognostic risk model and tumor immune environment modulation in lung adenocarcinoma based on single-cell and bulk RNA sequencing

Identification of prognostic risk model based on plasma cell markers in hepatocellular carcinoma through single-cell sequencing analysis

Construction and validation of transcription‐factor‐based prognostic signature for TACE non‐response and characterization of tumor microenvironment infiltration in hepatocellular carcinoma

Identification and validation of immune-related gene signature models for predicting prognosis and immunotherapy response in hepatocellular carcinoma

Unveiling efferocytosis-related signatures through the integration of single-cell analysis and machine learning: a predictive framework for prognosis and immunotherapy response in hepatocellular carcinoma

A Novel Prognostic Model for Predicting Patient Survival and Immunotherapy Responsiveness in Hepatocellular Carcinoma: Insights into the Involvement of T-cell Proliferation.

Integrative study reveals the prognostic and immunotherapeutic value of CTCs/CTMs -related prognostic signature in Hepatocellular carcinoma

Identification and validation of a tyrosine metabolism-related prognostic prediction model and characterization of the tumor microenvironment infiltration in hepatocellular carcinoma

Establishment and Validation of a Novel CD8+ T Cell-Associated Prognostic Signature for Predicting Clinical Outcomes and Immunotherapy Response in Hepatocellular Carcinoma Via Integrating Single-Cell RNA-seq and Bulk RNA-seq.