A signature based on m6A pattern and tumor microenvironment infiltration in clear cell renal cell carcinoma.

RNA N6-methyladenosine (m6A) has been found to have a critical impact on clear cell renal cell carcinoma (ccRCC) by affecting the tumor microenvironment (TME) and immune cell (IC) infiltration and is related to the treatment and survival rate of patients with ccRCC. However, the mechanism of m6A in TME and IC infiltration remained unclear. Nonnegative Matrix Factorization (NMF) clustering was performed on 650 ccRCC cases from the Cancer Genome Atlas (TCGA) and the Gene-Expression Omnibus (GEO) datasets. The immune infiltration was generated by the single-sample gene-set enrichment analysis (ssGSEA) algorithm. Survival analyses were performed using the Kaplan-Meier method, and the significance of the differences was determined using the log-rank test. The m6A score was constructed based on the expression of m6A regulators to quantify m6A modification. The package "survminer R" was employed to layer patients' low and high scores groups and predict the immunotherapy response. Three different patterns of m6A modification were established, and significant differences in TME and IC infiltration features were found in these three patterns. Survival analysis demonstrated that m6A cluster A and m6A gene cluster A experienced a longer survival time. Evaluation of m6A modification patterns in individual tumors was initiated by the m6A score. The low m6A score subtype was characterized by increased tumor mutation burden (TMB) and immune infiltration, whereas a high m6A score with a lack of immune cell infiltration showed significantly better overall survival. m6A score was also associated with the expression of programmed cell death protein 1 (PD-L1) and cytotoxic T lymphocyte antigen 4 (CTLA-4). Patients in the high m6A score group had high PD-L1 expression and low CTLA-4 expression. Significant differences in prognosis were identified among types of different TMB and m6A scores, where low TMB and high m6A score had longer survival time. This research indicated that m6A modification greatly affected TME and IC infiltration. Physicians can develop practical immunotherapy strategies for patients with ccRCC by evaluating m6A-associated genes.

Yang C ，Yu T ，Li Q ，Xie F ，Lin Q ... -  《American Journal of Translational Research》

N6-Methyladenosine Modification Patterns and Tumor Microenvironment Immune Characteristics Associated With Clinical Prognosis Analysis in Stomach Adenocarcinoma.

Background: N6-methyladenosine (m6A) modification is a part of epigenetic research that has gained increasing attention in recent years. m6A modification is widely involved in many biological behaviors of intracellular RNA by regulating mRNA, thus affecting disease progression and tumor occurrence. However, the effects of m6A modification on immune cell infiltration of the tumor microenvironment (TME) are uncertain in stomach adenocarcinoma (STAD). Methods: The Cancer Genome Map (TCGA) database was used to download transcriptome data, clinicopathological data, and survival data for m6A-regulated genes in 433 STAD tissues that meet the requirements of this study. GSE84437 data were obtained from the Gene Expression Omnibus (GEO) database. The correlation between 23 m6A regulated genes was analyzed using R software. Sample clustering analysis was carried out on the genes of the m6A regulatory factor, and survival analysis and differentiation comparison were made for patients in clustering grouping. Then, the Gene Set Enrichment Analysis (GSEA), the single-sample GSEA (ssGSEA), and other methods were conducted to assess the correlation among m6A modification patterns, TME cell infiltration characteristics, and immune infiltration markers. The m6A modification pattern of individual tumors was quantitatively evaluated using the m6A score scheme of the principal component analysis (PCA). Results: From the TCGA database, 94/433 (21.71%) samples were somatic cell mutations, and ZC3H13 mutations are the most common. Based on the consensus, matrix k-3 is an optimal clustering stability value to identify three different clusters. Three types of m6A methylation modification patterns were significantly different in immune infiltration. Thus, 1028 differentially expressed genes (DEGs) were identified. The survival analysis of the m6A score found that patients in the high m6A score group had a better prognosis than those in the low m6A score group. Further analysis of the survival curve combining tumor mutation burden (TMB) and m6A scores revealed that patients had a significantly lower prognosis in the low tumor mutant group and the low m6A score group (p = 0.003). The results showed that PD-L1 was significantly higher in the high m6A score group than in the low score group (p < 2.22e-16). The high-frequency microsatellite instability (MSI-H) subtype score was significantly different from the other two groups. Conclusions: This study systematically evaluated the modification patterns of 23 m6A regulatory factors in STAD. The m6A modification pattern may be a critical factor leading to inhibitory changes and heterogeneity in TME. This elucidated the TME infiltration characteristics in patients with STAD through the evaluation of the m6A modification pattern.

Meijing Z ，Tianhang L ，Biao Y 《Frontiers in Cell and Developmental Biology》

m6A Methylation Patterns and Tumor Microenvironment Infiltration Characterization in Clear-Cell Renal Cell Carcinoma.

