Comprehensive analysis of metabolism-related lncRNAs related to the progression and prognosis in osteosarcoma from TCGA.
Ontology highlight
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
Similar Datasets
Project description:Despite being significant in various diseases, including cancers, the impact of copper metabolism on osteosarcoma (OS) remains largely unexplored. This study aimed to use bioinformatics analyses to identify a reliable copper metabolism signature that could improve OS patient prognosis prediction, immune landscape understanding, and drug sensitivity. Through nonnegative matrix factorization (NMF) clustering, we revealed distinct prognosis-associated clusters of OS patients based on copper metabolism-related genes (CMRGs), showing differential gene expression linked to immune processes. The risk model, comprising 13 prognostic CMRGs, was established using least absolute shrinkage and selection operator (LASSO) Cox regression, closely associated with the OS microenvironment's immune situation and drug sensitivity. Furthermore, we developed an integrated nomogram, combining the risk score and clinical traits to quantitatively predict OS patient prognosis. The calibration plot, timeROC, and timeROC analyses demonstrated its predictable accuracy and clinical usefulness. Finally, we identified three independent prognostic signatures for OS patients: COX11, AP1B1, and ABCB6. This study confirmed the involvement of CMRGs in OS patient prognosis, immune processes, and drug sensitivity, suggesting their potential as promising prognostic signatures and therapeutic targets for OS.
| S-EPMC10497601 | biostudies-literature
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.
| S-EPMC9178207 | biostudies-literature
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.
| S-EPMC9478861 | biostudies-literature
Project description:5-Methyladenosine (m5C) is a type of epigenetic modification involved in the progression of various cancers. To investigate the role of m5C-related long non-coding RNAs (lncRNAs) in the prognosis and immune cell infiltration in hepatocellular carcinoma (HCC), we obtained patients' clinical information and transcriptome data of HCC from the Cancer Genome Atlas (TCGA) database. We applied Pearson correlation analysis to construct an m5C-related lncRNA-messenger RNA (mRNA) co-expression network. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analysis were employed to establish an m5C-related lncRNA prognostic risk model. We then verified the model using Kaplan-Meier analysis, principal component analysis, as well as univariate and multivariate Cox analyses. The expression of m5C-related lncRNAs was validated in HCC tissues and different cell lines. Combining the risk score and clinicopathological features, a nomogram was established for predicting the overall survival (OS) of HCC patients. Furthermore, gene set enrichment analysis (GSEA) revealed that some tumor-associated pathways were significantly enriched in the high-risk group. Immune cell infiltration analysis demonstrated that the levels of Treg cells, neutrophils, and M2 macrophages were higher in the high-risk group. In addition, patients with high tumor mutation burden (TMB) had worse OS than those with low TMB. We also assessed the immune checkpoint level and chemotherapeutic agent sensibility. Then in vitro experiments were performed to examine the biological function of MKLN1-AS in HCC cells and found that knockdown of MKLN1-AS suppressed the proliferation, migration, and invasion. In conclusion, m5C-related lncRNAs played a critical role in predicting the prognosis of patients with HCC and may serve as new therapeutic targets for HCC patients.
| S-EPMC9630339 | biostudies-literature
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.
| S-EPMC8943990 | biostudies-literature
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.
| S-EPMC8529014 | biostudies-literature
Project description:Dedifferentiation is the main concern associated with radioactive iodine (RAI) refractoriness in patients with papillary thyroid cancer (PTC), and the underlying mechanisms of PTC dedifferentiation remain unclear. The present work aimed to identify a useful signature to indicate dedifferentiation and further explore its role in prognosis and susceptibility to chemotherapy drugs. A total of five prognostic-related DR-lncRNAs were selected to establish a prognostic-predicting model, and corresponding risk scores were closely associated with the infiltration of immune cells and immune checkpoint blockade. Moreover, we built an integrated nomogram based on DR-lncRNAs and age that showed a strong ability to predict the 3- and 5-year overall survival. Interestingly, drug sensitivity analysis revealed that the low-risk group was more sensitive to Bendamustine and TAS-6417 than the high-risk group. In addition, knockdown of DR-lncRNAs (DPH6-DT) strongly promoted cell proliferation, invasion, and migration via PI3K-AKT signal pathway in vitro. Furthermore, DPH6-DT downregulation also increased the expression of vimentin and N-cadherin during epithelial-mesenchymal transition. This study firstly confirms that DR-lncRNAs play a vital role in the prognosis and immune cells infiltration in patients with PTC, as well as a predictor of the drugs' chemosensitivity. Based on our results, DR-lncRNAs can serve as a promising prognostic biomarkers and treatment targets.
| S-EPMC8909347 | biostudies-literature
Project description:Nowadays, primary liver cancer is still a major threat to human health. Anoikis is a particular form of programed cell death that has an inhibitory effect on neoplasm metastasis. Although several prognostic models based on anoikis-related genes for Hepatocellular carcinoma (HCC) have been established, signatures associated with anoikis-related lncRNAs have not been identified. To fill this blank space, the authors built up a prognostic signature and appraised its value in guiding immunotherapy. Eleven prognostic anoikis-related lncRNAs were identified through Least Absolute Shrinkage and Selection Operator Cox analysis. The accuracy of the risk signature in predicting prognosis was verified by K-M survival analysis and Receiver operating characteristic analysis. We further discovered that the high-risk group was often enriched in signal pathways related to cell growth and death and immune response; in addition, in the low-risk group, cells often undergo metabolic changes through gene set enrichment analysis. Finally, we realised that HCC patients in the high-risk group were upregulated in immune-checkpoint molecules and tend to have a higher tumour mutation burden level which indicated a higher sensitivity to immunotherapy. All in all, the anoikis-related lncRNAs risk signature showed excellent ability in predicting prognosis and may guide the application of immunotherapy in future clinical practice.
| S-EPMC10439496 | biostudies-literature
Project description:Osteosarcoma is a common malignant bone tumor in children and adolescents. The overall survival of osteosarcoma patients is remarkably poor. Herein, we sought to establish a reliable risk prognostic model to predict the prognosis of osteosarcoma patients. Patients ' RNA expression and corresponding clinical data were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and Gene Expression Omnibus databases. A consensus clustering was conducted to uncover novel molecular subgroups based on 200 hypoxia-linked genes. A hypoxia-risk models were established by Cox regression analysis coupled with LASSO regression. Functional enrichment analysis, including Gene Ontology annotation and KEGG pathway analysis, were conducted to determine the associated mechanisms. Moreover, we explored relationships between the risk scores and age, gender, tumor microenvironment, and drug sensitivity by correlation analysis. We identified two molecular subgroups with significantly different survival rates and developed a risk model based on 12 genes. Survival analysis indicated that the high-risk osteosarcoma patients likely have a poor prognosis. The area under the curve (AUC) value showed the validity of our risk scoring model, and the nomogram indicates the model's reliability. High-risk patients had lower Tfh cell infiltration and a lower stromal score. We determined the abnormal expression of three prognostic genes in osteosarcoma cells. Sunitinib can promote osteosarcoma cell apoptosis with down-regulation of KCNJ3 expression. In summary, the constructed hypoxia-related risk score model can assist clinicians during clinical practice for osteosarcoma prognosis management. Immune and drug sensitivity analysis can provide essential insights into subsequent mechanisms. KCNJ3 may be a valuable prognostic marker for osteosarcoma development.
| S-EPMC9853159 | biostudies-literature
Project description:N7-Methylguanosine (m7G) and long non-coding RNAs (lncRNAs) have been widely reported to play an important role in cancer. However, there is little known about the relationship between m7G-related lncRNAs and clear cell renal cell carcinoma (ccRCC). To find new potential biomarkers and construct an m7G-related lncRNA prognostic signature for ccRCC, we retrieved transcriptome data and clinical data from The Cancer Genome Atlas (TCGA), and divided the entire set into train set and test set with the ratio of 1:1 randomly. The m7G-related lncRNAs were identified by Pearson correlation analysis (|coefficients| > 0.4, and p < 0.001). Then we performed the univariate Cox regression and least absolute shrinkage and selection operator (LASSO) Cox regression analysis to construct a 12 m7G-related lncRNA prognostic signature. Next, principal component analysis (PCA), the Kaplan-Meier method, time-dependent receiver operating characteristics (ROC) were made to verify and evaluate the risk signature. A nomogram based on the risk signature and clinical parameters was developed and showed high accuracy and reliability for predicting the overall survival (OS). Functional enrichment analysis (GO, KEGG and GSEA) was used to investigate the potential biological pathways. We also performed the analysis of tumor mutation burden (TMB), immunological analysis including immune scores, immune cell infiltration (ICI), immune function, tumor immune escape (TIE) and immunotherapeutic drug in our study. In conclusion, using the 12 m7G-related lncRNA risk signature as a prognostic indicator may offer us insight into the oncogenesis and treatment response prediction of ccRCC.
| S-EPMC9043611 | biostudies-literature