Abstract
Recent evidence has implicated the transmembrane co-receptor neuropilin-1 (NRP1) in cancer progression. Primarily known as a regulator of neuronal guidance and angiogenesis, NRP1 is also expressed in multiple human malignancies, where it promotes tumor angiogenesis. However, non-angiogenic roles of NRP1 in tumor progression remain poorly characterized. In this study, we define NRP1 as an androgen-repressed gene whose expression is elevated during the adaptation of prostate tumors to androgen-targeted therapies (ATTs), and subsequent progression to metastatic castration-resistant prostate cancer (mCRPC). Using short hairpin RNA (shRNA)-mediated suppression of NRP1, we demonstrate that NRP1 regulates the mesenchymal phenotype of mCRPC cell models and the invasive and metastatic dissemination of tumor cells in vivo. In patients, immunohistochemical staining of tissue microarrays and mRNA expression analyses revealed a positive association between NRP1 expression and increasing Gleason grade, pathological T score, positive lymph node status and primary therapy failure. Furthermore, multivariate analysis of several large clinical prostate cancer (PCa) cohorts identified NRP1 expression at radical prostatectomy as an independent prognostic biomarker of biochemical recurrence after radiation therapy, metastasis and cancer-specific mortality. This study identifies NRP1 for the first time as a novel androgen-suppressed gene upregulated during the adaptive response of prostate tumors to ATTs and a prognostic biomarker of clinical metastasis and lethal PCa.
Similar content being viewed by others
Introduction
Prostate cancer (PCa)-associated mortality is due to therapy-resistant metastatic tumor burden. Although prostate-confined tumors are often treatable by surgery and/or radiation, distal metastatic disease remains incurable. Local recurrence following radical prostatectomy (RP) can be treated by salvage radiation therapy. However, this approach fails in some patients if tumor cells acquire radiation resistance, or if undetected (occult) distal metastases have already formed at the time of radiotherapy. Once metastatic disease is established, treatment relies on androgen-targeted therapies (ATTs), which exploit the androgen dependence of PCa cells. Although ATTs provide a temporary remission (usually 2–3 years), they inevitably promote adaptation of tumor cells to low androgen conditions, giving rise to lethal castration-resistant PCa (CRPC), a highly aggressive and metastatic PCa variant. Second- and third-generation ATTs such as abiraterone (Zytiga) and enzalutamide (Xtandi) delay progression to CRPC, but are not curative.1 Standard cytotoxic chemotherapies such as docetaxel (Taxotere) and cabazitaxel (Jevtana) offer additional survival benefits of only 2–22 months before inevitable relapse,2, 3 whereby greater survival benefits (>10–22 months) have been reported when administered early with ADT.4, 5 Nevertheless, it is clear that sequential treatment failure drives the emergence of increasingly aggressive, therapy-resistant PCa, progressively shortening time to relapse.
Although the mechanisms underlying treatment failure and progression to CRPC are not completely understood, it is increasingly apparent that a population of tumor cells can undergo an adaptive response to ATTs. This adaptive tumor response leads to the activation of alternative tumor-promoting pathways, alterations to androgen receptor (AR) function, and scavenging of adrenal and intra-tumoral steroids.6, 7, 8 The identification of the molecular determinants of the prostate tumor adaptation to ATTs relevant for the progression to CRPC will likely identify new candidates for therapeutic intervention. Used alongside current and emerging ATT and cytotoxic treatment combinations, these new therapies may lead to the improved clinical management of metastatic CRPC (mCRPC) and patient outcomes.
In this study, we have utilized the genome-wide transcriptional profiling of cell lines, xenografts and clinical PCa samples to identify a gene signature of tumor adaptation to ATTs. Refinement of this gene set has implicated the transmembrane glycoprotein neuropilin-1 (NRP1) in the adaptive response to ATTs and progression to CRPC. Here, we report that NRP1 is an androgen-suppressed gene overexpressed in therapy-resistant tumors and mCRPC. In cell line models of mCRPC, the short hairpin RNA (shRNA)-mediated inhibition of NRP1 expression led to a significant reduction in their invasive and metastatic capacity. Moreover, analysis of several PCa patient cohorts identified NRP1 as an independent prognostic indicator of early biochemical recurrence (BCR) following radiation therapy, metastasis and PCa-specific mortality. This study not only provides new insights into the function of NRP1 expressed by tumor cells, but also supports the rational use of anti-NRP1 agents alongside current ATT and cytotoxic regimes in the treatment of men with advanced PCa.
Results
Identification of a transcriptional signature representing the adaptive tumor response to ATTs
To characterize the adaptive response of PCa cells to ATTs, a human PCa xenograft (LNCaP) model of CRPC9, 10 was used to identify a cluster of genes significantly upregulated with castration and remaining elevated in CRPC (Figure 1a). Of these genes, a subset was repressed by the androgen dihydrotestosterone (DHT) and upregulated by enzalutamide (ENZ) in LNCaP cells in vitro (Figure 1b). Overlay of this data set with genes upregulated in human mCRPC versus localized PCa11 revealed a core transcriptional cluster of 120 androgen-regulated genes overexpressed in CRPC (Figure 1b, Supplementary Table 1). In an effort to identify clinically robust gene candidates, we analyzed the 120 gene set in publicly available gene expression profiling of patient metastatic versus localized PCa12, 13, 14, 15 (Figure 1c). This revealed NRP1 to be significantly upregulated (P<0.05; ⩾1.5 fc) in five of five data sets analyzed, with 1.9–4.7 times higher expression in metastatic samples than localized PCa samples (Figure 1d). Although studies have reported NRP1 to facilitate cancer progression in multiple tissues,16, 17, 18, 19, 20, 21 the biological role of NRP1 when expressed by PCa cells remains poorly understood. As such, we chose to characterize NRP1 expression during the adaptive response to ATTs and its role in prostate tumor progression.
NRP1 expression is suppressed by active androgen signaling
In vitro validation studies confirmed NRP1 mRNA expression to be suppressed in LNCaP cells by 48-h treatment with DHT (10 nM) or the synthetic androgen metribolone (R1881; 1 nM) (Figure 2a). To determine whether the AR directly binds to the NRP1 gene loci to mediate its transcriptional repression, analysis of a chromatin immunoprecipitation (ChIP) sequencing data set22 demonstrated AR enrichment at five putative AR-binding sites in multiple human prostate tumors (Figure 2b). We conducted ChIP assays coupled with quantitative PCR using primers that amplified these five putative binding sites, detecting significant AR enrichment at distal (distal_1) and intron 12 genomic regions, as well as a positive control locus (KLK3/prostate specific antigen (PSA)) in LNCaP cells cultured in the presence of androgens (Figure 2c). The repressive effect of AR on NRP1 expression was relieved by AR knockdown using small interfering RNA (Figure 2d and Supplementary Figure 1) or treatment with the AR inhibitors bicalutamide (Bic) and Enz (Figure 2e). In concordance, NRP1 mRNA (Figure 2f) and protein (Figure 2g) levels rose in response to in vitro androgen deprivation with 5% charcoal-stripped serum (CSS) supplementation over 7 days. Notably, the re-introduction of DHT attenuated the CSS-mediated increase in NRP1 mRNA and protein levels (Figures 2f and g), however, NRP1 levels remained elevated compared with basal levels in androgen containing media (fetal bovine serum (FBS)). Consistent with this in vitro data, surgical castration of mice significantly elevated NRP1 mRNA levels in LuCAP35 tumor xenografts as compared with sham-castrated controls (Figure 2h). Taken together, these data indicate that active AR signaling may directly mediate the suppression of NRP1 expression.
NRP1 is upregulated during the adaptive response to ATTs and progression to CRPC
We have previously established the LNCaP xenograft model to mimic the in vivo progression to CRPC following castration.9, 10 The transcriptional profiling of LNCaP tumors following castration of host mice identified NRP1 levels to peak at post-castration PSA nadir and remain elevated with progression to castration resistance (CR; Figure 3a). Increased NRP1 protein expression was also observed in an androgen-independent LNCaP variant cell line generated following long-term culture in androgen-deprived conditions (LNCaP-AI),23 compared with their androgen-dependent parental cells (Figure 3b). Next, we assessed the expression levels of NRP1 in clinical samples of CRPC. Analysis of gene expression data obtained from laser capture microdissected samples from treatment naive or hormone refractory primary tumor samples24 identified significantly higher NRP1 mRNA expression in men with hormone refractory progressive disease (PD; Figure 3c). This was confirmed by transcriptional profiling of mCRPC samples by Grasso et al.,11 which showed a significant upregulation of NRP1 expression in heavily pre-treated samples of mCRPC compared with either localized PCa or benign tissue samples (Figure 3d). No significant difference in expression was observed between benign and localized PCa samples (Figure 3d). Moreover, RNA-sequencing data of 150 mCRPC bone or soft tumor biopsies25 identified the substantial expression of NRP1 in these tumors to be further increased in patients treated with abiraterone (Abi) or Enz versus neither treatment (Figure 3e). Collectively, these data establish NRP1 to be dynamically regulated during the prostate tumor adaptation to ATTs and progression to mCRPC.
NRP1 is required for the invasion and metastatic dissemination of mCRPC cell models
As few reports have described a functional role for NRP1 when expressed directly by PCa cells, we next investigated the effect of NRP1 on the invasive and metastatic phenotype. Flow cytometric profiling revealed cell surface NRP1 levels to be relatively low in the benign prostate epithelial cell lines BPH-1 and RWPE-1, but elevated in tumor cells with increasing invasive and metastatic potential, with highest expression seen in the castration-resistant and metastatic PC3 line (Figure 4a). To determine the effect of NRP1 suppression on the phenotype of PCa cells, two stable NRP1 knockdown PC3 cell lines (shNRP1(1) and shNRP1(2)) were generated using the pLKO.1 lentiviral vector (Figure 4b). No significant difference in proliferation was observed between NRP1 knockdown and control cell lines as measured by either cell confluence or DNA content quantification (Figure 4c). However, we observed that NRP1 knockdown caused PC3 cells to grow as more compact and non-invasive colonies over a 10-day period in both monolayer (Figure 4d) and modified three-dimensional (3D) On-top Matrigel cultures (Figure 4e). A reduction in invasive capacity following NRP1 suppression in PC3 cells, as well as the additional metastatic and castrate resistant DU145 cell line, was confirmed using wound scratch invasion studies (Supplementary Figure 2).
As previous reports have implicated NRP1 in the epithelial-mesenchymal transition phenotype,26, 27 a developmental program capable of endowing tumor cells with invasive properties,28 we next assessed the mesenchymal properties of PC3 cells upon NRP1 knockdown. The expression of the prototypical mesenchymal and epithelial markers, vimentin and E-cadherin, respectively, were examined using the In-Cell Western technique LI-COR Biosciences (Lincoln, NE, USA) on monolayer cultures (Figure 4f) or immunofluorescence staining of cells grown in 3D On-top Matrigel cultures (Figure 4g). In both monolayer and 3D cultures, the suppression of NRP1 levels resulted in reduced vimentin and increased E-cadherin protein expression (Figures 4f and g). The reduction in NRP1 expression was also associated with increased cortical actin staining (Figure 4g), a characteristic of epithelial cell types. Next, we assessed the impact of NRP1 knockdown on metastatic dissemination. Control (shCntrol) and NRP1 knockdown PC3 cells were xenografted into the yolk sacs of dechorionated 2-day post-fertilization wild-type zebrafish embryos and metastasis, measured as cell dissemination outside the yolk sac, was imaged 5 days later. In comparison with control cells, which disseminated toward the head and tail of the fish, a significantly reduced proportion of zebrafish were positive for metastasis following injection of the NRP1 knockdown cell lines (Figures 4h and i). Collectively, these data provide evidence that NRP1 is required for the invasion and metastatic dissemination of mCRPC cell models, which may be mediated via its regulation of the mesenchymal phenotype.
Increased NRP1 expression is associated with primary tumor progression and therapy failure
Previous studies have reported a positive association between NRP1 expression and Gleason grade in immunohistochemical (IHC) analyses of small patient cohorts (5–17 specimens per Gleason grade).16, 18 We expanded upon these studies by performing IHC staining for NRP1 using a Gleason grade tissue microarray containing 176 patient specimens (cohort summarized in Supplementary Table 2). Results from this larger cohort supported previous observations, with strong NRP1 staining more frequently observed in higher Gleason grade tumors (Figure 5a). In addition, analysis of RNA-sequencing data from The Cancer Genome Atlas (TCGA) Prostate Adenocarcinoma (PRAD) cohort (http://xena.ucsc.edu) revealed increased NRP1 mRNA expression to be associated with increasing pathological stage and node status in a large cohort (n=498) of primary tumor samples (Figure 5b). RNA sequencing of the TCGA PRAD cohort also identified significantly elevated NRP1 expression in tumor samples from men with progressive disease (PD) following primary therapy versus those who had complete response (CR; Figure 5b). Kaplan–Meier analysis of the TCGA PRAD cohort showed significantly lower probability of relapse-free survival in men with higher than median NRP1 expression (Figure 5c).
NRP1 predicts biochemical failure in patients after postoperative radiation therapy
As TCGA analysis revealed elevated NRP1 expression in the primary tumors of patients who failed primary therapy, we next investigated the significance of NRP1 as a predictive biomarker for BCR after primary therapies. NRP1 expression at RP was analyzed with reference to patient outcome in a cohort of 130 patients who underwent post-RP adjuvant or salvage radiotherapy at the Kimmel Cancer Center, Thomas Jefferson University (TJU), Philadelphia, PA, USA29 (Figure 5d). All patients were diagnosed with pT3 or margin positive disease at the time of radiotherapy. Significantly higher NRP1 expression was observed in patients who developed post-radiotherapy BCR compared with patients who did not (Figures 5e and f). This was confirmed by Kaplan–Meier survival analysis revealing patients with high NRP1 expression to have a rapid progression to BCR, suggestive of a greater presence or outgrowth of occult metastases (Figure 5g). Finally, multivariate analysis identified high NRP1 expression as a significant independent predictor of BCR (P=0.019, Supplementary Table 3). Taken together, these data provide novel insight into the expression of NRP1 in clinical samples of primary PCa and its positive association with aggressive clinicopathological parameters.
NRP1 is an independent prognostic biomarker of metastatic progression and cancer-specific mortality
As a small minority of patients from the TJU cohort developed distal metastases during follow-up, we analyzed the ability of NRP1 to predict metastasis in a larger cohort. Hence, NRP1 expression at RP was analyzed for its ability to predict metastasis, defined by positive bone or computed tomography scan, in a cohort of 545 patients who had undergone RP at the Mayo Clinic30 (Figure 6a). NRP1 expression was found to be significantly higher in patients who were positive for metastasis (Figures 6b–c) and prostate cancer-specific mortality (PCSM; Figures 6d and e). Moreover, multivariate analysis identified high NRP1 expression at RP as a significant independent predictor of both metastatic progression (P=0.008) and PCSM (P=0.013; Table 1).
To further validate these findings, we performed a similar analysis in a Natural History cohort of 644 intermediate- and high-risk men who underwent RP at Johns Hopkins Medical Institutions (JHMI).31 In this unique cohort, patients received no adjuvant or salvage therapy following RP until metastatic progression. Again, high NRP1 expression was associated with lower probability of metastasis-free survival (Figure 7a) and was an independent predictor of metastatic progression (P=0.048, Table 1). Within a subgroup of 211 patients who were positive for BCR (Figure 7b), NRP1 was associated with a shorter time to PCSM (Figure 7b) and was a significant predictor of PCSM (P=0.034; Table 1). Overall, these results define NRP1 as a novel biomarker for identifying patients at risk of metastatic progression and death from PCa after RP.
Discussion
Despite the introduction of new androgen-targeting agents ENZ and Abi acetate into treatment regimens for CRPC, it remains incurable.1, 6 Androgens are potent growth and differentiation factors for prostate tissue, and loss of this differentiation pressure is an unintended consequence of ATTs. This can result in the initiation of de-differentiation and trans-differentiation cell plasticity programs, which enable survival under low androgen conditions.32 Hence, identifying biomarkers and molecular determinants of the adaptive progression to metastatic castrate-resistant disease is critical for the development of more effective treatments. Here, we identify NRP1 as one of a subset of androgen-suppressed genes persistently overexpressed in CRPC in both an in vivo human xenograft model of CRPC and clinical mCRPC.
Although ATTs initially suppress androgen signaling, this pathway is re-activated in CRPC.8, 33 A possible explanation of the sustained NRP1 overexpression in CRPC despite re-activated androgen signaling may involve the transcriptional activity of CRPC-associated AR splice variants.34 A recent report identified NRP1 to be upregulated by AR variants expressed in 22Rv1 cells, but not full-length AR, in the absence of androgens.35 In addition, NRP1 formed part of a 297 gene signature capable of distinguishing localized PCa from CRPC, and was 1 of 34 probes significantly associated with shorter time to post-RP BCR.35 In concordance, we report NRP1 to be associated with shorter time to relapse following primary therapy in the TCGA PRAD cohort and an independent predictor of PCa recurrence following post-RP radiation therapy in a 130-patient cohort from the TJU.
We report for the first time that high NRP1 expression in primary tumors at the time of RP is an independent predictor of clinical metastasis and cancer-specific mortality in two clinical cohorts from the Mayo Clinic and JHMI. Further studies will be required to decipher whether NRP1 is solely a biomarker for identifying patients at risk of metastatic progression and death from PCa or actively contributes to metastatic and therapy-resistant disease. We demonstrate NRP1 is required for the full metastatic potential of tumor cells using an in vivo Zebrafish model. Although studies using xeno-transplantation into mice and transgenic models will be necessary for validation, the Zebrafish model is becoming increasingly recognized as a rapid, robust and inexpensive assay that can faithfully recapitulate the metastatic potential of numerous human cancer cell models.36 Given the relatively short duration of studies (normally <1 week) that can be performed using the Zebrafish model, mouse studies will be critical for deciphering the role of NRP1 in the formation of occult metastases and their outgrowth over more extended periods.
Collectively, the data reported herein support the rational therapeutic targeting of NRP1 to inhibit prostate tumor progression. A recent Genentech-led Phase 1b clinical trial failed to demonstrate efficacy of the NRP1 antibody MNRP1685A (Vesencumab), targeted against the vascular endothelial growth factor (VEGF)-binding site of NRP1, as a mono drug therapy in a mixed cancer patient cohort, which did not include PCa patients.37 In combination with the VEGF inhibitor bevacizumab, MNRP1685A administration led to a high incidence of proteinuria, prompting withdrawal of the drug. As we report for the first time that NRP1 is expressed in a dynamic manner during the adaptive response of tumors to ATTs, NRP1-targeted compounds may prove more efficacious in a PCa context when administered as adjuvant therapies alongside agents targeting the androgen signaling axis. In the absence of additional VEGF inhibitors, renal side effects may be expected to be less pronounced.
Ligands outside the VEGF family, such as transforming growth factor-β,38 platelet-derived growth factor-BB,39 fibroblast growth factors-1, -2 and -4 and hepatocyte growth factor,40 also have a role in mediating the biological effects of NRP1. As such, agents aimed at more broadly inhibiting NRP1 activity, rather than specific blockade of the VEGF-binding site, may prove a more effective strategy in targeting NRP1-mediated PCa progression. Indeed, NRP1 can promote many aspects of tumorigenesis, such as angiogenesis, cell survival, migration, invasion and chemo-resistance.19, 41, 42 Further studies will be required to elucidate the contribution of various downstream signaling pathways to the aggressive phenotype mediated by NRP1 in the PCa setting. However, the therapeutic potential of targeting NRP1 has been highlighted by recent studies showing the blockade of NRP1 to inhibit the spread and growth of experimental models of human medulloblastoma, and non-small cell lung carcinoma xenografts.43, 44
In conclusion, this study provides the first comprehensive clinical evaluation of NRP1 expression in human PCa. We have identified NRP1 at RP as a prognostic biomarker of shorter time to BCR, metastasis and cancer-specific mortality, and report that NRP1 functions to enhance the metastatic potential of PCa cells. We provide the first evidence that NRP1 expression is upregulated by ATTs, the standard clinical treatment for recurrent and metastatic disease. Collectively, our findings provide the preclinical data to support the use of anti-NRP1-targeted therapies as novel co-targeted therapies to be used in an adjuvant setting alongside current ATT and cytotoxic regimes in the treatment of men with advanced disease.
Materials and methods
Cell culture
LNCaP, DU145 and PC3 cells were sourced from the American Type Culture Collection (Manassas, VA, USA). C4-2B and BPH-1 cells were obtained from Dr Leland Chung, Cedars-Sinai Medical Center, Los Angeles, CA, USA, and Dr Simon Hayward, Vanderbilt University Medical Center, Nashville, TN, USA, respectively. LNCaP-AI cells were obtained from Dr Ralph Buttyan, Vancouver Prostate Centre, Vancouver, BC, Canada. All cell lines undergo mycoplasma testing on a quarterly basis. LNCaP, DU145, PC3, C4-2B, BPH-1 and LNCaP-AI cells were grown in RPMI medium with 10% FBS. LAPC4 cells were cultured in Iscove’s modification of Dulbecco’s medium with 10% FBS. RWPE-1 cells were grown in keratinocyte medium containing recombinant human epidermal growth factor (5 ng/ml) and bovine pituitary extract (50 ng/ml). shRNA-transfected PC3 and DU145 cells were maintained in an additional 1 μg/ml puromycin.
Androgen and anti-androgen treatments
LNCaP cells were seeded into six-well dishes at 9x104 cells per well in growth medium before medium was replaced with RPMI containing 5% CSS 72 h later. After 48 h, medium was replaced with fresh RPMI/5% CSS and cells were treated with 10 μM ENZ or Bic (Selleck Chemicals, Houston, TX, USA) in the presence of 10 nM DHT or ethanol vehicle control, and RNA was harvested after 48 h. For AR knockdown, LNCaP cells were cultured in growth medium for 72 h before being transiently transfected with AR On-target Plus siRNA (target sequence 5′-CGAGAGAGCUGCAUCAGUU-3′, GE Dharmacon, GE Healthcare, Lafayette, CO, USA) or nonspecific control small interfering RNA (GE Dharmacon, GE Healthcare) at a final concentration of 50 nM using Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA, USA). Medium was replaced 48 h post-transfection with RPMI/5% CSS, with or without 10 nM DHT, and RNA collected 48 h later.
NRP1 shRNA knockdown
NRP1 shRNA pLKO.1 lentiviral vectors with mature anti-sense sequences 5′-AATACTAATGTCATCCACAGC-3′ (shNRP1(1)) and 5′-ATATAAGTGCATTCAAGGCTG-3′ (shNRP1(2)), as well as a control sequence targeting firefly luciferase (shCntrl), were obtained from Thermo Fisher Scientific. Viral particles were produced as previously described.45
Evaluation of NRP1 protein expression in clinical PCa tissues
A Gleason grade tissue microarray comprising 176 PCa patient samples was obtained from the Vancouver Prostate Centre Tissue Bank (see Supplementary Table 2 for a summary of clinicopathological features). Specimens were obtained from patients following informed consent using a protocol approved by the Clinical Research Ethics Board of the University of British Columbia and the BC Cancer Agency. Tissue microarray construction, IHC using a rabbit monoclonal antibody against human NRP1 (EPR3113, 2621-1, Epitomics, Burlingame, CA, USA) and evaluation of staining intensity were performed as previously described.46
Quantitative PCR
Total RNA was extracted with the RNeasy Mini Kit (Qiagen, Hilden, Germany) before reverse transcription with SuperScript III Reverse Transcriptase (Thermo Fisher Scientific). Quantitative PCR was performed using SYBR Green (Applied Biosystems, Foster City, CA, USA) using ViiA 7 or 7900 HT Fast Real Time PCR systems (Applied Biosystems). Gene expression was determined by the comparative Ct method, and normalized to the housekeeping gene RPL32. Primer sequences are shown in Supplementary Table 4.
Analysis of Mayo Clinic, JHMI and TJU cohorts
Affymetrix Human Exon 1.0 ST arrays were used to analyze NRP1 expression in RP samples derived from Mayo Clinic (GSE46691), JHMI31 and TJU (GSE72291) repositories (545, 644 and 130 patients, respectively). Within the JHMI repository, 361 patients who were positive for BCR post-RP were analyzed. Patient clinical characteristics and sample preparation methods have been described previously.30, 31 Data were normalized and summarized using the SCAN algorithm. For grouped analysis, samples were split by the median expression of NRP1 into groups of low and high expression. The prognostic value of NRP1 was evaluated using multivariable odds ratios, and area under the receiver operating characteristics curve for BCR, metastasis and PCSM endpoints. Metastatic progression was defined by a positive bone or computed tomography scan. Kaplan–Meier survival analysis curves were generated for the JHMI cohorts but not for the Mayo cohort because of its nested case–control study design. All studies analyzed adhered to the PRoBE and REMARK guidelines for blinded evaluation and analysis of prognostic biomarkers.47, 48
Zebrafish metastasis assays
Research was carried in accordance with protocols compliant with the Canadian Council on Animal Care and with the approval of the Animal Care Committee at the University of British Columbia. Wild-type zebrafish were maintained in aquaria according to standard protocols.49 Embryos were generated by natural pair-wise matings and raised at 28.5 °C on a 14-h light/10-h dark cycle in a 100 mm2 Petri dish containing aquarium water. In all, 0.2 mM phenylthiourea was added to the embryos at 10 h after fertilization to prevent pigment formation.50
PC3 cells were fluorescently labeled 24 h before microinjection with 1.5 μM CellTracker CM-Dil dye (Thermo Fisher Scientific) according to the manufacturer’s instructions. Wild-type embryos were dechorionated at 2 days after fertilization. Following tricane anesthetization, 50–70 cancer cells were microinjected into the yolk sacs of 50 animals per treatment group. Embryos were transferred to 100 mm2 dishes that contained aquaria water with added phenylthiourea. Embryos were visually assessed for presence of xenografts and only successfully xenografted embryos were included in the experiment. Embryos were kept at 35 °C for the duration of the experiment. Metastatic dissemination outside the yolk sac was assessed 5 days later by observation using the Zeiss (Oberkochen, Germany) Axio Observer microscope with Zen 2012 software (Carl Zeiss Microscopy GmbH, Jena, Germany). Fixed cells were used as a control to ensure that observed metastasis was not due to yolk sac absorption.
Flow cytometry
Cells were grown to 80% confluence, washed twice in phosphate-buffered saline and detached using non-enzymatic cell dissociation buffer (Sigma Aldrich, St Louis, MO, USA). After washing and resuspension (106 cells/ml) in phosphate-buffered saline/5% FBS, 100 μl of cell suspension was incubated with anti-human NRP1-APC (#446921) or mouse IgG2a Isotype Control-APC (#20102) (R&D Systems, Minneapolis, MN, USA) for 1 h on ice, then washed three times in phosphate-buffered saline/5% FBS. Propidium iodide (3 μl, 100 μg/ml) was added to the cells immediately before loading on a FACS Canto (BD Bioscience, Franklin Lakes, NJ, USA) to allow for viable cell gating. Data analysis was performed using Kaluza (Beckman Coulter, Brea, CA, USA) or Flowjo (Flowjo LLC, Ashland, OR, USA) software.
Western blotting and immunofluorescence
Western blotting was performed as described previously.45 Primary antibodies used were NRP1 (#sc-7239, Santa Cruz, Dallas, TX, USA), AR (#D6F11, Cell Signaling Technology, Danvers, MA, USA) and GAPDH (#14C10, Cell Signaling Technology).
Invasion and proliferation assays
3D Laminin-rich Extracellular Matrix (lrECM) On-Top Cultures (referred to as 3D On-top Matrigel assays) were performed with an initial seeding density of 1x103 PC3 cells per well of a 96-well plate and conducted as previously.45 For invasion assays, 15 000 PC3 or 25 000 DU145 cells were seeded overnight into Matrigel-coated (100 μg/ml in growth media) wells in a 96-well Image-lock plate (Essen BioScience Inc., Ann Arbor, MI, USA). Wounds were made through the monolayer of confluent cells using the 96-pin WoundMaker (Essen BioScience Inc.) according to the manufacturer’s instructions. Wells were washed twice with phosphate-buffered saline and matrigel (50 μl, 1 mg/ml in growth media) was added to each well and allowed to solidify before the initiation of imaging. Images were captured every 2 h for up to 48 h by the IncuCyte FLR live cell imaging system (Essen BioSciences Inc.). Wound closure kinetics were determined using the CellPlayer software module (Essen BioScience Inc.). For proliferation assays, cells were seeded as described for invasion assays and cell confluence determined using the same software. Proliferation was also assessed using the Quant-iT PicoGreen dsDNA Assay Kit (Thermo Fisher Scientific).
Statistical analysis
Data analysis was performed by one-way analysis of variance with Tukey’s post hoc test for multiple comparisons, unless otherwise stated. Zebrafish metastasis assay results were analyzed by chi-square test. Statistical significance was defined as P<0.05. For publicly available microarray expression data sets, the normalized expression data for NRP1 was downloaded from the Oncomine database. Kaplan–Meier survival curve and log rank tests were performed using GraphPad Prism v6 software (GraphPad Software, Inc., San Diego, CA, USA).
References
Heidenreich A, Bastian PJ, Bellmunt J, Bolla M, Joniau S, van der Kwast T et al. EAU guidelines on prostate cancer. Part II: treatment of advanced, relapsing, and castration-resistant prostate cancer. Eur Urol 2014; 65: 467–479.
de Bono JS, Oudard S, Ozguroglu M, Hansen S, Machiels JP, Kocak I et al. Prednisone plus cabazitaxel or mitoxantrone for metastatic castration-resistant prostate cancer progressing after docetaxel treatment: a randomised open-label trial. Lancet 2010; 376: 1147–1154.
Tannock IF, de Wit R, Berry WR, Horti J, Pluzanska A, Chi KN et al. Docetaxel plus prednisone or mitoxantrone plus prednisone for advanced prostate cancer. N Engl J Med 2004; 351: 1502–1512.
Sweeney CJ, Chen YH, Carducci M, Liu G, Jarrard DF, Eisenberger M et al. Chemohormonal therapy in metastatic hormone-sensitive prostate cancer. N Engl J Med 2015; 373: 737–746.
James ND, Sydes MR, Clarke NW, Mason MD, Dearnaley DP, Spears MR et al. Addition of docetaxel, zoledronic acid, or both to first-line long-term hormone therapy in prostate cancer (STAMPEDE): survival results from an adaptive, multiarm, multistage, platform randomised controlled trial. Lancet 2016; 387: 1163–1177.
Karantanos T, Corn PG, Thompson TC . Prostate cancer progression after androgen deprivation therapy: mechanisms of castrate resistance and novel therapeutic approaches. Oncogene 2013; 32: 5501–5511.
Locke JA, Guns ES, Lubik AA, Adomat HH, Hendy SC, Wood CA et al. Androgen levels increase by intratumoral de novo steroidogenesis during progression of castration-resistant prostate cancer. Cancer Res 2008; 68: 6407–6415.
Yuan X, Cai C, Chen S, Chen S, Yu Z, Balk SP . Androgen receptor functions in castration-resistant prostate cancer and mechanisms of resistance to new agents targeting the androgen axis. Oncogene 2014; 33: 2815–2825.
Qi J, Tripathi M, Mishra R, Sahgal N, Fazli L, Ettinger S et al. The E3 ubiquitin ligase Siah2 contributes to castration-resistant prostate cancer by regulation of androgen receptor transcriptional activity. Cancer Cell 2013; 23: 332–346.
Ettinger SL, Sobel R, Whitmore TG, Akbari M, Bradley DR, Gleave ME et al. Dysregulation of sterol response element-binding proteins and downstream effectors in prostate cancer during progression to androgen independence. Cancer Res 2004; 64: 2212–2221.
Grasso CS, Wu YM, Robinson DR, Cao X, Dhanasekaran SM, Khan AP et al. The mutational landscape of lethal castration-resistant prostate cancer. Nature 2012; 487: 239–243.
Chandran UR, Ma C, Dhir R, Bisceglia M, Lyons-Weiler M, Liang W et al. Gene expression profiles of prostate cancer reveal involvement of multiple molecular pathways in the metastatic process. BMC Cancer 2007; 7: 64.
Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS et al. Integrative genomic profiling of human prostate cancer. Cancer Cell 2010; 18: 11–22.
Vanaja DK, Cheville JC, Iturria SJ, Young CY . Transcriptional silencing of zinc finger protein 185 identified by expression profiling is associated with prostate cancer progression. Cancer Res 2003; 63: 3877–3882.
Varambally S, Yu J, Laxman B, Rhodes DR, Mehra R, Tomlins SA et al. Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression. Cancer Cell 2005; 8: 393–406.
Hansel DE, Wilentz RE, Yeo CJ, Schulick RD, Montgomery E, Maitra A . Expression of neuropilin-1 in high-grade dysplasia, invasive cancer, and metastases of the human gastrointestinal tract. Am J Surg Pathol 2004; 28: 347–356.
Hong TM, Chen YL, Wu YY, Yuan A, Chao YC, Chung YC et al. Targeting neuropilin 1 as an antitumor strategy in lung cancer. Clin Cancer Res 2007; 13: 4759–4768.
Latil A, Bieche I, Pesche S, Valeri A, Fournier G, Cussenot O et al. VEGF overexpression in clinically localized prostate tumors and neuropilin-1 overexpression in metastatic forms. Int J Cancer 2000; 89: 167–171.
Miao HQ, Lee P, Lin H, Soker S, Klagsbrun M . Neuropilin-1 expression by tumor cells promotes tumor angiogenesis and progression. FASEB J 2000; 14: 2532–2539.
Stephenson JM, Banerjee S, Saxena NK, Cherian R, Banerjee SK . Neuropilin-1 is differentially expressed in myoepithelial cells and vascular smooth muscle cells in preneoplastic and neoplastic human breast: a possible marker for the progression of breast cancer. Int J Cancer 2002; 101: 409–414.
Zhang S, Zhau HE, Osunkoya AO, Iqbal S, Yang X, Fan S et al. Vascular endothelial growth factor regulates myeloid cell leukemia-1 expression through neuropilin-1-dependent activation of c-MET signaling in human prostate cancer cells. Mol Cancer 2010; 9: 9.
Pomerantz MM, Li F, Takeda DY, Lenci R, Chonkar A, Chabot M et al. The androgen receptor cistrome is extensively reprogrammed in human prostate tumorigenesis. Nat Genet 2015; 47: 1346–1351.
Chen M, Feuerstein MA, Levina E, Baghel PS, Carkner RD, Tanner MJ et al. Hedgehog/Gli supports androgen signaling in androgen deprived and androgen independent prostate cancer cells. Mol Cancer 2010; 9: 89.
Best CJ, Gillespie JW, Yi Y, Chandramouli GV, Perlmutter MA, Gathright Y et al. Molecular alterations in primary prostate cancer after androgen ablation therapy. Clin Cancer Res 2005; 11 (19 Pt 1): 6823–6834.
Robinson D, Van Allen EM, Wu YM, Schultz N, Lonigro RJ, Mosquera JM et al. Integrative clinical genomics of advanced prostate cancer. Cell 2015; 161: 1215–1228.
Peng Y, Liu YM, Li LC, Wang LL, Wu XL . MicroRNA-338 inhibits growth, invasion and metastasis of gastric cancer by targeting NRP1 expression. PLoS One 2014; 9: e94422.
Mak P, Leav I, Pursell B, Bae D, Yang X, Taglienti CA et al. ERbeta impedes prostate cancer EMT by destabilizing HIF-1alpha and inhibiting VEGF-mediated snail nuclear localization: implications for Gleason grading. Cancer Cell 2010; 17: 319–332.
Tsai JH, Yang J . Epithelial-mesenchymal plasticity in carcinoma metastasis. Genes Dev 2013; 27: 2192–2206.
Den RB, Feng FY, Showalter TN, Mishra MV, Trabulsi EJ, Lallas CD et al. Genomic prostate cancer classifier predicts biochemical failure and metastases in patients after postoperative radiation therapy. Int J Radiat Oncol Biol Phys 2014; 89: 1038–1046.
Erho N, Crisan A, Vergara IA, Mitra AP, Ghadessi M, Buerki C et al. Discovery and validation of a prostate cancer genomic classifier that predicts early metastasis following radical prostatectomy. PLoS One 2013; 8: e66855.
Ross AE, Johnson MH, Yousefi K, Davicioni E, Netto GJ, Marchionni L et al. Tissue-based genomics augments post-prostatectomy risk stratification in a natural history cohort of intermediate- and high-risk men. Eur Urol 2015.
Nouri M, Ratther E, Stylianou N, Nelson CC, Hollier BG, Williams ED . Androgen-targeted therapy-induced epithelial mesenchymal plasticity and neuroendocrine transdifferentiation in prostate cancer: an opportunity for intervention. Front Oncol 2014; 4: 370.
Yuan X, Balk SP . Mechanisms mediating androgen receptor reactivation after castration. Urol Oncol 2009; 27: 36–41.
Sprenger CC, Plymate SR . The link between androgen receptor splice variants and castration-resistant prostate cancer. Horm Cancer 2014; 5: 207–217.
Lu J, Lonergan PE, Nacusi LP, Wang L, Schmidt LJ, Sun Z et al. The cistrome and gene signature of androgen receptor splice variants in castration resistant prostate cancer cells. J Urol 2015; 193: 690–698.
Teng Y, Xie X, Walker S, White DT, Mumm JS, Cowell JK . Evaluating human cancer cell metastasis in zebrafish. BMC Cancer 2013; 13: 453.
Patnaik A, LoRusso PM, Messersmith WA, Papadopoulos KP, Gore L, Beeram M et al. A phase Ib study evaluating MNRP1685A, a fully human anti-NRP1 monoclonal antibody, in combination with bevacizumab and paclitaxel in patients with advanced solid tumors. Cancer Chemother Pharmacol 2014; 73: 951–960.
Glinka Y, Stoilova S, Mohammed N, Prud'homme GJ . Neuropilin-1 exerts co-receptor function for TGF-beta-1 on the membrane of cancer cells and enhances responses to both latent and active TGF-beta. Carcinogenesis 2011; 32: 613–621.
Patel P, West-Mays J, Kolb M, Rodrigues JC, Hoff CM, Margetts PJ . Platelet derived growth factor B and epithelial mesenchymal transition of peritoneal mesothelial cells. Matrix Biol 2010; 29: 97–106.
West DC, Rees CG, Duchesne L, Patey SJ, Terry CJ, Turnbull JE et al. Interactions of multiple heparin binding growth factors with neuropilin-1 and potentiation of the activity of fibroblast growth factor-2. J Biol Chem 2005; 280: 13457–13464.
Jia H, Cheng L, Tickner M, Bagherzadeh A, Selwood D, Zachary I . Neuropilin-1 antagonism in human carcinoma cells inhibits migration and enhances chemosensitivity. Br J Cancer 2010; 102: 541–552.
Soker S, Miao HQ, Nomi M, Takashima S, Klagsbrun M . VEGF165 mediates formation of complexes containing VEGFR-2 and neuropilin-1 that enhance VEGF165-receptor binding. J Cell Biochem 2002; 85: 357–368.
Pan Q, Chanthery Y, Liang WC, Stawicki S, Mak J, Rathore N et al. Blocking neuropilin-1 function has an additive effect with anti-VEGF to inhibit tumor growth. Cancer Cell 2007; 11: 53–67.
Snuderl M, Batista A, Kirkpatrick ND, Ruiz de Almodovar C, Riedemann L, Walsh EC et al. Targeting placental growth factor/neuropilin 1 pathway inhibits growth and spread of medulloblastoma. Cell 2013; 152: 1065–1076.
Hollier BG, Tinnirello AA, Werden SJ, Evans KW, Taube JH, Sarkar TR et al. FOXC2 expression links epithelial-mesenchymal transition and stem cell properties in breast cancer. Cancer Res 2013; 73: 1981–1992.
Lin D, Dong X, Wang K, Wyatt AW, Crea F, Xue H et al. Identification of DEK as a potential therapeutic target for neuroendocrine prostate cancer. Oncotarget 2015; 6: 1806–1820.
McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM et al. REporting recommendations for tumour MARKer prognostic studies (REMARK). Eur J Cancer 2005; 41: 1690–1696.
Pepe MS, Feng Z, Janes H, Bossuyt PM, Potter JD . Pivotal evaluation of the accuracy of a biomarker used for classification or prediction: standards for study design. J Natl Cancer Inst 2008; 100: 1432–1438.
Westerfield M . The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Danio Rerio). University of Oregon Press: Eugene, OR, USA, 1995.
Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF . Stages of embryonic development of the zebrafish. Dev Dyn 1995; 203: 253–310.
Sun Y, Wang BE, Leong KG, Yue P, Li L, Jhunjhunwala S et al. Androgen deprivation causes epithelial-mesenchymal transition in the prostate: implications for androgen-deprivation therapy. Cancer Res 2012; 72: 527–536.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov 2012; 2: 401–404.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Science Signal 2013; 6: p11.
Acknowledgements
We thank Dr Leland Chung and Dr Simon Hayward for kindly providing the C4-2B and BPH-1 cell lines, respectively, used in this study. This research was supported by Cure Cancer Australia Foundation, Prostate Cancer Foundation of Australia and Cancer Australia PdCCRS grant 108878 and the Australian Government Department of Health, Queensland Government National and International Research Alliance Program funding to the Australian-Canadian Prostate Cancer Research Alliance, as well as by the Movember Foundation and the Prostate Cancer Foundation of Australia through a Movember Revolutionary Team Award.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing interests
The authors declare no conflict of interest.
Additional information
Supplementary Information accompanies this paper on the Oncogene website
Supplementary information
Rights and permissions
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/4.0/
About this article
Cite this article
Tse, B., Volpert, M., Ratther, E. et al. Neuropilin-1 is upregulated in the adaptive response of prostate tumors to androgen-targeted therapies and is prognostic of metastatic progression and patient mortality. Oncogene 36, 3417–3427 (2017). https://doi.org/10.1038/onc.2016.482
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/onc.2016.482
This article is cited by
sBioSITe enables sensitive identification of the cell surface proteome through direct enrichment of biotinylated peptides
Clinical Proteomics (2023)
NRP1 promotes prostate cancer progression via modulating EGFR-dependent AKT pathway activation
Cell Death & Disease (2023)
The Neuropilin-1/PKC axis promotes neuroendocrine differentiation and drug resistance of prostate cancer
British Journal of Cancer (2023)
Prostate cancer genetic risk and associated aggressive disease in men of African ancestry
Nature Communications (2023)
The novel transcriptomic signature of angiogenesis predicts clinical outcome, tumor microenvironment and treatment response for prostate adenocarcinoma
Molecular Medicine (2022)