Abstract:Postmenopausal osteoporosis (PMOP) is a common metabolic inflammatory disease. In conditions of estrogen deficiency, chronic activation of the immune system leads to a hypo-inflammatory phenotype and alterations in its cytokine and immune cell profile, although immune cells play an important role in the pathology of osteoporosis, studies on this have been rare. Therefore, it is important to investigate the role of immune cell-related genes in PMOP. PMOP-related datasets were downloaded from the Gene Expression Omnibus database. Immune cells scores between high bone mineral density (BMD) and low BMD samples were assessed based on the single sample gene set enrichment analysis method. Subsequently, weighted gene co-expression network analysis was performed to identify modules highly associated with immune cells and obtain module genes. Differential analysis between high BMD and low BMD was also performed to obtain differentially expressed genes. Module genes are intersected with differentially expressed genes to obtain candidate genes, and functional enrichment analysis was performed. Machine learning methods were used to filter out the signature genes. The receiver operating characteristic (ROC) curves of the signature genes and the nomogram were plotted to determine whether the signature genes can be used as a molecular marker. Gene set enrichment analysis was also performed to explore the potential mechanism of the signature genes. Finally, RNA expression of signature genes was validated in blood samples from PMOP patients and normal control by real-time quantitative polymerase chain reaction. Our study of PMOP patients identified differences in immune cells (activated dendritic cell, CD56 bright natural killer cell, Central memory CD4 T cell, Effector memory CD4 T cell, Mast cell, Natural killer T cell, T follicular helper cell, Type 1 T-helper cell, and Type 17 T-helper cell) between high and low BMD patients. We obtained a total of 73 candidate genes based on modular genes and differential genes, and obtained 5 signature genes by least absolute shrinkage and selection operator and random forest model screening. ROC, principal component analysis, and t-distributed stochastic neighbor embedding down scaling analysis revealed that the 5 signature genes had good discriminatory ability between high and low BMD samples. A logistic regression model was constructed based on 5 signature genes, and both ROC and column line plots indicated that the model accuracy and applicability were good. Five signature genes were found to be associated with proteasome, mitochondria, and lysosome by gene set enrichment analysis. The real-time quantitative polymerase chain reaction results showed that the expression of the signature genes was significantly different between the 2 groups. HIST1H2AG, PYGM, NCKAP1, POMP, and LYPLA1 might play key roles in PMOP and be served as the biomarkers of PMOP.

Identification of key immune genes of osteoporosis based on bioinformatics and machine learning

Multi-omics Analysis to Identify Key Immune Genes for Osteoporosis based on Machine Learning and Single-cell Analysis

Identification of mitophagy-related biomarkers in human osteoporosis based on a machine learning model

Screening and Analysis of Core Genes for Osteoporosis Based on Bioinformatics Analysis and Machine Learning Algorithms

Prognostic analysis and validation of diagnostic marker genes in patients with osteoporosis

Identification of differentially expressed genes, signaling pathways and immune infiltration in postmenopausal osteoporosis by integrated bioinformatics analysis

Identification of biomarkers associated with diagnosis of postmenopausal osteoporosis patients based on bioinformatics and machine learning

Screening of Important Markers in Peripheral Blood Mononuclear Cells to Predict Female Osteoporosis Risk Using LASSO Regression Algorithm and SVM Method

Identification of Potential Therapeutic Targets and Molecular Regulatory Mechanisms for Osteoporosis by Bioinformatics Methods

Identification of shared gene signatures and pathways for diagnosing osteoporosis with sarcopenia through integrated bioinformatics analysis and machine learning

Machine learning analysis of gene expression profile reveals a novel diagnostic signature for osteoporosis

Identification of Potential Osteoporosis miRNA Biomarkers Using Bioinformatics Approaches

A Machine Learning Framework for Screening Plasma Cell-Associated Feature Genes to Estimate Osteoporosis Risk and Treatment Vulnerability

Identification of combined biomarkers for predicting the risk of osteoporosis using machine learning

Consensus Gene Modules Related to Levels of Bone Mineral Density (BMD) among Smokers and Nonsmokers

Analysis and validation of biomarkers of immune cell-related genes in postmenopausal osteoporosis: An observational study

Identification of hub genes associated with osteoporosis development by comprehensive bioinformatics analysis

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

Identification of Biomarkers Associated with Diagnosis of Osteoarthritis Patients Based on Bioinformatics and Machine Learning

Identifying Key Genes and Functionally Enriched Pathways in Osteoporotic Patients by Weighted Gene Co-Expression Network Analysis

Integrated Single-cell RNA-seq and Bulk RNA-seq Identify Diagnostic Biomarkers for Postmenopausal Osteoporosis