N6-Methylandenosine-Related lncRNAs in Tumor Microenvironment Are Potential Prognostic Biomarkers in Colon Cancer.
Ontology highlight
ABSTRACT: Background
LncRNA dysregulation and the tumor microenvironment (TME) have been shown to play a vital role in the progression and prognosis of colon cancer (CC). We aim to reveal the potential molecular mechanism from the perspective of lncRNA in the TME and provide the candidate biomarkers for CC prognosis.Methods
ESTIMATE analysis was used to divide the CC patients into high and low immune or stromal score groups. The expression array of lncRNA was re-annotated by Seqmap. Microenvironment-associated lncRNAs were filtered through differential analysis. The m6A-associated lncRNAs were screened by Pearson correlation analysis. Lasso Cox regression analyses were performed to construct the m6A- and tumor microenvironment-related lncRNA prognostic model (m6A-TME-LM). Survival analysis was used to assess the prognostic efficacy of candidate lncRNAs. Enrichment analyses annotated the candidate genes' functions.Results
We obtained 25 common differentially expressed lncRNAs (DELs) associated with immune microenvironment and m6A-related genes for subsequent lasso analysis. Four out of these DELs were selected for the m6A-TME-LM. All the four lncRNAs were related to overall survival, and a test set testified the result. Further stratification analysis of the m6A-TME-LM retained its ability to predict OS for male and chemotherapy adjuvant patients and performed an excellent prognostic efficacy in the TNM stage III and IV subgroups. Network analysis also found the four lncRNAs mediated co-expression network was associated with tumor development.Conclusion
We constructed the m6A-TME-LM, which could provide a better prognostic prediction of CC.
SUBMITTER: Zhang H
PROVIDER: S-EPMC8231021 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play vital roles in the prognostic value and immune microenvironment of malignant tumors. Here, we constructed a m6A-related lncRNA signature in osteosarcoma samples from TCGA dataset and analyzed the association of the signature with tumor immune microenvironment. m6A-related lncRNAs were identified by performing Pearson's correlation analysis and were used to construct a novel m6A-related lncRNA signature in osteosarcoma. Validation in testing and entire cohorts confirmed the satisfactory accuracy of the risk signature. Principal-component analysis verifies the grouping ability of the risk signature. Functional enrichment analyses connected immune with the risk signature based on the six m6A-related lncRNAs. When patients were separated into high- and low-risk group based on their risk scores, we found that patients in the high-risk group had lower stromal scores, immune scores, and ESTIMATE scores, while the tumor purity was higher in the high-risk group than that in the low-risk group. As for immune cell infiltration, the proportion of monocytes was significantly higher in the low-risk group than that in the high-risk group. Of the six lncRNAs, AC004812.2 was a protective factor in osteosarcoma and low expression of AC004812.2 predicted worse overall survival. Overexpression of AC004812.2 inhibited 143B cell proliferation and increased the expression levels of IGF2BP1 and YTHDF1. In all, our m6A-related lncRNA signature was a potential prognostic biomarker and correlated with tumor immune microenvironment and immune cell infiltration, and AC004812.2 might be an important regulator of m6A modification and a promising therapeutic target in osteosarcoma.
| S-EPMC8847230 | biostudies-literature
Project description:The prognostic value of N6-methylandenosine-related long non-coding RNAs (m6A-related lncRNAs) was investigated in 646 lower-grade glioma (LGG) samples from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) datasets. We implemented Pearson correlation analysis to explore the m6A-related lncRNAs, and then univariate Cox regression analysis was performed to screen their prognostic roles in LGG patients. Twenty-four prognostic m6A-related lncRNAs were identified as prognostic lncRNAs and they were inputted in a least absolute shrinkage and selection operator (LASSO) Cox regression to establish a m6A-related lncRNA prognostic signature (m6A-LPS, including 9 m6A-related prognostic lncRNAs) in the TCGA dataset. Corresponding risk scores of patients were calculated and divided LGG patients into low- and high-risk subgroups by the median value of risk scores in each dataset. The m6A-LPS was validated in the CGGA dataset and it showed a robust prognostic ability in the stratification analysis. Principal component analysis showed that the low- and high-risk subgroups had distinct m6A status. Enrichment analysis indicated that malignancy-associated biological processes, pathways and hallmarks were more common in the high-risk subgroup. Moreover, we constructed a nomogram (based on m6A-LPS, age and World Health Organization grade) that had a strong ability to forecast the overall survival (OS) of the LGG patients in both datasets. We also establish a competing endogenous RNA (ceRNA) network based on seven of the twenty-four m6A-related lncRNAs. Besides, we also detected five m6A-related lncRNA expression levels in 22 clinical samples using quantitative real-time polymerase chain reaction assay.
| S-EPMC7390977 | biostudies-literature
Project description:N6-Adenosine methylation, yielding N6-methyladenosine (m6A), is a reversible epigenetic modification found in messenger RNAs and long non-coding RNAs (lncRNAs), which affects the fate of modified RNA molecules and is essential for the development and differentiation of immune cells in the tumor microenvironment (TME). Osteosarcoma (OS) is the most common primary bone tumor in children and adolescents, and is characterized by high mortality. Currently, the possible role of m6A modifications in the prognosis of OS is unclear. In the present study, we investigated the correlation between m6A-related lncRNA expression and the clinical outcomes of OS patients via a comprehensive analysis. Clinical and workflow-type data were obtained from the Genotype-Tissue Expression Program and The Cancer Genome Atlas. We examined the relationship between m6A modifications and lncRNA expression, conducted Kyoto Encyclopedia of Genes analysis and also gene set enrichment analysis (GSEA), implemented survival analysis to investigate the association of clinical survival data with the expression of m6A-related lncRNAs, and utilized Lasso regression to model the prognosis of OS. Furthermore, we performed immune correlation analysis and TME differential analysis to investigate the infiltration levels of immune cells and their relationship with clinical prognosis. LncRNA expression and m6A levels were closely associated in co-expression analysis. The expression of m6A-related lncRNAs was quite low in tumor tissues; this appeared to be a predicting factor of OS in a prognostic model, independent of other clinical features. The NOD-like receptor signaling pathway was the most significantly enriched pathway in GSEA. In tumor tissues, SPAG4 was overexpressed while ZBTB32 and DEPTOR were downregulated. Tissues in cluster 2 were highly infiltrated by plasma cells. Cluster 2 presented higher ESTIMATE scores and stromal scores, showing a lower tumor cell purity in the TME. In conclusion, m6A-related lncRNA expression is strongly associated with the occurrence and development of OS, and can be used to as a prognostic factor of OS. Moreover, m6A-related lncRNAs and infiltrating immune cells in the TME could serve as new therapeutic targets and prognostic biomarkers for OS.
| S-EPMC8790065 | biostudies-literature
Project description:HCC is one of the most common types of malignancies worldwide and the fourth-leading cause of cancer deaths. Thus, there is an urgent need to search for novel targeted therapies in HCC. 186 m6a-related lncRNAs were screened for subsequent analysis. Two distinct m6A modification clusters were identified to be associated with the overall prognosis in TCGA-LIHC based on the m6A-related lncRNAs profiling, followed by univariate Cox regression analysis. In addition, four m6A-related lncRNAs prognostic signatures were developed and validated that could predict the OS of HCC patients, followed by univariate Cox regression, LASSO regression, and multivariate Cox regression analysis. Moreover, four m6A-related lncRNAs were identified to be related to HCC prognosis. ESTIMATE was used to evaluate the stromal score, immune score, ESTIMATE score, and tumor purity of each HCC sample. ssGSEA was performed to identify the enrichment levels of 29 immune signatures in each sample. Finally, quantitative real-time polymerase chain reaction shown that KDM4A-AS1, BACE1-AS, and NRAV expressions were upregulated in HCC patients. We proved that our m6A-related lncRNAs signature had powerful and robust ability for predicting OS of different HCC subgroups.
| S-EPMC8854774 | biostudies-literature
Project description:Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive, fatal tumor. N6-methylandenosine (m6A) methylation is the major epigenetic modification of RNA including lncRNAs. The roles of m6A-related lncRNAs in PDAC have not been fully clarified. This study aims to assess gene signatures and prognostic value of m6A-related lncRNAs in PDAC. The Cancer Genome Atlas (TCGA) dataset and the International Cancer Genome Consortium (ICGC) dataset were explored to identify m6A-related lncRNAs. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were performed to construct the m6A-related lncRNAs prognostic riskscore (m6A-LPR) model to predict the overall survival (OS) in the TCGA training cohort. Kaplan-Meier curve with log-rank test and receiver operating characteristic (ROC) curve were used to evaluate the prognostic value of the m6A-LPR. Furthermore, the robustness of the m6A-LPR was further validated in the ICGC cohort. Tumor immunity was evaluated using ESTIMATE and CIBERSORT algorithms. A total of 262 m6A-related lncRNAs were identified in two datasets. In the TCGA training cohort, 28 prognostic m6A-related lncRNAs were identified and the m6A-LPR including four m6A-related lncRNAs was constructed. The m6A-LPR was able to identify high-risk patients with significantly poorer OS and accurately predict OS in both the TCGA training cohort and the ICGC validation cohort. Analysis of tumor immunity revealed that high-risk groups had remarkably lower stromal, immune, and ESTIMATE scores. Moreover, high-risk groups were associated with significantly higher levels of plasma B cells and resting NK cells infiltration, and lower levels of infiltrating resting memory CD4 T cells, monocytes, and resting mast cells. Our study proposed a robust m6A-related prognostic signature of lncRNAs for predicting OS in PDAC, which provides some clues for further studies focusing on the mechanism process underlying m6A modification of lncRNAs.
| S-EPMC8426490 | biostudies-literature
Project description:Both lncRNAs and the N6-methyladenosine (m6A) modification are key regulators of tumorigenesis and innate immunity. However, little is known about the m6A modification of lncRNAs and their clinical and immune relevance in bladder cancer. In this study, we identified m6A-related lncRNAs using Pearson correlation analysis in The Cancer Genome Atlas (TCGA) and the IMvigor210 datasets. Next, univariate Cox regression was performed using the TCGA dataset to filter prognostic m6A-related lncRNAs, which were further subjected to the least absolute shrinkage and selection operator (LASSO) Cox regression to establish a 12 m6A-related lncRNA prognostic score (m6A-LRS). The m6A-LRS was validated in the IMvigor210 dataset. In addition, high m6A-LRS tumors, characterized by decreased tumor mutation load and neoantigen load, showed poorer response to immunotherapy than those with low m6A-LRS in the IMvigor210 dataset. Further, we constructed an m6A-LRS-based nomogram that demonstrated a strong ability to predict overall survival in patients with bladder cancer. Moreover, enrichment analysis revealed that tumor-associated biological processes, oncogenic signaling, and tumor hallmarks were commonly associated with a high m6A-LRS. Gene set variation analysis also indicated that high m6A-LRS was associated with activation of canonical oncogenic signatures, such as the epithelial-to-mesenchymal transition, cell cycle regulators, and DNA replication, as well as activation of immunosuppressive signatures, such as the T-cell exhaustion and pan-fibroblast-TGF-β response signatures. Furthermore, we observed distinct tumor microenvironment cell infiltration characteristics between high- and low-risk tumors. High m6A-LRS tumors showed reduced infiltration of CD8+ T-cells and enhanced infiltration of macrophages and fibroblasts. Additionally, we established a competing endogenous RNA network based on the12 m6A-related lncRNAs. Finally, three lncRNAs (SNHG16, SBF2-AS1, and BDNF-AS) were selected for further validation. The qualitative PCR assay on 10 pairs of bladder cancer and adjacent normal control samples validated the differential expression, and methylated RNA immunoprecipitation (MeRIP) analysis demonstrated a robust m6A enrichment in T24 bladder cancer cells compared with normal uroepithelial cells (SVHUC-1). In conclusion, this study introduced an m6A-related lncRNA signature that identified a subgroup of patients with poor prognoses and suboptimal immune responses, thus providing novel approaches for treatment response prediction and patient stratification in bladder cancer.
| S-EPMC8387102 | biostudies-literature
Project description:BackgroundColon adenocarcinoma (COAD) is the most common type of colon cancer. To date, however, the prognostic values of m6A RNA methylation-related long non-coding RNAs (lncRNAs) in COAD are largely unknown.Materials and methodsThe m6A-related lncRNAs were identified from The Cancer Genome Atlas (TCGA) data set. Univariate and multivariate Cox regression analyses were performed to explore the prognostic m6A-related lncRNAs. Consistent clustering analysis was performed to classify the COAD patients into different subgroups based on the expression of m6A-related lncRNAs. The potential biological functions as well as differences in the stemness index and tumor immune microenvironment between different subgroups were analyzed. The prognostic m6A-related lncRNAs were used to establish an m6A-related lncRNA risk model to predict prognosis and survival status.ResultsWe identified 31 m6A-associated lncRNAs with prognostic values from the TCGA data set. Based on the expression of prognostic m6A-associated lncRNAs, TCGA-COAD patients were classified into three clusters using consistent clustering analysis. There was a low correlation of tumor stemness between the three clusters but a significant correlation with the tumor immune microenvironment as well as the tumor mutational load. Thirty-one prognostic-related m6A-associated lncRNAs were used to construct a risk model, which was further determined by survival analysis, receiver operating characteristic (ROC) curve, and univariate and multifactor Cox analysis. The m6A-related risk model demonstrates good performance in predicting prognosis and survival status. The model-based high-risk group exhibited poorer overall survival (OS) compared with the low-risk group.ConclusionIn this study, we construct a risk model that consists of 31 m6A-related lncRNAs with independent prognostic values in COAD. Our study shows the critical roles of these 31 m6A-related lncRNAs in the tumor immune microenvironment, indicating the prospect of informing prognostic stratification and the development of immunotherapeutic strategies for COAD patients.
| S-EPMC8321625 | biostudies-literature
Project description:N6-methyladenosine (m6A) is one of the most prevalent RNA modifications in mRNA and non-coding RNA. In this study, we identified 10 upregulated m6A regulators at both mRNA and protein levels, and 2,479 m6A-related lncRNAs. Moreover, the m6A-related long noncoding RNAs (lncRNAs) could clearly stratify the colon adenocarcinoma (COAD) samples into three subtypes. The subtype 2 had nearly 40% of samples with microsatellite instability (MSI), significantly higher than the two other subtypes. In accordance with this finding, the inflammatory response-related pathways were highly activated in this subtype. The subtype-3 had a shorter overall survival and a higher proportion of patients with advanced stage than subtypes 1 and 2 (p-value < 0.05). Pathway analysis suggested that the energy metabolism-related pathways might be aberrantly activated in subtype 3. In addition, we observed that most of the m6A readers and m6A-related lncRNAs were upregulated in subtype 3, suggesting that the m6A readers and the m6A-related lncRNAs might be associated with metabolic reprogramming and unfavorable outcome in COAD. Among the m6A-related lncRNAs in subtype 3, four were predicted as prognostically relevant. Functional inference suggested that CTD-3184A7.4, RP11-458F8.4, and RP11-108L7.15 were positively correlated with the energy metabolism-related pathways, further suggesting that these lncRNAs might be involved in energy metabolism-related pathways. In summary, we conducted a systematic data analysis to identify the key m6A regulators and m6A-related lncRNAs, and evaluated their clinical and functional importance in COAD, which may provide important evidences for further m6A-related researches.
| S-EPMC8484874 | biostudies-literature
Project description:To investigate the role of N6-methyladenosine (m6A)- related long non-coding RNAs (lncRNAs) in bladder cancer (BC). 50 m6A-related lncRNAs were screened out and were correlated with prognosis from BC samples in The Cancer Genome Atlas (TCGA). The lncRNAs were subdivided into cluster 1 and cluster 2 with consensus cluster analysis, and it was found that lncRNAs in cluster 2 were associated with poor prognosis and increased PD-L1 expression. Gene set enrichment analysis (GSEA) revealed tumor-related pathways in cluster 2. Through least absolute shrinkage and selection operator (LASSO) Cox regression analysis, univariate and multivariate Cox regression, and ROC analyses, 14 prognostic lncRNAs were selected and used to construct the m6A-related lncRNA prognostic signature (m6A-LPS), furthermore, that m6A-LPS was as a valuable independent prognostic factor. Interestingly, the m6A-LPS risk score was positively correlated with the immune score, PD-L1 expression, and the infiltration of immune cell subtypes in BC. SNHG16, a member of the high-risk group based on m6A-LPS, was highly expressed in BC tissues and cell lines and interfered with siRNA resulted in suppressed proliferation, migration, and invasion in vitro. Our study illustrates the role of m6A-related lncRNAs in BC. The m6A-LPS may be an important regulatory target of the tumor microenvironment (TME) in BC.
| S-EPMC8818872 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common type of lung cancer and is a severe threat to human health. Although many therapies have been applied to LUAD, the long-term survival rate of patients remains unsatisfactory. We aim to find reliable immune microenvironment-related lncRNA biomarkers to improve LUAD prognosis.MethodsESTIMATE analysis was performed to evaluate the degree of immune infiltration of each patient in TAGA LUAD cohort. Correlation analysis was used to identify the immune microenvironment-related lncRNAs. Univariate cox regression analysis, LASSO analysis, and Kaplan Meier analysis were used to construct and validate the prognostic model based on microenvironment-related lncRNAs.ResultsWe obtained 1,178 immune microenvironment-related lncRNAs after correlation analysis. One hundred and eighty of them are independent prognostic lncRNAs. Sixteen key lncRNAs were selected by LASSO method. This lncRNA-based model successfully predicted patients' prognosis in validation cohort, and the risk score was related to pathological stage. Besides, we also found that TP53 had the highest frequency mutation in LUAD, and the mutation of TP53 in the high-risk group, which was identified by our survival model, has a poor prognosis. lncRNA-mRNA co-expression network further suggested that these lncRNAs play a vital role in the prognosis of LUAD.ConclusionHere, we filtered 16 key lncRNAs, which could predict the survival of LUAD and may be potential biomarkers and therapeutic targets.
| S-EPMC8366027 | biostudies-literature