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.

Potential crosstalk between Naïve CD4+ T cells and SPP1+ Macrophages is associated with clinical outcome and therapeutic response in hepatocellular carcinoma

Identification of Molecular Subtypes for Integrated Multi-Omics Analysis for Use in Guiding Precision Medicine in Hepatocellular Carcinoma.

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

IDDF2019-ABS-0174 Targeting Monocyte-Intrinsic Enhancer Reprogramming Improves Immunotherapy Efficacy in Hepatocellular Carcinoma

Single-cell transcriptomics reveals the role of Macrophage-Naïve CD4 + T cell interaction in the immunosuppressive microenvironment of primary liver carcinoma

Identification of SPP1+ macrophages as an immune suppressor in hepatocellular carcinoma using single-cell and bulk transcriptomics

Multi-omics analysis of macrophage-associated receptor and ligand reveals a strong prognostic signature and subtypes in hepatocellular carcinoma

Characterizing ligand-receptor interactions and unveiling the pro-tumorigenic role of CCL16-CCR1 axis in the microenvironment of hepatocellular carcinoma

Construction of TME and Identification of crosstalk between malignant cells and macrophages by SPP1 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.

Integrated Analysis Reveals the Microenvironment of Non-Small Cell Lung Cancer and a Macrophage-Related Prognostic Model

Single cell analyses reveal the PD-1 blockade response-related immune features in hepatocellular carcinoma

Identification of M2-like macrophage-related signature for predicting the prognosis, ecosystem and immunotherapy response in 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

Identification of a Tumour Immune Barrier in the HCC Microenvironment That Determines the Efficacy of Immunotherapy

Comprehensive scRNA-seq Analysis and Identification of CD8_+T Cell Related Gene Markers for Predicting Prognosis and Drug Resistance of Hepatocellular Carcinoma

The immune-metabolic crosstalk between CD3+C1q+TAM and CD8+T cells associated with relapse-free survival in HCC

Single-cell RNA sequencing reveals cell subpopulations in the tumor microenvironment contributing to hepatocellular carcinoma

Single-cell RNA sequencing revealed potential targets for immunotherapy studies in hepatocellular carcinoma

Cellular crosstalk of macrophages and therapeutic implications in non-small cell lung cancer revealed by integrative inference of single-cell transcriptomics

POSTN + cancer-associated fibroblasts determine the efficacy of immunotherapy in hepatocellular carcinoma