Prognostic Signature of Osteosarcoma Based on 14 Autophagy-Related Genes

From simple factors to artificial intelligence: evolution of prognosis prediction in childhood cancer: a systematic review and meta-analysis

Background

Current paediatric cancer care requires innovative approaches to predict prognosis that facilitates personalised stratification, yet studies on the performance, composition and limitations of contemporary prognostic models are lacking. We aimed to compare the accuracy of traditional and advanced prognostic models.

Methods

A systematic search for this systematic review and meta-analysis (CRTN42022370251) was conducted in PubMed, Embase, Scopus, and the Cochrane Library databases on 28 June 2024. Studies on the accuracy of prognostic markers or models used in paediatric haematological malignancies, central nervous system (CNS), or non-CNS solid tumours (NCNSST) were included. Three model categories were defined using: 1-clinical parameters, 2-genomic-transcriptomic data, and 3-artificial intelligence (AI). Primary outcomes were area under the receiver operating characteristic curve with a 95% confidence interval (CI) for various overall survival intervals and event-free survival. Two independent groups performed selection and data extraction. We used data published by the authors and publicly available databases.

Findings

Of 12,982 studies, 358 were included in the meta-analysis and 27 in the systematic review, with limited data on AI-approaches. Most data were reported on NCNSST at 5-year OS, where a statistically significant difference was observed between Category-1 (0.75 CI: 0.72-0.79) and Category-2 (0.85 CI: 0.82-0.88) (p < 0.001), but not between Categories-2 and -3 (p = 0.2834) (0.82 CI: 0.77-0.88). Internal validation studies showed significantly better performance compared to those using external validation, highlighting the high risk of bias (ROB) inherent in internal validation. High ROB was most commonly experienced in the outcomes and statistical analysis domains, assessed using PROBAST and QUIPS.

Interpretation

It is advisable to introduce Category-2 and -3 models in a clinical setting, especially for NCNSST prognostic for aiding risk-stratification. Although AI-supported predictions in paediatric oncology are at an early stage of development, it is imperative to further explore their potential. This requires structured data collection and ethical sharing from paediatric oncology patients in sufficient quantity and quality.

Funding

None.

Development and validation of a pyroptosis-related prognostic signature associated with osteosarcoma metastasis and immune infiltration

Pyroptosis is a programmed cell death, which has garnered increasing attention because it relates to the immune and therapy response. However, few studies focus on the application of pyroptosis-related genes (PRGs) in predicting osteosarcoma (OS) patients' prognoses. In this study, the gene expression and clinical information of OS patients were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. Based on these PRGs and unsupervised clustering analysis, all OS samples can be classified into 2 clusters. The 8 key differential expressions for PRGs (LAG3, ITGAM, CCL2, TLR4, IL2RA, PTPRC, FCGR2B, and CD5) were established through the univariate Cox regression and utilized to calculate the risk score of all samples. According to the 8-gene signature, OS samples can be divided into high and low-risk groups and correlation analysis can be performed using immune cell infiltration and immune checkpoints. Finally, we developed a nomogram to improve the PRG-predictive model in clinical application. We verified the predictive performance using receiver operating characteristic (ROC) and calibration curves. There were significant differences in survival, immune cell infiltration and immune checkpoints between the low and high-risk groups. A nomogram was developed with clinical indicators and the risk scores were effective in predicting the prognosis of patients with OS. In this study, a prognostic model was constructed based on 8 PRGs were proved to be independent prognostic factors of OS and associated with tumor immune microenvironment. These 8 prognostic genes were involved in OS development and may serve as new targets for developing therapeutic drugs.

Develop a Novel Signature to Predict the Survival and Affect the Immune Microenvironment of Osteosarcoma Patients: Anoikis-Related Genes

Objective

Osteosarcoma (OS) represents a prevalent primary bone neoplasm predominantly affecting the pediatric and adolescent populations, presenting a considerable challenge to human health. The objective of this investigation is to develop a prognostic model centered on anoikis-related genes (ARGs), with the aim of accurately forecasting the survival outcomes of individuals diagnosed with OS and offering insights into modulating the immune microenvironment.

Methods

