Dataset Information

Comprehensive analysis of metabolism-related lncRNAs related to the progression and prognosis in osteosarcoma from TCGA.

ABSTRACT:

Background

Osteosarcoma is one of the most common malignant neoplasms in children and adolescents. Studies have shown that metabolism-related pathways are vital for the development and metastasis of osteosarcoma. Long non-coding RNA (lncRNA) plays a key role in the occurrence and progression of cancer in a variety of ways. However, the detailed molecular mechanisms of metabolism-related lncRNA in osteosarcoma remain to be deeply elucidated.

Methods

In this study, all metabolism-related mRNAs and lncRNAs in osteosarcoma were extracted and identified based on transcriptomic data from the TCGA database. Usingsurvival analysis, univariate and multivariate independent prognostic analysis, gene set enrichment analysis, and nomogram, a prognostic signature with metabolic lncRNAs as prognostic factors was constructed.

Results

Nine prognostic factors included lncRNA AC009779.2, lncRNA AL591895.1, lncRNA AC026271.3, lncRNA LPP-AS2, lncRNA LINC01857, lncRNA AP005264.1, lncRNA LINC02454, lncRNA AL133338.1, and lncRNA AC135178.5, respectively. Survival analysis indicated that alterations of specific lncRNA expression were strongly correlated with poor prognosis in osteosarcoma. Univariate and multivariate independent prognostic analysis showed that the prognostic signature had a good independent predictive ability for patient survival. The results of GSEA suggested that these predictors may be involved in the metabolism of certain substances or energy in cancer. The nomogram was further drawn for clinical guidance and assistance in clinical decision-making.

Conclusions

This study identified multiple metabolism-related lncRNAs, which may be novel therapeutic targets for osteosarcoma, and contributed to better explore the specific metabolic regulatory mechanisms of lncRNA in osteosarcoma.

SUBMITTER: Chen X

