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.

Integrating bulk and single-cell RNA sequencing reveals promotes hepatocellular carcinoma progression by activating the PI3K/AKT/mTOR pathway

Oncogene TUBA1C Promotes Migration and Proliferation in Hepatocellular Carcinoma and Predicts a Poor Prognosis.

Single-cell transcriptome analysis reveals the metabolic changes and the prognostic value of malignant hepatocyte subpopulations and predict new therapeutic agents for hepatocellular carcinoma

Single-cell transcriptomics reveals tumor microenvironment changes and prognostic gene signatures in hepatocellular carcinoma

Study on the prediction model of liver cancer based on chronic liver disease and the related molecular mechanism

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

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

Integrating Bulk and Single‐cell RNA Sequencing Reveals Cellular Heterogeneity and Immune Infiltration in Hepatocellular Carcinoma

Integrative Analysis Identifies Cell-Type-Specific Genes Within Tumor Microenvironment as Prognostic Indicators in Hepatocellular Carcinoma

Clinical and Prognostic Value of PPIA, SQSTM1, and CCL20 in Hepatocellular Carcinoma Patients by Single-Cell Transcriptome Analysis

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

Integrating TCGA and Single-Cell Sequencing Data for Hepatocellular Carcinoma: A Novel Glycosylation (GLY)/Tumor Microenvironment (TME) Classifier to Predict Prognosis and Immunotherapy Response

Integrated clinical and prognostic analyses of mTOR/Hippo pathway core genes in hepatocellular carcinoma

Large-scale analyses identify a cluster of novel long noncoding RNAs as potential competitive endogenous RNAs in progression of hepatocellular carcinoma

Development of an Autophagy-Based and Stemness-Correlated Prognostic Model for Hepatocellular Carcinoma Using Bulk and Single-Cell RNA-Sequencing

Prognostic value and potential mechanism of cellular senescence and tumor microenvironment in hepatocellular carcinoma: Insights from bulk transcriptomics and single‐cell sequencing analysis

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

Single-cell analysis reveals the intra-tumor heterogeneity and identifies MLXIPL as a biomarker in the cellular trajectory of hepatocellular carcinoma

Identifying epithelial-mesenchymal transition-related genes as prognostic biomarkers and therapeutic targets of hepatocellular carcinoma by integrated analysis of single-cell and bulk-RNA sequencing data

A Novel Approach for the Analysis of Single-Cell RNA Sequencing Identifies TMEM14B As a Novel Poor Prognostic Marker in Hepatocellular Carcinoma.