Genome-wide identification and annotation of HIF-1α binding sites in two cell lines using massively parallel sequencing

Transcriptional Response to Hypoxia: The Role of HIF-1-Associated Co-Regulators

The Hypoxia Inducible Factor 1 (HIF-1) plays a major role in the cellular response to hypoxia by regulating the expression of many genes involved in adaptive processes that allow cell survival under low oxygen conditions. Adaptation to the hypoxic tumor micro-environment is also critical for cancer cell proliferation and therefore HIF-1 is also considered a valid therapeutical target. Despite the huge progress in understanding regulation of HIF-1 expression and activity by oxygen levels or oncogenic pathways, the way HIF-1 interacts with chromatin and the transcriptional machinery in order to activate its target genes is still a matter of intense investigation. Recent studies have identified several different HIF-1- and chromatin-associated co-regulators that play important roles in the general transcriptional activity of HIF-1, independent of its expression levels, as well as in the selection of binding sites, promoters and target genes, which, however, often depends on cellular context. We review here these co-regulators and examine their effect on the expression of a compilation of well-characterized HIF-1 direct target genes in order to assess the range of their involvement in the transcriptional response to hypoxia. Delineating the mode and the significance of the interaction between HIF-1 and its associated co-regulators may offer new attractive and specific targets for anticancer therapy.

Systems approaches to understand oxygen sensing: how multi-omics has driven advances in understanding oxygen-based signalling

Hypoxia is a common denominator in the pathophysiology of a variety of human disease states. Insight into how cells detect, and respond to low oxygen is crucial to understanding the role of hypoxia in disease. Central to the hypoxic response is rapid changes in the expression of genes essential to carry out a wide range of functions to adapt the cell/tissue to decreased oxygen availability. These changes in gene expression are co-ordinated by specialised transcription factors, changes to chromatin architecture and intricate balances between protein synthesis and destruction that together establish changes to the cellular proteome. In this article, we will discuss the advances of our understanding of the cellular oxygen sensing machinery achieved through the application of 'omics-based experimental approaches.

A functional map of genomic HIF1α-DNA complexes in the eye lens revealed through multiomics analysis

BACKGROUND

During eye lens development the embryonic vasculature regresses leaving the lens without a direct oxygen source. Both embryonically and throughout adult life, the lens contains a decreasing oxygen gradient from the surface to the core that parallels the natural differentiation of immature surface epithelial cells into mature core transparent fiber cells. These properties of the lens suggest a potential role for hypoxia and the master regulator of the hypoxic response, hypoxia-inducible transcription factor 1 (HIF1), in the regulation of genes required for lens fiber cell differentiation, structure and transparency. Here, we employed a multiomics approach combining CUT&RUN, RNA-seq and ATACseq analysis to establish the genomic complement of lens HIF1α binding sites, genes activated or repressed by HIF1α and the chromatin states of HIF1α-regulated genes.

RESULTS

CUT&RUN analysis revealed 8375 HIF1α-DNA binding complexes in the chick lens genome. One thousand one hundred ninety HIF1α-DNA binding complexes were significantly clustered within chromatin accessible regions (χ2 test p < 1 × 10- 55) identified by ATACseq. Formation of the identified HIF1α-DNA complexes paralleled the activation or repression of 526 genes, 116 of which contained HIF1α binding sites within 10kB of the transcription start sites. Some of the identified HIF1α genes have previously established lens functions while others have novel functions never before examined in the lens. GO and pathway analysis of these genes implicate HIF1α in the control of a wide-variety of cellular pathways potentially critical for lens fiber cell formation, structure and function including glycolysis, cell cycle regulation, chromatin remodeling, Notch and Wnt signaling, differentiation, development, and transparency.

CONCLUSIONS

These data establish the first functional map of genomic HIF1α-DNA complexes in the eye lens. They identify HIF1α as an important regulator of a wide-variety of genes previously shown to be critical for lens formation and function and they reveal a requirement for HIF1α in the regulation of a wide-variety of genes not yet examined for lens function. They support a requirement for HIF1α in lens fiber cell formation, structure and function and they provide a basis for understanding the potential roles and requirements for HIF1α in the development, structure and function of more complex tissues.

Now a Nobel gas: oxygen

The recent bestowal of the Nobel Prize 2019 in Physiology or Medicine to Gregg L. Semenza, Sir Peter J. Ratcliffe, and William G. Kaelin Jr. celebrates a series of remarkable discoveries that span from the physiological research question on how oxygen deficiency (hypoxia) induces the red blood cell forming hormone erythropoietin (Epo) to the first clinical application of a novel family of Epo-inducing drugs to treat patients suffering from renal anemia. This review looks back at the most important findings made by the three Nobel laureates, highlights current research trends, and sheds an eye on future perspectives of hypoxia research, including emerging and potential clinical applications.

