Scientific Report
Extensive regulation of the non-coding transcriptome by hypoxia: role of HIF in releasing paused RNApol2 €del3,4,†, Spyros Oikonomopoulos1, Carme Camps1, Steffen Grampp4, Hani Choudhry1,2,†, Johannes Scho Adrian L Harris5, Peter J Ratcliffe3,*, Jiannis Ragoussis1,6,7,‡,** & David R Mole3,‡,***
Abstract Hypoxia is central to both ischaemic and neoplastic diseases. However, the non-coding transcriptional response to hypoxia is largely uncharacterized. We undertook integrated genomic analyses of both non-coding and coding transcripts using massively parallel sequencing and interfaced this data with pan-genomic analyses of hypoxia-inducible factor (HIF) and RNApol2 binding in hypoxic cells. These analyses revealed that all classes of RNA are profoundly regulated by hypoxia and implicated HIF as a major direct regulator of both the non-coding and coding transcriptome, acting predominantly through release of pre-bound promoter-paused RNApol2. These findings indicate that the transcriptional response to hypoxia is substantially more extensive than previously considered. Keywords HIF; hypoxia; non-coding; RNApol2; transcription Subject Categories Metabolism; Computational Biology DOI 10.1002/embr.201337642 | Received 17 June 2013 | Revised 30 September 2013 | Accepted 22 October 2013
Introduction Cells respond to physiological insults by altering their transcriptional output. Hypoxia (low oxygen level) is an important environmental stress with a central role in many physiological responses including adaptation to altitude, exercise, growth and development, as well as major pathophysiological processes such as ischaemic
1 2 3 4 5 6 7
70
vascular disease, inflammation, wound healing and cancer [1]. The transcription factor, hypoxia-inducible factor (HIF) orchestrates many of these responses to limit the disturbance of oxygen homeostasis or promote repair processes, through the transactivation of protein-coding genes with key roles in pathways such as apoptosis, differentiation, proliferation, energy metabolism and growth factor production [2]. Given the fundamental role of disturbed oxygen homeostasis in human disease and the potential for therapeutic manipulation of the HIF pathway, much interest has focused on understanding the extent and architecture of these pathways. To date, pan-genomic analyses of the transcriptional response to hypoxia using microarrays have focussed on protein-coding RNAs (mRNAs) and some microRNAs (miRNAs) [3, 4]. However, recent genomic analyses have transformed our perspective of regulatory transcriptional networks. It is now recognised that mRNA forms only part of a much broader transcriptional output that includes even greater numbers of RNAs that are not translated to protein [5]. These ‘non-coding’ RNAs have important regulatory roles that further shape the transcriptional output of the cell. They include short non-coding RNA (200 nucleotides). Recent work has demonstrated altered noncoding RNA expression in many types of cancer and revealed functions in cell cycle regulation, apoptosis, carcinogenesis and metastasis [5, 6]. However, little is known about the mechanisms regulating the expression of the non-coding transcriptome. Here, we provide the first pan-genomic analysis of both non-coding and coding transcriptional responses to hypoxia, including lncRNAs, miRNAs, piwiRNAs, snRNAs, and tRNAs. The analysis
The Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia The Henry Wellcome Building for Molecular Physiology, University of Oxford, Oxford, UK Department of Nephrology and Hypertension, Friedrich-Alexander-University Erlangen- Nuremberg, Erlangen, Germany The Weatherall Institute of Molecular Medicine, University of Oxford, Oxford, UK McGill University and Genome Quebec Innovation Centre, Montreal, Canada BSRC Alexander Fleming, Athens, Greece * Corresponding author: Tel: +44 1865 287990; Fax: +44 1865 287787; E-mail:
[email protected] ** Corresponding author: Tel: +1 514 398 6508; Fax: +1 514 398 1785; E-mail:
[email protected] *** Corresponding author: Tel: +44 1865 287788; Fax: +44 1865 287787; E-mail:
[email protected] † These authors contributed equally to this work. ‡ These authors contributed equally to this work.
EMBO reports
Vol 15 | No 1 | 2014
ª 2014 The Authors. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Hani Choudhry et al
reveals marked bidirectional hypoxia-inducible changes in transcript abundance across all classes of RNA and demonstrates the importance of HIF in regulating non-coding as well as coding transcriptional responses to hypoxia. Correlation with changes in histone marks and RNA polymerase 2 (RNApol2) positioning, indicates that activation of gene expression by HIF occurs commonly through long-range interactions that cause the release of pre-bound, promoter-paused RNApol2.
Results and Discussion Extensive regulation of the coding and non-coding transcriptome by hypoxia To gain a global overview of the transcriptional response to hypoxia, we integrated expression profile and chromatin maps of human breast cancer MCF-7 cells incubated in either 21% ambient oxygen (normoxia) or 1% ambient oxygen (hypoxia) (Fig 1A). Independent directional and nondirectional cDNA libraries were generated from polyadenylated RNA (polyA+) and total RNA depleted of ribosomal RNAs (ribo-) and analysed by high-throughput sequencing. For transcripts present in both data sets, expression levels determined by the two methods showed a high degree of correlation (supplementary Fig S1), although a subset of transcripts showed incomplete adenylation [7]. To determine transcript abundance across all classes of RNA, the ribosome depleted, directional RNA-seq dataset was mapped sequentially to databases of different RNA classes (Fig 1B). 42989 non-ribosomal transcripts were detected in normoxia and 43169 in hypoxia. The raw data for each class is presented in Fig 1C, whilst Fig 1D shows the median fold-change in hypoxia, together with the interquartile and full ranges, after exclusion of low-abundance transcripts where this ratio might be unreliable. This demonstrates that all classes of RNA are regulated by hypoxia, but reveals major classspecific differences. Several classes of transcript (snRNAs, tRNAs and piwiRNAs) demonstrated overall downregulation, whereas mRNAs, lncRNAs and miRNAs demonstrated overall upregulation. To confirm these findings, representative samples of snRNA, mRNA and lncRNA transcripts exhibiting regulation close to the median values were tested by quantitative (q)PCR. This analysis recapitulated differences between the classes, consistent with the finding of overall class-specific regulation (supplementary Fig S2). In addition to these global changes in hypoxia, a number of transcripts in each category show strong up- or downregulation, ranging from 15-fold up to 8-fold down (piwiRNA), >100-fold up to 30-fold down (miRNA) and 4-fold up to 12-fold down (lncRNA) (supplementary Table S1). Several microarray-based pan-genomic analyses have described heterogeneous responses in the miRNA response to hypoxia depending on cell type, hypoxic stimulus and profiling platform – reviewed in McCormick et al [4]. A number of these miRNAs (including mir-210) were also increased in our study (supplementary Table S2). In addition, upregulated miRNAs from our analysis were also induced across a panel of breast cancer cell lines (supplementary Fig S3). As few lncRNAs have been included in previous pan-genomic analyses our data identifies lncRNAs as a new class of transcript that are commonly regulated by hypoxia. The majority of these transcripts remain functionally uncharacterized. However, two lncRNAs
ª 2014 The Authors
EMBO reports
Hypoxic regulation of the non-coding transcriptome
A
B
C
D
Figure 1. Extensive regulation of transcript abundance by hypoxia. A RNA-seq was performed following 24 h culture in 21% (normoxia) or 1% (hypoxia) ambient oxygen. B Outline of the pipeline for mapping ribosomal depleted directional RNA-seq reads. C Raw data plotted as log2 fold-change by hypoxia on the vertical axis versus expression level on the horizontal axis. Vertical lines denote thresholds (higher for classes of longer RNAs) to select transcripts for further analysis. D Box-and-whisker plots of log2 fold-change by hypoxia for filtered transcripts in each class of RNA. The vertical dotted line denotes no foldchange.
(NEAT1 and MALAT1) included in a previous microarray-based analysis [8] were the most strongly upregulated lncRNAs in our dataset and were also upregulated in a panel of breast cancer cell lines (supplementary Fig S4). NEAT1 is a component of interchromatin-paraspeckles, which regulate nuclear retention of Adenosine to Inosine (A-to-I) edited RNAs [9]. MALAT1 promotes cellular proliferation, is associated with tumour growth, metastasis and recurrence and is frequently upregulated or mutated in solid tumours [10]. In addition, H19, another hypoxia-inducible oncogenic lncRNA [11] was also induced in our analysis but did not have a HIF-binding site close by. Hypoxic induction of non-annotated intergenic and antisense transcripts A proportion of reads did not map to an annotated transcript of any class. To define whether these represent previously non-annotated
EMBO reports
Vol 15 | No 1 | 2014
71
EMBO reports
Hypoxic regulation of the non-coding transcriptome
transcripts, they were first examined for contiguity, using the Cufflinks de novo transcript assembler. This identified 91 novel RNA transcripts with a median length of 9988 bp (range 269–38 2549 bp) many of which were upregulated in hypoxia (supplementary Tables S3 and S4) and were confirmed by qPCR (supplementary Fig S5). In keeping with assignment as bona fide non-annotated transcripts, reference to UCSC data (http://genome.ucsc.edu/cgi-bin/ hgTrackUi?hgsid=336849981&g=cpgIslandExt) revealed that 50 (55%) have a CpG Island within 1000 bp of the putative promoter. All 91 transcripts were also detected in the independently derived polyA+ RNA-seq dataset. Of the 91 non-annotated transcripts, 37 did not overlap with any annotated RefSeq gene, and were therefore classified as non-annotated intergenic transcripts (supplementary Table S3). All had a low coding potential [12], indicating that they belong to the lncRNA class of RNA. The remaining 54 overlapped with a previously documented RNA, but were expressed from the opposite DNA strand indicating that they are non-annotated antisense transcripts (NATS) (supplementary Table S4). Anti-sense transcripts may act in cis to regulate expression of the overlapping sense transcript [13]. Interestingly, the majority were upregulated by hypoxia, induction being commonly associated either with counterregulation (e.g. HIF-1a, TBX2) or with co-regulation of the overlapping sense transcript (e.g. SPAG4, CELSR2) (supplementary Fig S6). Thus, hypoxia-induced antisense transcripts are allied to both induction and repression in cis adding further tiers to regulation of the transcriptome by hypoxia. Transcriptional regulation of long non-coding RNAs by HIF As HIF is a major transcriptional regulator of the mRNA response to hypoxia [14], we next examined the role of HIF in the regulation of other classes of RNA. As a first step we examined for spatial association between each class of transcript and HIF-binding sites, using well validated, previously published HIF-1- and HIF-2-binding sites in MCF-7 cells [14]. That work considered only protein coding genes and reported a non-random distribution of HIF binding across the genome. Although HIF-binding sites were strongly enriched at the promoters of coding genes, many HIF-binding sites were observed to be remote from these promoters. In the light of extensive regulation of non-coding RNAs by hypoxia, we reanalysed the distribution of HIF-binding in relation to the location of all transcripts. The number of HIF-binding transcripts (defined as the closest promoter to each HIF-binding site) in each category broadly mirrored the total number in each class. Although the numbers in many classes were too small to permit firm conclusions, no enrichment of HIF binding was observed in the vicinity of classes of RNA showing global downregulation (piwiRNAs, snRNAs and tRNAs) suggesting that they are unlikely to be directly regulated by HIF. In contrast, HIF binding was enriched in the vicinity of both mRNAs and lncRNAs (Fig 2A) with HIF-2 binding a slightly higher proportion of lncRNAs than HIF-1. Nevertheless, despite promoter enrichment, the majority of HIF-1- and HIF-2-binding sites (50% and 70% respectively) were still found to lie in excess of 2.5-kb from the promoter of an expressed transcript of any class (Fig 2B and C). To assess the functional role of HIF in the regulation of different classes of RNA by hypoxia, we next examined for association of HIF binding with regulation by hypoxia across all classes of RNA using gene set enrichment analysis (GSEA) [15]. All expressed transcripts were pooled and ranked according to their hypoxic induction. HIF
72
EMBO reports
Vol 15 | No 1 | 2014
Hani Choudhry et al
binding was strongly enriched amongst transcripts upregulated (but not downregulated) by hypoxia (Fig 2D). Only mRNA, lncRNA, miRNA and non-annotated transcripts were present in the coreenrichment group of HIF-binding hypoxia-inducible transcripts. GSEA analysis by subgroup confirmed enrichment for mRNAs (Fig 2E) and lncRNAs (Fig 2F), but the number of miRNAs and non-annotated RNAs was too small to permit statistical analysis. Taken together, this establishes a non-random spatial association between HIF-binding sites and hypoxia-inducible transcripts for both mRNAs and lncRNAs, pointing to a major new role for HIF in the transcriptional regulation of lncRNAs. To pursue this further, we undertook polyA+ RNA-seq analysis of mRNA and lncRNA in hypoxic MCF-7 cells transfected with siRNAs targeting HIF-1a, HIF-2a or both isoforms simultaneously. Approximately 25% of all lncRNA transcripts detected in the ribosomal-depleted RNA-seq were detected in this analysis. Fig 3A portrays the results of HIF-a siRNA interventions on all HIF-binding lncRNA transcripts detected in this dataset. As expected, the heatmap confirms that many of these RNAs were upregulated by hypoxia and reveals a strong correlation with down regulation by siRNAs targeting one or both HIF-a isoforms, a result that was confirmed by GSEA (supplementary Fig S7). Individual datasets are illustrated for MALAT1 (downregulated by both HIF-1a and HIF-2a siRNA, Fig 3B) and NEAT1 (downregulated by HIF-2a but not by HIF-1a siRNA, Fig 3C). Thus, the presence of HIF binding, hypoxic upregulation and subsequent downregulation by the suppression of HIF-a defines lncRNAs as a new group of HIF target genes. HIF-dependent recruitment and release of promoter-paused RNApol2 in hypoxia To pursue the mechanism of transcriptional activation by hypoxia further, we next examined RNApol2 binding and histone H3K4me3 modification in normoxia and hypoxia. As expected, the promoterassociated abundance of each mark correlated with the basal transcript level. Strong peaks of RNApol2 and H3K4me3 signal were also seen at the putative transcriptional start sites (TSS) of the previously non-annotated transcripts consistent with assignment as bona fide genes (supplementary Fig S8) [16, 17]. Consideration of changes in RNApol2 binding in hypoxia revealed that for the majority of hypoxia-inducible genes, RNApol2 was already bound at the promoter in normoxic cells. Although there were exceptions, at the majority of promoter sites hypoxia did not increase promoter-bound RNApol2 substantially. This is illustrated in Fig 4A which depicts changes in RNApol2 ChIP signal in the vicinity of the promoters of the hundred most strongly hypoxia upregulated transcripts of any class. In hypoxia, little increase in loading of RNApol2 at the transcriptional start was seen; instead an increase in RNApol2 is observed across the body of the gene. Similar patterns of change in hypoxia were seen for both mRNAs and lncRNAs (supplementary Fig S9). Actively transcribing RNApol2 can recruit histone methyl transferase activity that trimethylates H3K4 either directly [18] or through p300 [19]. The distribution of H3K4me3 showed a very similar response to that of RNApol2, with little change in signal at the promoter in hypoxia, but an increase downstream of the transcriptional start site, consistent with the hypothesis that RNApol2 recruits histone methyl transferases that trimethylate H3K4 (Fig 4B). In contrast, consideration of DNAse hypersensitivity analyses in normoxic and hypoxic
ª 2014 The Authors
Hani Choudhry et al
EMBO reports
Hypoxic regulation of the non-coding transcriptome
D
A
E
B
F
C
G
Figure 2. HIF binding upregulates both the coding and non-coding transcriptome. A–C she proportion of HIF-1a and HIF-2a binding sites that are closest to transcribed loci of each RNA class (A). The distribution of (B) HIF-1a and (C) HIF-2a binding sites around the transcriptional start site of the nearest expressed gene irrespective of class (black bars). The grey bars show the distribution when only active mRNA genes are used (analogous to previous analyses). D–F GSEA analysis against fold-regulation by hypoxia for (D) all HIF-binding transcripts, (E) mRNA and (F) lncRNA. G Log2 fold-change of transcript abundance in hypoxia, ranked in order from highly upregulated to downregulated transcripts, from which the GSEA were derived.
MCF7 cells reported by the Encyclopedia of DNA Elements (ENCODE) consortium [17] revealed the expected hypersensitivity signals at these promoter regions, but did not reveal changes in hypoxia (Fig 4C). On close inspection the maximum RNApol2 signal was observed in the promoter-proximal region, 40–45 bp downstream of the transcriptional start. This suggests that promoter-proximal pausing is an important rate-limiting step for transcription of hypoxia-inducible genes and that reversal of pausing is a key step through which hypoxia regulates transcription [20]. To determine whether this is a feature of direct transcriptional regulation by HIF, as opposed to indirect effects of HIF or hypoxia itself on transcription, we performed similar analysis of RNApol2 distribution for genes that were downregulated by hypoxia and which (consistent with indirect regulation) do not show promoter enrichment of HIF binding. We also compared hypoxia-inducible changes in patterns of RNApol2 binding and H3K4
ª 2014 The Authors
trimethylation amongst groups of genes that were matched for the extent of upregulation by hypoxia, but were selected as being either the closest gene to a HIF-binding site (likely to represent direct transcriptional regulation) or were located >1 Mb from the nearest HIF-binding site (likely to represent indirect actions of HIF or hypoxia). Whilst a clear change in the ratio of promoter paused RNApol2 to that distributed along the gene body (travelling ratio) was observed for genes predicted to be direct HIF transcriptional targets this was not observed for genes predicted to be upregulated or downregulated indirectly by HIF or by hypoxia itself (supplementary Fig S10). Finally, ChIP-qPCR at a selection of these genes confirmed that the hypoxic induction of RNApol2 within the body of the gene was dependent on the presence of HIF-a subunits (Fig 4D–F). Amongst genes predicted to be direct transcriptional targets of HIF we observed a clear change in the distribution of travelling ratios in hypoxia, reflecting release of promoter-paused RNApol2
EMBO reports
Vol 15 | No 1 | 2014
73
EMBO reports
Hypoxic regulation of the non-coding transcriptome
A
Hani Choudhry et al
B
C
Figure 3. Direct transcriptional regulation oflncRNAs by HIF. A
Heat map showing fold regulation by hypoxia and by the indicated HIF siRNA for lncRNAs adjacent to HIF-binding site that are detected in the polyA RNA-seq analyses. B, C RNA-seq and HIF ChIP-seq genome browser tracks for the two most hypoxically upregulated HIF-binding lncRNAs: (B) MALAT1 and (C) NEAT1.
at most loci (Supplementary Fig S9E, F and Fig 4G). However, we also observed infrequent HIF target genes where RNApol2 was recruited de novo in hypoxia (Fig 4H). No association was detected between this type of behaviour and biological function. However, the promoters of these transcripts also show constitutive DNAse1 HSS signal indicating that other factors must limit normoxic RNApol2 binding at these sites. As DNA-seq cannot distinguish individual nucleosome positions, this might include mechanisms such as nucleosome positioning as well as DNA/histone methylation [21]. Taken together, our findings indicate that the regulation of gene expression by hypoxia is substantially more complex than has previously been considered. All RNA classes were regulated bi-directionally, our analysis identifying both general responses characteristic of particular types of RNA and high-amplitude effects on individual
74
EMBO reports
Vol 15 | No 1 | 2014
species. Both the non-random spatial association between HIF-binding sites and hypoxia-inducible genes and direct functional intervention by siRNA implicated HIF in the regulation of extensive networks of non-coding RNAs. RNApol2 profiling indicates that this direct regulation by HIF is mediated through activation of pre-bound RNApol2. These findings confirm and extend the role of HIF in the regulation of miRNAs and identify both intergenic and antisense lncRNAs as a new target class for the HIF transcriptional response. The ability of non-coding RNAs to regulate gene expression both in trans and in cis suggests the hitherto unrecognised involvement of HIF in a far broader network of gene regulation. Given the fundamental role of hypoxia in human pathology and emerging data on the importance of non-coding RNAs in cancer and other diseases, these findings open new avenues to the better understanding of human disease processes.
ª 2014 The Authors
Hani Choudhry et al
EMBO reports
Hypoxic regulation of the non-coding transcriptome
A
D
B
E
G
H
C
F
Figure 4. HIF predominantly regulates transcription through release of promoter-paused RNA-pol2. A–C Mean distribution of normoxic (red) and hypoxic (blue) RNApol2 binding at the 100 transcripts most upregulated by hypoxia. Inset shows expanded view at the TSS (FPKM = fragments per kilobase per million reads). The same plots for (B) H3K4me3 and (C) DNAse1 hypersensitivity are shown. D,E ChIP-qPCR analysis of RNApol2 within the body of (D) ALDOA and (E) NDRG1 gene shows suppression of the hypoxic induction by HIF-1a & HIF-2a siRNA (biological duplicates). F Immunoblot analysis of HIF-a levels. G,H RNA-seq and ChIP-seq tracks illustrating (G) release of promoter-paused RNApol2 and (H) de novo recruitment of RNApol2, despite constitutive DNAse1 hypersensitivity.
Materials and Methods Cell culture & HIF siRNA treatment Human MCF-7 breast carcinoma cells from the American Type Culture Collection (ATCC) were incubated for 24 h in an In vivo2 Hypoxia Work Station (Ruskinn Technology Ltd, Bridgend, UK) in an atmosphere of either normoxia (21% oxygen) or hypoxia (1% oxygen). Expression of HIF-1a and/or HIF-2a subunits was suppressed as previously described [22]. Ribosomal depleted RNA-seq Total RNA was prepared using the mirVana miRNA Isolation Kit (Ambion; Life Technologies Ltd, Paisley, UK) and treated with DNaseI (TURBO DNA-freeTM, Ambion). Ribosomal RNAs were depleted using RiboMinusTM Locked Nucleic Acid probes (Invitrogen; Life Technologies Ltd, Paisley, UK). The resulting RNA was fragmented and end-repaired using the Illumina directional protocol. The final cDNA library was purified using Agencourt
ª 2014 The Authors
AMPure beads (Beckman Coulter UK Ltd, High Wycombe, UK) optimised to retain small RNAs. PolyA+ selected RNA-seq Directional PolyA+ RNA libraries were prepared using the ScriptSeqTM v2 RNA-Seq kit (Epicentre, Madison, WI, USA) and non-directional libraries using the TruSeqRNA Sample Prep Kit (Illumina, San Diego, CA, USA). ChIP-seq Chromatin immunoprecipitation was performed as previously described [14]. In addition to HIF ChIPs, antibodies to H3K4me3 (Cell Signalling, Danvers, MA, USA, #9751) and RNApol2 (Santa Cruz Biotechnology, Dallas, TX, USA, #sc-899) were used. Libraries were prepared using the Illumina ChIP-seq kit. High throughput sequencing All libraries were sequenced using the GAII or HiSeq platforms (Illumina) according to the manufacturers protocol. Raw and mapped
EMBO reports
Vol 15 | No 1 | 2014
75
EMBO reports
Hypoxic regulation of the non-coding transcriptome
data for RNA-seq and ChIP-seq is available at EMBL-EBI Array Express (E-MTAB-1994, E-MTAB-1995).
8.
Hani Choudhry et al
Wollenick K, Hu J, Kristiansen G, Schraml P, Rehrauer H, BerchnerPfannschmidt U, Fandrey J, Wenger RH, Stiehl DP (2012) Synthetic transactivation screening reveals ETV4 as broad coactivator of hypoxia-
Bioinformatic analysis After trimming adapter sequences using the FASTX-tool kit (http:// hannonlab.cshl.edu/fastx_toolkit/), RNA-seq reads were sequentially mapped to RNA databases and expression analysis performed using the edgeR (Bioconductor, Seattle, WA, USA) package. Remaining unmapped reads were examined for non-annotated reads using the Cufflinks de novo transcript assembler (http://cufflinks.cbcb. umd.edu). ChIP-seq databases were analysed as previously described [14].
inducible factor signaling. Nucleic Acids Res 40: 1928 – 1943 9.
Clemson CM, Hutchinson JN, Sara SA, Ensminger AW, Fox AH, Chess A, Lawrence JB (2009) An architectural role for a nuclear noncoding RNA: NEAT1 RNA is essential for the structure of paraspeckles. Mol Cell 33: 717 – 726
10.
Schmidt LH, Spieker T, Koschmieder S, Schaffers S, Humberg J, Jungen D, Bulk E, Hascher A, Wittmer D, Marra A, Hillejan L, wiebe K, Berdel WE, Wiewrodt R, Muller-Tidow C (2011) The long noncoding MALAT-1 RNA indicates a poor prognosis in non-small cell lung cancer and induces migration and tumor growth. J Thorac Oncol 6: 1984 – 1992
Statistical analysis All assays were performed in duplicate or triplicate. Statistical analyses using two-tailed t-tests were performed in SPSS (IBM Corp, New York, NY, USA) and R (http://www.R-project.org).
11.
Matouk IJ, Mezan S, Mizrahi A, Ohana P, Abu-Lail R, Fellig Y, Degroot N, Galun E, Hochberg A (2010) The oncofetal H19 RNA connection: hypoxia, p53 and cancer. Biochim Biophys Acta 1803: 443 – 451
12.
Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G (2007) CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res 35:
Supplementary information for this article is available online: http://embor.embopress.org
W345 – 349 13.
Acknowledgements
epigenetic complexes that regulate transcription in human cells. Epige-
This work was funded by the Wellcome Trust (088182/Z/09/Z, 078333/Z/05/Z and WT091857MA), the Higher Education Funding Council for England, Cancer
netics 4: 296 – 301 14.
Research UK (A16016), the Ludwig Institute for Cancer Research, the Interdisci-
Schodel J, Oikonomopoulos S, Ragoussis J, Pugh CW, Ratcliffe PJ, Mole DR (2011) High-resolution genome-wide mapping of HIF-binding sites by ChIP-seq. Blood 117: e207 – 217
plinary Center of Clinical Research at the University of Erlangen-Nuremberg (TP J31) and King Abdulaziz University, Ministry of High Education for Saudi
Morris KV (2009) Long antisense non-coding RNAs function to direct
15.
Arabia.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP (2005) Gene set enrichment analysis: a knowledge-based approach for inter-
Author contributions
preting genome-wide expression profiles. Proc Natl Acad Sci USA 102:
HC, JS, ALH, PJR, IR, DRM designed experiments; HC, JS, SG, CC performed experiments; HC, SO, DRM performed data analysis; HC, PJR, IR, DRM wrote
15545 – 15550 16.
the manuscript.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K (2007) High-resolution profiling of histone methylations
Conflict of interest The authors declare that they have no conflict of interest.
in the human genome. Cell 129: 823 – 837 17.
Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis CA, Doyle F, Epstein CB, Frietze S, Harrow J, Kaul R, Khatun J, Lajoie BR, Landt SG, Lee BK, Pauli F, Rosenbloom KR, Sabo P, Safi A, Sanyal A, Shoresh N, et al (2012)
References
An integrated encyclopedia of DNA elements in the human genome. Nature 489: 57 – 74
1.
Semenza GL (2012) Hypoxia-inducible factors in physiology and medi-
18.
cine. Cell 148: 399 – 408 2. 3.
of Set1 histone methylase by elongating Pol II provides a localized
Kaelin WG Jr, Ratcliffe PJ (2008) Oxygen sensing by metazoans: the cen-
mark and memory of recent transcriptional activity. Mol Cell 11:
tral role of the HIF hydroxylase pathway. Mol Cell 30: 393 – 402
709 – 719
Ortiz-Barahona A, Villar D, Pescador N, Amigo J, del Peso L (2010) Genome-
19.
McGinty RK, Fernandez JP, Muir TW, Roeder RG (2013) SET1 and p300
genes by a probabilistic model integrating transcription-profiling data
Act Synergistically, through Coupled Histone Modifications, in Transcrip-
McCormick R, Buffa FM, Ragoussis J, Harris AL (2010) The role of hypoxia
tional Activation by p53. Cell 154: 297 – 310 20.
Adelman K, Lis JT (2012) Promoter-proximal pausing of RNA polymerase
21.
Zentner GE, Henikoff S (2012) Surveying the epigenomic landscape, one
22.
Elvidge GP, Glenny L, Appelhoff RJ, Ratcliffe PJ, Ragoussis J, Gleadle JM
regulated microRNAs in cancer. Curr Top Microbiol Immunol 345: 47 – 70 5.
Esteller M (2011) Non-coding RNAs in human disease. Nat Rev Genet 12:
II: emerging roles in metazoans. Nat Rev Genet 13: 720 – 731
861 – 874
76
Tang Z, Chen WY, Shimada M, Nguyen UT, Kim J, Sun XJ, Sengoku T,
wide identification of hypoxia-inducible factor binding sites and target and in silico binding site prediction. Nucleic Acids Res 38: 2332 – 2345 4.
Ng HH, Robert F, Young RA, Struhl K (2003) Targeted recruitment
base at a time. Genome Biol 13: 250
6.
Prensner JR, Chinnaiyan AM (2011) The emergence of lncRNAs in cancer biology. Cancer Discov 1: 391 – 407
(2006) Concordant regulation of gene expression by hypoxia and 2-oxo-
7.
Yang L, Duff MO, Graveley BR, Carmichael GG, Chen LL (2011) Genome-
glutarate-dependent dioxygenase inhibition: the role of HIF-1alpha, HIF-
wide characterization of non-polyadenylated RNAs. Genome Biol 12: R16
2alpha, and other pathways. J Biol Chem 281: 15215 – 15226
EMBO reports
Vol 15 | No 1 | 2014
ª 2014 The Authors
