Abstract:Background: Liver cancer stem cells, characterized by self-renewal and initiating cancer cells, were decisive drivers of progression and therapeutic resistance in hepatocellular carcinoma (HCC). However, a comprehensive understanding of HCC stemness has not been identified. Methods: RNA sequencing information, corresponding clinical annotation, and mutation data of HCC were downloaded from The Cancer Genome Atlas-LIHC project. Two stemness indices, mRNA expression-based stemness index (mRNAsi) and epigenetically regulated-mRNAsi, were used to comprehensively analyze HCC stemness. Estimation of Stromal and Immune Cells in Malignant Tumors using Expression Data and single-sample gene-set enrichment analysis algorithm were performed to characterize the context of tumor immune microenvironment (TIME). Next, differentially expressed gene (DEG) analysis and weighted gene co-expression network analysis (WGCNA) were employed to identify significant mRNAsi-related modules with hub genes. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment pathways were analyzed to functionally annotate these key genes. The least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to establish a prognostic signature. Kaplan–Meier survival curves and receiver operating characteristic (ROC) analysis were applied for prognostic value validation. Seven algorithms (XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT, and CIBERSORT-ABS) were utilized to draw the landscape of TIME. Finally, the mutation data were analyzed by employing “maftools” package. Results: mRNAsi was significantly elevated in HCC samples. mRNAsi escalated as tumor grade increased, with poor prognosis presenting the higher stemness index. The stemness-related (greenyellow) modules with 175 hub genes were screened based on DEGs and WGCNA. A prognostic signature was established using LASSO analysis of prognostic hub genes to classify samples into two risk subgroups, which exhibited good prognostic performance. Additionally, prognostic risk-clinical nomogram was drawn to estimate risk quantitatively. Moreover, risk score was significantly associated with contexture of TIME and immunotherapeutic targets. Finally, potential interaction between risk score with tumor mutation burden (TMB) was elucidated. Conclusion: This work comprehensively elucidated that stemness characteristics served as a crucial player in clinical outcome, complexity of TIME, and immunotherapeutic prediction from both mRNAsi and mRNA level. Quantitative identification of stemness characteristics in individual tumor will contribute into predicting clinical outcome, mapping landscape of TIME further optimizing precision immunotherapy.

Identification of a 9-gene signature related to stemness for prognosis in hepatocellular carcinoma basing on weighted gene coexpression network analysis

A Novel Nine Gene Signature Integrates Stemness Characteristics Associated with Prognosis in Hepatocellular Carcinoma

Identification and Validation of a New Gene Signature Predicting Prognosis of Hepatocellular Carcinoma Patients by Network Analysis of Stemness Indices.

Identification and validation of a novel nine-gene prognostic signature of stem cell characteristic in hepatocellular carcinoma

Identification of Prognostic Genes in Hepatocellular Carcinoma

Immunological Value of Prognostic Signature Based on Cancer Stem Cell Characteristics in Hepatocellular Carcinoma

Identification of cancer stem cell characteristics in liver hepatocellular carcinoma by WGCNA analysis of transcriptome stemness index

Development and Validation of a Novel Prognostic Signature Based on M6a/m5c/m1a-Related Genes in Hepatocellular Carcinoma

Identification and Validation of a Potential Stemness-Associated Biomarker in Hepatocellular Carcinoma

Integration of scRNA-seq and bulk RNA-seq constructs a stemness-related signature for predicting prognosis and immunotherapy responses in hepatocellular carcinoma

Identification of a Prognostic Transcriptome Signature for Hepatocellular Carcinoma with Lymph Node Metastasis

Integrated analysis revealing a novel stemness-metabolism-related gene signature for predicting prognosis and immunotherapy response in hepatocellular carcinoma

Identification of Novel Gene Signature Associated with Cell Glycolysis to Predict Survival in Hepatocellular Carcinoma Patients

Identification of a five-gene signature in association with overall survival for hepatocellular carcinoma

Identification of a glycolysis-related gene signature for predicting prognosis in patients with hepatocellular carcinoma

Integration of single-cell RNA-seq and bulk RNA-seq to construct liver hepatocellular carcinoma stem cell signatures to explore their impact on patient prognosis and treatment

Identification and validation of real hub genes in hepatocellular carcinoma based on weighted gene co-expression network analysis

A novel immune-related gene signature predicts the prognosis of hepatocellular carcinoma

Five metastasis-related mRNAs signature predicting the survival of patients with liver hepatocellular carcinoma

An Integrated Model Based on a Six-Gene Signature Predicts Overall Survival in Patients With Hepatocellular Carcinoma

Combined Signature of Nine Immune-Related Genes: a Novel Risk Score for Predicting Prognosis in Hepatocellular Carcinoma.