Author:
ren Min,Fan Bei,Cao Guangcai,Zong Rongrong,Feng Liaoliao,Sun Huiru
Abstract
Abstract
Background
It is widely acknowledged that hypoxia and m6A/m5C/m1A RNA modifications promote the occurrence and development of tumors by regulating the tumor microenvironment. This study aimed to establish a novel liver cancer risk signature based on hypoxia and m6A/m5C/m1A modifications.
Methods
We collected data from The Cancer Genome Atlas (TCGA-LIHC), the National Omics Data Encyclopedia (NODE-HCC), the International Cancer Genome Consortium (ICGC), and the Gene Expression Omnibus (GEO) databases for our study (GSE59729, GSE41666). Using Cox regression and least absolute shrinkage and selection operator (LASSO) method, we developed a risk signature for liver cancer based on differentially expressed genes related to hypoxia and genes regulated by m6A/m5C/m1A modifications. We stratified patients into high- and low-risk groups and assessed differences between these groups in terms of gene mutations, copy number variations, pathway enrichment, stemness scores, immune infiltration, and predictive capabilities of the model for immunotherapy and chemotherapy efficacy.
Results
Our analysis revealed a significantly correlated between hypoxia and methylation as well as m6A/m5C/m1A RNA methylation. The three-gene prognosis signature (CEP55, DPH2, SMS) combining hypoxia and m6A/m5C/m1A regulated genes exhibited strong predictive performance in TCGA-LIHC, NODE-HCC, and ICGC-LIHC-JP cohorts. The low-risk group demonstrated a significantly better overall survival compared to the high-risk group (p < 0.0001 in TCGA, p = 0.0043 in NODE, p = 0.0015 in ICGC). The area under the curve (AUC) values for survival at 1, 2, and 3 years are all greater than 0.65 in the three cohorts. Univariate and Multivariate Cox regression analyses of the three datasets indicated that the signature could serve as an independent prognostic predictor (p < 0.001 in the three cohorts). The high-risk group exhibited more genome changes and higher homologous recombination deficiency scores and stemness scores. Analysis of immune infiltration and immune activation confirmed that the signature was associated with various immune microenvironment characteristics. Finally, patients in the high-risk group experienced a more favorable response to immunotherapy, and various common chemotherapy drugs.
Conclusion
Our prognostic signature which integrates hypoxia and m6A/m5C/m1A-regulated genes, provides valuable insights for clinical prediction and treatment guidance for liver cancer patients.
Funder
Natural Science Basic Research Project of Shaanxi Province
Publisher
Springer Science and Business Media LLC
Reference85 articles.
1. Cisneros-Garza LE, González-Huezo MS, Moctezuma-Velázquez C, Ladrón de Guevara-Cetina L, Vilatobá M, García-Juárez I, Alvarado-Reyes R, Álvarez-Treviño GA, Allende-Pérez S, Bornstein-Quevedo L, et al. The second Mexican consensus on hepatocellular carcinoma. Part I: Epidemiology and diagnosis. Rev Gastroenterol Mex (Engl Ed). 2022;87:216–34.
2. Yang JD, Heimbach JK. New advances in the diagnosis and management of hepatocellular carcinoma. BMJ. 2020;371:m3544.
3. Craig AJ, von Felden J, Garcia-Lezana T, Sarcognato S, Villanueva A. Tumour evolution in hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2020;17:139–52.
4. Li H, Li CW, Li X, Ding Q, Guo L, Liu S, Liu C, Lai CC, Hsu JM, Dong Q, et al. MET inhibitors promote liver tumor evasion of the immune response by stabilizing PDL1. Gastroenterology. 2019;156:1849–1861.e1813.
5. Al-Salama ZT, Syed YY, Scott LJ. Lenvatinib: a review in hepatocellular carcinoma. Drugs. 2019;79:665–74.