Mechanisms of hypoxia signalling: new implications for nephrology

Studies of the regulation of erythropoietin (EPO) production by the liver and kidneys, one of the classical physiological responses to hypoxia, led to the discovery of human oxygen-sensing mechanisms, which are now being targeted therapeutically. The oxygen-sensitive signal is generated by 2-oxoglutarate-dependent dioxygenases that deploy molecular oxygen as a co-substrate to catalyse the post-translational hydroxylation of specific prolyl and asparaginyl residues in hypoxia-inducible factor (HIF), a key transcription factor that regulates transcriptional responses to hypoxia. Hydroxylation of HIF at different sites promotes both its degradation and inactivation. Under hypoxic conditions, these processes are suppressed, enabling HIF to escape destruction and form active transcriptional complexes at thousands of loci across the human genome. Accordingly, HIF prolyl hydroxylase inhibitors stabilize HIF and stimulate expression of HIF target genes, including the EPO gene. These molecules activate endogenous EPO gene expression in diseased kidneys and are being developed, or are already in clinical use, for the treatment of renal anaemia. In this Review, we summarize information on the molecular circuitry of hypoxia signalling pathways underlying these new treatments and highlight some of the outstanding questions relevant to their clinical use.

The Circadian Clock Protein CRY1 Is a Negative Regulator of HIF-1α

The circadian clock and the hypoxia-signaling pathway are regulated by an integrated interplay of positive and negative feedback limbs that incorporate energy homeostasis and carcinogenesis. We show that the negative circadian regulator CRY1 is also a negative regulator of hypoxia-inducible factor (HIF). Mechanistically, CRY1 interacts with the basic-helix-loop-helix domain of HIF-1α via its tail region. Subsequently, CRY1 reduces HIF-1α half-life and binding of HIFs to target gene promoters. This appeared to be CRY1 specific because genetic disruption of CRY1, but not CRY2, affected the hypoxia response. Furthermore, CRY1 deficiency could induce cellular HIF levels, proliferation, and migration, which could be reversed by CRISPR/Cas9- or short hairpin RNA-mediated HIF knockout. Altogether, our study provides a mechanistic explanation for genetic association studies linking a disruption of the circadian clock with hypoxia-associated processes such as carcinogenesis.

Copper affects the binding of HIF-1α to the critical motifs of its target genes

Copper regulates the target gene selection of HIF-1α under hypoxic conditions by affecting HIF-1α-DNA binding patterns across the genome.

The Transcription Factor Hif-1 Enhances the Radio-Resistance of Mouse MSCs

Mesenchymal stromal cells (MSCs) are multipotent progenitors supporting bone marrow hematopoiesis. MSCs have an efficient DNA damage response (DDR) and are consequently relatively radio-resistant cells. Therefore, MSCs are key to hematopoietic reconstitution following total body irradiation (TBI) and bone marrow transplantation (BMT). The bone marrow niche is hypoxic and via the heterodimeric transcription factor Hypoxia-inducible factor-1 (Hif-1), hypoxia enhances the DDR. Using gene knock-down, we have previously shown that the Hif-1α subunit of Hif-1 is involved in mouse MSC radio-resistance, however its exact mechanism of action remains unknown. In order to dissect the involvement of Hif-1α in the DDR, we used CRISPR/Cas9 technology to generate a stable mutant of the mouse MSC cell line MS5 lacking Hif-1α expression. Herein, we show that it is the whole Hif-1 transcription factor, and not only the Hif-1α subunit, that modulates the DDR of mouse MSCs. This effect is dependent upon the presence of a Hif-1α protein capable of binding to both DNA and its heterodimeric partner Arnt (Hif-1β). Detailed transcriptomic and proteomic analysis of Hif1a KO MS5 cells leads us to conclude that Hif-1α may be acting indirectly on the DNA repair process. These findings have important implications for the modulation of MSC radio-resistance in the context of BMT and cancer.

Advances in Hypoxia-Inducible Factor Biology

