Stratification of Hepatocellular Carcinoma Using N⁶-Methyladenosine

Background: The N⁶-methyladenosine (m⁶A) modification of eukaryotic mRNA is the most prevalent of such epigenetic modifications and has recently been identified as a potential player in the pathogenesis and progression of hepatocellular carcinoma (HCC). With the increasing emergence of immunotherapy in the treatment of HCC, we have evaluated the potential of m⁶A-related genes in predicting overall survival and the therapeutic efficacy of immunotherapy in HCC patients. Methods: We employed transcriptomic data from TCGA-LIHC and GSE76427, comprising a total of 485 HCC patients, as the training set. Based on 23 recognized m⁶A regulators, we performed clustering analysis on HCC patients. The intersecting differentially expressed genes (DEGs) among subtypes were used in least absolute shrinkage and selection operator (LASSO) Cox and multivariate Cox regression analyses to construct the risk model. For the quantification of a risk model of HCC patients, a risk score was developed and correlated with clinical and immunological parameters. Furthermore, a single-cell transcriptomic atlas was used to analyze the relationship between model genes and immune cell subpopulations. Mechanistic studies included in vitro assays to validate the association between the m⁶A-related gene ANLN and the progression of HCC. Results: Internal (TCGA and GEO) and external validation (ICGC) suggested that an 8-gene risk score provides an accurate and stable prognostic assessment for HCC. Furthermore, the high-risk score, characterized by elevated TP53 mutation frequency, tumor mutation burden (TMB), and tumor stem cell characteristics indicated a poor prognosis. The prognostic signature was associated with immune cell infiltration in HCC. Those patients with a high-risk score had lower immune tolerance with a better prediction of the efficacy of immunotherapy. The risk model helps to assess and predict the response and prognosis of HCC patients to immune checkpoint inhibitors (ICIs). Additionally, single-cell RNA sequencing data revealed that the high-risk group had a higher proportion of T cells and fewer immunosuppressive T cells, potentially correlating with a better response to immunotherapy. Finally, in vitro experiments showed that ANLN, an m⁶A-related gene, promoted the proliferation and migration of HCC cells. Conclusions: In this study, we identified and validated an m⁶A gene signature consisting of eight genes that can be used to predict prognosis and immunotherapy efficacy in HCC patients.