The impact of transcriptional cycling on gene regulation

A hierarchical regulatory network analysis of the vitamin D induced transcriptome reveals novel regulators and complete VDR dependency in monocytes

The transcription factor vitamin D receptor (VDR) is the high affinity nuclear target of the biologically active form of vitamin D₃ (1,25(OH)₂D₃). In order to identify pure genomic transcriptional effects of 1,25(OH)₂D₃, we used VDR cistrome, transcriptome and open chromatin data, obtained from the human monocytic cell line THP-1, for a novel hierarchical analysis applying three bioinformatics approaches. We predicted 75.6% of all early 1,25(OH)₂D₃-responding (2.5 or 4 h) and 57.4% of the late differentially expressed genes (24 h) to be primary VDR target genes. VDR knockout led to a complete loss of 1,25(OH)₂D₃–induced genome-wide gene regulation. Thus, there was no indication of any VDR-independent non-genomic actions of 1,25(OH)₂D₃ modulating its transcriptional response. Among the predicted primary VDR target genes, 47 were coding for transcription factors and thus may mediate secondary 1,25(OH)₂D₃ responses. CEBPA and ETS1 ChIP-seq data and RNA-seq following CEBPA knockdown were used to validate the predicted regulation of secondary vitamin D target genes by both transcription factors. In conclusion, a directional network containing 47 partly novel primary VDR target transcription factors describes secondary responses in a highly complex vitamin D signaling cascade. The central transcription factor VDR is indispensable for all transcriptome-wide effects of the nuclear hormone.

An oxidative DNA "damage" and repair mechanism localized in the VEGF promoter is important for hypoxia-induced VEGF mRNA expression

In hypoxia, mitochondria-generated reactive oxygen species not only stimulate accumulation of the transcriptional regulator of hypoxic gene expression, hypoxia inducible factor-1 (Hif-1), but also cause oxidative base modifications in hypoxic response elements (HREs) of hypoxia-inducible genes. When the hypoxia-induced base modifications are suppressed, Hif-1 fails to associate with the HRE of the VEGF promoter, and VEGF mRNA accumulation is blunted. The mechanism linking base modifications to transcription is unknown. Here we determined whether recruitment of base excision DNA repair (BER) enzymes in response to hypoxia-induced promoter modifications was required for transcription complex assembly and VEGF mRNA expression. Using chromatin immunoprecipitation analyses in pulmonary artery endothelial cells, we found that hypoxia-mediated formation of the base oxidation product 8-oxoguanine (8-oxoG) in VEGF HREs was temporally associated with binding of Hif-1α and the BER enzymes 8-oxoguanine glycosylase 1 (Ogg1) and redox effector factor-1 (Ref-1)/apurinic/apyrimidinic endonuclease 1 (Ape1) and introduction of DNA strand breaks. Hif-1α colocalized with HRE sequences harboring Ref-1/Ape1, but not Ogg1. Inhibition of BER by small interfering RNA-mediated reduction in Ogg1 augmented hypoxia-induced 8-oxoG accumulation and attenuated Hif-1α and Ref-1/Ape1 binding to VEGF HRE sequences and blunted VEGF mRNA expression. Chromatin immunoprecipitation-sequence analysis of 8-oxoG distribution in hypoxic pulmonary artery endothelial cells showed that most of the oxidized base was localized to promoters with virtually no overlap between normoxic and hypoxic data sets. Transcription of genes whose promoters lost 8-oxoG during hypoxia was reduced, while those gaining 8-oxoG was elevated. Collectively, these findings suggest that the BER pathway links hypoxia-induced introduction of oxidative DNA modifications in promoters of hypoxia-inducible genes to transcriptional activation.

Chromatin acetylation at transcription start sites and vitamin D receptor binding regions relates to effects of 1α,25-dihydroxyvitamin D3 and histone deacetylase inhibitors on gene expression

The nuclear hormone 1α,25-dihydroxyvitamin D₃ (1α,25(OH)₂D₃ or 1,25D) regulates its target genes via activation of the transcription factor vitamin D receptor (VDR) far more specifically than the chromatin modifier trichostatin A (TsA) via its inhibitory action on histone deacetylases. We selected the thrombomodulin gene locus with its complex pattern of five VDR binding sites and multiple histone acetylation and open chromatin regions as an example to investigate together with a number of reference genes, the primary transcriptional responses to 1α,25(OH)₂D₃ and TsA. Transcriptome-wide, 18.4% of all expressed genes are either up-or down-regulated already after a 90 min TsA treatment; their response pattern to 1α,25(OH)₂D₃ and TsA sorts them into at least six classes. TsA stimulates a far higher number of genes than 1α,25(OH)₂D₃ and dominates the outcome of combined treatments. However, 200 TsA target genes can be modulated by 1α,25(OH)₂D₃ and more than 1000 genes respond only when treated with both compounds. The genomic view on the genes suggests that the degree of acetylation at transcription start sites and VDR binding regions may determine the effect of TsA on mRNA expression and its interference with 1α,25(OH)₂D₃. Our findings hold true also for other HDAC inhibitors and may have implications on dual therapies using chromatin modifiers and nuclear receptor ligands.

Cyclical regulation of the insulin-like growth factor binding protein 3 gene in response to 1alpha,25-dihydroxyvitamin D3

The nuclear receptor vitamin D receptor (VDR) is known to associate with two vitamin D response element (VDRE) containing chromatin regions of the insulin-like growth factor binding protein 3 (IGFBP3) gene. In non-malignant MCF-10A human mammary cells, we show that the natural VDR ligand 1α,25-dihydroxyvitamin D(3) (1α,25(OH)(2)D(3)) causes cyclical IGFBP3 mRNA accumulation with a periodicity of 60 min, while in the presence of the potent VDR agonist Gemini the mRNA is continuously accumulated. Accordingly, VDR also showed cyclical ligand-dependent association with the chromatin regions of both VDREs. Histone deacetylases (HDACs) play an important role both in VDR signalling and in transcriptional cycling. From the 11 HDAC gene family members, only HDAC4 and HDAC6 are up-regulated in a cyclical fashion in response to 1α,25(OH)(2)D(3), while even these two genes do not respond to Gemini. Interestingly, HDAC4 and HDAC6 proteins show cyclical VDR ligand-induced association with both VDRE regions of the IGFBP3 gene, which coincides with histone H4 deacetylation on these regions. Moreover, combined silencing of HDAC4 and HDAC6 abolishes the cycling of the IGFBP3 gene. We assume that due to more efficient VDR interaction, Gemini induces longer lasting chromatin activation and therefore no transcriptional cycling but monotonically increasing IGFBP3 mRNA. In conclusion, 1α,25(OH)(2)D(3) regulates IGFBP3 transcription through short-term cyclical association of VDR, HDAC4 and HDAC6 to both VDRE-containing chromatin regions.

Dynamics of nuclear receptor target gene regulation

Ligand-regulated nuclear receptors, such as estrogen receptors, glucocorticoid receptor, vitamin D receptor, and peroxisome proliferator-activated receptors, belong to the most widely studied and best understood transcription factors. Therefore, the dynamic nature of transcriptional regulation was observed first with different members of the nuclear receptor superfamily, but is now also extended to other transcription factors, such as nuclear factor κB. Dynamic and in part cyclical processes were observed on the level of translocation into the nucleus, association with genomic binding sites, exchange of co-regulators and chromatin modifiers, occurrence of chromatin marks, and activities of RNA polymerase II resulting in mRNA synthesis. In this review, we summarize recent findings on the dynamic regulation of nuclear receptor target genes in the chromatin context.