Increasing evidence suggests the essential regulation of RNA N6-methyladenosine (m6A) modification in carcinogenesis and immune response. Nevertheless, the potential impacts of these modifications on the tumor microenvironment (TME) immune cell infiltration characteristics in clear-cell renal cell carcinoma (ccRCC) remain unclear. Utilizing a consensus clustering algorithm, we determined three m6A modification patterns and identified three m6A-related gene clusters among 569 ccRCC samples, which were associated with different biological functions and clinical outcomes. Thereafter, the m6A score was constructed using m6A-associated signature genes to accurately exploit the m6A modification patterns within individual tumors. The m6A score was further demonstrated to be noticeably related to ccRCC prognosis. In addition, the m6A score was found to be strongly correlated with tumor mutational burden (TMB), microsatellite instability, immune infiltration, immune checkpoint expression, and immunotherapy response, which was also validated in the pan-cancer analyses. Our findings thoroughly elucidated that m6A modification contributes to tumor microenvironment immune-infiltrating characteristics and prognosis in ccRCC. Assessing the m6A modification patterns of individual patients with ccRCC will offer novel insights into TME infiltration and help develop more effective treatment strategies.

Ma T ，Wang J ，Liu X ，Zhang W ，Meng L ，Zhang Y ... -  《Frontiers in Genetics》

The significance of m6A RNA methylation modification in prognosis and tumor microenvironment immune infiltration of cervical cancer.

Recent studies have highlighted that N6-methyladenosine (m6A) plays a significant role in tumorigenicity and progression. However, the mechanism of m6A modifications in the tumor microenvironment (TME) immune cell infiltration in cervical cancer (CC) remains unclear. Clinical and RNA sequencing data of 25 m6A RNA methylation regulators were acquired from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. LASSO Cox regression analysis was used to generate a prognostic risk signature. m6A modification patterns were identified based on the expression of 25 m6A regulators, and their correlation with TME immune cell-infiltrating characterization was analyzed. Principal component analysis was used to construct an m6A-scoring signature (m6A score) to evaluate the m6A modification patterns of individual CC samples and guide the selection of more effective immunotherapeutic strategies. Genetic and expression alterations of 25 m6A regulators were highly heterogeneous between CC and normal tissues. METTL14 and IGF2BP1 were selected to conduct the prognostic risk signature. Three m6A modification patterns were identified in 659 CC samples, which were associated with distinct clinical outcomes and biological pathways. The TME immune cell-infiltrating characterization of the three m6A modification patterns was highly consistent with 3 tumor immune phenotypes, including immune-excluded, immune-inflamed, and immune-desert phenotypes. Due to the heterogeneity of m6A modification patterns, an m6A scoring signature was established to evaluate the m6A modification patterns of individual CC samples. Univariate and multivariate Cox regression analyses revealed that the m6A score is a robust and independent prognostic biomarker for assessing the prognosis of CC patients. A low m6A score, characterized by higher somatic mutation and higher expression of proliferation-related and DNA repair-related genes, indicated poor overall survival. Activation of immune infiltration was exhibited by the high m6A score, which was likely to have a good response and clinical benefits to antiPD-1/L1 immunotherapy. This study highlights the prognostic value of 25 m6A regulators in CC. The m6A modification is related to immune regulation and the formation of TME heterogeneity and complexity. An m6A scoring signature to clarify the individual m6A modification pattern could enhance our understanding of TME immune cell-infiltrating characterization and guide immunotherapy strategies.

Guo Y ，Bai Y ，Wang L ，Xu Z ，Wang X ，Wang W ... -  《-》

Identification of m6A modification patterns and development of m6A-hypoxia prognostic signature to characterize tumor microenvironment in triple-negative breast cancer.

N6-methylation (m6A) modification of RNA has been found to have essential effects on aspects of the tumor microenvironment (TME) including hypoxia status and mobilization of immune cells. However, there are no studies to explore the combined effect of m6A modification and hypoxia on molecular heterogeneity and TME of triple-negative breast cancer (TNBC). We collected The Cancer Genome Atlas (TCGA-TNBC, N=139), the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC-TNBC, N=297), the GSE103091, GSE21653, and GSE135565 series from the Gene Expression Omnibus (GEO-TNBC, N=247), and FUSCCTNBC (N=245) for our study. The non-negative matrix factorization algorithm was used to cluster TNBC samples. Immune cell infiltration was analyzed by the CIBERSORT algorithm. The enrichment scores were calculated by single-sample gene set enrichment analysis(ssGSEA) to characterize TME in TNBC samples. Immunohistochemistry (IHC) and qRT-PCR were performed to detect the gene expression. Based on the expression of m6A-related genes, we identified three distinct m6A clusters (denoted A, B, and C) in TNBC samples. Comparing the TME characteristics among the three clusters, we observed that cluster C was strongly related to hypoxia status and immune suppression, whereas clusters A and B displayed more immune cell infiltration. Therefore, we combine m6A and hypoxia related genes to classify two m6A-hypoxia clusters of TNBC and screened six prognostic genes by LASSO-Cox regression to construct a m6A-hypoxia signature(MHPS), which divided TNBC samples into high- and low-risk groups. We identified different TME features, immune cell infiltration between the two groups, and a better immunotherapy response was observed in the low-risk group. A nomogram was constructed with tumor size, lymph node, and risk score to improve clinical application of MHPS. We identified distinct TME characteristics of TNBC based on three different m6A modification patterns. Then, we constructed a specific m6A-hypoxia signature for TNBC to evaluate risk and predict immunotherapy response of patients, to enable more accurate treatment in the future.

Shen X ，Zhong J ，He J ，Han J ，Chen N ... -  《Frontiers in Immunology》