Hypoxia-inducible factor (HIF), a central regulator for detecting and adapting to cellular oxygen levels, transcriptionally activates genes modulating oxygen homeostasis and metabolic activation. Beyond this, HIF influences many other processes. Hypoxia, in part through HIF-dependent mechanisms, influences epigenetic factors, including DNA methylation and histone acetylation, which modulate hypoxia-responsive gene expression in cells. Hypoxia profoundly affects expression of many noncoding RNAs classes that have clinicopathological implications in cancer. HIF can regulate noncoding RNAs production, while, conversely, noncoding RNAs can modulate HIF expression. There is recent evidence for crosstalk between circadian rhythms and hypoxia-induced signaling, suggesting involvement of molecular clocks in adaptation to fluxes in nutrient and oxygen sensing. HIF induces increased production of cellular vesicles facilitating intercellular communication at a distance-for example, promoting angiogenesis in hypoxic tumors. Understanding the complex networks underlying cellular and genomic regulation in response to hypoxia via HIF may identify novel and specific therapeutic targets.

Sexual dimorphism of the pulmonary transcriptome in neonatal hyperoxic lung injury: identification of angiogenesis as a key pathway

Bronchopulmonary dysplasia (BPD) is characterized by impaired alveolar secondary septation and vascular growth. Exposure to high concentrations of oxygen (hyperoxia) contributes to the development of BPD. The male sex is considered an independent risk factor for the development of BPD. The reasons underlying sexually dimorphic outcomes in premature neonates are not known. We hypothesized that sex-specific modulation of biological processes in the lung under hyperoxic conditions contributes to sex-based differences. Neonatal male and female mice (C57BL/6) were exposed to hyperoxia [95% [Formula: see text], postnatal day (PND) 1-5: saccular stage of lung development] and euthanized on PND 7 or 21. Pulmonary gene expression was studied using RNA-Seq on the Illumina HiSeq 2500 platform. Analysis of the pulmonary transcriptome revealed differential sex-specific modulation of crucial pathways such as angiogenesis, response to hypoxia, inflammatory response, and p53 pathway. Candidate genes from these pathways were validated at the mRNA level by qPCR. Analysis also revealed sex-specific differences in the modulation of crucial transcription factors. Focusing on the differential modulation of the angiogenesis pathway, we also showed sex-specific differential activation of Hif-1α-regulated genes using ChIP-qPCR and differences in expression of crucial genes (Vegf, VegfR2, and Phd2) modulating angiogenesis. We demonstrate the translational relevance of our findings by showing that our murine sex-specific differences in gene expression correlate with those from a preexisting human BPD data set. In conclusion, we provide novel molecular insights into differential sex-specific modulation of the pulmonary transcriptome in neonatal hyperoxic lung injury and highlight angiogenesis as one of the crucial differentially modulated pathways.

Hypoxia inducible factors regulate the transcription of the sprouty2 gene and expression of the sprouty2 protein

Receptor Tyrosine Kinase (RTK) signaling plays a major role in tumorigenesis and normal development. Sprouty2 (Spry2) attenuates RTK signaling and inhibits processes such as angiogenesis, cell proliferation, migration and survival, which are all upregulated in tumors. Indeed in cancers of the liver, lung, prostate and breast, Spry2 protein levels are markedly decreased correlating with poor patient prognosis and shorter survival. Thus, it is important to understand how expression of Spry2 is regulated. While prior studies have focused on the post-translation regulation of Spry2, very few studies have focused on the transcriptional regulation of SPRY2 gene. Here, we demonstrate that in the human hepatoma cell line, Hep3B, the transcription of SPRY2 is inhibited by the transcription regulating hypoxia inducible factors (HIFs). HIFs are composed of an oxygen regulated alpha subunit (HIF1α or HIF2α) and a beta subunit (HIF1β). Intriguingly, silencing of HIF1α and HIF2α elevates SPRY2 mRNA and protein levels suggesting HIFs reduce the transcription of the SPRY2 promoter. In silico analysis identified ten hypoxia response elements (HREs) in the proximal promoter and first intron of SPRY2. Using chromatin immunoprecipitation (ChIP), we show that HIF1α/2α bind near the putative HREs in the proximal promoter and intron of SPRY2. Our studies demonstrated that not only is the SPRY2 promoter methylated, but silencing HIF1α/2α reduced the methylation. ChIP assays also showed DNA methyltransferase1 (DNMT1) binding to the proximal promoter and first intron of SPRY2 and silencing HIF1α/2α decreased this association. Additionally, silencing of DNMT1 mimicked the HIF1α/2α silencing-mediated increase in SPRY2 mRNA and protein. While simultaneous silencing of HIF1α/2α and DNMT1 increased SPRY2 mRNA a little more, the increase was not additive suggesting a common mechanism by which DNMT1 and HIF1α/2α regulate SPRY2 transcription. Together these data suggest that the transcription of SPRY2 is inhibited by HIFs, in part, via DNMT1- mediated methylation.

