Dataset Information

Comprehensive Analysis of N6-Methyladenosine-Related lncRNA Signature for Predicting Prognosis and Immune Cell Infiltration in Patients with Colorectal Cancer

ABSTRACT:

Background

Colorectal cancer (CRC) is the third most common tumor worldwide. Aberrant N6-methyladenosine (m6A) modification can influence the progress of the CRC. Additionally, long noncoding RNA (lncRNA) plays a critical role in CRC and has a close relationship with m6A modification. However, the prognostic potential of m6A-related lncRNAs in CRC patients still remains to be clarified.

Methods

We use “limma” R package, “glmnet” R package, and “survival” R package to screen m6A-related-lncRNAs with prognostic potential. Then, we comprehensively analysed and integrated the related lncRNAs in different TNM stages from TCGA database using the LASSO Cox regression. Meanwhile, the relationship between functional enrichment of m6A-related lncRNAs and immune microenvironment in CRC was also investigated using the TCGA database. A prognostic model was constructed and validated to determine the association between m6A-related lncRNAs in different TNM stages and the prognosis of CRC.

Result

We demonstrated that three related m6A lncRNAs in different TNM stages were associated with the prognosis of CRC patients. Patients from the TCGA database were classified into the low-risk and the high-risk groups based on the expression of these lncRNAs. The patients in the low-risk group had longer overall survival than the patients in the high-risk group (P < 0.001). We further constructed and validated a prognostic nomogram based on these genes with a C-index of 0.80. The receiver operating characteristic curves confirmed the predictive capacity of the model. Meanwhile, we also found that the low-risk group has the correlation with the dendritic cell (DC). Finally, we discovered the relationship between the m6A regulators and the three lncRNAs.

Conclusion

The prognostic model based on three m6A-related lncRNAs exhibits superior predictive performance, providing a novel prognostic model for the clinical evaluation of CRC patients.

SUBMITTER: Zhao Z

PROVIDER: S-EPMC8568524 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The N6-methyladenosine (m6A) RNA modification can modify long non-coding RNAs (lncRNAs), thereby affecting the tumorigenesis and progression of tumors. However, the underlying role of m6A-modified lncRNAs in colorectal cancer (CRC) remains largely unknown. Therefore, our aim was to assess the prognostic value of m6A-modified lncRNAs in CRC patients. The gene expression and clinicopathological data of CRC were extracted from The Cancer Genome Atlas (TCGA) database. Pearson correlation analysis was used to investigate the m6A-modified lncRNAs. Consensus clustering was conducted to identify molecular subtypes of CRC, and the clinical significance of molecular subtypes was identified. The least absolute shrinkage and selection operator analysis (LASSO) was applied to establish a risk signature. Finally, a prognostic nomogram with risk score and clinicopathological variables was established. In total, 29 m6A-modified lncRNAs were identified as prognostic lncRNAs. Two molecular clusters were identified and significant differences were found with respect to clinicopathological features and prognosis. Cluster1 is associated with poor overall survival (OS), down-regulation of Programmed cell death ligand-1 (PD-L1) expression, lower immune score, and less immune cell infiltration. Then, an m6A-modified lncRNA signature for predicting OS was constructed in the TCGA training cohort. The signature demonstrated favorable prediction performance in both training and validation sets. Compared with low-risk patients, patients with high risk showed worse clinical outcomes, lower immune scores, and downregulated PD-L1 expression. Further analysis indicated that the signature was an independent prognostic indicator, and then a prognostic nomogram based on risk score, tumor location, and tumor stage was established. Our study identified a seven m6A-modified lncRNA signature and established a prognostic nomogram that reliably predicts OS in CRC. These findings may improve the understanding of m6A modifications in CRC and provide insights into the prognosis and treatment strategy of CRC.

Project description:BackgroundKidney renal clear cell carcinoma (KIRC) lacks effective prognostic biomarkers and the role and mechanism of N6-methyladenosine (m6A) modification of long noncoding RNAs (lncRNAs) in KIRC remain unclear.MethodsWe extracted standard mRNA-sequencing and clinical data from the TCGA database. The prognostic risk model was obtained by Lasso regression and Cox regression. We randomly divided the samples into training and test sets, each taking half of the cases. Based on Lasso regression and Cox regression for training set, the prognostic risk signature was constructed; risk scores were calculated with the R package "glmnet." Based on the median value of the prognostic risk score, risk scores were calculated for each patient and we divided all KIRC samples into high-risk and low-risk groups. Then, high- and low-risk subtypes were established and their prognosis, clinical features, and immune infiltration microenvironment were evaluated in test set and the entire sampled data set. The reliability of the prognostic model was confirmed by receiver operating characteristic curve analysis.ResultsWe found 28 prognostic m6A-related lncRNAs and established a m6A-related lncRNAs prognostic signature. Risk score=AC015813.1∗(0.0086)+EMX2OS∗(-0.0101)+LINC00173∗(0.0309)+PWAR5∗(-0.0146)+SNHG1∗(0.0043). The signature showed a better predictive ability than other clinical indicators, including tumor node metastasis classification (TNM), histological, and pathological stages. In the high-risk group, M0 macrophages, CD8+ T cells, and regulatory T cells had significantly higher scores. Contrarily, in the low-risk group, activated dendritic cells, M1 macrophages, mast resting cells, and monocytes had significantly higher scores. In the high-risk group, LSECtin was overexpressed. In the low-risk group, PD-L1 was overexpressed. Moreover, high-risk patients may benefit more from AZ628.ConclusionsIn conclusion, prognosis prediction of patients with KIRC and new insights for immunotherapy are provided by the m6A-related lncRNA prognostic signature.

