Functional and HRE motifs count analysis of induction of selected hypoxia-responsive genes by HIF-1 and HIF-2 in human umbilical endothelial cells.

We analyzed the effects of selective knockdown of either HIF-1α or HIF-2α on the transcriptional response to hypoxia of human umbilical endothelial cells at two time-points (2 h and 8 h) of hypoxia. We focused on 13 previously identified hypoxia-responsive genes, pre-selected to have different activation kinetics and different proportions of HRE motifs annotated to either HIF-1 or HIF-2 in open promoters - open chromatin DNase-hypersensitive sites (DHS) regions within ±1 kb of the gene start. We report that genes activated by both HIF-1 and 2 tend to be activated earlier than genes activated by HIF-1 only, which, in turn, tend to be activated earlier than genes activated by HIF-2 only. Moreover, for the 13 analyzed genes, we found that the effect of silencing HIF1A on the gene induction by hypoxia is greater for the genes with more HRE motifs annotated to HIF-1 in their promoter open chromatin DHS regions within ±1 kb and also within ±10 kb of the gene start. We corroborated and extended this finding by showing that among 232 genes previously identified as activated by hypoxia, the genes with ChIP-seq peak(s) for HIF-1α within a ±10 kb flank of the gene start contain more HRE motifs annotated to HIF-1 in the DHS regions within this flank than the genes with no ChIP-seq peaks. Also in the whole genome, the DHS regions intersecting ChIP-seq peaks for HIF-1α contain more HRE motifs annotated to HIF-1 than the DHS regions not intersecting the ChIP-seq peaks. This suggests a mechanism, by which higher promoter content of HRE motifs in DHS regions increases HIF-1 binding, which in turn increases gene induction by hypoxia.