A cellular senescence-related classifier based on a tumorigenesis- and immune infiltration-guided strategy can predict prognosis, immunotherapy response, and candidate drugs in hepatocellular carcinoma
Authors: Luo Yi , Liu Hao , Fu Hong , Ding Guo-Shan , Teng Fei 
Published: November 2022

Background

Cellular senescence plays an irreplaceable role in tumorigenesis, progression, and tumor microenvironment (TME) remodeling. However, to date, there is limited research delineating the landscape of cellular senescence in hepatocellular carcinoma (HCC), and an improved understanding on the interaction of tumor-associated cellular senescence with HCC prognosis, TME, and response to immunotherapy is warrant.

Methods

Tumorigenic and immune infiltration-associated senescence genes were determined by weighted gene co-expression network analysis (WGCNA) and the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm, and subsequently, a prognostic scoring model (named TIS) was constructed using multiple survival analysis algorithms to classify the senescence-related subtypes of HCC. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were conducted to identify the distinct hallmark pathways between high- and low-risk subtypes. Additionally, we carried out correlation analyses for TIS and clinical traits, senescence-associated secretory phenotype (SASP), immune infiltration and evasion, immune checkpoint factors, drug response, and immunotherapeutic efficacy. External experimental validation was conducted to delineate the association of CPEP3 (a TIS gene) with HCC phenotypes through assays of proliferation, colony formation, and invasion.

Results

A five-gene TIS, composed of NET1, ATP6V0B, MMP1, GTDC1, and CPEB3, was constructed and validated using TCGA and ICGC datasets, respectively, and showed a highly robust and plausible signature for overall survival (OS) prediction of HCC in both training and validation cohorts. Patients in the TIS-high group were accompanied by worse OS, activation of carcinogenetic pathways, infiltration of immunosuppressive cells, exclusion of effector killing cells, overexpression of immunomodulatory genes and SASP, and unsatisfied response to immunotherapy. In response to anticancer drugs, patients in the TIS-high group exhibited enhanced susceptibility to several conventional chemotherapeutic agents (5-fluorouracil, docetaxel, doxorubicin, gemcitabine, and etoposide), as well as several inhibitors of pathways involved in cellular senescence (cell-cycle inhibitors, bromodomain and extraterminal domain family (BET) inhibitors, PI3K-AKT pathway inhibitors, and multikinase inhibitors). Additionally, four putative drugs (palbociclib, JAK3 inhibitor VI, floxuridine, and lestaurtinib) were identified as potential compounds for patients in the TIS-high group. Notably, in vitro functional validation showed that CPEB3 knockdown boosted the phenotypes of proliferation, clonogenicity, and invasion in HCC cells, whereas CPEB3 overexpression attenuated these phenotypes.

Conclusions

Our study provides comprehensive clues demonstrating the role of novel TIS in predicting HCC prognosis, immunotherapeutic response, and candidate drugs. This work highlights the significance of tumorigenesis- and immune infiltration-related cellular senescence in cancer therapy.