Genomic determinants of gene regulation by 1,25-dihydroxyvitamin D3 during osteoblast-lineage cell differentiation

The biological effects of 1α,25-dihydroxyvitamin D3 (1,25 (OH)2D3) on osteoblast differentiation and function differ significantly depending upon the cellular state of maturation. To explore this phenomenon mechanistically, we examined the impact of 1,25(OH)2D3 on the transcriptomes of both pre-osteoblastic (POBs) and differentiated osteoblastic (OBs) MC3T3-E1 cells, and assessed localization of the vitamin D receptor (VDR) at sites of action on a genome-scale using ChIP sequence analysis. We observed that the 1,25(OH)2D3-induced transcriptomes of POBs and OBs were quantitatively and qualitatively different, supporting not only the altered biology observed but the potential for a change in VDR interaction at the genome as well. This idea was confirmed through discovery that VDR cistromes in POBs and OBs were also strikingly different. Depletion of VDR-binding sites in OBs, due in part to reduced VDR expression, was the likely cause of the loss of VDR-target gene interaction. Continued novel regulation by 1,25(OH)2D3, however, suggested that factors in addition to the VDR might also be involved. Accordingly, we show that transcriptomic modifications are also accompanied by changes in genome binding of the master osteoblast regulator RUNX2 and the chromatin remodeler CCAAT/enhancer-binding protein β. Importantly, genome occupancy was also highlighted by the presence of epigenetic enhancer signatures that were selectively changed in response to both differentiation and 1,25(OH)2D3. The impact of VDR, RUNX2, and C/EBPβ on osteoblast differentiation is exemplified by their actions at the Runx2 and Sp7 gene loci. We conclude that each of these mechanisms may contribute to the diverse actions of 1,25(OH)2D3 on differentiating osteoblasts.