Purpose:

HIF2α is a key driver of kidney cancer. Using a belzutifan analogue (PT2399), we previously showed in tumorgrafts (TG) that ∼50% of clear cell renal cell carcinomas (ccRCC) are HIF2α dependent. However, prolonged treatment induced resistance mutations, which we also identified in humans. Here, we evaluated a tumor-directed, systemically delivered, siRNA drug (siHIF2) active against wild-type and resistant-mutant HIF2α.

Experimental Design:

Using our credentialed TG platform, we performed pharmacokinetic and pharmacodynamic analyses evaluating uptake, HIF2α silencing, target gene inactivation, and antitumor activity. Orthogonal RNA-sequencing studies of siHIF2 and PT2399 were pursued to define the HIF2 transcriptome. Analyses were extended to a TG line generated from a study biopsy of a siHIF2 phase I clinical trial (NCT04169711) participant and the corresponding patient, an extensively pretreated individual with rapidly progressive ccRCC and paraneoplastic polycythemia likely evidencing a HIF2 dependency.

Results:

siHIF2 was taken up by ccRCC TGs, effectively depleted HIF2α, deactivated orthogonally defined effector pathways (including Myc and novel E2F pathways), downregulated cell cycle genes, and inhibited tumor growth. Effects on the study subject TG mimicked those in the patient, where HIF2α was silenced in tumor biopsies, circulating erythropoietin was downregulated, polycythemia was suppressed, and a partial response was induced.

Conclusions:

To our knowledge, this is the first example of functional inactivation of an oncoprotein and tumor suppression with a systemic, tumor-directed, RNA-silencing drug. These studies provide a proof-of-principle of HIF2α inhibition by RNA-targeting drugs in ccRCC and establish a paradigm for tumor-directed RNA-based therapeutics in cancer.

Translational Relevance

In this report, we address the need for new therapeutic strategies in clear cell renal cell carcinoma (ccRCC), the most common type of kidney cancer. We evaluate an innovative, systemically delivered, tumor-directed HIF2α siRNA drug (siHIF2) targeting both wild-type and resistant mutant HIF2α and provide proof-of-principle of activity in mouse models and a patient. Further, using orthogonal RNA-sequencing studies of siHIF2 and a belzutifan-related inhibitor in tumorgrafts, we define the HIF2α transcriptome in ccRCC with unprecedented detail implicating not only traditional (i.e., Myc) but also novel pathways such as E2F2 and E2F8 as well as an 11-gene multicancer network (PBK, BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK) that was previously identified but is enigmatically controlled. These studies provide insight into HIF2α-mediated tumorigenesis and establish a paradigm for RNA-based therapies in cancer.

Clear cell renal cell carcinoma (ccRCC), the most common type of kidney cancer, is characterized by inactivation of the von Hippel–Lindau (VHL) tumor suppressor gene (1). Inactivated either through mutation (80%) or epigenetically (10%), VHL loss is regarded as the signature event of ccRCC (2–5). The VHL protein normally functions to target hypoxia-inducible factor α (HIFα, both HIF1α and HIF2α) for degradation. When VHL is inactivated, HIFα subunits accumulate, bind their HIF1β partner, translocate to the nucleus, and activate gene expression (6). Among HIFα subunits, HIF2α is specifically regarded as the major oncogenic driver of ccRCC (7–11). HIF2 (HIF2α/HIF1β) regulates a plethora of genes promoting cell survival, stemness, proliferation, and angiogenesis (10, 12–14). Particularly prominent among the targets is vascular endothelial growth factor (VEGF), a secreted ligand that binds VEGF receptor 2 (VEGFR2) in endothelial cells promoting angiogenesis (10, 15–18). The VEGF/VEGFR2 axis is of such importance that it is the target of eight FDA-approved drugs for ccRCC to date (19). However, as a more proximal and broader effector, HIF2 would be a more attractive candidate for drug targeting (20).

Historically regarded as undruggable (21, 22), a vulnerability was discovered in the HIF2α structure leading to the identification of small molecules that bound HIF2α, dissociating it from its obligatory partner HIF1β and inhibiting its function (23–25). These chemicals were licensed to Peloton Therapeutics, Inc., which developed several analogue drugs generally referred to as PT drugs (PT2385, PT2399, or the recently FDA-approved PT2977, also called belzutifan). In our previous study, we evaluated PT2399 for the treatment of ccRCC (26). We showed that PT2399 inhibited tumor growth in approximately 50% of ccRCC tumorgrafts (TG). However, when sensitive TG-bearing mice were subjected to prolonged treatment, resistance developed, and we identified an acquired gatekeeper HIF2α mutation (G323E; ref. 26). In vitro, HIF2α (G323E) was sufficient to prevent PT2399 from dissociating HIF2 complexes (26). Subsequently, we identified the same mutation in patients treated with PT2385 in a phase I clinical trial (27). A second resistance mutation was also discovered in HIF1β at the interface with HIF2α (26). Resistance mutations are a well-known phenomenon with targeted therapies (28–30) and complementary therapeutic approaches are therefore needed.

There has been substantial progress with RNA-based and RNA-targeting therapeutics. mRNA-based vaccines have been developed for COVID-19 and the FDA has approved several RNA-targeting drugs for metabolic and other diseases (31, 32). However, to date, no RNA-targeting drug has been approved for cancer therapy. Recently, Arrowhead Pharmaceuticals reported the development of a HIF2α-targeting siRNA drug (33). Consisting of a specific HIF2α-targeting double-stranded (ds)RNA, a dynamic polyconjugate (DPC), and a tumor-directed ligand (an RGD peptide mimetic that interacts with integrins αvβ3 and αvβ5 expressed in ccRCC; refs. 34, 35), their first-generation siHIF2 drug (referred herein as A1HIF2) was taken up by ccRCC cell lines and suppressed tumor growth (33).

Expanding upon these results, we show that A1HIF2 silences HIF2α, suppressing HIF2 target genes, and inhibits the growth of multiple patient-derived HIF2-dependent TG lines from the same platform that previously credentialed PT drugs (26). Using orthogonal approaches integrating siHIF2 and PT2399, we characterize the HIF2 effector transcriptome in unprecedented detail. We show that A1HIF2 has activity against both wild-type and resistant mutant HIF2α. We present studies of a second-generation inhibitor and clinical candidate (A2HIF2), which shares the same HIF2α targeting sequence, but of a significantly smaller size. We provide proof-of-principle of activity in a patient from whom a TG line was generated who participated in the phase I clinical trial (36), where A2HIF2 suppressed not only HIF2α protein in the tumor but also HIF2-dependent paraneoplastic polycythemia and induced a partial response. To our knowledge, this is the first example of functional inactivation of an oncoprotein and tumor suppression with a systemic tumor-directed RNA-silencing drug.

Nomenclature

XP refers to established tumorgraft lines (37, 38). V and A represent, respectively, vehicle- and A1/A2HIF2-treated mice. P refers to PT drug-treated mice. An .x and .y were added to distinguish two vehicle-treated samples from PT2399 and A1HIF2 trials that otherwise had the same name (XP374V1).

Tumorgrafts

Studies were conducted in accordance with the U.S. Common Rule, and informed written consent was obtained from each subject. Investigations were performed after approval by an institutional review board and in accordance with an assurance filed with and approved by the U.S. Department of Health and Human Services. Patient-derived TGs were generated and maintained as previously described (38, 39). Briefly, tumor fragments were implanted orthotopically into 4- to 6-week-old male or female nonobese diabetic severe combined immunodeficient (NOD-SCID) mice according to a UT Southwestern Institutional Animal Care and Use Committee (IACUC)-approved protocol. Mice were monitored weekly for tumor growth by palpation and by magnetic resonance imaging (MRI) as indicated.

A1HIF2 and A2HIF2 targeting sequence

The sequence of A1HIF2 is: sense strand (5′ to 3′): ACGUAACGAUUUCAUGAAAT; antisense strand (5′ to 3′): TUUCAUGAAAUCGUUACGUTT.

The sequence of A2HIF2 is: sense strand (5′ to 3′): CAACGUAACGAUUUCAUGAAA; antisense strand (5′ to 3′): UUUCAUGAAAUCGUUACGUUG.

Drug uptake experiments

NOD/SCID mice were implanted orthotopically with tumor fragments. Once tumors reached ∼100 to 300 mm3, two mice per TG line (n = 4 total) were administered RGD-DPC-Cy3 and one mouse per line (n = 2 total) was administered DPC-Cy3 via intravenous injection. Tumors and normal kidneys were excised 4 hours after injection, and frozen tissue sections were prepared. Tissue sections were counterstained with anti-Na+/K+-ATPase α1 (CST, 23565) for membrane (green) and DAPI for cell nuclei (blue). Cover glass–mounted slides were analyzed using an LSM710 confocal microscope.

Drug trials

