Comprehensive Analysis of the Prognostic Value and Immune Infiltrates of the Three-m5C Signature in Colon Carcinoma
Ontology highlight
ABSTRACT: Background
The 5-methylcytosine (m5C) is one of the important forms of RNA post modification, and its regulatory mechanism in tumors has received increasing attention. However, its potential role in colorectal cancer remains unclear. Materials and Methods
Here, we systematically investigated the genetic variation and prognostic value of the 14 m5c RNA methylation regulators in colon cancer. The prognostic risk score was constructed using three m5C regulators, which was verified in the GSE17536 (N=177), GSE41258 (N=248) and GSE38832 (N=122) datasets. Results
The risk score developed from the three-m5C signature represents an independent prognostic factor, which can accurately predict the prognosis of patients with colon cancer in multiple datasets. The cytokine–cytokine receptor interaction and chemokine signaling pathway were significantly enriched in the low-risk score group. Further analysis showed that the three-m5C signature was related to tumor immune microenvironment (TIME), affecting the abundance of tumor-infiltrating immune cells. Especially, patients with low risk score had higher immune score than those with high risk score. In addition, gene set enrichment analysis (GSEA) confirmed that all three regulatory factors are associated with the MAPK/p38 signaling pathway. Conclusion
In conclusion, our study illustrates that the three-m5C signature may be involved in the regulation of colon cancer immune microenvironment in synergy with the MAPK signaling pathway. Therefore, further studying the three-m5C signature regulatory mechanisms might provide promising targets for improving the responsiveness of colon cancer to immunotherapy.
SUBMITTER: Geng Q
PROVIDER: S-EPMC8542737 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:5-methylcytosine (m5C) modification is involved in tumor progression. However, the lncRNAs associated with m5C in lung squamous cell carcinoma (LUSC) have not been elucidated. The Cancer Genome Atlas database was used to get the open-accessed transcriptional profiling and clinical information of LUSC patients. All the statistical analyses were performed based on R software v 4.0.0 and SPSS13.0. First, there were 614 m5C-related lncRNAs identified under the criterion of |R|>0.4 and p < 0.001 with m5C genes. Next, a prognosis model based on ERICD, AL021068.1, LINC01341, AC254562.3, and AP002360.1 was established, which showed good prediction efficiency in both the training and validation cohorts. Next, a nomogram plot was established by combining the risk score and clinical features for a better application in clinical settings. Pathway enrichment analysis showed that the pathways of angiogenesis, TGF-β signaling, IL6-JAK-STAT3 signaling, protein secretion, androgen response, interferon-α response, and unfolded protein response were significantly enriched in the high-risk patients. Immune infiltration analysis showed that the risk score was positively correlated with neutrophils, resting CD4+ memory T cells, and M2 macrophages, yet negatively correlated with follicular helper T cells, CD8+ T cells, and activated NK cells. Moreover, we found that high-risk patients might be more sensitive to immunotherapy, imatinib, yet resistant to erlotinib, gefitinib, and vinorelbine. In summary, our prognosis model is an effective tool that could robustly predict LUSC patient prognosis, which had the potential for clinical guidance.
| S-EPMC9355133 | biostudies-literature
Project description:Background: N6-methyladenosine (m6A), 5-methylcytosine (m5C) and N1-methyladenosine (m1A) are the main RNA methylation modifications involved in the progression of cancer. However, it is still unclear whether m6A/m5C/m1A-related long non-coding RNAs (lncRNAs) affect the prognosis of head and neck squamous cell carcinoma (HNSCC). Methods: We summarized 52 m6A/m5C/m1A-related genes, downloaded 44 normal samples and 501 HNSCC tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas (TCGA) database, and then searched for m6A/m5C/m1A-related genes co-expressed lncRNAs. We adopt the least absolute shrinkage and selection operator (LASSO) Cox regression to obtain m6A/m5C/m1A-related lncRNAs to construct a prognostic signature of HNSCC. Results: This prognostic signature is based on six m6A/m5C/m1A-related lncRNAs (AL035587.1, AC009121.3, AF131215.5, FMR1-IT1, AC106820.5, PTOV1-AS2). It was found that the high-risk subgroup has worse overall survival (OS) than the low-risk subgroup. Moreover, the results showed that most immune checkpoint genes were significantly different between the two risk groups (p < 0.05). Immunity microenvironment analysis showed that the contents of NK cell resting, macrophages M2, and neutrophils in samples of low-risk group were significantly lower than those of high-risk group (p < 0.05), while the contents of B cells navie, plasma cells, and T cells regulatory (Tregs) were on the contrary (p < 0.05). In addition, patients with high tumor mutational burden (TMB) had the worse overall survival than those with low tumor mutational burden. Conclusion: Our study elucidated how m6A/m5C/m1A-related lncRNAs are related to the prognosis, immune microenvironment, and TMB of HNSCC. In the future, these m6A/m5C/m1A-related lncRNAs may become a new choice for immunotherapy of HNSCC.
| S-EPMC8670092 | biostudies-literature
Project description:BackgroundN6-methyladenosine (m6A) and 5-methylcytosine (m5C) are the main RNA methylation modifications involved in the oncogenesis of cancer. However, it remains obscure whether m6A/m5C-related long non-coding RNAs (lncRNAs) affect the development and progression of low grade gliomas (LGG).MethodsWe summarized 926 LGG tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. 105 normal brain samples with RNA-seq data from the Genotype Tissue Expression project were collected for control. We obtained a molecular classification cluster from the expression pattern of sreened lncRNAs. The least absolute shrinkage and selection operator Cox regression was employed to construct a m6A/m5C-related lncRNAs prognostic signature of LGG. In vitro experiments were employed to validate the biological functions of lncRNAs in our risk model.ResultsThe expression pattern of 14 sreened highly correlated lncRNAs could cluster samples into two groups, in which various clinicopathological features and the tumor immune microenvironment were significantly distinct. The survival time of cluster 1 was significantly reduced compared with cluster 2. This prognostic signature is based on 8 m6A/m5C-related lncRNAs (GDNF-AS1, HOXA-AS3, LINC00346, LINC00664, LINC00665, MIR155HG, NEAT1, RHPN1-AS1). Patients in the high-risk group harbored shorter survival times. Immunity microenvironment analysis showed B cells, CD4 + T cells, macrophages, and myeloid-derived DC cells were significantly increased in the high-risk group. Patients in high-risk group had the worse overall survival time regardless of followed TMZ therapy or radiotherapy. All observed results from the TCGA-LGG cohort could be validated in CGGA cohort. Afterwards, LINC00664 was found to promote cell viability, invasion and migration ability of glioma cells in vitro.ConclusionOur study elucidated a prognostic prediction model of LGG by 8 m6A/m5C methylated lncRNAs and a critical lncRNA regulation function involved in LGG progression. High-risk patients have shorter survival times and a pro-tumor immune microenvironment.
| S-EPMC10320943 | biostudies-literature
Project description:MEX3A is a critical RNA-binding ubiquitin ligase that is upregulated in various types of cancer. However, the correlations of MEX3A with prognosis and its molecular mechanism in ovarian cancer (OC) remain unclear. The expression level, prognostic values, and the genetic variations of MEX3A were analyzed via Gene Expression Profiling Interactive Analysis (GEPIA) Oncomine, Kaplan-Meier plotter, and cBioPortal. We used the LinkedOmics database to investigate the functions of MEX3A coexpressed genes and performed visualizing gene interaction network analysis on the GeneMANIA website. The correlations between MEX3A and cancer immune infiltration were analyzed by the Tumor Immune Estimation Resource (TIMER) site and the TISIDB database. Furthermore, in vitro analysis was performed to evaluate the biological functions of MEX3A in OC cells. Our study showed that the expression of the MEX3A in OC was higher than in normal tissues; it had the greatest prognostic value in OC, and strong physical interaction with PABPC1, LAMTOR2, KHDRBS2, and IGF2BP2, which indicated the association between MEX3A and immune infiltration. We also found that MEX3A was negatively related to infiltrating levels of several types of immune cells, including macrophages, neutrophils, dendritic cells (DCs), B cells, and CD8+ T cells. Additionally, in vitro experiments demonstrated that MEX3A promotes proliferation and migration in OC cells. Taken together, MEX3A might influence the biological functions of OC cells by regulating the immune infiltration in the microenvironment as a prognostic biomarker and a potential therapeutic target.
| S-EPMC8195639 | biostudies-literature
Project description:BackgroundSingle nucleotide polymorphism (SNP) of complement component 2 (C2) has been found to be significantly associated with hepatocellular carcinoma (HCC). However, little is known about the role and mechanism of C2 in HCC. In the present study, we aimed to explore the prognostic value of C2 and its correlation with tumor-infiltrating immune cells in HCC.Materials and methodsmRNA expression was downloaded from TCGA (365 HCC patients and 50 healthy controls), GSE14520 (220 HCC patients and 220 adjacent normal tissues), and ICGC HCC (232 HCC patients) cohorts. Unpaired Student's t-tests or ANOVA tests were used to evaluate differences of C2 expression. Univariate and multivariate analyses were used to analyze the prognostic value of C2. CIBERSORT was used to calculate the proportion of 22 kinds of tumor-infiltrating immune cells.ResultsSignificantly lower C2 expression was found at HCC compared to healthy controls, and C2 was associated with TNM stages. Higher C2 expression was significantly associated with better prognosis, and multivariate analysis showed that C2 was also an independent factor for the prognosis of HCC. Moreover, elevated CD4 T cells were found at HCC patients with higher C2 expression while the higher proportion of macrophage M0 cells was found in HCC patients with lower C2 expression. KEGG analysis showed that "cell cycle," "AMPK signaling pathway," and "PPAR signaling pathway" were enriched in HCC patients with higher C2 expression.ConclusionC2 is a prognostic factor for HCC and may be used as a therapeutic target for future treatment of HCC.
| S-EPMC7312969 | biostudies-literature
Project description:RNA modifications implicate pathological and prognosis significance in cancer development and progression, of which, m6A and m5C are representative regulators. These RNA modifications could produce effects on the function of other RNA by regulating gene expression. Thus, in this study, we aimed to explore the correlation between m6A/m5C regulators and early-stage lung adenocarcinoma (LUAD). Only the early-stage LUAD samples were included in this investigation, and the RNA-seq dataset of The Cancer Genome Atlas (TCGA) cohort was utilized to evaluate the expression of 37 m6A/m5C regulated genes. Based on the expression level of these 37 genes, early-stage LUAD patients were divided into 2 clusters, which were performed by consensus clustering, and the m6A/m5C subtypes had significantly different prognostic outcomes (p < 0.001). Cluster1, which has a better prognosis, was characterized by the C3 (inflammatory) immune subtype, low immune infiltration, chemokine expression, major histocompatibility complex (MHC) expression, and immune checkpoint molecule expression. Furthermore, compared with cluster1, cluster2 showed a T cell exhaustion state, characterized by a high expression of immune checkpoint genes, and immune cells, such as T cells, CD8+ T cells, cytotoxic lymphocytes, NK cells, and so on. In addition, patients in cluster2 were with high tumor mutational burden (TMB) and numerous significant mutated oncogene and tumor suppressor genes, such as WNT10B, ERBB4, SMARCA4, TP53, and CDKN2A (p < 0.001). A total of 19 genes were mostly related to the prognosis of LUAD and were upregulated in cluster2 (p < 0.05), showing a positive correlation with the mRNA expression of 37 m6A/m5C regulated genes. The predictive risk model was constructed using Cox and LASSO (least absolute shrinkage and selection operator) regression analysis. Finally, a seven-gene m6A/m5C risk model, comprising of METTL3, NPLOC4, RBM15, YTHDF1, IGF2BP1, NSUN3, and NSUN7, was constructed to stratify the prognosis of early-stage LUAD (p = 0.0049, AUC = 0.791). The high-risk score was associated with a poorer prognosis. This model was also validated using two additional GEO datasets: GSE72094 (p = 0.011, AUC = 0.736) and GSE50081 (p = 0.012, AUC = 0.628). In summary, it was established that the m6A/m5C-regulated genes performed a crosstalk function in the mRNA expression of early-stage LUAD. By interacting with other mRNA genes, m6A/m5C modification disturbs DNA replication and the tumor immune microenvironment (TIME). The seven-gene risk model may be a critical tool for the prognostic assessment of early-stage LUAD.
| S-EPMC10095361 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) is the most common primary liver malignancy with no curative treatments. Plexin D1 (PLXND1) is a cellular receptor whose functions have been explored in several human cancers; however, the critical roles of PLXND1 in HCC have rarely been probed. Therefore, the present study attempted to elucidate the expression pattern, prognostic significance, and potential roles of PLXND1 in HCC. We found that PLXND1 expression in HCC tissues was significantly higher compared with normal liver tissue from Gene Expression Profiling Interactive Analysis (GEPIA) and Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB) databases. This result was further validated by immunohistochemistry staining (IHC) using tissue microarrays, which contained 216 HCC cases collected from our hospital. Additionally, PLXND1 expression showed a significant correlation with several clinical characteristics, including tumor grade and tumor hemorrhage (TH). Moreover, TISIDB and GEPIA databases were used to investigate the roles of PLXND1 in tumor-immune system interactions in HCC. As an immunoinhibitor, transforming growth factor-beta (TGF-β1) displayed the greatest correlations with PLXND1 in HCC. Finally, Kaplan-Meier curves and Cox analysis were conducted to further examine the potential clinical value of PLXND1 in HCC. We described a subclassification of HCC based on PLXND1 and TGF-β1 expression, which could be used to predict clinical outcomes and patient prognosis. Taken together, the results of this study indicate that PLXND1 might be a promising prognostic biomarker and potential therapeutic target in HCC.
| S-EPMC7820679 | biostudies-literature
Project description:The prognosis and immunotherapy response rates are unfavorable in patients with oral squamous cell carcinoma (OSCC). The tumor microenvironment is associated with tumor prognosis and progression, and the underlying mechanisms remain unclear. We obtained differentially expressed immune-related genes from OSCC mRNA data in The Cancer Genome Atlas (TCGA) database. Overall survival-related risk signature was constructed by univariate Cox regression analysis and LASSO Cox regression analysis. The prognostic performance was validated with receiver operating characteristic (ROC) analysis and Kaplan-Meier survival curves in the TCGA and Gene Expression Omnibus (GEO) datasets. The risk score was confirmed to be an independent prognostic factor and a nomogram was built to quantify the risk of outcome for each patient. Furthermore, a negative correlation was observed between the risk score and the infiltration rate of immune cells, as well as the expression of immunostimulatory and immunosuppressive molecules. Functional enrichment analysis between different risk score subtypes detected multiple immune-related biological processes, metabolic pathways, and cancer-related pathways. Thus, the immune-related gene signature can predict overall survival and contribute to the personalized management of OSCC patients.
| S-EPMC8724436 | biostudies-literature
Project description:BackgroundPancreatic cancer (PC) has a high mortality and dismal prognosis, predicting to be the second most lethal malignancy. 5-Methylcytosine (m5C) and long noncoding RNAs (lncRNAs) are both crucial in the prognostic outcome and immunotherapeutic effect for PC patients. Therefore, we aimed to create an m5C-related lncRNA signature (m5C-LS) for PC patients' prognosis and treatment.MethodsClinicopathological information and RNAseq data were acquired from The Cancer Genome Atlas (TCGA) database. Pearson's correlation analysis was used to extract m5C-related lncRNAs in PC. Univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analyses were adopted to build an m5C-LS. Kaplan-Meier (K-M), principal component analysis (PCA), and nomogram were utilized to assess model accuracy. In addition, we explored the model's possible immunotherapeutic responses and drug sensitivity targets.ResultsThree m5C-related lncRNAs were finally established to construct the risk signature, which has a good and independent predictive ability for PC patients. Based on the m5C-LS, patients were classified into the low- and high-m5C-LS group, with the latter having a worse prognosis. Furthermore, the m5C-LS allowed us to better discriminate the immunotherapeutic responses of PC patients in different subgroups.ConclusionsOur study constructed an m5C-LS and established a nomogram model that accurately predicted the prognosis of PC patients, as well as provides promising immunotherapeutic strategies in the future.
| S-EPMC8863480 | biostudies-literature
Project description:The development, incidence, and metastasis of tumors are all intimately correlated with 5-methylcytosine (m5C). However, uncertainty surrounds the function of m5C in acute myeloid leukemia (AML). In this study, multicenter AML data were collected and analyzed comprehensively to grasp the gene expression level, clinicopathological characteristics, prognostic significance of m5C in AML and its relationship with the tumor microenvironment (TME). The m5C gene-mediated scoring system (m5CSS) was created using principal component analysis, and multiple cox regression analyses were utilized to determine the prognostic relevance of the m5C score. The investigation of the correlation among m5C, immune characteristics, clinical characteristics, immune infiltration level, as well as drug reaction at immune checkpoints, and immunotherapy efficacy confirmed that the change of the characteristics of immune cell infiltration and patient prognosis are linked with the m5C gene. Moreover, the m5CSS was employed to assess the pattern of m5C modification. Further analyses showed that the m5C score can served as a reliable indicator of AML prognosis. Crucially, the prognostic value of the m5C score was validated in terms of drug resistance and immunotherapy. This work reveals that AML diversity and the generation of complex TMEs are both impacted by m5C modifications. Therefore, understanding the m5C modification pattern will improve grasp of TME infiltration characteristics and assist exploring more efficient immunotherapeutic approaches.
| S-EPMC10564437 | biostudies-literature