Dataset Information

An effective N6-methyladenosine-related long non-coding RNA prognostic signature for predicting the prognosis of patients with bladder cancer

ABSTRACT:

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.

SUBMITTER: Ma T

PROVIDER: S-EPMC8607649 | 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 (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 (m6A) 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.Methodsm6A-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 m6A-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, m6A regulators, m6A-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.ResultsWe identified 183 m6A-related ncRNAs, of which 14 were related to OS. LASSO regression analysis was further used to develop a prognostic risk model that included 10 m6A-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 m6A regulators and 12 m6A-related ncRNAs were significantly different between the two groups. Finally, this risk model could be used to predict immunotherapeutic responses.ConclusionOur study is the first to explore the potential application value of m6A-related ncRNAs in BC. The m6A-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 m6A-related ncRNAs could elucidate the molecular mechanisms of BC and lead to a new direction for the improvement of PPPM for BC.Supplementary informationThe online version contains supplementary material available at 10.1007/s13167-021-00259-w.

Project description:Long non-coding RNAs (lncRNAs) were reported to potentially play a regulatory role in the process of myocardial regeneration in the neonatal mouse. N6-methyladenosine (m6A) modification may play a key role in myocardial regeneration in mice and regulates a variety of biological processes through affecting the stability of lncRNAs. However, the map of m6A modification of lncRNAs in mouse cardiac development still remains unknown. We aimed to investigate the differences in the m6A status of lncRNAs during mouse cardiac development and reveal a potential role of m6A modification modulating lncRNAs in cardiac development and myocardial regeneration during cardiac development in mice. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) of the heart tissue in C57BL/6 J mice at postnatal day 1 (P1), P7 and P28 were performed to produce stagewise cardiac lncRNA m6A-methylomes in a parallel timeframe with the established loss of an intrinsic cardiac regeneration capacity and early postnatal development. There were significant differences in the distribution and abundance of m6A modifications in lncRNAs in the P7 vs P1 mice. In addition, the functional role of m6A in regulating lncRNA levels was established for selected transcripts with METTL3 silencing in neonatal cardiomyocytes in vitro. Based on our MeRIP-qPCR experiment data, both lncGm15328 and lncRNA Zfp597, that were not previously associated with cardiac regeneration, were found to be the most differently methylated at P1-P7. These two lncRNAs sponged several miRNAs which further regulated multiple mRNAs, including some of which have previously been linked with cardiac regeneration ability. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis revealed that differential m6A modifications were more enriched in functions and cellular signalling pathways related to cardiomyocyte proliferation. Our data suggested that the m6A modification on lncRNAs may play an important role in the regeneration of myocardium and cardiac development.

Project description:Background: Lung adenocarcinoma (LUAD) is the major subtype of lung cancer and is associated with very high mortality. Emerging studies have shown that N6-methyladenosine (m6A)-related long non-coding (lnc) RNAs play crucial roles in tumor prognosis and the tumor immune microenvironment (TME). We aimed to explore the expression patterns of different m6A-related lncRNAs concerning patient prognosis and construct an m6A-related lncRNA prognostic model for LUAD. Methods: The prognostic value of m6A-related lncRNAs was investigated in LUAD samples from The Cancer Genome Atlas (TCGA). Potential prognostic m6A-related lncRNAs were selected by Pearson's correlation and univariate Cox regression analysis. Patients were divided into clusters using principal component analysis and the m6A-related lncRNA prognostic signature was calculated using least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Results: Based on 91 prognostic m6A-related lncRNAs, we identified two m6A-related-lncRNA pattern clusters with different overall survival (OS) and different TMEs. We subsequently verified our findings multidimensionally by constructing a 13 m6A-related lncRNA prognostic signature (m6A-LPS) to calculate the risk score, which was robust in different subgroups. The receiver operating characteristic (ROC) curves and concordance index demonstrated that m6A-LPS harbored a promising ability to predict OS in TCGA data set and independent GSE11969 cohort. The risk score was also related to OS, TME, and clinical stage, and the risk score calculated by our model was also identified as independent prognostic predictive factors for LUAD patients after adjustment for age, smoking, gender, and stage. Enrichment analysis indicated that malignancy and drug resistance-associated pathways were more common in cluster2 (LUAD-unfavorable m6A-LPS). Furthermore, the results indicated that the signaling pathway enriched by the target gene of 13 m6A-related lncRNAs may be associated with metastasis and progression of cancer according to current studies. Conclusion: The current results indicated that different m6A-related-lncRNA patterns could affect OS and TME in patients with LUAD, and the prognostic signature based on 13 m6A-related lncRNAs may help to predict the prognosis in LUAD patients.