Bioinformatics analysis of the key genes in osteosarcoma metastasis and immune invasion

Machine learning survival prediction using tumor lipid metabolism genes for osteosarcoma

Osteosarcoma is a primary malignant tumor that commonly affects children and adolescents, with a poor prognosis. The existence of tumor heterogeneity leads to different molecular subtypes and survival outcomes. Recently, lipid metabolism has been identified as a critical characteristic of cancer. Therefore, our study aims to identify osteosarcoma's lipid metabolism molecular subtype and develop a signature for survival outcome prediction. Four multicenter cohorts-TARGET-OS, GSE21257, GSE39058, and GSE16091-were amalgamated into a unified Meta-Cohort. Through consensus clustering, novel molecular subtypes within Meta-Cohort patients were delineated. Subsequent feature selection processes, encompassing analyses of differentially expressed genes between subtypes, univariate Cox analysis, and StepAIC, were employed to pinpoint biomarkers related to lipid metabolism in TARGET-OS. We selected the most effective algorithm for constructing a Lipid Metabolism-Related Signature (LMRS) by utilizing four machine-learning algorithms reconfigured into ten unique combinations. This selection was based on achieving the highest concordance index (C-index) in the test cohort of GSE21257, GSE39058, and GSE16091. We identified two distinct lipid metabolism molecular subtypes in osteosarcoma patients, C1 and C2, with significantly different survival rates. C1 is characterized by increased cholesterol, fatty acid synthesis, and ketone metabolism. In contrast, C2 focuses on steroid hormone biosynthesis, arachidonic acid, and glycerolipid and linoleic acid metabolism. Feature selection in the TARGET-OS identified 12 lipid metabolism genes, leading to a model predicting osteosarcoma patient survival. The LMRS, based on the 12 identified genes, consistently accurately predicted prognosis across TARGET-OS, testing cohorts, and Meta-Cohort. Incorporating 12 published signatures, LMRS showed robust and significantly superior predictive capability. Our results offer a promising tool to enhance the clinical management of osteosarcoma, potentially leading to improved clinical outcomes.

Using machine learning methods to study the tumour microenvironment and its biomarkers in osteosarcoma metastasis

Background

The long-term prognosis for patients with osteosarcoma (OS) metastasis remains unfavourable, highlighting the urgent need for research that explores potential biomarkers using innovative methodologies.

Methods

This study explored potential biomarkers for OS metastasis by analysing data from the Cancer Genome Atlas Program (TCGA) and Gene Expression Omnibus (GEO) databases. The synthetic minority oversampling technique (SMOTE) was employed to tackle class imbalances, while genes were selected using four feature selection algorithms (Monte Carlo feature selection [MCFS], Borota, minimum-redundancy maximum-relevance [mRMR], and light gradient-boosting machine [LightGBM]) based on the gene expression matrix. Four machine learning (ML) algorithms (support vector machine [SVM], extreme gradient boosting [XGBoost], random forest [RF], and k-nearest neighbours [kNN]) were utilized to determine the optimal number of genes for building the model. Interpretable machine learning (IML) was applied to construct prediction networks, revealing potential relationships among the selected genes. Additionally, enrichment analysis, survival analysis, and immune infiltration were performed on the featured genes.

Results

In DS1, DS2, and DS3, the IML algorithm identified 53, 45, and 46 features, respectively. Using the merged gene set, we obtained a total of 79 interpretable prediction rules for OS metastasis. We subsequently conducted an in-depth investigation on 39 crucial molecules associated with predicting OS metastasis, elucidating their roles within the tumour microenvironment. Importantly, we found that certain genes act as both predictors and differentially expressed genes. Finally, our study unveiled statistically significant differences in survival between the high and low expression groups of TRIP4, S100A9, SELL and SLC11A1, and there was a certain correlation between these genes and 22 various immune cells.

Conclusions

The biomarkers discovered in this study hold significant implications for personalized therapies, potentially enhancing the clinical prognosis of patients with OS.

Osteosarcoma-Associated Immune Genes as Potential Immunotherapy and Prognosis Biomarkers

Immune-modulating therapies exhibit abundant promise compared to traditional treatment for osteosarcoma. We aim to establish an immune-related prognostic signatures in osteosarcoma. We identified the differentially expressed genes in osteosarcoma compared with normal controls using the GEO dataset. The intersection with immune-related genes was considered as differentially expressed immune genes. Potential prognosis-related differentially expressed genes were first analyzed with the multifactor Cox regression and then the step function performed the iteration. The best model was finally chosen as the immunological risk score signature model. And finally, we evaluated the correlation of genes in the prognostic model with immune cells, common immune checkpoints, and immune checkpoint blockade responses. We identified 1527 significantly upregulated and 2407 significantly downregulated genes in osteosarcoma compared to normal samples. In the 258 differentially expressed immune genes, 20 genes with independent prognostic significance were included in the step function. Finally, we constructed a prognostic signature for overall survival based on five immune genes (JAG2, PLXNB1, CMKLR1, NUDT6, and PSMC4) in osteosarcoma. These five genes are closely related to immune infiltration. Osteosarcoma with high JAG2 expression or low CMKLR1 expression may be associated with better immune checkpoint blockade response. JAG2 overexpression or CMKLR1 inhibition induced sensitivity to PD-1 antibody in osteosarcoma cells. We constructed a prognosis prediction signature containing five immune-related genes (JAG2, PLXNB1, CMKLR1, NUDT6, and PSMC4) with a significant prognostic value in osteosarcoma. Significant differences in immune infiltration and immunotherapy responses were identified between groups with different levels of these genes.

The immune cell infiltration-associated molecular subtypes and gene signature predict prognosis for osteosarcoma patients

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.

Single-cell RNA sequencing reveals the immune microenvironment landscape of osteosarcoma before and after chemotherapy

