Investigation of transcriptome mechanism associated with osteoporosis explored by microarray analysis

Microarray data of osteoporosis (OP) were analyzed based on prediction of transcription factors (TFs) or their targets as well as influences of TFs or TF network to uncover key TFs in OP. The microarray data E-GEOD-35956 was downloaded from the GPL570 platform. Differentially expressed genes (DEGs) with logarithm of fold change (|logFC|) >2 and P-value <0.05 were identified between OP samples and normal controls. TF genes were screened from the DEGs based on ITFP, Marbach 2016, TRRUST databases. TF targets were enriched from DEGs using Fisher's exact test. TF targets were selected based on their impact factors. TF targets were chosen from TF network analysis. Finally, key TFs were identified by based on TFs coverage. A total of 300 DEGs were obtained. There were no TF genes screened from the DEGs. In total 165, 87 and 178 TF targets were screened from DEGs respectively based on Fisher's exact test, influence of TFs or TF network analysis. According to the optimal TF set with TFs having maximum coverage of DEGs, 178 TF targets was the most. Thus, the optimal sets of TFs were FOXO1, KLF16, RXRA, RARA, HNF4A, CEBPB, ESR1, SOX8, ZNF219, and SP1. Altogether, these results suggested identified crucial TFs in OP might play a significant role in OP development, showing these key TFs probably would aid in unveiling the underlying molecular mechanisms and may be therapeutic targets, diagnostic or prognostic biomarkers for OP.