Abstract
Loss of the retinoblastoma (RB) tumor suppressor protein is a critical step in reprogramming biological networks that drive cancer progression, although mechanistic insight has been largely limited to the impact of RB loss on cell-cycle regulation. Here, isogenic modeling of RB loss identified disease stage–specific rewiring of E2F1 function, providing the first-in-field mapping of the E2F1 cistrome and transcriptome after RB loss across disease progression. Biochemical and functional assessment using both in vitro and in vivo models identified an unexpected, prominent role for E2F1 in regulation of redox metabolism after RB loss, driving an increase in the synthesis of the antioxidant glutathione, specific to advanced disease. These E2F1-dependent events resulted in protection from reactive oxygen species in response to therapeutic intervention. On balance, these findings reveal novel pathways through which RB loss promotes cancer progression and highlight potentially new nodes of intervention for treating RB-deficient cancers.
This study identifies stage-specific consequences of RB loss across cancer progression that have a direct impact on tumor response to clinically utilized therapeutics. The study herein is the first to investigate the effect of RB loss on global metabolic regulation and link RB/E2F1 to redox control in multiple advanced diseases.
This article is highlighted in the In This Issue feature, p. 2113
Introduction
The retinoblastoma (RB) tumor suppressor protein serves as a transcriptional corepressor that is frequently altered in cancer and has been linked to both prevention of tumor development and progression. Canonical RB function in cell-cycle control is tightly regulated by cyclin-dependent kinase (CDK) and cyclin phosphorylation as well as upstream CDK inhibitor (CDKI) proteins. In response to mitogenic stimuli, CDK/cyclin complexes hyperphosphorylate RB, promoting a conformational change and allowing for derepression of activator E2F transcription factors (E2F1, E2F2, and E2F3; refs. 1, 2). Loss of function (LOF) of RB alters gene networks, including but not limited to those associated with E2Fs, which have been shown to drive tumorigenesis and disease progression. Common alterations driving RB LOF in human malignancy include RB1 gene mutation, promoter methylation, and deep deletion, in addition to upstream pathway modulation such as upregulation of cyclins and CDKs or loss of CDKI function (3, 4). These alterations have been identified across cancers and proven to play a significant role in driving tumor phenotypes. Despite current understanding of RB function in preventing cellular transformation, significant gaps in understanding exist in discerning the role of RB in tumor progression.
Elucidating the role of RB loss in disease progression is critical, as RB LOF has been attributed to poor outcome in a number of tumor types. Specifically, RB loss has been associated with shorter overall survival (OS) in multiple myeloma, bladder cancer, and early-stage non–small cell lung cancer (NSCLC; refs. 5–8). In addition, RB loss has been associated with higher tumor grade and stage in bladder cancer and shorter event-free survival (EFS) in osteosarcoma (9, 10). Beyond the loss of RB protein, RB-loss gene signatures have been examined in breast and prostate cancers. Elevated expression of the RB-loss gene signature has been associated with shorter OS and relapse-free survival (RFS) in breast cancers and shorter progression-free survival (PFS), disease-specific survival (DSS), and OS in prostate cancer (11, 12). Other studies have linked RB loss to therapeutic relapse. In hormone receptor–positive breast cancers and prostate cancer, RB loss has been associated with tumor relapse and hormone therapy resistance, highlighting the clinical importance of RB loss on disease progression (13–15). Thus, although loss of RB function clearly can promote tumor initiation, the mechanisms by which RB functions to induce the aggressive tumor phenotypes that have been observed in human malignancie remain largely unstudied.
While the tumor-suppressive functions of RB in cancer affect cell-cycle control, clinical data indicate that functions beyond mitotic regulation likely occur in human disease. Strikingly, in cancers wherein RB loss is strongly associated with poor outcome, such as lung and prostate cancers, there is no correlation between RB loss and hyperproliferative activity (16, 17). Studies examining RB loss and Ki-67 positivity in lung adenocarcinoma and castration-resistant prostate cancer (CRPC) found no correlation between RB protein negativity and Ki-67 positivity, indicating that RB regulatory control expands beyond the cell cycle (16, 17). On the basis of these findings and the importance of RB in tumor progression, it is critical to determine the cellular and molecular mechanisms by which RB loss promotes disease progression in advanced cancers.
To address this important gap in understanding, studies herein have discovered unrealized functions of RB in controlling tumor metabolism. Isogenic models of RB loss in early-stage and advanced disease were developed, and subsequent molecular assessment identified distinct, stage-specific activities of the RB tumor suppressor protein. Key findings revealed that RB loss in advanced disease redirected E2F1 to serve as a critical regulator of glutathione synthesis, as evident through cistrome, transcriptome, and metabolome profiling. Further investigation revealed that these activities significantly increase glutathione both in vitro and in vivo, and conferred a resulting protection from reactive oxygen species (ROS) generation in response to cytotoxic chemotherapy. These studies shift paradigms in thinking with regard to metabolic control and reveal entirely new activities through which RB and E2F1 function in tumor suppression.
Results
Transcriptional Networks Governed by RB Depletion Are Disease Stage–Dependent
The biological significance of RB loss in driving progression to CRPC is well established (15, 18). RB loss is infrequent in localized disease, yet is found in 19% to 35% of CRPC and is associated with poor patient outcome (18–22). Importantly, the molecular consequence of RB loss on disease progression extends beyond cell-cycle control, highlighting a critical role for RB in regulating a myriad of protumorigenic pathways (16). Analyses of currently available prostate cancer patient cohorts revealed that 72% of RB1 gene alterations are deletions (Fig. 1A, top), and examination across hormone therapy–sensitive prostate cancer (HSPC) and CRPC cohorts revealed the frequency of this alteration is increased to an average of 10.9% in CRPC compared with 2.4% in HSPC, indicating a distinct function for RB1 depletion in CRPC (Fig. 1A, bottom). While RB loss induces the CRPC transition, castration resistance can be activated through multiple mechanisms, including alternative androgen receptor (AR) splicing, posttranslational modifications, cofactor perturbation, mutations, and intracrine androgen synthesis (23). Here, isogenic modeling was used to assess the relative impact of RB loss in HSPC versus existing CRPC from the identical lineage (Fig. 1B, top). Thus, HSPC models mimic the “RB1 depletion–induced CRPC” stage of disease (RBD-induced CRPC) as described previously (15, 16), while CRPC isogenic models represent the “RB1 depletion-post CRPC transition” (RBD-post CRPC), a more advanced stage of disease.
As RB function is frequently disrupted in cancer, the specific CDK4/6 inhibitor palbociclib was used to validate functional RB status in parental models utilized herein (24–28). Flow cytometry was performed to measure BrdU incorporation to rigorously assess the number of cells in S-phase in the absence or presence of palbociclib and revealed that both HSPC and CRPC control models showed a significant and similar decrease in S-phase incorporation following CDK4/6 inhibition (Supplementary Fig. S1A). In contrast, RB1-deficient variants proved resistant to palbociclib, confirming loss of RB function. Furthermore, Western blot analyses of total and phosphorylated RB protein revealed similar levels of phosphorylated RB present across both models that are both reduced following CDK4/6 inhibition (Supplementary Fig. S1B), thus indicating similar RB functional status.
As AR and E2F1 have both been shown to be sensitive to RB modulation in previous studies (15, 16), the effects of RB1 depletion on AR and E2F1 protein expression were examined across all models under castrate- and androgen-stimulated conditions. In castrate conditions, RB1 depletion promoted a greater than 1.2-fold increase in both AR and E2F1 protein levels in RBD-induced CRPC models, consistent with previous publications (15, 16). Conversely, no significant changes in AR and E2F1 protein expression were found in RBD-post CRPC models (Fig. 1B, bottom). Similar changes in AR and E2F1 expression were found in androgen-stimulated conditions showing a slight increase in AR and E2F1 protein expression in RBD-induced CRPC models and no change in RBD-post CRPC models (Fig. 1B, bottom). These data highlight a distinct consequence of RB1 depletion in the transition from HSPC to CRPC; however, they do not clarify the significance of RB1 depletion in CRPC, particularly the impact of RB loss on the protumorigenic functions of E2F1 and AR where no significant change in protein expression was observed. As expected, RBD-induced CRPC models exhibited a growth advantage in castrate conditions (ref. 15; Fig. 1C, top). However, there was no change in growth in the RBD-post CRPC models (Fig. 1C, bottom), further underscoring a distinct role for RB loss pre- and post-CRPC transition.
RB is a well-described transcriptional regulator; as such, RNA sequencing (RNA-seq) was performed in both RBD-induced CRPC and RBD-post CRPC isogenic models after androgen stimulation, so as to understand stage-specific differences. Consistency of biological replicates was confirmed by principal component analysis (PCA; Supplementary Fig. S1C). As shown, while RBD-induced CRPC exhibited significant changes in 4,313 transcripts (Fig. 1D, left), RBD-post CRPC resulted in 7,480 significantly altered transcripts (Fig. 1D, right;Padj < 0.05), suggesting a distinct role for RB loss in transcriptional regulation at this later stage of disease. To identify distinct transcriptional alterations regulated by RB between stages of disease, transcriptomic overlay was performed (Fig. 1E). Specifically, differentially expressed transcripts were identified compared with stage-specific controls in each model (FC > 1.5, Padj < 0.05). Expression changes were stratified as described in Fig. 1E: exclusively altered in RBD-induced CRPC (1,081 transcripts, light blue), exclusively altered in RBD-post CRPC (1,645 transcripts, light green), commonly altered in both stages (34 transcripts, gray), and two inversely regulated categories: RBD-induced CRPC upregulated, RBD-post CRPC downregulated (795 transcripts, dark blue) and RBD-post CRPC upregulated, RBD-induced CRPC downregulated (1,104 transcripts, dark green). Significantly, the clear differences in transcriptional regulation, including the notable number of inversely regulated transcripts, imply both distinct as well as potentially opposing functions of RB across stages of disease, which has not been previously described. Finally, to identify the putative transcriptional networks driven by these distinct gene-expression changes across stages of disease, gene set enrichment analysis (GSEA) was performed, revealing that pathways significantly affected by RBD-induced CRPC included immune response and cell-cycle control–related pathways, whereas those driven by RBD-post CRPC were primarily related to metabolic and growth factor signaling (Fig. 1F). Together these data highlight, for the first time, both the divergent transcriptional networks driven by RB1 depletion between stages of disease and the biologically relevant pathways associated with these changes.
RB1 Depletion–Induced E2F1 Rewiring is Disease Stage–Specific
It has been previously reported that RB loss in HSPC results in rewiring of the E2F1 cistrome, and the resulting transcriptional output is associated with poor outcome (16). To discern the relative impact of RB loss in distinct disease states, E2F1 function was assessed via chromatin immunoprecipitation sequencing (ChIP-seq) across disease progression (Fig. 2A, top). Whereas 32,235 E2F1 binding sites were identified in the stage-specific control, RBD-induced CRPC cells exhibited a reduction in total E2F1 binding sites with 23,950 sites identified after RB1 depletion, with 3,597 of those sites exclusive to the RB1-depleted condition (Fig. 2A, middle). Conversely, RBD-post CPRC induced a marked gain of total E2F1 binding compared with the stage-specific control, increasing the number of E2F1 sites from 13,594 to 31,647, with 18,624 of these sites being exclusive to the RB1-depleted condition (Fig. 2A, bottom). This increase in E2F1 binding suggests a putative gain in E2F1 function when RB is depleted after the CRPC transition, likely contributing to the transcriptional changes seen in this stage as reported in Fig. 1.
As RB loss has previously been shown to alter E2F1 preference for binding at canonically described promoter regions, genomic annotation of E2F1 binding sites was performed for each disease stage. As shown in Fig. 2B, RBD-induced CRPC resulted in a slight increase in the percent of E2F1 binding at promoter regions compared with stage-specific controls (Fig. 2B, top). Interestingly, RBD-post CRPC resulted in a significant shift in binding at promoter regions, decreasing from 40% to 23% of total binding, and at intronic and distal intergenic regions, increasing from a combined 40% to 59% of total binding (Fig. 2B, bottom), suggesting that E2F1 binding shifts toward intronic and distal intergenic regions in CRPC following RB loss. These differential shifts in E2F1 binding further underscore the disparate effects of RB loss pre- and post-transition to CRPC. Given the shift in the genomic elements associated with E2F1 binding after RB1 depletion, de novo motif analysis was performed to gain insight into the potential mechanism governing these changes. Interestingly, the most significant and highly enriched motif closely resembled BORIS, a known transcription factor that binds at promoters of genes that are frequently elevated in cancer (ref. 29; Fig. 2C, left). Furthermore, RB1-depleted exclusive E2F1 binding in CRPC was also enriched for motifs resembling BORIS in addition to the pioneer factor FOXM1 and the NF1-halfsite motif (Fig. 2C; Supplementary Fig. S2), suggesting this gained E2F1 binding may be associated with the function of these additional DNA binding proteins. To further identify potential novel cofactors associated with E2F1, specifically those gained exclusively after RB1 depletion in CRPC (compared with stage-specific controls), known motif analysis was performed. Analysis within a 1 kb window from the center of gained E2F1 binding revealed the most significantly enriched motifs included Foxa3, Foxa2, and Fox:Ebox forkhead motifs (Fig. 2D), supporting the potential for E2F1 association with forkhead proteins at these novel binding sites. These data provide insight into the mechanism and putative cofactors driving E2F1 to novel DNA binding sites after RB1 depletion in CRPC.
Given that RB depletion in existing CRPC resulted in distinct expansion of the E2F1 cistrome and concomitant transcriptional rewiring, the RBD-post CRPC–associated E2F1 cistrome was mapped to putative E2F1 target genes. Genes with transcriptional start sites (TSS) within 30 kb of the center of E2F1 binding after RB1 depletion were identified (17,393 genes). To prioritize genes for downstream investigation, these binding-associated genes were intersected with those exhibiting an increase in gene expression, specifically those seen to be exclusively or inversely regulated (1,981 genes) compared with RBD-induced CRPC. This overlay identified 1,599 genes that had E2F1 binding within 30 kb of the TSS in addition to increased expression after RB1 depletion distinct from RBD-induced CRPC. GSEA revealed that 27% of all pathways enriched were related to cell metabolism, with the most significantly enriched pathways involved in sphingolipid and amino acid metabolism (Fig. 3A). These data indicate that E2F1 binding in RBD-post CRPC likely drives a specific increase of E2F1-driven expression of genes involved in metabolic pathways and implicates RB as a regulator of cellular metabolism in late-stage disease.
RB Loss in Late-Stage Disease Results in E2F1-Dependent Metabolic Reprogramming
As E2F1 was identified as a potential, major metabolic regulator in RBD-post CRPC, the assortment of putative metabolic pathways regulated were investigated further. Target genes across pathways including GMPR (purine metabolism), PSAT1 (glycine serine and threonine metabolism), CHPT1 (glycerphospholipid metabolism), CERK (sphingolipid metabolism), and GSTA1 (glutathione metabolism) were confirmed to have a change in mRNA expression and E2F1 binding after RB1 depletion. The mRNA expression levels of GMPR, CHPT1, CERK, and GSTA1 were increased 1.2- to 2.0-fold after RB depletion (Fig. 3B, top). When validated in an additional RBD-post CRPC model in which RB1 was downregulated via miRNA, all selected target genes displayed 2.0- to 3.5-fold increases in expression (Fig. 3B, middle). Furthermore, expanding beyond prostate cancer, mRNA expression of these targets was confirmed in an invasive breast cancer model (Fig. 3B, bottom), showing a 1.5- to 3.0-fold increase in mRNA expression after RB1 depletion. Validation of RB, AR, and E2F1 protein levels within these additional models is also shown (Fig. 3B, right). These mRNA expression data highlight a significant role for RB in regulating the expression of a diverse range of metabolic genes functioning in numerous metabolic pathways across cancer types. To further confirm that these expression changes are likely to be regulated by E2F1 after RB1 depletion, the E2F1 binding sites identified within 30 kb of these gene TSSs were confirmed with a second antibody recognizing a distinct epitope of E2F1 (30). E2F1 binding was validated across targets, showing a greater than 2-fold increase in binding after RB1 depletion compared with controls (Fig. 3C). These data define a distinct role for RB and E2F1 in regulating a collection of metabolic pathways in advanced disease.
Because these data demonstrate a novel role for RB in regulating cell metabolism in CRPC, to assess biological outcome, steady-state metabolomic analysis was performed in isogenic RBD-post CRPC models. Metabolomics identified 89 significantly altered metabolites after RB1 depletion compared with stage-specific controls. KEGG pathway analysis of significantly altered metabolites identified the top pathways modulated after RB1 depletion as lipid metabolism (49%), amino acid metabolism (28%), and peptide metabolism (12%; Fig. 4A; Supplementary Fig. S3A). These data support a role for RB in regulating distinct lipid and amino acid metabolic pathways. Analysis of significantly altered metabolites within their metabolic pathways (Fig. 4B) showed a range of upregulated and downregulated lipid metabolites and pathways after RB1 depletion (Supplementary Fig. S3B). By contrast, there was a distinct upregulation of amino acid metabolites, with the majority exhibiting a 1.5-fold or greater increase in abundance (increased metabolites, red; decreased metabolites, blue; Fig. 4C). These data highlight a clear role for RB in regulation of metabolic function in RBD-post CRPC, with a particular effect in amino acid and lipid metabolism.
To prioritize specific metabolic pathways most likely to be directly regulated by E2F1 after RB1 depletion, pathway analysis from ChIP-seq, RNA-seq, and metabolomics were integrated as shown in Fig. 4D. This comparison included the pathways enriched from E2F1 binding–associated genes (31,647 binding sites, green circle), the pathways identified from the transcriptional alterations either exclusively or inversely regulated compared with RBD-induced CRPC models (1,981 genes, blue circle), and the pathways identified via metabolomic analyses after RB1 depletion in CRPC (89 metabolites, red circle). Integration of these complex datasets revealed seven metabolic pathways commonly altered after RB loss in late-stage disease as indicated by novel E2F1 regulation, transcriptional changes, and altered metabolites (Fig. 4D). These pathways included five amino acid synthesis pathways including glutamate metabolism and glutathione metabolism, both of which are of strong cancer significance (31). Together, these data identify an unexpected consequence of RB1 depletion on cancer cell metabolic control, reprogramming lipid and amino acid metabolism.
RB1 Depletion in Late-Stage Disease Protects against ROS through Increased Glutathione Synthesis
Data presented herein indicate E2F1 gain of function after RB1 depletion post-CRPC significantly alters metabolic control, resulting in increased amino acid and glutathione synthesis. Glutathione, an antioxidant that functions by neutralizing intracellular ROS, has been shown to protect cells from DNA-damaging agents and promote therapy resistance (32); thus, this pathway was prioritized for further investigation. The synthesis of the tripeptide glutathione requires the amino acids glutamate, cysteine, and glycine, which can be either transported into the cell by membrane transporters or metabolized from other sources such as amino acids. Analysis of three common cancer-associated amino acid transporters (SLC1A5/ASCT2, alanine-serine-cysteine transporter 2; LAT1/SLC7A5, l-type amino acid transporter 1; and xCT/SLC7A11, cystine/glutamate transporter), showed a significant increase in expression after RBD-post CRPC (Fig. 5A, left, red). ASCT2 imports glutamine, which is rapidly converted to glutamate by intracellular glutaminase for both glutamate supply and as a glutamate amine donor in the serine/glycine pathway. xCT exchanges glutamate to facilitate import of cystine, which can subsequently be reduced to cysteine within the cell. LAT1 transports methionine into the cell, which can ultimately be converted to cysteine through the trans-sulfuration pathway (33–36). Validation of the mRNA expression of SLC1A5/ASCT2 confirmed a greater than 1.6-fold increase in CRPC, breast cancer, NSCLC, and bladder cancer models after RB1 depletion. Furthermore, the mRNA of multiple crucial enzymes within the pathway were also found to be increased following RB1 depletion in CRPC. Cystathionase (CTH), the enzyme that converts cystathionine, the output of the transsulfuration pathway, to cysteine was confirmed to have a 1.6-fold mRNA increase in CRPC and NSCLC, and a 1.2-fold increase in bladder cancer models. Finally, gamma-glutamylcysteine synthetase (GCLC), the rate-limiting enzyme required for the first step of glutathione synthesis, was confirmed to have an average 4.1-fold increase in CRPC and a 2-fold increase in breast cancer, NSCLC, and bladder cancer models. The protein expression for each of these validated pathway components was confirmed with a 1.3-fold increase in GCLC and ASCT2 and a 2.2-fold increase in CTH following RB1 depletion in CRPC (Fig. 5B). Furthermore, metabolites within the pathway were also changed following RB1 depletion. Cystathionine, the precursor for the reaction catalyzed by CTH, was decreased, whereas glutamine and gamma-glutamylcysteine were increased (Fig. 5A, left, blue), highlighting a change in enzymatic activity within the pathway. To further examine the impact of RB loss on glutathione metabolism, the endpoint of the pathway, total glutathione was measured in CRPC, breast cancer, NSCLC, and bladder cancer disease models. RB1 depletion resulted in an average of a 1.5-fold increase in total glutathione across all models (Fig. 5C), further supporting an increase in glutathione synthesis following RB1 depletion in advanced disease and expanding the impact of these findings. These data implicate a significant role for RB in regulating the glutathione synthesis pathway at several key enzymatic points and across multiple advanced disease types.
Given that RB1 depletion results in an increase in the antioxidant glutathione, the effect of RB1 depletion on intracellular ROS was examined. Basal levels of intracellular ROS were decreased by 30% after RB1 depletion in CRPC models and 60% in breast cancer. Furthermore, when treated with 50 μmol/L menadione, a known inducer of ROS, intracellular ROS continued to display a reduction in the RB1-depleted condition (Fig. 5D), suggesting that the increase in glutathione after RB1 depletion lowers intracellular ROS, potentially acting to protect cells from full ROS induction and downstream effects. To examine the response to clinically relevant therapeutics, intracellular ROS was measured in a time-dependent manner following treatment with the cytotoxic agent doxorubicin. In response to doxorubicin, intracellular ROS peaked between 6 and 8 hours after treatment, while this gradual increase and ROS induction was significantly reduced in the absence of RB (Fig. 5E), further indicative of a protective role for RB LOF through regulation of glutathione production. In sum, these data indicate that RB1 depletion in CRPC and invasive breast cancer drives altered metabolic pathways, thereby promoting an increase in glutathione synthesis. This increase in glutathione significantly reduced basal intracellular ROS and protected from ROS induction following treatment with ROS inducers, including doxorubicin, highlighting a protective mechanism in cancer cells against cytotoxic agents currently used to target advanced disease.
Glutathione Synthesis is Directly Regulated by RB Loss–Induced E2F1 Function and is a Candidate for Therapeutic Intervention
To validate that these biological changes are a direct consequence of RB1 loss, PSM.7-LP, an established, constitutively activate RB, was utilized to rescue RB function in RB1-depleted models (37). Use of such a strategy is required, as it is well appreciated that the ectopic expression of wild-type RB is rapidly inactivated by endogenous CDKs and does not rescue RB loss (38, 39). As such, shRNA-resistant PSM.7-LP variants were generated and subsequently transfected into the RB1-depleted models. Expression of endogenous RB and PSM.7-LP was verified by immunoblot showing depleted full-length RB and the presence of PSM.7-LP (Fig. 6A, left). Expression of E2F1 was decreased by 1.4-fold in RB1 depletion models following transfection with PSM.7-LP. This anticipated decrease is a consequence of the presence of constitutively active RB, further validating a functional protein (Fig. 6A, left). To examine the effects of rescuing active RB on glutathione synthesis, mRNA expression of E2F1 and glutathione synthesis genes was measured following transfection of PSM.7-LP. Supporting the changes observed at the protein level, expression of PSM.7-LP led to a 2.1-fold reduction in E2F1 mRNA. Furthermore, mRNA expression of CTH, GCLC, and SLC1A5/ASCT2 was also significantly reduced by 1.5- to 2.5-fold following rescue of active RB (Fig. 6A, right). These data validate that glutathione synthesis genes are regulated directly by RB and indicate a positive correlation between E2F1 and the expression of CTH, GCLC, and SLC1A5/ASCT2. To further confirm that RB directly regulates glutathione synthesis, total glutathione was also measured following the expression of PSM.7-LP (Fig. 6B). Supporting the transcript data, total glutathione was significantly reduced by 1.5-fold with the rescue of active RB, confirming that RB is directly regulating glutathione synthesis in CRPC.
To delineate the underlying mechanism by which RB1 depletion drives increased glutathione synthesis and thus protection from ROS-inducing therapeutics, the potential for E2F1 regulation of glutathione synthesis was examined. ChIP-seq analyses indicated that E2F1 binding was enhanced after RB1 depletion at genes encoding for glutathione synthesis pathway components (Supplementary Fig. S4A). This trend was validated via ChIP-qPCR at the promoters of CTH, GCLC, and SLC1A5/ASCT2 where a 5-fold, 2.2-fold, and 5-fold increase in E2F1 binding was observed after RB1 depletion, respectively. This increase in E2F1 binding at glutathione synthesis–related gene promoters was abolished following the knockdown of E2F1 (Fig. 6C). Enhanced E2F1 binding was accompanied by an increase in relative mRNA (>2-fold change) and protein expression (>1.2-fold change) of these targets (Fig. 6D, left and right, respectively). Importantly, these alterations were abrogated by knockdown of E2F1 (Fig. 6D; Supplementary Fig. S4B), suggesting a direct role for E2F1 regulation of glutathione synthesis pathway components by E2F1 after RB1 depletion in advanced CRPC. Further supporting this concept, the marked 2-fold increase in total glutathione following RB1 depletion was largely abrogated after knockdown of E2F1 (Fig. 6E). Overall, these data provide a novel link between E2F1 and glutathione synthesis in advanced disease. These data reveal that RB loss confers an E2F1-dependent metabolic advantage in CRPC, with potential implications for therapeutic efficacy.
To further investigate the impact of RB regulation on glutathione metabolism, biological function was assessed throughout the pathway. To assess the initial step of glutamine transport, [3H]-l-glutamine uptake was measured before and after RB1 depletion. RB1 depletion alone resulted in a significant 1.6-fold increase in glutamine uptake, which could be inhibited using the ASCT2 inhibitor l-γ-glutamyl-p-nitroanilide (GPNA; Fig. 6F). Levels of leucine uptake, which is not an ASCT2 substrate, were unchanged after RB1 depletion (Supplementary Fig. S5), indicating that there was a specific increase in ASCT2-mediated glutamine uptake after RB1 depletion and that this increase can be therapeutically targeted. To target this pathway further, total glutathione and cytotoxicity of each isogenic model was assessed in response to the xCT (SLC7A11) inhibitor Erastin and the GCLC inhibitor buthionine sulfoximine (BSO; Fig. 6G, left). Erastin and BSO significantly reduced total glutathione in both RB1-depleted and control models by an average of 11-fold and 2.7-fold, respectively, verifying that both inhibitors directly target glutathione synthesis in the models utilized (Fig. 6G, right). Critically, inhibition of xCT (SLC7A11) or GCLC resulted in a significant decrease in cell viability at varying concentrations in RB1-depleted models of CRPC and breast cancer when compared with the RB1-intact controls (Fig. 6H). These data indicate that RB1 depletion leads to increased sensitivity to inhibition of glutathione synthesis, suggesting that RB1-depleted tumors are more reliant on this pathway for survival. In addition, these data reveal novel avenues to clinically target RB1-deficient tumors.
The consequences of RBD-post CRPC include a rearrangement of transcriptional networks, rewiring of E2F1 function, and reprogramming of cellular metabolism. More specifically, RB1 depletion has been shown to promote a gain in E2F1 binding at the promoters of genes that encode principal components of the glutathione synthesis pathway, with expression and protein changes in these genes dependent on E2F1. Moreover, increased expression of the glutamine transporter SLC1A5/ASCT2 was demonstrated to drive an increase in glutamine uptake in RBD-post CRPC cells. Further downstream in the pathway, gained E2F1 function after RB1 depletion induced an increase in total glutathione within cancer cells, and inhibition at various points within the pathway showed increased sensitivity in the RB1-depleted models. To validate these observations in vivo, tumor xenograft models were utilized. Specifically, isogenic RBD-post CRPC cell models were subcutaneously injected into athymic nude mice. When tumors reached 500 mm3, tumors were harvested and assessed for glutathione, mRNA, and protein. RB1 depletion within the tumors was validated via IHC (Fig. 7A). In agreement with in vitro studies, there was no change in tumor doubling time nor time to tumor take after RB1 depletion in CRPC in vivo (Supplementary Fig. S6). However, total glutathione in the RB1-depleted tumors was significantly increased (2.1-fold) compared with control (Fig. 7B). Further validating in vitro findings, mRNA expression of GCLC and SLC1A5/ASCT2 was significantly increased (>1.2 FC; Fig. 7C), while protein expression examined via IHC was also increased following RB1 depletion in vivo (Fig. 7D; Supplementary Fig. S7). Thus, these findings further validate the role of RB LOF in glutathione metabolic control in CRPC.
Importantly, the impact of RB on genes controlling glutathione metabolism was assessed in a novel and unselected metastatic CRPC tumor cohort of 98 patient samples from the Royal Marsden Hospital (London, United Kingdom). RNA-seq was performed to identify AR-positive tumors with and without intact RB. The expression of GCLC, CTH, SLC1A5, and SLC7A11 was examined across the cohort and compared with E2F1. There was no correlation between the expression of E2F1 and any of the glutathione genes investigated (GCLC, CTH, SLC1A5, and SLC7A11) in tumors with intact RB (Fig. 7E, left). However, significant positive correlations were observed between E2F1 and expression of all four glutathione genes in the RB-loss tumors, GCLC (Spearman correlationP = <0.001), CTH (P = 0.02), SLC1A5 (P = <0.001), and SLC7A11 (P = 0.02; Fig. 7E, right). These correlations validate our findings that RB loss drives an increase in glutathione synthesis through E2F1 in late-stage CRPC tumors.
Together, these data are the first to reveal distinct biological and molecular consequences of RB1 depletion in HSPC and CRPC. This study identifies a novel function for E2F1 in rewiring tumor metabolism, driving a redox advantage and protection from cytotoxic therapeutics (Fig. 7F), overall implicating RB as a biomarker for therapy response in advanced disease.
Discussion
Although significant clinical evidence links RB tumor suppressor loss to poor outcome, the underlying mechanisms remain loosely defined. For the first time, the studies herein reveal the disease stage–specific, protumorigenic impact of RB loss on disease progression, supported by the following key findings: (i) RB loss drives distinct biological networks contingent on stage of disease; (ii) rewiring of the E2F1 cistrome and transcriptome in response to RB loss is stage-dependent; (iii) the RB/E2F1 axis shifts regulation of metabolic pathways exclusively in advanced disease, including lipid and amino acid synthesis; (iv) E2F1-mediated networks after RB loss drive increased glutathione synthesis as observed using multiple in vitro and in vivo models of disease, and result in altered ROS production; (v) the relationship between RB and glutathione was confirmed in clinical specimens; and (vi) RB loss–induced ROS reduction facilitates bypass of cytotoxic therapy. Taken together, this study identifies the RB/E2F1 pathway as harboring stage-specific functions that regulate critical metabolic pathways in advanced cancers, which affects therapeutic efficacy.
The data herein reveal striking differences in the biological and molecular outputs governed by RB loss across prostate cancer progression. This evolution of RB regulation was illustrated by transcriptome and E2F1 cistrome analysis, where distinct gene regulation and E2F1 function were observed following RB1 loss in isogenic models of HSPC compared with models of CRPC derived from the same cell lineage(Figs. 1D–F and 2A and B). The study herein is the first to examine the consequence of RB1 depletion under physiologically relevant, AR-active conditions across stages of disease progression, providing the only comparison of the E2F1 cistrome in HSPC and CRPC after RB loss. RB1 depletion in CRPC resulted in 18,624 gained E2F1 binding sites compared with 3,597 gained sites after RB1 depletion in HSPC confirming stage-specific control of E2F1. This expansive shift in E2F1 binding after RB1 depletion in CRPC led to more than 30,000 E2F1 binding sites, largely distinct from those identified in early-stage disease. This expansion of E2F1 binding after RB loss is supported by previous studies in other prostate cancer models (16), implicating RB loss as a direct modulator of E2F1 function. De novo motif analysis revealed that the expanded E2F1 binding is enriched for additional cancer-associated motifs including BORIS, FOXM1, and NF1 motifs (refs. 40–42; Fig. 2C). BORIS and FOXM1 have been shown to be upregulated in prostate cancer, with BORIS positively correlating with Gleason score and AR protein levels and FOXM1 shown to play a role in increasing migratory, invasive, and proliferative abilities (43, 44). Furthermore, NF1 is amplified in neuroendocrine prostate cancer and suggested to regulate AR expression through activation of several kinases including mTOR (45, 46). In breast cancer, alternative splicing of BORIS has been shown to contribute to the Warburg effect, altering cell metabolism and driving tumorigenesis (47), while overexpression of FOXM1 has been shown to play significant roles in regulating cell-cycle progression, invasion, metastasis, and drug resistance (48). This shift in E2F1 binding after RB1 depletion suggests that E2F1 is rewired to novel sites to drive an expanded set of cancer-related biological networks. In addition, via known motif analysis, E2F1 is shown to bind in close proximity to a number of forkhead motifs, which are known coregulators of nuclear receptors, suggesting that these forkhead proteins may also contribute to the recruitment of E2F1 to these novel sites (Fig. 2D). Together, these data are the first to reveal a change in RB function across disease progression driven by altered control of E2F1 and provide molecular insight into the basis of E2F1-associated poor outcomes.
As RB loss is associated with poor outcome, these clinically significant changes in RB and E2F1 function emphasize the importance of understanding the timing of RB loss across disease progression. Previous studies have identified multiple RB-loss gene signatures that are associated with poor outcome including shorter OS, PFS, and DSS (12, 49). Furthermore, it was shown that these signatures can be used to predict RB loss in clinical samples, revealing the clinical significance of RB loss–induced gene expression changes (16). The data herein revealed 1,599 genes with increased expression and E2F1 binding within 30 kb of the gene promoter after RB1 depletion. When compared with previously published RB-loss gene signatures, this list of 1,599 genes included 23% of genes identified in breast cancer and 25% of genes identified in prostate cancer signatures associated with shorter OS (12, 49). Furthermore, when examined beyond breast and prostate cancers, more than 35% of genes identified in multiple NSCLC gene signatures that predicted for poor therapeutic response and shorter OS were also present in this 1,599-gene list (50, 51), indicating that a significant number of novel RB/E2F1–regulated genes in CRPC are likely to have clinical impact. After identifying differential functions of RB and E2F1 across disease progression, the difference from the previously developed gene signatures is unsurprising; however, it highlights the importance of examining RB loss in advanced disease. The RB/E2F1 axis drives networks involved in metabolic control specifically in CPRC, which is not seen in HSPC (Fig. 3A). The shift to metabolic control without the advantage of increased cell growth (Fig. 1C) indicates an evolution away from tumor growth and metastasis toward tumor survival and resistance to therapeutics. As the studies herein reveal that RB loss in late-stage disease render tumors less responsive to therapy, the potential differences in the early- and late-stage signatures are of clinical significance.
Although RB loss in advanced disease induced widespread unique gene networks, a critical consequence of RB loss is the observed metabolic shift in lipid and amino acid synthesis. Analysis of the 1,599 genes with increased expression and E2F1 binding after RB loss revealed 16 enriched metabolic pathways, suggesting a role for E2F1 and RB in regulating metabolism in advanced disease (Fig. 3A). Metabolomic analysis further highlighted a significant change in metabolites involved in lipid and amino acid metabolism (Fig. 4). Beyond this, previous studies have linked RB depletion to altered mitochondrial function. RB loss induces a decrease in mitochondrial protein expression in lung and colon cancers (52). Interestingly, the change in protein expression did not correlate with a corresponding change in mRNA level, suggesting that mRNA expression of mitochondrial genes does not necessarily correlate with protein expression. Although mRNA expression of these mitochondria proteins was also not altered in models herein, further assessment revealed a reduction in oxygen consumption rate (OCR), a measurement of mitochondrial function, in C4–2-shRB cells compared with C4–2-shCon (Supplementary Fig. S8). These data support this previous study, revealing first that mitochondrial function can be altered without changes in mitochondrial gene expression at the mRNA level, and second that RB loss does induce a shift in mitochondrial function in CRPC in addition to lung and colon cancer. Overall, these studies support a role for RB in regulating a multitude of metabolic pathways through mechanisms beyond transcriptional control.
In addition to the studies herein showing that RB and E2F1 function to regulate metabolic control through transcription, E2F1 has also been shown to regulate nucleotide synthesis through transcriptional regulation of thymidine kinase (TK1) and dihydrofolate reductase (DHFR), and has been shown to regulate glucose oxidation through transcription of pyruvate dehydrogenase (PDK1; refs. 53–55). Interestingly, blocking amino acid pathways via inhibition of ASCT2 or LAT1 has been shown to directly shut down E2F transcriptional activation, suggesting a reciprocal feedback between amino acid metabolism and E2F regulation that should be further examined (56, 57). Furthermore, the RB/E2F1 axis has been associated with metabolic control through overexpression of c-Myc. While c-Myc is a known regulator of numerous metabolic transcripts (58), E2F1 occupancy at the c-Myc promoter and mRNA expression of c-Myc is unchanged after RB loss in the models utilized herein. Even with this lack of change in regulation, c-Myc may still have a significant role in metabolism and should be investigated further in the context of RB loss. The RB/E2F1-regulated metabolic genes identified within this study expand on previous RB-loss gene signatures. While genes such as TK1 and RRM2 have already been recognized as altered after RB loss (49), 87 of the 92 increased metabolic genes identified here are exclusive to RB loss in advanced disease. Thus, the studies herein reveal a novel RB-loss/E2F1-regulated metabolic gene signature for advanced disease. This study is the first to explore the genome-wide contribution of E2F1 to metabolic control and highlights that this metabolic control is specific to RB loss in advanced disease.
Further investigation of the RB loss–induced metabolic shift revealed coordinated alteration in glutathione synthesis and downstream ROS generation. Under normal physiologic conditions, cells maintain a redox balance by controlling intracellular ROS with ROS scavengers, such as glutathione. Genes involved in the glutathione synthesis pathway, such as GCLC, CTH, and SLC1A5/ASCT2, and the increase in glutathione observed after RB1 depletion were found to be directly regulated by E2F1 (Figs. 5A and C and 6C–E). This study supports previous data in multiple model systems that RB depletion drives an increase in glutathione. In mouse embryonic fibroblasts (MEF), triple knockout (TKO) of all RB family members (Rb1, Rbl1, and Rbl2) revealed a change in glutamine uptake and increased glutathione production controlled by E2F3. This increase in glutamine uptake and utilization was required for increased proliferation, and modulation of E2F1 had no effect on glutamine uptake or cell proliferation (59). Although the study herein does not examine the effects of RB loss on E2F3 regulation of metabolism, this may be of interest for future studies. Furthermore, increased glutathione was seen in a Drosophila model system with mutant RBF1 through a flux of glutamine, suggesting that this mechanism is conserved across species (60). The elevated glutathione observed herein coincided with lower intracellular basal ROS levels and a lower ROS induction following treatment with cytotoxic agents (Fig. 5D and E). As cancer cells display elevated basal levels of ROS due to the increased proliferation and altered metabolism, this shift in redox homeostasis is often exploited therapeutically through the application of ROS-inducing agents (61, 62). These cytotoxic agents, such as cisplatin and doxorubicin, generate high levels of ROS to induce oxidative stress that leads to damage of protein, lipids, and DNA, promoting apoptosis, necrosis, and autophagic cell death (63). Thus, alterations within this redox control can have significant clinical implications. These studies are the first to confirm that the glutathione increase after RB1 depletion is directly controlled by E2F1 function and that this increase coincides with protection from ROS induction following treatment with clinically approved cytotoxic therapy.
Taken together, the studies herein are the first to demonstrate stage-specific consequences of RB loss in advanced disease. Through transcriptome, cistrome, and metabolome analysis, RB1 depletion rewired E2F1 to drive glutathione synthesis exclusively in advanced disease. This increase in glutathione synthesis drives an acquired protection from ROS generation in response to cytotoxic chemotherapy. These studies reveal a novel role for RB/E2F1 in regulating cell metabolism and implicate RB loss as a stage-specific driver of clinical response to treatment.
Methods
Cell Lines, Cell Culture, and Reagents
LNCaP, C4–2, MCF7, H1299, and UMUC3 cells were obtained from ATCC. LNCaP95 cells were a kind gift from Dr. Jun Luo at John Hopkins University (Baltimore, MD). C4–2-, MCF7-, and LNCaP-derived cell lines were maintained in improved minimum essential medium (IMEM; Thermo Fisher Scientific, 10024CV) supplemented with 5% heat-inactivated FBS, 1% l-glutamine (2 mmol/L), and 1% penicillin–streptomycin (100 U/mL). H1299 and UMUC3 cells were maintained in IMEM supplemented with 10% heat-inactivated FBS, 1% l-glutamine (2 mmol/L), and 1% penicillin–streptomycin (100 U/mL). LNCaP95-derived cell lines were maintained in modified IMEM (Thermo Fisher Scientific, A1048801) supplemented with 10% charcoal dextran-treated (CDT) media and 1% penicillin–streptomycin. All cells were maintained at 37°C with 5% CO2. LNCaP, C4–2, MCF7, H1299, and UMUC3 cell lines were authenticated by ATCC. All cell lines were checked for Mycoplasma upon thawing.
RB1 Knockdown Cell Lines
LNCaP, C4–2, H1299, and MCF7 isogenic models were previously developed as described previously (64, 65). Briefly, cells were transfected with either a shRNA plasmid targeting RB (sequence: 5′-CGCATACTCCGGTTAGGACTGTTATGAA-3′) or a control plasmid (MSCV donor) using FuGENE Transfection Reagent. Stable clones were isolated after 3–4 days of selection with 2.5 μg/mL of puromycin (13). LNCaP95 isogenic models were developed as described previously (65). Cells were transfected with either a plasmid targeting RB (sequence: 5′-GCAGTTCGATATCTACTGAAA-3′) or a nonspecific control (sequence: 5′-GCTGAGGTGATAAACAGTTACA-3′). Newly generated LNCaP95 RB-depleted (miCon and miRB) cell lines underwent at least three rounds of antibiotic selection with 2.5 μg/mL puromycin. These models were labeled LN95-miCon for the control and LN95-miRB for the RB-depleted line.
In Vitro RNA-seq
RNA was extracted and purified utilizing TRIzol reagent and RNeasy Mini Kit (Qiagen) following manufacturer's instructions. RNA-seq libraries were constructed using TruSeq RNA Sample Prep Kit v2 (Illumina, RS-122–2001) and sequenced on the NExtSeq 500 at the Dana-Farber Cancer Institute Molecular Biology Core Facility. RNA-seq was aligned against the hg19 human genome using STAR v2.5.2a (66). Differential gene expression was generated using DESeq2 v1.12.4 (67). GSEA was performed using gene sets from the Molecular Signature Database (MSigDB; ref. 68). RNA-seq data has been deposited in the Gene Expression Omnibus (GEO) under accession number GSE154190.
Proliferation Assays
LNCaP-shCon, LNCaP-shRB, C4–2-shCon, and C4–2-shRB cells were plated at equal densities in 5% CDT media supplemented with 1% penicillin and 1% l-glutamine, to measure relative cell growth in castrate conditions. To measure cell growth in response to androgen stimulation, cells were treated with 0.1 nmol/L DHT after 72 hours, which was replenished every 48 hours. Day 0 timepoints represent cell number 24 hours after plating or 24 hours after DHT treatment. Cell number was quantified at day 4 and day 6 using the Quanti-IT Pico Green dsDNA Assay Kit (Thermo Fisher Scientific).
RB Rescue
PSM.7-LP, an established, constitutively activated RB, was utilized as described previously (37). C4–2 isogenic cells were seeded at equal densities in 5% CDT and transfected using Lipofectamine 2000 (Thermo Fisher Scientific 11668019) with 1 μg of either empty vector (EV) or PSM.7-LP DNA for 24 hours. Cells were treated with 10 nmol/L DHT either 16 hours (mRNA) or 24 hours (protein or glutathione) before harvest. Use of such a strategy is required, as it is well appreciated that the ectopic expression of wild-type RB is rapidly inactivated by endogenous CDKs and does not rescue RB loss (38, 39). As such, shRNA-resistant PSM.7-LP variants were generated and subsequently transfected into the RB knockdown models.
ChIP-seq
E2F1 ChIP-seq was performed as described previously (16). Briefly, cells were cross-linked with 1% fresh formaldehyde for 10 minutes at room temperature. Chromatin was sheared to approximately 200 bp using a Diaganode Ultrasonicator for 30 cycles (30 seconds on, 30 seconds off). E2F1 antibody utilized for ChIP-seq was purchased from Santa Cruz Biotechnology (sc-193). Validation ChIP-qPCR was performed with E2F1 antibody purchased from Bethyl Laboratories (A300–766A), recognizing epitopes distinct from the antibody purchased from Santa Cruz Biotechnology to avoid any accessibility concerns. E2F1 ChIP-seq libraries were constructed using ThruPLEX-FD Prep Kit (Rubicon Genomics). Sequencing was performed on NextSeq 500 at the Dana-Farber Cancer Institute Molecular Biology Core Facility. E2F1 ChIP-seq was deposited in GEO under accession number GSE154191. ChIP-seq was analyzed as described previously (16). Cis-regulatory element analysis was performed using CEAS v1.0.2. Motif analyses were performed through Homer v4.8.3 (69).
Gene Expression
Cells were plated at equal densities in hormone-deficient media for 72 hours followed by 16 hours of 10 nmol/L DHT treatment. RNA was isolated using TRIzol (Invitrogen). RT-PCR was performed using SuperScript VILO cDNA Synthesis Kit (Thermo Fisher Scientific). ABI StepOne Real-Time PCR using PowerSybr (Thermo Fisher Scientific 43–676–59) was utilized to perform quantitative PCR (qPCR) analyses. Primer sequences: GMPR (Fw: 5-′CTCAAGCTCTTCTACGGGA-3′; Rev: 5′-CTCAGAGGCTCTGTACTCAG-3′); PSAT1 (Fw: 5′-AAATGAGTCACAGGTCATCAG-3′; Rev: 5′-CTGGAACAGCTAGCAATTCC-3′);CHPT (Fw: 5′-AGAAGAACAAACTCTTGTTCCC-3′; Rev: 5′-CACTGCCATAAATACTGTGGA-3′); CERK (Fw: 5′-CTCATCCGGAAATGCTCCA-3′; Rev: 5′-AGTGAAGTCAAACTGGTCCT-3′); GSTA1 (Fw: 5′-AAGCTCCACTACTTCAATGC-3′; Rev: 5′-CTCTTCAAACTCTACTCCAGC-3′); CTH (Fw: 5′-ATAGCCGTTCTGGAAATCC-3′; Rev: 5′-AGGCCAAACAGTACTTAGC-3′); GCLC (Fw: 5′-AACCCAAACCATCCTACCC-3′; Rev: 5′-ACATTGTTCCTCCGTAGGG-3′); SLC1A5 (Fw: 5′-TCCTGGATCATGTGGTACG-3′; Rev: 5′-TAAACCCACATCCTCCATCTC-3′); 18S (Fw: 5′-CGGCGACGACCCATTCGAAC-3′; Rev: 5′-GAATCGAACCCTGATTCCCCGTC-3′).
Immunoblotting
LNCaP- and C4–2 derived cells were plated in equal densities in 5% CDT media for 72 hours followed by 24 hours of 10 nmol/L DHT. LN95-derived cells were plated in equal densities in 10% CDT media for 24 hours and treated with 10 nmol/L DHT for 24 hours. MCF7-derived cells were plated in equal densities in 5% FBS media and harvested after 48 hours. UMUC3 cells were transfected with siCtl or siRB as described below in 10% FBS media. Cells were harvested after 72 hours. Cell lysates were generated as described previously (70). Lysate was resolved by SDS-PAGE and transferred to polyvinylidene difluoride membrane. Proteins were analyzed using the following antibodies at 1:1,000 dilution: RB (BD Pharmingen 554136), phospho-RB (Cell Signaling Technology 9307S), AR (N-20, directed against amino acids 1–20 by Bethyl Laboratories), E2F1 (Cell Signaling Technology, 3742S), GAPDH (Santa Cruz Biotechnology, sc-25778), CTH (Santa Cruz Biotechnology, sc-365382), GCLC (Abcam, ab53179), SLC1A5/ASCT2 (Cell Signaling Technology, 8057S), and Vinculin (Sigma-Aldrich, V9264).
Flow Cytometry
LNCaP and C4–2 isogenic cell models were plated in triplicate and treated for 24 hours with 0.5 μmol/L of the CDK4/6 inhibitor palbociclib (Selleck Chemicals S1116). Adherent and nonadherent cells were harvested together after 2 hours of 1:1,000 BrdU labeling (Life Technologies 00–0103) and fixed with 100% ethanol. After additional propidium iodide (PI) staining, flow cytometry was performed using a Guava easyCyte flow cytometer. At least 10,000 events were analyzed per sample. Analysis was performed using the Guava InCyte software.
ChIP-qPCR
qPCR analyses were performed on the ABI StepOne Real-Time PCR using PowerUPsybr (Fisher Scientific, A25742). Primer sequences are listed below. For ChIP-qPCR, C4–2-shCon and C4–2-shRB cells were transfected as described below. Cells were transferred to 5% CDT media following transfection for 72 hours and treated with 10 nmol/L DHT for 3 hours. The following primer sequences were used: GMPR (Fw: 5′-TGTCCCTAACCAGTCTGCCT-3′; Rev: 5′-ACGCCGCGGTCTTAAACA-3′); PSAT1 (Fw: 5′-CACCTTCTTCTGGTTTGGGCT-3′; Rev: 5′-CGAGTGAGCTGGAAACCTGC-3′); CHPT (Fw: 5′-TGTTCCGACCCCTAGGAGAG-3′; Rev: 5′-CGCCAACGGCAACTTGTTAT-3′); CERK (Fw: 5′-CAGCCAGTCTGGAAGAGCAG-3′; Rev: 5′-AGCACTGTGCAAGCCAAATG-3′); GSTA1 (Fw: 5′-CTCAACGGAGTCCCACAGAC-3′; Rev: 5′-TTCCAGGAGAACACCCAAGC-3′); CTH (Fw: 5′-GGCTTACGTCTTTTAGCACCG-3′; Rev: 5′-GGTCCAAGTGGTGTTAGTCCC-3′); GCLC (Fw: 5′- AGGCCGATTCCCAGTAGTGT -3′; Rev: 5′-GGGGGCTGTATCACAACAAC -3′); SLC1A5 (Fw: 5′-CTGCCAGGGATAAGGGTAGG-3′; Rev: 5′-AACCCCTGTGGTTTAAGGGC -3′); Desert (Fw: 5′-CTAGGAGGGTGGAGGTAGGG-3′; Rev: 5′-GCCCCAAACAGGAGTAATGA-3′).
Metabolomics
Metabolomics were performed by Metabolon, Inc. C4–2-shCon and C4–2-shRB cells were seeded at 1.5 million cells per plate in 5% CDT media for 72 hours followed by 10 nmol/L DHT treatment for 16 hours in quadruplicate. Cell pellets were collected following Metabolon sample guidelines, flash-frozen, and stored at −80°C until shipment. Metabolite identification was performed through ultra-high-performance liquid-phase chromatography and gas chromatography separation, coupled with mass spectrometry. Compounds were identified by comparison to libraries maintained by Metabolon. Data visualization was performed using Metabolon's custom-built Cytoscape plugin.
Amino Acid Uptake Assay
C4–2–derived cells were assessed using a [3H]-labeled amino acid uptake assay, as described previously (71, 72). Briefly, cells (1 × 105/well) were incubated at 37°C in 96-well plates with [3H]-l-glutamine (400 nmol/L; Perkin Elmer) or [3H]-l-leucine (200 nmol/L; Perkin Elmer) in glutamine-free RPMI or leucine-free RPMI media (Invitrogen) for 15 minutes in the presence or absence of 2 mmol/L GPNA (MP Biochemicals; diluted in H2O). Cells were collected and transferred to filter paper using a 96-well plate harvester (Wallac, PerkinElmer), dried, exposed to scintillation fluid (PerkinElmer), and [3H] levels analyzed using a liquid scintillation counter (PerkinElmer).
Glutathione Assay
C4–2-, LN95-, MCF7-, and H1299-derived cells were seeded at equal densities in 5% CDT media for 72 hours followed by 24 hours 10 nmol/L DHT treatment. Relative glutathione was measured using a Glutathione Assay Kit (Sigma-Aldrich, CS0260) following manufacturer's instructions. UMUC3 cells were transfected with siCtl or siRB as described below and harvested after 72 hours in 10% FBS IMEM. For assays with siE2F1, transfection was performed as described below in 5% CDT media for 72 hours following transfection and treated with 10 nmol/L DHT for 24 hours. For assays with drug treatments, C42 isogenic cells were seeded at equal densities in 5% CDT for 72 hours. Cells were treated with 10 nmol/L DHT for 24 hours before an additional 24-hour treatment with either 0.5 mmol/L BSO (Sigma-Aldrich B2515) or 1 μmol/L Erastin (Sigma-Aldrich E7781) prior to harvest. RB rescue with PSM.7-LP assays was performed after 72-hour transfection as described above in 5% CDT and 24 hours of 10 nmol/L DHT treatment. Relative fluorescence was measured on a BioTek Synergy HT microplate reader and analysis using BioTek Gen5 2.09 software.
ROS Assay
C4–2-, LN95-, and MCF7-derived cell lines were seeded at equal densities. C4–2 cells were plated in 5% CDT media, LN95 cells were plated in 10% CDT media, and MCF7 cells were plated in 5% FBS media. C4–2 and LN95 isogenic cells were treated with 10 nmol/L DHT for 24 hours following 72 hours of CDT. Menadione (Sigma-Aldrich, M5625) treatment was utilized as a positive control. Relative ROS levels were measured by fluorescence using the ROS-Glo Assay (Promega, G8821). Fluorescence was measured on a BioTek Synergy HT microplate reader and analysis was done using BioTek Gen5 2.09 software.
RNA Interference
C4–2-shCon and C4–2-shRB cell lines were seeded on poly-l-lysine–coated plates in 5% FBS media for 24 hours. Cells were transfected with either scramble control (C4–2-shCon siCtl, C4–2-shRB siCtl) or E2F1 (C4–2-shCon siE2F1, C4–2-shRB siE2F1) siRNA pools (Dharmacon) for 8 hours according to the manufacturer's protocol. After transfection, the cells were moved to 5% CDT media for 72 hours followed by appropriate 10 nmol/L DHT treatment. UMUC3 cells were seeded on poly-l-lysine–coated plates in 10% FBS media for 24 hours. Cells were transfected with either scramble control (UMUC3 siCtl) or RB (UMUC3 siRB) siRNA pools (Dharmacon) for 8 hours according to the manufacturer's protocol. Cells were harvested after 72 hours.
MTT Cell Viability Assay
C4–2 and MCF7 isogenic pairs were seeded at equal densities in a 96-well plate. C4–2–derived models were seeded in 5% CDT for 72 hours followed by 24 hours of 10 nmol/L DHT. MCF7-derived models were seeded in 5% FBS for 48 hours. Cells were treated with varying concentrations of Erastin (0.1 μmol/L to 500 μmol/L) or BSO (0.01 mmol/L to 50 mmol/L) for 24 hours. MTT assay was performed according to the manufacturer's instructions (Sigma-Aldrich 11465007001). Absorbance was measured at 625 nm on a BioTek Synergy HT microplate reader and analysis using BioTek Gen5 2.09 software.
Xenografts
C4–2-shCon and C4–2-shRB cells were resuspended in 100 μL of saline with 50% Matrigel (BD Biosciences 354234). Cell mixture (3 × 106 cells per injection) was injected subcutaneously into the flank of 5- to 6-week-old athymic nude male mice following protocols approved by the Institutional Animal Care and Use Committee at Thomas Jefferson University (Philadelphia, PA). Tumor volume was measured with digital calipers multiple times per week. Mice were sacrificed, and tumors were harvested once tumors reached an approximate size of 500 mm3. Tumors were cut into four pieces for glutathione, mRNA, and IHC. Glutathione analysis was performed following manufacturer's guidelines using a Glutathione Assay Kit (Sigma-Aldrich, CS0260). For RNA extraction, tumor tissue was homogenized, and RNA was extracted using TRIzol (Invitrogen). IHC staining was performed by the Sidney Kimmel Cancer Center Shared Resource Translational Pathology Core. Antibodies utilized included RB (Abcam ab226979), GCLC (Abcam ab 53179), and SLC1A5/ASCT2 (Cell Signaling Technology 8057S).
Royal Marsden Hospital Cohort
RMH cohort patient samples (n = 98) were collected and analyzed with written informed consent under the CCR2472 protocol approved by the Ethics Committee at the Royal Marsden NHS Foundation Trust Hospital (London, United Kingdom).
RMH Sample RNA-seq
RNA, from fresh tissue biopsy, quality was analyzed using Agilent RNA ScreenTape Assay (Agilent). Five hundred nanograms of total RNA was used for library preparation using the NEBNext rRNA Depletion Kit followed by NEBNext Ulta II directional RNA assay kit as per manufacturer's protocol (New England Biolabs). Library quality was confirmed using the Agilent High sensitivity D1000 ScreenTape Assay (Agilent). The libraries were quantified and normalized by qPCR using Generead Library Quant Kit (Qiagen). Library clustering and sequencing were performed on the Illumina NovaSeq 6000. The libraries were run across two lanes of an Illumina NovaSeq S2 flowcell using 150 base pair pair-end v1 kit and eight base pair dual indexes. Base calling and quality scoring were performed using Real-Time Analysis (version v3.4.4) and FASTQ file generation and demultiplexing using Illumina bcl2fastq2 (version 2.20). CRPC transcriptome reads were aligned to the human reference genome (GRCh37/hg19) using TopHat2 (version 2.0.7). Gene expression, fragments per kilobase of transcript per million mapped reads (FPKM), was calculated using Cufflinks. AR activity score was an accumulation measurement of AR pathway activity based on 43 genes regulated by AR in prostate cancer cell lines and metastatic prostate cancer as described previously (73).
Seahorse Assay
Cellular oxygen consumption rate (OCR) was measured with a Seahorse XF96e bioanalyzer (Seahorse Bioscience) using a Cell Mito Stress Test Kit according to manufacturer's instructions. Briefly, C4–2-shCon or C4–2-shRB cells (1.5 × 104 cells/well) were seeded in 4 wells (technical replicates) each of a Seahorse XF 96-well assay plate in full growth medium. After overnight attachment, the medium was carefully washed and replaced with prewarmed running medium (consisting of nonbuffered DMEM, Agilent Technologies) supplemented with 1 mmol/L sodium pyruvate, 2 mmol/L glutamine, and 10 mmol/L glucose, pH 7.4). Plates were incubated for 60 minutes in a non-CO2 incubator at 37°C before four basal measurements were undertaken, determining oxygen and proton concentration in the medium. The ATP synthase inhibitor oligomycin (1 μmol/L) was immediately added followed by four further measurements; carbonyl cyanide-4-(trifluoromethoxy) phenylhydrazone (FCCP; 1 μmol/L) was then added with another four further measurements, before addition of the complex I and III inhibitors rotenone/antimycin A (0.5 μmol/L) and a final four further measurements taken.
Statistical Analysis
All experiments were performed in technical triplicate with at least three biological replicates per condition. Data are displayed as mean ± SEM. Statistical significance (P < 0.05) was determined using Student t test, one-way ANOVA, and two-way ANOVA on GraphPad Prism Software as appropriate.
Authors' Disclosures
W. Yuan reports other support from Jilin Huarui Gene Technology Ltd outside the submitted work. J.S. de Bono reports personal fees and other support from Amgen, Bioxcel Therapeutics, Boehringer Ingelheim, Eisai, Menarini Silicon Biosystems, Qiagen, and Terumo; grants, personal fees, and other support from Astellas, AstraZeneca, Bayer, Cellcentric, Daiichi, Genentech Roche, Genmab, GlaxoSmithKline, Harpoon, Janssen, Merck Serono, Merck Sharp & Dohme, Orion Pharma, Pfizer, Sanofi Aventis, Sierra Oncology, and Taiho, and grants, personal fees, and other support from Vertex Pharmaceuticals outside the submitted work; in addition, J.S. de Bono has a patent for DNA damage repair inhibitors for treatment of cancer (patent no. WO 2005 053662 - no personal income) issued and licensed to AstraZeneca and a patent for 17-substituted steroids useful in cancer treatment (patent no. US5604213 - no personal income) issued and licensed to Janssen. K.E. Knudsen reports grants from NCI during the conduct of the study; other support from Janssen and Genentech, and other support from CellCentric outside the submitted work. No disclosures were reported by the other authors.
Authors' Contributions
A.C. Mandigo: Conceptualization, data curation, formal analysis, funding acquisition, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. W. Yuan: Formal analysis. K. Xu: Investigation. P. Gallagher: Investigation. A. Pang: Investigation. Y. Guan: Investigation. A.A. Shafi: Writing–review and editing. C. Thangavel: Investigation. B. Sheehan: Formal analysis. D. Bogdan: Formal analysis. A. Paschalis: Formal analysis. J.J. McCann: Funding acquisition, writing–review and editing. T.S. Laufer: Investigation. N. Gordon: Writing–review and editing. I.A. Vasilevskaya: Writing–review and editing. E. Dylgjeri: Writing–review and editing. S.N. Chand: Writing–review and editing. M.J. Schiewer: Conceptualization, writing–review and editing. J. Domingo-Domenech: Writing–review and editing. R.B. Den: Resources. J. Holst: Conceptualization, resources, writing–review and editing. P.A. McCue: Formal analysis. J.S. de Bono: Resources. C. McNair: Software, formal analysis, methodology, writing–original draft, writing–review and editing. K.E. Knudsen: Conceptualization, resources, supervision, funding acquisition, writing–original draft, project administration, writing–review and editing.
Acknowledgments
This work was supported by NIH T32 grant (GM100836; to A.C. Mandigo); NIH R01 grants (5R01CA17640105, 5R01CA18256905, 5R01CA21732903; to K.E. Knudsen); NCI F99 grant (F99CA212225; to J. McCann); and the Sidney Kimmel Cancer Center (SKCC) Support Grant (5P30CA056036). Additional support was provided by the SKCC Cancer Genomics, Translational Research/Pathology and Biostatistics core services, specifically Dr. Benjamin Leiby. We would like to thank Dr. Jun Lou (John Hopkins University, Baltimore, MD) for providing the LNCaP95 cells utilized in this study. The sponsors were critical in study design, data collection, analysis, interpretation, and review of the manuscript.