Dataset Information

M6A Classification Combined With Tumor Microenvironment Immune Characteristics Analysis of Bladder Cancer.

ABSTRACT:

Background

Studies have shown that N6-methyl adenosine (m6A) plays an important role in cancer progression; however, the underlying mechanism of m6A modification in tumor microenvironment (TME) cell infiltration of bladder cancer remains unclear. This study aimed to investigate the role of m6A modification in TME cell infiltration of bladder cancer.

Methods

The RNA expression profile and clinical data of bladder cancer were obtained from The Cancer Genome Atlas and Gene Expression Omnibus. We assessed the m6A modification patterns of 664 bladder cancer samples based on 20 m6A regulators through unsupervised clustering analysis and systematically linked m6A modification patterns to TME cell infiltration characteristics. Gene ontology and gene set variation analyses were conducted to analyze the underlying mechanism based on the assessment of m6A methylation regulators. Principal component analysis was used to construct the m6A score to quantify m6A modification patterns of bladder cancer.

Results

The genetic and expression alterations in m6A regulators were highly heterogeneous between normal and bladder tissues. Three m6A modification patterns were identified. The cell infiltration characteristics were highly consistent with the three immune phenotypes, including immune rejection, immune inflammation, and immune desert. The biological functions of three m6A modification patterns were different. Cox regression analyses revealed that the m6A score was an independent signature with patient prognosis (HR = 1.198, 95% CI: 1.031-1.390). Patients with a low-m6A score were characterized by increased tumor mutation burden, PD-L1 expression, and poorer survival. Patients in the low-m6A score group also showed significant immune responses and clinical benefits in the CTLA-4 immunotherapy cohort (p =0.0069).

Conclusions

The m6A methylation modification was related to the formation of TME heterogeneity and complexity. Assessing the m6A modification pattern of individual bladder cancer will improve the understanding of TME infiltration characteristics.

SUBMITTER: Zhu H

PROVIDER: S-EPMC8479184 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The emergence of immunotherapy has provided an option of treatment methods for bladder cancer (BC). However, the beneficiaries of immunotherapy are still limited to small-scale patients, and immunotherapy-related adverse events often occur. It is a major challenge for clinical work to study the immune subtypes of BC and the molecular mechanism of immune escape, and identify the immune responders accurately. Here, we explore the immune molecular subtypes of bladder cancer and potential escape mechanisms. First, we screened the expression profiles of 303 differentially expressed immune-related genes in BC patients from the Cancer Genome Atlas (TCGA) database, and successfully identified 4 molecular subtypes of BC. By comparing the clinical characteristics, immune cells infiltration, the expression of checkpoint genes, human leukocyte antigen (HLA) genes, and gene mutation status of different subtypes, we identified different clinical and immunological characteristics of 4 subtypes. Among 4 subtypes, Cluster 2 met the general characteristics of immunotherapy responders and responded well to immunotherapy, while Cluster 4 had the highest expression of immune characteristics, and is similar to the immune environment of normal bladder tissue. Then, the weighted gene co-expression network analysis (WGCNA) of immune-related genes revealed that brown module was positively correlated with subtypes. Pathway enrichment analysis explored the major pathways associated with subtypes, which are also associated with immune escape mechanisms. Moreover, the decision tree model, which was constructed by the principle of random forest screening factors, was also validated in internal validation set and external validation set from the Gene Expression Omnibus (GEO) cohort (GSE133624), and could achieve accurate subtypes prediction for BC patients with high-throughput sequencing. Taken together, we explored the immune molecular subtypes and their mechanisms of BC, and these results may provide guidance for the development of new BC immunotherapy strategies.

