Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Research Article
Transparent Process
Open Access

Multi-level interaction between HIF and AHR transcriptional pathways in kidney carcinoma

Véronique N Lafleur, Silvia Halim, Hani Choudhry, View ORCID ProfilePeter J Ratcliffe, View ORCID ProfileDavid R Mole  Correspondence email
Véronique N Lafleur
1https://ror.org/052gg0110NDM Research Building, University of Oxford, Old Road Campus, Oxford, UK
Roles: Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Silvia Halim
1https://ror.org/052gg0110NDM Research Building, University of Oxford, Old Road Campus, Oxford, UK
Roles: Data curation, Software, Formal analysis, Validation, Visualization, Methodology, Writing—original draft, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hani Choudhry
2Department of Biochemistry, Faculty of Science, Center of Innovation in Personalized Medicine, King Fahd Center for Medical Research, King Abdulaziz University, Jeddah, Saudi Arabia
Roles: Conceptualization, Supervision, Funding acquisition, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter J Ratcliffe
3https://ror.org/052gg0110Ludwig Institute for Cancer Research, University of Oxford, Old Road Campus, Oxford, UK
4The Francis Crick Institute, London, UK
Roles: Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Peter J Ratcliffe
David R Mole
1https://ror.org/052gg0110NDM Research Building, University of Oxford, Old Road Campus, Oxford, UK
Roles: Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David R Mole
  • For correspondence: david.mole@ndm.ox.ac.uk
Published 1 February 2023. DOI: 10.26508/lsa.202201756
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

Hypoxia-inducible factor (HIF) and aryl hydrocarbon receptor (AHR) are members of the bHLH-PAS family of transcription factors that underpin cellular responses to oxygen and to endogenous and exogenous ligands, respectively, and have central roles in the pathogenesis of renal cancer. Composed of heterodimers, they share a common HIF-1β/ARNT subunit and similar DNA-binding motifs, raising the possibility of crosstalk between the two transcriptional pathways. Here, we identify both general and locus-specific mechanisms of interaction between HIF and AHR that act both antagonistically and cooperatively. Specifically, we observe competition for the common HIF-1β/ARNT subunit, in cis synergy for chromatin binding, and overlap in their transcriptional targets. Recently, both HIF and AHR inhibitors have been developed for the treatment of solid tumours. However, inhibition of one pathway may promote the oncogenic effects of the other. Therefore, our work raises important questions as to whether combination therapy targeting both of these pro-tumourigenic pathways might show greater efficacy than targeting each system independently.

Introduction

The basic helix–loop–helix (bHLH) PER-ARNT-SIM (PAS) family of proteins are a closely related group of transcription factors that have important roles in both the physiological adaptation to environmental stimuli and in the pathogenesis of cancer (1, 2, 3). Although each has its own repertoire of transcriptional targets, their operation as heterodimers that use common dimerization partners within distinct subfamilies raises the possibility of complex interactions.

Notably, the hypoxia-inducible factor (HIF) responds to the hypoxic tumour microenvironment, which is a common feature of many solid cancers, in addition to being directly activated as a consequence of mutation in oncogenic and tumour suppressor pathways (4, 5). Dysregulation of HIF is linked to transcriptional programs with key roles in tumourigenesis and has direct implications for patient prognosis (6, 7). The most direct and profound up-regulation of the HIF pathway occurs in clear cell renal cell carcinoma (ccRCC), the most prevalent form of kidney cancer (8). This is associated with bi-allelic inactivation of the von Hippel Lindau (VHL) tumour suppressor, which forms part of an E3 ubiquitin ligase complex that targets HIF-α subunits (HIF-1α and HIF-2α; also known as EPAS-1) for proteasomal degradation in the presence of oxygen (9). As a result of VHL loss, physiological HIF-α degradation is disrupted, leading to constitutive activation of HIF. After stabilisation, the HIF-α proteins interact with HIF-1β (also known as ARNT; aryl hydrocarbon receptor nuclear translocator) to form functional HIF-1 and HIF-2 heterodimers that bind to chromatin (10). Despite a common consensus DNA-binding motif, HIF-1 and HIF-2 manifest distinct and overlapping binding patterns and transcriptional outputs and have opposing functions in ccRCC biology (11, 12, 13). Indeed, although HIF-2 is pro-tumourigenic, with its targeting recently approved as a treatment for VHL syndrome-associated tumours (14), HIF-1 is generally considered a tumour suppressor in ccRCC (15). Consequently, the overall effect of HIF activation on ccRCC tumourigenesis results from fine-tuning this balance between the HIF-1 and HIF-2 transcriptional outputs (16, 17, 18, 19).

The aryl hydrocarbon receptor (AHR) is a ligand-activated bHLH-PAS protein that acts as a mediator of metabolic responses to environmental contaminants, dietary compounds, and endogenous metabolites (20). Many of these contaminants, including cigarette smoke, polycyclic aromatic hydrocarbons, and halogenated aromatic hydrocarbons, are known carcinogens, with AHR being required for polycyclic aromatic hydrocarbon– and halogenated aromatic hydrocarbon–induced toxicities (2, 21, 22, 23, 24, 25). Like HIF, AHR can also be modulated by the tumour microenvironment and has critical roles in cancer progression, tumour aggressiveness, patient prognosis, and in tumour immunity (26, 27). AHR also dimerizes with HIF-1β/ARNT, and current evidence suggests that HIF and AHR compete directly for binding to HIF-1β/ARNT (28). However, observations also point to a more complex and context-dependent crosstalk that includes gene-specific cooperative events (29).

As an excretory organ, the kidneys are highly exposed to carcinogens and AHR has been found to be up-regulated in ccRCC tumours, in which it is associated with pathological tumour stage, histological grade, and poor prognosis (30). This raises the possibility of a consequential interplay between the AHR and HIF pathways in ccRCC. However, it remains unclear how the context-dependent competition between HIF and AHR directly impacts their chromatin binding abilities genome wide and whether it can involve locus-specific regulation or cooperation. To better understand these events and their integration for pathway outputs in the context of ccRCC, we have used chromatin immunoprecipitation coupled to next-generation DNA sequencing (ChIP-seq) to examine pan-genomic patterns of HIF-1α, HIF-2α, AHR, and HIF-1β/ARNT binding under co-activating stimuli and have linked these findings to gene expression changes analysed by RNA-sequencing. Our work reveals direct competition between HIF-α and AHR proteins for interaction with HIF-1β/ARNT, but that the symmetry of this competition is sensitive to the relative stoichiometry of the transcription factors. In ccRCC cells, we observed bi-directional competition between HIF and AHR that leads to a widespread antagonism of chromatin binding and to global downstream changes in gene expression. Strikingly, this competition is diminished at regions bound by both HIF and AHR, where binding of one transcription factor tends to reinforce binding of the other, and at a small set of genes that are direct transcriptional targets of both transcription factors. Overall, this reveals an unexpected complexity in the overlap and interaction between the HIF and AHR transcriptional programs. Finally, we provide evidence of this antagonistic relationship in ccRCC tumours, which display an overall negative correlation between HIF- and AHR-target gene expression. These findings shed light on the complex and direct crosstalk between the HIF and AHR transcriptional pathways, which encompasses both a widespread inhibitory interaction and specific cooperative events in ccRCC. Given the important role of HIF and AHR as regulators of cancer progression, therapeutic inhibitors have now been developed to target each pathway (14, 26). Our work, therefore, has important therapeutic implications because inhibition of one pathway may promote the oncogenic effects of the other. This raises the possibility that targeting both pathways simultaneously may lead to enhanced synthetic lethality.

Results

Uni- and bi-directional competition between HIF-α and AHR in different cell lines

We first surveyed a panel of cell lines for crosstalk between the HIF and AHR transcriptional pathways, using the expression of canonical target genes for each transcription factor. To start, we determined optimal conditions for target gene induction using RCC4 renal cancer cells re-constituted with WT VHL (RCC4+VHL) to restore oxygen-dependent HIF-α regulation. Firstly, we found maximal or near maximal induction of HIF-target genes after incubation in 0.5% hypoxia for 6 h that was maintained at 18 h, when compared with normoxia (ambient oxygen) (Fig S1A). Secondly, to activate AHR, we used two endogenous AHR ligands—ITE, an indole-based dietary ligand, and FICZ, a tryptophan photo-metabolite. Incubation of RCC4+VHL cells with 100 nM ITE resulted in a time course of AHR-target gene induction which was similar to that of HIF, whereas incubation of these cells with 10 nM FICZ resulted in a more transient AHR-target gene induction that was maximal at 6 h and waning by 18 h, when compared with DMSO (vehicle control) (Fig S1B). Therefore, we chose to examine the 6-h time point for these and subsequent experiments.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Induction of HIF- and AHR-target gene expression.

(A, B) RT-qPCR analysis of HIF-target genes in RCC4+VHL cells exposed to normoxia or 0.5% hypoxia for 6 or 18 h (B) RT-qPCR analysis of AHR-target genes in RCC4+VHL cells exposed to DMSO, 100 nM ITE or 10 nM FICZ for 6 or 18 h. (C) Immunoblots showing comparative levels of AHR, HIF-1α, HIF-2α, and HIF-1β/ARNT (CST-5537) protein in cell lines exposed to 0.5% hypoxia for 18 h (β-Actin and Coomassie blue stain used as loading controls). (D) RT-qPCR analysis of AHR-target genes in RCC4+VHL cells exposed to control conditions, 10 nM FICZ, 0.5% hypoxia or both for 6 h. (E) Immunoblots showing comparative levels of AHR, HIF-1α, HIF-2α and HIF-1β/ARNT protein in RCC4+VHL cells exposed to either 0.5% hypoxia or 50 μM FG4592 for the specified duration. (F, G) RT-qPCR analysis of AHR-target genes in (F) RCC4+VHL cells exposed to control conditions, 100 nM ITE, 50 μM FG4592, or both for 6 h, and (G) RCC4+VHL cells exposed to control conditions, 10 nM FICZ, 50 μM FG4592, or both for 6 h. (H) Immunoblots showing HIF-1α, HIF-2α, and HIF-1β (CST-5537) protein in HKC8 in parental and HIF-a KO cell lines exposed to normoxia or 0.5% hypoxia for 6 h (β-actin used as loading control). (I, J, K, L, M, N) RT-qPCR analysis of (I) NDRG1, (J) PFKFB4, and (K) INSIG2 in WT, HIF-1α KO, HIF-2α KO, or HIF-1α/2α double KO HKC8 cells exposed to normoxia or 0.5% hypoxia for 6 h RT-qPCR analysis of HIF-target genes in (L) RCC4+VHL cells exposed to control conditions, 10 nM FICZ, 0.5% hypoxia, or both for 6 h, (M) RCC4+VHL cells exposed to control conditions, 100 nM ITE, 50 μM FG4592, or both for 6 h, and (N) RCC4+VHL cells exposed to control conditions, 10 nM FICZ, 50 μM FG4592, or both for 6 h. Data information: data in RT-qPCR analyses (A, B, D, F, G, I, J, K, L, M, N) are presented as mean ± SD of three biological replicates. *P ≤ 0.05, **P ≤ 0.01 (paired t-test).

In four of the cell lines analysed for crosstalk—RCC4+VHL, HKC8, HepG2, and PC3—incubation in hypoxia resulted in inhibition of the canonical AHR-target gene CYP1B1 (Fig 1A). No significant effect was seen in 786O+VHL cells, which express only the HIF-2α isoform of HIF. Conversely, hypoxic induction of the canonical HIF-target gene NDRG1 was largely unaffected by treatment with ITE in HKC8, HepG2, and PC3 cells (Fig 1B). However, modest suppression of hypoxic NDRG1 was observed in RCC4+VHL and 786O+VHL cells incubated with ITE. These cell lines have particularly high levels of AHR compared with HIF-α, suggesting that the stoichiometry may favour bi-directional crosstalk (Fig S1C).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Context-specific crosstalk between HIF and AHR.

