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.

Liver hepatocellular carcinoma (LIHC) is a prevalent form of primary liver cancer. Research has demonstrated the contribution of tumor stem cells in facilitating tumor recurrence, metastasis, and treatment resistance. Despite this, there remains a lack of established cancer stem cells (CSCs)-associated genes signatures for effectively predicting the prognosis and guiding the treatment strategies for patients diagnosed with LIHC. The single-cell RNA sequencing (scRNA-seq) and bulk RNA transcriptome data were obtained based on public datasets and computerized firstly using CytoTRACE package and One Class Linear Regression (OCLR) algorithm to evaluate stemness level, respectively. Then, we explored the association of stemness indicators (CytoTRACE score and stemness index, mRNAsi) with survival outcomes and clinical characteristics by combining clinical information and survival analyses. Subsequently, weighted co-expression network analysis (WGCNA) and Cox were applied to assess mRNAsi-related genes in bulk LIHC data and construct a prognostic model for LIHC patients. Single-sample gene-set enrichment analysis (ssGSEA), Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) and Tumor Immune Estimation Resource (TIMER) analysis were employed for immune infiltration assessment. Finally, the potential immunotherapeutic response was predicted by the Tumor Immune Dysfunction and Exclusion (TIDE), and the tumor mutation burden (TMB). Additionally, pRRophetic package was applied to evaluate the sensitivity of high and low-risk groups to common chemotherapeutic drugs. A total of four genes (including STIP1, H2AFZ, BRIX1, and TUBB) associated with stemness score (CytoTRACE score and mRNAsi) were identified and constructed a risk model that could predict prognosis in LIHC patients. It was observed that high stemness cells occurred predominantly in the late stages of LIHC and that poor overall survival in LIHC patients was also associated with high mRNAsi scores. In addition, pathway analysis confirmed the biological uniqueness of the two risk groups. Personalized treatment predictions suggest that patients with a low risk benefited more from immunotherapy, while those with a high risk group may be conducive to chemotherapeutic drugs. The current study developed a novel prognostic risk signature with genes related to CSCs, which provides novel ideas for the diagnosis, prognosis and treatment of LIHC.

Liu L ，Zhang M ，Cui N ，Liu W ，Di G ，Wang Y ，Xi X ，Li H ，Shen Z ，Gu M ，Wang Z ，Jiang S ，Liu B ... -  《PLoS One》

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

Cancer stem cells (CSCs) are characterized by self-renewal and -differential potential as compared to common cancer cells and play an important role in the development and therapeutic resistance of liver hepatocellular carcinoma (LIHC). However, the specific pathogenesis of LIHC stem cells is still unclear, and the genes involved in the stemness of LIHC stem cells are currently unknown. In this study, we investigated novel biomarkers associated with LIHC and explored the expression characteristics of stem cell-related genes in LIHC. We found that mRNA expression-based stemness index (mRNAsi) was significantly overexpressed in liver cancer tissues. Further, mRNAsi expression in LIHC increased with the tumor pathological grade, with grade 4 tumors harboring the greatest stem cell features. Upon establishing mRNAsi scores based on mRNA expression of every gene, we found an association with poor overall survival in LIHC. Moreover, modules of interest were determined based on weighted gene co-expression network analysis (WGCNA) inclusion criteria, and three significant modules (red, green, and brown) and 21 key genes (DCN, ECM1, HAND2, PTGIS, SFRP1, SRPX, COLEC10, GRP182, ADAMTS7, CD200, CDH11, COL8A1, FAP, LZTS1, MAP1B, NAV1, NOTCH3, OLFML2A, PRR16, TMEM119, and VCAN) were identified. Functional analysis of these 21 genes demonstrated their enrichment in pathways involved in angiogenesis, negative regulation of DNA-binding transcription factor activity, apoptosis, and autophagy. Causal relationship with proteins indicated that the Wnt, Notch, and Hypoxia pathways are closely related to LIHC tumorigenesis. To our knowledge, this is the first report of a novel CSC biomarker, mRNAsi, to predict the prognosis of LIHC. Further, we identified 21 key genes through mRNA expression network analysis, which could be potential therapeutic targets to inhibit the stemness of cancer cells in LIHC.

Bai KH ，He SY ，Shu LL ，Wang WD ，Lin SY ，Zhang QY ，Li L ，Cheng L ，Dai YJ ... -  《Cancer Medicine》

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

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.

Xu Q ，Xu H ，Chen S ，Huang W ... -  《Frontiers in Cell and Developmental Biology》

Predictive value of a stemness-based classifier for prognosis and immunotherapy response of hepatocellular carcinoma based on bioinformatics and machine-learning strategies.

Significant advancements have been made in hepatocellular carcinoma (HCC) therapeutics, such as immunotherapy for treating patients with HCC. However, there is a lack of reliable biomarkers for predicting the response of patients to therapy, which continues to be challenging. Cancer stem cells (CSCs) are involved in the oncogenesis, drug resistance, and invasion, as well as metastasis of HCC cells. Therefore, in this study, we aimed to create an mRNA expression-based stemness index (mRNAsi) model to predict the response of patients with HCC to immunotherapy. We retrieved gene expression and clinical data of patients with HCC from the GSE14520 dataset and the Cancer Genome Atlas (TCGA) database. Next, we used the "one-class logistic regression (OCLR)" algorithm to obtain the mRNAsi of patients with HCC. We performed "unsupervised consensus clustering" to classify patients with HCC based on the mRNAsi scores and stemness subtypes. The relationships between the mRNAsi model, clinicopathological features, and genetic profiles of patients were compared using various bioinformatic methods. We screened for differentially expressed genes to establish a stemness-based classifier for predicting the patient's prognosis. Next, we determined the effect of risk scores on the tumor immune microenvironment (TIME) and the response of patients to immune checkpoint blockade (ICB). Finally, we used qRT-PCR to investigate gene expression in patients with HCC. We screened CSC-related genes using various bioinformatics tools in patients from the TCGA-LIHC cohort. We constructed a stemness classifier based on a nine-gene (PPARGC1A, FTCD, CFHR3, MAGEA6, CXCL8, CABYR, EPO, HMMR, and UCK2) signature for predicting the patient's prognosis and response to ICBs. Further, the model was validated in an independent GSE14520 dataset and performed well. Our model could predict the status of TIME, immunogenomic expressions, congenic pathway, and response to chemotherapy drugs. Furthermore, a significant increase in the proportion of infiltrating macrophages, Treg cells, and immune checkpoints was observed in patients in the high-risk group. In addition, tumor cells in patients with high mRNAsi scores could escape immune surveillance. Finally, we observed that the constructed model had a good expression in the clinical samples. The HCC tumor size and UCK2 genes expression were significantly alleviated and decreased, respectively, by treatments of anti-PD1 antibody. We also found knockdown UCK2 changed expressions of immune genes in HCC cell lines. The novel stemness-related model could predict the prognosis of patients and aid in creating personalized immuno- and targeted therapy for patients in HCC.

Chen E ，Zou Z ，Wang R ，Liu J ，Peng Z ，Gan Z ，Lin Z ，Liu J ... -  《Frontiers in Immunology》

Constructing immune and prognostic features associated with ADCP in hepatocellular carcinoma and pan-cancer based on scRNA-seq and bulk RNA-seq.

Zhang Z ，Li Y ，Quan Z ，Li Y ，Zhu L ，Sun S ，Chen X ... -  《Frontiers in Immunology》