|Home | About | Journals | Submit | Contact Us | Français|
The Krebs cycle enzyme fumarate hydratase (FH) is a human tumor suppressor whose inactivation is associated with the development of leiomyomata, renal cysts, and tumors. It has been proposed that activation of hypoxia inducible factor (HIF) by fumarate-mediated inhibition of HIF prolyl hydroxylases drives oncogenesis. Using a mouse model, we provide genetic evidence that Fh1-associated cyst formation is Hif independent, as is striking upregulation of antioxidant signaling pathways revealed by gene expression profiling. Mechanistic analysis revealed that fumarate modifies cysteine residues within the Kelch-like ECH-associated protein 1 (KEAP1), abrogating its ability to repress the Nuclear factor (erythroid-derived 2)-like 2 (Nrf2)-mediated antioxidant response pathway, suggesting a role for Nrf2 dysregulation in FH-associated cysts and tumors.
► Fh1-associated renal cyst formation is independent of Hif/Phd pathway ► Nrf2-mediated antioxidant response pathway is upregulated following Fh1 loss ► Fumarate modifies cysteine residues in Keap1 by succination
Activation of hypoxia pathways is strongly associated with poor prognosis in cancer. Inactivation of the tumor suppressor gene encoding fumarate hydratase (FH) causes activation of hypoxia signaling, and it has been proposed that this plays a causal role in renal cyst and tumor development. Analyses of mice defective in both FH and HIF signaling reveal that inactivation of Hif-1α but not Hif-2α actually exacerbates renal cyst development. A HIF-independent association has been identified between NRF2 pathway activation and FH deficiency, which correlates with fumarate-mediated modification of cysteine residues in KEAP1, the negative regulator of NRF2. Our data introduce NRF2 dysregulation, rather than HIF activation, as a candidate oncogenic pathway in FH-associated disease.
Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an inherited cancer syndrome in which affected individuals are at risk of developing benign cutaneous and uterine leiomyomas, renal cysts, and aggressive collecting duct and Type 2 papillary renal cell carcinomas (pRCC) (Kiuru et al., 2001; Tomlinson et al., 2002; Toro et al., 2003). The early onset of cysts in hereditary renal cancer syndromes including HLRCC and von Hippel-Lindau (VHL) disease (Kaelin, 2008; Lehtonen et al., 2006) together with observation of dysplastic changes in cystic epithelium and solid tumors arising from cyst walls, strongly suggests that cystic change represents an early stage in carcinogenesis. Patients with HLRCC carry heterozygous germline mutations in the gene encoding the Krebs cycle enzyme fumarate hydratase (FH) and tumor formation is associated with loss of heterozygosity at this locus (Tomlinson et al., 2002). Although FH has a central role in cellular energy metabolism, the mechanisms underlying FH-associated tumorigenesis remain to be defined (Frezza et al., 2011).
Much interest has focused on the possibility that upregulation of one or more transcriptional pathways mediated by hypoxia inducible factor (HIF) may underlie the oncogenic process (Ratcliffe, 2007). In normal cells, oxygen dependent prolyl hydroxylation of HIF-alpha (HIF-α) subunits promotes their degradation by the VHL ubiquitin E3 ligase-proteasome pathway (Kaelin, 2008). In FH-deficient cells and tumors, fumarate accumulates to very high levels, and competitively inhibits the 2-oxoglutarate (2OG) dependent dioxygenases that catalyze HIF prolyl hydroxylation, thus mimicking hypoxia and stabilizing the HIF complex (Isaacs et al., 2005; O'Flaherty et al., 2010; Pollard and Ratcliffe, 2009). This process, termed pseudohypoxia, has been proposed to drive tumor development by activation of oncogenic HIF target genes. In favor of this hypothesis, activation of the HIF transcriptional cascade is strongly associated with aggressive behavior and poor prognosis across a wide range of common human cancers and HIF target genes include many with potentially oncogenic actions in angiogenesis, energy metabolism, invasion and dedifferentiation (Kaelin and Ratcliffe, 2008). HIF is strongly upregulated, both in HLRCC-associated human pRCC and in the hyperplastic renal cysts that develop in mice following targeted inactivation of Fh1 (the murine homolog of FH), and studies of gene expression patterns in these tissues have revealed strong signatures of HIF activation (Ashrafian et al., 2010; Isaacs et al., 2005; Pollard et al., 2005; Pollard et al., 2007). Based on these findings, it has been proposed that pharmacological downregulation of HIF pathways by agents that promote HIF hydroxylation in the face of high fumarate levels might provide an effective treatment for FH-associated neoplasia (Mackenzie et al., 2007; Tennant and Gottlieb, 2010).
It has been shown (Alderson et al., 2006; Frizzell et al., 2011; Nagai et al., 2007) that fumarate, in addition to its role as an allosteric regulator of 2OG-dioxygenase, also modifies cysteine residues in many proteins, forming S-(2-succinyl)-cysteine (2SC). Critically, these modifications have functional consequences as exemplified by inactivation of glyceraldehyde-3-phosphate dehydrogenase (Blatnik et al., 2008a, 2008b). Recently, we have reported that FH-deficient cells and tumors accumulate high levels of 2SC (Bardella et al., 2011). Furthermore, this modification is highly specific and absent in normal tissues and other tumor types and therefore a candidate mechanism for tumorigenesis.
To define the role, if any, of HIF activation in FH-associated neoplasia, we combined inactivation of Fh1 with Hif-1α, Hif-2α, or both Hif-α isoforms, measured the frequency of renal cyst formation in a mouse model recapitulating the cystic phenotype of the human disease, and compared the outcome with that of genetic inactivation of the Hif prolyl hydroxylases (Phds). To extend our analyses and understanding of events underpinning cyst formation following the loss of FH, and to identify potential HIF-independent oncogenic pathways, we compared gene expression patterns in Fh1- and Fh1; Hif-1α-deficient kidneys, where Fh1-associated profiles are not confounded by Hif activation. We provide evidence for an alternative mechanism by which fumarate may activate oncogenic pathways.
To assess the role of HIF activation in FH-associated renal cystic disease, we determined if parallel inactivation of Hif-1α or Hif-2α would ameliorate the hyperplastic renal cystic phenotype in mice with renal tubule specific inactivation of Fh1 (Pollard et al., 2007). Accordingly, mice bearing conditionally inactivated alleles of Hif-1α (Cramer et al., 2003; Higgins et al., 2004), Hif-2α (Gruber et al., 2007), and Fh1 (Pollard et al., 2007) were intercrossed with transgenic animals expressing Cre recombinase under the control of a kidney specific cadherin (Ksp) promoter (Shao et al., 2002) to generate mice that were transgenic for Ksp-Cre and homozygous for one or more conditionally inactivated alleles. A total of seven lines were generated as follows: Fh1fl/fl Ksp-Cre+/− (Fh1−/−), Hif-1αfl/fl Ksp-Cre+/− (Hif-1α−/−), Hif-2αfl/fl Ksp-Cre+/− (Hif-2α−/−), Hif-1αfl/fl Hif-2αfl/fl Ksp-Cre+/− (Hif-1α−/−Hif-2α−/−), Fh1fl/fl Hif-1αfl/fl Ksp-Cre+/− (Fh1−/−Hif-1α-/-), Fh1fl/fHif-2αfl/fl Ksp-Cre+/− (Fh1−/−Hif-2α−/−), and Fh1fl/fl Hif-1αfl/fl Hif-2αfl/fl Ksp-Cre+/− (Fh1−/−Hif-1α−/−Hif-2α−/−). Control animals (Fh1+/+) were littermates bearing conditional alleles without Ksp-Cre or carrying a wild-type allele with Ksp-Cre. For each allele, PCR amplification of genomic DNA from the kidneys was used to verify efficient tissue-specific recombination (Figure 1A).
To determine whether loss of the Hif genes alone in kidney tubules would generate major abnormalities that might confound assessment of cyst development in the combined Fh1−/−Hif-α−/− genotypes, we first analyzed kidneys from control, Hif-1α−/−, Hif-2α−/−, and Hif-1α−/−Hif-2α−/− mice. No major anatomical abnormalities and in particular no cysts were observed by 40 weeks of age in any of these animals (Figure 1B). By comparison, cyst development in Fh1−/− mice is observed from 13 weeks of age (Figure 2A) and is followed by ill health or death from renal failure by 50–65 weeks (Pollard et al., 2007). We therefore conclude that inactivation of Hif-1α or Hif-2α, either alone, or in combination, is not sufficient to initiate cyst formation or to disrupt the renal tubule architecture.
Next, we analyzed kidneys from mice in which Hif-1α, Hif-2α, or both had been deleted in renal tubules in parallel with Fh1. Histological analysis was performed at 13, 17, and 24 weeks of age (Figure 2A). Combined deletion of Fh1 and Hif-1α in Fh1−/−Hif-1α−/− mice did not ameliorate the development of cystic disease as had been postulated. In contrast, parallel inactivation of Hif-1α strikingly accelerated both the initiation and progression of cystic disease compared with mice lacking Fh1 alone, evidenced by increased numbers of dilated tubules and microcysts that progressed to larger and more frequent cysts. By comparison, cysts were never observed in control mice (Figures (Figures1B1B and and2B).2B). Quantification of the numbers of microcysts (>0.1 mm) in kidneys from control, Fh1−/− and Fh1−/−Hif-1α−/− mice at three time points revealed a significant increase in Fh1−/−Hif-1α−/− mice at 13 and particularly at 17 weeks of age as compared with Fh1−/− mice (p < 0.01). This marked exacerbation of the cystic phenotype is also seen in the numbers of macrocysts (>0.5 mm) observed at 17 and 24 weeks of age (Figure 2B).
The development of renal cysts in mice deficient for the Vhl tumor suppressor gene can be repressed by genetic inactivation of the Hif dimerization partner Hif-1β (Arnt) but not Hif-1α; thus suggesting a causal role for Hif-2α in renal cyst development (Rankin et al., 2006). This is in keeping with a potential oncogenic role for HIF-2α in the pathogenesis of clear cell renal cancer (Kaelin, 2007). Since inactivation of Fh1, like Vhl, could potentially result in stabilization of both Hif-α isoforms, we investigated the role of Hif-2α in cyst development in our model. Combined inactivation of Hif-2α and Fh1 failed to ameliorate either the initiation or development of cysts at the three time points analyzed (Figure 2A) and quantification of microcysts and macrocysts showed no differences from Fh1 deficient mice (Figure 2B). These data suggest that distinct mechanisms underlie Vhl- and Fh1-dependent cyst formation.
Although mice subjected to parallel inactivation of Fh1 and Hif-2α were not phenotypically different from animals in which Fh1 alone was inactivated, an inverse relationship between Hif-1α and Hif-2α activity has been proposed previously in some settings (Carroll and Ashcroft, 2006). Therefore, to assess any compensatory or synergistic effects between the two Hif-α isoforms we inactivated both Hif-1α and Hif-2α in combination with Fh1 in the kidney tubules (Figures 2A and 2B). These mice had an apparently identical phenotype to Fh1−/−Hif-1α−/− mice and therefore, at least in the context of this model, Hif-2α does not appear to have a causal role in the cystic phenotype.
HIF prolyl hydroxylation is catalyzed by three closely related enzymes, PHD-1, -2, and -3, (also termed Egln2, -1 and -3) (Epstein et al., 2001). Fumarate, which accumulates in cells lacking FH, mimics hypoxia by allosteric inhibition of the PHDs, allowing HIF to escape destruction and to activate transcription (Isaacs et al., 2005; O'Flaherty et al., 2010). Furthermore, PHDs have been proposed to hydroxylate other substrates in addition to HIF and it is plausible that fumarate-mediated PHD inhibition drives oncogenesis via HIF-independent pathways (Luo et al., 2011; Xie et al., 2009). A number of conditional mouse models exist for the Phd enzymes (Aragonés et al., 2008; Fong and Takeda, 2008; Mazzone et al., 2009). Mice lacking both Phd1 and -3 are viable and develop mild erythrocytosis, although levels of Hif-1α or Hif-2α in the kidney remain normal (Takeda et al., 2008), whereas deletion of Phd2 results in an embryonic lethal phenotype (Takeda et al., 2006). Therefore, to investigate the role of Phds in the regulation of Hif and possibly other relevant pathways affecting cyst formation, we generated mice lacking both Phd1 and −3, (Phd1−/−Phd3−/−) or Phd1, -2, and -3 (Phd1−/−Phd2fl/fl Ksp-Cre+/−Phd3−/− from herein termed Phd1−/−Phd2−/−Phd3−/−). PCR amplification of genomic DNA from the kidneys verified recombination (Figure 1A). Kidney sections were analyzed from Phd knockout mice at two time points, 17 and 24 weeks of age and compared with littermate controls. Whereas kidneys from control animals had no abnormalities, kidneys from Phd1−/−Phd3−/− double knockout mice and Phd1−/−Phd2−/−Phd3−/− triple knockout mice developed subtle tubular cell vacuolization, but no dilation or cyst formation (Figure 3A). Since loss of all three Phds in the kidney failed to induce cyst formation, despite the increased levels of Hif-1α (Figures 3B and 3C), we conclude that cystogenesis in Fh1-deficient mice is both Hif and Phd independent.
To gain a better understanding of the molecular mechanisms causing renal cyst development, comparative genome-wide transcript profiling was performed. Since we observed the largest phenotypic differences in cyst development between control mice and Fh1−/− or Fh1−/−Hif-1α−/− mice these were chosen as the experimental groups between which we compared gene expression profiles. All mice were aged 15 weeks and therefore in the early stages of cystic disease as determined by histological analysis. Pairwise comparisons between the above groups revealed 489 mRNAs to be differentially regulated (see Table S1available online), which were then subjected to IPA pathway analysis (Ingenuity Systems, http://www.ingenuity.com) as in Table S2. The most significant differentially regulated canonical pathway was the Nrf2 -mediated antioxidant pathway which was upregulated in both Fh1−/− and Fh1−/−Hif-1α−/− mice (Figure 4A).
To validate the results obtained from the microarray analysis and to quantify the differential expression of some genes in the Nrf2- mediated gene pathway we analyzed a subset of genes by quantitative reverse transcription polymerase chain reaction (Q-PCR) in control, Fh1−/−, Fh1−/−Hif-1α−/− and Fh1−/−Hif-2α−/− mouse kidneys. First, we quantified the expression of Fh1 and Hif-1α (Figure 4B) and confirmed that Fh1 is significantly reduced (p < 0.05) in the Fh1−/−, Fh1−/−Hif-1α−/− and Fh1−/−Hif-2α−/− mouse kidneys compared with controls and that Hif-1α was significantly reduced (p < 0.05) in the Fh1−/−Hif-1α−/− mouse kidneys compared with controls, Fh1−/− and Fh1−/− Hif-2α−/−. Next, gene expression levels were determined for three of the Nrf2 target genes highlighted in the microarray analysis; glutathione S-transferase alpha 1 (Gsta1), heme oxygenase 1 (Hmox1), and NAD(P)H dehydrogenase [quinone]1 (Nqo1) (Figure 4B). All three transcripts were increased very significantly in Fh1−/−, Fh1−/−Hif-1α−/− and Fh1−/−Hif-2α−/− mouse kidneys compared with controls (p < 0.02). Also, we confirmed that expression of the Hif target gene pyruvate dehydrogenase kinase, isozyme 1 (Pdk1) was increased in Fh1−/− kidneys, but was reduced by parallel Hif-1α or Hif-2α deletion (Figure 4B). To confirm changes in the expression of Fh1, Hif-1α, Nrf2, and Nqo1 at a cellular level, we undertook immunohistochemistry (IHC) on kidney sections from Fh1−/− and Fh1−/−Hif-1α−/− mice at 17 weeks of age. IHC staining confirmed that Fh1 (Figure 4C) is deleted in the cysts, but in both groups continues to be expressed in the interstitium and in a proportion of the renal tubules. Hif-1α is stabilized in the nuclei of cells lining the cysts in Fh1−/− kidneys while this staining is absent in the Fh1−/−Hif-1α−/− kidneys (Figure 4C). Renal cysts from both groups manifest increased nuclear expression of Nrf2; elevated expression in cysts was particularly striking for Nqo1 as compared with the interstitium and some noncystic tubules (Figure 4C). Similarly, strong staining was observed for NRF2 and NQO1 exclusively in the tumor cells and not the stroma in FH-associated Type 2 pRCC (Figure 4D).
Thus, we conclude that expression of Nrf2, and some Nrf2 target genes, are elevated in those kidney tubule cells that have lost Fh1 and line renal cysts, and that, at least in the mouse model, this occurs independently of Hif-1α and Phds and at an early stage of renal cyst formation. Furthermore, elevated levels of these genes can be detected specifically in the tumor cells of HLRCC-associated Type 2 pRCC.
The findings of increased expression of Nrf2 and Nrf2 target genes in FH mutant material from both hyperplastic cysts of knockout mice and human tumors indicate that the Nrf2-mediated antioxidant response is activated during the development of FH-associated disease. To determine if this occurs as a direct consequence of Fh1 inactivation, rather than indirectly in response to subsequent pathology, we analyzed changes in the expression of Nrf2 target genes and other relevant proteins in vitro in wild-type and Fh1−/− mouse embryonic fibroblasts (MEFs) and isogenic stable Fh1−/− transfectants re-expressing human FH (Fh1-/-+FH) (O'Flaherty et al., 2010). In normal cells, the activity of the NRF2 pathway is controlled by KEAP1; a component of an E3 ubiquitin ligase complex that targets NRF2 for degradation (Zhang, 2006). Mutation, deletion, or oxidation of KEAP1 leads to accumulation of nuclear NRF2, enhanced binding to antioxidant response elements (AREs) and activation of downstream target genes (Hayes et al., 2010; Taguchi et al., 2011; Zhang, 2010). To assay the activity of this pathway in MEFs we prepared nuclear and cytoplasmic fractions and analyzed Nrf2, Keap1 and Nqo1 by immunoblotting (Figure 5A). Increased levels of Nrf2 were observed in both the nuclear and cytosolic fractions of Fh1−/− cells compared with either wild-type cells, or to two independent Fh1−/−+FH clones. Whereas the level of Keap1 was unaffected by FH status, Nqo1 was clearly increased in Fh1−/− MEFs relative to the other cell lines.
Since, in addition to NRF2, a complex network of transcription factors, including NRF1, can be recruited to the AREs to modulate transcriptional activity (Biswas and Chan, 2010; Nguyen et al., 2003) we asked whether activation of the antioxidant response observed in Fh1-deficient cells was directly mediated by Nrf2. We therefore depleted Nrf1 or Nrf2 by siRNA in Fh1−/− or wild-type MEFs. The cells were transfected at 0h and 24h and harvested at 48h for analysis. Efficient and specific knockdown of both Nrf1 and Nrf2 was confirmed at the protein level (Figure 5B). In both cell types, the expression levels of Gsta1, Hmox1, and Nqo1 were reduced by depletion of Nrf2, but not by either Nrf1 knockdown, or a non-targeting control (Figure 5C). To determine whether NRF2 also mediates upregulation of the antioxidant pathway in human FH-associated cancer we performed siRNA knockdown of NRF2 in UOK 262 cells, derived from lymph node metastases in an HLRCC patient with aggressive recurring kidney cancer (Yang et al., 2010). Efficient depletion of NRF2 was confirmed by immunoblotting (Figure 5D) and associated with striking reduction of HMOX1 and NQO1 (Figure 5E). Taken together, the data indicate that activation of the Nrf2 pathway occurs as a direct consequence of inactivation of FH in mouse and human cells.
As shown in Figure 4, upregulation of Nrf2 and its target genes Gsta1, Hmox1, and Nqo1 were observed in kidneys from Fh1−/−, Fh1−/− Hif-1α−/−, and Fh1−/− Hif-2α−/− mice, suggesting that this occurs independently of Hif. Similarly, Nrf2 and Nqo1 are both increased in cells lining Fh1-associated renal cysts. To further test any relationship to the Hif/Phd pathway we generated PhdΔ123 MEFs that lack all three of the Phd enzymes. Validation of the genotype and confirmation of upregulation of Hif-1α and Hif target genes in these cells is provided in Figure S1. We then compared expression of Gsta1, Hmox1, and Nqo1 and the Hif target genes Hexokinase 2 (Hk2), Pdk1 and Glucose transporter 1 (Slc2a1) in wild-type, Fh1−/−, Fh1−/−+FH, and PhdΔ123 MEFs. Q-PCR analysis (Figure 5F) demonstrated that while Fh1−/−MEFs express significantly elevated levels of both antioxidant response and Hif-target genes, PhdΔ123 MEFs upregulated Hif-target genes, but not antioxidant response genes. Thus we conclude that upregulation of the antioxidant response pathway in Fh1 deficient cells occurs entirely independently of Hif/Phd dysregulation.
In previous work, we have demonstrated that the intracellular accumulation of fumarate in Fh1−/− cells can be corrected in the face of a persistent defect in mitochondrial oxidative metabolism by the re-expression of an extramitochondrial form of human FH (Fh1−/−+FHΔMTS) (O'Flaherty et al., 2010). Therefore, to distinguish whether Nrf2 activation was a consequence of accumulated fumarate or the mitochondrial defect, we compared Keap1, Nrf2, and Nqo1 protein levels between wild-type, knockout, and Fh1−/−+FHΔMTS MEFs. Only Fh1−/− MEFs showed elevated levels of Nrf2 and Nqo1, whereas the Fh1−/−+FHΔMTS MEFs restored the Nrf2 and Nqo1 levels to those of wild-type MEFs (Figure 5G). Similarly, expression of Gsta1, Hmox1, and Nqo1 were elevated only in the Fh1−/− MEFs as compared with wild-type and Fh1−/−+FHΔMTS MEFs (Figure 5H). Elevation of the Hif target gene Pdk1 was observed in only the Fh1 knockout MEFs (Figure 5H), consistent with previous work (O'Flaherty et al., 2010). Taken together, these findings suggest that activation of Nrf2 signaling occurs as a direct consequence of fumarate accumulation in Fh1−/− cells, rather than as a consequence of defective oxidative metabolism.
Since the NRF2-KEAP1 pathway is regulated by critical cysteine residues in KEAP1 and proteins containing succinated cysteine residues are detected readily in FH-deficient cells and tumors (Bardella et al., 2011), we postulated that the succination of cysteine residues in Keap1 might account for the increased levels of Nrf2 and its downstream targets observed in Fh1 deficient cells and tissues.
We were unable to immunoprecipitate endogenous Keap1 with available antibodies; therefore to determine if Keap1 was indeed subject to succination in Fh1 defective cells we generated Fh1+/+ and Fh1−/− MEFs stably expressing KEAP1-V5 and confirmed cytoplasmic localization of KEAP1 by immunofluorescence (IF) (Figure 6A). To test the function of KEAP1, we analyzed Nrf2 and Nrf2 target gene expression in these cells and demonstrated that as expected, stable expression of KEAP1-V5 reduced levels of Nrf2 protein, and Gsta1, Hmox1, and Nqo1 transcript levels in Fh1+/+ MEFs. In contrast, expression of KEAP1-V5 in Fh1−/− MEFs increased levels of these transcripts indicating that KEAP1 function was strikingly affected by Fh1 status. Furthermore, only KEAP1-V5 immunoprecipitated from Fh1−/− MEFs exhibited strong immunoreactivity for 2SC (Figure 6B, top panel) implying that at least some cysteine residues within KEAP1 were succinated specifically by high levels of fumarate within Fh1−/− cells (Figure 6B).
To identify the site(s) of modification in KEAP1 precisely, V5 immunoprecipitates from the KEAP1-V5 transfectants were subject to enzymatic digestion with trypsin, chymotrypsin or elastase and subsequent analysis by tandem mass spectrometry (UPLC-MS/MS). This revealed mass increments of 116.01 Da corresponding precisely to the predicted mass of the succinyl modification on 17 peptides derived from Fh1−/− but not Fh1+/+ transfectants (Table 1). MS/MS analysis identified succination modification at residues Cys38, Cys151, Cys241, Cys288, Cys319, and Cys613 (Table 1 and Figure 6D). Cys151 and Cys288 have previously been implicated in the regulation of KEAP1 activity by oxidant stress (Nguyen et al., 2009).
Precise quantification of the extent of modification by MS/MS is difficult since the ionization efficiency of modified and unmodified peptide species may differ. However, the addition of two carboxyl residues in the succinated peptide would be predicted to impair ionization, suggesting that succination might be underestimated. The abundance of modified peptides therefore supports the existence of very high levels of succination on KEAP1. Susceptibility of cysteine residues to succination varies, and previous studies have indicated that high pKa thiols such as those found in glutathione (GSH) or N-acetylcysteine are not efficiently targeted for succination (Alderson et al., 2006; Blatnik et al., 2008a). Consistent with this, we observed normal or somewhat elevated levels of reduced GSH in Fh1−/− kidneys, despite high levels of fumarate (Figure S2), therefore suggesting that fumarate does not react directly with GSH.
Taken together, our results suggest that increased levels of fumarate in Fh1−/− cells promote activation of the Nrf2 antioxidant response pathway by succination of specific redox sensitive cysteine residues in KEAP1 (Figure 7) and not by general oxidant stress.
The unexpected demonstration that the gene encoding the Krebs cycle enzyme fumarate hydratase conforms to the classical genetic model of a tumor suppressor, predisposing individuals carrying germline mutations to cancers bearing somatic inactivation of the second allele (Tomlinson et al., 2002), has raised great interest in defining the associated oncogenic pathway(s).
Structural, biochemical, and biological analyses have established that fumarate, which accumulates in FH-defective cells, binds to PHDs and inhibits their catalytic activity leading to upregulation of HIF transcriptional pathways, as occurs in hypoxia (Hewitson et al., 2007; Koivunen et al., 2007). Hypoxia and activation of the HIF system are commonly associated with aggressive cancer (Harris, 2002), but despite intense investigation, cause and effect have remained difficult to distinguish. Since HIF activation is a direct consequence of inactivation of the FH tumor suppressor, irrespective of hypoxia, this link might indicate causality. Indeed, striking activation of HIF was observed in the mouse model described above and in FH-associated human cancer, as well as in tumors linked to inactivation of the succinate dehydrogenase enzyme complex and mutations in genes encoding isocitrate dehydrogenases 1 and 2, which have also been defined as tumor suppressors or oncogenes (Pollard and Ratcliffe, 2009)).
Our findings clearly demonstrate that despite striking activation of Hif and a number of Hif-target genes in Fh1-deficient cells, neither upregulation of Hif nor inactivation of the Phds is required or responsible for the hyperplastic cystic phenotype observed in the mouse model. Surprisingly we found that, rather than ameliorating cyst development, combined inactivation of Hif-1α (but not Hif-2α) and Fh1, greatly exacerbated cystic hyperplasia. Thus, in this setting, upregulation of Hif-1α appears to exert an antiproliferative effect. While this is apparently at odds with the frequently observed upregulation of HIF-1α in cancer, it is consistent with emerging evidence for differential effects of HIF-1α and HIF-2α in tumor biology. HIF-1α antagonizes MYC function, whereas HIF-2α promotes MYC activity (Gordan et al., 2007), and overexpression of HIF-1α and HIF-2α have contrasting effects on the growth of experimental tumors from VHL-defective RCC lines (Raval et al., 2005). Furthermore, mutational analyses reveal a modest but significant prevalence of HIF-1α inactivating mutations in VHL-associated clear cell RCC (Dalgliesh et al., 2010; Morris et al., 2009; Shen et al., 2011). Nevertheless, the finding that inactivation of Hif-2α had no effect on Fh1-associated cystic disease either alone, or in combination with Hif-1α inactivation, differs from findings reported in a similar mouse model of VHL-associated renal neoplasia. In this latter model, combined inactivation of Arnt, but not Hif-1α, ameliorated Vhl-associated renal cystic disease, implying that Hif-2α might be responsible for the cyst development associated with Vhl loss (Rankin et al., 2006). Hence, we conclude that despite the common activation of HIF pathways in VHL- and FH-associated neoplasia, the oncogenic mechanisms are likely to be different.
Combined inactivation of Fh1 and Hif-1α in our mouse model enabled Fh1 dependent changes in transcript profiles to be interrogated without confounding influences from activation of extensive HIF-dependent transcriptional cascades and revealed striking activation of the Nrf2-mediated antioxidant signaling pathway. Further analyses in cell lines derived from Fh1−/− MEFs, mouse cystic tissues and FH-associated human cancer demonstrated that activation of the canonical NRF2 antioxidant pathway arose as a direct consequence of FH inactivation. Though we cannot exclude other influences on Nrf2 dysregulation, the demonstration of high levels of succination on critical cysteine residues in KEAP1, the abnormal activity of transfected KEAP1 in Fh1−/− cells and the maintenance of GSH levels in Fh1−/− cells all argue that Nrf2 activation results from succination of KEAP1 rather than general oxidant stress, at least under the conditions of these experiments (Figure 7). NRF2 acts as a master regulator controlling the ability of mammalian cells to adapt rapidly to stress caused by oxidants and electrophiles, through the induction of ARE containing genes (Nguyen et al., 2003). KEAP1 complexes with Cullin 3 (CUL3) forming an ubiquitin E3 ligase that degrades NRF2. Although not fully understood, the interactions by which KEAP1 controls the levels of NRF2, its cellular localization and transcriptional activity are complex. However, all current models propose that cysteine residues of KEAP1 are modified in response to oxidative stress, resulting in compromised function of the ubiquitin E3 ligase complex that effects proteasomal degradation of NRF2 and enhanced NRF2 stability (Nguyen et al., 2009).
Given that succination of multiple proteins occurs in FH-defective cells and tumors (Bardella et al., 2011), we enquired whether defective KEAP1 function in Fh1 deficient cells might be associated with succination of critical regulatory cysteine residues. MS/MS analysis provided clear evidence of succination on KEAP1 residues Cys151 and Cys288 in association with defective regulation of Nrf2 in Fh1−/− cells. These two residues are among those cysteines, Cys23, 151, 273, 288, and 613, that are conserved between mouse and human and which have been identified as having functional roles in the activity of Keap1 (Hayes et al., 2010; McMahon et al., 2010; Taguchi et al., 2011). Transgenic complementation studies have shown that Cys273 and 288 are essential for Keap1 to repress Nrf2 activity in vivo, while Cys151 is important in facilitating Nrf2 activation in studies with MEFs from a Keap1 (C151S) transgenic mouse model (Yamamoto et al., 2008). Although Cys613 is part of a zinc sensor system (McMahon et al., 2010) and might be modified by fumarate leading to zinc signaling (Cousins et al., 2006) our microarray data and pathway analyses do not suggest that this pathway is dysregulated in Fh1-associated cystic disease.
Interestingly, a recent study has described activation of Nrf2 by exogenous fumarate both in vitro and in vivo (Linker et al., 2011); similar to our findings, they provide direct evidence that KEAP1 is modified at Cys151, though not at Cys288. Whereas this study utilized cell permeable fumaric acid esters (mono- and dimethylfumarate), we have demonstrated that pathophysiological levels of fumarate associated with cancer are sufficient to succinate KEAP1 and activate Nrf2 signaling.
To our knowledge, mutations of KEAP1, NRF2, or downstream target genes, which might shed further light on the role of this pathway in FH-associated oncogenesis, have not yet been described in Type 2 pRCC. Given the extensive transcriptional cascade regulated by NRF2, whether and how dysregulation of KEAP1/NRF2 signaling drives oncogenesis requires further investigation. However, both KEAP1 and NRF2 somatic missense mutations have been identified in a variety of tumors. Moreover, functional assays and the clustering of mutations at sites that disrupt KEAP1/NRF2 regulation have suggested that dysfunction of KEAP1 contributes in some way to oncogenesis in these settings (Hayes and McMahon, 2009; Taguchi et al., 2011). Though KEAP1/NRF2 dysregulation has been considered as an adaptive response that might particularly affect later stages of oncogenesis, recent data in mouse models of pancreatic and lung cancer, where Nrf2 ablation was associated with reduced cellular proliferation, have suggested an early effect (DeNicola et al., 2011). Our data indicating that Nrf2 dysregulation occurs early in the course of hyperplastic cyst development, as a direct consequence of Fh1 inactivation, are consistent with this possibility.
In summary, our investigations have revealed that despite the striking upregulation of the HIF transcriptional cascade in FH-associated neoplasia, these pathways do not appear to contribute to hyperplastic renal cyst formation, at least in a mouse model that recapitulates many features of the human disease; rather, our findings have raised the possibility of an alternative oncogenic action of fumarate through the activation of antioxidant response pathways by succination of KEAP1, and possibly other proteins with tumor suppressor functions.
All experimental procedures were in line with AACR guidelines and passed ethical review by Oxford University's Medical Sciences divisional Local Ethical Review panel. Experiments were performed under UK Home Office regulations, as required by the terms of the Animal (Scientific Procedures) Act 1986. The Fh1, Hif-1α, and Hif-2α conditional knockout and Ksp-cre mice have all been described previously (Gruber et al., 2007; Higgins et al., 2004; Pollard et al., 2007; Shao et al., 2002) as have the Phd1, -2, and -3 knockout mice (Aragonés et al., 2008; Mazzone et al., 2009). Each line had been backcrossed with C57/BL6J for at least five generations and was intercrossed to obtain littermates for the appropriate genotypes. Genotyping was determined by PCR (details on request).
Total RNA was extracted from each sample using a miVana kit (Ambion). The RNA quantity and quality were determined using the Agilent BioAnalyzer 2100 (Agilent Technologies). Gene expression data were obtained by hybridizing a total of 12 mouse samples from three experimental groups: control, Fh1−/− and Fh1−/−Hif-1α−/− (n = 4 per group) to Illumina MouseWG-6 BeadChips. Chips were scanned with Illumina BeadArray Reader; GenomeStudioV2010.1 (Illumina Inc) was used for data extraction. Data was imported to GeneSpring GX 11.0.2 (Agilent Technologies, Inc., Santa Clara, CA) normalized with Shift to 75 percentile and baseline transformed to median of all samples to identify significantly differentially expressed genes which were then subject to IPA pathway analysis (Ingenuity Systems, http://www.ingneuity.com).
The authors are grateful for financial support to UCARE (Oxford), Cancer Research UK; [C6079/A9485] to P.J.R. and P.J.P. and [C10843/A12027] to P.J.P.; The Wellcome Trust; [WT091112MA] to P.J.P., [WT091857MA] to P.J.R. and [090532/Z/09/Z] for core infrastructure support to the Wellcome Trust Centre for Human Genetics, Oxford. P.J.P. is a Beit Memorial Fellow. B.M.K. is supported by the Biomedical Research Centre (NIHR), Oxford, UK. P.C. is grateful for funding from Longterm structural Methusalem funding (Flemish Government); Fondation Leducq - Transatlantic Network Artemis; Belgian Federation against Cancer (196-2008) and Federal Government IUAP Program (P06/30). The authors thank Randall Johnson, Celeste Simon, and Peter Igarashi for providing the Hif-1α, Hif-2α, and Ksp-Cre mice, respectively, and W. Marston Linehan (NCI/NIH) for the UOK 262 cells. The authors are grateful to the staff at the Oxford University Biomedical Services facilities, David Trudgian, Mikael Altun, Rebecca Konietzny, Melroy Miranda, Kaori Igarashi, Kaori Saito, Maki Ohishi, and Sana Ota for technical assistance and to Prof. John Baynes for critical reading of the manuscript. P.J.P. and N.F. have filed for a patent application [USC-268-P(849)] covering immunohistochemical screening with a 2SC antibody for determination of FH mutations. P.J.R., P.H.M., and C.W.P. are cofounders and shareholders in ReOx, Ltd.
Published: October 17, 2011
Supplemental Information includes two tables, two figures, and Supplemental Experimental Procedures and can be found online at doi:10.1016/j.ccr.2011.09.006.
Microarray data analysis was deposited in the GEO public database (GSE29988) (http://www.ncbi.nlm.nih.gov/geo/).