Similarity-based metric analysis approach for predicting osteogenic differentiation correlation coefficients and discovering the novel osteogenic-related gene FOXA1 in BMSCs.

Background: As a powerful tool, bioinformatics analysis is playing an increasingly important role in many fields. Osteogenic differentiation is a complex biological process involving the fine regulation of numerous genes and signaling pathways.

Method: Osteogenic differentiation-related genes are collected from the online databases. Then, we proposed two indexes Jaccard similarity and Sorensen-Dice similarity to measure the topological relevance of genes in the human PPI network. Furthermore, we selected three pathways involving osteoblast-related transcription factors, osteoblast differentiation, and RUNX2 regulation of osteoblast differentiation for investigation. Subsequently, we performed functional a enrichment analysis of these top-ranked genes to check whether these candidate genes identified by similarity-based metrics are enriched in some specific biological functions and states. we performed a permutation test to investigate the similarity score with four well-known osteogenic differentiation-related pathways including hedgehog signaling pathway, BMP signaling, ERK pathway, and Wnt signaling pathway to check whether these osteogenic differentiation-related pathways can be regulated by FOXA1. Lentiviral transfection was used to knockdown and overexpress gene FOXA1 in human bone mesenchymal stem cells (hBMSCs). Alkaline phosphatase (ALP) staining and Alizarin Red staining (ARS) were employed to investigate osteogenic differentiation of hBMSCs.

Result: After data collection, human PPI network involving 19,344 genes is included in our analysis. After simplifying, we used Jaccard and Sorensen-Dice similarity to identify osteogenic differentiation-related genes and integrated into a final similarity matrix. Furthermore, we calculated the sum of similarity scores with these osteogenic differentiation-related genes for each gene and found 337 osteogenic differentiation-related genes are involved in our analysis. We selected three pathways involving osteoblast-related transcription factors, osteoblast differentiation, and RUNX2 regulation of osteoblast differentiation for investigation and performed functional enrichment analysis of these top-ranked 50 genes. The results collectively demonstrate that these candidate genes can indeed capture osteogenic differentiation-related features of hBSMCs. According to the novel analyzing method, we found that these four pathways have significantly higher similarity with FOXA1 than random noise. Moreover, knockdown FOXA1 significantly increased the ALP activity and mineral deposits. Furthermore, overexpression of FOXA1 dramatically decreased the ALP activity and mineral deposits.

Conclusion: In summary, this study showed that FOXA1 is a novel significant osteogenic differentiation-related transcription factor. Moreover, our study has tightly integrated bioinformatics analysis with biological knowledge, and developed a novel method for analyzing the osteogenic differentiation regulatory network.