Identification and Development of a Novel 4-Gene Immune-Related Signature to Predict Osteosarcoma Prognosis
Ontology highlight
ABSTRACT:
SUBMITTER: Cao M
PROVIDER: S-EPMC7785859 | biostudies-literature | 2020 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundIn childhood and adolescence, the prevailing bone tumor is osteosarcoma associated with frequent recurrence and lung metastasis. This research focused on predicting the survival and immune landscape of osteosarcoma by developing a prognostic signature and establishing aging-related genes (ARGs) subtypes.MethodsThe training group comprised of the transcriptomic and associated clinical data of 84 patients with osteosarcoma accessed at the TARGET database and the validation group consisted of 53 patients from GSE21257. The aging-related subtypes were identified using unsupervised consensus clustering analysis. The ARG signature was developed utilizing multivariate Cox analysis and LASSO regression. The prognostic value was assessed using the univariate and multivariate Cox analyses, Kaplan-Meier plotter, time-dependent ROC curve, and nomogram. The functional enrichment analyses were performed by GSEA, GO, and KEGG analysis, while the ssGSEA, ESTIMATE, and CIBERSORT analyses were conducted to reveal the immune landscape in osteosarcoma.ResultsThe two clusters of osteosarcoma patients formed based on 543 ARGs, depicted a considerable difference in the tumor microenvironment, and the overall survival and immune cell infiltration rate varied as well. Among these, the selected 23 ARGs were utilized for the construction of an efficient predictive prognostic signature for the overall survival prediction. The testing in the validation group of osteosarcoma patients confirmed the status of the high-risk score as an independent indicator for poor prognosis, which was already identified as such using the univariate and multivariate Cox analyses. Furthermore, the ARG signature could distinguish different immune-related functions, infiltration status of immune cells, and tumor microenvironment, as well as predict the immunotherapy response of osteosarcoma patients.ConclusionThe aging-related subtypes were identified and a prognostic signature was developed in this research, which determined different prognoses and allowed for treatment of osteosarcoma patients to be tailored. Additionally, the immunotherapeutic response of individuals with osteosarcoma could also be predicted by the ARG signature.
| S-EPMC9581613 | biostudies-literature
Project description:BackgroundAnoikis is a specialized form of programmed apoptosis that occurs in two model epithelial cell lines and plays an important role in tumors. However, the prognostic value of anoikis-related lncRNA (ARLncs) in osteosarcoma (OS) has not been reported.MethodsBased on GTEx and TARGET RNA sequencing data, we carried out a thorough bioinformatics analysis. The 27 anoikis-related genes were obtained from the Gene Set Enrichment Analysis (GSEA). Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analysis were successively used to screen for prognostic-related ARLncs. To create the prognostic signature of ARLncs, we performed multivariate Cox regression analysis. We calculated the risk score based on the risk coefficient, dividing OS patients into high- and low-risk subgroups. Additionally, the relationship between the OS immune microenvironment and risk prognostic models was investigated using function enrichment, including Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), single-sample gene set enrichment analysis (ssGSEA), and GSEA analysis. Finally, the potential effective drugs in OS were found by immune checkpoint and drug sensitivity screening.ResultsA prognostic signature consisting of four ARLncs (AC079612.1, MEF2C-AS1, SNHG6, and TBX2-AS1) was constructed. To assess the regulation patterns of anoikis-related lncRNA genes, we created a risk score model. According to a survival analysis, high-risk patients have a poor prognosis as they progress. By using immune functional analysis, the lower-risk group demonstrated the opposite effects compared with the higher-risk group. GO and KEGG analysis showed that the ARLncs pathways and immune-related pathways were enriched. Immune checkpoints and drug sensitivity analysis might be used to determine the better effects of the higher group.ConclusionWe identified a novel prognostic model based on a four-ARLncs signature that might serve as potential prognostic indicators that can be used to predict the prognosis of OS patients, and immunotherapy and drugs that may contribute to improving the overall survival of OS patients and advance our understanding of OS.
| S-EPMC10076677 | biostudies-literature
Project description:Background: Gastric cancer (GC) is the fifth most common malignancy and the third leading cause of tumor-related deaths globally. Herein, we attempted to build a novel immune-related gene (IRG) signature that could predict the prognosis and immunotherapeutic efficiency for GC patients. Methods: The mRNA transcription data and corresponding clinical data of GC were downloaded from The Cancer Genome Atlas (TCGA) database as the training group and the GSE84437 data set as the testing cohort, followed by acquisition of IRGs from the InnateDB resource and ImmPort database. Using the univariate Cox regression analysis, an IRG signature was developed. Several immunogenomic analyses were performed to illustrate the associations between the immune risk score and tumor mutational burden, immune cell infiltrations, function of immune infiltration, clinical characteristics, immune subtype, and immunotherapeutic response. Results: The analysis of 343 GC samples and 30 normal samples from the TCGA database gave rise to 8,713 differentially expressed genes (DEGs) and 513 differentially expressed immune-related genes (DEIRGs) were extracted. The novel IRG signature contained eight DEIRGs (FABP4, PI15, RNASE2, CGB5, INHBE, RLN2, DUSP1, and CD36) and was found to serve as an independent predictive and prognostic factor for GC. Then, the GC patients were separated into the high- and low-risk groups based on the median risk score, wherein the low-risk group presented a better prognosis and was more sensitive to immunotherapy than did the high-risk group. According to the time-dependent ROC curves and AUCs, the immunotherapeutic value of the signature was better than the Tumor Immune Dysfunction and Exclusion (TIDE) and T-cell inflammatory signature (TIS) scores. In addition, the AUCs of the risk score for predicting 1-, 2-, and 3-year OS were 0.675, 0.682, and 0.710, respectively, which indicated that the signature had great predictive power. Conclusion: This study presents a novel IRG signature based on the tumor immune microenvironment, which could improve the prediction of the prognosis and immunotherapeutic efficiency for GC patients. The powerful signature may serve as novel biomarkers and provide therapeutic targets for precision oncology in clinical practice.
| S-EPMC9186121 | biostudies-literature
Project description:BackgroundCellular senescence, a novel hallmark of cancer, is associated with patient outcomes and tumor immunotherapy. However, at present, there is no systematic study on the use of cellular senescence-related long non-coding RNAs (CSR-lncRNAs) to predict survival in patients with osteosarcoma. In this study, we aimed to identify a CSR-lncRNAs signature and to evaluate its potential use as a survival prognostic marker and predictive tool for immune response of osteosarcoma.MethodsWe downloaded a cohort of patients with osteosarcoma from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. We performed differential expression and co-expression analyses to identify CSR-lncRNAs. We performed univariate and multivariate Cox regression analyses along with the random forest algorithm to identify lncRNAs significantly correlated with senescence. Subsequently, we assessed the predictive models using survival curves, receiver operating characteristic curves, nomograms, C-index, and decision curve analysis. Based on this model, patients with osteosarcoma were divided into two groups according to their risk scores. Then, using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, we compared their clinical characteristics to uncover functional differences. We further conducted immune infiltration analyses using estimation of stromal and immune cells in malignant tumor tissues using expression data (ESTIMATE), cell-type identification by estimating relative subsets of rna transcripts (CIBERSORT), and single-sample gene set enrichment analysis for the two groups. We also evaluated the expression of the target genes of immune checkpoint inhibitors (ICIs).ResultsWe identified six lncRNAs that were significantly correlated with senescence and accordingly established a novel cellular senescence-related lncRNA prognostic signature incorporating these lncRNAs. The nomogram indicated that the risk model was an independent prognostic factor that could predict the survival of patients with osteosarcoma. This model demonstrated high accuracy upon validation. Further analysis revealed that patients with osteosarcoma in the low-risk group exhibited better clinical outcomes and enhanced immune infiltration.ConclusionsThe six-CSR-lncRNA prognostic signature effectively predicted survival outcomes and patients in the low-risk group might have improved immune infiltration.
| S-EPMC11319968 | biostudies-literature
Project description:Glycolysis ensures energy supply to cancer cells, thereby facilitating tumor progression. Here, we identified glycolysis-related genes that could predict the prognosis of patients with osteosarcoma. We examined 198 glycolysis-related genes that showed differential expression in metastatic and non-metastatic osteosarcoma samples in the TARGET database, and identified three genes (P4HA1, ABCB6, and STC2) for the establishment of a risk signature. Based on the signature, patients in the high-risk group had poor outcomes. An independent Gene Expression Omnibus database GSE21257 was selected as the validation cohort. Receiver operating characteristic curve analysis was performed and the accuracy of predicting the 1- and 3-year survival rates was shown by the areas under the curve. The results were 0.884 and 0.790 in the TARGET database, and 0.740 and 0.759 in the GSE21257, respectively. Furthermore, we applied ESTIMATE algorithm and performed single sample gene set enrichment analysis to compare tumor immunity between high- and low-risk groups. We found that the low-risk group had higher immune scores and immune infiltration levels than the high-risk group. Finally, we chose P4HA1 as a representative gene to verify the function of risk genes in vitro and in vivo and found that P4HA1 could promote the metastasis of osteosarcoma cells. Our study established a novel glycolysis-related risk signature that could predict the prognosis of patients with osteosarcoma.
| S-EPMC8148463 | biostudies-literature
Project description:BackgroundInflammatory response took part in the progression of tumor and was regarded as the hallmark of cancer. However, the prognostic relationship between osteosarcoma and inflammatory response-associated genes (IRGs) was unclear. This research aimed to explore the correlations between osteosarcoma prognosis and IRG signature.MethodsThe inflammatory response-associated differentially expressed messenger RNAs (DEmRNAs) were screened out through Gene Expression Omnibus (GEO) and Molecular Signature Database (MSigDB) databases. Univariate and multivariate cox regression analyses were utilized to construct the IRG signature. The prognostic value of signature was investigated through Kaplan-Meier (KM) survival curve and nomogram. DEmRNAs among high and low inflammatory response-associated risks were identified and functional enrichment analyses were conducted. ESTIMATE, CIBERSORT and single-sample gene set enrichment analyses (ssGSEA) were implied to reveal the alterations in immune infiltration. All the above results were validated in Target database. The expression of IRGs was also validated in different cell lines by quantitative real-time PCR (qRT-PCR) and osteosarcoma patient samples by immunohistochemistry.ResultsThe IRG signature that consisted of two genes (MYC, CLEC5A) was established. In training and validation datasets, patients with lower risk scores survived longer and the IRG signature was confirmed as the independent prognostic factor in osteosarcoma. The nomogram was constructed and the calibration curves demonstrated the reliability of this model. Functional analysis of risk score-associated DEmRNAs indicated that immune-related pathways and functions were significantly enriched. ssGSEA revealed that 14 immune cells and 11 immune functions were significantly dysregulated. The qRT-PCR results indicated IRGs were significantly differently expressed in osteosarcoma and osteoblast cell lines. The immunohistochemistry analyses of patients' samples revealed the same result.ConclusionThe novel osteosarcoma inflammatory response-associated prognostic signature was established and validated in this study. This model could serve as the biomarker and therapeutic target for osteosarcoma in the future.
| S-EPMC8668229 | biostudies-literature
Project description:Immune-related genes (IRGs) are responsible for osteosarcoma (OS) initiation and development. We aimed to develop an optimal IRGs-based signature to assess of OS prognosis. Sample gene expression profiles and clinical information were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) and Genotype-Tissue Expression (GTEx) databases. IRGs were obtained from the ImmPort database. R software was used to screen differentially expressed IRGs (DEIRGs) and functional correlation analysis. DEIRGs were analyzed by univariate Cox regression and iterative LASSO Cox regression analysis to develop an optimal prognostic signature, and the signature was further verified by independent cohort (GSE39055) and clinical correlation analysis. The analyses yielded 604 DEIRGs and 10 hub IRGs. A prognostic signature consisting of 13 IRGs was constructed, which strikingly correlated with OS overall survival and distant metastasis (p?<?0.05, p?<?0.01), and clinical subgroup showed that the signature's prognostic ability was independent of clinicopathological factors. Univariate and multivariate Cox regression analyses also supported its prognostic value. In conclusion, we developed an IRGs signature that is a prognostic indicator in OS patients, and the signature might serve as potential prognostic indicator to identify outcome of OS and facilitate personalized management of the high-risk patients.
| S-EPMC7591524 | biostudies-literature
Project description:ObjectiveOsteosarcoma (OS) is a common bone malignancy with poor prognosis. We aimed to investigate the relationship between cuproptosis-related lncRNAs (CRLncs) and the survival outcomes of patients with OS.MethodsTranscriptome and clinical data of 86 patients with OS were downloaded from The Cancer Genome Atlas (TCGA). The GSE16088 dataset was downloaded from the Gene Expression Omnibus (GEO) database. The 10 cuproptosis-related genes (CRGs) were obtained from a recently published article on cuproptosis in Science. Combined analysis of OS transcriptome data and the GSE16088 dataset identified differentially expressed CRGs related to OS. Next, pathway enrichment analysis was performed. Co-expression analysis obtained CRLncs related to OS. Univariate COX regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were used to construct the risk prognostic model of CRLncs. The samples were divided evenly into training and test groups to verify the accuracy of the model. Risk curve, survival, receiver operating characteristic (ROC) curve, and independent prognostic analyses were performed. Next, principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE) analysis were performed. Single-sample gene set enrichment analysis (ssGSEA) was used to explore the correlation between the risk prognostic models and OS immune microenvironment. Drug sensitivity analysis identified drugs with potential efficacy in OS. Real-time quantitative PCR, Western blotting, and immunohistochemistry analyses verified the expression of CRGs in OS. Real-time quantitative PCR was used to verify the expression of CRLncs in OS.ResultsSix CRLncs that can guide OS prognosis and immune microenvironment were obtained, including three high-risk CRLncs (AL645608.6, AL591767.1, and UNC5B-AS1) and three low-risk CRLncs (CARD8-AS1, AC098487.1, and AC005041.3). Immune cells such as B cells, macrophages, T-helper type 2 (Th2) cells, regulatory T cells (Treg), and immune functions such as APC co-inhibition, checkpoint, and T-cell co-inhibition were significantly downregulated in high-risk groups. In addition, we obtained four drugs with potential efficacy for OS: AUY922, bortezomib, lenalidomide, and Z.LLNle.CHO. The expression of LIPT1, DLAT, and FDX1 at both mRNA and protein levels was significantly elevated in OS cell lines compared with normal osteoblast hFOB1.19. The mRNA expression level of AL591767.1 was decreased in OS, and that of AL645608.6, CARD8-AS1, AC005041.3, AC098487.1, and UNC5B-AS1 was upregulated in OS.ConclusionCRLncs that can guide OS prognosis and the immune microenvironment and drugs that may have a potential curative effect on OS obtained in this study provide a theoretical basis for OS survival research and clinical decision-making.
| S-EPMC9373797 | biostudies-literature
Project description:PurposeBladder cancer (BCa) is one of the most common malignant tumors in the urogenital system, characterized by the high recurrence rate, mortality rate and poor prognosis. Based on cuproptosis-related long noncoding RNAs (CRLs), this study set out to create a prediction signature to evaluate the prognosis of patients with BCa.MethodsRNA-seq data including CRLs and related clinicopathological data were gathered from The Cancer Genome Atlas (TCGA) database (n = 428). The predictive signature was constructed after correlation analysis. Subsequently, relying on the analyzed data from the TCGA database and our sample collection, we examined and verified the connections between CRLs model and important indexes included prognosis, route and functional enrichment, tumor immune evasion, tumor mutation, and treatment sensitivity.ResultsPatients in the high-risk group had lower overall survival (OS) than that of low-risk group. Compared with clinicopathological variables, CRLs features have better predictive value according to receiver operating characteristic (ROC) curve. The expression level of CRLs was highly associated with the tumor progress, tumor microenvironment and tumor immune escape. Additionally, we identified that the mutation of TP53, TTN, KMT2D and MUC16 gene were founded in patients with BCa. Lapatinib, pazopanib, saracatinib, gemcitabine, paclitaxel and palenolactone had good antitumor effects for BCa patients in the high-risk group (all P < 0.001).ConclusionThis study revealed the effects of CRLs on BCa and further established CRLs model, which can be used in clinic for predicting prognosis, immunological response and treatment sensitivity inpatient with BCa.
| S-EPMC9715909 | biostudies-literature
Project description:Host immune dysregulation involves in the initiation and development of osteosarcoma (OS). However, the exact role of immune cells in OS remains unknown. We aimed to distinguish the molecular subtypes and establish a prognostic model in OS patients based on immunocyte infiltration. The gene expression profile and corresponding clinical feature of OS patients were obtained from TARGET and GSE21257 datasets. MCP-counter and univariate Cox regression analyses were applied to identify immune cell infiltration-related molecular subgroups. Functional enrichment analysis and immunocyte infiltration analysis were performed between two subgroups. Furthermore, Cox regression and LASSO analyses were performed to establish the prognostic model for the prediction of prognosis and metastasis in OS patients. The subgroup with low infiltration of monocytic lineage (ML) was related to bad prognosis in OS patients. 435 DEGs were screened between the two subgroups. Functional enrichment analysis revealed these DEGs were involved in immune- and inflammation-related pathways. Three important genes (including TERT, CCDC26, and IL2RA) were identified to establish the prognostic model. The risk model had good prognostic performance for the prediction of metastasis and overall survival in OS patients. A novel stratification system was established based on ML-related signature. The risk model could predict the metastasis and prognosis in OS patients. Our findings offered a novel sight for the prognosis and development of OS.
| S-EPMC10908810 | biostudies-literature