(A, B) RT-qPCR analysis of (A) CYP1B1 and (B) NDRG1 mRNA expression in RCC4+VHL, 786O+VHL, HKC8, HepG2, and PC3 cells subjected to incubation for 6 h in control conditions, 0.5% hypoxia, 100 nM ITE, or both. (C, D) RT-qPCR analysis of (C) AHR-target genes and (D) HIF-target genes in RCC4+VHL cells subjected to the same conditions as above. (E) RT-qPCR analysis of CYP1B1 mRNA expression in control HKC8 cells and HKC8 cells deficient in HIF-1α, HIF-2α, or both. Data information: data are presented as mean ± SD of four (A, B) or three (C, D, E) biological replicates. *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001 (paired t-test).

To explore this further, we examined additional HIF- and AHR-target genes, and other stimuli, in RCC4+VHL cells. In addition to CYP1B1, the ITE-mediated induction of ALDH1A3, BMF, and ASB2 was suppressed in hypoxia (Fig 1C), whereas the induction of CYP27A1 was unaffected. Similarly, FICZ-mediated induction of CYP1B1, BMF, and ASB2 was also suppressed in hypoxia (Fig S1D). Furthermore, the HIF-prolyl hydroxylase inhibitor FG4592 induced HIF-1α and HIF-2α independently of hypoxia (Fig S1E) and resulted in down-regulation of all five AHR-target genes in ITE-treated cells and 4/5 genes in FICZ-treated cells (Fig S1F and G). This strongly implicates HIF rather than non-HIF pathways in the hypoxic suppression of AHR-target genes. However, to test this more formally, we examined the effect of hypoxia on the ITE-mediated induction of CYP1B1 in HKC8 cells deficient in HIF-1α, HIF-2α, or both isoforms (Fig 1E). Effective knockout of these transcription factors was confirmed by abolition of hypoxic induction of HIF protein and HIF-target genes (Fig S1H–K). Cells lacking both HIF-1α and HIF-2α isoforms did not show any effect of hypoxia on ITE-mediated CYP1B1 expression, confirming its dependence on HIF. Furthermore, in HKC8 cells, the hypoxic suppression of CYP1B1 expression appeared largely mediated by HIF-1α because cells lacking this isoform exhibited no hypoxic regulation of CYP1B1, whereas loss of HIF-2α resulted in a more modest diminution in the effect of hypoxia.

Conversely, in addition to NDRG1, modest ITE-mediated suppression of the HIF-target genes PFKFB4 and INSIG2 was observed in hypoxia (Fig 1D), whereas FICZ led to down-regulation of NDRG1 alone (Fig S1L). Furthermore, when HIF was stimulated with FG4592, neither ITE nor FICZ had a significant impact upon the induction of HIF-target genes (Fig S1M and N). Thus, although the crosstalk between HIF and AHR in RCC4+VHL cells operates in both directions, under the conditions used, HIF has a larger effect on AHR-target genes than AHR does on HIF-target genes.

Direct competition between HIF-α and AHR for HIF-1β/ARNT dimerization in ccRCC cells

We next explored the mechanisms underlying the crosstalk between the HIF and AHR transcriptional pathways in RCC4+VHL cells. The antagonistic relationship between HIF and AHR has previously been linked to competition for the common dimerization partner HIF-1β/ARNT, which may be limiting under specific circumstances (28). Therefore, we first assessed the competition between the HIF-α and AHR proteins for binding to HIF-1β/ARNT in RCC4+VHL cells. We initially determined optimal conditions for HIF-α and AHR induction in these cells. Firstly, exposure of RCC4+VHL cells to 6 h of 0.5% hypoxia resulted in maximal stabilisation of both HIF-1α and HIF-2α proteins (Fig S2A). Secondly, upon ligand binding, cytoplasmic AHR rapidly translocates to the nucleus and binds chromatin in complex with HIF-1β/ARNT. It subsequently undergoes protease-mediated degradation as part of a negative feedback mechanism (31). We therefore used translocation of AHR to the nucleus as a marker of its activation and observed rapid and transient increases in nuclear levels in response to 100 nM ITE and 10 nM FICZ that were maximal at 30 min (Fig S2B). Cells were therefore incubated for 6 h in either 0.5% hypoxia or normoxia, followed by treatment with ITE, FICZ, or DMSO (control) for the final 30 min (Fig 2A and B). These conditions resulted in strong induction of binding of each transcription factor to chromatin (Fig S2C and D).

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Analysis of AHR, HIF-1α, HIF-2α, and HIF-1β/ARNT protein levels in RCC4 cells.

(A) Immunoblots showing comparative levels of AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (CST-5537), and NDRG1 protein in RCC4+VHL and RCC4+VA cells exposed to 0.5% hypoxia for the specified durations (β-actin and Coomassie blue stain used as loading controls). (B) Immunoblots showing comparative levels of AHR, HIF-1β/ARNT (CST-5537) and CYP1B1 protein in nuclear and cytoplasmic extracts from RCC4+VHL cells exposed 100 nM ITE or 10 nM FICZ for the specified durations (α-tubulin, PolII and Coomassie blue stain used as loading controls). (C, D) ChIP-qPCR analysis of HIF-1β/ARNT chromatin binding (C) at the AHR binding site at the CYP1B1 gene locus and (D) at the HIF binding site at the NDRG1 locus in RCC4+VHL cells exposed to the specified conditions.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. HIF and AHR compete for HIF-1β/ARNT.

RCC4+VHL cells were subjected to the indicated conditions and nuclear extracts were prepared and immunoprecipitated using antibodies directed against HIF-1β/ARNT (rabbit pAb NB100-110). (A) Input samples and immunoprecipitates were then immunoblotted for AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (mouse mAb NB100-124), and HDAC2. (B) Quantitation of immunoprecipitation results from three biological replicates presented as mean ± SD *P ≤ 0.05, **P ≤ 0.01 (paired t-test). *ns, non-specific band.

Nuclear and cytoplasmic extracts were then prepared from these cells and immunoprecipitated using antibodies directed against HIF-1β/ARNT (Figs 2A and B and S3A). Although AHR, HIF-1α, and HIF-2α were all present in both the nucleus and the cytoplasm, the fraction bound to HIF-1β/ARNT was almost exclusively nuclear (Fig 2A and B). This is consistent with the previously described role for HIF-1β/ARNT as a nuclear translocator and a chromatin binding partner. Dimerization of HIF-1α and HIF-2α with HIF-1β/ARNT was observed under hypoxic incubation, whereas dimerization of AHR with HIF-1β/ARNT was strongly induced by both ITE and FICZ treatments (Fig 2A and B). Strikingly, ITE and FICZ treatments both reduced HIF-1α-HIF-1β/ARNT and HIF-2α-HIF-1β/ARNT dimer formation in hypoxic cells. Conversely, hypoxic incubation repressed AHR-HIF-1β/ARNT dimer formation in both ITE- and FICZ-treated cells even more strongly. Comparable inhibition of the AHR-HIF-1β/ARNT dimer was observed after stabilisation of HIF-α proteins by the HIF-prolyl hydroxylase inhibitor FG4592 (Fig S3B). Interestingly, relocation of AHR from nuclear to cytoplasmic cell fractions was observed under HIF-stabilising conditions, even in the presence of ITE and FICZ (Figs 2A and B and S3A and B). Thus, AHR nuclear localisation appears to correlate with its ability to interact with HIF-1β/ARNT and therefore to be repressed by competition with HIF-α. This suggests that the competition between HIF-α and AHR not only alters the ability of AHR to form an active transcription factor in the nucleus but also increases its cytoplasmic localisation, which may indirectly affect non-canonical AHR functions. Furthermore, re-oxygenation of cells after combined hypoxic and ITE treatments led to rapid degradation of the HIF-α proteins and to complete restoration of the AHR-HIF-1β/ARNT dimer (Fig S3C). Thus, the HIF-α and AHR proteins are able to compete reciprocally for HIF-1β/ARNT in RCC4+VHL cells, with HIF-α stabilisation having a greater impact on AHR binding than vice versa. Conversely, in HepG2 cells, in which crosstalk at the level of gene expression was found to be uni-directional, with HIF induction suppressing AHR-target gene expression but not vice versa, induction of HIF-α in hypoxia inhibited binding of AHR to HIF-1β/ARNT but induction of AHR did not inhibit HIF dimerization (Fig S3D). These cell type differences suggest that the bi-directionality or uni-directionality of the crosstalk between HIF and AHR is unlikely to be a consequence of inherent properties of HIF-α and AHR and more likely to result from differences in stoichiometry of the two transcription factors.

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. IP-IB analysis of HIF-1β/ARNT interacting proteins.

(A) RCC4+VHL cells were subjected to the indicated conditions, and cytoplasmic extracts were prepared and immunoprecipitated using antibodies directed against HIF-1β/ARNT (rabbit pAb NB100-110). Input samples and immunoprecipitates were then immunoblotted for AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (mouse mAb NB100-124), and α-tubulin. (B) RCC4+VHL cells were subjected to the indicated conditions and nuclear extracts were prepared and immunoprecipitated using antibodies directed against HIF-1β/ARNT (NB100-110). Cytoplasmic extracts, input nuclear samples, and immunoprecipitates were then immunoblotted for AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (NB100-124), and HDAC2. (C) RCC4+VHL cells were exposed to hypoxia and/or ITE followed by re-oxygenation in 21% oxygen for 30 min. Nuclear extracts were prepared and immunoprecipitated using antibodies directed against HIF-1β/ARNT (NB100-110). Input samples and immunoprecipitates were then immunoblotted for AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (NB100-124), and HDAC2. (D) HepG2 cells were subjected to the indicated conditions, and nuclear extracts were prepared and immunoprecipitated using antibodies directed against HIF-1β/ARNT (NB100-110). Input samples and immunoprecipitates were then immunoblotted for AHR, HIF-1α, HIF-2α, HIF-1β/ARNT (NB100-124), and β-actin. *ns, non-specific band.

Pan-genomic antagonism between HIF and AHR chromatin binding

Although competition between HIF-α and AHR for binding to HIF-1β/ARNT would likely result in widespread suppression of signaling in response to activation of the opposing factor, our analysis of target gene expression has shown heterogenous responses at individual target genes (Fig 1). Therefore, to further assess locus-specific crosstalk between HIF and AHR across the genome, we first performed ChIP-seq analysis of HIF-1α, HIF-2α, HIF-1β/ARNT, and AHR chromatin binding. Experiments were performed on RCC4+VHL cells incubated in either 0.5% hypoxia or normoxia for 6 h and treated with ITE, FICZ, or DMSO (control) for the final 30 min. All experiments were performed in duplicate and compared with input chromatin control samples. Only peaks identified by both the MACS and TPIC peak callers in both replicates of each condition were used for analysis.

Preliminary analysis was performed to examine the pairwise correlation between samples. The signal intensity at each peak (averaged between the two replicates) was plotted as a heat map (r2 = 1, dark green, r2 = 0, white), and samples were grouped using unsupervised hierarchical clustering (Fig S4). Three distinct clusters were identified, composed of samples from each of the HIF-1α, HIF-2α, and AHR ChIP-seq analyses. Samples from the HIF-1β/ARNT analyses clustered with the AHR analyses when cells were stimulated with ITE or FICZ alone but clustered with the HIF-1α analyses when cells were treated with hypoxia, even in the presence of ITE or FICZ. This indicates that HIF-1β/ARNT binding correlates more closely with HIF-1α binding than with AHR binding when both factors are present. Clustering of the hypoxic HIF-1β/ARNT analyses with HIF-1α rather than HIF-2α likely reflects the greater number of HIF-1 sites identified (n = 828) compared with HIF-2 sites (n = 401).

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Preliminary analysis of ChIP-seq data.

Canonical AHR, HIF-1, and HIF-2 binding sites (those binding HIF-1β/ARNT and either AHR, HIF-1α or HIF-2α in both ChIP-seq replicates) were identified from ChIP-seq analyses of RCC4+VHL cells exposed to combinations of either normoxia or hypoxia and either DMSO (control), ITE or FICZ. For each ChIP-seq analysis, the combined signal for both replicates was determined at each of the defined sites. The signal intensity for each ChIP-seq analysis was then correlated pairwise with each of the other ChIP-seq analyses across all these sites. R2 values are displayed as a heat map, and the ChIP-seq analyses are ordered according to unsupervised hierarchical clustering based on these correlations.

