Identification of Hepatocellular Carcinoma Subtypes Based on Global Gene Expression Profiling to Predict the...

Background: Hepatocellular carcinoma (HCC) is a highly heterogeneous tumor, and distinguishing its subtypes holds significant value for diagnosis, treatment, and the prognosis. Methods: Unsupervised clustering analysis was conducted to classify HCC subtypes. Subtype signature genes were identified using LASSO, SVM, and logistic regression. Survival-related genes were identified using Cox regression, and their expression and function were validated via qPCR and gene interference. GO, KEGG, GSVA, and GSEA were used to determine enriched signaling pathways. ESTIMATE and CIBERSORT were used to calculate the stromal score, tumor purity, and immune cell infiltration. TIDE was employed to predict the patient response to immunotherapy. Finally, drug sensitivity was analyzed using the oncoPredict algorithm. Results: Two HCC subtypes with different gene expression profiles were identified, where subtype S1 exhibited a significantly shorter survival time. A subtype scoring formula and a nomogram were constructed, both of which showed an excellent predictive performance. COL11A1 and ACTL8 were identified as survival-related genes among the signature genes, and the downregulation of COL11A1 could suppress the invasion and migration of HepG2 cells. Subtype S1 was characterized by the upregulation of pathways related to collagen and the extracellular matrix, as well as downregulation associated with the xenobiotic metabolic process and fatty acid degradation. Subtype S1 showed higher stromal scores, immune scores, and ESTIMATE scores and infiltration of macrophages M0 and plasma cells, as well as lower tumor purity and infiltration of NK cells (resting/activated) and resting mast cells. Subtype S2 was more likely to benefit from immunotherapy. Subtype S1 appeared to be more sensitive to BMS-754807, JQ1, and Axitinib, while subtype S2 was more sensitive to SB505124, Pevonedistat, and Tamoxifen. Conclusions: HCC patients can be classified into two subtypes based on their gene expression profiles, which exhibit distinctions in terms of signaling pathways, the immune microenvironment, and drug sensitivity.