The study's training cohort comprised 86 OS patients sourced from The Cancer Genome Atlas database, while the validation cohort consisted of 53 OS patients extracted from the Gene Expression Omnibus database. Differential analysis utilized the GSE33382 dataset, encompassing three normal samples and 84 OS samples. Subsequently, the study executed gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses. Identification of differentially expressed ARGs associated with OS prognosis was carried out through univariate COX regression analysis, followed by LASSO regression analysis to mitigate overfitting risks and construct a robust prognostic model. Model accuracy was assessed via risk curves, survival curves, receiver operating characteristic curves, independent prognostic analysis, principal component analysis, and t-distributed stochastic neighbor embedding (t-SNE) analysis. Additionally, a nomogram model was devised, exhibiting promising potential in predicting OS patient prognosis. Further investigations incorporated gene set enrichment analysis to delineate active pathways in high- and low-risk groups. Furthermore, the impact of the risk prognostic model on the immune microenvironment of OS was evaluated through tumor microenvironment analysis, single-sample gene set enrichment analysis (ssGSEA), and immune infiltration cell correlation analysis. Drug sensitivity analysis was conducted to identify potentially effective drugs for OS treatment. Ultimately, the verification of the implicated ARGs in the model construction was conducted through the utilization of real-time quantitative polymerase chain reaction (RT-qPCR).

Results

The ARGs risk prognostic model was developed, comprising seven high-risk ARGs (CBS, MYC, MMP3, CD36, SCD, COL13A1, and HSP90B1) and four low-risk ARGs (VASH1, TNFRSF1A, PIP5K1C, and CTNNBIP1). This prognostic model demonstrates a robust capability in predicting overall survival among patients. Analysis of immune correlations revealed that the high-risk group exhibited lower immune scores compared to the low-risk group within our prognostic model. Specifically, CD8+ T cells, neutrophils, and tumor-infiltrating lymphocytes were notably downregulated in the high-risk group, alongside significant downregulation of checkpoint and T cell coinhibition mechanisms. Additionally, three immune checkpoint-related genes (CD200R1, HAVCR2, and LAIR1) displayed significant differences between the high- and low-risk groups. The utilization of a nomogram model demonstrated significant efficacy in prognosticating the outcomes of OS patients. Furthermore, tumor metastasis emerged as an independent prognostic factor, suggesting a potential association between ARGs and OS metastasis. Notably, our study identified eight drugs-Bortezomib, Midostaurin, CHIR.99021, JNK.Inhibitor.VIII, Lenalidomide, Sunitinib, GDC0941, and GW.441756-as exhibiting sensitivity toward OS. The RT-qPCR findings indicate diminished expression levels of CBS, MYC, MMP3, and PIP5K1C within the context of OS. Conversely, elevated expression levels were observed for CD36, SCD, COL13A1, HSP90B1, VASH1, and CTNNBIP1 in OS.

Conclusion

The outcomes of this investigation present an opportunity to predict the survival outcomes among individuals diagnosed with OS. Furthermore, these findings hold promise for progressing research endeavors focused on prognostic evaluation and therapeutic interventions pertaining to this particular ailment.

Novel pyroptosis-related lncRNAs and ceRNAs predict osteosarcoma prognosis and indicate immune microenvironment signatures

Objective

To study pyroptosis-related biomarkers that are associated with the prognosis and immune microenvironment characteristics of osteosarcoma (OS). The goal is to establish a foundation for the prognosis and treatment of OS.

Methods

We retrieved transcriptome and clinical data from The Cancer Genome Atlas (TCGA) database for 88 OS patients. Using this data, we constructed a prognostic model to identify pyroptosis-related genes (PRGs) associated with OS prognosis. To further explore the biological function of these PRGs, we performed enrichment analysis. To identify pyroptosis-related long non-coding RNAs (PRLncs) associated with the prognosis of OS, we performed co-expression analysis. Subsequently, a risk prognostic model was constructed using these PRLncs to generate a risk score, termed as PRLncs-score, thereby obtaining PRLncs associated with the prognosis of OS. The accuracy of the prognostic model was verified through survival analysis, risk curve, independent prognostic analysis, receiver operating characteristic (ROC) curve, difference analysis between high- and low-risk groups, and clinical correlation analysis. And to determine whether PRLncs-score is independent prognostic factor for OS. In addition, we further conducted external and internal validation for the risk prognosis model. Further analyses of immune cell infiltration and tumor microenvironment were performed. A pyroptosis-related competitive endogenous RNA (PRceRNA) network was constructed to obtain PRceRNAs associated with the prognosis of OS and performed gene set enrichment analysis (GSEA) on PRceRNA genes.