Next, canonical AHR binding sites (those with AHR and HIF-1β/ARNT signal in both replicates of a condition) were identified in normoxic cells treated with either DMSO, ITE or FICZ (Fig 3A). In total, 1,176 sites were identified in one or more condition. As expected, binding of AHR was strongly induced across these sites by both ITE and FICZ, with concomitant increases also seen for HIF-1β/ARNT signal (Fig S5A–K). Signal intensity for both AHR and HIF-1β/ARNT correlated very strongly between the two treatments (r2 = 0.96 and 0.94, respectively) with all sites having comparable signal with either ITE or FICZ (Fig S5L and M). However, under the conditions used, there was, on average, slightly stronger binding with ITE than with FICZ.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. HIF and AHR inhibit binding of each other to chromatin in trans.

(A) Canonical AHR binding sites (those binding both AHR and HIF-1β/ARNT in both ChIP-seq replicates) were identified in RCC4+VHL cells treated with DMSO (control), ITE, or FICZ. The Venn diagram shows sites detected under the different stimuli. (B) Canonical HIF-1 binding sites in hypoxic RCC4+VHL cells (binding both HIF-1α and HIF-1β/ARNT in both biological replicates) and canonical HIF-2 binding sites (binding both HIF-2α and HIF-1β/ARNT in both biological replicates) were identified. The Venn diagram shows the number of sites bound by one or both factors. (C, D, E) Most significantly enriched sequence logos in DREME analysis (32) of (C) canonical AHR unique sites, (D) canonical HIF-1 unique sites, and (E) canonical HIF-2 unique sites. (F, G) Scatter plot and (G) box-and-whisker plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + ITE versus normoxia + ITE. (H, I) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (J, K) Scatter plot and (K) box-and-whisker plot showing normalised HIF-1α signal intensity at canonical HIF-1 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus hypoxia + ITE. (L, M) The same analysis for HIF-1β/ARNT signal at canonical HIF-1 sites. (N, O) Scatter plot and (O) box-and-whisker plot showing normalised HIF-2α signal intensity at canonical HIF-2 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus hypoxia + ITE. (P, Q) The same analysis for HIF-1β/ARNT signal at canonical HIF-2 sites. Data information: P-values in box-and-whisker plots (G, I, K, M, O, Q) by Wald test with correction using Benjamini–Hochberg procedure.

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. Induction of AHR and HIF chromatin binding in response to single stimuli.

(A, B) Scatter plot and (B) box-and-whisker plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells treated with normoxia + ITE versus normoxia + DMSO. (C, D) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (E, F) Scatter plot and (F) box-and-whisker plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells treated with normoxia + FICZ versus normoxia + DMSO. (G, H) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (I, J, K) IGV tracks at the (I) CYP1B1, (J) ALDH1A3, and (K) BMF gene loci showing superimposed ChIP-seq signal for AHR and HIF-1β/ARNT signal in control (light colours) and ITE-stimulated (dark colours) RCC4+VHL cells. (L) Scatter plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells treated with normoxia + ITE versus normoxia + FICZ. (M) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (N, O) Scatter plot and (O) box-and-whisker plot showing normalised HIF-1α signal intensity at canonical HIF-1 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus normoxia + DMSO. (P, Q) The same analysis for HIF-1β/ARNT signal at canonical HIF-1 sites. (R, S) Scatter plot and (S) box-and-whisker plot showing normalised HIF-2α signal intensity at canonical HIF-2 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus normoxia + DMSO. (T, U) The same analysis for HIF-1β/ARNT signal at canonical HIF-2 sites. (V, W, X) IGV tracks at the (V) NDRG1, (W) PFKFB4, and (X) INSIG2 gene loci showing superimposed ChIP-seq signal for HIF-1α, HIF-2α, and HIF-1β/ARNT in normoxic (light colours) and hypoxic (dark colours) RCC4+VHL cells. Data information: P-values in box-and-whisker plots (B, D, F, H, O, Q, S, U) by Wald test with correction using Benjamini–Hochberg procedure.

Canonical HIF binding sites were identified in a similar way, using peaks that had either HIF-1α or HIF-2α signal together with HIF-1β/ARNT signal, in both replicates of a condition. In total, 983 canonical HIF binding sites were identified in hypoxic cells treated with DMSO (828 HIF-1 sites and 401 HIF-2 sites, of which 246 were present in both datasets) (Fig 3B). 101 of these sites were also detected in normoxia, suggesting low-level HIF activation under this condition. Consistent with this, no sites were detected in normoxia alone. As expected, HIF-1α and HIF-2α binding was strongly induced by hypoxia at canonical HIF-1 and HIF-2 sites, respectively, with concomitant increases also observed in HIF-1β/ARNT binding (Fig S5N–X). Structural analysis has shown that the HIF-1β/ARNT subunit of HIF and AHR binds a GTG DNA sequence (33). Consistent with this, both AHR and HIF binding sites were strongly enriched for a consensus NCGTG motif (Fig 3C–E). However, HIF-1 and HIF-2 strongly favoured the ACGTG variant (the hypoxia response element), whereas AHR preferentially recognised the GCGTG motif (the xenobiotic/dioxin response element).

We next investigated the interaction between HIF and AHR activation on the binding of the opposing factor to chromatin. Hypoxia significantly reduced binding of both AHR and HIF-1β/ARNT at canonical AHR sites in ITE-treated cells (Figs 3F–I and S6A–C). This was independent of signal intensity at the site. Hypoxia also had a weaker effect on AHR and HIF-1β/ARNT binding in FICZ-treated cells (Fig S6D–G). Taken together, these results indicate that hypoxia reduces AHR DNA binding in RCC4+VHL cells. Conversely, ITE treatment significantly reduced HIF-1α binding at canonical HIF-1 sites and HIF-2α binding at canonical HIF-2 sites with concomitant reductions in HIF-1β/ARNT binding (Figs 3J–Q and S6H–J). Again, weaker effects were seen when the cells were treated with FICZ (Fig S6K–R). Overall, AHR activation reduces both HIF-1 and HIF-2 binding in RCC4+VHL cells. Thus, we observe bi-directional competition between HIF-α and AHR at the level of chromatin binding in RCC4+VHL cells that reflects their bi-directional competition for HIF-1β/ARNT.

Figure S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S6. AHR and HIF chromatin binding in RCC4+VHL cells in response to combined FICZ and hypoxic stimuli.

(A, B, C) IGV tracks at the (A) CYP1B1, (B) ALDH1A3, and (C) BMF gene loci showing superimposed ChIP-seq signals for AHR and HIF-1β/ARNT in ITE-treated RCC4+VHL cells in normoxia (dark colours) and in ITE-treated RCC4+VHL cells in hypoxia (light colours). (D, E) Scatter plot and (E) box-and-whisker plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + FICZ versus normoxia + FICZ. (F, G) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (H, I, J) IGV tracks at the (H) NDRG1, (I) PFKFB4, and (J) INSIG2 gene loci showing superimposed ChIP-seq signals for HIF-1α, HIF-2α, and HIF-1β/ARNT in hypoxic RCC4+VHL cells treated with DMSO (dark colours) and hypoxic RCC4+VHL cells treated with ITE (light colours). (K, L) Scatter plot and (L) box-and-whisker plot showing normalised HIF-1α signal intensity at canonical HIF-1 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus hypoxia + FICZ. (M, N) The same analysis for HIF-1β/ARNT signal at canonical HIF-1 sites. (O, P) Scatter plot and (P) box-and-whisker plot showing normalised HIF-2α signal intensity at canonical HIF-2 binding sites (averaged across the two replicates) in RCC4+VHL cells treated with hypoxia + DMSO versus hypoxia + FICZ. (Q, R) The same analysis for HIF-1β/ARNT signal at canonical HIF-2 sites. Data information: P-values in box-and-whisker plots (E, G, L, N, P, R) by Wald test with correction using Benjamini–Hochberg procedure.

In clear cell renal cancer, HIF-α is constitutively stabilised irrespective of oxygen levels because of VHL inactivation. Because HIF activated in this way may differ from HIF stabilised by hypoxia, we examined the effect of this normoxic HIF activation on binding of AHR and HIF-1β/ARNT to canonical AHR sites by comparing signals in VHL WT RCC4+VHL cells and VHL-deficient RCC4+VA cells in normoxia (Fig S7). VHL loss led to a significant reduction in both AHR and HIF-1β/ARNT binding at canonical AHR sites in cells treated with ITE, comparable to that observed under hypoxia (Fig S7A–D). Conversely, in VHL-defective RCC4+VA cells, ITE significantly reduced binding of HIF-1α, HIF-2α, and HIF-1β/ARNT at canonical HIF binding sites (Fig S7E–L). Taken together, this indicates that the crosstalk between AHR and HIF on chromatin binding is a direct consequence of HIF activation rather than an off-target effect of hypoxia. It also confirms the relevance of our findings to kidney cancer, indicating that HIF activation by VHL loss can alter AHR binding and that AHR stimulation can reciprocally modulate the HIF response.

Figure S7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S7. AHR and HIF chromatin binding in RCC4+VHL and RCC4+VA cells in response to ITE stimuli.

(A, B) Scatter plot and (B) box-and-whisker plot showing normalised AHR signal intensity at canonical AHR binding sites (averaged across the two replicates) in RCC4+VHL cells versus RCC4+VA cells both treated with normoxia + DMSO. (C, D) The same analysis for HIF-1β/ARNT signal at canonical AHR sites. (E, F) Scatter plot and (F) box-and-whisker plot showing normalised HIF-1α signal intensity at canonical HIF-1 binding sites (averaged across the two replicates) in RCC4+VA cells treated with normoxia + DMSO versus normoxia + ITE. (G, H) The same analysis for HIF-1β/ARNT signal at canonical HIF-1 sites. (I, J) Scatter plot and (J) box-and-whisker plot showing normalised HIF-2α signal intensity at canonical HIF-2 binding sites (averaged across the two replicates) in RCC4+VA cells treated with normoxia + DMSO versus normoxia + ITE. (K, L) The same analysis for HIF-1β/ARNT signal at canonical HIF-2 sites. Data information: P-values in box-and-whisker plots (B, D, F, H, J, L) by Wald test with correction using Benjamini–Hochberg procedure.

Reduced competition between HIF and AHR at shared chromatin binding sites

Of the 1,933 canonical HIF and AHR binding sites identified, 221 (11%) were overlapping HIF and AHR binding sites. This overlap occurred at 19% of all AHR sites and was higher at HIF-2 sites (42% of sites) than at HIF-1 sites (18% of sites) (Fig 4A). Of note, canonical AHR binding sites also showed a similar promoter-distal distribution to canonical HIF-2 binding sites in contrast to the more promoter-proximal distribution observed for canonical HIF-1 binding sites (Fig 4B–D). This similarity in distribution may in part explain the greater overlap between AHR and HIF-2 binding sites.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. HIF and AHR promote binding of each other to chromatin in cis.

(A) Venn diagram showing the number of individual and overlapping canonical HIF-1, HIF-2, and AHR binding sites. (B, C, D) Bar charts showing the frequency distribution for (B) canonical AHR sites, (C) canonical HIF-1 sites, and (D) canonical HIF-2 sites with respect to the distance upstream or downstream from the closest transcriptional start site. (E) Box-and-whisker plot showing the log2 fold change in AHR signal intensity between hypoxia + ITE and normoxia + ITE for canonical AHR binding sites without HIF binding (Indiv Sites) and with HIF-1 binding (Shared Sites). (F) The same analysis for canonical AHR sites without HIF binding (Indiv Sites) and with HIF-2 binding (Shared Sites). (G) Box-and-whisker plot showing the log2 fold change in HIF-α signal intensity between hypoxia + ITE and hypoxia + DMSO for canonical HIF-1 binding sites without AHR binding (Indiv Sites) and with AHR binding (Shared Sites). (H) The same analysis for canonical HIF-2 sites without AHR binding (Indiv Sites) and with AHR binding (Shared Sites). (I, J) Box-and-whisker plots showing the distance between (I) AHR and HIF-1α and (J) AHR and HIF-2α peak summits (grey boxes). For comparison, the distance between AHR peak summits when stimulated with either ITE or FICZ is also shown (red boxes). Data information: P-values in box-and-whisker plots (E, F, G, H) by Wald test with correction using Benjamini–Hochberg procedure.

