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.

Association Between SEMA3A Signaling Pathway Genes and BMD/OP Risk: an Epidemiological and Experimental Study

SERUM C-C MOTIF LIGAND 11/EOTAXIN-1 MAY SERVE AS A CANDIDATE BIOMARKER FOR POSTMENOPAUSAL OSTEOPOROSIS

Evaluation of Relationship Between SPON1 Gene and Genetic Susceptibility of Postmenopausal Osteoporosis

Serum Sema3A is in a weak positive association with bone formation marker osteocalcin but not related to bone mineral densities in postmenopausal women.

Association Between NF-κB Signal Pathway-Related Gene Polymorphisms and Response to Alendronate Treatment in Postmenopausal Chinese Women with Low Bone Mineral Density

Estrogen receptor α and vitamin D receptor gene polymorphisms and bone mineral density: association study of healthy pre- and postmenopausal Chinese women

SAA1 Gene Polymorphisms in Osteoporosis Patients

Endochondral Ossification Pathway Genes and Postmenopausal Osteoporosis: Association and Specific Allele Related Serum Bone Sialoprotein Levels in Han Chinese.

The Influence of Lys3Asn Polymorphism in the Osteoprotegerin Gene on Bone Mineral Density in Chinese Postmenopausal Women

Evaluation of common variants in the CNR2 gene and its interaction with abdominal obesity for osteoporosis susceptibility in Chinese post-menopausal females

In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis.

Relationship of COL9A1 and SOX9 Genes with Genetic Susceptibility of Postmenopausal Osteoporosis

Influences of multiple gene interactions on bone mineral density and osteoporotic fractures in postmenopausal women

Integrating Genome- Wide Association Study with Regulatory SNP Annotations Identified Novel Candidate Genes for Osteoporosis

Relationship between polymorphisms and mutations at rs7125942 and rs3736228 of LRP5 gene and bone metabolism in postmenopausal women

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

Analysis of Recently Identified Osteoporosis Susceptibility Genes in Han Chinese Women.

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

Identification of a role for the ARHGEF3 gene in postmenopausal osteoporosis

Integrative Analysis of LGR5/6 Gene Variants, Gut Microbiota Composition and Osteoporosis Risk in Elderly Population

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