RESUMEN
The major liver cancer subtype is hepatocellular carcinoma (HCC). Studies have indicated that a better prognosis is related to the presence of tumor-infiltrating lymphocytes (TILs) in HCC. However, the molecular pathways that drive immune cell variation in the tumor microenvironment (TME) remain poorly understood. Glycosylation (GLY)-related genes have a vital function in the pathogenesis of numerous tumors, including HCC. This study aimed to develop a GLY/TME classifier based on glycosylation-related gene scores and tumor microenvironment scores to provide a novel prognostic model to improve the prediction of clinical outcomes. The reliability of the signatures was assessed using receiver operating characteristic (ROC) and survival analyses and was verified with external datasets. Furthermore, the correlation between glycosylation-related genes and other cells in the immune environment, the immune signature of the GLY/TME classifier, and the efficacy of immunotherapy were also investigated. The GLY score low/TME score high subgroup showed a favorable prognosis and therapeutic response based on significant differences in immune-related molecules and cancer cell signaling mechanisms. We evaluated the prognostic role of the GLY/TME classifier that demonstrated overall prognostic significance for prognosis and therapeutic response before treatment, which may provide new options for creating the best possible therapeutic approaches for patients.
RESUMEN
Background: The development of HCC is often associated with extensive metabolic disturbances. Single cell RNA sequencing (scRNA-seq) provides a better understanding of cellular behavior in the context of complex tumor microenvironments by analyzing individual cell populations. Methods: The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) data was employed to investigate the metabolic pathways in HCC. Principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) analysis were applied to identify six cell subpopulations, namely, T/NK cells, hepatocytes, macrophages, endothelial cells, fibroblasts, and B cells. The gene set enrichment analysis (GSEA) was performed to explore the existence of pathway heterogeneity across different cell subpopulations. Univariate Cox analysis was used to screen genes differentially related to The Overall Survival in TCGA-LIHC patients based on scRNA-seq and bulk RNA-seq datasets, and LASSO analysis was used to select significant predictors for incorporation into multivariate Cox regression. Connectivity Map (CMap) was applied to analysis drug sensitivity of risk models and targeting of potential compounds in high risk groups. Results: Analysis of TCGA-LIHC survival data revealed the molecular markers associated with HCC prognosis, including MARCKSL1, SPP1, BSG, CCT3, LAGE3, KPNA2, SF3B4, GTPBP4, PON1, CFHR3, and CYP2C9. The RNA expression of 11 prognosis-related differentially expressed genes (DEGs) in normal human hepatocyte cell line MIHA and HCC cell lines HCC-LM3 and HepG2 were compared by qPCR. Higher KPNA2, LAGE3, SF3B4, CCT3 and GTPBP4 protein expression and lower CYP2C9 and PON1 protein expression in HCC tissues from Gene Expression Profiling Interactive Analysis (GEPIA) and Human Protein Atlas (HPA) databases. The results of target compound screening of risk model showed that mercaptopurine is a potential anti-HCC drug. Conclusion: The prognostic genes associated with glucose and lipid metabolic changes in a hepatocyte subpopulation and comparison of liver malignancy cells to normal liver cells may provide insight into the metabolic characteristics of HCC and the potential prognostic biomarkers of tumor-related genes and contribute to developing new treatment strategies for individuals.
RESUMEN
[This corrects the article DOI: 10.3389/fonc.2023.1104262.].