Establishment and validation of a novel lysosome-related gene signature for predicting prognosis and immune landscape in hepatocellular carcinoma.

Recent studies have shown that lysosomes not only provide energy for tumor cell growth, but also participate in the occurrence and development of malignant tumors by regulating various ways of tumor cell death. However, the role of lysosome associated genes (LSAGs) in hepatocellular carcinoma (HCC) remains unclear. Transcriptome data and clinical data of HCC were downloaded from the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) databases. We identified differential expression of LSAGs by comparing tumor tissue with normal liver tissue. Subsequently, we used univariate COX analysis and least absolute shrinkage and selection operator (LASSO) COX regression to construct the prognostic feature of LSAGs. Kaplan-Meier survival curve and receiver operating characteristic curve were used to evaluate the predictive ability of LSAGs feature. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used for functional enrichment analysis of risk differential genes. The relationship between LSAGs score and tumor microenvironment and chemotherapy drug sensitivity was analyzed. Finally, the cellular communication of tumor cells with high and low expression of model LSAGs was explored. We identified sixteen prognostic associated LSAGs, four of which were selected to construct prognostic feature of LSAGs. Patients in the low LSAGs group had a better prognosis than those in the high LSAGs group. GO and KEGG analyses showed that risk differential genes were enriched in leukocyte migration, cytokine-cytokine receptor interaction and PI3K-Akt signaling pathway. The group with low LSAGs score had lower immune score. Patients in the high LSAGs group were more sensitive to drugs for chemotherapy. In addition, tumor cells with high expression of model LSAGs showed stronger association with immune cells through the interleukin-2 (IL2), fibroblast growth factor (FGF), adiponectin, and bone morphogenetic proteins (BMP) signaling pathways. We established a LSAGs signature that had the ability to predict clinical prognosis and immune landscape, proposing potential therapeutic targets for HCC.

Li H ，Li J ，Qu X ，Dai H ，Liu J ，Ma M ，Wang J ，Dong W ，Wang W ... -  《-》

Characterizing the key genes of COVID-19 that regulate tumor immune microenvironment and prognosis in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC), a highly heterogeneous malignant tumor associated with a poor prognosis, is a common cause of cancer-related deaths worldwide, with a limited survival benefit for patients despite ongoing therapeutic breakthroughs. Coronavirus disease 2019 (COVID-19), a severe infectious disease caused by severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), is a global pandemic and a serious threat to human health. The increased susceptibility to SARS-CoV-2 infection and a poor prognosis in patients with cancer necessitate the exploration of the potential link between the two. No studies have investigated the relationship of COVID-19 genes with the prognosis and tumor development in patients with HCC. We screened prognosis-related COVID-19 genes in HCC, performed molecular typing, developed a stable and reliable COVID-19 genes signature for predicting survival, characterized the immune microenvironment in HCC patients, and explored new molecular therapeutic targets. Datasets of HCC patients, including RNA sequencing data and clinical information, were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. Prognosis-related COVID-19 genes were identified by univariate Cox analysis. Molecular typing of HCC was performed using the consensus non-negative matrix factorization method (cNMF), followed by the analysis of survival, tumor microenvironment, and pathway enrichment for each subtype. Prognostic signatures were constructed using LASSO-Cox regression models, and receiver operating characteristic (ROC) curves were used to validate the predictive performance of the signature. The same approach was used for the test and external validation sets. Seven software packages were applied to determine the abundance of immune infiltration in HCC patients and investigate its relationship with the risk scores. Gene set enrichment analysis (GSEA) was used to explore the potential mechanisms by which the COVID-19 genes affect hepatocarcinogenesis and prognosis. Three types of machine learning methods were combined to identify the most critical genes in the signature and localize their expression at the single cell level. We identified 53 prognosis-related COVID-19 genes and classified HCC into two molecular subtypes (C1, C2) by using the NMF method. The prognosis of C2 was significantly better than that of C1, and the two subtypes differed remarkably in terms of the tumor immune microenvironment and biological functions. The 17 COVID-19 genes were screened using the LASSO regression method to develop a 17 COVID-19 genes signature, which demonstrated a good predictive performance for 1-, 2- and 3-year OS of patients with HCC. The risk score as an independent prognostic factor for HCC has better predictive accuracy than traditional clinical variables. Patients in the TCGA cohort were categorized by risk score into the high- and low-risk groups, with the high-risk group mainly enriched in the immune modulation-related pathways and the low-risk group mainly enriched in the metabolism-related pathways, suggesting that the COVID-19 genes may affect disease progression and prognosis by regulating the tumor immune microenvironment and metabolism in HCC. NOL10 was identified as the most critical gene in the signature and hypothesized to be a potential therapeutic target for HCC. Objectively, the COVID-19 genes signature developed in this study, as an independent prognostic factor in HCC patients, is closely associated with the prognosis and tumor immune microenvironment of HCC patients and indicates that they may regulate the development of HCC in multiple ways, providing us with new perspectives for understanding the molecular mechanisms of HCC and finding effective therapeutic targets.

Gao S ，Zhang L ，Wang H 《-》