The hypoxia response element HIF-recognition motif (RCGTG) and the xenobiotic/dioxin response element AHR-recognition motif (GCGTG) are very similar. This raises the possibility that in addition to a global competition for HIF-1β/ARNT, HIF and AHR may compete for access to the same binding motif on chromatin at these overlapping sites. To dissect this further, we examined whether shared AHR and HIF sites were subjected to a different level of competition than sites that bound either transcription factor alone. Contrary to our hypothesis, hypoxia suppressed AHR binding at shared sites less than it did at sites that bound AHR alone. Specifically, the hypoxic repression of AHR binding at sites shared with HIF-1 (median log2FC: −1.199) or HIF-2 (median log2FC: −1.264) was significantly weaker than at individual sites (median log2FC: −1.458) (Fig 4E and F). Conversely, the ITE-dependent reduction in HIF-1α binding was also weaker at sites shared with AHR (median log2FC: −0.490) than at individual sites (median log2FC: −0.591) (Fig 4G). The difference in repression of HIF-2α binding at shared sites (median log2FC: −0.697) versus individual sites (median log2FC: −0.703) did not reach statistical significance (Fig 4H). Taken together, this suggests that rather than competing for the same binding motif, HIF and AHR positively reinforce chromatin binding of each other in cis when bound at overlapping sites to an extent that partially mitigates the competition for HIF-1β/ARNT.

To gain better insight into the potential mechanisms underlying this, we attempted to resolve the precise binding overlap between these factors at shared sites by comparing the peak summits. As a control, we established that the median distance between summits for AHR ChIP-seq signal when stimulated with either ITE or FICZ was just 19 bp. However, the median distance between AHR and HIF-1α summits was significantly larger at 39 and 38 bp (for ITE and FICZ, respectively), whereas the distance between AHR and HIF-2α summits was 62 and 59 bp (Fig 4I and J). Taken together, this suggests that even where AHR and HIF-α ChIP-seq peaks overlap, the two transcription factors generally bind to closely associated but distinct sites (Fig S8, e.g., loci), consistent with HIF and AHR binding to the two closely related but largely distinct recognition motifs ACGTG and GCGTG, respectively.

Figure S8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S8. AHR and HIF binding signal at co-bound sites.

(A, B, C) IGV tracks at the (A) PAX2, (B) LUCAT1 and (C) DDIT4 gene loci showing superimposed ChIP-seq signals for AHR and HIF-1β/ARNT in ITE-treated RCC4+VHL cells and for HIF-1α, HIF-2α, and HIF-1β/ARNT in hypoxic cells. The distance between the XRE and HRE motifs most proximal to peak summits is represented by dashed lines.

Contrasting antagonistic and cooperative activities of HIF and AHR in gene regulation

The observation of both general and locus-specific effects in the interaction between HIF and AHR transcription pathways suggested that there might be heterogenous consequences on target gene expression that go beyond general antagonism. Therefore, to relate changes in the DNA binding of HIF and AHR across the genome to downstream changes in gene expression, we performed RNA-sequencing on RCC4+VHL cells incubated in 0.5% hypoxia and/or treated with ITE for 6 h.

Gene set enrichment analysis (GSEA) revealed that AHR-bound genes—those closest to canonical AHR-only binding sites—were enriched among ITE-up-regulated genes (enrichment score [ES]: 0.98; P-value = 0.03) and also enriched among hypoxia-down-regulated genes in cells treated with ITE (ES: −0.41; P-value = 1 × 10−04) (Fig 5A and B). This is consistent with transcriptional activation by AHR that is inhibited in hypoxia and consistent with the overall changes in AHR chromatin binding seen by ChIP-seq. Conversely, HIF-bound genes—those closest to canonical HIF-only binding sites—were enriched among hypoxia-up-regulated genes (ES: 0.97; P-value = 1 × 10−04), although they were not significantly enriched among ITE-down-regulated genes under hypoxia (ES: −0.21; P-value = 0.25) (Fig 5C and D). Furthermore, AHR-specific target genes that were up-regulated by ITE and bound by AHR, but not by HIF, were down-regulated in hypoxia (median log2FC: −0.235; P-value = 3 × 10−06, compared with normoxia), whereas HIF-specific target genes that were up-regulated by hypoxia and bound by HIF, but not by AHR, were down-regulated by ITE (median log2FC: −0.076; P-value = 0.009, compared with DMSO) (Fig 5E). This bi-directional nature of the crosstalk between HIF and AHR in gene regulation is consistent with that observed for chromatin binding.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Heterogenous effects of crosstalk between HIF and AHR on gene regulation.

(A) Gene set enrichment analysis (GSEA) showing enrichment of genes adjacent to canonical AHR binding sites among genes induced in normoxia + ITE compared with normoxia + DMSO (ranked on the x-axis). (B) GSEA showing enrichment of the same genes among genes down-regulated in hypoxia + ITE compared with normoxia + ITE. (C) GSEA showing enrichment of genes adjacent to canonical HIF binding sites among genes up-regulated in hypoxia + DMSO compared with normoxia + DMSO. (D) GSEA showing enrichment of the same genes among genes down-regulated in hypoxia + ITE compared with hypoxia + DMSO. (E) Box-and-whisker plot showing down-regulation of AHR-bound only and ITE-induced genes in hypoxia + ITE compared with normoxia + ITE (blue; P-value = 3 × 10−06; paired t-test) and of HIF-bound only and hypoxia-induced genes in hypoxia + ITE compared with hypoxia + DMSO (red; P-value = 0.009; paired t-test). (F) GSEA showing enrichment of HIF and AHR co-bound genes among genes up-regulated in hypoxia + DMSO compared with normoxia + DMSO. (G) GSEA showing enrichment of HIF and AHR co-bound genes among genes up-regulated in normoxia + ITE compared with normoxia + DMSO. (H) Box-and-whisker plot showing the regulation by hypoxia (blue), ITE (red), or the combination of both stimuli (grey) of common HIF- and AHR-target genes.

In contrast, co-bound genes—those closest to both an AHR and HIF binding site—were significantly enriched among both hypoxia-up-regulated (ES: 0.97; P-value = 0.03) and ITE-up-regulated genes (ES: 0.96; P-value = 0.07), consistent with these genes being direct transcriptional targets of both transcription factors (Fig 5F and G). This analysis defined a small set of 14 genes (PAX2, MOAP1, WSB1, FAM65C, MFNG, CARD10, LIMCH1, NR3C1, PADI1, SLC2A1, ARNTL, INSIG1, PAG1, VLDLR) that are bound by both HIF and AHR and up-regulated by both hypoxia and ITE. For these genes, there was an additive effect on gene regulation when the two stimuli were combined (Fig 5H). This characterises a third point of interaction between the two transcriptional pathways, in which a small number of genes are directly regulated by both transcription factors.

Evidence of HIF and AHR antagonism in ccRCC tumours

To further evaluate the relevance of our findings to kidney cancer, we next examined the crosstalk between HIF and AHR in RNA-sequencing analyses of ccRCC tumours from The Cancer Genome Atlas (TCGA-KIRC) cohort. Because HIF-α and AHR mRNA levels do not represent activation of these pathways, we first used our ChIP-seq and RNA-seq analyses from RCC4+VHL cells to define an AHR gene signature comprising 67 AHR-bound and ITE-up-regulated genes and an HIF gene signature comprising 236 HIF-bound and hypoxia-induced genes (genes bound and up-regulated by both factors were excluded). After normalisation, individual gene expression levels were combined to give composite scores reflecting overall AHR or HIF activation. We then compared these composite scores for clear cell kidney tumours and paired surrounding kidney tissue. As expected, HIF composite scores were significantly elevated in ccRCC tumours compared with patient-matched normal kidney (P = 6 × 10−13), confirming constitutive HIF activation in these tumours as a result of VHL loss (Fig 6A). In contrast, AHR composite scores were significantly lower in tumours compared with patient-matched normal kidney (P = 1 × 10−5) (Fig 6B). Importantly, when we repeated the analysis using a composite score based on common HIF- and AHR-target genes, no significant difference was observed, indicating that the increase in HIF activity is balanced by the reduction in AHR activity at these common loci (Fig 6C). In ccRCC, VHL inactivation frequently results (at least in part) from copy number loss. In the TCGA-KIRC tumour samples, increasing VHL copy number loss was associated with a significant increase in the HIF composite score (P = 7 × 10−6) (Fig 6D), whereas an inverse association was seen for the AHR composite score (P = 5 × 10−11) (Fig 6E). Conversely, no significant association was observed for a composite score based on the common HIF- and AHR-target genes (Fig 6F). Finally, we compared the HIF and AHR composite scores for each individual tumour and found a significant negative correlation between the two composite scores, further supporting an antagonism between the pathways (Fig 6G). Taken together, we observe a negative correlation between overall HIF and AHR transcriptional activation in ccRCC tumours from the TCGA-KIRC cohort. This is consistent with the antagonism seen between the two transcription factors in ccRCC cell lines and indicates that endogenous AHR activity in ccRCC interacts with HIF activation.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Crosstalk between HIF and AHR in ccRCC tumours and in response to chemotherapeutic agents.

(A) Box-and-whisker plot showing the composite score for an HIF metagene (based on 236 genes that were bound only by HIF and up-regulated by hypoxia in DMSO-treated RCC4+VHL cells) in RNA-seq analysis of normal and tumour samples from the TCGA-KIRC cohort of clear cell renal cancer. (B, C) The same analysis for (B) an AHR composite score based on 67 genes that were bound only by AHR and up-regulated by ITE in normoxic RCC4+VHL cells and (C) a composite score based on common HIF- and AHR-bound and regulated genes. (D, E, F) The same analyses performed on tumour samples stratified according to copy number alteration at the VHL gene locus. (G) Scatter plot showing the correlation between AHR and HIF composite scores for the TCGA-KIRC tumour samples. P-values (A, B, C, D, E, F, G) by paired Wilcoxon rank sum test. (H) RT-qPCR analysis of NDRG1, CYP1B1, and ALDH1A3 mRNA levels in 786O cells exposed to DMSO (control), PT2385, ITE, or both, presented as mean ± SD of five biological replicates. *P ≤ 0.05, **P ≤ 0.01 (paired t-test).

Finally, the HIF pathway has long been considered a target for the treatment of ccRCC tumours. In this tumour type, HIF-2 has a role as a key ccRCC oncogene, whereas HIF-1 appears to act more as a tumour suppressor and is often repressed by the tumour. Thus, effort has focussed on the development of small molecule HIF-2 dimerization inhibitors and their use as therapeutic agents (34). The first HIF-2 dimerization inhibitor (belzutifan/PT2977) has recently been approved by the US Food and Drug Administration for the treatment of ccRCC in patients with the familial VHL syndrome (14, 35). Because blocking HIF-2 dimerization will liberate HIF-1β/ARNT, we next evaluated the effect of the first-generation HIF-2 inhibitor PT2385 on AHR activity. 786O parental cells (defective for both WT VHL and HIF-1α) were pre-treated with PT2385 for 18 h (to ensure complete inhibition of constitutive HIF-2 activity and downstream gene expression) followed by ITE for 18 h. PT2385 treatment completely abolished expression of the HIF-specific gene NDRG1, which was also significantly repressed by ITE treatments (Fig 6H). Conversely, PT2385 treatment significantly increased the ITE-dependent induction of the AHR-specific genes CYP1B1 and ALDH1A3. Taken together, this indicates that in addition to direct effects on HIF-2 signaling in ccRCC tumours, belzutifan may have indirect effects on AHR activation. Because AHR activation may promote tumourigenesis, this raises the possibility that there may be therapeutic synergy in targeting both pathways simultaneously.