Results

We obtained five PRGs (CHMP4C, BAK1, GSDMA, CASP1, and CASP6) that predicted OS prognosis and seven PRLncs (AC090559.1, AP003119.2, CARD8-AS1, AL390728.4, SATB2-AS1, AL133215.2, and AC009495.3) and one PRceRNA (CARD8-AS1-hsa-miR-21-5p-IL1B) that predicted OS prognosis and indicated characteristics of the OS immune microenvironment. The PRLncs-score, in combination with other clinical features, was established as an independent prognostic factor for OS patients. Subsequent scrutiny of the tumor microenvironment and immune infiltration indicated that patients with low-PRLncs-scores were associated with reduced metastatic risk, improved survival rates, heightened levels of immune cells and stroma, and increased immune activity compared to those with high-PRLncs-scores.

Conclusion

The study's findings offer insight into the prognosis of OS and its immune microenvironment, and hold promise for improving early diagnosis and immunotherapy.

A glycometabolic gene signature associating with immune infiltration and chemosensitivity and predicting the prognosis of patients with osteosarcoma

Background

Accumulating evidence has suggested that glycometabolism plays an important role in the pathogenesis of tumorigenesis. However, few studies have investigated the prognostic values of glycometabolic genes in patients with osteosarcoma (OS). This study aimed to recognize and establish a glycometabolic gene signature to forecast the prognosis, and provide therapeutic options for patients with OS.

Methods

Univariate and multivariate Cox regression, LASSO Cox regression, overall survival analysis, receiver operating characteristic curve, and nomogram were adopted to develop the glycometabolic gene signature, and further evaluate the prognostic values of this signature. Functional analyses including Gene Ontology (GO), kyoto encyclopedia of genes and genomes analyses (KEGG), gene set enrichment analysis, single-sample gene set enrichment analysis (ssGSEA), and competing endogenous RNA (ceRNA) network, were used to explore the molecular mechanisms of OS and the correlation between immune infiltration and gene signature. Moreover, these prognostic genes were further validated by immunohistochemical staining.

Results

A total of four genes including PRKACB, SEPHS2, GPX7, and PFKFB3 were identified for constructing a glycometabolic gene signature which had a favorable performance in predicting the prognosis of patients with OS. Univariate and multivariate Cox regression analyses revealed that the risk score was an independent prognostic factor. Functional analyses indicated that multiple immune associated biological processes and pathways were enriched in the low-risk group, while 26 immunocytes were down-regulated in the high-risk group. The patients in high-risk group showed elevated sensitivity to doxorubicin. Furthermore, these prognostic genes could directly or indirectly interact with other 50 genes. A ceRNA regulatory network based on these prognostic genes was also constructed. The results of immunohistochemical staining showed that SEPHS2, GPX7, and PFKFB3 were differentially expressed between OS tissues and adjacent normal tissues.

Conclusion

The preset study constructed and validated a novel glycometabolic gene signature which could predict the prognosis of patients with OS, identify the degree of immune infiltration in tumor microenvironment, and provide guidance for the selection of chemotherapeutic drugs. These findings may shed new light on the investigation of molecular mechanisms and comprehensive treatments for OS.

Molecular Subtyping and Survival Analysis of Osteosarcoma Reveals Prognostic Biomarkers and Key Canonical Pathways

Osteosarcoma (OS) is a common bone malignancy in children and adolescents. Although histological subtyping followed by improved OS treatment regimens have helped achieve favorable outcomes, a lack of understanding of the molecular subtypes remains a challenge to characterize its genetic heterogeneity and subsequently to identify diagnostic and prognostic biomarkers for developing effective treatments. In the present study, global analysis of DNA methylation, and mRNA and miRNA gene expression in OS patient samples were correlated with their clinical characteristics. The mucin family of genes, MUC6, MUC12, and MUC4, were found to be highly mutated in the OS patients. Results revealed the enrichment of molecular pathways including Wnt signaling, Calcium signaling, and PI3K-Akt signaling in the OS tumors. Survival analyses showed that the expression levels of several genes such as RAMP1, CRIP1, CORT, CHST13, and DDX60L, miRNAs and lncRNAs were associated with survival of OS patients. Molecular subtyping using Cluster-Of-Clusters Analysis (COCA) for mRNA, lncRNA, and miRNA expression; DNA methylation; and mutation data from the TARGET dataset revealed two distinct molecular subtypes, each with a distinctive gene expression profile. Between the two subtypes, three upregulated genes, POP4, HEY1, CERKL, and seven downregulated genes, CEACAM1, ABLIM1, LTBP2, ISLR, LRRC32, PTPRF, and GPX3, associated with OS metastasis were found to be differentially regulated. Thus, the molecular subtyping results provide a strong basis for classification of OS patients that could be used to develop better prognostic treatment strategies.