Mitochondrial-Related Transcriptome Feature Correlates with Prognosis, Vascular Invasion, Tumor Microenvironment, and Treatment Response in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is the most common subtype of primary liver cancer, which was highly correlated with metabolic dysfunction. Nevertheless, the association between nuclear mitochondrial-related transcriptome and HCC remained unclear. A total of 147 nuclear mitochondrial-related genes (NMRGs) were downloaded from the MITOMAP: A Human Mitochondrial Genome Database. The training dataset was downloaded from The Cancer Genome Atlas (TCGA), while validation datasets were retrieved from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). The univariate and multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses were applied to construct a NMRG signature, and the value of area under receiver operating characteristic curve (AUC) was utilized to assess the signature and nomogram. Then, data from the Genomics of Drug Sensitivity in Cancer (GDSC) were used for the evaluation of chemotherapy response in HCC. Functional enrichment of differentially expressed genes (DEGs) between HCC and paired normal tissue samples demonstrated that mitochondrial dysfunction was significantly associated with HCC development. Survival analysis showed a total of 35 NMRGs were significantly correlated with overall survival (OS) of HCC, and the LASSO Cox regression analysis further identified a 25-NMRG signature and corresponding prognosis score based on their transcriptional profiling. HCC patients were divided into high- and low-risk groups according to the median prognosis score, and high-risk patients had significantly worse OS (median OS: 27.50 vs. 83.18 months, P < 0.0001). The AUC values for OS at 1, 3, and 5 years were 0.79, 0.77, and 0.77, respectively. The prognostic capacity of NMRG signature was verified in the GSE14520 dataset and ICGC-HCC cohort. Besides, the NMRG signature outperformed each NMRG and clinical features in prognosis prediction and could also differentiate whether patients presented with vascular invasions (VIs) or not. Subsequently, a prognostic nomogram (C-index: 0.753, 95% CI: 0.703~0.804) by the integration of age, tumor metastasis, and NMRG prognosis score was constructed with the AUC values for OS at 1, 3, and 5 years were 0.82, 0.81, and 0.82, respectively. Notably, significant enrichment of regulatory and follicular helper T cells in high-risk group indicated the potential treatment of immune checkpoint inhibitors for these patients. Interestingly, the NMRG signature could also identify the potential responders of sorafenib or transcatheter arterial chemoembolization (TACE) treatment. Additionally, HCC patients in high-risk group appeared to be more sensitive to cisplatin, vorinostat, and methotrexate, reversely, patients in low-risk group had significantly higher sensitivity to paclitaxel and bleomycin instead. In summary, the development of NMRG signature provided a more comprehensive understanding of mitochondrial dysfunction in HCC, helped predict prognosis and tumor microenvironment, and provided potential targeted therapies for HCC patients with different NMRG prognosis scores.

Wang Y ，Song F ，Zhang X ，Yang C ... -  《-》

A Novel Oxidative Phosphorylation-Associated Gene Signature for Prognosis Prediction in Patients with Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a common type of malignant tumor with high morbidity and mortality. The oxidative phosphorylation (OXPHOS) metabolic pathway produces adenosine triphosphate (ATP) by delivering electrons to transmembrane protein complexes in the mitochondria. This research was dedicated to identifying an OXPHOS-associated signature for the assessment of prognosis of HCC patients. A total of 371 HCC patients from the Cancer Genome Atlas (TCGA) and 231 HCC patients from the International Cancer Genome Consortium (ICGC) with RNA expression data and clinical data were employed as construction and validation cohorts, respectively. The least absolute shrinkage and selection operator (LASSO) Cox regression was applied to establish a multigene signature in the TCGA cohort, and the ICGC cohort was used for validation. The prognostic value of the risk signature was evaluated using univariate and multivariate Cox regression, Kaplan-Meier curves, and receiver operating characteristic (ROC) curves. The potential enrichment of biological functions was investigated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Meanwhile, we analyzed the correlation between the risk score and the tumor microenvironment (TME). A five-gene signature including ATP6V0B, ATP6V1C1, ATP6V1E1, TIMM9, and UQCRH was identified by LASSO Cox regression to classify patients into low- and high-risk groups. ROC curve analysis indicated that the five-gene signature is a prospective prognostic factor in HCC patients. Univariate and multivariate Cox regression analyses demonstrated that the risk score was an independent prognostic factor for overall survival (OS). Functional analysis showed that differentially expressed genes (DEGs) between the low- and high-risk groups were enriched in mitosis and the cell cycle pathway. In addition, the five-gene signature was associated with innate immune cell infiltration, immune subtypes, and tumor stemness. A novel OXPHOS-associated gene signature can be used for prognostic prediction for patients with HCC.

Chen W ，Yang Z ，Chen Y 《-》

Construction and validation of a novel lysosomal signature for hepatocellular carcinoma prognosis, diagnosis, and therapeutic decision-making.

Chen J ，Gao G ，He Y ，Zhang Y ，Wu H ，Dai P ，Zheng Q ，Huang H ，Weng J ，Zheng Y ，Huang Y ... -  《Scientific Reports》