Dataset Information

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

ABSTRACT:

Background

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).

Methods

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.

Results

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.

Conclusion

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.

SUBMITTER: Guo W

PROVIDER: S-EPMC8149800 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Gastric cancer is the fourth most frequent cancer and has a high death rate. Immunotherapy represented by PD-1 has brought hope for the treatment of advanced gastric cancer. Methylation of the m6A genes is linked to the onset and progression of numerous cancers, but there are few studies on gastric cancer. The main purpose of this study aims to analyze the relationship between m6A RNA methylation regulators, PD-L1, prognosis and tumor immune microenvironment (TIME) in gastric cancer. The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEx) databases were used to acquire transcriptomic data and clinical information from gastric cancer patients. The changes in m6A regulator expression levels in gastric cancer tissues and normal tissues were studied. Consensus clustering analysis was used to separate gastric cancer samples into two categories. We employed Least Absolute Shrinkage, Selection Operator (LASSO) Cox regression analysis, Gene Set Enrichment Analysis (GSEA), and cBioPortal to analyze the m6A regulators, PD-L1 and TIME in gastric cancer. In gastric cancer tissues, the majority of m6A regulatory factors are considerably overexpressed. Two gastric cancer subgroups (Cluster1/2) based on consensus clustering of 21 m6A regulators. PD-L1 and PD-1 expression levels were significantly higher in gastric cancer tissues, and they were significantly linked with METTL3, WTAP, HNRNPD, ZC3H7B, METTL14, FTO, PCIF1, HNRNPC, YTHDF1 and YTDHF2. Cluster1 showed a large increase in resting memory CD4+ T cells, regulatory T cells, naïve B cells, active NK cells, and resting Mast cells. Cluster1 and Cluster2 were shown to be involved in numerous critical signaling pathways, including base excision repair, cell cycle, nucleotide excision repair, RNA degradation, and spliceosome pathways. Gastric cancer RiskScores based on prognostic factors have been found as independent prognostic indicators. The amount of tumor-infiltrating immune cells is dynamically affected by changes in the copy number of m6A methylation regulators associated with TIME.

Project description:Recently, N 6-methyladenosine (m6A) RNA methylation in eukaryotic mRNA has become increasingly obvious in the pathogenesis and prognosis of cancer. Moreover, tumor microenvironment is involved in the regulation of tumorigenesis. In our research, the clinical data, including 374 tumor and 50 normal patients, were obtained from The Cancer Genome Atlas (TCGA). Then 19 m6A regulators were selected from other studies. Hepatocellular carcinoma (HCC) patients were clustered in cluster1/2, according to the consensus clustering for the m6A RNA regulators. We found that m6A regulators were upregulated in cluster1. The cluster1 was associated with higher programmed death ligand 1 (PD-L1) expression level, higher immunoscore, worse prognosis, and distinct immune cell infiltration compared with cluster2. Five risk signatures were identified, including YTH N6-methyladenosine RNA-binding protein 1, YTHDF2, heterogeneous nuclear ribonucleoprotein C, WT1-associated protein, and methyltransferase-like 3, based on univariate Cox and least absolute shrinkage and selection operator regression analysis. High-risk group and low-risk group HCC patients were selected based on the risk score. Similarly, the high-risk group was extremely associated with higher PD-L1 expression level, higher grade, and worse overall survival (OS). Also, cluster1 was mainly enriched in high-risk group. Receiver operating characteristic (ROC) and a nomogram were used to predict the ability and the probability of 3- and 5-year OS of HCC patients. The time-dependent ROC curve (AUC) reached 0.77, 0.67, and 0.68 at 1, 3, and 5 years in the training dataset. Also, AUC areas of 1, 3, and 5 years were 0.7, 0.63, and 0.55 in the validation dataset. The gene set enrichment analysis showed that MTOR signaling pathway and WNT signaling pathway were correlated with cluster1 and high-risk group. Collectively, the research showed that the m6A regulators were significantly associated with tumor immune microenvironment in HCC. Risk characteristics based on m6A regulators may predict prognosis in patients with HCC and provide a new therapeutic target for improving the efficacy of immunotherapy.

Project description:Background: Hexokinase 2 not only plays a role in physiological function of human normal tissues and organs, but also plays a vital role in the process of glycolysis of tumor cells. However, there are few comprehensive studies on HK2 in esophageal carcinoma (ESCA) needs further study. Methods: Oncomine, Tumor Immune Estimation Resource (TIMER), The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database were used to analyze the expression differences of HK2 in Pan-cancer and ESCA cohort, and to analyze the correlation between HK2 expression level and clinicopathological features of TCGA ESCA samples. GO/KEGG, GGI, and PPI analysis of HK2 was performed using R software, LinkedOmics, GeneMANIA and STRING online tools. The correlation between HK2 and ESCA immune infiltration was analyzed TIMER and TCGA ESCA cohort. The correlation between HK2 expression level and m6A modification of ESCA was analyzed by utilizing TCGA ESCA cohort. Results: HK2 is highly expressed in a variety of tumors, and its high expression level in ESCA is closely related to the weight, cancer stages, tumor histology and tumor grade of ESCA. The analysis results of GO/KEGG showed that HK2 was closely related to cell adhesion molecule binding, cell-cell junction, ameboidal-type cell migration, insulin signaling pathway, hif-1 signaling pathway, and insulin resistance. GGI showed that HK2 associated genes were mainly involved in the glycolytic pathway. PPI showed that HK2 was closely related to HK1, GPI, and HK3, all of which played an important role in tumor proliferation. The analysis results of TIMER and TCGA ESCA cohort indicated that the HK2 expression level was related to the infiltration of various immune cells. TCGA ESCA cohort analyze indicated that the HK2 expression level was correlated with m6A modification genes. Conclusion: HK2 is associated with tumor immune infiltration and m6A modification of ESCA, and can be used as a potential biological target for diagnosis and therapy of ESCA.