Genome-wide analysis of HIF-2α chromatin binding sites under normoxia in human bronchial epithelial cells (BEAS-2B) suggests its diverse functions

Constitutive functional HIF-2α was recently identified in cancer and stem cell lines under normoxia. In this study, BEAS-2B, a bronchial epithelial cell line, was shown to constitutively express active HIF-2α under normoxia and exhibit markers of pluripotency including Oct-4, Nanog, and sphere formation. Oct-4 expression was reduced after knockdown of HIF-2α under normoxia. Global enrichment analysis of HIF-2α demonstrated the diverse functions of HIF-2α under normoxia. Bioinformatics analysis of the enriched loci revealed an enhancer role of HIF-2α binding sites, involvement of HIF-2α interacting proteins, and enriched de novo motifs which suggest the diverse role of HIF-2α in pseudohypoxia. The low ratio of the discovered loci overlapping with those revealed in cancer cell lines 786-O (16.1%) and MCF-7 (15.9%) under hypoxia indicated a prevailing non-canonical mechanism. Hypoxia had positive, marginal or adverse effects on the enrichment of the selected loci in ChIP-PCR assays. Deletion of the N-terminal activation domain (N-TAD) of HIF-2α disrupted the reporting activity of two of the loci annotated to ELN and ANKRD31. Hypoxia incurring abundance variation of HIF-2α may misrepresent the N-TAD functions as canonical hypoxia inducible features via C-TAD activation. Elucidation of the pseudohypoxia functions of constitutive HIF-2α is useful for resolving its role in malignancy and pluripotency.

Hypoxia, Hypoxia-inducible Transcription Factors, and Renal Cancer

CONTEXT

Renal cancer is a common urologic malignancy, and therapeutic options for metastatic disease are limited. Most clear cell renal cell carcinomas (ccRCC) are associated with loss of von Hippel-Lindau tumor suppressor (pVHL) function and deregulation of hypoxia pathways.

OBJECTIVE

This review summarizes recent evidence from genetic and biological studies showing that hypoxia and hypoxia-related pathways play critical roles in the development and progress of renal cancer.

EVIDENCE ACQUISITION

We used a systematic search for articles using the keywords hypoxia, HIF, renal cancer, and VHL.

EVIDENCE SYNTHESIS

Identification of the tumor suppressor pVHL has allowed the characterization of important ccRCC-associated pathways. pVHL targets α-subunits of hypoxia-inducible transcription factors (HIF) for proteasomal degradation. The two main HIF-α isoforms have opposing effects on RCC biology, possibly through distinct interactions with additional oncogenes. Furthermore, HIF-1α activity is commonly diminished by chromosomal deletion in ccRCCs, and increased HIF-1 activity reduces tumor burden in xenograft tumor models. Conversely, polymorphisms at the HIF-2α gene locus predispose to the development of ccRCCs, and HIF-2α promotes tumor growth. Genetic studies have revealed a prominent role for chromatin-modifying enzyme genes in ccRCC, and these may further modulate specific aspects of the HIF response. This suggests that, rather than global activation of HIF, specific components of the response are important in promoting kidney cancer. Some of these processes are already targets for current therapeutic strategies, and further dissection of this pathway might yield novel methods of treating RCC.

CONCLUSIONS

In contrast to many tumor types, HIF-1α and HIF-2α have opposing effects in ccRCC biology, with HIF-1α acting as a tumor suppressor and HIF-2α acting as an oncogene. The overall effect of VHL inactivation will depend on fine-tuning of the HIF response.

PATIENT SUMMARY

High levels of hypoxia-inducible transcription factors (HIF) are particularly important in the clear cell type of kidney cancer, in which they are no longer properly regulated by the von Hippel-Lindau protein. The two HIF-α proteins have opposing effects on tumor evolution.

Integrative Analysis of DCE-MRI and Gene Expression Profiles in Construction of a Gene Classifier for Assessment of Hypoxia-Related Risk of Chemoradiotherapy Failure in Cervical Cancer

PURPOSE

A 31-gene expression signature reflected in dynamic contrast enhanced (DCE)-MR images and correlated with hypoxia-related aggressiveness in cervical cancer was identified in previous work. We here aimed to construct a dichotomous classifier with key signature genes and a predefined classification threshold that separated cervical cancer patients into a more and less hypoxic group with different outcome to chemoradiotherapy.

EXPERIMENTAL DESIGN

