The TMPRSS2:ERG gene fusion is common in androgen receptor (AR) positive prostate cancers, yet its function remains poorly understood. From a screen for functionally relevant ERG interactors, we identify the arginine methyltransferase PRMT5. ERG recruits PRMT5 to AR-target genes, where PRMT5 methylates AR on arginine 761. This attenuates AR recruitment and transcription of genes expressed in differentiated prostate epithelium. The AR-inhibitory function of PRMT5 is restricted to TMPRSS2:ERG-positive prostate cancer cells. Mutation of this methylation site on AR results in a transcriptionally hyperactive AR, suggesting that the proliferative effects of ERG and PRMT5 are mediated through attenuating AR’s ability to induce genes normally involved in lineage differentiation. This provides a rationale for targeting PRMT5 in TMPRSS2:ERG positive prostate cancers. Moreover, methylation of AR at arginine 761 highlights a mechanism for how the ERG oncogene may coax AR towards inducing proliferation versus differentiation.https://doi.org/10.7554/eLife.13964.001
Prostate cancers are among the most common types of cancer in men, which, like other cancers, are driven by genetic mutations. Roughly half of all prostate cancers contain a genetic change that incorrectly fuses two genes together, causing the cells to produce abnormally high levels of a protein called ERG.
ERG is a transcription factor, a protein that binds to specific sequences of DNA to influence the activity of nearby genes. ERG represses genes that help to prevent prostate cancers from growing, and so promotes prostate cancer development. Like most other transcription factors, ERG is difficult to target with drugs and no therapies that directly prevent the activity of ERG currently exist.
Mounir et al. wanted to find out whether ERG cooperates with other proteins to cause prostate cancer cells to grow, with the hope that these proteins could be more easily targeted with a drug. By using various biochemical techniques in human prostate cancer cell lines, Mounir et al. found that ERG interacts with an enzyme called PRMT5. This interaction enables PRMT5 to chemically modify other proteins to change their activity. In the case of prostate cancer cells, PRMT5 inappropriately modifies the androgen receptor, a protein that regulates the growth of normal prostate cells. This abnormal modification contributes to the excessive growth of the cancer cells.
Although PRMT5 will be easier to target with drugs than ERG, it also has many other roles besides those described by Mounir et al. Much more work is therefore needed to investigate whether PRMT5 could be safely targeted to treat patients with prostate cancer.https://doi.org/10.7554/eLife.13964.002
Prostate cancer (PC) is highly prevalent and lethal (Siegel et al., 2015). Drugs targeting the Androgen Receptor (AR), a 'lineage driver' of PC (Garraway and Sellers, 2006), are an important therapeutic approach. AR is an androgen (i.e. testosterone)-activated nuclear hormone receptor that regulates normal prostate gland growth and differentiation. In PC however, AR facilitates unregulated proliferation (Mills, 2014). While it is unclear how AR and other lineage factors switch between promoting normal lineage differentiation vs. tumor growth, it is hypothesized that somatic mutations in additional genes may facilitate such changes (Garraway and Sellers, 2006). Many PCs bear chromosomal translocations resulting in aberrant expression of the ETS transcription factor ERG, most commonly through the TMPRSS2:ERG fusion (Shah and Chinnaiyan, 2009). TMPRSS2:ERG, alone or in combination with additional genetic alterations, promotes prostate tumor formation in mice (Baena et al., 2013; Chen et al., 2013; King et al., 2009; Klezovitch et al., 2008; Mounir et al., 2015; Tomlins et al., 2008). ERG is recruited to many AR target genes and represses AR-dependent transcription (Yu et al., 2010), suggesting ERG functions at least in part through attenuating AR target gene expression. However, ERG also regulates the expression of AR-independent genes thought to drive oncogenic function (Klezovitch et al., 2008; Mounir et al., 2015; Tomlins et al., 2008; Wang et al., 2008). We explored whether a deeper mechanistic understanding of ERG proliferative function could yield therapeutic insights into targeting this key PC oncogene.
To identify genes that selectively facilitate the growth of TMPRSS2:ERG positive PC cells, we performed a pooled short hairpin RNA (shRNA) screen in TMPRSS2:ERG and AR-positive VCaP prostate cancer cells, using ERG-negative 22Rv1 cells as a control (Materials and methods). The shRNA pool targets 648 genes involved in transcriptional and epigenetic regulation (Supplementary file 1). While ERG shRNAs were not in the pool, AR shRNAs were preferentially depleted from VCaP cells, underscoring AR dependence in this cell line. Thirty two (32) genes showing VCaP-selective shRNA depletion (Materials and methods) were considered for further study (Figure 1A; Supplementary file 1).
We next narrowed the shRNA screen hit list by focusing on candidates more likely to be ERG interacting proteins. We immunoprecipitated ERG from VCaP cells, and then identified co-immunoprecipitated proteins by mass spectrometry (Materials and methods). Identified proteins (Supplementary file 2) included AR and DNA-PKcs, previously known ERG interactors (Brenner et al., 2011; Yu et al., 2010). Eight of the VCaP-selective shRNA screen hits that also co-immunoprecipitated with ERG were further validated by directed ERG co-immunoprecipitation experiments in VCaP cells. Of these, AR and PRMT5 were the only proteins that co-immunoprecipitated with ERG but not IgG control; these interactions were not overtly influenced by exposure to an androgen analog (R1881, Figure 1B; Figure 1—figure supplement 1A). We next tested whether the ERG/PRMT5 interaction is observed in other models. PRMT5 co-immunoprecipitated with ERG in 22Rv1 cells ectopically expressing ERG. This interaction was still observed upon expression of ERG bearing mutations in the DNA binding domain ('Dx', Figure 1C), suggesting DNA binding is not required for the ERG/PRMT5 interaction. Reciprocal co-immunoprecipitation experiments using overexpressed ERG and PRMT5 in AR-negative 293 and PC3 cells suggest the ERG/PRMT5 interaction can occur in the absence of AR (Figure 1—figure supplement 1B). Further work in 293 cells using truncated ERG constructs suggested that the conserved ETS DNA binding domain of ERG was necessary for the observed co-immunoprecipitation with PRMT5 (Figure 1—figure supplement 1C). Given this evidence that ERG and PRMT5 co-exist in a protein complex, we focused further efforts on PRMT5, as to our knowledge it has not been previously linked to ERG biology.
To validate the growth effects of PRMT5 knockdown, we transduced ERG-positive VCaP cells, and ERG-negative 22Rv1 and LNCaP PC cells, with three independent doxycycline (Dox)-inducible shRNA vectors targeting PRMT5 and a non-targeting control shRNA (NTC). PRMT5 knockdown was robust in all cell lines (Figure 1—figure supplement 1D). Robust growth inhibition was observed in VCaP cells; in contrast PRMT5 knockdown had no growth effects in ERG-negative 22Rv1 cells, and only minor effects in ERG negative LNCaP cells (Figure 1D–F). Deletion of methylthioadenosine phosphorylase (MTAP), which is common across cancers, is a major determinant of sensitivity to PRMT5 inhibition (Kryukov et al., 2016; Mavrakis et al., 2016); as VCaP, LNCaP, and 22Rv1 cells are all MTAP intact, the observed sensitivity of VCaP to PRMT5 shRNA is not due to MTAP deletion. The project Achilles shRNA screen dataset (Kryukov et al., 2016) contains three prostate cancer cell lines (VCaP, 22Rv1 and TMPRSS2:ERG positive NCI-H660) and one PRMT5 hairpin likely to have minimal off-target effects. This shRNA shows a trend of sensitivity in ERG-positive lines, in agreement with our findings (Figure 1—figure supplement 1E).
PRMT5 is a protein arginine methyltransferase that regulates multiple signaling pathways through the mono- and symmetric di-methylation of arginines on its target proteins (Yang and Bedford, 2013).To determine whether the antiproliferative effects of PRMT5 knockdown in ERG positive VCaP cells were mediated through methyltransferase activity, we expressed shRNA-resistant wild-type PRMT5, or a catalytically inactive G365A/R368A double mutant (Materials and methods) (Antonysamy et al., 2012) along with PRMT5 shRNA in VCaP cells. WT PRMT5, but not the G365A/R368A mutant, rescued the effects of PRMT5 knockdown on VCaP cell proliferation (Figure 1G), indicating a requirement for PRMT5 catalytic function to support VCaP proliferation.
To understand pathways affected by PRMT5, we performed transcriptional profiling of PRMT5 knockdown in VCaP cells, followed by the identification of significantly altered pathways (Figure 2A; Supplementary file 3; see Materials and methods). Among these, AR activation was the second most significantly affected pathway, and is a key pathway in common with previous reports of ERG knockdown in VCaP cells (Chen et al., 2013; Mounir et al., 2015; Yu et al., 2010). AR pathway upregulation was apparent using multiple published AR gene signatures (Figure 2—figure supplement 1). Using quantitative PCR of reverse transcribed RNA (qRT-PCR), we confirmed that knockdown of either PRMT5 or ERG increased the expression of the AR target genes PSA, NKX3-1 and SLC45A3 (Figure 2B). Expression of shRNA-resistant WT PRMT5, but not the G365A/R368A mutant, rescued the effects of PRMT5 knockdown on AR target gene expression (Figure 2C; Figure 2—figure supplement 2A), demonstrating that PRMT5 methyltransferase activity is required for repression of AR target genes. The effect of PRMT5 knockdown was restricted to genes co-regulated by both AR and ERG, as PRMT5 knockdown did not affect previously published (Mounir et al., 2015) AR-independent ERG target genes in VCaP (Figure 2—figure supplement 2B). In addition, PRMT5 knockdown did not induce AR target gene expression in ERG negative 22Rv1 or LNCaP PC cells (Figure 2—figure supplement 2C). In 22Rv1 cells, exogenous ERG expression is sufficient to attenuate PSA, NKX3-1, and SLC45A3 expression (Mounir et al., 2015). Under these conditions, PRMT5 knockdown restored the expression of these genes to baseline levels, yet had no effect on their expression in the background of the DNA-binding defective (i.e. inactive) ERG mutant ('Dx', Figure 2—figure supplement 2D). These data indicate that PRMT5’s ability to repress AR function is dependent on ERG.
We next used chromatin immunoprecipitation (ChIP) to investigate ERG, AR, and PRMT5 recruitment to the AR targets PSA and NKX3-1 following modulation of ERG or PRMT5 expression. ERG knockdown in VCaP reduced its recruitment to previously characterized (Wei et al., 2010) binding sites on both genes (Figure 2D). We also observed PRMT5 recruitment to these same sites, which was dramatically reduced upon ERG knockdown. Conversely, ERG expression in 22Rv1 cells induced PRMT5 recruitment to these same sites (Figure 2—figure supplement 3A). In VCaP, PRMT5 knockdown reduced its own recruitment to both genes but had virtually no effect on ERG recruitment (Figure 2D). As expected (Yu et al., 2010), ERG expression in 22Rv1 reduced AR recruitment to both genes (Figure 2—figure supplement 3A), and ERG knockdown in VCaP increased AR recruitment (Figure 2D). Like ERG, PRMT5 knockdown in VCaP strongly induced AR recruitment (Figure 2D).
To extend these findings to a genome wide scale, we performed AR and ERG ChIP-seq experiments in androgen (DHT or R1881) stimulated VCaP cells upon PRMT5 knockdown (Materials and methods; available PRMT5 antibodies did not work in our hands for ChIP-seq). ERG and AR recruitment were robust in these experiments, as judged by recovery of the canonical DNA binding sites of these proteins (Figure 2—figure supplement 3B–D). In agreement with directed ChIP experiments, AR ChIP-seq demonstrated that PRMT5 knockdown increased the recruitment of AR at a subset of peaks (6% p-val 0.01, q-val 0.151), (Figure 2E), but did not significantly affect ERG binding on the same set of sites. Collectively, these data support a model where ERG recruits PRMT5 to AR targets, and PRMT5 is required for ERG-dependent attenuation of AR binding to specific regulatory regions.
In many cell types, PRMT5 represses gene expression through symmetric di-methylation of histone H4 at arginine 3 (H4R3me2s) (Yang and Bedford, 2013), suggesting this as a mechanism for PRMT5 function at AR target genes. However, ERG expression in 22Rv1 cells did not alter H4R3me2s levels at the PSA or NKX3-1 loci (Figure 3—figure supplement 1). We therefore hypothesized that PRMT5 directly methylates AR. To test this, we immunoprecipitated AR from ERG-positive VCaP cells, and blotted for AR, mono-methyl arginine (MMA), or symmetric di-methyl arginine (SDMA). MMA and SDMA signals were indeed detected in AR immunoprecipitates from VCaP cells, and were reduced following knockdown of either ERG or PRMT5 (Figure 3A). To confirm these findings in an additional ERG-dependent model, we also immunoprecipitated AR from the prostates of wild-type and TMPRSS2:ERG transgenic mice, the latter of which show ERG-dependent hyperproliferative phenotypes (Mounir et al., 2015). Like VCaP, immunoprecipitated AR from TMPRSS2:ERG mouse prostates showed high mono- (MMA) and symmetric di-methylation (SDMA) levels compared to wild-type controls (Figure 3B). On the other hand, LNCaP cells, which bear a translocation in the ETS factor ETV1, do not show any SDMA signal on AR, and show reduced levels of MMA versus VCaP (Figure 3C), suggesting that AR symmetric dimethylation is unique to ERG versus ETV1. We next immunoprecipitated AR from 22Rv1 cells, which express wild-type AR as well as a roughly 80 kDa truncated AR variant that lacks its ligand binding domain (LBD) (Dehm et al., 2008). In these cells, ERG expression increased MMA and SDMA signals on wild-type AR but not the truncated variant, and these ERG-dependent signals were reduced upon PRMT5 knockdown (Figure 3D), suggesting the AR LBD is mono- and symmetrically di-methylated in an ERG- and PRMT5-dependent manner.
To further understand if the observed AR methylation was directly dependent on PRMT5, we performed biochemical assays using purified PRMT5 (complexed with its requisite binding partner MEP50/WDR77) and AR LBD. PRMT5 activity, as judged by production of SAH (the by-product of SAM-dependent substrate methylation), was observed in the presence of AR LBD as substrate, but not in the presence of ERG ETS DNA binding domain or pointed (PNT) domain. PRMT5 activity in the presence of AR LBD further increased with the addition of ERG ETS domain protein to the reaction, but not with PNT domain (Figure 3—figure supplement 2A). Unlike PRMT5, purified PRMT1 showed no activity in the presence of AR LBD (Figure 3—figure supplement 2B). PRMT5 activity on AR, in the absence or presence of ERG ETS domain, was reduced with the addition of the tool PRMT inhibitor AMI-1 (Cheng et al., 2004). PRMT5 dependent methylation of the AR-LBD was also observed in western blots, where increased signal for MMA (but not SDMA, likely due to low overall PRMT5 activity in vitro and the distributive nature of this enzyme [Wang et al., 2014]) was observed at the correct size for the AR-LBD protein in the presence of PRMT5. This activity was increased by ERG ETS domain, and inhibited by AMI-1 (Figure 3—figure supplement 2C–D). Together this data indicates that AR LBD is a substrate of the PRMT5 enzyme in vitro, and that the ERG ETS domain can facilitate greater PRMT5 activity on AR.
To identify the arginine methylation site(s) on AR, we cloned AR cDNA and mutated all arginines in the LBD to lysine (see Materials and methods). We expressed each construct in AR-negative, ERG-negative RWPE-1 prostate cells. R1881 stimulates exogenous AR to induce target gene expression in these cells. ERG expression represses this effect (Figure 3—figure supplement 3A–B) and induces PRMT5-dependent mono- and symmetric di-methylation of AR (Figure 3E). RWPE-1 cells co-expressing ERG with each of the AR constructs were assessed for SDMA and MMA modification of AR following immunoprecipitation. All mutants were expressed at equivalent levels to wild-type AR (Figure 3—figure supplement 3C). Only one AR mutant, R761K, completely lacked MMA and SDMA of AR in the presence of ERG (Figure 3F; Figure 3—figure supplement 3D). We confirmed these results in RWPE-1 cells using a proximity ligation assay to detect symmetric di-methylation of AR (Materials and methods). We detected strong proximity signals in RWPE-1 nuclei upon expression of wild-type AR and exposure to R1881 (which stimulates nuclear translocation of AR) that were dependent upon ERG expression. These signals were not detected in the AR R761K mutant cell line (Figure 3G). This demonstrates that arginine 761 is the likely target of ERG- and PRMT5-dependent AR methylation.
Available structural information (Chandra V et al., 2008; Emsley et al., 2010; Helsen et al., 2012) (see Materials and methods) suggests R761 methylation on AR may affect LBD interactions with the DNA binding domain (DBD, Figure 4A), resulting in altered DNA binding and gene activation. To test this, we evaluated the recruitment of WT versus R761K mutant AR in RWPE-1 cells treated with R1881, using the PSA locus as a model AR-regulated gene. Relative to WT AR, the AR R761K mutant showed enhanced recruitment to PSA and increased PSA expression. Moreover, R761K mutation prevented ERG-dependent attenuation of AR recruitment to PSA and PSA expression (Figure 4B,C). In contrast to these effects, the AR R761K mutation did not affect overall AR protein levels, the ability of AR to shuttle into the nucleus upon R1881 induction (Figure 4—figure supplement 1), or the levels of ERG or PRMT5 recruitment to PSA (Figure 4—figure supplement 1C). These results indicate R761 methylation mediated through ERG and PRMT5 attenuates ligand-dependent AR activation, likely through modulating interactions between the LBD and DBD.
Increased AR activity drives PC proliferation (Chen et al., 2004). However, we noted that RWPE-1 cells expressing AR R761K, despite increased AR activity, were more prone to R1881-induced growth arrest (Figure 4D). This suggested that heightened AR activity via R761K mutation (or loss of R761 methylation) may block proliferation in certain contexts. To further explore this, we expressed R761K mutant AR in VCaP cells. Despite transcriptional hyperactivity, as judged by increased PSA expression relative to either WT or an irrelevant R789K mutant AR, VCaP cells expressing AR R761K proliferated poorly in reduced-androgen media (Figure 4E and F).
Detailed mechanisms of how the tumor context affects AR and other lineage oncogenes to switch their function from lineage differentiation to proliferation in cancer has remained elusive (Garraway and Sellers, 2006). Our results indicate that a primary effect of ERG in facilitating PC proliferation is modulating AR function. We propose that ERG- and PRMT5-dependent methylation of R761 on AR reduces AR recruitment to genes that would otherwise induce differentiation, yet R761 methylation allows sufficient AR function to promote proliferation. While our results raise the question of how R761-methylated AR may still facilitate proliferation, R761 methylation could prove a relevant biomarker for AR-mediated proliferation versus arrest in TMPRSS2:ERG positive cells. Finally, as PRMT5 is an enzyme required downstream of ERG in facilitating AR proliferative function, exploring therapeutic PRMT5 inhibition in TMPRSS2:ERG positive prostate cancers may be warranted.
VCaP, 22Rv1, LNCaP, PC3, RWPE-1, and 293T cells were obtained from ATCC and maintained in recommended media unless otherwise specified. Cell identities were verified by SNP analyses using ABI TaqMan SNP genotyping assays (Asuragen) and tested for mycoplasma contamination using the MycoAlert Mycoplasma Detection kit (Lonza).
22Rv1 and VCaP cells were used for a screen using a custom shRNA library (Cellecta, Inc.) targeting transcriptional and epigenetic regulators similar to a previously reported library (Hoffman et al., 2014). shRNA library design and construction and viral packing were performed as previously described (Hoffman et al., 2014). To obtain an MOI of 0.3, the required volume of virus was determined using a 10 point dose curve ranging from 0 to 1ml of viral supernatant in the presence of 10 ug/ml polybrene. Infection efficiency was determined by the percentage of RFP positive cells measured by FACS analysis. Screens were run in duplicate. For the 22Rv1 screen, 14.4 million cells were plated 24 hr prior to infection in T-225 flasks. On the day of infection, the culture media was replaced with fresh media containing 10 ug/mL polybrene and sufficient virus was added for an MOI of 0.3. 24 hr after infection, the culture media was replaced with fresh media containing puromycin. 72 hr following puromycin addition, cells were trypsinized, and 14.4 million cells were plated into new flasks. For the VCaP screen, 23.7 million cells were plated 24 hr prior to infection in T-225 flasks. On the day of infection, the culture media was replaced with fresh media containing 10 ug/mL polybrene and sufficient virus was added for an MOI of 0.3. 24 hr after infection, the culture media was replaced with fresh media and cells were allowed to recover for 72 hr prior to puromycin selection. 5 days following puromycin addition, cells were trypsinized, and 23.7 million cells were plated into new flasks. At each passage, an aliquot of cells was used to measure transduction efficiency determined by measuring the% RFP positive cells and was typically > 90%. Cells were maintained in culture and split when confluence reached 90% and at each passage, 14.4 million cells (22Rv1 screen) and 23.7 million cells (VCaP screen) were passaged into new flasks, ensuring a representation of >1000 cells/shRNA in the library and the% RFP positive cells was measured to ensure stability of the transduced population over time. When the cells reached 5-population doublings, 40 million cells were harvested by centrifugation and stored at −20°C. Purification of genomic DNA and PCR for library production were performed as previously described (Hoffman et al., 2014). shRNA screen data analysis was performed as previously described (Hoffman et al., 2014). For gene based hit calling, the The Redundant siRNA Activity or RSA metric was used as described (Hoffman et al., 2014). Briefly, the RSA down p-value is the statistical score that models the probability of a gene 'hit' based on the collective activities of multiple shRNAs per gene. The RSA down p-value reports the statistically significant genes causing a loss in viability. All hits showing an RSA score >10-5 in the 22Rv1 screen and <10-5 in the VCaP screen (total of 32 genes) were used for further analysis (see Supplementary file 1 for full list).
Nuclear extracts of VCaP cells (~600 million cells) were pre-cleared using agarose beads (Trueblot rabbit kit, eBioscience) and used for pulldown with either an ERG antibody (Santa Cruz#353) or an anti-rabbit IgG antibody. All antibodies were pre-coupled to beads (AminoLink Plus kit, Thermo Scientific), washed then used for pulldown. Immunoprecipitations were eluted at pH2.5 followed by TCA protein precipitation, alkylated with iodoacetamide, and separated on a NuPage 4–12% Bis-Tris gradient gel (Invitrogen). Complete gel lanes were excised using a LEAP 2DiD robot and in-gel digested with trypsin (Tecan Freedom EVO 20). Peptide sequencing for the resulting resulting 16 digest samples was performed by liquid chromatography-tandem mass spectrometry using an Eksigent 1D+ high-pressure liquid chromatography system coupled to a LTQ-Orbitrap XL mass spectrometer (Thermo Scientific). Peptide mass and fragmentation data were searched against a combined forward-reverse IPI database (v3.55) using Mascot 2.2 (Matrix Science). Peptide and protein validation were done using Transproteomic pipeline v3.3sqall (Institute for Systems Biology; [http://tools.proteomecenter.org/software.php]) using a false positive threshold of <1% for protein identifications (See Supplementary file 2 for full list; in red are the ERG interactors also identified as shRNA screen hits).
The Dox-inducible shRNA vector (pLKO-Tet-On) was previously described, as were the sequences of the nontargeting control and ERG shRNA inserts and stable cell line generation(Mounir et al., 2015). PRMT5 shRNA sequences are as follows:
All shRNA were expressed from a puromycin resistant vector. Lentiviral production and cell transduction was as previously described (Mounir et al., 2015).
The Dox-inducible ERG constructs named ERG or ERG DNAx used for stable and inducible expression in 22Rv1 cells were generated as previously described (Mounir et al., 2015). The cDNA expression rescue constructs were cloned into a pRETRO retroviral vector under a CMV promoter and containing a neomycin-IRES-YFP selection cassette. The HA-tagged PRMT5 sequence was generated synthetically to resist knockdown by all PRMT5 shRNAs used in this study and was cloned into a Gateway compatible entry vector. In order to resist knockdown by all PRMT5 shRNA sequences, 6–7 silent mutations (in red below) were introduced into the HA-PRMT5 cDNA sequence to produce a shRNA resistant version (HA-scPRMT5). For shPRMT5-1: the AGGGACTGGAATACGCTAATT target sequence was converted to CGCGATTGGAACACGTTGATT (underlines denote altered bases).
For shPRMT5-2: the GCGGATAAAGCTGTATGCTGT target sequence was converted to CAGAATCAAGCTCTACGCCGT (underlines denote altered bases).
Full sequence of scPRMT5 below:
Generation of the PRMT5 catalytic dead mutant was performed by site-directed mutagenesis (QuickChange II, Agilent) through mutation of G365 to A and R368 to A (Antonysamy et al., 2012).
VCaP, 22Rv1, LNCaP and RWPE-1 cell lines (ATCC) were grown in vendor-recommended media and maintained in a humidified 5% CO2 incubator at 37°C. Doxycycline (Dox, Sigma) was used at 100 ng/ml. Cell proliferation was measured in 6-well plates (Corning) using automated confluence readings (IncuCyte EX, Essen Bioscience). R1881 (Sigma) and charcoal-stripped serum (Omega Scientific) were used where indicated.
2–5 ug of AR, ERG or IgG control antibody was coupled to 1mg of magnetic beads according to manufacturer’s protocol (Invitrogen Dynabeads Antibody Coupling kit#143.11D). After coupling, 1 mg of the antibody/bead mixture was incubated with 1–5 mg of protein lysate overnight under rotation at 4°C. IP samples were then washed with RIPA buffer containing protease/phosphatase inhibitor cocktail for 3–4 washes and resuspended in non-reducing loading buffer, boiled and loaded on a gel for western blot analysis. Immunoprecipitations were performed using the following antibodies: anti-ERG antibody (Epitomics# 2805–1), anti-HA antibody (Roche#11815016001), anti-AR antibody (Thermo Scientific# MA5-13426) or anti-IgG antibody (Rockland Immunochemicals# RL011-0102).
Procedures were previously described (Mounir et al., 2015) and the following antibodies were used at 1:1000 dilutions and incubated overnight at 4°C: ERG (Epitomics# 2805–1), PRMT5 (CST# 2252; SIGMA#P0493), GAPDH (Millipore#MAB374), H4 (CST#2592), AR (Santa Cruz#sc-7305), HA (Roche#11815016001), Symmetric Di-methyl arginine (SDMA, CST#13222), Mono-methyl arginine (MMA, CST#8711), TRIP12 (Abcam#ab86220), EIF4E (CST#9742), CDC42 (BD Transduction Laboratories#610929), HDAC1 (CST#2062), SMARCB1/SNF5 (Bethyl Laboratories# A301-087A), SMARCE1/BAF57 (Bethyl Laboratories# A300–810A).
Generation of labeled cDNA, hybridization to Affymetrix U133plus2 human arrays, and data normalization were performed as described (Mounir et al., 2015). For the candidate signatures in Figure supplement 2A (Malik et al., 2015; Mendiratta et al., 2009; Nelson et al., 2002), a two-tailed fisher’s exact test was used to determine if probesets representing genes in those signatures were under- or over-represented in the set of probesets that were up- or down-regulated at least 1.5-fold compared to expressed but non-differentially-expressed probesets, with a nominal p-value of 0.05 or less. For an unbiased approach (Figure 2A), pathways derived from GO terms and transcription-factor networks were analyzed for overrepresentation via a one-tailed interpolated fisher’s exact test, using genes that varied 1.5-fold or more with a nominal p-value of 0.05 or less compared to all genes represented on the array; Benjamini-Hochberg (BH) correction was then applied to these p-values (Wiederschain et al., 2007). The VCaP microarray dataset (Figure 2A) is available at the NCBI Gene Expression Omnibus (accession number GSE65965).
Black line (Figure supplement 2A) represents expressed probe set position and is ranked by average fold-change. Blue, green, and red lines indicate where the probe sets mapping to genes in the androgen receptor activation signatures appear in our data set and show the cumulative sum of the probe sets in the androgen receptor activation signatures that overlap with our gene list (only the highest expressing probe set was used per gene). The dashed line represents the hypothetical cumulative sum for a random list of genes that are unenriched.
RNA isolation was performed as previously described (Mounir et al., 2015). Taqman reactions (Applied Biosystems) were performed using Gene Expression master mix, FAM-labeled probes for PSA, NKX3-1 and SLC45A3 and VIC-labeled probe for Beta-2-macroglobulin (B2M) as a normalization control. Samples were run on a 7900HT Real-Time PCR machine (Applied Biosystems) and data was analyzed and normalized according to manufacturer’s instructions (2-ΔCt method).
VCaP, 22Rv1 and RWPE-1 cells were treated as specified followed by cross-linking with 1% formaldehyde for 10 min. Cells were next lysed in 1% SDS and sonicated until DNA ladder is below 1 kb (Diagenode). Sheared chromatin was then used for IP with specific primary antibodies (2-4 ug; previously tested) pre-complexed with Protein A/G Dynabeads and incubated overnight under rotation at 4°C. The next day, ChIP samples were washed with RIPA buffer and TE followed by reverse crosslinking using 1%SDS and 30 ug/ml proteinase K (Invitrogen) at 65°C for 6 hr with beads. The eluates were then purified using the QIAquick PCR purification kit (Qiagen) and used for qPCR with the following primer sets:
PSA -4100: acctgctcagcctttgtctc AND ttgtttactgtcaaggacaatcg
PSA -3800: agaattgcctcccaacactg AND cagtcgatcgggacctagaa
PSA -100: cttccacagctctgggtgt AND aaaccttcattccccaggac
PSA +700: agccccagactcttcattca AND atgcagatttggggaatcag
NKX3-1 -2800: gagagcagctgttcctccac AND acgagccttttccacctttc
NKX3-1 -200: agggaggagagctggagaag AND tcctccctaggggattcct
NKX3-1 +2150: accaggatgaggatgtcacc AND cagggacagagagagccttg
NKX3-1 -+10800: tctctcgttggctcctgatt AND ccagcttttgttccttcctg
NKX3-1 +62100: cggtttattgcccatgaaga AND aacagggctcacagtgcttt
VCaP cells harboring Dox inducible shRNAs targeting PRMT5 where grown to 80% confluency and re-seeded (day 0) into full media containing 100 ng/ml Doxycycline. On day 3 media was replaced with ‘hormone reduced’ media, containing 100 ng/ml Doxycycline. Cells were stimulated with adding indicated ligands (DHT (Sigma), R1881 [Sigma]) or vehicle on day 4 and harvesting of cells was performed on day 5. Cells were harvested by fixation using 1% methanol free formaldedyde (Polysciences, Cat#18814) in PBS at room temperature. Fixation was stopped after 8 min by replacing fixation buffer with ice cold PBS containing 125 mM glycine and 5 mg/ml BSA. Cells were further washed once using ice-cold PBS and re-suspended into 500 ul of PBS containing Complete Protease Inhibitors (Roche). Cells where then pelleted and supernatant removed, and the resulting pellet either snap frozen in liquid Nitrogen or immediately re-suspended in lysis buffer for further processing.
For ChIP the resulting cell lysate was sonicated using a Covaris E210 instrument according to manufacturers recommendations. Each ChIP reaction was performed using soluble fraction chromatin corresponding to 7.5 ug purifed DNA and 4 ug of antibodies. Antibodies were allowed to bind overnight before capture on protein A magnetic beads (Invitrogen, Dynal). Bound beads where washed 4 times in RIPA buffer containing 500 mM LiCl, and 2 times with TE buffer before being re-suspended in containing 100 mM NaHCO3 and 1% (w/v) SDS. Crosslink reversal was done at 65C°C for 6 hr and ChIP DNA were isolated using DNA purification beads (MagBio). ChIP-seq libraries were generated using the KAPA HTP library preparation kit (Kapa biosciences). All handling of samples after sonication was done on using a Sciclone NGS Workstation (P/N SG3-31020-0300, PerkinElmer). Sequencing was performed on an Illumina NextSeq500 instrument.
Reads passing Illumina standard QC were mapped to genome version Hg19 using BWA, and binding sites (‘peaks’) were identified using MACS2, evolutionary conservation scores at peak locations was calculated using Phastcons and enriched DNA motifs using MDscan and Seqpos. These were performed using the ChiLin QC pipeline (liulab.dfci.harvard.edu/WEBSITE/software). Peaks from AR ChIP-sequencing samples with a MACS2 enrichment score higher than 10 were extended to a uniform 400 bp across all samples and overlapping peaks where collapsed to generate a union of all peaks. This resulted in a Cistrome of 25,593 peaks that were used in all genome wide ChIP analyses. Using the features of this Cistrome as GTF, read counts from BWA mapped Bam files were processed using the Qlucore 3.1.19 software. Heatmaps and statistical test (two-sided t-test using correction for multiple hypothesis testing) of differential binding scores on the 25,593 features were performed in Qlucore v3.1.19. The AR and ERG ChIPseq datasets are available at the NCBI Gene Expression Omnibus (Accession number GSE79128).
RWPE-1 cells were treated as specified, fixed with 4% paraformaldehyde for 45 min at room temperature (Electron Microscopy Sciences), blocked with 5% goat serum, 0.5% Triton X-100 in PBS for 2 hr and incubated with 1:50 dilutions of AR (LSBIO #LS-C87494) and symmetric di-methyl arginine antibodies (CST#13222) in 5% goat serum and 0.05% Triton X-100. Fixed samples were incubated overnight at 4°C in primary antibody before incubations with proximity ligation assay (PLA) secondary antibodies (Duolink, Sigma).The secondary antibody incubation, ligation, amplification and final wash steps were performed according to the manufacturer’s specifications. Confocal microscopy was performed using an LSM 510 META (Carl Zeiss, Inc.) with a 40x C-Apochromat objective, NA 1.2. Images were collected and processed using Zen software (Carl Zeiss, Inc.).
Full length Homo sapiens androgen receptor (AR) sequence (transcript variant 1; NM_000044.3) was synthesized to include a 5’ NotI and 3’ BamHI sites and used as a template for the mutagenesis of each arginine in the ligand binding domain of AR into lysine (QuikChange XL site-directed mutagenesis kit, Agilent). Following sequence verification to ensure mutation incorporation, each AR mutant sequence was cloned into the pLVX vector via the 5’ NotI and 3’ BamHI as previously described (Mounir et al., 2015).
Immunofluorescence of RWPE-1 cells was performed by fixing cells for 45 min at room temperature by adding 4% paraformaldehyde (Electron Microscopy Sciences) and incubated with 1:50 dilution of AR antibody (LSBIO #LS-C87494). Confocal microscopy was performed using an LSM 510 META (Carl Zeiss, Inc.) with a 40x C-Apochromat objective, NA 1.2. Images were collected and processed using Zen software (Carl Zeiss, Inc.).
The heterodimeric structure of PPARγ-RXRα (PDB Code: 3DZY) was used as a template to overlay individual domain structures of AR including the DBD in complex with DNA (PDB Code: 1R4I) and the LBD in complex with coactivator peptide TIF2(iii) and ligand R1881 (PDB Code: 2AO6). A superposition was achieved using secondary structure matching in COOT (Emsley et al., 2010). The AR DBD was superposed onto chain A of RXRα and the AR-LBD was superposed onto chain B of PPARγ, resulting in the final overlay shown in Figure 4A. A dimethylated Arg was superposed onto R761 (NCBI Reference Sequence: NP_000035.2; some publications may refer to it as R760 when using the previous Reference Sequence number) within the AR LBD using least squares fit and matching only mainchain atoms, after which the non-methylated Arg was removed from the resulting model. While structure determination of AR containing its intact DBD and LBD has remained elusive, structural data of each individual domain, as well as intact structures within the nuclear receptor family, lead to a valuable understanding of interdomain communication. A previous study utilized the heterodimeric structure of PPARγ and RXRα to apply in silico three-dimensional alignment and docking analysis, followed by mutational analysis, to propose a DBD-LBD interface within AR, including R761 (Figure 4A) (Chandra V et al., 2008; Helsen et al., 2012). We hypothesize that R761 is involved in key interactions at this interface and that its methylation would add hydrophobicity, eliminating any polar interactions, as well as steric bulk. This disruption at the DBD-LBD interface could result in destabilization of the quaternary structure, resulting in an inhibitory effect on AR activation. While AR R761K closely mimics the properties of wt AR, the substitution eliminates the possibility of methylation by PRMT5 and, therefore, eliminates the possibility of this type of disruption, accounting for the observed increase in activation.
The gene encoding human AR LBD (residues 663–919), was inserted into a pGEX-6P-1 vector and expressed as a GST-tagged fusion protein in BL21 Star (DE3) cells. Cells were grown in TB2 medium containing 10 µM dihydrotestosterone (DHT) and induced with 1 mM IPTG for 14 –16 hr at 16°C. Cells were resuspended in buffer A containing 50 mM Tris-HCl (pH 7.3), 150 mM NaCl, 10% glycerol, 0.25 mM TCEP, and 10 µM DHT, to which 50 µg/ml DNase I and protease inhibitor cocktail (Roche) were added. Cells were lysed using an M-110L Microfluidizer at 18,000 psi, followed by the addition of 0.5% CHAPS to the lysate prior to high speed centrifugation. For one-step batch purification, the soluble extract was incubated with 2 ml of glutathione sepharose 4 fast flow medium (GE Healthcare) for 1 hr at 4°C with rotational mixing. The sepharose medium was washed in buffer A with the addition of 0.5% CHAPS. Elution was accomplished by resuspending and incubating the media for 10 min in the wash buffer plus 10–20 mM reduced glutathione. The eluted fractions were then combined and concentrated to 0.3 mg/ml.
PRMT5 enzymatic activity was assessed by monitoring S-adenosyl-L-homocysteine (SAH) product formation utilizing liquid chromatography-tandem mass spectrometry (LC-MS/MS). 0.5–1 uM of PRMT5/MEP50 recombinant enzyme (BPS, cat#51045) was incubated with 50 uM SAM, 2 uM GST-AR LBD and/or 5 uM ETS or PNT ERG protein for 2 hr at 37°C. Reactions were quenched to 0.1% HCOOH followed by addition of [β,β,γ,γ-2H4]-SAH (SAH-D4) in 20% DMSO as an internal standard for MS quantification. Samples were sonicated with a Hendrix SM-100 sonicator (Microsonics Systems) and centrifuged. SAH was separated from the reaction mixture by reversed phase chromatography using polar endcapped C18 reversed phase columns (Synergi Hydro-RP, 2.5 μm, 100 Å, 20 x 2 mm, Phenomenex) and detected using a 4000 QTRAP Hybrid Triple Quadrupole/Linear Ion Trap LC-MS/MS system (AB Sciex).
Crystal structure of the human PRMT5:MEP50 complexProceedings of the National Academy of Sciences of the United States of America 109:17960–17965.https://doi.org/10.1073/pnas.1209814109
Molecular determinants of resistance to antiandrogen therapyNature Medicine 10:33–39.https://doi.org/10.1038/nm972
Small molecule regulators of protein arginine methyltransferasesThe Journal of Biological Chemistry 279:23892–23899.https://doi.org/10.1074/jbc.M401853200.
Lineage dependency and lineage-survival oncogenes in human cancerNature Reviews. Cancer 6:593–602.https://doi.org/10.1038/nrc1947
Evidence for DNA-binding domain--ligand-binding domain communications in the androgen receptorMolecular and Cellular Biology 32:3033–3043.https://doi.org/10.1128/MCB.00151-12
Functional epigenetics approach identifies BRM/SMARCA2 as a critical synthetic lethal target in brg1-deficient cancersProceedings of the National Academy of Sciences of the United States of America 111:3128–3133.https://doi.org/10.1073/pnas.1316793111
A causal role for ERG in neoplastic transformation of prostate epitheliumProceedings of the National Academy of Sciences of the United States of America 105:2105–2110.https://doi.org/10.1073/pnas.0711711105
Targeting the MLL complex in castration-resistant prostate cancerNature Medicine 21:344–352.https://doi.org/10.1038/nm.3830
Genomic strategy for targeting therapy in castration-resistant prostate cancerJournal of Clinical Oncology 27:2022–2029.https://doi.org/10.1200/JCO.2008.17.2882
Maintaining and reprogramming genomic androgen receptor activity in prostate cancerNature Reviews. Cancer 14:187–198.https://doi.org/10.1038/nrc3678
TMPRSS2:ERG blocks neuroendocrine and luminal cell differentiation to maintain prostate cancer proliferationOncogene, 34, 10.1038/onc.2014.308.
The program of androgen-responsive genes in neoplastic prostate epitheliumProceedings of the National Academy of Sciences of the United States of America 99:11890–11895.https://doi.org/10.1073/pnas.182376299
The discovery of common recurrent transmembrane protease serine 2 (tmprss2)-erythroblastosis virus E26 transforming sequence (ETS) gene fusions in prostate cancer: Significance and clinical implicationsAdvances in Anatomic Pathology 16:145–153.https://doi.org/10.1097/PAP.0b013e3181a12da7
Genome-wide analysis of ets-family dna-binding in vitro and in vivoThe EMBO Journal 29:2147–2160.https://doi.org/10.1038/emboj.2010.106
Contribution of polycomb homologues bmi-1 and mel-18 to medulloblastoma pathogenesisMolecular and Cellular Biology 27:4968–4979.https://doi.org/10.1128/MCB.02244-06
Scott A ArmstrongReviewing Editor; Memorial Sloan Kettering Cancer Center, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your work entitled "ERG signaling in prostate cancer is driven through PRMT5-dependent methylation of the Androgen Receptor" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Kevin Struhl as the Senior Editor.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission. As you will see below, the reviewers have some concerns that will require quite an additional effort to address. Ordinarily, we allow just two months to return a revised version but we are concerned that this required experiments may take longer to complete. To help us assess the likelihood of a successful completion of the required additional work, please send us a note indicating how you plan to proceed and an estimate of the time it will take to complete the work.
In this manuscript, the authors present two major findings: 1) PRMT5 methylates AR at R761 and these methylation results in decreased AR recruitment at target genes and decreased AR function and 2) ERG targets AR to be methylated. Overall, the findings address long-standing questions of ERG function in prostate tumorigenesis and possibly posit new therapies to ERG translocated prostate cancers. This would be certainly relevant to the prostate field. In addition, it is appropriate for a general audience for implications of a general enzyme that affects specific programs. The experiments are well done and logical. The manuscript is well written. However, at this point, there are several critical experiments that should be performed to solidify the conclusion.
1) Definition of the ERG and AR binding genome wide with and without PRMT5 and ERG knockdown and correlation with the gene expression is very important for the conclusion. Is AR binding specifically depleted in ERG co-bound sites? The changes genome wide may not reflect the selected sites, as highlighted by the effect of FOXA1 on AR. This does require time and effort. But AR and ERG ChIP-seq is well optimized and ChIP-seq today is cost effective so this is eminently doable.
2) One weakness of the manuscript, intrinsic to prostate cancer research, is the use of a single ERG positive cell line. So any hit of the screen could be ERG specific, or even more likely, VCAP specific. Since there are no other available ERG dependent cell lines, the authors should note in the manuscript this major caveat. It would be very helpful if the authors could examine xenografts or perform some validation of the mechanisms in primary human tissue samples. For example, authors can IP AR from several ERG positive and ERG negative prostatectomy specimens to assess for methylation.
3) Mass spec evidence of methylation at R761 would be helpful. Since you know what you're looking for, IP of AR in VCAP and another ERG negative line and looking for this should not be hard and would provide definitive evidence. The mutagenesis is highly suggestive, but could still change conformation in a way such that R761 is not the direct methylation site.https://doi.org/10.7554/eLife.13964.023
- Thomas Westerling
- Gilles Buchwalter
- Myles Brown
All authors except Thomas Westerling, Gilles Buchwalter, and Myles Brown received no external funding for this work.
We thank Huili Zhai, Veronica Sanz-Vash, Victor Lin, Yi Yang, Alyson Freeman, Konstantinos Mavrakis, Kenneth Crawford, Jessica Cherry, Aron Jaffe, Margaret McLaughlin, Nanxin Li, Frank Stegmeier, Pascal Fortin, Gyorgy Petrovics, Shiv Srivastava, and Wojciech Wrona for helpful discussions.
- Scott A Armstrong, Reviewing Editor, Memorial Sloan Kettering Cancer Center, United States
© 2016, Mounir et al
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.