Identification and validation of autophagy-related genes during osteogenic differentiation of bone marrow mesenchymal stem cells

Osteogenic differentiation of bone marrow mesenchymal stem cells (BMSCs) is an essential stage in bone formation. Autophagy plays a pivotal role in the self-renewal potential and pluripotency of stem cells. This study aimed to explore the function of autophagy-related genes during osteogenic differentiation of BMSCs.

The differentially expressed autophagy-related genes (ARGs) were obtained from the GEO and HADb databases. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using R software. The PPI and hub gene mining networks were constructed using the STRING database and Cytoscape. Finally, the RT-qPCR was conducted to validate the expression level of ARGs in BMSCs.

Thirty-seven differentially expressed ARGs were finally obtained, including 12 upregulated and 25 downregulated genes. GO and KEGG enrichment analysis showed that most of these genes were enriched in apoptosis and autophagy. The PPI network revealed strong interactions between differentially expressed ARGs. The expression level of differentially expressed ARGs tested by RT-qPCR showed 6 upregulated ARGs, including FOXO1, MAP1LC3C, CTSB, FOXO3, CALCOCO2, FKBP1A, and 4 downregulated  ARGs, including MAPK8IP1, NRG1, VEGFA, and ITGA6 were consistent with the expression of high-throughput sequencing data.

We identified 37 ARGs during osteogenic differentiation using bioinformatics analysis. FOXO1, MAP1LC3C, CTSB, FOXO3, CALCOCO2, FKBP1A, MAPK8IP1, NRG1, VEGFA, and ITGA6 may regulate osteogenic differentiation of hBMSCs by involving autophagy pathway. This study provides new insight into the osteogenic differentiation of hBMSCs and may be available in developing therapeutic strategies for maxillofacial bone defects.