Dataset Information

N6-methyladenosine-related non-coding RNAs are potential prognostic and immunotherapeutic responsiveness biomarkers for bladder cancer.

ABSTRACT:

Background

Bladder cancer (BC) is a commonly occurring malignant tumor of the urinary system, demonstrating high global morbidity and mortality rates. BC currently lacks widely accepted biomarkers and its predictive, preventive, and personalized medicine (PPPM) is still unsatisfactory. N6-methyladenosine (m⁶A) modification and non-coding RNAs (ncRNAs) have been shown to be effective prognostic and immunotherapeutic responsiveness biomarkers and contribute to PPPM for various tumors. However, their role in BC remains unclear.

Methods

m⁶A-related ncRNAs (lncRNAs and miRNAs) were identified through a comprehensive analysis of TCGA, starBase, and m6A2Target databases. Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m⁶A-related ncRNA-based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO). A nomogram containing independent prognostic factors was developed. Differences in BC clinical characteristics, m⁶A regulators, m⁶A-related ncRNAs, gene expression patterns, and differentially expressed genes (DEGs)-associated molecular networks between the high- and low-risk groups in TCGA dataset were also analyzed. Additionally, the potential applicability of the risk model in the prediction of immunotherapeutic responsiveness was evaluated based on the "IMvigor210CoreBiologies" data set.

Results

We identified 183 m⁶A-related ncRNAs, of which 14 were related to OS. LASSO regression analysis was further used to develop a prognostic risk model that included 10 m⁶A-related ncRNAs (BAALC-AS1, MIR324, MIR191, MIR25, AC023509.1, AL021707.1, AC026362.1, GATA2-AS1, AC012065.2, and HCP5). The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively. A nomogram containing 3 independent prognostic factors (risk score, age, and clinical stage) was developed and was found to demonstrate high prognostic prediction accuracy (AUC = 0.83). Moreover, the risk model could also predict BC progression. A higher risk score indicated a higher pathological grade and clinical stage. We identified 1058 DEGs between the high- and low-risk groups in TCGA dataset; these DEGs were involved in 3 molecular network systems, i.e., cellular immune response, cell adhesion, and cellular biological metabolism. Furthermore, the expression levels of 8 m⁶A regulators and 12 m⁶A-related ncRNAs were significantly different between the two groups. Finally, this risk model could be used to predict immunotherapeutic responses.

Conclusion

Our study is the first to explore the potential application value of m⁶A-related ncRNAs in BC. The m⁶A-related ncRNA-based risk model demonstrated excellent performance in predicting prognosis and immunotherapeutic responsiveness. Based on this model, in addition to identifying high-risk patients early to provide them with focused attention and targeted prevention, we can also select beneficiaries of immunotherapy to deliver personalized medical services. Furthermore, the m⁶A-related ncRNAs could elucidate the molecular mechanisms of BC and lead to a new direction for the improvement of PPPM for BC.

Supplementary information

The online version contains supplementary material available at 10.1007/s13167-021-00259-w.

SUBMITTER: Lu M

PROVIDER: S-EPMC8648947 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: The role of N6-methyladenosine (m6A)-related long non-coding RNAs (lncRNAs) in osteosarcoma (OS) has not been fully studied yet. We aimed to identify m6A-related lncRNAs that could act as prognostic biomarkers for OS. Methods: Pearson correlation was performed to identify m6A-related lncRNAs. Univariate and multivariate Cox regression analyses were performed to construct the risk model and assess whether the risk score was an independent prognostic factor for patients with OS. Gene Set Enrichment Analysis (GSEA) was performed to analyze the functions of genes in high-risk and low-risk groups. StarBase and Cytoscape were used to construct a competing endogenous RNA (ceRNA) network based on m6A-related prognostic lncRNA signature. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to analyze the function of genes involved in the ceRNA network. Results: We extracted 122 common lncRNAs from TCGA and Gene Expression Omnibus (GEO) databases. Pearson correlation results revealed 59 significant m6A-related lncRNAs in The Cancer Genome Atlas (TCGA) database, from which 2 were screened to construct a risk signature in TCGA dataset, which was then validated in the GEO dataset. A corresponding risk score was calculated and shown to be an independent prognostic factor for patients with OS. Enrichment analysis indicated that cell proliferation-related biological processes were more common in the high-risk group, while immune-related biological processes were more common in the low-risk group. Moreover, we established a nomogram that had a good ability to predict the overall survival of patients with OS. Additionally, a ceRNA network based on small nucleolar RNA host gene 7 (SNHG7) and small nucleolar RNA host gene 12 (SNHG12) was constructed, with genes that were enriched in hepatocellular carcinoma, gastric cancer, and non-small-cell lung cancer pathways. Conclusion: Our study revealed the prognostic role of m6A-related lncRNAs in OS and identified SNHG7 and SNHG12 as potential biomarkers for predicting the prognosis of patients with OS. These findings have enriched our understanding of the role of m6A modification in the dysregulation of lncRNAs in OS.

Project description:BackgroundBladder 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.MethodsConsensus 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.ResultsA 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.ConclusionsWe 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.

Project description:Background: N6-methyladenosine (m6A) is the most common form of mRNA- and long noncoding RNA (lncRNA)-specific internal modification encountered in eukaryotes, with important effects on mRNA stability, translation, and splicing. The role of m6A-modified lncRNAs (m6A-lncRNAs) in bladder cancer (BLCA) is rarely reported. This study aimed to evaluate an efficient prognostic model of BLCA in patients, based on m6A-lncRNAs, and to discover potential biological targets. Methods: Differentially expressed lncRNAs were investigated in 433 BLCA samples derived from The Cancer Genome Atlas (TCGA) database. Kaplan–Meier and univariate Cox regression analyses were performed to screen for m6A-lncRNAs with prognostic roles in BLCA. We implemented Pearson correlation analysis to analyze 18 potentially prognostic lncRNAs and 20 known m6A-associated genes. Next, the data were imputed using least absolute shrinkage and selection operator (LASSO) Cox regression to establish an m6A-lncRNA prognostic signature. Results: We established an integrated risk score (RS) containing five m6A-lncRNAs and constructed a nomogram that had the ability to forecast the overall survival (OS) of patients with BLCA. We showed that the predictive accuracy of the RS for BLCA prognosis was high, which was confirmed by the area under the receiver operating characteristic (ROC) curve. We analyzed the correlation between tumor immune infiltrating cells and RS in high- and low-risk patients with BLCA and used tumor immune dysfunction and exclusion to predict the effect of immunotherapy. We screened out the most relevant modules of RS through the weighted gene co-expression network analysis network and explored their potential biological functions using GO and KEGG analyses. Conclusion: Our findings demonstrate that, compared with nomograms constructed using a single prognostic factor, the integrated RS represents a superior model for predicting survival in patients with BLCA, which may improve the clinical management of BLCA.