An integrative genomics approach identifies Hypoxia Inducible Factor-1 (HIF-1)-target genes that form the core response to hypoxia

Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction

The transcriptional response driven by Hypoxia-inducible factor (HIF) is central to the adaptation to oxygen restriction. Hence, the complete identification of HIF targets is essential for understanding the cellular responses to hypoxia. Herein we describe a computational strategy based on the combination of phylogenetic footprinting and transcription profiling meta-analysis for the identification of HIF-target genes. Comparison of the resulting candidates with published HIF1a genome-wide chromatin immunoprecipitation indicates a high sensitivity (78%) and specificity (97.8%). To validate our strategy, we performed HIF1a chromatin immunoprecipitation on a set of putative targets. Our results confirm the robustness of the computational strategy in predicting HIF-binding sites and reveal several novel HIF targets, including RE1-silencing transcription factor co-repressor (RCOR2). In addition, mapping of described polymorphisms to the predicted HIF-binding sites identified several single-nucleotide polymorphisms (SNPs) that could alter HIF binding. As a proof of principle, we demonstrate that SNP rs17004038, mapping to a functional hypoxia response element in the macrophage migration inhibitory factor (MIF) locus, prevents induction of this gene by hypoxia. Altogether, our results show that the proposed strategy is a powerful tool for the identification of HIF direct targets that expands our knowledge of the cellular adaptation to hypoxia and provides cues on the inter-individual variation in this response.

Functional characterization of two CITED3 homologs (gcCITED3a and gcCITED3b) in the hypoxia-tolerant grass carp, Ctenopharyngodon idellus

Background

CITED proteins belong to a family of non-DNA-binding transcriptional co-regulators that are characterized by a conserved ED-rich domain at the C-terminus. This family of genes is involved in the regulation of a variety of transcriptional responses through interactions with the CBP/p300 integrators and various transcription factors. In fish, very little is known about the expression and functions of CITEDs.

Results

We have characterized two closely related but distinct CITED3 genes, gcCited3a and gcCited3b, from the hypoxia-tolerant grass carp. The deduced gcCITED3a and gcCITED3b proteins share 72% amino acid identity, and are highly similar to the CITED3 proteins of both chicken and Xenopus. Northern blot analysis indicates that the mRNA expression of gcCited3a and gcCited3b is strongly induced by hypoxia in the kidney and liver, respectively. Luciferase reporter assays demonstrated that both gene promoters are activated by gcHIF-1. Further, ChIP assays comparing normal and hypoxic conditions reveal differential in vivo binding of gcHIF-1 to both gene promoters in kidney and liver tissues. HRE-luciferase reporter assays demonstrated that both gcCITED3a and gcCITED3b proteins inhibit gcHIF-1 transcriptional activity, and GST pull-down assays confirmed that both proteins bind specifically to the CH1 domain of the grass carp p300 protein.

Conclusion

The grass carp gcCITED3a and gcCITED3b genes are differentially expressed and regulated in different fish organs in response to hypoxic stress. This is the first report demonstrating in vivo regulation of two closely-related CITED3 isogenes by HIF-1, as well as CITED3 regulation of HIF-1 transcriptional activity in fish. Overall, our findings suggest that unique molecular mechanisms operate through these two gcCITED3 isoforms that likely play an important regulatory role in the hypoxic response in the grass carp.

Generating specificity and diversity in the transcriptional response to hypoxia

The sensing of oxygen levels and maintenance of oxygen homeostasis is crucial for cells. The hypoxic-sensitive regulation of gene expression allows information about the oxygen status to be converted into appropriate cellular responses. Although there is a core transcriptional pathway, the signalling cascade can be modified to allow diversity and specificity in the transcriptional output. In this Review, we discuss recent advances in our understanding of the mechanisms and factors that contribute to the observed diversity and specificity. A deeper knowledge about how hypoxic signalling is tuned will further our understanding of the cellular hypoxic response in normal physiology and how it becomes derailed in disease.

Validation of suitable house keeping genes for hypoxia-cultured human chondrocytes

Background

Hypoxic culturing of chondrocytes is gaining increasing interest in cartilage research. Culturing of chondrocytes under low oxygen tension has shown several advantages, among them increased synthesis of extracellular matrix and increased redifferentiation of dedifferentiated chondrocytes. Quantitative gene expression analyses such as quantitative real-time PCR (qRT-PCR) are powerful tools in the investigation of underlying mechanisms of cell behavior and are used routinely for differentiation and phenotype assays. However, the genes used for normalization in normoxic cell-cultures might not be suitable in the hypoxic environment. The objective of this study was to determine hypoxia-stable housekeeping genes (HKG) for quantitative real-time PCR (qRT-PCR) in human chondrocytes cultured in 21%, 5% and 1% oxygen by geNorm and NormFinder analyses.

Results

The chondrocytic response to the hypoxic challange was validated by a significant increase in expression of the hypoxia-inducible gene ankyrin repeat 37 as well as SOX9 in hypoxia. When cultured on the 3-dimentional (3D) scaffold TATA-binding protein (TBP) exhibited the highest expression stability with NormFinder while Ribosomal protein L13a (RPL13A) and beta2-microglobulin (B2M) were the most stable using geNorm analysis. In monolayer RPL13A were the most stable gene using NormFinder, while geNorm assessed RPL13A and human RNA polymerase II (RPII) as most stable. When examining the combination of (3D) culturing and monolayer RPL13A and B2M showed the highest expression stability from geNorm analysis while RPL13A also showed the highest expression stability using NormFinder. Often used HKG such as beta actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were the most unstable genes investigated in all comparisons. The pairwise variations for the two most stable HKG in each group were all below the cut-off value of 0.15, suggesting that the two most stable HKG from geNorm analysis would be sufficient for qRT-PCR.

Conclusion

All data combined we recommend RPL13A, B2M and RPII as the best choice for qRT-PCR analyses when comparing normoxic and hypoxic cultured human chondrocytes although other genes might also be suitable. However, the matching of HKG to target genes by means of a thorough investigation of the stability in each study would always be preferable.