Comprehensive analysis of scRNA-Seq and bulk RNA-Seq data reveals dynamic changes in tumor-associated neutrophils in the tumor microenvironment of hepatocellular carcinoma and leads to the establishment of a neutrophil-related prognostic model.

Abstract

Analysis of hepatocellular carcinoma (HCC) single-cell sequencing data was conducted to explore the role of tumor-associated neutrophils in the tumor microenvironment. Analysis of single-cell sequencing data from 12 HCC tumor cores and five HCC paracancerous tissues identified cellular subpopulations and cellular marker genes. The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases were used to establish and validate prognostic models. xCELL, TIMER, QUANTISEQ, CIBERSORT, and CIBERSORT-abs analyses were performed to explore immune cell infiltration. Finally, the pattern of tumor-associated neutrophil roles in tumor microenvironmental components was explored. A total of 271 marker genes for tumor-associated neutrophils were identified based on single-cell sequencing data. Prognostic models incorporating eight genes were established based on TCGA data. Immune cell infiltration differed between the high- and low-risk groups. The low-risk group benefited more from immunotherapy. Single-cell analysis indicated that tumor-associated neutrophils were able to influence macrophage, NK cell, and T-cell functions through the IL16, IFN-II, and SPP1 signaling pathways. Tumor-associated neutrophils regulate immune functions by influencing macrophages and NK cells. Models incorporating tumor-associated neutrophil-related genes can be used to predict patient prognosis and immunotherapy responses.