A training cohort of 42 patients and two independent cohorts of 108 and 131 patients were included. Gene expression data were generated from tumor biopsies by two Illumina array generations (WG-6, HT-12). Technical transfer of the classifier to a reverse transcription quantitative PCR (RT-qPCR) platform was performed for 74 patients. The amplitude ABrix in the Brix pharmacokinetic model was extracted from DCE-MR images of 64 patients and used as an indicator of hypoxia.

RESULTS

Classifier candidates were constructed by integrative analysis of ABrix and gene expression profiles in the training cohort and evaluated by a leave-one-out cross-validation approach. On the basis of their ability to separate patients correctly according to hypoxia status, a 6-gene classifier was identified. The classifier separated the patients into two groups with different progression-free survival probability. The robustness of the classifier was demonstrated by successful validation of hypoxia association and prognostic value across cohorts, array generations, and assay platforms. The prognostic value was independent of existing clinical markers, regardless of clinical endpoints.

CONCLUSIONS

A robust DCE-MRI-associated gene classifier has been constructed that may be used to achieve an early indication of patients' risk of hypoxia-related chemoradiotherapy failure. Clin Cancer Res; 22(16); 4067-76. ©2016 AACR.

Spatial Interplay between Polycomb and Trithorax Complexes Controls Transcriptional Activity in T Lymphocytes

Trithorax group (TrxG) and Polycomb group (PcG) proteins are two mutually antagonistic chromatin modifying complexes, however, how they together mediate transcriptional counter-regulation remains unknown. Genome-wide analysis revealed that binding of Ezh2 and menin, central members of the PcG and TrxG complexes, respectively, were reciprocally correlated. Moreover, we identified a developmental change in the positioning of Ezh2 and menin in differentiated T lymphocytes compared to embryonic stem cells. Ezh2-binding upstream and menin-binding downstream of the transcription start site was frequently found at genes with higher transcriptional levels, and Ezh2-binding downstream and menin-binding upstream was found at genes with lower expression in T lymphocytes. Interestingly, of the Ezh2 and menin cooccupied genes, those exhibiting occupancy at the same position displayed greatly enhanced sensitivity to loss of Ezh2. Finally, we also found that different combinations of Ezh2 and menin occupancy were associated with expression of specific functional gene groups important for T cell development. Therefore, spatial cooperative gene regulation by the PcG and TrxG complexes may represent a novel mechanism regulating the transcriptional identity of differentiated cells.

Heterogeneous Effects of Direct Hypoxia Pathway Activation in Kidney Cancer

General activation of hypoxia-inducible factor (HIF) pathways is classically associated with adverse prognosis in cancer and has been proposed to contribute to oncogenic drive. In clear cell renal carcinoma (CCRC) HIF pathways are upregulated by inactivation of the von-Hippel-Lindau tumor suppressor. However HIF-1α and HIF-2α have contrasting effects on experimental tumor progression. To better understand this paradox we examined pan-genomic patterns of HIF DNA binding and associated gene expression in response to manipulation of HIF-1α and HIF-2α and related the findings to CCRC prognosis. Our findings reveal distinct pan-genomic organization of canonical and non-canonical HIF isoform-specific DNA binding at thousands of sites. Overall associations were observed between HIF-1α-specific binding, and genes associated with favorable prognosis and between HIF-2α-specific binding and adverse prognosis. However within each isoform-specific set, individual gene associations were heterogeneous in sign and magnitude, suggesting that activation of each HIF-α isoform contributes a highly complex mix of pro- and anti-tumorigenic effects.

Destruction of a distal hypoxia response element abolishes trans-activation of the PAG1 gene mediated by HIF-independent chromatin looping

A crucial step in the cellular adaptation to oxygen deficiency is the binding of hypoxia-inducible factors (HIFs) to hypoxia response elements (HREs) of oxygen-regulated genes. Genome-wide HIF-1α/2α/β DNA-binding studies revealed that the majority of HREs reside distant to the promoter regions, but the function of these distal HREs has only been marginally studied in the genomic context. We used chromatin immunoprecipitation (ChIP), gene editing (TALEN) and chromosome conformation capture (3C) to localize and functionally characterize a 82 kb upstream HRE that solely drives oxygen-regulated expression of the newly identified HIF target gene PAG1. PAG1, a transmembrane adaptor protein involved in Src signalling, was hypoxically induced in various cell lines and mouse tissues. ChIP and reporter gene assays demonstrated that the -82 kb HRE regulates PAG1, but not an equally distant gene further upstream, by direct interaction with HIF. Ablation of the consensus HRE motif abolished the hypoxic induction of PAG1 but not general oxygen signalling. 3C assays revealed that the -82 kb HRE physically associates with the PAG1 promoter region, independent of HIF-DNA interaction. These results demonstrate a constitutive interaction between the -82 kb HRE and the PAG1 promoter, suggesting a physiologically important rapid response to hypoxia.

