Overweight/obesity-related transcriptomic signature as a correlate of clinical outcome, immune microenvironment, and treatment response in hepatocellular carcinoma

The pandemic of overweight and obesity (quantified by body mass index (BMI) ≥ 25) has rapidly raised the patient number of non-alcoholic fatty hepatocellular carcinoma (HCC), and several clinical trials have shown that BMI is associated with the prognosis of HCC. However, whether overweight/obesity is an independent prognostic factor is arguable, and the role of overweight/obesity-related metabolisms in the progression of HCC is scarcely known. In the present study, clinical information, mRNA expression profile, and genomic data were downloaded from The Cancer Genome Atlas (TCGA) as a training cohort (TCGA-HCC) for the identification of overweight/obesity-related transcriptome. Machine learning and the Cox regression analysis were conducted for the construction of the overweight/obesity-associated gene (OAG) signature. The Kaplan-Meier curve, receiver operating characteristic (ROC) curve, and the Cox regression analysis were performed to assess the prognostic value of the OAG signature, which was further validated in two independent retrospective cohorts from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). Subsequently, functional enrichment, genomic profiling, and tumor microenvironment (TME) evaluation were utilized to characterize biological activities associated with the OAG signature. GSE109211 and GSE104580 were retrieved to evaluate the underlying response of sorafenib and transcatheter arterial chemoembolization (TACE) treatment, respectively. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic response. Overweight/obesity-associated transcriptome was mainly involved in metabolic processes and noticeably and markedly correlated with prognosis and TME of HCC. Afterward, a novel established OAG signature (including 17 genes, namely, , , , , , , , , , , , , , , , , and ) divided patients into high and low OAG score groups with distinct prognosis (median overall survival (OS): 24.87 83.51 months, p < 0.0001), and the values of area under ROC curve (AUC) in predicting 1-, 2-, 3-, and 4-year OS were 0.81, 0.80, 0.83, and 0.85, respectively. Moreover, the OAG score was independent of clinical features and also exhibited a good ability for prognosis prediction in the ICGC-LIHC-JP cohort and GSE54236 dataset. Expectedly, the OAG score was also highly correlated with metabolic processes, especially oxidative-related signaling pathways. Furthermore, abundant enrich