Abstract:Objective: Osteonecrosis of the femoral head (ONFH) is a common orthopedic condition that will prompt joint dysfunction, significantly impacting patients' quality of life. However, the specific pathogenic mechanisms underlying this disease remain elusive. The objective of this study is to examine the differentially expressed messenger RNAs (DE mRNAs) and key genes linked to ONFH, concurrently investigating the immune cell infiltration features in ONFH patients through the application of the CIBERSORT algorithm. Methods: Microarray was applied to scrutinize mRNA expression profiles in both ONFH patients and healthy controls, with data integration sourced from the GEO database. DE mRNAs were screened using the Limma method. The biological functions of DE mRNAs were explored through the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, Gene Ontology (GO) functional analysis, and Gene Set Enrichment Analysis (GSEA). Additionally, support vector machine-recursive feature elimination (SVM-RFE) and the least absolute shrinkage and selection operator (LASSO) were employed to discern diagnostic biomarkers associated with the disease. Receiver operating characteristic (ROC) analysis was utilized to assess the statistical performance of the feature genes. The validation of key genes was performed using qRT-PCR in bone tissues obtained from ONFH patients and healthy controls. Osteogenic differentiation of BMSC was then performed and detected by alkaline phosphatase staining (ALP) and qRT-PCR to verify the correlation between key genes and osteogenic differentiation. Finally, immune cell infiltration analysis was executed to evaluate immune cell dysregulation in ONFH, concurrently exploring the correlation between the infiltration of immune cells and key genes. Results: After consolidating the datasets, the Limma method revealed 107 DEGs, comprising 76 downregulated and 31 upregulated genes. Enrichment analysis revealed close associations of these DE mRNAs with functions such as cell migration, osteoblast differentiation, cartilage development and extracellular region. Machine learning algorithms further identified APOD, FBXO43 and LRP12 as key genes. ROC curves demonstrated the high diagnostic efficacy of these genes. The results of qRT-PCR showed that the expression levels of key genes were consistent with those of microarray analysis. In addition, the results of in vitro experiments showed that APOD was closely related to osteogenic differentiation of BMSC. Immune infiltration analysis suggested a close correlation between ONFH and imbalances in levels of Neutrophils, Monocytes, Macrophages M2, Dendritic cells activated and Dendritic cells resting. Conclusion: APOD is closely related to osteogenic differentiation of BMSCs and can be used as a diagnostic marker of ONFH. Immune cell infiltration significantly differs between controls and ONFH patients.

Identification of microtubule-associated biomarker using machine learning methods in osteonecrosis of the femoral head and osteosarcoma

Identification of hub genes and therapeutic drugs in osteonecrosis of the femoral head through integrated bioinformatics analysis and literature mining

Identification of a Programmed Cell Death-related Signature for Prediction of Osteosarcoma

Identification of Candidate Biomarkers Correlated With the Pathogenesis and Prognosis of Non-small Cell Lung Cancer via Integrated Bioinformatics Analysis

Comprehensive Analysis of Copper-Metabolism-related Genes about Prognosis and Immune Microenvironment in Osteosarcoma

Identification and Development of a Novel 4-Gene Immune-Related Signature to Predict Osteosarcoma Prognosis

Identification of Immune-Related Prognostic Genes and LncRNAs Biomarkers Associated with Osteosarcoma Microenvironment

Bioinformatic analysis of related immune cell infiltration and key genes in the progression of osteonecrosis of the femoral head

A ubiquitin-related gene signature for predicting prognosis and constructing molecular subtypes in osteosarcoma

A novel immune cell signature for predicting osteosarcoma prognosis and guiding therapy

A diagnostic signature developed based on the necroptosis-related genes and its association with immune infiltration in osteosarcoma

Transcriptome Analysis Identifies Novel Prognostic Genes in Osteosarcoma

Identification of Potential Therapeutic Targets and Immune Cell Infiltration Characteristics in Osteosarcoma Using Bioinformatics Strategy.

Exploring gene biomarkers and targeted drugs for ferroptosis and cuproptosis in osteosarcoma: A bioinformatic approach

TRIM68, PIKFYVE, and DYNLL2: The Possible Novel Autophagy- and Immunity-Associated Gene Biomarkers for Osteosarcoma Prognosis

Identification and Analysis of Three Hub Prognostic Genes Related to Osteosarcoma Metastasis

Identification of the Transcription Co-Factor–Related Gene Signature and Risk Score Model for Osteosarcoma

Construction of survival-related co-expression modules and identification of potential prognostic biomarkers of osteosarcoma using WGCNA

Comprehensive investigation of tumor immune microenvironment and prognostic biomarkers in osteosarcoma through integrated bulk and single-cell transcriptomic analysis

A novel molecular signature for predicting prognosis and immunotherapy response in osteosarcoma based on tumor-infiltrating cell marker genes

Screening candidate microRNA-mRNA network for predicting the response to chemoresistance in osteosarcoma by bioinformatics analysis.