Prognostic Role of M6A-Associated Immune Genes and Cluster-Related Tumor Microenvironment Analysis: A Multi-Omics Practice in Stomach Adenocarcinoma

Front Cell Dev Biol. 2022 Jun 24:10:935135. doi: 10.3389/fcell.2022.935135. eCollection 2022.

Abstract

N6-methylandrostenedione (m6A) methylation plays a very important role in the development of malignant tumors. The immune system is the key point in the progression of tumors, particularly in terms of tumor treatment and drug resistance. Tumor immunotherapy has now become a hot spot and a new approach for tumor treatment. However, as far as the stomach adenocarcinoma (STAD) is concerned, the in-depth research is still a gap in the m6A-associated immune markers. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases is extremely important for our research, where we obtained gene mutation, gene expression data and relevant clinical information of STAD patients. Firstly, the samples from GEO were used as external validation groups, while the TCGA samples were divided into a training group and an internal validation group randomly. Using the way of Single factor COX-LASSO- and multi-factor Cox to construct the prognostic model. Then, all samples were subjected to cluster analysis to generate high and low expression groups of immune gene. Meanwhile, we also collected the correlation between these types and tumor microenvironment. On this basis, a web version of the dynamic nomogram APP was developed. In addition, we performed microenvironmental correlation, copy number variation and mutation analyses for model genes. The prognostic model for STAD developed here demonstrated a very strong predictive ability. The results of cluster analysis manifested that the immune gene low expression group had lower survival rate and higher degree of immune infiltration. Therefore, the immune gene low expression group was associated with lower survival rates and a higher degree of immune infiltration. Gene set enrichment analysis suggested that the potential mechanism might be related to the activation of immunosuppressive functions and multiple signaling pathways. Correspondingly, the web version of the dynamic nomogram APP produced by the DynNom package has successfully achieved rapid and accurate calculation of patient survival rates. Finally, the multi-omics analysis of model genes further enriched the research content. Interference of RAB19 was confirmed to facilitate migration of STAD cells in vitro, while its overexpression inhibited these features. The prognostic model for STAD constructed in this study is accurate and efficient based on multi-omics analysis and experimental validation. Additionally, the results of the correlation analysis between the tumor microenvironment and m6Ascore are the basics of further exploration of the pathophysiological mechanism in STAD.

Keywords: M6A; bioinformatic analysis; prognosis; stomach adenocarcinoma; survival; tumor microenvironment.