Dataset Information

RNA N6-Methyladenosine Patterns in Hepatocellular Carcinoma Reveal a Distinct Immune Infiltration Landscape and Clinical Significance

ABSTRACT:

Background

RNA N6-methyladenosine (m6A) methylation, the most abundant and prominent form of epigenetic modification, is involved in hepatocellular carcinoma (HCC) initiation and progression. However, the role of m6A methylation in HCC tumor microenvironment (TME) formation is unexplored. This study aimed to reveal the TME features of HCC patients with distinct m6A expression patterns and establish a prognostic model based on m6A signatures for HCC cohorts.

Material/Methods

We classified the m6A methylation patterns in 365 HCC samples based on 21 m6A modulators using a consensus clustering algorithm. Single-sample gene set enrichment analysis algorithm was used to quantify the abundance of immune cell infiltration. Gene set variation analysis revealed the biological characteristics between the m6A modification patterns. The m6A-based prognostic model was constructed using a training set with least absolute shrinkage and selection operator regression and validated in internal and external datasets.

Results

Two distinct m6A modification patterns exhibiting different TME immune-infiltrating characteristics, heterogeneity, and prognostic variations were identified in the HCC cohort. After depicting the immune landscape of TME in HCC, we found patients with high LRPPRC m6A modulator expression had depletion of T cells, cytotoxic cells, dendritic cells, and cytolytic activity response. A high m6A score, characterized by suppression of immunity, indicated an immune-excluded TME phenotype, with poor survival. A nomogram was developed to facilitate HCC clinical decision making.

Conclusions

Our results highlight the nonnegligible role of m6A methylation in TME formation and reveal a potential clinical application of the m6A-associated prognostic model for patients with HCC.

SUBMITTER: Zhao H

PROVIDER: S-EPMC8555444 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundSeveral studies have revealed that N6-methyladenosine (m6A) regulation is involved in various biological processes and cancer progression. Nevertheless, the potential effects of m6A modifications in the tumor immune microenvironment (TIME) and on immune regulation in pancreatic adenocarcinoma (PAAD) remains unclear.MethodsA consensus clustering algorithm was used to identify different m6A modification patterns and construct an m6A-associated gene signature based on 23 m6A regulators in PAAD. The CIBERSORT and ssGSEA algorithms were used to estimate the components of the immune cells in each sample. The PCA algorithm was used to develop the m6Ascore system for the evaluation of m6A modification patterns in each sample.ResultsTwo m6A modification patterns with different biological properties and prognoses were identified in 176 PAAD patient samples. The features of TIME between the two patterns were similar, with two definite immune phenotypes: immune-inflamed and immune-excluded. Based on the m6A phenotype-associated signature genes, we constructed an m6Ascore system to investigate the m6A modification pattern of each sample, profile the dissection of physiological processes, immune infiltration, clinical prognosis, immunotherapy, and genetic variation. Patients with low m6Ascore scores had better clinical outcomes, enhanced immune infiltration, and lower expression of immunotherapeutic drug targets, such as CD274 and PDCD1LG2. Further research indicated that the m6Ascore and tumor mutation burden were significantly correlated, and patients with low m6Ascore had higher mutation rates in SMAD4 and TTN. Moreover, TNFRSF21 was significantly upregulated in PAAD tumor tissues and cell lines. Lower expression of TNFRSF21 had a prominent advantage in survival and was correlated with a low level of immune infiltration. PAAD samples with different TNFRSF21 expression levels showed significantly distinct sensitivities to chemotherapeutic agents.ConclusionsThis study revealed that m6A modification patterns could play an important role in the diversity and complexity of TIME, and the m6Ascore system could serve as an independent and powerful prognostic biomarker and is latently related to PAAD immunotherapies. Quantitative determination of m6A modification patterns in individual patients will be instrumental in mapping the TIME landscape and further optimizing precision immunotherapy.