The genome-wide analysis of the binding sites of the transcription factor vitamin D receptor (VDR) is essential for a global appreciation the physiological impact of the nuclear hormone 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃). Genome-wide analysis of lipopolysaccharide (LPS)-polarized THP-1 human monocytic leukemia cells via chromatin immunoprecipitation sequencing (ChIP-seq) resulted in 1,318 high-confidence VDR binding sites, of which 789 and 364 occurred uniquely with and without 1,25(OH)₂D₃ stimulation, while only 165 were common. We re-analyzed five public VDR ChIP-seq datasets with identical peak calling settings (MACS, version 2) and found, using a novel consensus summit identification strategy, in total 23,409 non-overlapping VDR binding sites, 75% of which are unique within the six analyzed cellular models. LPS-differentiated THP-1 cells have 22% more genomic VDR locations than undifferentiated cells and both cell types display more overlap in their VDR locations than the other investigated cell types. In general, the intersection of VDR binding profiles of ligand-stimulated cells is higher than those of unstimulated cells. De novo binding site searches and HOMER screening for binding motifs formed by direct repeats spaced by three nucleotides (DR3) suggest for all six VDR ChIP-seq datasets that these sequences are found preferentially at highly ligand responsive VDR loci. Importantly, all VDR ChIP-seq datasets display the same relationship between the VDR occupancy and the percentage of DR3-type sequences below the peak summits. The comparative analysis of six VDR ChIP-seq datasets demonstrated that the mechanistic basis for the action of the VDR is independent of the cell type. Only the minority of genome-wide VDR binding sites contains a DR3-type sequence. Moreover, the total number of identified VDR binding sites in each ligand-stimulated cell line inversely correlates with the percentage of peak summits with DR3 sites.

The physiological functions of vitamin D are mediated by its metabolite 1α,25-dihydroxyvitamin D₃ (1,25(OH)₂D₃) activating the transcription factor vitamin D receptor (VDR). In THP-1 human monocytes we demonstrated epigenome-wide effects of 1,25(OH)₂D₃ at 8979 loci with significantly modulated chromatin accessibility. Maximal chromatin opening was observed after 24 h, while after 48 h most sites closed again. The chromatin-organizing protein CTCF bound to 14% of the 1,25(OH)₂D₃-sensitive chromatin regions. Interestingly, 1,25(OH)₂D₃ affected the chromatin association of CTCF providing an additional mechanism for the epigenome-wide effects of the VDR ligand. The 1,25(OH)₂D₃-modulated transcriptome of THP-1 cells comprised 1284 genes, 77.5% of which responded only 24 h after stimulation. During the 1,25(OH)₂D₃ stimulation time course the proportion of down-regulated genes increased from 0% to 44.9% and the top-ranking physiological function of the respective genes shifted from anti-microbial response to connective tissue disorders. The integration of epigenomic and transcriptomic data identified 165 physiologically important 1,25(OH)₂D₃ target genes, including HTT and NOD2, whose expression can be predicted primarily from epigenomic data of their genomic loci. Taken together, a large number of 1,25(OH)₂D₃-triggered epigenome-wide events precede and accompany the transcriptional activation of target genes of the nuclear hormone.

Vitamin D₃ has transcriptome- and genome-wide effects and activates, via the binding of its metabolite 1α,25-dihydroxyvitamin D₃ to the transcription factor vitamin D receptor (VDR), several hundred target genes. Using samples from a 5-month vitamin D₃ intervention study (VitDmet), we recently reported that the expression of 12 VDR target genes in peripheral blood mononuclear cells (PBMCs) as well as 12 biochemical and clinical parameters of the study participants are significantly triggered by vitamin D₃. In this study, we performed a more focused selection of further 12 VDR target genes and demonstrated that changes of their mRNA expression in PBMCs of VitDmet subjects significantly correlate with alterations of 25-hydroxyvitamin D3 serum levels. Network and self-organizing map analysis of these datasets together with that of the other 24 parameters was followed by relevance calculations and identified changes in parathyroid hormone serum levels and the expression of the newly selected genes STS, BCL6, ITGAM, LRRC25, LPGAT1 and TREM1 as well as of the previously reported genes DUSP10 and CD14 as the most relevant parameters for describing vitamin D responsiveness in vivo. Moreover, parameter relevance ranking allowed the segregation of study subjects into high and low responders. Due to the long intervention period the vitamin D response was not too prominent on the level of transcriptional activation. Therefore, we performed in the separate VitDbol trial a short-term but high dose stimulation with a vitamin D3 bolus. In PBMCs of VitDbol subjects we observed direct transcriptional effects on the selected VDR target genes, such as an up to 2.1-fold increase already one day after supplementation onset. In conclusion, both long-term and short-term vitamin D₃ supplementation studies allow monitoring the vitamin D responsiveness of human individuals and represent new types of human in vivo vitamin D₃ investigations.

The transcription factor vitamin D receptor (VDR) is the exclusive nuclear target of the biologically active form of vitamin D (1,25(OH)₂D₃). In THP-1 human monocytes we obtained a highly accurate VDR cistrome after 2 and 24h ligand stimulation comprising >11,600 genomic loci, 78% of which were detected exclusively after 24h. In contrast, a group of 510 persistent VDR sites occurred at all conditions and some 2100 VDR loci were only transiently occupied. Machine learning and statistical analysis as well as a comparison with the re-analyzed B cell VDR cistrome indicated a subgroup of 339 highly conserved persistent VDR sites that were suited best for describing vitamin D-triggered gene regulatory scenarios. The 1,25(OH)₂D₃-dependent transcriptome of THP-1 cells comprised 587 genes, 311 of which were primary targets with main functions in the immune system. More than 97% of the latter genes were located within 1,25(OH)₂D₃-modulated topologically associated domains (TADs). The number of persistent and transient VDR sites was found to be the main discriminator for sorting these TADs into five classes carrying vitamin D target genes involved in distinct biological processes. In conclusion, specific regulation of biological processes by vitamin D depends on differences in time-dependent VDR binding.