Quantification of m6A RNA methylation modulators pattern was a potential biomarker for prognosis and associated with tumor immune microenvironment of pancreatic adenocarcinoma.

m6A is the most prevalent and abundant form of mRNA modifications and is closely related to tumor proliferation, differentiation, and tumorigenesis. In this study, we try to conduct an effective prediction model to investigated the function of m6A RNA methylation modulators in pancreatic adenocarcinoma and estimated the potential association between m6A RNA methylation modulators and tumor microenvironment infiltration for optimization of treatment. Expression of 28 m6A RNA methylation modulators and clinical data of patients with pancreatic adenocarcinoma and normal samples were obtained from TCGA and GTEx database. Differences in the expression of 28 m6A RNA methylation modulators between tumour (n = 40) and healthy (n = 167) samples were compared by Wilcoxon test. LASSO Cox regression was used to select m6A RNA methylation modulators to analyze the relationship between expression and clinical characteristics by univariate and multivariate regression. A risk score prognosis model was conducted based on the expression of select m6A RNA methylation modulators. Bioinformatics analysis was used to explore the association between the m6Ascore and the composition of infiltrating immune cells between high and low m6Ascore group by CIBERSORT algorithm. Evaluation of m6Ascore for immunotherapy was analyzed via the IPS and three immunotherapy cohort. Besides, the biological signaling pathways of the m6A RNA methylation modulators were examined by gene set enrichment analysis (GSEA). Expression of 28 m6A RNA methylation modulators were upregulated in patients with PAAD except for MTEEL3. An m6Ascore prognosis model was established, including KIAA1429, IGF2BP2, IGF2BP3, METTL3, EIF3H and LRPPRC was used to predict the prognosis of patients with PAAD, the high risk score was an independent prognostic indicator for pancreatic adenocarcinoma, and a high risk score presented a lower overall survival. In addition, m6Ascore was related with the immune cell infiltration of PAAD. Patients with a high m6Ascore had lower infiltration of Tregs and CD8+T cells but a higher resting CD4+ T infiltration. Patients with a low m6Ascore displayed a low abundance of PD-1, CTLA-4 and TIGIT, however, the IPS showed no difference between the two groups. The m6Ascore applied in three immunotherapy cohort (GSE78220, TCGA-SKCM, and IMvigor210) did not exhibit a good prediction for estimating the patients' response to immunotherapy, so it may need more researches to figure out whether the m6A modulator prognosis model would benefit the prediction of pancreatic patients' response to immunotherapy. Modulators involved in m6A RNA methylation were associated with the development of pancreatic cancer. An m6Ascore based on the expression of IGF2BP2, IGF2BP3, KIAA1429, METTL3, EIF3H and LRPPRC is proposed as an indicator of TME status and is instrumental in predicting the prognosis of pancreatic cancer patients.

Wang L ，Zhang S ，Li H ，Xu Y ，Wu Q ，Shen J ，Li T ，Xu Y ... -  《BMC CANCER》

Comprehensive Analysis of m6A RNA Methylation Regulators and the Immune Microenvironment to Aid Immunotherapy in Pancreatic Cancer.

Guo Y ，Wang R ，Li J ，Song Y ，Min J ，Zhao T ，Hua L ，Shi J ，Zhang C ，Ma P ，Yang C ，Zhu L ，Gan D ，Li S ，Liu X ，Su H ... -  《Frontiers in Immunology》

