Abstract:Backgrounds: As a systemic skeletal dysfunction, osteoporosis (OP) is characterized by low bone mass and bone microarchitectural damage. The global incidences of OP are high. Methods: Data were retrieved from databases like Gene Expression Omnibus (GEO), GeneCards, Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), Gene Expression Profiling Interactive Analysis (GEPIA2), and other databases. R software (version 4.1.1) was used to identify differentially expressed genes (DEGs) and perform functional analysis. The Least Absolute Shrinkage and Selection Operator (LASSO) logistic regression and random forest algorithm were combined and used for screening diagnostic markers for OP. The diagnostic value was assessed by the receiver operating characteristic (ROC) curve. Molecular signature subtypes were identified using a consensus clustering approach, and prognostic analysis was performed. The level of immune cell infiltration was assessed by the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm. The hub gene was identified using the CytoHubba algorithm. Real-time fluorescence quantitative PCR (RT-qPCR) was performed on the plasma of osteoporosis patients and control samples. The interaction network was constructed between the hub genes and miRNAs, transcription factors, RNA binding proteins, and drugs. Results: A total of 40 DEGs, eight OP-related differential genes, six OP diagnostic marker genes, four OP key diagnostic marker genes, and ten hub genes (TNF, RARRES2, FLNA, STXBP2, EGR2, MAP4K2, NFKBIA, JUNB, SPI1, CTSD) were identified. RT-qPCR results revealed a total of eight genes had significant differential expression between osteoporosis patients and control samples. Enrichment analysis showed these genes were mainly related to MAPK signaling pathways, TNF signaling pathway, apoptosis, and Salmonella infection. RT-qPCR also revealed that the MAPK signaling pathway (p38, TRAF6) and NF-kappa B signaling pathway (c-FLIP, MIP1β) were significantly different between osteoporosis patients and control samples. The analysis of immune cell infiltration revealed that monocytes, activated CD4 memory T cells, and memory and naïve B cells may be related to the occurrence and development of OP. Conclusions: We identified six novel OP diagnostic marker genes and ten OP-hub genes. These genes can be used to improve the prognostic of OP and to identify potential relationships between the immune microenvironment and OP. Our research will provide insights into the potential therapeutic targets and pathogenesis of osteoporosis.

Identification of 12 hub genes associated to the pathogenesis of osteoporosis based on microarray and single-cell RNA sequencing data

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

What do osteoporosis and osteoarthritis have in common? An integrated study of overlapping differentially expressed genes in bone mesenchymal stem cells of osteoporosis and osteoarthritis

Identification of potential pathogenic genes related to osteoporosis and osteoarthritis

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

Single-cell RNA sequencing analysis dissected the osteo-immunology microenvironment and revealed key regulators in osteoporosis

Study on Key Genes and Regulatory Networks Associated with Osteoporosis by Microarray Technology

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

Integration of single‐cell and RNA‐seq data to explore the role of focal adhesion‐related genes in osteoporosis

Identifying Crucial Biomarkers in Osteoporosis and Ulcerative Colitis Through Bioinformatics Analysis of Co-expressed Genes

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

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

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

Revealing the key role of cuproptosis in osteoporosis via the bioinformatic analysis and experimental validation of cuproptosis-related genes

Identification of Biomarkers Modulating Osteogenic Differentiation With Immune and Inflammation Landscapes in Mesenchymal Stem Cells: A Bioinformatics-based Comprehensive Study

RETRACTED ARTICLE: Key Genes Associated with Osteoporosis Revealed by Genome Wide Gene Expression Analysis

Identification of Key Osteoporosis Genes Through Comparative Analysis of Men's and Women's Osteoblast Transcriptomes

Transcriptomic Data Identified Key Transcription Factors for Osteoporosis in Caucasian Women

Transcriptome Sequencing Identifies Abnormal Lncrnas and Mrnas and Reveals Potentially Hub Immune-Related Mrna in Osteoporosis with Vertebral Fracture.

Unveiling the unexplored novel signatures for osteoporosis via a detailed bioinformatics and molecular experiments based approach

Employing single-cell RNA sequencing coupled with an array of bioinformatics approaches to ascertain the shared genetic characteristics between osteoporosis and obesity