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

Background: Postmenopausal osteoporosis is a systemic metabolic disorder typified by an imbalance in bone turnover, where bone resorption supersedes bone formation. This imbalance primarily arises from a decline in bone mass induced by estrogen deficiency, and an elevated risk of fractures resulting from degradation of bone microstructure. Despite recognizing these changes, the precise causative factors and potential molecular pathways remain elusive. In this study, we aimed to identify differentially expressed genes (DEGs), associated pathways, and the role of immune infiltration in osteoporosis, leveraging an integrated bioinformatics approach to shed light on potential underlying molecular mechanisms. Methods: We retrieved the expression profiles of GSE230665 from the Gene Expression Omnibus database, comprising 15 femur samples, including 12 postmenopausal osteoporosis samples and 3 normal controls. From the aggregated microarray datasets, we derived differentially expressed genes (DEGs) for further bioinformatics analysis. We used WGCNA, analyzed DEGs, PPI, and conducted GO analysis to identify pivotal genes. We then used the CIBERSORT method to explore the degree of immune cell infiltration within femur specimens affected by postmenopausal osteoporosis. To probe into the relationship between pivotal genes and infiltrating immune cells, we conducted correlation analysis. Results: We identified a total of 12,204 DEGs. Among these, 12,157 were up-regulated, and 47 were down-regulated. GO and KEGG pathway analyses indicated that these DEGs predominantly targeted cellular protein localization activity and associated signaling pathways. The protein-protein interaction network highlighted four central hub-genes: RPL31, RPL34, EEF1G, and BPTF. Principal component analysis indicated a positive correlation between the expression of these genes and resting NK cells (as per CIBERSORT). In contrast, the expression of RPL31, RPL34, and EEF1G showed a negative correlation with T cells (gamma delta per CIBERSORT). Conclusions: Immune infiltration plays a role in the development of osteoporosis.