PROVIDER: S-EPMC8381543 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background: Osteosarcoma (OS) is the most common primary tumor of bone in adolescents, and its survival rate is generally less than 20% when metastases occur. Necroptosis, a novel form of programmed necrotic cell death distinct from apoptosis, has been increasingly recognized as a promising therapeutic strategy. This study sought to identify long non-coding RNAs (lncRNAs) associated with necrotizing apoptosis to predict prognosis and target drug use to improve patient survival. Methods: Transcriptomic data and clinical data from 85 OS patients with survival time data and expression profiles from 85 random normal adipose tissue samples were extracted from the UCSC Xena website (http://xena.ucsc.edu/). Nine necroptosis-associated differential prognostic lncRNAs were then identified by analysis of variance, correlation analysis, univariate Cox (uni-Cox) regression, and Kaplan-Meier analysis. Then, patients were randomized into training or testing groups. According to uni-Cox, we obtained prognostic lncRNAs in the training group and intersected them with the abovementioned nine lncRNAs to obtain the final necrotizing apoptosis-related differential prognostic lncRNAs (NRlncRNAs). Next, we performed the least absolute shrinkage and selection operator (LASSO) to construct a risk model of NRlncRNAs. Kaplan-Meier analysis, ROC curves, nomograms, calibration curves, and PCA were used to validate and evaluate the models and grouping. We also analyzed the differences in tumor immunity and drugs between risk groups. Results: We constructed a model containing three NRlncRNAs (AL391121.1, AL354919.2, and AP000851.2) and validated its prognostic predictive power. The value of the AUC curve of 1-, 3-, and 5-year survival probability was 0.806, 0.728, and 0.731, respectively. Moreover, we found that the overall survival time of patients in the high-risk group was shorter than that in the low-risk group. GSEA and ssGSEA showed that immune-related pathways were mainly abundant in the low-risk group. We also validated the differential prediction of immune checkpoint expression, tumor immunity, and therapeutic compounds in the two risk groups. Conclusion: Overall, NRlncRNAs have important functions in OS, and these three NRlncRNAs can predict the prognosis of OS and provide guidance for immunotherapy in OS.

Project description:Osteosarcoma is the most frequent bone tumor. Notwithstanding that significant medical progress has been achieved in recent years, the 5-year overall survival of osteosarcoma patients is inferior. Regulation of fatty acids and lactate plays an essential role in cancer metabolism. Therefore, our study aimed to comprehensively assess the fatty acid and lactate metabolism pattern and construct a fatty acid and lactate metabolism-related risk score system to predict prognosis in osteosarcoma patients. Clinical data and RNA expression data were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and Gene Expression Omnibus (GEO) databases. We used the least absolute shrinkage and selection operator (LASSO) and Cox regression analyses to construct a prognostic risk score model. Relationships between the risk score model and age, gender, tumor microenvironment characteristics, and drug sensitivity were also explored by correlation analysis. We determined the expression levels of prognostic genes in osteosarcoma cells via Western blotting. We developed an unknown fatty acid and lactate metabolism-related risk score system based on three fatty acid and lactate metabolism-related genes (SLC7A7, MYC, and ACSS2). Survival analysis showed that osteosarcoma patients in the low-risk group were likely to have a better survival time than those in the high-risk group. The area under the curve (AUC) value shows that our risk score model performs well in predicting prognosis. Elevated fatty acids and lactate risk scores weaken immune function and the environment of the body, which causes osteosarcoma patients' poor survival outcomes. In general, the constructed fatty acid and lactate metabolism-related risk score model can offer essential insights into subsequent mechanisms in available research. In addition, our study may provide rational treatment strategies for clinicians based on immune correlation analysis and drug sensitivity in the future.

Project description:BackgroundN6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approaches of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between m6A and lncRNAs enrolling in the shaping of TIME remains unclear.MethodsThe RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering method was implemented to classify GC patients into two clusters. Survival analysis, the infiltration level of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. A competing endogenous RNA (ceRNA) network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of pathways in which m6A-related lncRNAs enriched. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal lncRNAs, and risk model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the response to immune checkpoint inhibitors (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect expression patterns of the selected lncRNAs in the 35 tumor tissues and their paired adjacent normal tissues, and validated the prognostic value of risk model in our cohort (N = 35).ResultsThe expression profiles of 15 lncRNAs were included to cluster patients into 2 subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and lower mutation rates of the genes. Different immune cell infiltration patterns were also displayed between the two clusters. GSEA showed that cluster1 preferentially enriched in tumor hallmarks and tumor-related biological pathways. KEGG pathway analysis found that the target mRNAs which m6A-related lncRNAs regulated by sponging miRNAs mainly enriched in vascular smooth muscle contraction, cAMP signaling pathway and cGMP-PKG signaling pathway. Next, eight lncRNAs were selected by LASSO regression algorithm to construct risk model. Patients in the high-risk group had poor prognoses, which were consistent in our cohort. As for predicting responses to ICIs therapy, patients from high-risk group were found to have lower tumor mutation burden (TMB) scores and account for large proportion in the Microsatellite Instability-Low (MSI-L) subtype. Moreover, patients had distinct immunophenoscores in different risk groups.ConclusionOur study revealed that the interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to ICIs therapy.

Project description:Background: Abnormal expression of lncRNA is closely related to the occurrence and metastasis of osteosarcoma. The tumor immune microenvironment (TIM) is considered to be an important factor affecting the prognosis and treatment of osteosarcoma. This study aims to explore the effect of immune-related lncRNAs (IRLs) on the prognosis of osteosarcoma and its relationship with the TIM. Methods: Ninety-five osteosarcoma samples from the TARGET database were included. Iterative LASSO regression and multivariate Cox regression analysis were used to screen the IRLs signature with the optimal AUC. The predict function was used to calculate the risk score and divide osteosarcoma into a high-risk group and low-risk group based on the optimal cut-off value of the risk score. The lncRNAs in IRLs signature that affect metastasis were screened for in vitro validation. Single sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithms were used to evaluate the role of TIM in the influence of IRLs on osteosarcoma prognosis. Results: Ten IRLs constituted the IRLs signature, with an AUC of 0.96. The recurrence and metastasis rates of osteosarcoma in the high-risk group were higher than those in the low-risk group. In vitro experiments showed that knockdown of lncRNA (AC006033.2) could increase the proliferation, migration, and invasion of osteosarcoma. ssGSEA and ESTIMATE results showed that the immune cell content and immune score in the low-risk group were generally higher than those in the high-risk group. In addition, the expression levels of immune escape-related genes were higher in the high-risk group. Conclusion: The IRLs signature is a reliable biomarker for the prognosis of osteosarcoma, and they alter the prognosis of osteosarcoma. In addition, IRLs signature and patient prognosis may be related to TIM in osteosarcoma. The higher the content of immune cells in the TIM of osteosarcoma, the lower the risk score of patients and the better the prognosis. The higher the expression of immune escape-related genes, the lower the risk score of patients and the better the prognosis.

Project description:Background: Cuproptosis is a novel form of cell death discovered in recent. A great quantity of researches has confirmed the close relationships and crucial roles between long non-coding RNAs (lncRNAs) with the progression of colorectal cancer (CRC). However, the relationship between cuproptosis and lncRNAs remains unclear in CRC. Methods: 1,111 co-expressed lncRNAs with 16 cuproptosis regulators were retrieved from CRC samples of The Cancer Genome Atlas (TCGA) database. Through univariate Cox and least absolute shrinkage and selection operator regression analysis, a prognosis model was constructed with 15 lncRNAs. The Kaplan-Meier, receiver operating characteristic curve, C-index and principal component analysis identified the prognostic power. Furthermore, a cuproptosis-related cluster was generated based on the 15 lncRNAs by unsupervised methods. The correlations between the cuproptosis-related signatures with immune cell infiltration and anti-tumor therapy were explored by multiple algorithms. Results: A risk score and nomogram with great prediction ability were constructed for CRC prognosis evaluation. The immune activate pathways, immune infiltration cells, immune functions, immune score and immune activation genes were remarkably enriched in the high risk group. The cuproptosis-related cluster was generated, of which the cluster 2 showed longer overall survival. The immune cell infiltration analysis indicated the similar results of cluster 2 with the high risk group, implying a significant marker for "hot tumor." The cluster 2 also presented high expression of immune checkpoint molecules, MSI-H status and higher susceptibility to multiple immunotherapy drugs. Conclusion: We appraised a novel cuproptosis-related prognosis model and molecular signature associated with prognosis, immune infiltration and immunotherapy. The identification of cuproptosis-related lncRNAs improved our understanding of immune infiltration and provided a significant marker for prognosis and immunotherapy in CRC.