Discussion

In this work, we demonstrate diverse and unexpected interactions between the HIF and AHR pathways that operate both generally and in a locus-specific manner to influence chromatin binding and gene expression. These interactions are both antagonistic and cooperative and include (1) competition between HIF-α and AHR subunits for binding to HIF-1β/ARNT, (2) synergistic binding, in cis, of the two transcription factors when bound close to each other on chromatin, and (3) partial overlap in the transcriptional targets of the two transcription factors. Our findings provide new insights into both the nature of the interplay between the HIF and AHR pathways and its implications for their therapeutic manipulation in cancer.

Mechanistic implications

Firstly, we have shown that AHR and HIF-α compete for binding to HIF-1β/ARNT, leading to global antagonism between the two pathways. Previous studies have shown that this can be reversed by overexpressing HIF-1β/ARNT, confirming that it is a limiting subunit (36). Also, consistent with previous reports (28), activation of HIF inhibited AHR-target gene expression and AHR-HIF-1β/ARNT dimerization in all settings, but there was no significant effect of AHR activation on HIF-target genes in three out of five cell lines studied. Only in two ccRCC cell lines, 786O + VHL and RCC4+VHL, did activation of AHR inhibit HIF-target gene expression and only in RCC4+VHL cells did the crosstalk between the two factors operate in both directions. Notably, these two cell lines had relatively high ratios of AHR to HIF-α subunits. Therefore, although the relative affinities of HIF-1α, HIF-2α, and AHR for HIF-1β/ARNT may influence the symmetry and strength of their competition (37), the variation in behaviour between cell lines and its broad correlation with the relative stoichiometry of the different subunits suggests that the latter also has a role in defining the directionality of this competition for HIF-1β/ARNT. In this respect, it is also important to consider other β-class bHLH-PAS paralogues of HIF-1β/ARNT: ARNT2, ARNTL, and ARNTL2 (1, 2, 3). Interaction of AHR and HIF-α with ARNTL and ARNTL2 is relatively weak and therefore unlikely to be of significance in vivo (3, 38, 39, 40). However, studies suggest that ARNT2 can substitute for HIF-1β/ARNT to dimerize with both AHR and HIF-α subunits and direct transcription (3, 38, 39, 40). Nevertheless, ARNT2 expression is highly tissue specific, being predominantly expressed in the brain and so is unlikely to contribute to the crosstalk between AHR and HIF-α outside of restricted cell types (40). Indeed, analysis of publicly available RNA-seq analyses indicates that ARNT mRNA is 7–35 times more abundant than ARNT2 mRNA in the cell lines studied.

The limited availability of HIF-1β/ARNT raises the possibility of even more complex interplay between HIF, AHR, and other bHLH-PAS family members that are also known to dimerize with HIF-1β/ARNT. Specifically, AHRR (aryl hydrocarbon receptor repressor), the SIM proteins (single-minded 1 and 2), and NPAS1 and NPAS3 (neuronal PAS domain protein 1 and 3) have all been reported to dimerize with HIF-1β/ARNT (41, 42, 43). Notably, a competitive relationship between AHR and AHRR has previous been described, whereby AHRR was found to sequester HIF-1β/ARNT and reduce AHR transcriptional activity (43, 44). As a consequence, AHR was redirected towards a non-canonical cytoplasmic E3 ligase function, targeting the estrogen receptor α for degradation (44). The availability of HIF-1β/ARNT was therefore suggested as the basis for a switch between AHR canonical and non-canonical functions, and notably, we have observed relocation of AHR from nuclear to cytoplasmic cell fractions in response to HIF-α stabilisation (Figs 2A and B and S3A and B). Whether reduced HIF-1β/ARNT availability under HIF-stabilising conditions is able to redirect AHR to other non-canonical cytoplasmic functions in RCC cells requires further investigation but suggests an additional level of interplay between HIF and the AHR pathway.

Secondly, analysis of AHR and HIF DNA-binding across the genome revealed several further levels of interaction. In RCC4+VHL cells, the competition between AHR and HIF-α subunits for HIF-1β/ARNT was reflected in binding to chromatin, with activation of either HIF or AHR leading to a reciprocal global reduction in binding of the other factor. This is likely a direct consequence of competition for HIF-1β/ARNT. However, where HIF and AHR bound close to each other, this mutual inhibition was reduced. This indicates that rather than competing for binding to shared chromatin motifs, HIF and AHR positively reinforce the binding of each other in cis. Indeed, when bound at overlapping sites, AHR and HIF binding do not coincide precisely, indicating that they bind to separate response elements. Furthermore, cooperation between transcription factors is commonly observed at enhancers and can result through a direct interaction between simultaneously co-bound factors or indirectly through effects on chromatin remodelling or nucleosome positioning (45, 46).

Thirdly, a small set of genes that are bound by both transcription factors are also induced by each factor independently and in combination, indicating that they are direct transcriptional targets of both AHR and HIF, which can act synergistically. Notably, this list of genes includes ARNTL, which is the dimerization partner for the regulator of circadian rhythms, CLOCK (38). This raises the possibility of still further crosstalk between AHR, HIF, and other bHLH-PAS transcription factors.

Relevance to kidney cancer

AHR mRNA levels are elevated in many types of cancer when compared with normal tissues (including ccRCC), and AHR has important roles in regulating cellular proliferation, invasion, migration, and immunity that help drive tumour initiation, promotion, progression, and metastasis (26, 47, 48, 49). However, AHR activation can have both pro- and anti-tumourigenic actions depending upon the cell line or tumour model used. In ccRCC, activation of AHR increases cancer cell invasion and is associated with poor patient survival (30). Conversely, pharmacological stimulation of AHR may also slow proliferation and migration of ccRCC cells in vitro and in xenograft models (50).

Similarly, HIF is activated in many cancer types, in which it is associated with intra-tumour hypoxia, poor patient prognosis, and resistance to therapy. However, like AHR, HIF can be both pro- and anti-tumourigenic depending upon the setting, and in ccRCC, HIF-2 promotes tumour growth, whereas HIF-1 is tumour suppressive. This has led to the development of selective HIF-2 inhibitors, which specifically block dimerization of HIF-2α with HIF-1β/ARNT to slow renal cancer progression (34). However, the new data reveal that in ccRCC tumours, competition for HIF-1β/ARNT results in a reciprocal relationship between HIF and AHR and constitutive activation of HIF is associated with suppression of the AHR pathway. Therefore, therapies that block HIF-2 dimerization can release more HIF-1β/ARNT to bind with AHR, thereby potentially promoting the oncogenic effects of the AHR pathway in kidney cancer. In summary, our data predict both that therapeutic inhibition of AHR in RCC has activating effects on components of the HIF pathway and that inhibition of HIF has activating effects on the AHR pathway. This raises important questions as to whether combination therapy targeting both pathways might show synergistic efficacy over and above targeting each system separately.

Materials and Methods

Cell culture

RCC4 renal cancer cells were a gift from CH Buys. RCC4+VHL and RCC4+VA (“vector alone”) clones were generated by stable transfection with pcDNA3-VHL or pcDNA3, respectively (9). 786O renal cancer cells, HepG2 liver cancer cells, and PC3 prostate cancer cells were purchased directly from ATCC. The 786O+VHL clone was a gift from WG Kaelin. HKC8 renal tubular cells were a gift from L.C. Racusen (51). HKC8 HIF-1α/KO, HIF-2α/KO, and 1α/2α dKO clones were generated through CRISPR-Cas9 disruption, as previously described (13). Unless used directly from a certified source, the identity of all cell lines was confirmed by STR genotyping. RNA-seq analysis of RCC4 and 786O cells also confirmed the presence of unique VHL gene coding mutations (chr3:10,183,725C > G and chr3:10,183,841 G > del, respectively) that are as previously described. All cell lines were grown in DMEM supplemented with 100 U/ml penicillin, 100 μg/ml streptomycin, 2 mM L-glutamine, and 10% foetal bovine serum (Sigma-Aldrich) and regularly tested for mycoplasma infection. RCC4+VHL, RCC4+VA and 786O+VHL cells were maintained in 0.5 mg/ml G418. Cells were maintained at 37°C under an atmosphere of 5% CO2 in air. Cells were allowed to grow in phenol red–free DMEM supplemented with 100 U/ml penicillin, 100 μg/ml streptomycin, 2 mM L-glutamine, and 10% dextran-coated charcoal-treated foetal calf serum for ∼48 h before treatments to minimise basal AHR activity. Treatments of sub-confluent cells were performed as indicated with 100 nM ITE, 10 nM FICZ (both from Tocris Bioscience), 50 μM FG4592 (Selleckchem), 1 μM PT2385 (MedChemExpress), or 0.1% DMSO for specified duration. Hypoxic incubations were conducted for specified duration at 0.5% ambient oxygen concentration (5% CO2; balance N2) within an InvivO2 400 workstation (Baker Ruskinn).

Whole cell lysate preparation

Cells were washed with PBS, and lysates were prepared in urea-SDS lysis buffer (6.7 M urea, 10% glycerol, 1% SDS, 10 mM Tris–HCl, pH 8.0) supplemented with freshly added cOmplete Protease Inhibitor Cocktail (Roche) and 1 mM DTT. Lysates were incubated 15 min at 4°C, then sonicated using a Bioruptor (Diagenode) at medium intensity for four cycles of 15 s on/30 s off, before being centrifuged for 5 min at 16,000g at 4°C. Protein concentrations were determined by BCA assay (Pierce).

Cytoplasmic and nuclear lysate preparation

Cytoplasmic and nuclear protein extraction was performed based on the published Method 2 from Klenova et al (52). Briefly, cells were washed with PBS and resuspended in Tween-20 lysis buffer (0.5% Tween-20; 25 mM Tris–HCl, pH 8.0; 2 mM EDTA, pH 8.0) containing 50 mM NaCl and supplemented with freshly added cOmplete Protease Inhibitor Cocktail (Roche). Lysates were incubated for 15 min at 4°C, then centrifuged for 1 min at 6,000g at 4°C. Supernatants were recovered as cytoplasmic fractions. Nuclear pellets were washed and resuspended in Tween-20 lysis buffer containing 500 mM NaCl. Lysates were incubated 15 min at 4°C before, then supplemented with Tween-20 lysis buffer without NaCl (for a final concentration of 250 mM) and centrifuged for 15 min at 10,000g at 4°C. Supernatants were kept as nuclear fractions. Protein concentrations were determined by BCA assay (Pierce).

Immunoprecipitation

Cytoplasmic or nuclear cell lysates (between 150–300 μg per experiment) were diluted to 1 ml in Tween-20 buffer (0.5% Tween-20; 25 mM Tris–HCl, pH 8.0; 2 mM EDTA, pH 8.0) containing 150 mM NaCl and supplemented with freshly added cOmplete Protease Inhibitor Cocktail (Roche). Samples were pre-cleared by incubating with protein A-agarose beads (Millipore) for 1 h at 4°C on a rotating mixer. Samples were then centrifuged for 3 min at 400g. Supernatants were incubated with non-immunized rabbit serum (X0902; Dako) or HIF-1β/ARNT antibody (rabbit polyclonal; cat no. NB100-110; Novus Biologicals) (between 2–4 μl) overnight at 4°C on a rotating mixer, before incubation with protein A-agarose beads (between 25–50 μl) for 1 h at 4°C. Beads were washed three times in Tween-20 buffer for 5 min at 4°C. Protein complexes were eluted from beads by adding Laemmli sample buffer with 50 mM DTT and incubating at 95°C for 5 min. Proteins of interest were detected by immunoblotting IP and input samples (HIF-1β/ARNT was detected using mouse monoclonal antibody NB100-124; see below).

Immunoblotting

