Background: N6-methyladenosine (m6A) is related to the progression of multiple cancers. However, the underlying influences of m6A-associated genes on the tumor immune microenvironment in hepatocellular carcinoma (HCC) remain poorly understood. Therefore, we sought to construct a survival prediction model using m6A-associated genes to clarify the molecular and immune characteristics of HCC. Methods: HCC case data were downloaded from The Cancer Genome Atlas (TCGA). Then, by applying consensus clustering, we identified two distinct HCC clusters. Next, four m6A-related genes were identified to construct a prognostic model, which we validated with Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) datasets. Additionally, the molecular and immune characteristics in different subgroups were analyzed. Results: m6A RNA methylation regulators were differentially expressed between HCC and normal samples and linked with immune checkpoint expression. Using consensus clustering, we divided HCC samples into two subtypes with distinct clinical features. Cluster 2 was associated with unfavorable prognosis, higher immune checkpoint expression and immune cell infiltration levels. In addition, the immune and carcinogenic signaling pathways were enriched in cluster 2. Furthermore, we constructed a risk model using four m6A-associated genes. Patients with different risk scores had distinct survival times, expression levels of immunotherapy biomarkers, TP53 mutation rates, and sensitivities to chemotherapy and targeted therapy. Similarly, the model exhibited an identical impact on overall survival in the validation cohorts. Conclusion: The constructed m6A-based signature may be promising as a biomarker for prognostics and to distinguish immune characteristics in HCC.
Keywords: M6A; hepatocellular carcinoma; immune microenvironment; prognosis; therapy.
Copyright © 2021 Du, Ma, Zhu, Liu, Jiao, Yuan and Wang.