Verification of Ferroptosis Subcluster-Associated Genes Related to Osteosarcoma and Exploration of Immune Targeted Therapy

Background

Despite tremendous advances in treating osteosarcoma (OS), the survival rates of patients have failed to improve dramatically over the past decades. Ferroptosis, a newly discovered iron-dependent type of regulated cell death, is implicated in tumors, and its features in OS remain unascertained. We designed to determine the involvement of ferroptosis subcluster-related modular genes in OS progression and prognosis.

Methods

The OS-related datasets retrieved from GEO and TARGET database were clustered for identifying molecular subclusters with different ferroptosis-related genes (FRGs) expression patterns. Weighted gene coexpression network analysis (WGCNA) was applied to identify modular genes from FRG subclusters. The least absolute shrinkage and selection operator (LASSO) algorithm and multivariable Cox regression analysis were adopted to develop the prognostic model. Potential mechanisms of development and prognosis in OS were explored by gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Then, a comprehensive analysis was conducted for immune checkpoint markers and assessment of predictive power to drug response. The protein expression levels of the three ferroptosis subcluster-related modular genes were verified by immunohistochemistry.

Results

Two independent subclusters presenting diverse expression profiles of FRGs were obtained, with significantly different survival states. Ferroptosis subcluster-related modular genes were screened with WGCNA, and the GESA results showed that ferroptosis subcluster-related modular genes could affect the cellular energy metabolism, thus influencing the development and prognosis of osteosarcoma. A prognostic model was established by incorporating three ferroptosis subcluster-related modular genes (LRRC1, ACO2, and CTNNBIP1) and a nomogram by integrating clinical features, and they were evaluated for the predictive power on OS prognosis. The 20 immune checkpoint-related genes confirmed the insensitivity to tumor immunotherapy in high-risk patients. IC50s of Axitinib and Cytarabine suggested a higher sensitivity to the targeted drug. Finally, the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and immunohistochemistry were consistent with bioinformatics analysis.

Conclusion

Ferroptosis are closely associated with the OS prognosis. The risk-scoring model incorporating three ferroptosis subcluster-related modular genes has shown outstanding advantages in predicting patient prognosis.

Metastatic Progression of Osteosarcomas: A Review of Current Knowledge of Environmental versus Oncogenic Drivers

Simple Summary

Osteosarcomas are heterogeneous bone tumors with complex genetic and chromosomic alterations. The numerous patients with metastatic osteosarcoma have a very poor prognosis, and only those who can have full surgical resection of the primary tumor and of all the macro metastasis can survive. Despite the recent improvements in prediction and early detection of metastasis, big efforts are still required to understand the specific mechanisms of osteosarcoma metastatic progression, in order to reveal novel therapeutic targets.

Abstract

Metastases of osteosarcomas are heterogeneous. They may grow simultaneously with the primary tumor, during treatment or shortly after, or a long time after the end of the treatment. They occur mainly in lungs but also in bone and various soft tissues. They can have the same histology as the primary tumor or show a shift towards a different differentiation path. However, the metastatic capacities of osteosarcoma cells can be predicted by gene and microRNA signatures. Despite the identification of numerous metastasis-promoting/predicting factors, there is no efficient therapeutic strategy to reduce the number of patients developing a metastatic disease or to cure these metastatic patients, except surgery. Indeed, these patients are generally resistant to the classical chemo- and to immuno-therapy. Hence, the knowledge of specific mechanisms should be extended to reveal novel therapeutic approaches. Recent studies that used DNA and RNA sequencing technologies highlighted complex relations between primary and secondary tumors. The reported results also supported a hierarchical organization of the tumor cell clones, suggesting that cancer stem cells are involved. Because of their chemoresistance, their plasticity, and their ability to modulate the immune environment, the osteosarcoma stem cells could be important players in the metastatic process.

Identification and Validation of Autophagy-Related Genes in Necrotizing Enterocolitis