Comprehensive Analysis of PD-L1 Expression, Immune Infiltrates, and m6A RNA Methylation Regulators in Esophageal Squamous Cell Carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most common cancer types and represents a threat to global public health. N6-Methyladenosine (m6A) methylation plays a key role in the occurrence and development of many tumors, but there are still few studies investigating ESCC. This study attempts to construct a prognostic signature of ESCC based on m6A RNA methylation regulators and to explore the potential association of these regulators with the tumor immune microenvironment (TIME). The transcriptome sequencing data and clinical information of 20 m6A RNA methylation regulators in 453 patients with ESCC (The Cancer Genome Atlas [TCGA] cohort, n = 95; Gene Expression Omnibus [GEO] cohort, n = 358) were obtained. The differing expression levels of m6A regulators between ESCC and normal tissue were evaluated. Based on the expression of these regulators, consensus clustering was performed to investigate different ESCC clusters. PD-L1 expression, immune score, immune cell infiltration and potential mechanisms among different clusters were examined. LASSO Cox regression analysis was utilized to obtain a prognostic signature based on m6A RNA methylation modulators. The relationship between the risk score based on the prognostic signature and the TIME of ESCC patients was studied in detail. Six m6A regulators (METTL3, WTAP, IGF2BP3, YTHDF1, HNRNPA2B1 and HNRNPC) were observed to be significantly highly expressed in ESCC tissues. Two molecular subtypes (clusters 1/2) were determined by consensus clustering of 20 m6A modulators. The expression level of PD-L1 in ESCC tissues increased significantly and was significantly negatively correlated with the expression levels of YTHDF2, METL14 and KIAA1429. The immune score, CD8 T cells, resting mast cells, and regulatory T cells (Tregs) in cluster 2 were significantly increased. Gene set enrichment analysis (GSEA) shows that this cluster involves multiple hallmark pathways. We constructed a five-gene prognostic signature based on m6A RNA methylation, and the risk score based on the prognostic signature was determined to be an independent prognostic indicator of ESCC. More importantly, the prognostic value of the prognostic signature was verified using another independent cohort. m6A regulators are related to TIME, and their copy-number alterations will dynamically affect the number of tumor-infiltrating immune cells. Our study established a strong prognostic signature based on m6A RNA methylation regulators; this signature was able to accurately predict the prognosis of ESCC patients. The m6A methylation regulator may be a key mediator of PD-L1 expression and immune cell infiltration and may strongly affect the TIME of ESCC.

Guo W ，Tan F ，Huai Q ，Wang Z ，Shao F ，Zhang G ，Yang Z ，Li R ，Xue Q ，Gao S ，He J ... -  《Frontiers in Immunology》

The N6-Methyladenosine- (m6A-) Associated Genes Act as Strong Key Biomarkers for the Prognosis of Pancreatic Adenocarcinoma.

Pancreatic adenocarcinoma (PAAD) has become the major cause of cancer-related deaths globally. The m6A (N6-methyladenosine) alteration plays a crucial function in carcinogenesis and tumor progression. The role of genes related to m6A and their expression level in pancreatic cancer is not identified yet. The objective of this research analysis is a demonstration of the m6A RNA methylation regulators based as biomarkers for the PAAD diagnosis. About 23 extensively reported m6A RNA methylation regulators were identified through the Cancer Genome Atlas (TCGA) database. This identification was based on consensus clustering analysis, protein-protein integration (PPI) analysis, risk prognostic model, Cox-regression analysis, String Spearman analysis, and LASSO Cox-regression. Herein, we conclude that 23 m6A methylation regulators have a strong link with the clinical and molecular characteristics of PAAD. The three subgroups (1/2) of pancreatic adenocarcinoma were identified using the clustering of 23 m6A regulators. Subgroup cluster 2 had a lower survival rate than the subgroup of cluster 1, and the difference in grades between the two groups was substantial. An assessment was performed using the 23 reported m6A methylation regulators. Eight of these can be used as independent PAAD prognostic markers. The consequences of variable IGF2BP3 expression in PAAD were then investigated further. The key finding of this study was that the m6A methylation regulator gene has the main role in pancreatic tumors, and it may be used as a biomarker in the prognosis of the PAAD and for therapy purposes.

Li F ，Zhang P 《-》

Characterization of the m6A-Associated Tumor Immune Microenvironment in Prostate Cancer to Aid Immunotherapy.

The aim of this study was to elucidate the correlation between m6A modification and the tumor immune microenvironment (TIME) in prostate cancer (PCa) and to identify the m6A regulation patterns suitable for immune checkpoint inhibitors (ICIs) therapy. We evaluated the m6A regulation patterns of PCa based on 24 m6A regulators and correlated these modification patterns with TIME characteristics. Three distinct m6A regulation patterns were determined in PCa. The m6A regulators cluster with the best prognosis had significantly increased METTL14 and ZC3H13 expression and was characterized by low mutation rate, tumor heterogeneity, and neoantigens. The m6A regulators cluster with a poor prognosis had markedly high KIAA1429 and HNRNPA2B1 expression and was characterized by high intratumor heterogeneity and Th2 cell infiltration, while low Th17 cell infiltration and Macrophages M1/M2. The m6Ascore was constructed to quantify the m6A modification pattern of individual PCa patients based on m6A-associated genes. We found that the low-m6Ascore group with poor prognosis had a higher immunotherapeutic response rate than the high-m6Ascore group. The low-m6Ascore group was more likely to benefit from ICIs therapy. This study was determined that immunotherapy is more effective in low-m6Ascore PCa patients with poor prognosis.

Liu Z ，Zhong J ，Zeng J ，Duan X ，Lu J ，Sun X ，Liu Q ，Liang Y ，Lin Z ，Zhong W ，Wu W ，Cai C ，Zeng G ... -  《Frontiers in Immunology》