Chemotherapy, a primary treatment for osteosarcoma (OS), has limited knowledge regarding its impact on tumor immune microenvironment (TIME). Here, tissues from 6 chemotherapy-naive OS patients underwent single-cell RNA sequencing (scRNA-seq) and were analyzed alongside public dataset (GSE152048) containing 7 post-chemotherapy OS tissues. CD45+ (PTPRC+) cells were used for cell clustering and annotation. Changes in immune cell composition pre- and post-chemotherapy were characterized. Totally, 28,636 high-quality CD45+ (PTPRC+) cells were extracted. Following chemotherapy, the proportions of regulatory T cells (Tregs) and activated CD8 T cells decreased, while CD8 effector T cells increased. GO analysis indicated that differentially expressed genes (DEGs) in T cells were associated with cell activation, adaptive immune response, and immune response to tumor cells. Furthermore, the proportions of plasma cells increased, while naive B cells decreased. B cell surface receptors expression was upregulated, and GO analysis revealed DEGs of B cells were mainly enriched in B cell-mediated immunity and B cell activation. Moreover, M2 polarization of macrophages was suppressed post-chemotherapy. Overall, this study elucidates chemotherapy remodels the OS TIME landscape, triggering immune heterogeneity and enhancing anti-tumor properties.

Plasma protein changes reflect colorectal cancer development and associated inflammation

Introduction

Colorectal cancer (CRC) is the third most common malignancy and the second leading cause of death worldwide. Efficient non-invasive blood-based biomarkers for CRC early detection and prognosis are urgently needed.

Methods

To identify novel potential plasma biomarkers, we applied a proximity extension assay (PEA), an antibody-based proteomics strategy to quantify the abundance of plasma proteins in CRC development and cancer-associated inflammation from few μL of plasma sample.

Results

Among the 690 quantified proteins, levels of 202 plasma proteins were significantly changed in CRC patients compared to age-and-sex-matched healthy subjects. We identified novel protein changes involved in Th17 activity, oncogenic pathways, and cancer-related inflammation with potential implications in the CRC diagnosis. Moreover, the interferon γ (IFNG), interleukin (IL) 32, and IL17C were identified as associated with the early stages of CRC, whereas lysophosphatidic acid phosphatase type 6 (ACP6), Fms-related tyrosine kinase 4 (FLT4), and MANSC domain-containing protein 1 (MANSC1) were correlated with the late-stages of CRC.

Discussion

Further study to characterize the newly identified plasma protein changes from larger cohorts will facilitate the identification of potential novel diagnostic, prognostic biomarkers for CRC.

Construction of Osteosarcoma Diagnosis Model by Random Forest and Artificial Neural Network

Osteosarcoma accounts for 28% of primary bone malignancies in adults and up to 56% in children and adolescents (<20 years). However, early diagnosis and treatment are still inadequate, and new improvements are still needed. Missed diagnoses exist due to fewer traditional diagnostic methods, and clinical symptoms are often already present before diagnosis. This study aimed to develop novel and efficient predictive models for the diagnosis of osteosarcoma and to identify potential targets for exploring osteosarcoma markers. First, osteosarcoma and normal tissue expression microarray datasets were downloaded from the Gene Expression Omnibus (GEO). Then we screened the differentially expressed genes (DEGs) in the osteosarcoma and normal groups in the training group. Next, in order to explore the biologically relevant role of DEGs, Metascape and enrichment analyses were also performed on DEGs. The “randomForest” and “neuralnet” packages in R software were used to select representative genes and construct diagnostic models for osteosarcoma. The next step is to validate the model of the artificial neural network. Then, we performed an immune infiltration analysis by using the training set data. Finally, we constructed a prognostic model using representative genes for prognostic analysis. The copy number of osteosarcoma was also analyzed. A random forest classifier identified nine representative genes (ANK1, TGFBR3, RSF21, HSPB8, ITGA7, RHD, AASS, GREM2, NFASC). HSPB8, RHD, AASS, and NFASC were genes we identified that have not been previously reported to be associated with osteosarcoma. The osteosarcoma diagnostic model we constructed has good performance with areas under the curves (AUCs) of 1 and 0.987 in the training and validation groups, respectively. This study opens new horizons for the early diagnosis of osteosarcoma and provides representative markers for the future treatment of osteosarcoma. This is the first study to pioneer the establishment of a genetic diagnosis model for osteosarcoma and advance the development of osteosarcoma diagnosis and treatment.

Screening and identification of susceptibility genes for osteosarcoma based on bioinformatics analysis

Background

Gene play an important role in malignant tumors. However, there is still insufficient research on genetic variations in osteosarcoma (OS) patients. Therefore, we aimed to analyze the gene expression profile of OS using bioinformatics and to explore the pathogenesis of OS at the molecular level.

Methods

The gene chip dataset of OS samples was downloaded from the Gene Expression Omnibus (GEO) database for screening differentially expressed genes (DEGs). The R language clusterProfiler software package was used to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs. The central node proteins of the protein interaction network were analyzed by the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, Cytoscape, and its plug-ins cytoHubba and NetworkAnalyzer to find the key genes.

Results

A total of 631 DEGs were obtained, including 362 upregulated genes and 269 downregulated genes. DEGs were mainly involved in the regulation of leukocyte chemotaxis and migration, vascular development, and other biological processes (BPs); mediation of receptor ligand activity, growth factor binding, growth factor activity, integrin binding, and other molecular functions (MFs); and were enriched in the extracellular matrix (ECM).

Conclusions

DEGs in the ECM and growth factors play a key role in the development of OS. The leukocyte transendothelial migration pathway and the PI3K-AKT pathway are closely related to OS, and the related molecular mechanism is worthy of further study.