Samples were prepared for SDS–PAGE by adding Laemmli sample buffer with 50 mM DTT and incubating at 95°C for 5 min. Proteins were then resolved on acrylamide gels using Mini-PROTEAN Tetra cells (Bio-Rad) and transferred to polyvinylidene difluoride membrane (Immobilon-P; Millipore) using mini-trans-blot cells (Bio-Rad). Membranes were blocked for 1 h in 5% milk solution prepared in PBS-T (phosphate-buffered saline containing 0.1% Tween-20), before being incubated in primary antibodies overnight at 4°C (in 5% milk/PBS-T). Primary antibodies used were anti-AHR (A-3, cat no. sc-133088; Santa Cruz Biotechnology), anti-HIF-1α (cat no. 610959; BD Biosciences), anti-HIF-2α (mouse monoclonal, 190 b), anti-HIF-1β/ARNT (D28F3, cat no. 5537; Cell Signaling), anti-HIF-1β/ARNT (H1β234, cat no. NB100-124; Novus Biologicals), anti-CYP1B1 (EPR14972, cat no. ab185954; Abcam), anti-NDRG1 (cat no. 5196; Cell Signaling), anti-HDAC2 (cat no. PA1-861; Invitrogen), anti-α-tubulin (cat no. ab4074; Abcam), anti-PolII (CTD4H8, cat no. 05-623; Millipore), and anti-β-actin-HRP (cat no. ab49900; Abcam). HRP-conjugated secondary antibodies (rabbit, cat no. P0048; mouse, cat no. P0047; Dako) were incubated for 1 h at room temperature. Membranes were developed using either SuperSignal West Dura Extended Duration Substrate or SuperSignal West Femto Maximum Sensitivity Substrate (Thermo Fisher Scientific) and imaged using a ChemiDoc XRS+ imaging system (Bio-Rad). After immunoblot analysis, membranes were stained with Coomassie brilliant blue to visualize separated proteins. Densitometric analysis was performed using Image Lab Software (Bio-Rad) and expressed as the mean ± SD of at least three independent experiments. Statistical significance was tested by paired t-test using GraphPad Prism software (v9).

Chromatin immunoprecipitation (ChIP) coupled to qPCR

ChIP experiments were performed as previously described (53) using antibodies directed against AHR (rabbit monoclonal; D5S6H, cat no. 83200; Cell Signaling), HIF-1α (rabbit polyclonal, PM14), HIF-2α (rabbit polyclonal, PM9), or HIF-1β/ARNT (rabbit polyclonal; cat no. NB100-110; Novus Biologicals) (54, 55). Non-immunized rabbit serum (X0902; Dako) was used as a negative control. One 15-cm plate of cells was used for each ChIP. qPCR analysis was performed using Fast SYBR Green Master Mix on a StepOne thermocycler (both from Applied Biosystems) with the following oligonucleotides: HIF binding site at the NDRG1 locus (forward: TCCCTCCCAATCTCTCTCTTCTT; reverse: CACCATCAGCACAGCAAACTAC) and AHR binding site at the CYP1B1 locus (forward: CGCTCCGGCGAACTTTAT; reverse: GATATGACTGGAGCCGACTTTC).

RNA extraction and RT–qPCR

Cells were lysed in TRIzol and total RNA was prepared using Direct-zol RNA MiniPrep Kit (Zymo Research), according to manufacturer’s instructions. Equal amounts of RNA were used for cDNA synthesis using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Expression analyses were performed using Fast SYBR Green Master Mix on a StepOne thermocycler (both from Applied Biosystems, now Thermo Fisher Scientific) and calculated using the Pfaffl method. The list of primers used for analysis is found in Table S1, with HPRT1 used as reference gene for normalisation. Data are expressed as the mean ± SD of at least three independent experiments. Statistical significance was tested by paired t-test using GraphPad Prism software (v9).

Table S1 List of primers for RT–qPCR analysis.

ChIP-sequencing

ChIP experiments were performed as above in duplicate, in accordance with ENCODE consortium guidelines (The_Encode_Consortium, 2017). Automated library preparation was performed at the Oxford Genomics Centre (University of Oxford) using the Apollo prep system (PrepX Comp ILMN 32i DNA Lib 96. cat no. 640102; Takara). Sequencing was performed on the NovaSeq 6000 platform, using 150 bp paired-end configuration, according to Illumina protocols (Illumina).

RNA-sequencing

RNA extraction was performed as above from six independent replicates, in accordance with ENCODE consortium guidelines (The_ENCODE_Consortium, 2016). PolyA+ RNA library preparation was conducted at GENEWIZ/Azenta Life Sciences using the NEBNext Ultra II Directional RNA Library Prep Kit for Illumina. Sequencing was performed on the NovaSeq platform, using 150 bp paired-end configuration, according to Illumina protocols (Illumina).

Bioinformatic analysis of ChIP-seq data

Preliminary analysis

Illumina adaptor sequences were trimmed using TrimGalore (0.3.3). Reads were aligned to GRCh37 using BWA (0.7.5a-r405). Low-quality mapping was removed (MapQ < 15) using SAMtools (0.1.19) (56). Reads mapping to Duke Encode black list regions https://hgwdev.gi.ucsc.edu/cgi-bin/hgFileUi?db=hg19&g=wgEncodeMapability were excluded using BEDTools (2.17.0) (57). Duplicate reads were excluded using Picard tools (2.0.1) (http://broadinstitute.github.io/picard/). Read densities were normalised and expressed as reads per kilobase per million reads (RPKM) (58). One million random non-overlapping regions selected from ENCODE DNase Cluster II peaks (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeRegDnaseClustered/) were used as a control.

Peak calling

ChIP-seq peaks were identified using T-PIC (Tree shape Peak Identification for ChIP-Seq) (59) and MACS (model-based analysis of ChIP-seq) (60) in control mode. Peaks detected by both peak callers were filtered quantitatively using the total count under the peak to include only peaks that were above the 99.99th percentile of random background regions selected from the ENCODE DNase II cluster (P-value < 0.0001). Only peaks from each independent replicate that overlapped by at least 1 base pair (BEDTools 2.17.0 (57)) were considered.

Quantitation of ChIP-seq signal

Manipulation of .bam files was performed using SAMtools (0.1.19) (56). Briefly, SAMtools merge was used to merge sorted alignment .bam files from each replicate into a single .bam output file, which was then indexed using SAMtools index. SAMtools bedcov was then used to determine the read depth per bed region, which was normalised to the total reads in each dataset as determined using SAMtools flagstat.

Normalisation of ChIP-seq signal

ChIP-seq normalisation was performed using MAnorm2 (v1.0.0) (61). Briefly, genomic intervals were created using profile_bins from MAnorm2 with a bin size of 1,000 and shift size of 150 to count the reads that fall within each genomic interval. MA normalisation was then performed using “normalize” function from MAnorm2 with default options.

Differential binding analysis

Differential binding analysis for each comparison was performed using ChIPComp (v1.16.0) (62), and the Wald test P-value was corrected for multiple testing using the Benjamini–Hochberg procedure.

Peak annotation

