Development and validation of multi-omic prognostic signature of anoikis-related genes in liver hepatocellular carcinoma.

Liver hepatocellular carcinoma (LIHC) is characterized by high morbidity, rapid progression and early metastasis. Although many efforts have been made to improve the prognosis of LIHC, the situation is still dismal. Inability to initiate anoikis process is closely associated with cancer proliferation and metastasis, affecting patients' prognosis. In this study, a corresponding gene signature was constructed to comprehensively assess the prognostic value of anoikis-related genes (ARGs) in LIHC. Using TCGA-LIHC dataset, the mRNA levels of the differentially expressed ARGs in LIHC and normal tissues were compared by Student t test. And prognostic ARGs were identified through Cox regression analysis. Prognostic signature was established and then externally verified by ICGC-LIRI-JP dataset and GES14520 dataset via LASSO Cox regression model. Potential functions and mechanisms of ARGs in LIHC were evaluated by functional enrichment analyses. And the immune infiltration status in prognostic signature was analyzed by ESTIMATE algorithm and ssGSEA algorithm. Furthermore, ARGs expression in LIHC tissues was validated via qRT-PCR and IHC staining from the HPA website. A total of 97 differentially expressed ARGs were detected in LIHC tissues. Functional enrichment analysis revealed these genes were mainly involved in MAP kinase activity, apoptotic signaling pathway, anoikis and PI3K-Akt signaling pathway. Afterward, the prognostic signature consisting of BSG, ETV4, EZH2, NQO1, PLK1, PBK, and SPP1 had a moderate to high predictive accuracy and served as an independent prognostic indicator for LIHC. The prognostic signature was also applicable to patients with distinct clinical parameters in subgroup survival analysis. And it could reflect the specific immune microenvironment in LIHC, which indicated high-risk group tended to profit from ICI treatment. Moreover, qRT-PCR and IHC staining showed increasing expression of BSG, ETV4, EZH2, NQO1, PLK1, PBK and SPP1in LIHC tissues, which were consistent to the results from TCGA database. The current study developed a novel prognostic signature comprising of 7 ARGs, which could stratify the risk and effectively predict the prognosis of LIHC patients. Furthermore, it also offered a potential indicator for immunotherapy of LIHC.

Ding D ，Wang D ，Qin Y 《-》

A novel signature incorporating lipid metabolism- and immune-related genes to predict the prognosis and immune landscape in hepatocellular carcinoma.

Liver hepatocellular carcinoma (LIHC) is a highly malignant tumor with high metastasis and recurrence rates. Due to the relation between lipid metabolism and the tumor immune microenvironment is constantly being elucidated, this work is carried out to produce a new prognostic gene signature that incorporates immune profiles and lipid metabolism of LIHC patients. We used the "DEseq2" R package and the "Venn" R package to identify differentially expressed genes related to lipid metabolism (LRDGs) in LIHC. Additionally, we performed unsupervised clustering of LIHC patients based on LRDGs to identify their subgroups and immuno-infiltration and Gene Ontology (GO) enrichment analysis on the subgroups. Next, we employed multivariate, LASSO and univariate Cox regression analyses to determine variables and to create a prognostic profile on the basis of immune- and lipid metabolism-related differential genes (IRDGs and LRDGs). We separated patients into low- and high-risk groups in accordance with the best cut-off value of risk score. We conducted Decision Curve Analysis (DCA), Receiver Operating Characteristic curve analysis as a function of time as well as Survival Analysis to evaluate this signature's prognostic value. We incorporated the clinical characteristics of patients into the risk model to obtain a nomogram prognostic model. GEO14520 and ICGC-LIRI JP datasets were employed to externally confirm the accuracy and robustness of signature. The gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied for investigating the underlying mechanisms. Immune infiltration analysis was implemented to examine the differences in immune between both risk groups. Single-cell RNA sequencing (scRNA-SEQ) was utilized to characterize the genes that were involved in the distribution of signature and expression characteristics of different LIHC cell types. The patients' sensitivity in both risk groups to commonly used chemotherapeutic agents and semi-inhibitory concentrations (IC50) of the drugs was assessed using the GDSC database. On the basis of the differentially expressed genes (DEGs) in the two groups, the CMAP database was adopted for the prediction of potential small-molecule compounds. Small-molecule compounds were molecularly docked with prognostic markers. Lastly, we investigated the prognostic gene expression levels in normal and LIHC tissues with immunohistochemistry (IHC) and quantitative reverse transcription polymerase chain reaction(qRT-PCR). We built and verified a prognostic signature with seven genes that incorporated immune profiles and lipid metabolism. Patients were classified as low- and high-risk groups depending on their prognostic profiles. The overall survival (OS) was markedly lower in the high-risk group as compared to low-risk group. Time-dependent ROC curves more precisely predicted patients' survival at 1, 3 and 5 years; the area under the ROC curve was 0.81 (1 year), 0.75 (3 years) and 0.77 (5 years). The DCA curves showed the value of the prognostic genes in this signature for clinical applications. We included the patients' clinical characteristics in the risk model for both multivariate and univariate Cox regression analyses, and the findings revealed that the risk model represents an independent factor that influences OS in LIHC patients. With immune analysis, GSVA and GSEA, we identified that there are remarkable differences between the two risk groups in immune pathways, lipid metabolism, tumor development, immune cell infiltration and immune microenvironment, response to immunotherapy, and sensitivity to chemotherapy. Moreover, those with higher risk scores presented greater sensitivity to the chemotherapeutic agents. Experiments in vitro further elucidated the roles of SPP1 and FLT3 in the LIHC immune microenvironment. Furthermore, four small-molecule drugs that could target LIHC were screened. In vitro qRT-PCR , IHC revealed that the SPP1,KIF18A expressions were raised in LIHC in tumor samples, whereas FLT3,SOCS2 showed the opposite trend. We developed and verified a new signature comprising immune- and lipid metabolism-associated markers and to assess the prognosis and the immune status of LIHC patients. This signature can be applied to survival prediction, individualized chemotherapy, and immunotherapeutic guidance for patients with liver cancer. This study also provides potential targeted therapeutics and novel ideas for the immune evasion and progression of LIHC.