Typically, 6 NOD/SCID mice were implanted orthotopically with tissues from the TG lines. Once tumors reached ∼6 mm in diameter, mice were allocated to vehicle or A1HIF2. A1HIF2 was administered at 10 mg/kg in D5W via intravenous injection every 2 weeks. Where applicable, A2HIF2 was administered at 20 mg/kg in D5W via intravenous injection once a week. D5W was used as a control. Tumor growth was evaluated by MRI. Tumor volume was measured by multiplying tumor length, width, and depth to reduce bias. Toxicity was evaluated by monitoring mouse weights once weekly. Mice were euthanized according to IACUC guidelines at the end of drug trials, when they became ill or tumor reached 20 mm in diameter. Tumor tissues, plasma, and serum were collected, processed, and stored as previously described (38, 39).

Immunohistochemistry

Immunohistochemistry (IHC) was performed using a Dako Autostainer Link 48. HIF2α IHC procedures and interpretation were standardized in our laboratory based on expression profiles in well-characterized cell lines and human ccRCC tissues with known expression by western blot analyses as described previously (26). Tissues were stained with HIF2α antibodies (1:1,800, clone A3, Bethyl Laboratories or 1:75, sc-46691, Santa Cruz) using citrate buffer antigen retrieval (K800521, Dako). Appropriate positive and negative controls were used with each run of immunostaining. The percentage of tumor cells staining in the entire section examined was recorded by a trained pathologist. Only nuclear reactivity was regarded as a positive expression.

Real-time PCR

Total RNA was isolated using TRIzol according to the manufacturer's instructions (15596026, Invitrogen). Total RNA (2 μg) for each sample was reverse transcribed using iScript Reverse Transcription Supermix (170-8841, Bio-Rad). RT-PCR was performed on an ABI 7500 Real-Time PCR machine using iTaq Universal SYBR Green SMX (1725124, Bio-Rad). Primers are available upon request.

Western blot

For western blot analyses, HIF2α antibody (A700-003, Bethyl) was diluted at 1:1,000 in 5% BSA in Tris-buffered saline and incubated overnight at 4°C. Tubulin antibody (T5168, Sigma) was diluted at 1:5,000. Anti-rabbit and anti-mouse secondary antibodies (31460 and 31430, Pierce) were used to detect primary antibodies. ECL substrate (1705061, Bio-Rad) was used for signal development.

ELISA

ELISA was performed using a human VEGF ELISA kit (DVE00, R&D Systems) as previously described (26).

RNA-sequencing analyses