Project description:BackgroundN6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approaches of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between m6A and lncRNAs enrolling in the shaping of TIME remains unclear.MethodsThe RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering method was implemented to classify GC patients into two clusters. Survival analysis, the infiltration level of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. A competing endogenous RNA (ceRNA) network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of pathways in which m6A-related lncRNAs enriched. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal lncRNAs, and risk model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the response to immune checkpoint inhibitors (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect expression patterns of the selected lncRNAs in the 35 tumor tissues and their paired adjacent normal tissues, and validated the prognostic value of risk model in our cohort (N = 35).ResultsThe expression profiles of 15 lncRNAs were included to cluster patients into 2 subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and lower mutation rates of the genes. Different immune cell infiltration patterns were also displayed between the two clusters. GSEA showed that cluster1 preferentially enriched in tumor hallmarks and tumor-related biological pathways. KEGG pathway analysis found that the target mRNAs which m6A-related lncRNAs regulated by sponging miRNAs mainly enriched in vascular smooth muscle contraction, cAMP signaling pathway and cGMP-PKG signaling pathway. Next, eight lncRNAs were selected by LASSO regression algorithm to construct risk model. Patients in the high-risk group had poor prognoses, which were consistent in our cohort. As for predicting responses to ICIs therapy, patients from high-risk group were found to have lower tumor mutation burden (TMB) scores and account for large proportion in the Microsatellite Instability-Low (MSI-L) subtype. Moreover, patients had distinct immunophenoscores in different risk groups.ConclusionOur study revealed that the interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to ICIs therapy.

Project description:Objectives: Bladder cancer (BLCA) is the most common malignant tumor in the urinary system, while the prognosis of muscle-invasive bladder cancer (MIBC) is poor. Cuproptosis might be a promising therapeutic approach to trigger tumor cell death. This study aimed to figure out the role of cuproptosis in BLCA and constructed a new cuproptosis scoring system to guide clinical diagnosis and individualize treatments. Methods: Consensus clustering was used to classify 490 patients with BLCA from TCGA and GEO cohorts. Survival outcomes and functional enrichment analyses were performed between the different subtypes. The cuproptosis scoring system was constructed by LASSO-Cox analysis. ESTIMATE, CIBERSORT, and ssGSEA were used to investigate the tumor microenvironment (TME). Drug sensitivity was evaluated with pRRophetic. An immunotherapy cohort was used to investigate the treatment response. The cuproptosis scoring system was verified in our own cohort with quantitative real-time PCR. Results: An overview of 12 cuproptosis genes (CuGs) in the TCGA database was depicted. Based on the mRNA expression profiles of CuGs, patients were classified into two cuproptosis molecular patterns. Based on the differential genes between the two cuproptosis patterns, the patients were classified into two cuproptosis gene clusters. There were distinct survival outcomes, signaling pathways, and TME between the two subtypes. A 7-gene cuproptosis scoring system was constructed. Patients with high cuproptosis scores showed worse OS and more immunosuppressing TME than those with low cuproptosis scores. The two cuproptosis score groups had distinct mutation profiles. Patients with high cuproptosis scores tended to be sensitive to chemotherapy drugs, but insensitive to immune checkpoint inhibitors (ICIs) treatment. Conclusion: This study depicted the landscape of cuproptosis in BLCA. We constructed a cuproptosis scoring system to predict the prognosis of BLCA patients. There were significant differences in survival outcomes, TME, mutation profiles, and drug sensitivities in high and low cuproptosis score patients. The cuproptosis scoring system could help oncologists comprehensively understand the tumor characteristic of BLCA and make individualized treatment strategies.

Project description:BackgroundNon-small-cell lung cancer (NSCLC) is a major component of lung cancer and is significantly correlated with poor prognosis. N6-methyladenosine (m6A) RNA methylation is closely related to the occurrence, progression, and prognosis of cancer. The potential biological functions and mechanisms of m6A RNA methylation in the immune microenvironment are still unclear.MethodsWe assessed m6A RNA methylation modification patterns in 1326 NSCLC patient samples based on 20 m6A regulators, linking these clusters to the tumor microenvironment and immune cell infiltration. The m6Ascore was created to quantify the m6A modification patterns of individual tumors. We then assessed the value of NSCLC patients in terms of clinical prognosis and immunotherapy response.ResultsAccording to different mRNA expression levels, two different m6A clusters were identified. m6A aggregation was significantly associated with clinical prognostic characteristics, the tumor microenvironment, and immune-related biological processes. Fifteen differential genes were screened based on these two m6A clusters, and to further investigate the mechanisms of action of these differential genes, they were subjected to unsupervised clustering analysis, which classified them into four different genomic isoforms. Prognostic analysis indicated that the survival advantage of the m6A gene cluster A modification mode was significantly prominent. We continued to construct the m6Ascore, which was used as a scoring tool to evaluate tumor typing, immunity, and prognosis. Patients with a low m6Ascore showed a significant survival advantage, and the group with a low m6Ascore had a better prognosis predicted by immunotherapy. The anti-PD-1/L1 immunotherapy cohort showed that a lower m6Ascore was associated with higher efficacy of immunotherapy.ConclusionsThe results suggest that m6A RNA methylation regulators make an important difference in the tumor immune microenvironment of patients with NSCLC. m6A gene characterization and the construction of the m6Ascore provide us with a richer understanding of m6A RNA methylation modification patterns in NSCLC patients and help to predict clinical prognosis and immunotherapeutic response.