Differentially Expressed Genes Reveal the Biomarkers and Molecular Mechanism of Osteonecrosis

Osteonecrosis is one of the most refractory orthopedic diseases, which seriously threatens the health of old patients. High-throughput sequencing (HTS) and microarray analysis have confirmed as an effective way for investigating the pathological mechanism of disease. In this study, GSE7716, GSE74089, and GSE123568 were obtained from Gene Expression Omnibus (GEO) database and used to identify differentially expressed genes (DEGs) by R language. Subsequently, the DEGs were analyzed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Moreover, the protein-protein interaction (PPI) network of DEGs was analyzed by STRING database and Cytoscape. The results showed that 318 downregulated genes and 58 upregulated genes were observed in GSE7116; 690 downregulated genes and 1148 upregulated genes were screened from 34183 genes in GSE74089. The DEGs involved in progression of osteonecrosis involved inflammation, immunological rejection, and bacterial infection-related pathways. The GO enrichment showed that osteonecrosis was related with extracellular matrix, external encapsulating structure organization, skeletal system development, immune response activity, cell apoptosis, mononuclear cell differentiation, and serine/threonine kinase activity. Moreover, PPI network showed that the progression of osteonecrosis of the femoral head was related with CCND1, CDH1, ESR1, SPP1, LOX, JUN, ITGA, ABL1, and VEGF, and osteonecrosis of the jaw is related with ACTB, CXCR4, PTPRC, IL1B, CXCL8, TNF, JUN, PTGS2, FOS, and RHOA. In conclusion, this study identified the hub factors and pathways which might play important roles in progression of osteonecrosis and could be used as potential biomarkers for diagnosis and treatment of osteonecrosis.