Cell type-dependent modulation of the gene encoding heat shock protein HSPA2 by hypoxia-inducible factor HIF-1: Down-regulation in keratinocytes and up-regulation in HeLa cells

HSPA2 belongs to the multigene HSPA family, whose members encode chaperone proteins. Although expression and function of HSPA2 is mainly associated with spermatogenesis, recent studies demonstrated that in humans, the gene is active in various cancers, as well as in normal tissues, albeit in a cell type-specific manner. In the epidermis, HSPA2 is expressed in keratinocytes in the basal layer. Currently, the mechanisms underlying the regulation of HSPA2 expression remain unknown. This study was aimed at determining whether HIF-1 and its binding site, the hypoxia-response element (HRE) located in the HSPA2 promoter, are involved in HSPA2 regulation. As a model system, we used an immortal human keratinocyte line (HaCaT) and cervical cancer cells (HeLa) grown under control or hypoxic conditions. Using an in vitro gene reporter assay, we demonstrated that in keratinocytes HSPA2 promoter activity is reduced under conditions that facilitate stabilization of HIF-1α, whereas HIF-1 inhibitors abrogated the suppressive effect of hypoxia on promoter activity. Chromatin immunoprecipitation revealed that HIF-1α binds to the HSPA2 promoter. In keratinocytes, hypoxia or overexpression of a stable form of HIF-1α attenuated the expression of endogenous HSPA2, whereas targeted repression of HIF-1α by RNAi increased transcription of HSPA2 under hypoxia. Conversely, in HeLa cells, HSPA2 expression increased under conditions that stimulated HIF-1α activity, whereas inhibition of HIF-1α abrogated hypoxia-induced up-regulation of HSPA2 expression. Taken together, our results demonstrate that HIF-1 can exert differential, cell context-dependent regulatory control of the HSPA2 gene. Additionally, we also showed that HSPA2 expression can be stimulated during hypoxia/reoxygenation stress.

HIF1α Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes

Transcriptional mediators of cell stress pathways, including HIF1α, ATF4, and p53, are key to normal development and play critical roles in disease, including ischemia and cancer. Despite their importance, mechanisms by which pathways mediated by these transcription factors interact with one another are not fully understood. In addressing the controversial role of HIF1α in cardiomyocytes (CMs) during heart development, we discovered a mid-gestational requirement for HIF1α for proliferation of hypoxic CMs, involving metabolic switching and a complex interplay among HIF1α, ATF4, and p53. Loss of HIF1α resulted in activation of ATF4 and p53, the latter inhibiting CM proliferation. Bioinformatic and biochemical analyses revealed unexpected mechanisms by which HIF1α intersects with ATF4 and p53 pathways. Our results highlight previously undescribed roles of HIF1α and interactions among major cell stress pathways that could be targeted to enhance proliferation of CMs in ischemia and may have relevance to other diseases, including cancer.

Analysis of RNA decay factor mediated RNA stability contributions on RNA abundance

Background

Histone epigenome data determined by chromatin immunoprecipitation sequencing (ChIP-seq) is used in identifying transcript regions and estimating expression levels. However, this estimation does not always correlate with eventual RNA expression levels measured by RNA sequencing (RNA-seq). Part of the inconsistency may arise from the variance in RNA stability, where the transcripts that are more or less abundant than predicted RNA expression from histone epigenome data are inferred to be more or less stable. However, there is little systematic analysis to validate this assumption. Here, we used stability data of whole transcriptome measured by 5′-bromouridine immunoprecipitation chase sequencing (BRIC-seq), which enabled us to determine the half-lives of whole transcripts including lincRNAs, and we integrated BRIC-seq with ChIP-seq to achieve better estimation of the eventual transcript levels and to understand the importance of post-transcriptional regulation that determine the eventual transcript levels.

Results

We identified discrepancies between the RNA abundance estimated by ChIP-seq and measured RNA expression from RNA-seq; for number of genes and estimated that the expression level of 865 genes was controlled at the level of RNA stability in HeLa cells. ENCODE data analysis supported the idea that RNA stability control aids to determine transcript levels in multiple cell types. We identified UPF1, EXOSC5 and STAU1, well-studied RNA degradation factors, as controlling factors for 8% of cases. Computational simulations reasonably explained the changes of eventual mRNA levels attributable to the changes in the rates of mRNA half-lives. In addition, we propose a feedback circuit that includes the regulated degradation of mRNAs encoding transcription factors to maintain the steady state level of RNA abundance. Intriguingly, these regulatory mechanisms were distinct between mRNAs and lincRNAs.

