Shao-Hua Ren1,2, Ya-Fei Qin1,2, Hong Qin1,2, Hong-da Wang1,2, Guang-Ming Li1,2, Yang-Lin Zhu1,2, Cheng-Lu Sun1,2, Bo Shao1,2, Jing-Yi Zhang1,2, Jing-Peng Hao3, Hao Wang1,2. 1. Department of General Surgery, Tianjin Medical University General Hospital, Tianjin, People's Republic of China. 2. Tianjin General Surgery Institute, Tianjin Medical University General Hospital, Tianjin, People's Republic of China. 3. Department of Anorectal Surgery, The Second Hospital of Tianjin Medical University, Tianjin, People's Republic of China.
Abstract
Background: The prognostic value of m6A-related genes in hepatocellular carcinoma (HCC) and its correlation with the immune microenvironment still requires further investigation. Methods: Consensus clustering by m6A related genes was used to classify 374 patients with HCC from The Cancer Genome Atlas (TCGA) database. Then we performed the least absolute shrinkage and selection operator (LASSO) to construct the m6A related genes model. The International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO) datasets were used to verify and evaluate the model. ESTIMATE, CIBERSORTx, the expression levels of immune checkpoint genes, and TIDE were used to investigate the tumor microenvironment (TME) and the response to immunotherapy. Gene set enrichment analyses (GSEA), tumor-associated macrophages (TAMs), and gene-drug sensitivity were also analyzed. Results: By expression value and regression coefficient of five m6A related genes, we constructed the risk score of each patient. The patients with a higher risk score had a considerably poorer prognosis in the primary and validated cohort. For further discussing TME and the response to immunotherapy, we divided the entire set into two groups based on the risk score. Our findings implied that the tumor-infiltrating lymphocytes (TILs) were proportional to the risk scores, which seemed to contradict that patients with higher scores had a poor prognosis. Further, we found that the high-risk group had higher expression of PD-L1, CTLA-4, and PDCD1, indicating immune dysfunction, which may be a fundamental reason for poor prognosis. This was further reinforced by the fact that the low-risk group responded better than the high-risk group to monotherapy and combination therapy. Conclusion: The m6A related risk score is a new independent prognostic factor that correlates with immunotherapy response. It can provide a new therapeutic strategy for improving individual immunotherapy in HCC.
Background: The prognostic value of m6A-related genes in hepatocellular carcinoma (HCC) and its correlation with the immune microenvironment still requires further investigation. Methods: Consensus clustering by m6A related genes was used to classify 374 patients with HCC from The Cancer Genome Atlas (TCGA) database. Then we performed the least absolute shrinkage and selection operator (LASSO) to construct the m6A related genes model. The International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO) datasets were used to verify and evaluate the model. ESTIMATE, CIBERSORTx, the expression levels of immune checkpoint genes, and TIDE were used to investigate the tumor microenvironment (TME) and the response to immunotherapy. Gene set enrichment analyses (GSEA), tumor-associated macrophages (TAMs), and gene-drug sensitivity were also analyzed. Results: By expression value and regression coefficient of five m6A related genes, we constructed the risk score of each patient. The patients with a higher risk score had a considerably poorer prognosis in the primary and validated cohort. For further discussing TME and the response to immunotherapy, we divided the entire set into two groups based on the risk score. Our findings implied that the tumor-infiltrating lymphocytes (TILs) were proportional to the risk scores, which seemed to contradict that patients with higher scores had a poor prognosis. Further, we found that the high-risk group had higher expression of PD-L1, CTLA-4, and PDCD1, indicating immune dysfunction, which may be a fundamental reason for poor prognosis. This was further reinforced by the fact that the low-risk group responded better than the high-risk group to monotherapy and combination therapy. Conclusion: The m6A related risk score is a new independent prognostic factor that correlates with immunotherapy response. It can provide a new therapeutic strategy for improving individual immunotherapy in HCC.
Authors: William C Reinhold; Margot Sunshine; Hongfang Liu; Sudhir Varma; Kurt W Kohn; Joel Morris; James Doroshow; Yves Pommier Journal: Cancer Res Date: 2012-07-15 Impact factor: 12.701
Authors: Li Zhong; Dan Liao; Meifang Zhang; Cuiling Zeng; Xinchun Li; Ruhua Zhang; Haiqing Ma; Tiebang Kang Journal: Cancer Lett Date: 2018-11-10 Impact factor: 8.679
Authors: Jose J G Marin; Marta R Romero; Elisa Herraez; Maitane Asensio; Sara Ortiz-Rivero; Anabel Sanchez-Martin; Luca Fabris; Oscar Briz Journal: Semin Liver Dis Date: 2021-09-20 Impact factor: 6.115
Authors: Cristina P Rodriguez; Qian Vicky Wu; Jenna Voutsinas; Jonathan R Fromm; Xiuyun Jiang; Venu G Pillarisetty; Sylvia M Lee; Rafael Santana-Davila; Bernardo Goulart; Christina S Baik; Laura Q M Chow; Keith Eaton; Renato Martins Journal: Clin Cancer Res Date: 2019-12-03 Impact factor: 12.531
Authors: Stephen Mok; Richard C Koya; Christopher Tsui; Jingying Xu; Lídia Robert; Lily Wu; Thomas Graeber; Brian L West; Gideon Bollag; Antoni Ribas Journal: Cancer Res Date: 2013-11-18 Impact factor: 12.701
Authors: Peng Jiang; Shengqing Gu; Deng Pan; Jingxin Fu; Avinash Sahu; Xihao Hu; Ziyi Li; Nicole Traugh; Xia Bu; Bo Li; Jun Liu; Gordon J Freeman; Myles A Brown; Kai W Wucherpfennig; X Shirley Liu Journal: Nat Med Date: 2018-08-20 Impact factor: 53.440
Authors: Aaron M Newman; Chloé B Steen; Chih Long Liu; Andrew J Gentles; Aadel A Chaudhuri; Florian Scherer; Michael S Khodadoust; Mohammad S Esfahani; Bogdan A Luca; David Steiner; Maximilian Diehn; Ash A Alizadeh Journal: Nat Biotechnol Date: 2019-05-06 Impact factor: 54.908