Osteoporosis is a common bone disease characterized by loss of bone mass, reduced bone strength, and deterioration of bone microstructure. ROS-induced oxidative stress plays an important role in osteoporosis. However, the biomarkers and molecular mechanisms of oxidative stress are still unclear. We obtained the datasets from the Gene Expression Omnibus (GEO) database, and performed differential analysis, Venn analysis, and weighted correlation network analysis (WGCNA) analysis out the hub genes. Then, the correlation between inflammatory factors and hub genes was analyzed, and a Mendelian randomization (MR) analysis was performed on cytokines and osteoporosis outcomes. In addition, "CIBERSORT" was used to analyze the infiltration of immune cells and single-cell RNA-seq data was used to analyze the expression distribution of hub genes and cell-cell communications. Finally, we collected human blood samples for RT-qPCR and Elisa experiments, the miRNA-mRNA network was constructed using the miRBase database, the 3D structure was predicted using the RNAfold, Vfold3D database, and the drug sensitivity analysis was performed using the RNAactDrug database. We obtained three differentially expressed genes associated with oxidative stress: DBH, TAF15, and STAT4 by differential, WGCNA clustering, and Venn screening analyses, and further analyzed the correlation of these 3 genes with inflammatory factors and immune cell infiltration and found that STAT4 was significantly and positively correlated with IL-2. Single-cell data analysis showed that the STAT4 gene was highly expressed mainly in dendritic cells and monocytes. In addition, the results of RT-qPCR and Elisa experiments verified that the expression of STAT4 was consistent with the previous analysis, and a significant causal relationship between IL-2 and STAT4 SNPs and osteoporosis was found by Mendelian randomization. Finally, through miRNA-mRNA network and drug sensitivity analysis, we analyzed to get Palbociclib/miR-141-3p/STAT4 axis, which can be used for the prevention and treatment of osteoporosis. In this study, we proposed the Palbociclib/miR-141-3p/STAT4 axis for the first time and provided new insights into the mechanism of oxidative stress in osteoporosis.

Osteoporosis (OP) is a chronic bone metabolic disease and a serious global public health problem. Several studies have shown that mitophagy plays an important role in bone metabolism disorders; however, its role in osteoporosis remains unclear. The Gene Expression Omnibus (GEO) database was used to download GSE56815, a dataset containing low and high BMD, and differentially expressed genes (DEGs) were analyzed. Mitochondrial autophagy-related genes (MRG) were downloaded from the existing literature, and highly correlated MRG were screened by bioinformatics methods. The results from both were taken as differentially expressed (DE)-MRG, and Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed. Protein-protein interaction network (PPI) analysis, support vector machine recursive feature elimination (SVM-RFE), and Boruta method were used to identify DE-MRG. A receiver operating characteristic curve (ROC) was drawn, a nomogram model was constructed to determine its diagnostic value, and a variety of bioinformatics methods were used to verify the relationship between these related genes and OP, including GO and KEGG analysis, IP pathway analysis, and single-sample Gene Set Enrichment Analysis (ssGSEA). In addition, a hub gene-related network was constructed and potential drugs for the treatment of OP were predicted. Finally, the specific genes were verified by real-time quantitative polymerase chain reaction (RT-qPCR). In total, 548 DEGs were identified in the GSE56815 dataset. The weighted gene co-expression network analysis(WGCNA) identified 2291 key module genes, and 91 DE-MRG were obtained by combining the two. The PPI network revealed that the target gene for AKT1 interacted with most proteins. Three MRG (NELFB, SFSWAP, and MAP3K3) were identified as hub genes, with areas under the curve (AUC) 0.75, 0.71, and 0.70, respectively. The nomogram model has high diagnostic value. GO and KEGG analysis showed that ribosome pathway and cellular ribosome pathway may be the pathways regulating the progression of OP. IPA showed that MAP3K3 was associated with six pathways, including GNRH Signaling. The ssGSEA indicated that NELFB was highly correlated with iDCs (cor = -0.390, < 0.001). The regulatory network showed a complex relationship between miRNA, transcription factor(TF) and hub genes. In addition, 4 drugs such as vinclozolin were predicted to be potential therapeutic drugs for OP. In RT-qPCR verification, the hub gene NELFB was consistent with the results of bioinformatics analysis. Mitophagy plays an important role in the development of osteoporosis. The identification of three mitophagy-related genes may contribute to the early diagnosis, mechanism research and treatment of OP.