For A1HIF2 RNA sequencing (RNA-seq), 11 vehicle- and 11 A1HIF2-treated tumor RNA samples from XP165, XP283, XP289, and XP374 underwent RNA-seq at the New York Genome Center. For PT2399 RNA-seq, 12 vehicle- and 12 PT2399-treated tumor RNA samples from XP144, XP164, XP373, XP374, and XP453 were accessed from Sequence Read Archive SRP073253. Downstream bioinformatic analyses were performed as follows. Briefly, quality control for raw data was accomplished using fastp to filter low-quality reads and trim adaptors (40). Passed reads were aligned to both mouse and human reference sequences with bbmap (41). Mouse reads were filtered out and human reads were remapped to the NCBI hg38 using Hisat2 (41). PCR duplicates were further removed by picard (https://broadinstitute.github.io/picard/). Quantification of genes was performed using HTSeq (42). Differential gene-expression analysis was measured using DESeq2 (43). A false discovery rate (FDR) cutoff of 0.05 was applied to identify statistically significant genes between comparison groups. Gene ontology (GO) enrichment analysis was performed using clusterProfiler in R (44). Gene set enrichment analysis (GSEA) was performed using GSEA software (https://www.gsea-msigdb.org/gsea/index.jsp; ref. 40).

ChIP-seq analyses

Chromatin immunoprecipitation and sequencing (ChIP-seq) data were retrieved from Gene-Expression Omnibus (accession code: GSE120887). Normalized bigwig files of called peaks were loaded into the Integrative Genomics Viewer. Peaks around the regulatory region of a gene locus were used to determine whether the gene may be a direct target of HIF2α.

Statistical analyses

Significance of relative mRNA expression values between A1HIF2 and vehicle groups was determined using Student t tests. Significance of changes in tumor volume or weight was determined using a Mann–Whitney U test. Overlap of gene signatures was evaluated using a hypergeometric distribution. Data were analyzed and plotted using GraphPad Prism.

Data and materials availability statement

The data and materials, including TG lines, generated in this study are available upon request from the corresponding author.

siHIF2 targeting of ccRCC tumorgrafts

We evaluated a first-generation siHIF2 (A1HIF2) developed by Arrowhead Pharmaceuticals (33) and investigated its efficacy using the same TG platform previously credentialing PT drugs (26). We evaluated four different patient-derived ccRCC TGs that were known (or suspected) to be HIF2α dependent and expressed the target integrins αvβ3 and αvβ5 (αvβ3/β5; XP374, XP165, XP283, and XP289; refs. 26, 37).

A summary of patient characteristics and treatment prior to TG line generation is shown in Table 1. XP374 was derived from a 48-year-old man who presented with a 14-cm ccRCC with sarcomatoid dedifferentiation that invaded beyond the kidney. XP165 was generated from an abdominal wall metastasis of a 42-year-old man with ccRCC, who had been previously treated with systemic therapy including an anti-VEGF antibody (bevacizumab). XP283 was derived from a high-grade locally invasive primary ccRCC from a 77-year-old woman who went on to develop metastatic disease. XP289 was derived from a 62-year-old female who presented with metastatic ccRCC and underwent a radical nephrectomy of a 6.5-cm locally advanced tumor.

Table 1.

Characteristics of patients corresponding to each tumorgraft line.

XPIDRaceAgeSexSourceHistologySarcomatoidRhabdoidNucleolar GradeStagePrior therapy
XP165 White 42 Male Abdominal mass ccRCC NA NA – HD-IL2, interferon, bevacizumab 
XP283 White 77 Female Kidney ccRCC Absent NA T3aN0M0 Treatment naïve 
XP289 White 62 Female Kidney ccRCC Absent NA T3aN0M1 Treatment naïve 
XP374 White 48 Male Kidney ccRCC Present Present T4N0M0 Treatment naïve 
XP1487a White 70 Female Lymph node ccRCC Absent Absent – Sunitinib, nivolumab, axitinib, ARO-HIF2 
XPIDRaceAgeSexSourceHistologySarcomatoidRhabdoidNucleolar GradeStagePrior therapy
XP165 White 42 Male Abdominal mass ccRCC NA NA – HD-IL2, interferon, bevacizumab 
XP283 White 77 Female Kidney ccRCC Absent NA T3aN0M0 Treatment naïve 
XP289 White 62 Female Kidney ccRCC Absent NA T3aN0M1 Treatment naïve 
XP374 White 48 Male Kidney ccRCC Present Present T4N0M0 Treatment naïve 
XP1487a White 70 Female Lymph node ccRCC Absent Absent – Sunitinib, nivolumab, axitinib, ARO-HIF2 

Abbreviation: NA, not available.

aTG line from tumor biopsy of phase I trial participant.

To evaluate A1HIF2 tumor targeting, we determined tumor uptake in two different TG lines in the presence and absence of the RGD integrin-targeting ligand. While an RGD-conjugated Cy3-labeled dynamic polyconjugate (RGD-DPC-Cy3) was taken up by the XP283 TG, this was not the case for DPC-Cy3 lacking the RGD peptide mimetic (Supplementary Fig. S1A). Similar results were observed in a second TG line, XP289 (Supplementary Fig. S1B). Of note, RGD-DPC-Cy3 uptake was specific for the tumor with little uptake in normal kidney cells. Thus, the RGD ligand is necessary and sufficient for ccRCC tumor targeting.

Figure 1.

Tumor growth inhibition by siHIF2 in ccRCC TGs. A, Representative T2-weighted MRI images at baseline and after administration of A1HIF2 with corresponding tumor volume quantitation. B, Tumor volume by MRI (n = 3 per arm) at baseline and at the end of treatment. C, Actual tumor volume measurements at the end of treatment. D, Tumor weight at the end of treatment. *, P < 0.05; **, P < 0.01; ****, P < 0.0001. Error bars represent standard deviation.

Figure 1.

Tumor growth inhibition by siHIF2 in ccRCC TGs. A, Representative T2-weighted MRI images at baseline and after administration of A1HIF2 with corresponding tumor volume quantitation. B, Tumor volume by MRI (n = 3 per arm) at baseline and at the end of treatment. C, Actual tumor volume measurements at the end of treatment. D, Tumor weight at the end of treatment. *, P < 0.05; **, P < 0.01; ****, P < 0.0001. Error bars represent standard deviation.

Close modal

Tumor growth inhibition by siHIF2

Next, we evaluated the impact of A1HIF2 on the growth of orthotopically implanted TGs. Cohorts of TG-bearing mice were established and subjected to MRI. Once tumors reached ∼100 to 300 mm3, mice were allocated to either A1HIF2 [10 mg/kg intravenously (i.v.) every two weeks] or vehicle. Tumor growth was monitored longitudinally using MRI and tumors were collected at the end of the study. To assess HIF2α downregulation, IHC analyses were performed. When sufficient tissue was available, the impact of A1HIF2 on HIF2α protein was also measured by western blot.

A1HIF2 inhibited tumor growth by >75% in XP374 (P = 0.0049; Fig. 1A and B; Supplementary Fig. S2A). Similar results were observed by evaluating tumor volumes and weights at the end of the study (Fig. 1C and D). IHC showed significant suppression of HIF2α by A1HIF2 with <5% of the tumor cells staining positively (Fig. 2A). HIF2α suppression was confirmed by western blot (Fig. 2B). Unlike studies we previously performed with other FDA-approved kidney cancer drugs, such as sunitinib (26, 38), A1HIF2 did not induce weight loss (Supplementary Fig. S2B).

Figure 2.

HIF2α silencing and VEGF downregulation by siHIF2. A, Representative IHC images illustrating HIF2α protein depletion by A1HIF2 in TGs. B, Western blot analyses of HIF2α in A1HIF2-treated TG-bearing mice. C, Human VEGF ELISA of circulating tumor-produced VEGF in A1HIF2- or vehicle-treated mice. *, P < 0.05. Error bars represent standard deviation.

Figure 2.

HIF2α silencing and VEGF downregulation by siHIF2. A, Representative IHC images illustrating HIF2α protein depletion by A1HIF2 in TGs. B, Western blot analyses of HIF2α in A1HIF2-treated TG-bearing mice. C, Human VEGF ELISA of circulating tumor-produced VEGF in A1HIF2- or vehicle-treated mice. *, P < 0.05. Error bars represent standard deviation.

Close modal

For XP165, an MRI on day 21 of A1HIF2 treatment showed significant tumor growth inhibition (P = 0.0094; Fig. 1A and B; Supplementary Fig. S2A). Physical analyses revealed a ∼60% reduction in tumor size compared with vehicle-treated mice and a similar reduction in tumor weight (Fig. 1C and D). IHC showed depletion of HIF2α by A1HIF2 with <5% of tumor cells staining positively (Fig. 2A), and results were confirmed by western blot (Fig. 2B). There was no evidence of toxicity (Supplementary Fig. S2B).

For XP283, MRI analyses after 3 doses of A1HIF2 showed >90% reduction in tumor growth (P = 0.017; Fig. 1A and B; Supplementary Fig. S2A), and this was confirmed at the end of the trial (Fig. 1C and D). A1HIF2 suppressed HIF2α expression by IHC (Fig. 2A), but due to the small tumor size, western blot analyses could not be performed on trial tumors. However, significant depletion of HIF2α was achieved by A1HIF2 in this TG line as well (Fig. 2B).

In XP289, A1HIF2 inhibited tumor growth (P = 0.0049) decreasing tumor size and weight by > 60% (Fig. 1; Supplementary Fig. S2A). A1HIF2 suppressed HIF2α expression by IHC (Fig. 2A), and suppression was also observed by western blot (Fig. 2B).

To further evaluate the impact of HIF2α suppression, we measured VEGF, which plays a key role in ccRCC tumorigenesis. We leveraged the species difference between the tumor and host using an ELISA assay specific for human VEGF (hVEGF) and measured the effect of A1HIF2 on circulating VEGF levels. A1HIF2 markedly suppressed tumor VEGF production (Fig. 2C). Although VEGF may also be regulated by HIF1α (15, 45), as expected, quantitative RT-PCR showed suppression of HIF2α (EPAS1) by A1HIF2 but not HIF1α (HIF1A; Supplementary Fig. S3. In addition, no significant change was observed in mRNA levels of HIF1α canonical target genes (LDHA, PGK1; Supplementary Fig. S3). Overall, these data show that A1HIF2 specifically depletes HIF2α in human ccRCC transplants in mice, where it inhibits VEGF production and tumor growth with no appreciable toxicity.

Figure 3.

Orthogonal analyses of siHIF2 and PT2399 in TGs define the HIF2 transcriptome in ccRCC. A, Venn diagram showing shared downregulated genes (FDR < 0.05 and LogFC < −1) by A1HIF2 and PT2399 in TGs (11 vehicle- and 11 A1HIF2-treated tumors from XP165, XP283, XP289, and XP374; and 12 vehicle- and 12 PT2399-treated tumors from XP144, XP164, XP373, XP374, and XP453). B, Heat map of the overlapping 147 significantly downregulated genes. C, Unsupervised GSEA of overlapping 147 genes showing convergence on MYC, E2F, and G2M targets.

Figure 3.

Orthogonal analyses of siHIF2 and PT2399 in TGs define the HIF2 transcriptome in ccRCC. A, Venn diagram showing shared downregulated genes (FDR < 0.05 and LogFC < −1) by A1HIF2 and PT2399 in TGs (11 vehicle- and 11 A1HIF2-treated tumors from XP165, XP283, XP289, and XP374; and 12 vehicle- and 12 PT2399-treated tumors from XP144, XP164, XP373, XP374, and XP453). B, Heat map of the overlapping 147 significantly downregulated genes. C, Unsupervised GSEA of overlapping 147 genes showing convergence on MYC, E2F, and G2M targets.

Close modal

Orthogonal analyses illuminate the HIF2 program in ccRCC

Previous studies of HIF2 in ccRCC tumorigenesis have largely been conducted using RCC tumor cell lines, and how these cell lines reflect the role of HIF2 in RCC tumorigenesis is unknown. We evaluated the impact of A1HIF2 on gene expression by RNA-seq. We performed RNA-seq on 22 TGs from mice implanted with the 4 TG lines (XP165, XP283, XP289, and XP374) treated with either vehicle or A1HIF2. We identified 294 genes that were downregulated at least 2-fold (logFC < −1) at an FDRq < 0.05.

To isolate the effect of HIF2, we performed orthogonal studies leveraging the same ccRCC TG platform and the HIF2 inhibitor, PT2399 (26). We reasoned that by examining the convergence between these two very different HIF2 inhibitors, the HIF2 program may be defined with unprecedented accuracy. Furthermore, because TGs (i) recapitulate the biology of human ccRCC (as shown for example by their coclustering with the corresponding patient tumors in unsupervised gene expression analyses; refs. 37, 38) and (ii) enable direct studies of tumor cells without contamination from the stroma (which is murine), how HIF2 supports ccRCC transformation could be precisely dissected.

We integrated data from A1HIF2 with PT2399 RNA-seq studies involving 12 PT2399-treated and 12 controls from XP144, XP164, XP373, XP374, and XP453 (26). Focusing again on genes downregulated at least 2-fold at an FDRq < 0.05, we identified 419 genes downregulated by PT2399 in TGs. Speaking to the specificity of the two drugs, unsupervised GO analyses of the individual gene lists showed that among the top 10 categories, 8 biological process, 7 cellular components, and 6 molecular functions were shared between A1HIF2 and PT2399 (Supplementary Fig. S4A). KEGG analysis further confirmed these results (Supplementary Fig. S4B).

Figure 4.

First-generation (A1HIF2) and second-generation siHIF2 (A2HIF2). A, Schematic illustration of A1HIF2 and A2HIF2, which share the same RNAi trigger. B, Comparative analyses of HIF2α depletion by A1HIF2 and A2HIF2. C, Depletion of ectopically expressed HIF2α with acquired resistance mutation (G323E) in HEK293T by siHIF2 (A1HIF2). D, Clinical trial schema of A2HIF2 (NCT04169711).

Figure 4.

First-generation (A1HIF2) and second-generation siHIF2 (A2HIF2). A, Schematic illustration of A1HIF2 and A2HIF2, which share the same RNAi trigger. B, Comparative analyses of HIF2α depletion by A1HIF2 and A2HIF2. C, Depletion of ectopically expressed HIF2α with acquired resistance mutation (G323E) in HEK293T by siHIF2 (A1HIF2). D, Clinical trial schema of A2HIF2 (NCT04169711).

Close modal

Remarkably, 50% of the genes downregulated by A1HIF2 (147 out of 294) overlapped with the list of PT2399 downregulated genes (Fig. 3A and B). The probability of such convergence at random was quite low (P < 0.0001). The significance is further highlighted by considering that (i) the studies with A1HIF2 and PT2399 were performed on TG lines that only partially overlapped; (ii) both drugs are administered to mice systemically [one intravenously (A1HIF2) and the other one by gavage (PT2399)]; (iii) a single time point was analyzed; and (iv) the studies were performed several years apart.

An investigation of the 147-gene list showed extensive cell cycle genes (Supplementary Fig. S5A). They included genes encoding cyclins (Cyclin D, Cyclin E, and Cyclin A); cyclin-dependent and other cell cycle kinases such as CDK2, CDK1, and PLK1; origin recognition complex (ORC) proteins (ORC1 and ORC6); mini-chromosome maintenance (MCM) complex proteins (MCM2, MCM3, MCM4, MCM5, and MCM6); ORC/MCM regulatory proteins (CDC6 and CDC45); anaphase promoting complex regulators (CDC20, BUBR1, MAD2, and MPS1); as well as proteins implicated in cytokinesis including ANLN, AURKA, KIF4A, KIF20A, KIF14, NUSAP1, KIF23, KIF20B, RACGAP1, AURKB, and PRC1 (Supplementary Fig. S5B).

Figure 5.

HIF2α inhibition and tumor suppression by siHIF2 in TG from a clinical trial participant (106-00C). A, Schematic illustrating orthotopic TG line generation from ultrasound-guided biopsy (yellow arrows) of a supraclavicular lymph node (white arrowheads) from 106-00C (day 16) along with representative H&E images of the patient core biopsy and corresponding TG (XP1487). B, MR images of orthotopic tumor-bearing TGs treated with A1HIF2 (or vehicle). C, Tumor volume and tumor weight analyses. D, Representative HIF2α IHC images from A1HIF2-treated or control TGs. E, Western blot analysis of HIF2α protein in TGs treated with A1HIF2 or vehicle. F, Integrated analyses of two A1HIF2 trials in TGs from 106-00C (see also Supplementary Fig. S8). *, P < 0.05. Error bars represent standard deviation.

Figure 5.

HIF2α inhibition and tumor suppression by siHIF2 in TG from a clinical trial participant (106-00C). A, Schematic illustrating orthotopic TG line generation from ultrasound-guided biopsy (yellow arrows) of a supraclavicular lymph node (white arrowheads) from 106-00C (day 16) along with representative H&E images of the patient core biopsy and corresponding TG (XP1487). B, MR images of orthotopic tumor-bearing TGs treated with A1HIF2 (or vehicle). C, Tumor volume and tumor weight analyses. D, Representative HIF2α IHC images from A1HIF2-treated or control TGs. E, Western blot analysis of HIF2α protein in TGs treated with A1HIF2 or vehicle. F, Integrated analyses of two A1HIF2 trials in TGs from 106-00C (see also Supplementary Fig. S8). *, P < 0.05. Error bars represent standard deviation.

Close modal

Notably, the 147-gene list included an independently defined and previously reported 11-gene network that was found to be upregulated across 32 cancer types (46). Every member of the list was present (PBK, BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK). Assuming 25,000 genes in the human genome, the total number of 11-gene combinations is ∼6 × 1040, and the probability that a particular such group be found among 147 genes is P = 1.98 × 10−25. Interestingly, all 11 genes contain HIF-binding sites (hypoxia-response elements) in the promoter region [−1,000 to 100 from the transcription start site (TSS)]. These data suggest that this previously defined pan-cancer network is regulated by HIF2.

GSEA of overlapping genes showed a convergence on Myc and E2F programs (Fig. 3C; Supplementary Fig. S4C). While Myc has been extensively implicated as an effector downstream of HIF2 (11, 47), our data implicate E2F in ccRCC development. Both E2F2 and E2F8 were downregulated by A1HIF2 and PT2399. In addition, among the 147 commonly downregulated genes, 44 are putative E2F target genes (Supplementary Table S1).

To further characterize the 147-gene list, we integrated our results with elegant HIF2α ChIP-seq experiments in human kidney proximal tubule (HKC8) cells (48–50). In these experiments, HKC8 cells were edited to disrupt HIF2α, HIF1α, or HIF1β and subjected to ChIP for all 3 proteins (GSE120887; Supplementary Fig. S6A and Supplementary Table S2). To identify putative direct HIF2α target genes, we focused on genes: (i) with predominant HIF2α (compared with HIF1α) binding (TSS ± 5 Kb); (ii) unaffected by HIF1α disruption; and (iii) dampened by HIF2α knockout (Supplementary Fig. S6B and S6C). We found 24 HIF2 target genes, including canonical targets such as SERPINE1 and VEGFA (Supplementary Fig. S6B and C and Supplementary Table S2). However, many genes not previously implicated as putative direct HIF2 target genes in ccRCC were discovered, including cell cycle genes, such as ASPM, AURKA, BLM, CDCA2, CDCA3, CHAF1A, FEN1, GTSE1, KIF15, KNL1, MAD2L1, MCM2, NUF2, and PCLAF.

Figure 6.

HIF2α downregulation, Epo suppression, and tumor growth inhibition by siHIF2 in a clinical trial participant (106-00C). A, Axial (top) and coronal (bottom) CT images of subcarinal lymph node target lesion showing rapid progression during the 2-week washout period and deep prolonged response after two doses of A2HIF2. While first CT scan was performed without iodinated contrast, marked decrease in lymph node enhancement between baseline (second) CT and subsequent CTs illustrates profound antiangiogenic effect. B, Treatment timeline, tumor change (LD, longest diameter), and erythropoietin (EPO) levels prior to and following A2HIF2 administration (red lines). C, H&E and HIF2α IHC of lymph node biopsy at baseline and 16 days after treatment onset with percent of tumor cells positively staining for HIF2α.

Figure 6.

HIF2α downregulation, Epo suppression, and tumor growth inhibition by siHIF2 in a clinical trial participant (106-00C). A, Axial (top) and coronal (bottom) CT images of subcarinal lymph node target lesion showing rapid progression during the 2-week washout period and deep prolonged response after two doses of A2HIF2. While first CT scan was performed without iodinated contrast, marked decrease in lymph node enhancement between baseline (second) CT and subsequent CTs illustrates profound antiangiogenic effect. B, Treatment timeline, tumor change (LD, longest diameter), and erythropoietin (EPO) levels prior to and following A2HIF2 administration (red lines). C, H&E and HIF2α IHC of lymph node biopsy at baseline and 16 days after treatment onset with percent of tumor cells positively staining for HIF2α.

Close modal

To assess the value of studying HIF2 function directly in TGs, we contextualized our findings with previous studies in ccRCC cell lines. We compared our RNA-seq results using the A1HIF2 siRNA to three previously published studies in which HIF2α was depleted in cell lines using sgRNA or siRNA. The first data set (GSE72959), which involves microarray analysis of 786O ccRCC cells after CRISPR-based HIF2α gene (EPAS1) editing (51) identified 929 downregulated genes (Supplementary Fig. S7A). Only 20 genes overlapped with the 294 A1HIF2-downregulated genes (Supplementary Fig. S7A). A second similar data set where EPAS1 was edited in 786O cells and output was evaluated by RNA-seq (GSE149005; ref. 52) identified 355 downregulated genes, but only 18 genes overlapped (Supplementary Fig. S7B). Finally, microarray analysis of A498 cells treated with a HIF2α siRNA (GSE16622; ref. 13) identified 9 downregulated genes and only one that overlapped (Supplementary Fig. S7C). Thus, a very small percentage of genes deregulated by HIF2α targeting in cell lines were validated in TGs. Furthermore, fewer than 10% of genes regulated by HIF2 in ccRCC TGs were identified in experiments using cell lines.

One cell line data set involved orthogonal analyses similar to ours. In the GSE72959 data set, 786O cells were treated with either a HIF2α sgRNA or PT2399. Nine hundred and twenty-nine genes were downregulated by the EPAS1 sgRNA, 48 by PT2399, and 28 overlapped (Supplementary Fig. S7D). Thus, whereas 50% of the genes downregulated by A1HIF2 in TGs were also downregulated by PT2399, only 3% of the genes downregulated by EPAS1 sgRNA were also downregulated by PT2399 in cell lines. When comparing the 28 genes with the 147 overlapping genes in TGs, only one gene overlapped, DEPP1. DEPP1 (also known as DEPP or C10orf10) is a hypoxia-responsive gene (53) implicated in autophagy and is likely a target of HIF2 in ccRCC (54).

HIF2α silencing by siHIF2 clinical candidate (A2HIF2)

Next, we evaluated a second-generation siHIF2, which preserves the RGD targeting moiety and the RNAi trigger sequence (see Materials and Methods), but is otherwise significantly smaller and is the clinical candidate, ARO-HIF2 (referred to here as A2HIF2; Fig. 4A). XP165 TG-bearing mice were treated with A1HIF2 (10 mg/kg in D5W via i.v. injection every 2 weeks), A2HIF2 (20 mg/kg in D5W via i.v. injection once a week), or vehicle. Western blot of tumor samples showed suppression of HIF2α by both A1HIF2 and A2HIF2 (Fig. 4B). However, despite A2HIF2 administration in molar excess compared with A1HIF2, HIF2α appeared to be downregulated to a lesser extent (Fig. 4B). Similar results were observed in other TG lines.

Resistant mutant HIF2α is inhibited by siHIF2

Our previous studies showed that TGs undergoing prolonged treatment with PT2399 acquired resistance and we identified a gatekeeper HIF2α mutation (G323E; ref. 26). The same mutation was subsequently found in patients who developed acquired resistance to the related drug (PT2385) in the phase I clinical trial (27). This mutation is outside the HIF2α siRNA targeting sequence and would not be expected to affect siHIF2 activity. Indeed, A1HIF2 effectively depleted resistant mutant HIF2α (Fig. 4C).

siHIF2 activity in TG from clinical trial participant

A2HIF2 advanced to a phase I clinical trial (Study of ARO-HIF2 in patients with advanced clear cell renal cell carcinoma; NCT04169711), which opened at UT Southwestern Medical Center (UTSW) in 2020. A2HIF2 was evaluated in three sequential dose-escalation cohorts (Fig. 4D) of patients with advanced ccRCC previously treated with conventional antiangiogenic therapy and immune checkpoint inhibitors. The study included pretreatment and on-treatment (after two doses) tumor biopsies for pharmacodynamic analyses.

Biopsy samples from UTSW patients participating in the study were implanted orthotopically in NOD/SCID mice using previously published protocols (38, 39). A post-treatment biopsy specimen from a supraclavicular lymph node of a participant (106-00C) successfully engrafted, and 69 days after implantation, the tumor was passaged into a second cohort of mice, where it grew again. Histologic analyses showed that the TG line (XP1487) retained the features of the parent tumor (Fig. 5A).

XP1487 TG samples were expanded for siHIF2 studies. We observed a >5-fold reduction in tumor size with A1HIF2 (Fig. 5B and C). HIF2α depletion was noted by IHC and confirmed by western blot analyses (Fig. 5D and E). Similar results were obtained with a second trial (Supplementary Fig. S8A and S8B), and a statistically significant difference in tumor growth was observed in combined analyses (Fig. 5F). As expected, A1HIF2 was well tolerated and the mice did not lose weight (Supplementary Fig. S8C).

HIF2α depletion, suppression of paraneoplastic polycythemia, and tumor growth inhibition in clinical trial participant

Patient 106-00C, a 70-year-old female, was diagnosed with metastatic ccRCC in 2017. She had been treated with sunitinib (dose reduced due to toxicity) for 7 months and had received nivolumab for 2 months, when she developed hyperprogression with new onset of brain, bone, and scalp metastases. She received stereotactic radiation to brain and bone metastases and was switched to axitinib, on which she had rapid improvement. Her disease slowly progressed, however, necessitating a gradual increase in axitinib, on which she remained for 2 years. She developed paraneoplastic polycythemia, a syndrome resulting from tumor secretion of erythropoietin (Epo), which induces supranormal red blood cell (RBC) counts (55, 56). Paraneoplastic polycythemia is most often observed with ccRCC, where it is likely driven by HIF2. To reduce RBC counts, phlebotomies were needed as often as every 2 weeks. In January 2021, the patient enrolled in cohort 2 (525 mg weekly) of the ARO-HIF2 phase I trial (NCT04169711). Her disease was particularly aggressive with clinical progression during the 2-week washout period mandated by the study necessitating repeat scans, which showed >30% increase by RECIST v1.1 criteria (Fig. 6A and B). She received two doses of weekly A2HIF2 and a biopsy on day 16 showed downregulation of HIF2α (compared with a pretreatment biopsy; Fig. 6C). Unfortunately, she bled from an intestinal metastasis, which precipitated an acute coronary syndrome, and she came off study. Despite stopping siHIF2 therapy, she had a deep partial response (>60%; Fig. 6A and B). This was accompanied by a substantial decrease in her Epo, which normalized (Fig. 6B). During this time, she enjoyed an outstanding and unprecedented quality of life since her diagnosis. However, Epo levels began to rise over time, and scans showed progression after 4 months. She was started on tivozanib, which reduced her Epo and controlled her disease, at the cost, however, of fatigue, diarrhea, loss of appetite, and hypertension, which are common side effects of tyrosine kinase inhibitors (Fig. 6B).

Here, we show that systemic delivery of a tumor-directed HIF2α siRNA drug resulted in target inactivation and antitumor activity in ccRCC. siHIF2 effectively depleted HIF2α protein in human ccRCC tumors in mice inhibiting target gene expression and VEGF production, and this was associated with substantial tumor growth inhibition. Proof-of-principle of activity is also presented in a patient that participated in the phase I trial.

siHIF2 inhibited not only wild-type HIF2α but also drug-resistant HIF2α. Previously, we reported that belzutifan analogues are associated with resistance mutations, including a gatekeeper mutation, which we first identified in preclinical models, and subsequently reported in patients (26, 27). Resistance mutations are commonplace for targeted therapeutics. With the recent FDA approval of belzutifan, resistance mutations are likely to be more widely identified. Resistance mutations occur not only in HIF2α but also in the HIF1β partner protein that is dislodged by belzutifan (27). HIF1β mutations stabilize the HIF2α/HIF1β complex (57). These mutations illustrate a fundamental dependency on HIF2 in ccRCC, and regardless of whether HIF2α or HIF1β are mutated, HIF2 should be inactivated by siHIF2.

We present proof-of-principle of activity in a patient (106-00C) who participated in the phase I clinical trial of A2HIF2 (36). The patient enrolled in the second planned dose-escalation cohort. Her disease was characterized by sensitivity to antiangiogenic therapy and resistance to immune checkpoint inhibitors and was particularly aggressive with objective progression in the 2-week washout period prior to A2HIF2. She had paraneoplastic polycythemia, which likely indicates a HIF2-driven ccRCC. She experienced a deep partial response, which was accompanied by suppression of Epo, normalization of her hemoglobin, and discontinuation of phlebotomies. IHC studies of pre- and on-treatment tumor biopsies showed HIF2α depletion by siHIF2. That tumor control was due to siHIF2 is supported by the following. First, there was substantial tumor growth preceding siHIF2 but the trend reversed following A2HIF2 administration. Second, tumor regression was accompanied by suppression of Epo. Third, there was tumor growth over time, which was accompanied by an Epo increase, and both responded to tivozanib, a potent FDA-approved ccRCC drug. Finally, her tumor was similarly responsive to siHIF2 in a TG line derived from a biopsy sample. Taken together, these data provide compelling evidence for siHIF2 activity against ccRCC.

In the phase I trial, A2HIF2 was administered weekly. However, in patient 106-00C, the disease remained controlled, and Epo was suppressed for 4 months after just two A2HIF2 infusions. The prolonged activity is consistent with other systemically administered FDA-approved siRNA therapeutics in nononcologic conditions, which can be administered every 3 to 6 months (58). Although we cannot exclude idiosyncrasies about the particular patient, this case study illustrates the complexity of establishing administration schedules for RNA-silencing therapeutics for oncologic indications. In this particular patient, erythropoietin could be used as a pharmacodynamic marker, but paraneoplastic polycythemia is uncommon, and in the majority of patients where tumors do not secrete Epo, Epo could not be used as a biomarker of HIF2α depletion in tumor cells. This emphasizes the need for approaches to evaluate target depletion by siRNA therapeutics. In the case of HIF2α, this could be accomplished by a HIF2α probe, and one has been generated by substituting a fluor atom in PT2385 for 18F (26, 59). Preclinical dual TG models revealed selective uptake in HIF2α-expressing ccRCCs by positron emission tomography, and a clinical trial is ongoing (NCT04989959). Regardless, if generalizable, the data in 106-00C suggest that the A2HIF2 dosing interval could be significantly prolonged.

How many ccRCC patients may benefit from siHIF2 is unclear. In principle, all HIF2-dependent tumors should be responsive to siHIF2. Among these, a minority may have paraneoplastic polycythemia. However, tumors that are sensitive to belzutifan are likely HIF2 dependent. Furthermore, our discovery of HIF2-resistance mutations (26, 27) suggests that HIF2 represents a core dependency. This is unsurprising as ccRCC is thought to be initiated by VHL loss and HIF2 is possibly the major effector downstream. Thus, much like EGFR, where successive generations of inhibitors have been developed targeting resistant mutant EGFR, disabling mutant HIF2, which could be accomplished with siHIF2, should provide additional tumor control. However, the extent to which HIF2-dependent tumors will be responsive to siHIF2 will depend on siHIF2 uptake, which depends in turn on integrin expression (36).

We speculate that, as for targeted therapies, resistance will also arise to siHIF2 drugs. Resistance may involve mutations in the siRNA target sequence. However, RNA-targeting drugs could be redesigned to target other sequences in HIF2α. In addition, multiple sequences could be targeted simultaneously, thereby reducing the chances of such resistance mutations. Nevertheless, other mechanisms of resistance may emerge, leading for example to reduced drug uptake (i.e., downregulating integrin levels).

Given their specificity, siRNA-based therapeutics have the potential to be particularly well tolerated. Currently, FDA-approved drugs targeting angiogenesis downstream of HIF2 cause a plethora of adverse effects including hypertension, diarrhea, thyroid dysfunction, and loss of appetite (19). Many patients lose weight, which is also seen in mice (38). Belzutifan and its analogues are significantly better tolerated, which is likely related to their inhibiting VEGF signaling just in the tumor (in other sites, VEGF production is uncoupled from HIF2; ref. 26). However, belzutifan and its analogues cause anemia and hypoxia (27, 60, 61, 62), which reflect two physiologic processes that engage HIF2, RBC production (63), and respiration (64). Anemia is frequently observed with chronic belzutifan administration and results from suppressing physiologic Epo production (63). However, hypoxia develops infrequently and appears to be associated with particular circumstances such as body habitus or insults such as anesthesia. By directing siRNA to tumor cells, siRNA drugs have the potential to be particularly well tolerated. Nevertheless, toxicities may be observed, as in the case of A2HIF2, that may pertain to other aspects of the drug (36).

This study provides unique insight into HIF2 effector pathways in ccRCC. We leverage a state-of-the-art experimental system characterized by (i) two highly specific inhibitors acting through different mechanisms (PT2399 and siHIF2); (ii) TGs, which faithfully reproduce human ccRCC (37, 38); and (iii) the ability to eliminate confounding effects from the stroma. The importance of the approach is highlighted by comparative studies of tumor cell lines. Whereas in TGs, 50% of the genes downregulated by siHIF2 were downregulated by PT2399, this was the case for only 3% of the genes in cell lines. Furthermore, fewer than 10% of HIF2-regulated genes in ccRCC TGs were previously identified in cell line studies.

Several interesting observations emerged from gene-expression analyses, which defined the HIF2 transcriptome in ccRCC with unprecedented accuracy. First, our data suggest that E2F is an important effector downstream of HIF2. Interestingly, E2F was previously linked to HIF2 in a genetically engineered mouse model (65). Although previous studies were performed in mice and did not characterize the process further, we found that HIF2 regulates the levels of E2F2, which promotes cell cycle entry (66, 67). Our data suggest that HIF2 may be implicated in cell cycle regulation. One unresolved question in ccRCC is where mitogenic signaling comes from. Most tumors have activating mutations in growth factor receptors or downstream signal transduction proteins, but these are rarely mutated in ccRCC. Our data suggest that perhaps HIF2 itself is responsible for inducing cell proliferation, and many cell cycle genes we discovered were previously implicated as putative direct HIF2 targets in normal kidney tubular cells. It is noteworthy that HIF2α was originally identified in endothelial cells where hypoxia, unlike in most other tissues, induces a proliferative response (18). We also found that HIF2 regulates E2F8, which has been implicated in angiogenesis (68). Previously, E2F8 was shown to form a transcriptional complex with HIF1 and induce VEGF expression (68, 69). Given the high degree of similarity between HIF1 and HIF2, HIF2 may similarly interact with E2F8 to promote VEGF expression. Finally, we identified all 11 members of a previously reported multicancer gene network and given our finding of hypoxia-response elements in the promoters of all genes, they may be directly regulated by HIF2.

Another aspect of the gene-expression analyses is the paucity of HIF2-independent effects. Previously, we showed that PT2399 had minimal effects on the transcriptome of ccRCC TGs devoid of HIF2α (26). The low number of genes deregulated by siHIF2 suggests that despite engaging the RNAi machinery to generate the final active siRNA drug, siHIF2 is similarly highly specific.

To our knowledge, this is the first example of target depletion and antitumor activity using a tumor-directed, systemically administered siRNA. This provides a proof-of-concept for what may become a new generation of RNA-silencing cancer therapeutics. The same approach could be deployed to target noncanonical VHL substrates including ZHX2 and SFMBT1 (70, 71) as well as VHL synthetic lethal partners such TBK1, CDK4, and GLS (72–74). In addition, the same paradigm could find application in other tumor types.

A. Christie reports grants from the NIH/NCI during the conduct of the study, as well as grants from CPRIT outside the submitted work. O. Fatunde reports other support from Arrowhead Pharmaceuticals during the conduct of the study. S.C. Wong reports other support from Arrowhead Pharmaceuticals, Inc during the conduct of the study, as well as other support from Arrowhead Pharmaceuticals, Inc outside the submitted work. I. Pedrosa reports personal fees from Merck outside the submitted work. P. Kapur reports personal fees from ClearNano Inc outside the submitted work. T. Schluep reports personal fees from Arrowhead Pharmaceuticals during the conduct of the study, as well as personal fees from Arrowhead Pharmaceuticals outside the submitted work. S.B. Kanner reports a patent for US20160348107 issued. J. Hamilton reports other support from Arrowhead Pharmaceuticals during the conduct of the study, as well as other support from Arrowhead Pharmaceuticals outside the submitted work. J. Brugarolas reports grants and personal fees from Arrowhead Pharmaceuticals during the conduct of the study, as well as personal fees from Eisai, Johnson and Johnson, Exelixis, and Calithera outside the submitted work; in addition, J. Brugarolas has patents for HIF2 biomarkers and mechanisms of resistance pending. No disclosures were reported by the other authors.

Y. Ma: Conceptualization, data curation, formal analysis, validation, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing. A. Joyce: Data curation, investigation, visualization. O. Brandenburg: Data curation, investigation, visualization. F. Saatchi: Data curation, formal analysis, investigation, visualization, writing–original draft, writing–review and editing. C. Stevens: Data curation, investigation, visualization. V. Toffessi Tcheuyap: Data curation, investigation, visualization, project administration. A. Christie: Formal analysis, investigation, visualization. Q.N. Do: Data curation, investigation, visualization. O. Fatunde: Data curation. A. Macchiaroli: Data curation. S.C. Wong: Investigation. L. Woolford: Data curation, investigation, visualization. Q. Yousuf: Data curation. J. Miyata: Data curation. D. Carrillo: Data curation. O. Onabolu: Data curation. T. McKenzie: Data curation. A. Mishra: Formal analysis. T. Hardy: Data curation. W. He: Data curation. D. Li: Data curation. A. Ivanishev: Data curation. Q. Zhang: Formal analysis. I. Pedrosa: Data curation, formal analysis, investigation. P. Kapur: Data curation, formal analysis, investigation, visualization. T. Schluep: Investigation. S.B. Kanner: Investigation. J. Hamilton: Investigation. J. Brugarolas: Conceptualization, resources, formal analysis, supervision, funding acquisition, validation, investigation, methodology, writing–original draft, project administration, writing–review and editing.

We acknowledge the patients who participated in this study and donated samples for research. This work was supported by funding from Arrowhead Pharmaceuticals as well as the UTSW SPORE grant (P50CA196516). Y. Ma, C. Stevens, V. Toffessi Tcheuyap, A. Christie, J. Miyata, D. Carrillo, O. Onabolu, T. McKenzie, I. Pedrosa, P. Kapur, and J. Brugarolas are supported by P50CA196516. F. Saatchi is supported by DoD CDMRP Kidney Cancer Research Program (Award No. W81XWH2110631). Q.N. Do is supported by NICHD K25 award (K25HD104004). Imaging was facilitated by the Small Animal Imaging Research Program supported in part by U24CA126608 and Cancer Prevention and Research Institute of Texas (RP210099). We acknowledge the UTSW Tissue Resource, which is supported in part by the Harold C. Simmons Cancer Center through an NCI Cancer Center Support Grant, 1P30CA142543.

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

1.
Shen
C
,
Kaelin
WG
Jr.
The VHL/HIF axis in clear cell renal carcinoma
.
Semin Cancer Biol
2013
;
23
:
18
25
.
2.
Latif
F
,
Tory
K
,
Gnarra
J
,
Yao
M
,
Duh
F-M
,
Orcutt
ML
, et al
.
Identification of the von Hippel-Lindau disease tumor suppressor gene
.
Science
1993
;
260
:
1317
20
.
3.
Gnarra
JR
,
Tory
K
,
Weng
Y
,
Schmidt
L
,
Wei
MH
,
Li
H
, et al
.
Mutations of the VHL tumour suppressor gene in renal carcinoma
.
Nat Genet
1994
;
7
:
85
90
.
4.
Herman
JG
,
Latif
F
,
Weng
Y
,
Lerman
MI
,
Zbar
B
,
Liu
S
, et al
.
Silencing of the VHL tumor-suppressor gene by DNA methylation in renal carcinoma
.
Proc Natl Acad Sci U S A
1994
;
91
:
9700
4
.
5.
Nickerson
ML
,
Jaeger
E
,
Shi
Y
,
Durocher
JA
,
Mahurkar
S
,
Zaridze
D
, et al
.
Improved identification of von Hippel-Lindau gene alterations in clear cell renal tumors
.
Clin Cancer Res
2008
;
14
:
4726
34
.
6.
Kaelin
WG
Jr
,
Ratcliffe
PJ
.
Oxygen sensing by metazoans: the central role of the HIF hydroxylase pathway
.
Mol Cell
2008
;
30
:
393
402
.
7.
Kondo
K
,
Kim
WY
,
Lechpammer
M
,
Kaelin
WG
.
Inhibition of HIF2alpha is sufficient to suppress pVHL-defective tumor growth
.
PLoS Biol
2003
;
1
:
E83
.
8.
Zimmer
M
,
Doucette
D
,
Siddiqui
N
,
Iliopoulos
O
.
Inhibition of hypoxia-inducible factor is sufficient for growth suppression of VHL−/− tumors
.
Mol Cancer Res
2004
;
2
:
89
95
.
9.
Shen
C
,
Beroukhim
R
,
Schumacher
SE
,
Zhou
J
,
Chang
M
,
Signoretti
S
, et al
.
Genetic and functional studies implicate HIF1alpha as a 14q kidney cancer suppressor gene
.
Cancer Discov
2011
;
1
:
222
35
.
10.
Keith
B
,
Johnson
RS
,
Simon
MC
.
HIF1alpha and HIF2alpha: sibling rivalry in hypoxic tumour growth and progression
.
Nat Rev Cancer
2011
;
12
:
9
22
.
11.
Gordan
JD
,
Lal
P
,
Dondeti
VR
,
Letrero
R
,
Parekh
KN
,
Oquendo
CE
, et al
.
HIF-alpha effects on c-Myc distinguish two subtypes of sporadic VHL-deficient clear cell renal carcinoma
.
Cancer Cell
2008
;
14
:
435
46
.
12.
Covello
KL
,
Kehler
J
,
Yu
H
,
Gordan
JD
,
Arsham
AM
,
Hu
C-J
, et al
.
HIF-2alpha regulates Oct-4: effects of hypoxia on stem cell function, embryonic development, and tumor growth
.
Genes Dev
2006
;
20
:
557
70
.
13.
Bertout
JA
,
Majmundar
AJ
,
Gordan
JD
,
Lam
JC
,
Ditsworth
D
,
Keith
B
, et al
.
HIF2alpha inhibition promotes p53 pathway activity, tumor cell death, and radiation responses
.
Proc Natl Acad Sci U S A
2009
;
106
:
14391
6
.
14.
Roberts
AM
,
Watson
IR
,
Evans
AJ
,
Foster
DA
,
Irwin
MS
,
Ohh
M
.
Suppression of hypoxia-inducible factor 2alpha restores p53 activity via Hdm2 and reverses chemoresistance of renal carcinoma cells
.
Cancer Res
2009
;
69
:
9056
64
.
15.
Ema
M
,
Taya
S
,
Yokotani
N
,
Sogawa
K
,
Matsuda
Y
,
Fujii-Kuriyama
Y
.
A novel bHLH-PAS factor with close sequence similarity to hypoxia-inducible factor 1alpha regulates the VEGF expression and is potentially involved in lung and vascular development
.
Proc Natl Acad Sci U S A
1997
;
94
:
4273
8
.
16.
Compernolle
V
,
Brusselmans
K
,
Acker
T
,
Hoet
P
,
Tjwa
M
,
Beck
H
, et al
.
Loss of HIF-2alpha and inhibition of VEGF impair fetal lung maturation, whereas treatment with VEGF prevents fatal respiratory distress in premature mice
.
Nat Med
2002
;
8
:
702
10
.
17.
Turner
KJ
,
Moore
JW
,
Jones
A
,
Taylor
CF
,
Cuthbert-Heavens
D
,
Han
C
, et al
.
Expression of hypoxia-inducible factors in human renal cancer: relationship to angiogenesis and to the von Hippel-Lindau gene mutation
.
Cancer Res
2002
;
62
:
2957
61
.
18.
Tian
H
,
McKnight
SL
,
Russell
DW
.
Endothelial PAS domain protein 1 (EPAS1), a transcription factor selectively expressed in endothelial cells
.
Genes Dev
1997
;
11
:
72
82
.
19.
Motzer
RJ
,
Jonasch
E
,
Boyle
S
,
Carlo
MI
,
Manley
B
,
Agarwal
N
, et al
.
NCCN guidelines insights: kidney cancer, version 1.2021
.
J Natl Compr Canc Netw
2020
;
18
:
1160
70
.
20.
Brugarolas
J
.
Renal-cell carcinoma–molecular pathways and therapies
.
N Engl J Med
2007
;
356
:
185
7
.
21.
Koehler
AN
.
A complex task? Direct modulation of transcription factors with small molecules
.
Curr Opin Chem Biol
2010
;
14
:
331
40
.
22.
Lazo
JS
,
Sharlow
ER
.
Drugging undruggable molecular cancer targets
.
Annu Rev Pharmacol Toxicol
2016
;
56
:
23
40
.
23.
Key
J
,
Scheuermann
TH
,
Anderson
PC
,
Daggett
V
,
Gardner
KH
.
Principles of ligand binding within a completely buried cavity in HIF2alpha PAS-B
.
J Am Chem Soc
2009
;
131
:
17647
54
.
24.
Rogers
JL
,
Bayeh
L
,
Scheuermann
TH
,
Longgood
J
,
Key
J
,
Naidoo
J
, et al
.
Development of inhibitors of the PAS-B domain of the HIF-2alpha transcription factor
.
J Med Chem
2013
;
56
:
1739
47
.
25.
Scheuermann
TH
,
Tomchick
DR
,
Machius
M
,
Guo
Y
,
Bruick
RK
,
Gardner
KH
.
Artificial ligand binding within the HIF2alpha PAS-B domain of the HIF2 transcription factor
.
Proc Natl Acad Sci U S A
2009
;
106
:
450
5
.
26.
Chen
W
,
Hill
H
,
Christie
A
,
Kim
MS
,
Holloman
E
,
Pavia-Jimenez
A
, et al
.
Targeting renal cell carcinoma with a HIF-2 antagonist
.
Nature
2016
;
539
:
112
7
.
27.
Courtney
KD
,
Ma
Y
,
Diaz de Leon
A
,
Christie
A
,
Xie
Z
,
Woolford
L
, et al
.
HIF-2 complex dissociation, target inhibition, and acquired resistance with PT2385, a first-in-class HIF-2 inhibitor, in patients with clear cell renal cell carcinoma
.
Clin Cancer Res
2020
;
26
:
793
803
.
28.
Fisher
R
,
Pusztai
L
,
Swanton
C
.
Cancer heterogeneity: implications for targeted therapeutics
.
Br J Cancer
2013
;
108
:
479
85
.
29.
Housman
G
,
Byler
S
,
Heerboth
S
,
Lapinska
K
,
Longacre
M
,
Snyder
N
, et al
.
Drug resistance in cancer: an overview
.
Cancers (Basel)
2014
;
6
:
1769
92
.
30.
Dagogo-Jack
I
,
Shaw
AT
.
Tumour heterogeneity and resistance to cancer therapies
.
Nat Rev Clin Oncol
2018
;
15
:
81
94
.
31.
Crooke
ST
,
Baker
BF
,
Crooke
RM
,
Liang
X-H
.
Antisense technology: an overview and prospectus
.
Nat Rev Drug Discov
2021
;
20
:
427
53
.
32.
Paunovska
K
,
Loughrey
D
,
Dahlman
JE
.
Drug delivery systems for RNA therapeutics
.
Nat Rev Genet
2022
;
23
:
265
80
.
33.
Wong
SC
,
Cheng
W
,
Hamilton
H
,
Nicholas
AL
,
Wakefield
DH
,
Almeida
A
, et al
.
HIF2alpha-targeted RNAi therapeutic inhibits clear cell renal cell carcinoma
.
Mol Cancer Ther
2018
;
17
:
140
9
.
34.
Vogetseder
A
,
Thies
S
,
Ingold
B
,
Roth
P
,
Weller
M
,
Schraml
P
, et al
.
αv-Integrin isoform expression in primary human tumors and brain metastases
.
Int J Cancer
2013
;
133
:
2362
71
.
35.
Wechsel
HW
,
Petri
E
,
Feil
G
,
Nelde
HJ
,
Bichler
KH
,
Loesr
W
, et al
.
Renal cell carcinoma: immunohistological investigation of expression of the integrin alpha v beta 3
.
Anticancer Res
1999
;
19
:
1529
32
.
36.
Brugarolas
J
,
Beckermann
K
,
Rini
BI
,
Vogelzang
NJ
,
Lam
ET
,
Hamilton
JC
, et al
.
Initial results from the phase 1 study of ARO-HIF2 to silence HIF2-alpha in patients with advanced ccRCC (AROHIF21001)
.
J Clin Oncol
2022
;
40
(
6_suppl
):
339
.
37.
Elias
R
,
Tcheuyap
VT
,
Kaushik
AK
,
Singla
N
,
Gao
M
,
Reig Torras
O
, et al
.
A renal cell carcinoma tumorgraft platform to advance precision medicine
.
Cell Rep
2021
;
37
:
110055
.
38.
Sivanand
S
,
Pena-Llopis
S
,
Zhao
H
,
Kucejova
B
,
Spence
P
,
Pavia-Jimenez
A
, et al
.
A validated tumorgraft model reveals activity of dovitinib against renal cell carcinoma
.
Sci Transl Med
2012
;
4
:
137ra75
.
39.
Pavía-Jiménez
A
,
Tcheuyap
VT
,
Brugarolas
J
.
Establishing a human renal cell carcinoma tumorgraft platform for preclinical drug testing
.
Nat Protoc
2014
;
9
:
1848
59
.
40.
Chen
S
,
Zhou
Y
,
Chen
Y
,
Gu
J
.
fastp: an ultra-fast all-in-one FASTQ preprocessor
.
Bioinformatics
2018
;
34
:
i884
-
i890
.
41.
Bushnell
B
.
BBMap: a fast, accurate, splice-aware aligner
.
In: 9th Annual Genomics of Energy and Environment Meeting
,
Walnut Creek, CA
,
March 17–20, 2014
.
2014
:
United States. p. Medium: ED
.
42.
Anders
S
,
Pyl
PT
,
Huber
W
.
HTSeq–a Python framework to work with high-throughput sequencing data
.
Bioinformatics
2015
;
31
:
166
9
.
43.
Love
MI
,
Huber
W
,
Anders
S
.
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
.
Genome Biol
2014
;
15
:
550
.
44.
Yu
G
,
Wang
L-G
,
Han
Y
,
He
Q-Y
.
clusterProfiler: an R package for comparing biological themes among gene clusters
.
OMICS
2012
;
16
:
284
7
.
45.
Hu
C-J
,
Wang
L-Y
,
Chodosh
LA
,
Keith
B
,
Simon
MC
.
Differential roles of hypoxia-inducible factor 1alpha (HIF-1alpha) and HIF-2alpha in hypoxic gene regulation
.
Mol Cell Biol
2003
;
23
:
9361
74
.
46.
Liu
Y
,
Xiang
J
,
Peng
G
,
Shen
C
.
Omics- and pharmacogenomic evidence for the prognostic, regulatory, and immune-related roles of PBK in a pan-cancer cohort
.
Front Mol Biosci
2021
;
8
:
785370
.
47.
Gordan
JD
,
Bertout
JA
,
Hu
C-J
,
Diehl
JA
,
Simon
MC
.
HIF-2alpha promotes hypoxic cell proliferation by enhancing c-myc transcriptional activity
.
Cancer Cell
2007
;
11
:
335
47
.
48.
Smythies
JA
,
Sun
M
,
Masson
N
,
Salama
R
,
Simpson
PD
,
Murray
E
, et al
.
Inherent DNA-binding specificities of the HIF-1alpha and HIF-2alpha transcription factors in chromatin
.
EMBO Rep
2019
;
20
.
49.
Schmid
V
,
Lafleur
VN
,
Lombardi
O
,
Li
R
,
Salama
R
,
Colli
L
, et al
.
Co-incidence of RCC-susceptibility polymorphisms with HIF cis-acting sequences supports a pathway tuning model of cancer
.
Sci Rep
2019
;
9
:
18768
.
50.
Lauer
V
,
Grampp
S
,
Platt
J
,
Lafleur
V
,
Lombardi
O
,
Choudhry
H
, et al
.
Hypoxia drives glucose transporter 3 expression through hypoxia-inducible transcription factor (HIF)-mediated induction of the long noncoding RNA NICI
.
J Biol Chem
2020
;
295
:
4065
78
.
51.
Cho
H
,
Du
X
,
Rizzi
JP
,
Liberzon
E
,
Chakraborty
AA
,
Gao
W
, et al
.
On-target efficacy of a HIF-2alpha antagonist in preclinical kidney cancer models
.
Nature
2016
;
539
:
107
11
.
52.
Hong
K
,
Hu
L
,
Liu
X
,
Simon
JM
,
Ptacek
TS
,
Zheng
X
, et al
.
USP37 promotes deubiquitination of HIF2alpha in kidney cancer
.
Proc Natl Acad Sci U S A, 2020.
117
:
13023
32
.
53.
Klee
K
,
Storti
F
,
Maggi
J
,
Todorova
V
,
Karademir
D
,
Berger
W
, et al
.
The expression of decidual protein induced by progesterone (DEPP) is controlled by three distal consensus hypoxia responsive element (HRE) in hypoxic retinal epithelial cells
.
Genes (Basel)
2020
;
11
:
111
.
54.
Stepp
MW
,
Folz
RJ
,
Yu
J
,
Zelko
IN
.
The c10orf10 gene product is a new link between oxidative stress and autophagy
.
Biochim Biophys Acta
2014
;
1843
:
1076
88
.
55.
Morais
C
,
Johnson
DW
,
Vesey
DA
,
Gobe
GC
.
Functional significance of erythropoietin in renal cell carcinoma
.
BMC Cancer
2013
;
13
:
14
.
56.
Palapattu
GS
,
Kristo
B
,
Rajfer
J
.
Paraneoplastic syndromes in urologic malignancy: the many faces of renal cell carcinoma
.
Rev Urol
2002
;
4
:
163
70
.
57.
Wu
D
,
Su
X
,
Lu
J
,
Li
S
,
Hood
BL
,
Vasile
S
, et al
.
Bidirectional modulation of HIF-2 activity through chemical ligands
.
Nat Chem Biol
2019
;
15
:
367
76
.
58.
Ray
KK
,
Landmesser
U
,
Leiter
LA
,
Kallend
D
,
Dufour
R
,
Karakas
M
, et al
.
Inclisiran in patients at high cardiovascular risk with elevated LDL cholesterol
.
N Engl J Med
2017
;
376
:
1430
40
.
59.
Courtney
KD
,
Infante
JR
,
Lam
ET
,
Figlin
RA
,
Rini
BI
,
Brugarolas
J
, et al
.
Phase I dose-escalation trial of PT2385, a first-in-class hypoxia-inducible factor-2α antagonist in patients with previously treated advanced clear cell renal cell carcinoma
.
J Clin Oncol
2018
;
36
:
867
74
.
60.
Choueiri
TK
,
Bauer
TM
,
Papadopoulos
KP
,
Plimack
ER
,
Merchan
JR
,
McDermott
DF
, et al
.
Inhibition of hypoxia-inducible factor-2alpha in renal cell carcinoma with belzutifan: a phase 1 trial and biomarker analysis
.
Nat Med
2021
;
27
:
802
5
.
61.
Choueiri
TK
,
Plimack
ER
,
Bauer
TM
,
Merchan
JR
,
Papadopoulos
KP
,
McDermott
DF
, et al
.
Phase I/II study of the oral HIF-2 α inhibitor MK-6482 in patients with advanced clear cell renal cell carcinoma (RCC)
.
J Clin Oncol
2020
;
38
(
6_suppl
):
611
.
62.
Jonasch
E
,
Donskov
F
,
Iliopoulos
O
,
Rathmell
WK
,
Narayan
VK
,
Maughan
BL
, et al
.
Belzutifan for renal cell carcinoma in von Hippel-Lindau disease
.
N Engl J Med
2021
;
385
:
2036
46
.
63.
Haase
VH
.
Regulation of erythropoiesis by hypoxia-inducible factors
.
Blood Rev
2013
;
27
:
41
53
.
64.
Cheng
X
,
Prange-Barczynska
M
,
Fielding
JW
,
Zhang
M
,
Burrell
AL
,
Lima
JDCC
, et al
.
Marked and rapid effects of pharmacological HIF-2alpha antagonism on hypoxic ventilatory control
.
J Clin Invest
2020
;
130
:
2237
51
.
65.
Hoefflin
R
,
Harlander
S
,
Schäfer
S
,
Metzger
P
,
Kuo
F
,
Schönenberger
DSE
, et al
.
HIF-1alpha and HIF-2alpha differently regulate tumour development and inflammation of clear cell renal cell carcinoma in mice
.
Nat Commun
2020
;
11
:
4111
.
66.
Beier
R
,
Bürgin
A
,
Kiermaier
A
,
Fero
M
,
Karsunky
H
,
Saffrich
R
, et al
.
Induction of cyclin E-cdk2 kinase activity, E2F-dependent transcription and cell growth by Myc are genetically separable events
.
Embo J
2000
;
19
:
5813
23
.
67.
Joyce
NC
,
Harris
DL
,
McAlister
JC
,
Ali
RR
,
Larkin
DFP
.
Effect of overexpressing the transcription factor E2F2 on cell cycle progression in rabbit corneal endothelial cells
.
Invest Ophthalmol Vis Sci
2004
;
45
:
1340
8
.
68.
Weijts
BGMW
,
Bakker
WJ
,
Cornelissen
PWA
,
Liang
K-H
,
Schaftenaar
FH
,
Westendorp
B
, et al
.
E2F7 and E2F8 promote angiogenesis through transcriptional activation of VEGFA in cooperation with HIF1
.
Embo J
2012
;
31
:
3871
84
.
69.
Bakker
WJ
,
Weijts
BGMW
,
Westendorp
B
,
de Bruin
A
.
HIF proteins connect the RB-E2F factors to angiogenesis
.
Transcription
2013
;
4
:
62
6
.
70.
Liu
X
,
Simon
JM
,
Xie
H
,
Hu
L
,
Wang
J
,
Zurlo
G
, et al
.
Genome-wide screening identifies SFMBT1 as an oncogenic driver in cancer with VHL loss
.
Mol Cell
2020
;
77
:
1294
306
.
71.
Zhang
J
,
Wu
T
,
Simon
J
,
Takada
M
,
Saito
R
,
Fan
C
, et al
.
VHL substrate transcription factor ZHX2 as an oncogenic driver in clear cell renal cell carcinoma
.
Science
2018
;
361
:
290
5
.
72.
Hu
L
,
Xie
H
,
Liu
X
,
Potjewyd
F
,
James
LI
,
Wilkerson
EM
, et al
.
TBK1 is a synthetic lethal target in cancer with VHL loss
.
Cancer Discov
2020
;
10
:
460
75
.
73.
Nicholson
HE
,
Tariq
Z
,
Housden
BE
,
Jennings
RB
,
Stransky
LA
,
Perrimon
N
, et al
.
HIF-independent synthetic lethality between CDK4/6 inhibition and VHL loss across species
.
Sci Signal
2019
;
12
:
eaay0482
.
74.
Okazaki
A
,
Gameiro
PA
,
Christodoulou
D
,
Laviollette
L
,
Schneider
M
,
Chaves
F
, et al
.
Glutaminase and poly(ADP-ribose) polymerase inhibitors suppress pyrimidine synthesis and VHL-deficient renal cancers
.
J Clin Invest
2017
;
127
:
1631
45
.