Project description:Background Bladder cancer (BLCA) typically has a poor prognosis due to high relapse and metastasis rates. A growing body of evidence indicates that N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) play crucial roles in the progression of BLCA and the treatment response of patients with BLCA. Therefore, we conducted a comprehensive RNA-seq analysis of BLCA using data from The Cancer Genome Atlas (TCGA) to establish an m6A-related lncRNA prognostic signature (m6A-RLPS) for BLCA. Methods Consensus clustering analysis was used to investigate clusters of BLCA patients with varying prognoses. The least absolute shrinkage and selection operator Cox regression were used to develop the m6A-RLPS. The ESTIMATE and CIBERSORT algorithms were used to evaluate the immune composition. Results A total of 745 m6A-related lncRNAs were identified using Pearson correlation analysis (|R| > 0.4, p < 0.001). Fifty-one prognostic m6A-related lncRNAs were screened using univariate Cox regression analysis. Through consensus clustering analysis, patients were divided into two clusters (clusters 1 and 2) with different overall survival rates and tumor stages based on the differential expression of the lncRNAs. Enrichment analysis demonstrated that terms related to tumor biological processes and immune-related activities were increased in patient cluster 2, which was more likely to exhibit low survival rates. Nine m6A-related prognostic lncRNAs were finally determined and subsequently used to construct the m6A-RLPS, which was verified to be an independent predictor of prognosis using univariate and multivariate Cox regression analyses. Further, a nomogram based on age, tumor stage, and the m6A-RLPS was generated and showed high accuracy and reliability with respect to predicting the survival outcomes of BLCA patients. The prognostic signature was found to be strongly correlated to tumor-infiltrating immune cells and immune checkpoint expression. Conclusions We established a novel m6A-RLPS with a favorable prognostic value for patients with BLCA. We believe that this prognostic signature can provide new insights into the tumorigenesis of BLCA and predict the treatment response in patients with BLCA. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-021-08981-4.

Project description:BackgroundPatients with advanced clear cell renal cell carcinoma (ccRCC) have a poor prognosis and lack effective prognostic biomarkers. N6-methyladenosine-related lncRNAs (m6A-related long noncoding RNAs [lncRNAs]) have been confirmed to be associated with the development of multiple tumors, but its role in ccRCC is not clear.MethodsGene expression data and clinical information of ccRCC patients were extracted from The Cancer Genome Atlas Database. The prognostic m6A-related lncRNAs were obtained by Pearson's correlation analysis and univariate Cox regression analysis. Afterward, the cluster classification and its correlation with prognosis, clinical characteristics, and immunity were analyzed. LASSO regression was used to establish the prognostic risk model. The predictive performance of the prognostic model was evaluated and validated by survival analysis and receiver operating characteristic curve analysis, et al. The expression of immune checkpoints and immune cell infiltration in patients with different risks were systematically analyzed.ResultsA total of 27 prognostic m6A-related lncRNAs were identified. These m6A-related lncRNAs were differentially expressed between tumor and normal tissues. Among them, 24 high-risk m6A-related lncRNAs were overexpressed in Cluster 2 and correlated with poor prognosis, low stromal score, high expression of immune checkpoints, and immunosuppressive cells infiltration. Based upon, a prognostic risk model composed of seven m6A-related lncRNAs was constructed. After a series of analyses, it was proved that this model had good sensitivity and specificity, and could predict the prognosis of patients with different clinical stratification. The expression of PD-1, PD-L1, CTLA-4, LAG-3, TIM-3, and TIGIT were significantly increased in the high-risk patients, and there was a correlation between the risk score and immune cell infiltration.ConclusionsThe seven m6A-related lncRNAs prognostic risk signature showed reliable prognostic predictive power for ccRCC and was associated with the expression of immune checkpoints and immune cell infiltration. This seven m6A-related lncRNAs signature will be helpful in managing ccRCC and guiding individualized immunotherapy.