Screening of key candidate genes and pathways for osteocytes involved in the differential response to different types of mechanical stimulation using a bioinformatics analysis

This study aimed to predict the key genes and pathways that are activated when different types of mechanical loading are applied to osteocytes. mRNA expression datasets (series number of GSE62128 and GSE42874) were obtained from Gene Expression Omnibus database (GEO). High gravity-treated osteocytic MLO-Y4 cell-line samples from GSE62128 (Set1), and fluid flow-treated MLO-Y4 samples from GSE42874 (Set2) were employed. After identifying the differentially expressed genes (DEGs), functional enrichment was performed. The common DEGs between Set1 and Set2 were considered as key DEGs, then a protein–protein interaction (PPI) network was constructed using the minimal nodes from all of the DEGs in Set1 and Set2, which linked most of the key DEGs. Several open source software programs were employed to process and analyze the original data. The bioinformatic results and the biological meaning were validated by in vitro experiments. High gravity and fluid flow induced opposite expression trends in the key DEGs. The hypoxia-related biological process and signaling pathway were the common functional enrichment terms among the DEGs from Set1, Set2 and the PPI network. The expression of almost all the key DEGs (Pdk1, Ccng2, Eno2, Egln1, Higd1a, Slc5a3 and Mxi1) were mechano-sensitive. Eno2 was identified as the hub gene in the PPI network. Eno2 knockdown results in expression changes of some other key DEGs (Pdk1, Mxi1 and Higd1a). Our findings indicated that the hypoxia response might have an important role in the differential responses of osteocytes to the different types of mechanical force.