Aging is an influential risk factor for progression of both degenerative and oncological diseases of the bone. Osteosarcoma, considered the most common primary mesenchymal tumor of the bone, is a worldwide disease with poor 5-year survival. This study investigated the role of aging-/senescence-induced genes (ASIGs) in contributing to osteosarcoma diagnosis, prognosis, and therapeutic agent prediction.
Therapeutically Applicable Research to Generate Effective Treatments (TARGET), Gene Expression Omnibus (GEO), and The Cancer Genome Atlas (TCGA) were used to collect relevant gene expression and clinical data of osteosarcoma and paracancerous tissues. Patients were clustered by consensus using prognosis-related ASIGs. ssGSEA, ESTIMATE, and TIMER were used to determine the tumor immune microenvironment (TIME) of subgroups. Functional analysis of differentially expressed genes between subgroups, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set variation analyses (GSVAs), was performed to clarify functional status. Prognostic risk models were constructed by univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression. SCISSOR was used to identify relevant cells in osteosarcoma single-cell data for different risk groups. The effect of immunotherapy was predicted based on TIDE scores and chemotherapy drug sensitivity using CTRP and PRISM.
Three molecular subgroups were identified based on prognostic differentially expressed ASIGs. Immunological infiltration levels of the three groups differed significantly. Based on GO and KEGG analyses, differentially expressed genes between the three subgroups mainly relate to immune and aging regulation pathways; GSVA showed substantial variations in multiple Hallmark pathways among the subgroups. The ASIG risk score built based on differentially expressed genes can predict patient survival and immune status. We also developed a nomogram graph to accurately predict prognosis in combination with clinical characteristics. The correlation between the immune activation profile of patients and the risk score is discussed. Through single-cell analysis of the tumor microenvironment, we identified distinct risk-group-associated cells with significant differences in immune signaling pathways. Immunotherapeutic efficacy and chemotherapeutic agent screening were evaluated based on risk score.
Aging-related prognostic genes can distinguish osteosarcoma molecular subgroups. Our novel aging-associated gene signature risk score can be used to predict the osteosarcoma immune landscape and prognosis. Moreover, the risk score correlates with the TIME and provides a reference for immunotherapy and chemotherapy in terms of osteosarcoma.