The distance from each binding site to the nearest TSS in the Ensembl hg19 database was determined using PeakAnnotator (v1.4, https://www.bioinformatics.org/peakanalyzer/wiki/Main/Download) and plotted as a histogram in R (63).

Data visualization

Locus-specific data were visualized using Integrative Genomics Viewer (IGV; 2.14.13) (64). Chromosomal coordinates and gene annotation are from the RefSeq hg19 (GRCh37) build.

Bioinformatic analysis of RNA-seq data

Preliminary analysis

Adaptor sequences were trimmed using TrimGalore (0.3.3). Reads were then aligned to GRCh37 using HISAT2 (2.0.5) (65), and non-uniquely mapped fragments were excluded using Picard tools (2.0.1) (http://broadinstitute.github.io/picard/). Total read counts for each UCSC-defined gene were extracted using HTSeq (0.5.4p3) (66) with “intersection-strict” mode, and significantly regulated genes were identified using DESeq2 (67).

GSEA

GSEA used 10,000 permutations, weighted ES, and pre-ranking of genes (68). Both differential expression significance according to DESeq2 and fold difference between the two conditions were used to rank genes according to the equation:πi=φi−log10Pυi,where ϕi is the log2 fold change and Ρυi is the P-value for gene i (69).

TCGA RNA-seq data analysis

To compare expression in ccRCC with normal kidney, gene expression levels were extracted from TCGA-KIRC level 3 RNA-seq data for 72 paired renal tumours and surrounding normal renal tissue and compared using a paired Wilcoxon rank sum test.

Bioinformatic statistical analyses

All statistical analyses were performed using R, version 4.0.5.

Data Availability

High-throughput sequencing data are available at Gene Expression Omnibus, accession number GSE206029.

Acknowledgements

This study was supported by the National Institute for Health Research (DR Mole; NIHR-RP-2016-06-004), Cancer Research UK (DR Mole; A416016), the Deanship of Scientific Research, King Abdulaziz University, Ministry of High Education for Saudi Arabia (H Choudhry, DR Mole, and PJ Ratcliffe), the Ludwig Institute for Cancer Research (PJ Ratcliffe), the Wellcome Trust (PJ Ratcliffe; 106241/Z/14/Z), and the Francis Crick Institute (PJ Ratcliffe), which receives its core funding from Cancer Research UK (FC001501), the UK Medical Research Council (FC001501), and the Wellcome Trust (FC001501). The computational aspects of this research were supported by the Wellcome Trust Core Award Grant Number 203141/Z/16/Z and the NIHR Oxford BRC. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health.

Author Contributions

  • VN Lafleur: conceptualization, data curation, formal analysis, validation, investigation, visualization, methodology, and writing—original draft, review, and editing.

  • S Halim: data curation, software, formal analysis, validation, visualization, methodology, and writing—original draft, review, and editing.

  • H Choudhry: conceptualization, supervision, funding acquisition, and writing—review and editing.

  • PJ Ratcliffe: conceptualization, formal analysis, supervision, funding acquisition, investigation, and writing—review and editing.

  • DR Mole: conceptualization, resources, data curation, formal analysis, supervision, funding acquisition, validation, investigation, visualization, methodology, project administration, and writing—original draft, review, and editing.

Conflict of Interest Statement

PJ Ratcliffe is a scientific co-founder of ReOx Ltd., a company, which is developing inhibitors of the HIF hydroxylase enzymes. The other authors declare no competing interest.

  • Received October 6, 2022.
  • Revision received January 11, 2023.
  • Accepted January 12, 2023.
  • © 2023 Lafleur et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Bersten DC,
    2. Sullivan AE,
    3. Peet DJ,
    4. Whitelaw ML
    (2013) bHLH-PAS proteins in cancer. Nat Rev Cancer 13: 827–841. doi:10.1038/nrc3621
    OpenUrlCrossRefPubMed
  2. ↵
    1. Vazquez-Rivera E,
    2. Rojas BL,
    3. Parrott JC,
    4. Shen AL,
    5. Xing Y,
    6. Carney PR,
    7. Bradfield CA
    (2022) The aryl hydrocarbon receptor as a model PAS sensor. Toxicol Rep 9: 1–11. doi:10.1016/j.toxrep.2021.11.017
    OpenUrlCrossRef
  3. ↵
    1. Gu YZ,
    2. Hogenesch JB,
    3. Bradfield CA
    (2000) The PAS superfamily: Sensors of environmental and developmental signals. Annu Rev Pharmacol Toxicol 40: 519–561. doi:10.1146/annurev.pharmtox.40.1.519
    OpenUrlCrossRefPubMed
  4. ↵
    1. Keith B,
    2. Johnson RS,
    3. Simon MC
    (2011) HIF1α and HIF2α: Sibling rivalry in hypoxic tumour growth and progression. Nat Rev Cancer 12: 9–22. doi:10.1038/nrc3183
    OpenUrlCrossRefPubMed
  5. ↵
    1. Xiong Q,
    2. Liu B,
    3. Ding M,
    4. Zhou J,
    5. Yang C,
    6. Chen Y
    (2020) Hypoxia and cancer related pathology. Cancer Lett 486: 1–7. doi:10.1016/j.canlet.2020.05.002
    OpenUrlCrossRef
  6. ↵
    1. Semenza GL
    (2012) Hypoxia-inducible factors: Mediators of cancer progression and targets for cancer therapy. Trends Pharmacological Sciences 33: 207–214. doi:10.1016/j.tips.2012.01.005
    OpenUrlCrossRefPubMed
  7. ↵
    1. Ratcliffe PJ
    (2013) Oxygen sensing and hypoxia signalling pathways in animals: The implications of physiology for cancer. J Physiol 591: 2027–2042. doi:10.1113/jphysiol.2013.251470
    OpenUrlCrossRefPubMed
  8. ↵
    1. Lombardi O,
    2. Mole DR
    (2021) HIF pathways in clear cell renal cancer Innovative Bioanalysis [Working Title].
  9. ↵
    1. Maxwell PH,
    2. Wiesener MS,
    3. Chang G-W,
    4. Clifford SC,
    5. Vaux EC,
    6. Cockman ME,
    7. Wykoff CC,
    8. Pugh CW,
    9. Maher ER,
    10. Ratcliffe PJ
    (1999) The tumour suppressor protein VHL targets hypoxia-inducible factors for oxygen-dependent proteolysis. Nature 399: 271–275. doi:10.1038/20459
    OpenUrlCrossRefPubMed
  10. ↵
    1. Wu D,
    2. Rastinejad F
    (2017) Structural characterization of mammalian bHLH-PAS transcription factors. Curr Opin Struct Biol 43: 1–9. doi:10.1016/j.sbi.2016.09.011
    OpenUrlCrossRef
  11. ↵
    1. Salama R,
    2. Masson N,
    3. Simpson P,
    4. Sciesielski LK,
    5. Sun M,
    6. Tian Y-M,
    7. Ratcliffe PJ,
    8. Mole DR
    (2015) Heterogeneous effects of direct hypoxia pathway activation in kidney cancer. PLoS One 10: e0134645. doi:10.1371/journal.pone.0134645
    OpenUrlCrossRefPubMed
  12. ↵
    1. Raval RR,
    2. Lau KW,
    3. Tran MGB,
    4. Sowter HM,
    5. Mandriota SJ,
    6. Li JL,
    7. Pugh CW,
    8. Maxwell PH,
    9. Harris AL,
    10. Ratcliffe PJ
    (2005) Contrasting properties of hypoxia-inducible factor 1 (HIF-1) and HIF-2 in von Hippel-Lindau-associated renal cell carcinoma. Mol Cell Biol 25: 5675–5686. doi:10.1128/mcb.25.13.5675-5686.2005
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Smythies JA,
    2. Sun M,
    3. Masson N,
    4. Salama R,
    5. Simpson PD,
    6. Murray E,
    7. Neumann V,
    8. Cockman ME,
    9. Choudhry H,
    10. Ratcliffe PJ, et al.
    (2019) Inherent DNA-binding specificities of the HIF-1α and HIF-2α transcription factors in chromatin. EMBO Rep 20: e46401. doi:10.15252/embr.201846401
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Deeks ED
    (2021) Belzutifan: First approval. Drugs 81: 1921–1927. doi:10.1007/s40265-021-01606-x
    OpenUrlCrossRef
  15. ↵
    1. Schödel J,
    2. Grampp S,
    3. Maher ER,
    4. Moch H,
    5. Ratcliffe PJ,
    6. Russo P,
    7. Mole DR
    (2016) Hypoxia, hypoxia-inducible transcription factors, and renal cancer. Eur Urol 69: 646–657. doi:10.1016/j.eururo.2015.08.007
    OpenUrlCrossRefPubMed
  16. ↵
    1. Shen C,
    2. Beroukhim R,
    3. Schumacher SE,
    4. Zhou J,
    5. Chang M,
    6. Signoretti S,
    7. Kaelin WG Jr.
    (2011) Genetic and functional studies implicate HIF1α as a 14q kidney cancer suppressor gene. Cancer Discov 1: 222–235. doi:10.1158/2159-8290.cd-11-0098
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Xu J,
    2. Wang B,
    3. Xu Y,
    4. Sun L,
    5. Tian W,
    6. Shukla D,
    7. Barod R,
    8. Grillari J,
    9. Grillari-Voglauer R,
    10. Maxwell PH, et al.
    (2012) Epigenetic regulation of HIF-1α in renal cancer cells involves HIF-1α/2α binding to a reverse hypoxia-response element. Oncogene 31: 1065–1072. doi:10.1038/onc.2011.305
    OpenUrlCrossRefPubMed
  18. ↵
    1. Schmid V,
    2. Lafleur VN,
    3. Lombardi O,
    4. Li R,
    5. Salama R,
    6. Colli L,
    7. Choudhry H,
    8. Chanock S,
    9. Ratcliffe PJ,
    10. Mole DR
    (2019) Co-incidence of RCC-susceptibility polymorphisms with HIF cis-acting sequences supports a pathway tuning model of cancer. Sci Rep 9: 18768. doi:10.1038/s41598-019-55098-7
    OpenUrlCrossRef
  19. ↵
    1. Hoefflin R,
    2. Harlander S,
    3. Schäfer S,
    4. Metzger P,
    5. Kuo F,
    6. Schönenberger D,
    7. Adlesic M,
    8. Peighambari A,
    9. Seidel P,
    10. Chen CY, et al.
    (2020) HIF-1α and HIF-2α differently regulate tumour development and inflammation of clear cell renal cell carcinoma in mice. Nat Commun 11: 4111. doi:10.1038/s41467-020-17873-3
    OpenUrlCrossRefPubMed
  20. ↵
    1. Larigot L,
    2. Juricek L,
    3. Dairou J,
    4. Coumoul X
    (2018) AhR signaling pathways and regulatory functions. Biochimie open 7: 1–9. doi:10.1016/j.biopen.2018.05.001
    OpenUrlCrossRefPubMed
  21. ↵
    1. Nebert DW
    (1989) The Ah locus: Genetic differences in toxicity, cancer, mutation, and birth defects. Crit Rev Toxicol 20: 153–174. doi:10.3109/10408448909017908
    OpenUrlCrossRefPubMed
  22. ↵
    1. Boffetta P,
    2. Jourenkova N,
    3. Gustavsson P
    (1997) Cancer risk from occupational and environmental exposure to polycyclic aromatic hydrocarbons. Cancer Causes Control 8: 444–472. doi:10.1023/a:1018465507029
    OpenUrlCrossRefPubMed
  23. ↵
    1. Deutsch-Wenzel RP,
    2. Brune H,
    3. Grimmer G,
    4. Dettbarn G,
    5. Misfeld J
    (1983) Experimental studies in rat lungs on the carcinogenicity and dose-response relationships of eight frequently occurring environmental polycyclic aromatic hydrocarbons. J Natl Cancer Inst 71: 539–544. doi:10.1093/jnci/71.3.539
    OpenUrlCrossRefPubMed
  24. ↵
    1. Leclerc D,
    2. Staats Pires AC,
    3. Guillemin GJ,
    4. Gilot D
    (2021) Detrimental activation of AhR pathway in cancer: An overview of therapeutic strategies. Curr Opin Immunol 70: 15–26. doi:10.1016/j.coi.2020.12.003
    OpenUrlCrossRef
  25. ↵
    1. Beischlag TV,
    2. Morales JL,
    3. Hollingshead BD,
    4. Perdew GH
    (2008) The aryl hydrocarbon receptor complex and the control of gene expression. Crit Rev Eukaryot Gene Expr 18: 207–250. doi:10.1615/critreveukargeneexpr.v18.i3.20
    OpenUrlCrossRefPubMed
  26. ↵
    1. Paris A,
    2. Tardif N,
    3. Galibert MD,
    4. Corre S
    (2021) AhR and Cancer: From gene profiling to targeted therapy. Int J Mol Sci 22: 752. doi:10.3390/ijms22020752
    OpenUrlCrossRef
  27. ↵
    1. Zhao H,
    2. Chen L,
    3. Yang T,
    4. Feng YL,
    5. Vaziri ND,
    6. Liu BL,
    7. Liu QQ,
    8. Guo Y,
    9. Zhao YY
    (2019) Aryl hydrocarbon receptor activation mediates kidney disease and renal cell carcinoma. J Transl Med 17: 302. doi:10.1186/s12967-019-2054-5
    OpenUrlCrossRefPubMed
  28. ↵
    1. Vorrink SU,
    2. Domann FE
    (2014) Regulatory crosstalk and interference between the xenobiotic and hypoxia sensing pathways at the AhR-ARNT-HIF1α signaling node. Chem Biol Interact 218: 82–88. doi:10.1016/j.cbi.2014.05.001
    OpenUrlCrossRefPubMed
  29. ↵
    1. Lee K,
    2. Burgoon LD,
    3. Lamb L,
    4. Dere E,
    5. Zacharewski TR,
    6. Hogenesch JB,
    7. LaPres JJ
    (2006) Identification and characterization of genes susceptible to transcriptional cross-talk between the hypoxia and dioxin signaling cascades. Chem Res Toxicol 19: 1284–1293. doi:10.1021/tx060068d
    OpenUrlCrossRefPubMed
  30. ↵
    1. Ishida M,
    2. Mikami S,
    3. Shinojima T,
    4. Kosaka T,
    5. Mizuno R,
    6. Kikuchi E,
    7. Miyajima A,
    8. Okada Y,
    9. Oya M
    (2015) Activation of aryl hydrocarbon receptor promotes invasion of clear cell renal cell carcinoma and is associated with poor prognosis and cigarette smoke. Int J Cancer 137: 299–310. doi:10.1002/ijc.29398
    OpenUrlCrossRef
  31. ↵
    1. Pollenz RS
    (2002) The mechanism of AH receptor protein down-regulation (degradation) and its impact on AH receptor-mediated gene regulation. Chemico-Biological Interactions 141: 41–61. doi:10.1016/s0009-2797(02)00065-0
    OpenUrlCrossRefPubMed
  32. ↵
    1. Bailey TL
    (2011) DREME: motif discovery in transcription factor ChIP-seq data. Bioinformatics 27: 1653–1659. doi:10.1093/bioinformatics/btr261
    OpenUrlCrossRefPubMed
  33. ↵
    1. Seok SH,
    2. Lee W,
    3. Jiang L,
    4. Molugu K,
    5. Zheng A,
    6. Li Y,
    7. Park S,
    8. Bradfield CA,
    9. Xing Y
    (2017) Structural hierarchy controlling dimerization and target DNA recognition in the AHR transcriptional complex. Proc Natl Acad Sci U S A 114: 5431–5436. doi:10.1073/pnas.1617035114
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. W Choi WS,
    2. Boland J,
    3. Lin J
    (2021) Hypoxia-inducible factor-2α as a novel target in renal cell carcinoma. J Kidney Cancer VHL 8: 1–7. doi:10.15586/jkcvhl.v8i2.170
    OpenUrlCrossRef
  35. ↵
    1. Jonasch E,
    2. Donskov F,
    3. Iliopoulos O,
    4. Rathmell WK,
    5. Narayan VK,
    6. Maughan BL,
    7. Oudard S,
    8. Else T,
    9. Maranchie JK,
    10. Welsh SJ, et al.
    (2021) Belzutifan for renal cell carcinoma in von Hippel–Lindau disease. N Engl J Med 385: 2036–2046. doi:10.1056/NEJMoa2103425
    OpenUrlCrossRefPubMed
  36. ↵
    1. Vorrink SU,
    2. Severson PL,
    3. Kulak MV,
    4. Futscher BW,
    5. Domann FE
    (2014) Hypoxia perturbs aryl hydrocarbon receptor signaling and CYP1A1 expression induced by PCB 126 in human skin and liver-derived cell lines. Toxicol Appl Pharmacol 274: 408–416. doi:10.1016/j.taap.2013.12.002
    OpenUrlCrossRefPubMed
  37. ↵
    1. Gradin K,
    2. McGuire J,
    3. Wenger RH,
    4. Kvietikova I,
    5. fhitelaw ML,
    6. Toftgård R,
    7. Tora L,
    8. Gassmann M,
    9. Poellinger L
    (1996) Functional interference between hypoxia and dioxin signal transduction pathways: Competition for recruitment of the arnt transcription factor. Mol Cell Biol 16: 5221–5231. doi:10.1128/mcb.16.10.5221
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Hogenesch JB,
    2. Gu YZ,
    3. Jain S,
    4. Bradfield CA
    (1998) The basic-helix-loop-helix-PAS orphan MOP3 forms transcriptionally active complexes with circadian and hypoxia factors. Proc Natl Acad Sci U S A 95: 5474–5479. doi:10.1073/pnas.95.10.5474
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Hogenesch JB,
    2. Chan WK,
    3. Jackiw VH,
    4. Brown RC,
    5. Gu YZ,
    6. Pray-Grant M,
    7. Perdew GH,
    8. Bradfield CA
    (1997) Characterization of a subset of the basic-helix-loop-helix-PAS superfamily that interacts with components of the dioxin signaling pathway. J Biol Chem 272: 8581–8593. doi:10.1074/jbc.272.13.8581
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Hirose K,
    2. Morita M,
    3. Ema M,
    4. Mimura J,
    5. Hamada H,
    6. Fujii H,
    7. Saijo Y,
    8. Gotoh O,
    9. Sogawa K,
    10. Fujii-Kuriyama Y
    (1996) cDNA cloning and tissue-specific expression of a novel basic helix-loop-helix/PAS factor (Arnt2) with close sequence similarity to the aryl hydrocarbon receptor nuclear translocator (Arnt). Mol Cell Biol 16: 1706–1713. doi:10.1128/mcb.16.4.1706
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Probst MR,
    2. Fan CM,
    3. Tessier-Lavigne M,
    4. Hankinson O
    (1997) Two murine homologs of the Drosophila single-minded protein that interact with the mouse aryl hydrocarbon receptor nuclear translocator protein. J Biol Chem 272: 4451–4457. doi:10.1074/jbc.272.7.4451
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Teh CH,
    2. Lam KK,
    3. Loh CC,
    4. Loo JM,
    5. Yan T,
    6. Lim TM
    (2006) Neuronal PAS domain protein 1 is a transcriptional repressor and requires arylhydrocarbon nuclear translocator for its nuclear localization. J Biol Chem 281: 34617–34629. doi:10.1074/jbc.m604409200
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Mimura J,
    2. Ema M,
    3. Sogawa K,
    4. Fujii-Kuriyama Y
    (1999) Identification of a novel mechanism of regulation of Ah (dioxin) receptor function. Genes Dev 13: 20–25. doi:10.1101/gad.13.1.20
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Luecke-Johansson S,
    2. Gralla M,
    3. Rundqvist H,
    4. Ho JC,
    5. Johnson RS,
    6. Gradin K,
    7. Poellinger L
    (2017) A molecular mechanism to switch the aryl hydrocarbon receptor from a transcription factor to an E3 ubiquitin ligase. Mol Cell Biol 37: e00630-16. doi:10.1128/MCB.00630-16
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Long HK,
    2. Prescott SL,
    3. Wysocka J
    (2016) Ever-changing landscapes: Transcriptional enhancers in development and evolution. Cell 167: 1170–1187. doi:10.1016/j.cell.2016.09.018
    OpenUrlCrossRefPubMed
  46. ↵
    1. Spitz F,
    2. Furlong EEM
    (2012) Transcription factors: From enhancer binding to developmental control. Nat Rev Genet 13: 613–626. doi:10.1038/nrg3207
    OpenUrlCrossRefPubMed
  47. ↵
    1. Murray IA,
    2. Patterson AD,
    3. Perdew GH
    (2014) Aryl hydrocarbon receptor ligands in cancer: Friend and foe. Nat Rev Cancer 14: 801–814. doi:10.1038/nrc3846
    OpenUrlCrossRefPubMed
  48. ↵
    1. Button EL,
    2. Bersten DC,
    3. Whitelaw ML
    (2017) HIF has Biff - crosstalk between HIF1a and the family of bHLH/PAS proteins. Exp Cell Res 356: 141–145. doi:10.1016/j.yexcr.2017.03.055
    OpenUrlCrossRef
  49. ↵
    1. Mascanfroni ID,
    2. Takenaka MC,
    3. Yeste A,
    4. Patel B,
    5. Wu Y,
    6. Kenison JE,
    7. Siddiqui S,
    8. Basso AS,
    9. Otterbein LE,
    10. Pardoll DM, et al.
    (2015) Metabolic control of type 1 regulatory T cell differentiation by AHR and HIF1-α. Nat Med 21: 638–646. doi:10.1038/nm.3868
    OpenUrlCrossRefPubMed
  50. ↵
    1. Itkin B,
    2. Breen A,
    3. Turyanska L,
    4. Sandes EO,
    5. Bradshaw TD,
    6. Loaiza-Perez AI
    (2020) New treatments in renal cancer: The AhR ligands. Int J Mol Sci 21: 3551. doi:10.3390/ijms21103551
    OpenUrlCrossRef
  51. ↵
    1. Racusen LC,
    2. Monteil C,
    3. Sgrignoli A,
    4. Lucskay M,
    5. Marouillat S,
    6. Rhim JG,
    7. Morin JP
    (1997) Cell lines with extended in vitro growth potential from human renal proximal tubule: Characterization, response to inducers, and comparison with established cell lines. J Lab Clin Med 129: 318–329. doi:10.1016/s0022-2143(97)90180-3
    OpenUrlCrossRefPubMed
  52. ↵
    1. Klenova E,
    2. Chernukhin I,
    3. Inoue T,
    4. Shamsuddin S,
    5. Norton J
    (2002) Immunoprecipitation techniques for the analysis of transcription factor complexes. Methods 26: 254–259. doi:10.1016/s1046-2023(02)00029-4
    OpenUrlCrossRefPubMed
  53. ↵
    1. Mole DR,
    2. Blancher C,
    3. Copley RR,
    4. Pollard PJ,
    5. Gleadle JM,
    6. Ragoussis J,
    7. Ratcliffe PJ
    (2009) Genome-wide association of hypoxia-inducible factor (HIF)-1alpha and HIF-2alpha DNA binding with expression profiling of hypoxia-inducible transcripts. J Biol Chem 284: 16767–16775. doi:10.1074/jbc.M901790200
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Lau KW,
    2. Tian Y-M,
    3. Raval RR,
    4. Ratcliffe PJ,
    5. Pugh CW
    (2007) Target gene selectivity of hypoxia-inducible factor-alpha in renal cancer cells is conveyed by post-DNA-binding mechanisms. Br J Cancer 96: 1284–1292. doi:10.1038/sj.bjc.6603675
    OpenUrlCrossRefPubMed
  55. ↵
    1. Wiesener MS,
    2. Jürgensen JS,
    3. Rosenberger C,
    4. Scholze CK,
    5. Hörstrup JH,
    6. Warnecke C,
    7. Mandriota S,
    8. Bechmann I,
    9. Frei UA,
    10. Pugh CW, et al.
    (2003) Widespread hypoxia-inducible expression of HIF-2alpha in distinct cell populations of different organs. FASEB J 17: 271–273. doi:10.1096/fj.02-0445fje
    OpenUrlCrossRefPubMed
  56. ↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R
    (2009) The sequence alignment/map format and SAMtools. Bioinformatics 25: 2078–2079. doi:10.1093/bioinformatics/btp352
    OpenUrlCrossRefPubMed
  57. ↵
    1. Dale RK,
    2. Pedersen BS,
    3. Quinlan AR
    (2011) Pybedtools: A flexible Python library for manipulating genomic datasets and annotations. Bioinformatics 27: 3423–3424. doi:10.1093/bioinformatics/btr539
    OpenUrlCrossRefPubMed
  58. ↵
    1. Mortazavi A,
    2. Williams BA,
    3. McCue K,
    4. Schaeffer L,
    5. Wold B
    (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5: 621–628. doi:10.1038/nmeth.1226
    OpenUrlCrossRefPubMed
  59. ↵
    1. Hower V,
    2. Evans SN,
    3. Pachter L
    (2011) Shape-based peak identification for ChIP-Seq. BMC Bioinformatics 12: 15. doi:10.1186/1471-2105-12-15
    OpenUrlCrossRefPubMed
  60. ↵
    1. Zhang Y,
    2. Liu T,
    3. Meyer CA,
    4. Eeckhoute J,
    5. Johnson DS,
    6. Bernstein BE,
    7. Nusbaum C,
    8. Myers RM,
    9. Brown M,
    10. Li W, et al.
    (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9: R137. doi:10.1186/gb-2008-9-9-r137
    OpenUrlCrossRefPubMed
  61. ↵
    1. Tu S,
    2. Li M,
    3. Chen H,
    4. Tan F,
    5. Xu J,
    6. Waxman DJ,
    7. Zhang Y,
    8. Shao Z
    (2021) MAnorm2 for quantitatively comparing groups of ChIP-seq samples. Genome Res 31: 131–145. doi:10.1101/gr.262675.120
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Chen L,
    2. Wang C,
    3. Qin ZS,
    4. Wu H
    (2015) A novel statistical method for quantitative comparison of multiple ChIP-seq datasets. Bioinformatics 31: 1889–1896. doi:10.1093/bioinformatics/btv094
    OpenUrlCrossRefPubMed
  63. ↵
    1. Salmon-Divon M,
    2. Dvinge H,
    3. Tammoja K,
    4. Bertone P
    (2010) PeakAnalyzer: Genome-wide annotation of chromatin binding and modification loci. BMC Bioinformatics 11: 415. doi:10.1186/1471-2105-11-415
    OpenUrlCrossRefPubMed
  64. ↵
    1. Robinson JT,
    2. Thorvaldsdottir H,
    3. Winckler W,
    4. Guttman M,
    5. Lander ES,
    6. Getz G,
    7. Mesirov JP
    (2011) Integrative genomics viewer. Nat Biotechnol 29: 24–26. doi:10.1038/nbt.1754
    OpenUrlCrossRefPubMed
  65. ↵
    1. Kim D,
    2. Paggi JM,
    3. Park C,
    4. Bennett C,
    5. Salzberg SL
    (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37: 907–915. doi:10.1038/s41587-019-0201-4
    OpenUrlCrossRefPubMed
  66. ↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    (2014) HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 31: 166–169. doi:10.1093/bioinformatics/btu638
    OpenUrlCrossRefPubMed
  67. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  68. ↵
    1. Subramanian A,
    2. Tamayo P,
    3. Mootha VK,
    4. Mukherjee S,
    5. Ebert BL,
    6. Gillette MA,
    7. Paulovich A,
    8. Pomeroy SL,
    9. Golub TR,
    10. Lander ES, et al.
    (2005) Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550. doi:10.1073/pnas.0506580102
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Xiao Y,
    2. Hsiao TH,
    3. Suresh U,
    4. Chen HIH,
    5. Wu X,
    6. Wolf SE,
    7. Chen Y
    (2014) A novel significance score for gene selection and ranking. Bioinformatics 30: 801–807. doi:10.1093/bioinformatics/btr671
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Multi-level interaction between HIF and AHR transcriptional pathways in kidney carcinoma
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Crosstalk between HIF and AHR
Véronique N Lafleur, Silvia Halim, Hani Choudhry, Peter J Ratcliffe, David R Mole
Life Science Alliance Feb 2023, 6 (4) e202201756; DOI: 10.26508/lsa.202201756

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Crosstalk between HIF and AHR
Véronique N Lafleur, Silvia Halim, Hani Choudhry, Peter J Ratcliffe, David R Mole
Life Science Alliance Feb 2023, 6 (4) e202201756; DOI: 10.26508/lsa.202201756
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 6, No. 4
April 2023
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Data Availability
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • Design of Parkin activating mutations
  • ppTPP reveals rapid proteome remodeling
  • MicroRNA biomarkers in AD tears
Show more Research Article

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2023 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.