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.

CD69 + Vδ1γδ T cells are anti-tumor subpopulations in hepatocellular carcinoma

Hepatocellular Carcinoma‐infiltrating Γδ T Cells Are Functionally Defected and Allogenic Vδ2+ Γδ T Cell Can Be a Promising Complement

Limited Predictive or Prognostic Role of Tumor-Infiltrating Tissue-Resident Memory CD8 T Cells in Patients with Hepatocellular Carcinoma Receiving Immunotherapy

Increased Cd4(+)Cd69(+)Cd25(-) T Cells In Patients With Hepatocellular Carcinoma Are Associated With Tumor Progression

CD8+CXCR5+T cells infiltrating hepatocellular carcinomas are activated and predictive of a better prognosis.

Expression of TIGIT/CD155 and Correlations with Clinical Pathological Features in Human Hepatocellular Carcinoma

The Unique Characteristic in Peripheral Immune Cells in Patients with Advanced Hepatocellular Carcinoma.

The Identification and Functional Analysis of CD8+PD-1+CD161+ T Cells in Hepatocellular Carcinoma.

The Functional Impairment of HCC-infiltrating Γδ T Cells, Partially Mediated by Regulatory T Cells in a TGFβ- and IL-10-dependent Manner

Intrahepatic CD69+Vδ1 T cells re-circulate in the blood of patients with metastatic colorectal cancer and limit tumor progression

Decreased granzyme-B expression in CD11c+CD8+ T cells associated with disease progression in patients with HBV-related hepatocellular carcinoma

Increased Frequency of CD4+CD25highFoxP3+ Regulatory T Cells in Patients with Hepatocellular Carcinoma

Low Counts of Γδ T Cells in Peritumoral Liver Tissue Are Related to More Frequent Recurrence in Patients with Hepatocellular Carcinoma after Curative Resection.

Genomic and Transcriptional Profiling of Tumor Infiltrated CD8+ T Cells Revealed Functional Heterogeneity of Antitumor Immunity in Hepatocellular Carcinoma.

Identification of a FOXP3(+)CD3(+)CD56(+) Population with Immunosuppressive Function in Cancer Tissues of Human Hepatocellular Carcinoma.

The Mechanisms of Gamma Delta T Cells in the Occurrence and Development of Hepatocellular Carcinoma

CD69 is a Promising Immunotherapy and Prognosis Prediction Target in Cancer

Single-cell RNA sequencing highlights the role of PVR/PVRL2 in the immunosuppressive tumour microenvironment in hepatocellular carcinoma

PD1Hi CD8+ T Cells Correlate with Exhausted Signature and Poor Clinical Outcome in Hepatocellular Carcinoma

Phenotype and Function of Cxcr5(+)Cd45ra(-)Cd4(+) T Cells Were Altered in Hbv-Related Hepatocellular Carcinoma and Elevated Serum Cxcl13 Predicted Better Prognosis

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