In this study, we screened the hub genes of osteoporosis by bioinformatics, and further analyzed their relationship with osteoporosis and their mechanism of action.

Firstly, the gene expression profile data set was downloaded from the public gene expression database, and the differentially expressed genes were screened by R software, and their functions and pathways were analyzed. Then, the online tool String was used to construct the protein interaction network, and the software cytoscape was imported to screen the core genes and construct the clustering module.

A total of 1334 differentially expressed genes were screened, including 722 up-regulated genes and 612down-regulated genes. GO analysis shows that the functions are mainly concentrated in extracellular matrix structural constituent, signal receptor activator activity, transmembrane transporter binding and cytokine binding. The enrichment of KEGG pathway shows that the differential genes are mainly involved in PI3K-Akt signaling pathway, MAPK signaling pathway, Rap1 signaling pathway and Ras signaling pathway. Ten key genes, namely AKT1, EGF, VEGFA, PROM1, TP53, NES, CD21, SNAL1, FGF13, LIF, and one clustering module were selected according to protein interaction network.

The functions and roles of key genes and clustering modules and their possible relationship with osteoporosis were screened and analyzed, which provided new ideas for revealing the potential molecular mechanism and drug targets of osteoporosis.