BACKGROUND

Autophagy plays an essential role in the occurrence and progression of necrotizing enterocolitis (NEC). We intend to carry out the identification and validation of the probable autophagy-related genes of NEC via bioinformatics methods and experiment trials.

METHODS

The autophagy-related differentially expressed genes (arDEGs) of NEC were identified by analyzing the RNA sequencing data of the experiment neonatal mouse model and dataset GSE46619. Protein-protein interactions (PPIs), Gene Ontology (GO) enrichment analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used for the arDEGs. Then, co-expressed autophagy-related genes in two datasets were identified by Venn analysis and verified by qRT-PCR in experimental NEC.

RESULTS

Autophagy increased in experimental NEC and 47 arDEGs were identified in experimental NEC by RNA-sequencing. The PPI results proclaimed those genes interplayed with each other. The GO and KEGG enrichment results of arDEGs reported certain enriched pathways related to autophagy and macroautophagy. Furthermore, 22 arDEGs were identified in human NEC from dataset GSE46619. The GO and KEGG enrichment analysis of these genes showed similar enriched terms with the results of experimental NEC. Finally, HIF-1a, VEGFA, ITGA3, ITGA6, ITGB4, and NAMPT were identified as co-expressed autophagy-related genes by Venn analysis in human NEC from dataset GSE46619 and experimental NEC. The result of quantified real-time PCR (qRT-PCR) revealed that the expression levels of HIF-1a and ITGA3 were upregulated, while VEGFA and ITGB4 were downregulated in experimental NEC.

CONCLUSION

We identified 47 arDEGs in experimental NEC and 22 arDEGs in human NEC via bioinformatics analysis. HIF-1a, ITGA3, VEGFA, and ITGB4 may have effects on the progression of NEC through modulating autophagy.

Large-Scale Whole Genome Sequencing Study Reveals Genetic Architecture and Key Variants for Breast Muscle Weight in Native Chickens

Breast muscle weight (BrW) is one of the most important economic traits in chicken, and directional breeding for that results in both phenotypic and genetic changes. The Jingxing yellow chicken, including an original (without human-driven selection) line and a selected line (based on selection for increased intramuscular fat content), were used to dissect the genetic architecture and key variants associated with BrW. We detected 1069 high-impact single nucleotide polymorphisms (SNPs) with high conserved score and significant frequency difference between two lines. Based on the annotation result, the ECM-receptor interaction and fatty acid biosynthesis were enriched, and muscle-related genes, including MYOD1, were detected. By performing genome-wide association study for the BrW trait, we defined a major haplotype and two conserved SNPs that affected BrW. By integrated genomic and transcriptomic analysis, IGF2BP1 was identified as the crucial gene associated with BrW. In conclusion, these results offer a new insight into chicken directional selection and provide target genetic markers by which to improve chicken BrW.

A Novel Pyroptosis-Related Signature for Predicting Prognosis and Indicating Immune Microenvironment Features in Osteosarcoma

Osteosarcoma is a common malignant bone tumor with a propensity for drug resistance, recurrence, and metastasis. A growing number of studies have elucidated the dual role of pyroptosis in the development of cancer, which is a gasdermin-regulated novel inflammatory programmed cell death. However, the interaction between pyroptosis and the overall survival (OS) of osteosarcoma patients is poorly understood. This study aimed to construct a prognostic model based on pyroptosis-related genes to provide new insights into the prognosis of osteosarcoma patients. We identified 46 differentially expressed pyroptosis-associated genes between osteosarcoma tissues and normal control tissues. A total of six risk genes affecting the prognosis of osteosarcoma patients were screened to form a pyroptosis-related signature by univariate and LASSO regression analysis and verified using GSE21257 as a validation cohort. Combined with other clinical characteristics, including age, gender, and metastatic status, we found that the pyroptosis-related signature score, which we named “PRS-score,” was an independent prognostic factor for patients with osteosarcoma and that a low PRS-score indicated better OS and a lower risk of metastasis. The result of ssGSEA and ESTIMATE algorithms showed that a lower PRS-score indicated higher immune scores, higher levels of tumor infiltration by immune cells, more active immune function, and lower tumor purity. In summary, we developed and validated a pyroptosis-related signature for predicting the prognosis of osteosarcoma, which may contribute to early diagnosis and immunotherapy of osteosarcoma.