Conclusions

Integrative analysis of ChIP-seq, RNA-seq and our BRIC-seq showed that transcriptional regulation and RNA degradation are independently regulated. In addition, RNA stability is an important determinant of eventual transcript levels. RNA binding proteins, such as UPF1, STAU1 and EXOSC5 may play active roles in such controls.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1358-y) contains supplementary material, which is available to authorized users.

Aberrant transcriptional regulations in cancers: genome, transcriptome and epigenome analysis of lung adenocarcinoma cell lines

Here we conducted an integrative multi-omics analysis to understand how cancers harbor various types of aberrations at the genomic, epigenomic and transcriptional levels. In order to elucidate biological relevance of the aberrations and their mutual relations, we performed whole-genome sequencing, RNA-Seq, bisulfite sequencing and ChIP-Seq of 26 lung adenocarcinoma cell lines. The collected multi-omics data allowed us to associate an average of 536 coding mutations and 13,573 mutations in promoter or enhancer regions with aberrant transcriptional regulations. We detected the 385 splice site mutations and 552 chromosomal rearrangements, representative cases of which were validated to cause aberrant transcripts. Averages of 61, 217, 3687 and 3112 mutations are located in the regulatory regions which showed differential DNA methylation, H3K4me3, H3K4me1 and H3K27ac marks, respectively. We detected distinct patterns of aberrations in transcriptional regulations depending on genes. We found that the irregular histone marks were characteristic to EGFR and CDKN1A, while a large genomic deletion and hyper-DNA methylation were most frequent for CDKN2A. We also used the multi-omics data to classify the cell lines regarding their hallmarks of carcinogenesis. Our datasets should provide a valuable foundation for biological interpretations of interlaced genomic and epigenomic aberrations.

Genome-wide identification of hypoxia-inducible factor-1 and -2 binding sites in hypoxic human macrophages alternatively activated by IL-10

Macrophages (MΦ) often accumulate in hypoxic areas, where they significantly influence disease progression. Anti-inflammatory cytokines, such as IL-10, generate alternatively activated macrophages that support tumor growth. To understand how alternative activation affects the transcriptional profile of hypoxic macrophages, we globally mapped binding sites of hypoxia-inducible factor (HIF)-1α and HIF-2α in primary human monocyte-derived macrophages prestimulated with IL-10. 713 HIF-1 and 795 HIF-2 binding sites were identified under hypoxia. Pretreatment with IL-10 altered the binding pattern, with 120 new HIF-1 and 188 new HIF-2 binding sites emerging. HIF-1 binding was most prominent in promoters, while HIF-2 binding was more abundant in enhancer regions. Comparison of ChIP-seq data obtained in other cells revealed a highly cell type specific binding of HIF. In MΦ HIF binding occurred preferentially in already active enhancers or promoters. To assess the roles of HIF on gene expression, primary human macrophages were treated with siRNA against HIF-1α or HIF-2α, followed by genome-wide gene expression analysis. Comparing mRNA expression to the HIF binding profile revealed a significant enrichment of hypoxia-inducible genes previously identified by ChIP-seq. Analysis of gene expression under hypoxia alone and hypoxia/IL-10 showed the enhanced induction of a set of genes including PLOD2 and SLC2A3, while another group including KDM3A and ADM remained unaffected or was reduced by IL-10. Taken together IL-10 influences the DNA binding pattern of HIF and the level of gene induction.

Construction of mate pair full-length cDNAs libraries and characterization of transcriptional start sites and termination sites

To identify and characterize transcript structures ranging from transcriptional start sites (TSSs) to poly(A)-addition sites (PASs), we constructed and analyzed human TSS/PAS mate pair full-length cDNA libraries from 14 tissue types and four cell lines. The collected information enabled us to define TSS cluster (TSC) and PAS cluster (PAC) relationships for a total of 8530/9400 RefSeq genes, as well as 4251/5618 of their putative alternative promoters/terminators and 4619/4605 intervening transcripts, respectively. Analyses of the putative alternative TSCs and alternative PACs revealed that their selection appeared to be mostly independent, with rare exceptions. In those exceptional cases, pairs of transcript units rarely overlapped one another and were occasionally separated by Rad21/CTCF. We also identified a total of 172 similar cases in which TSCs and PACs spanned adjacent but distinct genes. In these cases, different transcripts may utilize different functional units of a particular gene or of adjacent genes. This approach was also useful for identifying fusion gene transcripts in cancerous cells. Furthermore, we could construct cDNA libraries in which 3′-end mate pairs were distributed randomly over the transcripts. These libraries were useful for assembling the internal structure of previously uncharacterized alternative promoter products, as well as intervening transcripts.