Yang T ，Luo Y ，Liu J ，Liu F ，Ma Z ，Liu G ，Li H ，Wen J ，Chen C ，Zeng X ... -  《Frontiers in Oncology》

Identification of a Necroptosis-Related Prognostic Signature and Associated Regulatory Axis in Liver Hepatocellular Carcinoma.

Liver hepatocellular carcinoma (LIHC) ranks the sixth in global cancer incidence with poor prognosis. Necroptosis is a kind of regulated cell death and has been proved to be of significance in cancer occurrence and progression. However, few studies comprehensively discuss the potential applications of necroptosis-related genes (NRGs) in the prognostic evaluation and immunotherapy of LIHC. The prognostic signature in the present study was built up using LASSO Cox regression analysis. Integrated bioinformatics tools were utilized to explore the potential mRNA-miRNA-lncRNA regulatory axis in LIHC. Furthermore, qRT-PCR method was used to verify the EZH2 expression in LIHC tissues. Furthermore, prognostic performance of EZH2 in LIHC was assessed by Kaplan-Meier method. A total of 14 NRGs were differentially expressed in LIHC tissues. The overall genetic mutation status of these NRGs in LIHC was also shown. NRGs were significantly correlated with programmed necrotic cell death, as well as Toll-like receptor signaling pathway in GO and KEGG pathway analysis. Kaplan-Meier analysis revealed that ALDH2, EZH2, NDRG2, PGAM5, RIPK1, and TRAF2 were related to the prognosis. A prognostic signature was constructed by these six genes and showed medium to high accuracy in the prediction of LIHC patients' prognosis. Further analysis revealed that NRGs were correlated with pathological stage, immune infiltration, and drug resistance in LIHC. Moreover, we identified a potential lncRNA TUG1/miR-26b-5p/EZH2 regulatory axis in LIHC, which might affect the progression of LIHC. qRT-PCR suggested a higher mRNA level of EZH2 in LIHC tissues. And a poor overall survival rate was detected in LIHC patients with high EZH2 expression. Moreover, EZH2 expression and cancer stage were identified as the independent risk factors affecting LIHC patients' prognosis. In the present study, we conducted comprehensive bioinformatic analyses and built up a necroptosis-related prognostic signature containing four genes (ALDH2, EZH2, NDRG2, and PGAM5) for patients with LIHC, and this prognostic signature showed a medium to high predictive accuracy. And our study also identified a lncRNA TUG1/miR-26b-5p/EZH2 regulatory axis, which might be of great significance in LIHC progression. In addition, based on the data from our center, the result of qRT-PCR and survival analysis showed a higher mRNA level of EZH2 in LIHC tissues and an unfavorable prognosis in high EZH2 expression group, respectively.

He A ，Huang Z ，Wang J ，Lu H ，Zhang R ，Wu L ，Feng Q ... -  《-》

Identification of anoikis related subtypes and construction of prognostic model in hepatocellular carcinoma.

Fang W ，Chen Z 《-》

Construction of an anoikis‒related prognostic signature to predict immunotherapeutic response and prognosis in hepatocellular carcinoma.

Anoikis resistance is an important inducer of tumor metastasis. The role of anoikis-related genes (ARGs) in hepatocellular carcinoma (HCC) remains unclear. A list of ARGs was obtained and regression analysis was employed to assemble an anoikis-related prognostic signature and risk score formula from mRNA data and clinical prognostic data downloaded from The Cancer Genome Atlas database. Quantitative real-time PCR (qRT-PCR) was performed on clinical samples to validate the selected ARGs expressions. Kaplan‒Meier curves, ROC curves, and Cox regression analyses were used to demonstrated the prognostic value of this signature. Biological functional enrichment analysis and immune infiltration analysis were utilized to analyze the differences in potential biological functions, immune cell infiltration, immune functions, and immunotherapeutic responses. A prognostic signature based on 6 ARGs and corresponding prognostic nomogram were established. Our qRT-PCR results showed a higher expression of 6 ARGs in HCC tissues (p value < 0.05). Kaplan‒Meier curves, ROC curves, and Cox regression analyses demonstrated good prognostic value of the signature in HCC (p value < 0.05). Significant differences between the enriched biological functions and immune landscapes of patients in different risk groups were discovered (p value < 0.05). In addition, patients with higher risk scores possibly had poor therapeutic response to transhepatic arterial chemotherapy and embolization or sorafenib, but their responses to immunotherapy might be more effective. A successful anoikis-related prognostic signature and corresponding clinical nomogram were established, which might facilitate better predictions of prognosis and therapeutic responses for HCC patients.

Xiong C ，Pan G ，Wang H ，Meng G ，Yan L ，Li R ，Yan Y ，Yang Y ，Zhang X ，Yang C ，Dong Z ，Li T ... -  《-》