Molecular responses to hypoxia-inducible factor 1α and beyond

Cellular response to changes in oxygen tension during normal development or pathologic processes is, in part, regulated by hypoxia-inducible factor (HIF), an oxygen-sensitive transcription factor. HIF activity is primarily controlled through post-translational modifications and stabilization of HIF-1α and HIF-2α proteins and is regulated by a number of cellular pathways involving both oxygen-dependent and -independent mechanisms. Stabilization of HIF-1α activates transcription of genes that participate in key pathways in carcinogenesis, such as angiogenesis, dedifferentiation, and invasion. Since its discovery more than two decades ago, HIF-1α has become a hot topic in molecular research and has been implicated not only in disease pathology but also in prognosis. In this review, we will focus on recent insights into HIF-1α regulation, function, and gene expression. We will also discuss emerging data on the involvement of HIF in cancer prognosis and therapeutic interventions.

A comparison of the rest complex binding patterns in embryonic stem cells and epiblast stem cells

We detected and characterized the binding sites of the representative Rest complex components Rest, Sin3A, and Lsd1. We compared their binding patterns in mouse embryonic stem (ES) cells and epiblast stem (EpiS) cells. We found few Rest sites unique to the EpiS cells. The ES-unique site features were distinct from those of the common sites, namely, the signal intensities were weaker, and the characteristic gene function categories differed. Our analyses showed that the Rest binding sites do not always overlap with the Sin3A and Lsd1 binding sites. The Sin3A binding pattern differed remarkably between the ES and EpiS cells and was accompanied by significant changes in acetylated-histone patterns in the surrounding regions. A series of transcriptome analyses in the same cell types unexpectedly showed that the putative target gene transcript levels were not dramatically different despite dynamic changes in the Rest complex binding patterns and chromatin statuses, which suggests that Rest is not the sole determinant of repression at its targets. Nevertheless, we identified putative Rest targets with explicitly enhanced transcription upon Rest knock-down in 143 and 60 common and ES-unique Rest target genes, respectively. Among such sites, several genes are involved in ES cell proliferation. In addition, we also found that long, intergenic non-coding RNAs were apparent Rest targets and shared similar features with the protein-coding target genes. Interestingly, such non-coding target genes showed less conservation through evolution than protein-coding targets. As a result of differences in the components and targets of the Rest complex, its functional roles may differ in ES and EpiS cells.

Hypoxia inducible factor pathway inhibitors as anticancer therapeutics

Hypoxia is a significant feature of solid tumor cancers. Hypoxia leads to a more malignant phenotype that is resistant to chemotherapy and radiation, is more invasive and has greater metastatic potential. Hypoxia activates the hypoxia inducible factor (HIF) pathway, which mediates the biological effects of hypoxia in tissues. The HIF complex acts as a transcription factor for many genes that increase tumor survival and proliferation. To date, many HIF pathway inhibitors indirectly affect HIF but there have been no clinically approved direct HIF inhibitors. This can be attributed to the complexity of the HIF pathway, as well as to the challenges of inhibiting protein-protein interactions.

Genome-wide analysis of FOXO3 mediated transcription regulation through RNA polymerase II profiling

Forkhead box O (FOXO) transcription factors are key players in diverse cellular processes affecting tumorigenesis, stem cell maintenance and lifespan. To gain insight into the mechanisms of FOXO-regulated target gene expression, we studied genome-wide effects of FOXO3 activation. Profiling RNA polymerase II changes shows that FOXO3 regulates gene expression through transcription initiation. Correlative analysis of FOXO3 and RNA polymerase II ChIP-seq profiles demonstrates FOXO3 to act as a transcriptional activator. Furthermore, this analysis reveals a significant part of FOXO3 gene regulation proceeds through enhancer regions. FOXO3 binds to pre-existing enhancers and further activates these enhancers as shown by changes in histone acetylation and RNA polymerase II recruitment. In addition, FOXO3-mediated enhancer activation correlates with regulation of adjacent genes and pre-existence of chromatin loops between FOXO3 bound enhancers and target genes. Combined, our data elucidate how FOXOs regulate gene transcription and provide insight into mechanisms by which FOXOs can induce different gene expression programs depending on chromatin architecture.