Next Article in Journal
Dietary Leucine Supplement Ameliorates Hepatic Steatosis and Diabetic Nephropathy in db/db Mice
Next Article in Special Issue
BRCA1 Interacting Protein COBRA1 Facilitates Adaptation to Castrate-Resistant Growth Conditions
Previous Article in Journal
Single Nucleotide Polymorphisms in the Vitamin D Receptor Gene (VDR) May Have an Impact on Acute Pancreatitis (AP) Development: A Prospective Study in Populations of AP Patients and Alcohol-Abuse Controls
Previous Article in Special Issue
Expression of Phospho-ELK1 and Its Prognostic Significance in Urothelial Carcinoma of the Upper Urinary Tract
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 7
10.3390/ijms19071920
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A hnRNP K–AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer

by
Matteo Capaia
1,†,
Ilaria Granata
2,†,
Mario Guarracino
2,
Andrea Petretto
3,
Elvira Inglese
3,
Carlo Cattrini
1,4,
Nicoletta Ferrari
5,
Francesco Boccardo
1,4 and
Paola Barboro
1,*
1
Academic Unit of Medical Oncology, Ospedale Policlinico San Martino-IRCCS, L.go R. Benzi 10, 16132 Genova, Italy
2
Institute for High Performance Computing and Networking (ICAR), National Research Council (CNR), Via Pietro Castellino 111, 80131 Napoli, Italy
3
Core Facilities-Proteomics Laboratory, Giannina Gaslini Institute, L.go G. Gaslini 5, 16147 Genova, Italy
4
Department of Internal Medicine and Medical Specialties, School of Medicine, University of Genova, L.go R. Benzi 10, 16132 Genova, Italy
5
Molecular Oncology and Angiogenesis, Ospedale Policlinico San Martino-IRCCS, L.go R. Benzi 10, 16132 Genova, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2018, 19(7), 1920; https://doi.org/10.3390/ijms19071920
Submission received: 15 June 2018 / Revised: 27 June 2018 / Accepted: 29 June 2018 / Published: 30 June 2018
(This article belongs to the Special Issue Sex Hormone Receptor Signals in Human Malignancies)
Browse Figures
Versions Notes

Abstract

:
The major challenge in castration-resistant prostate cancer (CRPC) remains the ability to predict the clinical responses to improve patient selection for appropriate treatments. The finding that androgen deprivation therapy (ADT) induces alterations in the androgen receptor (AR) transcriptional program by AR coregulators activity in a context-dependent manner, offers the opportunity for identifying signatures discriminating different clinical states of prostate cancer (PCa) progression. Gel electrophoretic analyses combined with western blot showed that, in androgen-dependent PCa and CRPC in vitro models, the subcellular distribution of spliced and serine-phosphorylated heterogeneous nuclear ribonucleoprotein K (hnRNP K) isoforms can be associated with different AR activities. Using mass spectrometry and bioinformatic analyses, we showed that the protein sets of androgen-dependent (LNCaP) and ADT-resistant cell lines (PDB and MDB) co-immunoprecipitated with hnRNP K varied depending on the cell type, unravelling a dynamic relationship between hnRNP K and AR during PCa progression to CRPC. By comparing the interactome of LNCaP, PDB, and MDB cell lines, we identified 51 proteins differentially interacting with hnRNP K, among which KLK3, SORD, SPON2, IMPDH2, ACTN4, ATP1B1, HSPB1, and KHDRBS1 were associated with AR and differentially expressed in normal and tumor human prostate tissues. This hnRNP K–AR-related signature, associated with androgen sensitivity and PCa progression, may help clinicians to better manage patients with CRPC.

Graphical Abstract

1. Introduction

In the last years, experimental evidence has supported the role of the androgen receptor (AR) in castration-resistant prostate cancer (CRPC) development. Almost all patients with metastatic prostate cancer (PCa) initially treated with androgen deprivation therapy (ADT) progress to CRPC, in many cases following the reactivation of the AR pathway. Several mechanisms, which are not necessarily mutually exclusive, have been proposed to explain CRPC development [1]. They include AR amplification or overexpression, AR mutations that can modify the ligand specificity, AR gain of function, changes in the expression levels of AR coregulators, and involvement of alternative pathways that could be completely independent of AR signaling [2,3]. However, to date, the molecular mechanisms through which hormone-sensitive PCa cells acquire the ability to resist to hormone deprivation need further investigation.
Multi-omics studies showed that CRPC is a heterogeneous group of diseases characterized by different genotypes and phenotypes [4,5,6]. Emerging evidence suggests that phenotypic plasticity, driving the adaptation to ADT stress, further amplifies cellular heterogeneity and can contribute to ADT resistance [7,8]. These findings indicate that an adequate response to cytotoxic or targeted therapies cannot disregard the broad spectrum of cellular sub-clones with different clinical behaviors. Since multiple pathways are involved in CRPC development and progression, it is evident that therapies, to be useful, should selectively target driving molecular alterations at a specific stage of PCa evolution. For this purpose, we explored the possibility of developing a signature for identifying PCa and CRPC subtypes with different androgen responsiveness.
A new scenario to develop alternative targeted therapies was recently proposed by Liu et al. [9], reporting the context dependency by which AR coregulators control selectively an AR target gene set reflecting PCa biology and evolution. This finding also provides a proof of principle for the identification of PCa subtypes associated with an AR coregulator and the corresponding subset of AR-related genes.
The heterogeneous nuclear ribonucleoprotein K (hnRNP K) is a multifunctional protein playing a pivotal role in regulating numerous cellular functions, such as transcription, signal transduction, alternative splicing, and chromatin remodeling [10]. This functional versatility depends on post-translation modifications (PTM) that modulate its interactions with nucleic acids and proteins [11]. Increasing evidence for the involvement of hnRNP K in cancer progression was reported [12]. In PCa patients, we have demonstrated that its overexpression positively correlates with Gleason score and poor patients prognosis [13] and that the concomitant expression of both AR and cytoplasmic hnRNP K has a potential prognostic value [14]. In PCa cell lines, hnRNP K regulates AR activity by inhibiting its translation [15], and, in nucleoplasm, hnRNP K phosphorylation shapes the AR–DNA complex after anti-androgen treatments [16]. HnRNP K appears to be also able to regulate neuroendocrine differentiation [17].
Kelly et al. [18] underlined the importance of the adaptive phenotype acquired during ADT leading to cellular reprogramming that ultimately resulted in tumor heterogeneity and different AR status. To reproduce this behavior, we obtained the androgen-resistant cell lines PDB and MDB by treating LNCaP cell line for over one year with the anti-androgen bicalutamide (BIC) in the presence or absence of 5-α-dihydrotestosterone (DHT), respectively. Transcriptomic and proteomic analyses highlighted the high degree of phenotypic plasticity that characterizes our CRPC models and allows their adaptation under stress conditions. BIC-resistant cell lines represent two sub-populations with AR levels and subcellular localization similar to the parental LNCaP cell line, but with reduced functionality depending on AR phosphorylation status. Interestingly, partial (PDB) or minimal (MDB) AR transcriptional activity correlated with enhanced tumorigenicity and decreased sensitivity to treatment with a novel anti-androgen, enzalutamide, compared to the parental cell line [19].
Given the above findings, in this study, we investigated the role of hnRNP K in androgen-resistance using in vitro models of androgen-dependent or castration-resistant PCa. We hypothesized that hnRNP K, working like an AR transcriptional collaborator, could participate in regulating the different AR transcriptional programs during PCa development and progression. Consequently, a hnRNP K signature associated to AR activity could be useful to identify clinically distinct PCa subgroups.

2. Results

2.1. Role of hnRNP K in ADT Resistance

We previously reported that in the androgen-dependent cell line LNCaP, changes in the AR binding property of hnRNP K was associated with cell growth and AR activity [16,20]. Here, using our in vitro resistant models, we investigated the role of hnRNP K in the resistant cell lines PDB and MDB which can be considered models mimicking two CRPC subpopulations [19].
HnRNP K silencing decreased both AR and PSA expression in LNCaP cells (Figure 1a,b), while it was less effective in BIC-resistant cells lines, in particular MDB. Furthermore, hnRNP K silencing in LNCaP induced a 68% reduction of AR activity, as evaluated by the luciferase assay (Figure 1c). Because of their reduced AR activity [19], it was not possible to determine luciferase activity in PDB and MDB cells.
BIC removal and restored availability of androgen in the culture media determined comparable features in prostate-specific antigen (PSA) expression in both our ADT-resistant cell lines, with a maximum effect after 3–4 weeks (Figure 1d,e). However, in MDB cells, the substantial increase in PSA synthesis was associated with the overexpression of both hnRNP K and AR, while, in PDB cells, it was independent of it, probably due to AR hypersensitivity developed during ADT.
Overall, these results suggest that hnRNP K in prostate cancer cell lines is likely to act as an AR transcriptional collaborator that regulates AR activity through different molecular mechanisms depending on cell differentiation, supporting the hypothesis of the involvement of hnRNP K in ADT resistance.

2.2. Role of hnRNP K Phosphorylation in ADT Resistance

The mechanistic role of phosphorylation in regulating hnRNP K transcriptional activity [21,22] as well as its increased expression in several neoplasms with an aggressive phenotype [12] have been described. Here, we evaluated both its expression level and phosphorylation status in LNCaP and BIC-resistant PDB and MDB cell lines.
Unexpectedly, using quantitative western blotting (WB) analysis (Figure 2a), we detected a significant hnRNP K decrease in PBD and MDB cells compared to the parental cell line LNCaP, suggesting its distinct functional role in androgen-resistance compared to PCa where hnRNP K overexpression has shown diagnostic and prognostic value [13].
By evaluating the total hnRNP K phosphorylation status by means of the monodimensional phosphate affinity gel electrophoresis (1D Phos-tag), it was possible to identify a non-phosphorylated isoform (phK0) and two species characterized by intermediate (phK1) and maximal (phK2) phosphorylation (Figure 2b). As shown in the histogram of Figure 2b (right panel), significant differences were found only for PDB, showing phK0 decrease, compared to LNCaP and MDB, and phK1 increase with respect to MDB.
Phosphorylation and cellular compartmentalization of hnRNP K isoforms generate a regulatory system involved in cell growth [23] and translation regulation [24]. Aberrant hnRNP K hyperphosphorylation and cytoplasmic accumulation are peculiar features of several human tumors, often associated with a worse prognosis [12]. As serine residues phosphorylation and hnRNP K functions involved in regulating its intracellular distribution, cell growth, and transcription (PhosphoSitePlus database: www.phosphosite.org) are closely related, we evaluated the role of the subcellular distribution of hnRNP K phosphorylated isoforms in PCa evolution. Using Phos-tag bidimensional gel electrophoresis (2D Phos-tag) and WB (Figure 3), we analyzed nuclear and cytoplasmic extracts from the androgen-dependent PCa cell line LNCaP and from cell lines with different CRPC phenotypes with respect to AR status and transcriptional activity: three AR-positive cell lines (PDB, hypersensitive; MDB, inactive; 22Rv1, ARv7, androgen-independent) and the AR-negative PC3 cell line [25]. Each spot detected with the anti-hnRNP K antibody was attributed, according to Kimura el al. [26], to the alternatively spliced isoforms 1 and 2 and to the phosphorylated isoforms at Ser116, Ser284, Ser353 residues (pS116, pS284, pS353), as schematically represented in Figure S1. Quantitative analysis of each spot was carried out using PDQuest software and is reported in Table 1 and Figure S2.
The level of nuclear hnRNP K phosphorylated isoform 1 was higher than that of isoform 2 in resistant cell lines compared to androgen-dependent LNCaP, while, in the cytoplasm, alternatively spliced isoforms phosphorylation was regulated in an opposite manner. Interestingly, pS353, localized in the nuclear shuttling (KNS) domains regulating hnRNP K intracellular localization, was detected exclusively in isoform 2 (Figure S2).
Quantitative analysis of nuclear hnRNP K phosphorylated forms in all cell lines revealed that: (i) pS353 decreased in resistant cell lines compared to LNCaP; (ii) an over 36% pS284 increase was observed in cell lines with active AR axis, while pS284 decreased to about 23% in MDB and PC3 cells; (iii) both high percentages of non-phosphorylated forms and low pS116 percentages were detected in androgen-independent cell lines MDB (AR-positive), 22Rv1 (ARv7), and PC3 (AR-negative), independently of the AR status.
The cytoplasmic hnRNP K isoforms distribution in the different cell lines showed that pS116 could discriminate between LNCaP and CRPC lines, regardless of AR functionality. Cytoplasmic pS353 showed an opposite trend in PDB and 22Rv1 cell lines with aberrant AR activation, in comparison to LNCaP. The pS284 isoforms were overrepresented in androgen-independent MDB and 22Rv1 cell lines.
These findings suggest that differential phosphorylation at specific serine residues and compartmentalization of hnRNP K splicing isoforms 1 and 2 correlate with different resistant phenotypes, providing further evidence for hnRNP K involvement in CRPC evolution.

2.3. Characterization of the hnRNP K Interactome in LNCaP, PDB, and MDB Cell Lines

HnRNP K may regulate several cellular functions, such as transcription and signal transduction, by PTMs that modify its binding partners [11,27]. Using a co-immunoprecipitation assay coupled with mass spectrometry (MS), we identified the proteins directly or indirectly interacting with hnRNP K in LNCaP, PDB, and MDB total extracts. A total of 254 proteins were identified (Table S1).
As shown in Figure 4, the number of proteins interacting with hnRNP K gradually decreased from 221 in LNCaP, to 111 in PDB, and 55 in MDB cell lines. Protein–protein interaction networks for LNCaP, PDB, and MDB cell lines were created through Intact and Reactome databases and then clusterized through the Cytoscape app clustermaker2 [28]. Three significant clusters were obtained for each cell line, and the proteins belonging to each of them were enriched through JEPETTO Cytoscape plug-in [29], according to the Gene Ontology (GO) Molecular Function annotation restricted to prostate tissue. In cluster 1, hnRNP K was the primary hub connecting about 50% of proteins both in LNCaP and in PDB cells, while cluster 2 was predominant in the resistant cell line MDB. Interestingly, some proteins involved in AR binding (Figure 5, red arrow) were significantly clustered only in the androgen-dependent cell line LNCaP. In all cell lines, enriched GO terms in cluster 1 and 2 were different, while cluster 3, referred mainly to mRNA binding, represented about 20% of hnRNP K interacting proteins. Since, in our experiments, MS analysis was not able to efficiently detect AR, we verified hnRNP K–AR association in all cell lines using WB analysis of hnRNP K co-immunoprecipitates (co-IPs) obtained from LNCaP, PDB, and MDB (Figure S3).
By comparing hnRNP K co-IP abundance values obtained from LNCaP, PDB, and MDB, we identified three sets of differentially interacting proteins. The data were linearly modeled using Limma R package, and from the three contrasts generated (PDBvsLNCaP, MDBvsLNCaP, MDBvsPDB) we identified 51 proteins differentially interacting (p < 0.05) with hnRNP K (DIhKP), as reported in Table S2. Decreased interaction with hnRNP K was observed in all comparisons, except for PDBvsLNCaP, where eight proteins increased their binding capacity.
Enrichment analysis was performed to identify the biological processes and molecular functions associated to DIhKPs, according to Gene Ontology annotation (Table S3). Among the significant terms, the highest percentages corresponded to proteins involved in translation (17% for GO Biological Processes) and RNA binding (35% for GO Molecular Function), indicating that in PCa resistance, the central role of hnRNP K is linked to protein synthesis regulation.
In PCa and CRPC tissues, the altered profile of AR-interacting transcription factors (TFs) determines AR transcriptional program modifications [30]. By performing the enrichment of transcription factor binding sites (TFBS) through Database for Annotation Visualization and Integrated Discovery, we identified six TFs potentially regulating the expression of DIhKPs (hK-TFs) (Table S3). It is significant that five of these hK-TFs revealed a direct association with AR and PCa. ZIC2 and HOXB13 are AR-controlled genes [9]; the latter, in the absence of androgen, promotes androgen-independent growth in PCa cell line [31], while HOXA13 overexpression is significantly associated to PCa poor prognosis [32]. CMYB shares with AR common genes that correlate with advanced stages of PCa [33]. The TF XBP1, mediating cellular stress responses, is proposed as a predictor marker for PCa [34]. The bidirectional crosstalk between PPARG and AR regulates growth and development both in normal and in tumor prostate, as recently highlighted [35]. Using Cytoscape, we generated a network visualization of the predicted regulatory interaction between hK-TFs and 51 DIhKPs (Figure S4). In this network, PPARG and PAX5 were the main hK-TFs, regulating about 25% and 18% of DIhKPs, respectively. In addition, all DIhKPs and hK-TFs appeared to be potentially regulated by one or more hK-TF. Interestingly, hnRNP K expression resulted to be potentially regulated by four hK-TFs, among which PPARG, HOX13, and ZIC2 are controlled by AR.
As LNCaP, PDB and MDB cell lines can mimic distinct stages of PCa progression [19], we can infer that the Cytoscape network representation of the hK-TFs and the controlled DIhKPs grouped by comparing PDB versus LNCaP, MDB versus LNCaP, MDB and PDB versus LNCaP, and MDB versus PDB could represent proteins associated with different stages of PCa progression (Figure 5). Although PPARG represents the main hK-TF controlling DIhKPs in all networks, the global regulatory feature of hnRNP K partners is different in all comparisons. In PDB versus LNCaP, the eight proteins with increased hnRNP K interaction (Table S2) were under transcriptional control of PPARG and XBP1 (Figure 5a). Conversely, in MDB versus LNCaP, excluding PPARG, the number of DIhKPs controlled by others hK-TF was equally distributed (Figure 5b). Interestingly, despite the exiguous DIhKP number, the networks that featured the differences between androgen-dependent and -resistant cell lines or between MDB and PDB cell lines (Figure 5c,d) showed, respectively, an enhanced role of XBP1 and HOX13 in regulating hnRNP K partners transcription.
Overall, these findings indicate that, in cell lines with different androgen-resistant phenotypes, decreased hnRNP K expression and modified subcellular distribution of its isoforms were associated with dynamic remodeling of the hnRNP K interacting protein network depending on cellular differentiation, thus providing a rationale for DIhKPs analysis in human prostate cancer tissues.

2.4. Signature Identification for Potential PCa and CRPC Patient Stratification

From the above results, we hypothesized that some DIhKPs could differentiate distinct stages of PCa progression to CRPC, thus conferring clinical relevance to the findings obtained in our CRPC cell lines PDB and MDB. In this section, we describe the analytical strategies used to identify, among the 51 DIhKPs, a protein set that might putatively discriminate among clinical stages of PCa progression.
To exclude any potential bias due to the different expression of DIhKPs in the three cell lines, we verified the differential abundance calculated in our previous work [19] using total LNCaP, PDB, and MDB extracts by MS. For all proteins, the concordance between differential expression and altered hnRNP K binding capacity (respectively, differentially expressed protein—DEP—and DIhKP columns in Table S4) was evaluated, excluding 16 proteins (highlighted in light grey in Table S4) from further analysis. Moreover, referring to the expression determined by immunohistochemistry (IHC) reported in “The Human Protein Atlas” (http://www. proteinatlas.org/), we verified that only three proteins (highlighted in heavy grey in Table S4) out of the 51 DIhKPs were not expressed in normal and tumor prostate tissues. According to these findings, 32 DIhKPs were retained and considered for further evaluations.
By analyzing AR-related profiles reported in the literature evaluating AR signature [36], androgen-responsive genes [37], AR-controlled genes, and co-regulators [9,38], we verified that six out of the 32 DIhKPs retained for our analysis correlated with AR activity (Table S5).
Next, we evaluated differential expression, both at protein and transcript level, in human prostate tissues for the selected 32 DIhKPs (Table S5). Using gene expression data from four published datasets (expression profiling by array) available on Gene Expression Omnibus (GEO) portal (GSE32269, GSE101607, GSE29650, GSE66532), we performed the differential expression analyses to identify transcript alterations in tissues at different stages of PCa development. Interestingly, 53% of the 32 DIhKPs selected in our in vitro models were differentially expressed also in the first two datasets comparing non-tumoral tissue, PCa, and two subgroups of CRPC (AR-driven and non-AR-driven) bone metastases. By analyzing recent and exhaustive literature of proteomic studies in PCa tissues reported by Larkin et al. [39], we observed different expression of seven DIhKPs in PCa compared with non-tumoral or hyperplastic prostate tissues.
Combining all data reported in Table S5, we defined a hnRNP K-AR-related signature comprising eight DIhKPs selected on the basis of their AR association as coregulators and/or interactors and their differential expression in both mCRPC subgroups and PCa (Table 2).
For all proteins, evidence in the literature, reported in Table 2, shows a clear association with PCa progression. KLK3, SORD, and SPON2 expression differentiates hormone-dependent PCa from the two mCRPC subgroups. IMPDH2 tissue expression can differentiate PCa and mCRPC non-AR driven subgroups. Although there was no evidence for differential expression in androgen-dependent PCa and mCRPC tissues, the other four proteins (ACTN4, ATP1B1, HSPB1, KHDRBS1) were associated with AR. Of note, three AR co-regulators showed decreased hnRNP K interaction with respect to the parental cell line, depending on the cellular context (ACTN4 and HSPB1 for PDB and KHDRBS1 for MDB cell line). All the hK-TFs are potentially involved in regulating the transcription of the eight selected proteins. Exclusively the expression of ACTN4, an AR coregulator, is regulated by all hK-TFs, while the other proteins are controlled by one or more hK-TF.
Globally, these observations suggest that the hnRNP K–AR-related signature identified in our in vitro models might be associated to advanced stages of PCa and could identify molecular mCRPC subgroups.

3. Discussion

In the last years, experimental evidence has demonstrated that hnRNP K overexpression and aberrant cytoplasmic localization have an important role in PCa progression [12]. Recent findings suggest a relationship between AR and hnRNP K, as the latter regulates AR expression [57,58] and mRNA translation [15,24] and shapes a complex with AR and DNA modulated by anti-androgens [20]. In PCa tissues, the association of AR with cytoplasmic hnRNP K has a potential prognostic value [14], while the deregulation of the AKT–hnRNPK–AR–β-catenin pathway is involved in neuroendocrine differentiation [17]. Despite these evidence, the role of hnRNP K in androgen resistance has not been established yet.
HnRNP K is a modular protein that regulates a wide range of biological processes, among which transcription, splicing, and translation involved in the regulation of gene expression. By its ability to interact with both nucleic acids and proteins, hnRNP K acts like a docking platform coordinating the cross-talk between cell signaling pathways [10,21]. Consequently, the hnRNP K interactome is dynamic and changes depending on cellular compartments and external stimuli [27]. Moreover, hnRNP K functional versatility is regulated by several PTMs that modify its binding with nucleic acid and proteins [11,26].
On the basis of these premises, we explored the role of hnRNP K in androgen-resistance acquisition using three PCa cell lines: the androgen-responsive LNCaP and two resistant cell lines representing in vitro CRPC models, namely, the androgen-hypersensitive PDB and the androgen-insensitive MDB cell lines [19].
In LNCaP cell line, hnRNP K was found to regulate both AR expression and AR transcriptional activity (Figure 1a,b), while these functions were compromised in resistant cell lines, maximally in the most aggressive MDB. The adaptive phenotype characterizing PDB and MDB cells [19] under less stressful conditions (BIC removal) was found to drive a temporary recovery of AR activity which appeared to depend on the cellular context (Figure 1d,e). In the AR-inactive MDB cell line, restored androgen levels temporary reactivated the AR axis, determining a PSA rise correlated with AR and hnRNP K overexpression. Conversely, in PDB cell line, BIC removal increased PSA expression independently of AR and hnRNP K but associated with increased hnRNP K phosphorylation (Figure 2b) that could drive the AR hypersensitivity acquired during prolonged BIC exposure in androgen-containing medium.
The functional flexibility of hnRNP K arises, as for all the others hnRNPs, from the expression of alternatively spliced isoforms, PTMs, and subcellular distribution, which modify its binding activity to nucleic acids and proteins [59]. In response to mitogenic stimuli, the phosphorylated isoforms 1 and 2 generated by hnRNP K alternative splicing were differentially distributed in the nucleus and the cytoplasm [26], whereas the isoform 2 was downregulated in the colon carcinoma cell line HCT116 after p53 activation [60]. In the present study, we demonstrate a role for hnRNP K in ADT resistance acquisition depending on its alternative splicing, phosphorylation, and subcellular localization (Table 1 and Figure S2). Regardless of the AR status, the resistant cell lines showed higher levels of nuclear isoforms 1 and cytoplasmic isoform 2 than the androgen-responsive LNCaP. A quantitative comparison of the nuclear profiles of the pS116, pS284, and pS353 isoforms in the androgen-responsive LNCaP cell lines and in the resistant cell lines with different AR status and androgen responsiveness (Figure 3 and Table 1) showed a relationship between hnRNP K phosphorylation and AR transcriptional activity. In particular, to an active AR corresponded a low percentage of non-phosphorylated hnRNP K nuclear isoforms and vice versa (24% in LNCaP versus 66% in PC3). Interestingly, Ser116-phosphorylated isoforms could discriminate between resistant cell lines with different AR activity (nuclear fraction) and between LNCaP and all resistant cell lines (cytoplasmic fraction). It has been shown that the nuclear–cytoplasmic trafficking of hnRNP K may depend on Ser284 and Ser353 phosphorylation, determining inhibition of translation [24], and that Ser353 was involved in hnRNP K transactivation activity [61]. In line with these observations, our analysis showed an altered subcellular distribution of pS284 and pS353 in all CRPC cell lines in comparison to LNCaP.
Overall, these results suggest that the decreased phosphorylation of hnRNP K isoform 2, harboring the pS353 residue, and the subcellular distribution of pS116 and pS284 could have a significant role in ADT resistance.
It is increasingly evident that the identification of a molecular signature for specific stages of PCa would be useful not only for the selection of an appropriate and more effective therapy, but also for understanding the molecular basis of ADT resistance with relevant therapeutic implications [62]. The complex genomic and transcriptomic landscape of PCa has been explored through molecular analyses leading to the evidence that AR is involved in the development of drug resistance [61]. The selection pressure exerted by ADT induces alterations both in AR transcriptional program and in androgen-responsive gene (ARG) expression [37], triggering different downstream signaling pathways in PCa tissues that promote resistance [30,63]. Among the mechanisms proposed for AR reactivation during ADT, coregulators can play a critical role in influencing AR transcriptional activity and can represent new potential therapeutic targets for CRPC [64,65,66,67]. Recently, it was reported that AR coregulators control ARGs associated with different cellular processes in a context-dependent manner [9]. The dysregulated expression of AR coregulators in androgen-dependent PCa and in CRPC was reported to correlate with a poor prognosis and a more aggressive disease [65]. Therefore, it is reasonable to assume that the altered hnRNP K expression and function, by modifying the interaction with its partners involved in protein expression regulation (Figure 4 and Table S3), could alter the AR transcriptional program and androgen responsiveness in all stages of PCa.
In addition to AR coregulators, Obinata et al. [68] suggested a role for TFs in the control of the altered AR transcriptional program during anti-androgen treatment. The transcriptional regulatory networks of hK-TFs and DIhKPs subsets showed in Figure 5 confirm this hypothesis, since the composition of the hnRNP K interactome changed dynamically from LNCaP to the resistant cell lines PDB (AR-hypersensitive) and MDB (AR-inactive). Among the six hK-TFs, AR-controlled PPARG and HOX13 can be associated to ADT resistance. A recent study [35] has reported that PPARG controls prostate cell growth and differentiation and is involved in androgen-dependent PCa and CRPC evolution, acting as a tumor suppressor through its interplay with AR. The TF HOXA13, whose expression is associated with unfavorable survival [32], and the HOXB13, which controls PCa cell growth by inhibiting AR signaling [31], are also involved in PCa evolution.
The findings that emerge from our study demonstrate, for the first time, a clear-cut relationship between hnRNP K and AR, since, among the DIhKPs and the six hK-TFs, we found three AR coregulators (ACTN4, HSB1, KHDRBS1), four AR-controlled proteins (KLK3, SORD, SPON2, ATP1B1), and three TFs associated with AR activity (PPARG, ZIC2, HOX13) (Figure 6 and Table 2). In both androgen-resistant cell lines, the main hK-TF was PPARG which controls two distinct DIhKPs sets: in PDB (Figure 5a), the expression of eight upregulated DIhKPs was associated with the downregulation of the AR coregulators ACTN4 and HSB1, while, in MDB cell line (Figure 5b), all proteins under PPARG control were downregulated, except for PDIA4, including the AR coregulator KHDRBS1. Thus, we can infer that hnRNP K acts as a coordinator between AR and PPARG pathways and, depending on the AR coregulator involved, differently regulates the AR axis in androgen-dependent or -independent cell lines. As a confirmation, Moss et al. [69] reported different AR activity regulation by the PPARG ligand ciglitazone depending on androgen sensitivity, while Yang et al. [70] described that hnRNP K and KHDRBS1 interaction can regulate their activities in signal transduction pathways.
Understanding the context and the functional role of key proteins in ADT resistance could allow the identification of a molecular signature, rather than single proteins, that could overcome the main limit of prostate cancer treatments: tumor heterogeneity. Previously, we demonstrated that, in our CRPC model, the AR transcriptional program was modified by ADT [19]. In this study, we found that hnRNP K expression, alternatively spliced isoforms localization, and impaired phosphorylation promote ADT resistance depending on the cellular context and external stimuli. Using different analytical strategies, we identified eight proteins with dissimilar roles in the AR axis and correlation with clinical PCa progression. Of note, KLK3, SORD, and SPON2 can discriminate between non-AR-driven and AR-driven mCRPC and are characterized by common features: (i) decreased interaction with hnRNP K; (ii) control by AR; (iii) association with a single hK-TF (PAX5, PPARG, HOX13).
The eight proteins highlighted in this study (KLK3, SORD, SPON2, IMPDH2, ACTN4, ATP1B1, HSPB1, and KHDRBS1) represent a signature associated with AR modulation during the transition from androgen-dependent to castration-resistant prostate cancers. These proteins might discriminate specific stages of PCa and thus predict patients’ prognosis and response to treatments. Further investigations on PCa tissues or liquid biopsies from CRPC patients are needed to validate the clinical usefulness of the hnRNP K–AR-related signature and to evaluate whether some of the proteins included in our signature might become potential targets for new drugs development.

4. Materials and Methods

4.1. Cell Culture

The human prostate cancer cell lines LNCaP, 22Rv1, and PC3 were purchased from the American Type Culture Collection (CRL-1740, CRL-2505, CRL1435) and cultured in RPMI (Celbio, Milan, Italy) medium containing 10% fetal bovine serum, 2 mM glutamine, 1% penicillin, and 1% streptomycin. The LNCaP and 22Rv1 medium was also supplemented with 10 mM HEPES, 10 mM sodium pyruvate, and 0.25% glucose. The passage numbers at which LNCaP cells were used ranged from 22 to 30. The BIC-resistant cell lines PDB and MDB, obtained in our laboratory [19], were cultured in RPMI medium without phenol red, containing 10% charcoal-stripped fetal bovine serum, 2 mM glutamine, 10 mM HEPES, 10 mM sodium pyruvate, 0.25% glucose, 1% penicillin, 1% streptomycin, and 10 μM BIC. The PDB cell line was supplemented with 0.1 nM DHT.

4.2. Cell Fractionation

All buffers were supplemented with a protease inhibitors cocktail (5 mM Na2S2O5, 1 mM phenylmethylsulfonyl fluoride, 0.5 mM benzamidine, 20 mg/mL leupeptin, 10 mg/mL pepstatin A, 25 mg/mL aprotinin), 1 mM Na3VO4, and 1 mM dithiothreitol. The fractionation procedures were performed at 4 °C. The cells were mechanically harvested with a sterile cell scraper and washed three times in 30 mL of PBS. Total protein extracts, obtained using RIPA buffer, and cytoplasmic and nuclear extracts were prepared as reported [16,20]. After cell fractionation, the pellet was delipidated, cleaned, and solubilized for 1D or 2D gel electrophoresis as previously described [20]. The protein concentrations were determined using the Bio-Rad (Hercules, CA, USA) protein microassay with bovine serum albumin as a standard.

4.3. 1D and 2D Gel Electrophoresis and WB Analyses

Eight micrograms of proteins extracted from different cell fractions were loaded onto an 8–14% linear gradient (1D) or 7.5% (1D Phos-tag) polyacrylamide gels, containing 25 mM polyacrylamide -bound Mn2+-Phos-tagTM ligand (Phos-tag Consortium, Hiroshima, Japan) and separated at 5 mA/gel for 16 h at a constant temperature of 12 °C. 2D Phos-tag was performed as described [20]. After electrophoresis separation, the proteins were transferred to a Hybond-P membrane (GE Healthcare, Piscataway, NJ, USA) and probed with anti-hnRNP K (Santa Cruz Biotechnology Inc., Santa Cruz, CA, USA; 1:5000), anti-AR (DAKO, Carpinteria, CA, USA; 1:500), and anti-PSA (Cell Signaling Technology, Danvers, MA, USA; 1:2000) antibodies at 4 °C overnight. HRP-conjugated secondary antibodies (Cell Signaling Technologies, Danvers, MA, USA; 1:2000) were used, and protein bands were detected by chemiluminescent HRP substrate (Immobilon Western, Millipore, Darmstadt, Germany) and Hyper film-ECL (GE-healthcare, Lafayette, CO, USA), according to the manufacturer’s instructions. As PDB and MDB cell lines showed a global downregulation of protein expression compared to LNCaP cell line [19], the conventional normalization approach was inadequate. Different results were in fact obtained using β-actin- or Sypro Ruby- (Invitrogen, Eugene, OR, USA) normalized data to evaluate the relative amounts of AR and hnRNP K in total extracts from LNCaP, PDB, and MDB (Figure S5). For 1D analysis, the relative amount of immunoreactive bands on WB films scanned with a GS-800 densitometer (BioRad, Hercules, CA, USA), was assessed using Bio-Rad Quantity One software and normalizing the optical density of each band to the total optical density of the corresponding lane in the Sypro Ruby-stained gel providing the total protein content [16]. In Figure S6, the Sypro Ruby-stained gels used for quantitative analysis of WB reported in Figure 1a,d and Figure 2a are shown. The quantitative analysis of 1D Phos-tag WB was carried out comparing the relative percentage of phK0, phK1, and phK2 hnRNP K isoforms in each cell lines. The 2D Phos-tag analyses, including detection, alignment, and matching for spots present in all 2D maps, were done using the software package PDQuest (ver. 8.6, BioRad, Hercules, CA, USA).
Statistical significance was evaluated by two-tailed Student’s t-test within OriginPro 7.5 software (OriginLab, Northampton, MA, USA), and the differences were considered statistically significant when p < 0.05.

4.4. hnRNP K Silencing and Reporter Assay

ON-TARGET plus SMART pool for human hnRNP K (Dharmacon, GE- healthcare, Lafayette, CO, USA) was used to knockdown the expression of hnRNP K; siCONTROL non-targeting (NT) siRNA pool (Dharmacon, Lafayette, CO, USA) was used as a negative control. The cells were grown in their proper medium and transfected with Lipofectamine 2000 (Thermo Fisher, Waltham, MA, USA) following the manufacturer’s instructions. Protein expression was analyzed by WB 72 h after transfection.
To measure AR function, the cells were transfected with Lipofectamine 2000 utilizing the Cignal androgen receptor dual-luciferase reporter kit (Qiagen, Hilden, Germany), following the manufacturer’s instructions. Luciferase activity was assayed in triplicate 48 h after transfection using the dual-luciferase reporter assay kit (Promega, Madison, WI, USA). To control for general effects on transcription, Renilla luciferase was co-transfected in all reporter assays, and luciferase values represent the ratio luciferase/Renilla.

4.5. Co-Immunoprecipitation

Co-immunoprecipitation, employed to identify proteins directly or indirectly interacting with hnRNP K, was carried out using the Pierce Direct IP kit according to the manufacturer’s instructions (Pierce Biotechnology, Rockford, IL, USA). Briefly, total cellular extracts, obtained with IP Lysis/Wash buffer containing a protease inhibitor cocktail, were pre-cleared with agarose resin and passed through an AminoLink plus resin with covalently immobilized anti-hnRNP K antibodies (Santa Cruz Biotechnology, Santa Cruz, CA, USA). The eluted proteins were subjected to MS for identification of the hnRNP K interactome. Control resin and quenched antibody-coupling resin without anti-hnRNP K antibodies were utilized to identify non-specific interactions. Common proteins between these groups and the co-IPs were excluded from bioinformatic analyses. For MS analysis, three independent experiments with LNCaP, PDB, and MDB were performed.
The hnRNP K co-IPs were analyzed by WB for hnRNP K and AR expression evaluation in LNCaP, PDB, and MDB cell lines (Figure S3).

4.6. MS Analysis

The samples were processed by the FASP Protein Digestion Kit (Expedeon, Harston, UK). Briefly, the co-IP samples were mixed with 0.3 mL of 8 M urea in 0.1 M Tris/HCl pH 8.5 (UA solution), loaded into the filtration devices, and alkylated in 0.1 mL of 50 mM iodoacetamide in UA solution for 1 h in darkness at room temperature. The samples were digested using, sequentially, 1 µg of LysC and 3 µg of Trypsin in 50 mM NaHCO3 solution at 37 °C overnight. Peptides were collected by centrifugation of the filter units for 10 min, and the filter devices were rinsed with two 40 µL washes of 50 mM NaHCO3 and 50 µL 0.5 M NaCl to eliminate the hydrophobic interactions. Each sample digest was desalted on StageTips and analyzed by liquid chromatography–tandem mass spectrometry (LC–MS/MS).
All MS experiments were performed on a nanoscale high-performance liquid chromatography system connected to a hybrid linear trap quadrupole (LTQ) Orbitrap mass spectrometer. The MS instrument was operated in data-dependent mode to automatically switch between full-scan MS and MS/MS acquisition. Survey full-scan MS spectra were acquired in the Orbitrap analyzer with resolution R = 60,000. The 10 most intense peptide ions with charge states ≥2 were sequentially isolated and fragmented by collision-induced dissociation in the LTQ mass spectrometer. The raw mass spectrometric data were analyzed using the MaxQuant pipeline. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies, and proteome-wide protein quantification [71]. The statistical and the pathways analyses were done with the freely available Perseus [72] and Cytoscape software [73] respectively.

4.7. Network Analysis

Cytoscape (v 3.6.0) was used to visualize the protein–protein interaction network obtained from the above proteomic data (co-IP) and to highlight functional annotations [28]. Statistically significant proteins were used as the queries to interrogate Intact and Reactome database in order to obtain a protein–protein interaction network. The network was treated as undirected, and all duplicated edges and self-loop were removed. To highlight the most important proteins, a network analysis was performed. In particular, the degrees and closeness were assessed. The degree is the simplest topological index, corresponding to the number of nodes adjacent to a given node, where “adjacent” means directly connected. The degree allows an immediate evaluation of the regulatory relevance of the node. The closeness is a node centrality index. The closeness of a specific node is calculated by computing the shortest path between the node and all other nodes in the graph and then calculating the summa. From the biological point of view, the closeness can be interpreted as the “probability” of a protein to be functionally relevant for several other proteins. Moreover, the MaxLFQ expression data were loaded into the network as attribute node. Subsequently, a clusterization of significant proteins related to the expression data and to the topology of the network was performed as previously described, by clustermaker2 (Cytoscape App) [28] with the “autoSOME Clustering” algorithm (node attribute = expression; distance metric: Euclidean, ignore node with no data). Each obtained cluster was functionally tagged by an enrichment, performed by JEPETTO [29] (network: PSICQUIC, reference database: GO MF/prostate tissue) if statistically significant results were obtained.

4.8. Differential Expression and Enrichment Analyses

Differential expression analyses were carried out both on protein detected through MS spectroscopy for their binding to hnRNP K and on gene expression data extracted from four datasets published and available on the GEO portal (GSE32269, GSE101607, GSE29650, GSE66532).
In the first case, the log-transformed values of protein binding were imported and analyzed through linear model fitting and empirical Bayes approach provided by Limma R package [74]. Three different contrast matrices were generated: PDB versus LNCaP, MDB versus LNCaP, MDB versus PDB. Proteins showing a p-value ≤ 0.05 and a |log2FC| ≥ 1 were considered differentially interacting (DIhKPs).
Microarray data from GEO datasets were imported through GEOquery R package [75], “GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor”, and comparison analysis of genes was performed using Limma.
For each dataset, pairwise comparisons between sample groups, as indicated in Table S5, were performed. Genes having a |log2FC| ≥ 1 and the corresponding p-value ≤ 0.05 were considered statistically significant.
Enrichment analysis was performed on the 51 DIhKPs to identify the over-represented GO annotations regarding Biological Processes and Molecular Function, as well as the transcription factor binding sites through the category UCSC-TFBS provided by the functional annotation tool DAVID v 6.8 [76].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/7/1920/s1.

Author Contributions

Conceptualization, P.B.; Data curation, M.C. and I.G.; Formal analysis, I.G. and E.I.; Funding acquisition, M.R.G., A.P., N.F., F.B., and P.B.; Investigation, M.C., I.G., E.I., N.F., and P.B.; Project administration, P.B.; Resources, M.C. and N.F.; Supervision, M.R.G., A.P., and P.B.; Writing–original draft, P.B.; Writing–review & editing, C.C., N.F., and F.B.

Funding

This work was supported by Compagnia di San Paolo (2012.1588), IRCCS AOU San Martino-IST (Project 5x1000), MIUR PON02-00619, and National Research University Higher School of Economics RSF Grant 14-41-00039.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ADTAndrogen deprivation therapy
ARAndrogen receptor
ARGAndrogen-responsive gene
BICBicalutamide
Co-IPCo-immunoprecipitate
CRPCCastration-resistant prostate cancer
DAVIDDatabase for Annotation Visualization and Integrated Discovery
DEPDifferentially expressed protein
DIhKPDifferentially interacting hnRNP K protein
DHT5-α-dihydrotestosterone
GOGene ontology
hK-TFTranscription factor potentially regulating the expression of differentially interacting hnRNP K protein
LTQLinear trap quadrupole
MSMass spectrometry
PTMPost-translation modifications
TFTranscription factor
TFBSTranscription factor binding sites
PAGEPolyacrylamide gel electrophoresis
PCaProstate cancer
WBWestern blotting

References

  1. Katsogiannou, M.; Ziouziou, H.; Karaki, S.; Andrieu, C.; Henry de Villeneuve, M.; Rocchi, P. The hallmarks of castration-resistant prostate cancers. Cancer Treat. Rev. 2015, 41, 588–597. [Google Scholar] [CrossRef] [PubMed]
  2. Cattrini, C.; Zanardi, E.; Vallome, G.; Cavo, A.; Cerbone, L.; Di Meglio, A.; Fabbroni, C.; Latocca, M.M.; Rizzo, F.; Messina, C.; et al. Targeting androgen-independent pathways: New chances for patients with prostate cancer? Crit. Rev. Oncol. Hematol. 2017, 118, 42–53. [Google Scholar] [CrossRef] [PubMed]
  3. Zarif, J.C.; Miranti, C.K. The importance of non-nuclear AR signaling in prostate cancer progression and therapeutic resistance. Cell Signal. 2016, 28, 348–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Shah, R.B.; Mehra, R.; Chinnaiyan, A.M.; Shen, R.; Ghosh, D.; Zhou, M.; Macvicar, G.R.; Varambally, S.; Harwood, J.; Bismar, T.A.; et al. Androgen-independent prostate cancer is a heterogeneous group of diseases: Lessons from a rapid autopsy program. Cancer Res. 2004, 64, 9209–9216. [Google Scholar] [CrossRef] [PubMed]
  5. Beltran, H.; Yelensky, R.; Frampton, G.M.; Park, K.; Downing, S.R.; MacDonald, T.Y.; Jarosz, M.; Lipson, D.; Tagawa, S.T.; Nanus, D.M.; et al. Targeted next-generation sequencing of advanced prostate cancer identifies potential therapeutic targets and disease heterogeneity. Eur. Urol. 2012, 63, 920–926. [Google Scholar] [CrossRef] [PubMed]
  6. Ciccarese, C.; Massari, F.; Iacovelli, R.; Fiorentino, M.; Montironi, R.; Di Nunno, V.; Giunchi, F.; Brunelli, M.; Tortora, G. Prostate cancer heterogeneity: Discovering novel molecular targets for therapy. Cancer Treat. Rev. 2017, 54, 68–73. [Google Scholar] [CrossRef] [PubMed]
  7. Wyatt, A.W.; Gleave, M.E. Targeting the adaptive molecular landscape of castration-resistant prostate cancer. EMBO Mol. Med. 2015, 7, 878–894. [Google Scholar] [CrossRef] [PubMed]
  8. Emmons, M.F.; Faiao-Flores, F.; Smalley, K.S. The role of phenotypic plasticity in the escape of cancer cells from targeted therapy. Biochem. Pharmacol. 2016, 122, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Liu, S.; Kumari, S.; Hu, Q.; Senapati, D.; Venkadakrishnan, V.B.; Wang, D.; DePriest, A.D.; Schlanger, S.E.; Ben-Salem, S.; Valenzuela, M.M.; et al. A comprehensive analysis of coregulator recruitment, androgen receptor function and gene expression in prostate cancer. Elife 2017, 6, e28482. [Google Scholar] [CrossRef] [PubMed]
  10. Bomsztyk, K.; Denisenko, O.; Ostrowski, J. hnRNP K: One protein multiple processes. Bioessays 2004, 26, 629–638. [Google Scholar] [CrossRef] [PubMed]
  11. Mikula, M.; Ostrowski, J. Heterogeneous Nuclear Ribonucleoprotein K—Hnrnp K. The Transcription Factor Encyclopedia. Genome Biol. 2012, 13, R24. Available online: http://www.cisreg.ca/cgi-bin/tfe/articles.pl?tfid=551 (accessed on 27 November 2017).
  12. Barboro, P.; Ferrari, N.; Balbi, C. Emerging roles of heterogeneous nuclear ribonucleoprotein K (hnRNP K) in cancer progression. Cancer Lett. 2014, 352, 152–159. [Google Scholar] [CrossRef] [PubMed]
  13. Barboro, P.; Repaci, E.; Rubagotti, A.; Salvi, S.; Boccardo, S.; Spina, B.; Truini, M.; Introini, C.; Puppo, P.; Ferrari, N.; et al. Heterogeneous nuclear ribonucleoprotein K: Altered pattern of expression associated with diagnosis and prognosis of prostate cancer. Br. J. Cancer 2009, 100, 1608–1616. [Google Scholar] [CrossRef] [PubMed]
  14. Barboro, P.; Salvi, S.; Rubagotti, A.; Boccardo, S.; Spina, B.; Truini, M.; Carmignani, G.; Introini, C.; Ferrari, N.; Boccardo, F.; et al. Prostate cancer: Prognostic significance of the association of heterogeneous nuclear ribonucleoprotein K and androgen receptor expression. Int. J. Oncol. 2014, 44, 1589–1598. [Google Scholar] [CrossRef] [PubMed]
  15. Mukhopadhyay, N.K.; Kim, J.; Cinar, B.; Ramachandran, A.; Hager, M.H.; Di Vizio, D.; Adam, R.M.; Rubin, M.A.; Raychaudhuri, P.; De Benedetti, A.; et al. Heterogeneous nuclear ribonucleoprotein K is a novel regulator of androgen receptor translation. Cancer Res. 2009, 69, 2210–2218. [Google Scholar] [CrossRef] [PubMed]
  16. Barboro, P.; Repaci, E.; Ferrari, N.; Rubagotti, A.; Boccardo, F.; Balbi, C. Androgen receptor and heterogeneous nuclear ribonucleoprotein K colocalize in the nucleoplasm and are modulated by bicalutamide and 4-hydroxy-tamoxifen in prostatic cancer cell lines. Prostate 2011, 71, 1466–1479. [Google Scholar] [CrossRef] [PubMed]
  17. Ciarlo, M.; Benelli, R.; Barbieri, O.; Minghelli, S.; Barboro, P.; Balbi, C.; Ferrari, N. Regulation of neuroendocrine differentiation by AKT/hnRNPK/AR/beta-catenin signaling in prostate cancer cells. Int. J. Cancer 2012, 131, 582–590. [Google Scholar] [CrossRef] [PubMed]
  18. Kelly, K.; Balk, S.P. Reprogramming to resist. Science 2017, 355, 29–30. [Google Scholar] [CrossRef] [PubMed]
  19. Ferrari, N.; Granata, I.; Capaia, M.; Piccirillo, M.; Guarracino, M.R.; Vene, R.; Brizzolara, A.; Petretto, A.; Inglese, E.; Morini, M.; et al. Adaptive phenotype drives resistance to androgen deprivation therapy in prostate cancer. Cell Commun. Signal. 2017, 15, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Barboro, P.; Borzi, L.; Repaci, E.; Ferrari, N.; Balbi, C. Androgen receptor activity is affected by both nuclear matrix localization and the phosphorylation status of the heterogeneous nuclear ribonucleoprotein K in anti-androgen-treated LNCaP cells. PLoS ONE 2013, 8, e79212. [Google Scholar] [CrossRef] [PubMed]
  21. Moumen, A.; Magill, C.; Dry, K.L.; Jackson, S.P. ATM-dependent phosphorylation of heterogeneous nuclear ribonucleoprotein K promotes p53 transcriptional activation in response to DNA damage. Cell Cycle 2013, 12, 698–704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Habelhah, H.; Shah, K.; Huang, L.; Burlingame, A.L.; Shokat, K.M.; Ronai, Z. Identification of new JNK substrate using ATP pocket mutant JNK and a corresponding ATP analogue. J. Biol. Chem. 2001, 276, 18090–18095. [Google Scholar] [CrossRef] [PubMed]
  23. Laury-Kleintop, L.D.; Tresini, M. Hammond, O. Compartmentalization of hnRNP-K during cell cycle progression and its interaction with calponin in the cytoplasm. J. Cell. Biochem. 2005, 95, 1042–1056. [Google Scholar] [CrossRef] [PubMed]
  24. Habelhah, H.; Shah, K.; Huang, L.; Ostareck-Lederer, A.; Burlingame, A.L.; Shokat, K.M.; Hentze, M.W.; Ronai, Z. ERK phosphorylation drives cytoplasmic accumulation of hnRNP-K and inhibition of mRNA translation. Nat. Cell Biol. 2001, 3, 325–330. [Google Scholar] [CrossRef] [PubMed]
  25. Sampson, N.; Neuwirt, H.; Puhr, M.; Klocker, H.; Eder, I.E. In vitro model systems to study androgen receptor signaling in prostate cancer. Endocr. Relat. Cancer 2013, 20, R49–R64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kimura, Y.; Nagata, K.; Suzuki, N.; Yokoyama, R.; Yamanaka, Y.; Kitamura, H.; Hirano, H.; Ohara, O. Characterization of multiple alternative forms of heterogeneous nuclear ribonucleoprotein K by phosphate-affinity electrophoresis. Proteomics 2010, 10, 3884–3895. [Google Scholar] [CrossRef] [PubMed]
  27. Mikula, M.; Dzwonek, A.; Karczmarski, J.; Rubel, T.; Dadlez, M.; Wyrwicz, L.S.; Bomsztyk, K.; Ostrowski, J. Landscape of the hnRNP K protein-protein interactome. Proteomics 2006, 6, 2395–2406. [Google Scholar] [CrossRef] [PubMed]
  28. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  29. Winterhalter, C.; Widera, P.; Krasnogor, N. JEPETTO: A Cytoscape plugin for gene set enrichment and topological analysis based on interaction networks. Bioinformatics 2013, 30, 1029–1030. [Google Scholar] [CrossRef] [PubMed]
  30. Sharma, N.L.; Massie, C.E.; Ramos-Montoya, A.; Zecchini, V.; Scott, H.E.; Lamb, A.D.; MacArthur, S.; Stark, R.; Warren, A.Y.; Mills, I.G.; et al. The androgen receptor induces a distinct transcriptional program in castration-resistant prostate cancer in man. Cancer Cell 2013, 23, 35–47. [Google Scholar] [CrossRef] [PubMed]
  31. Kim, Y.R.; Oh, K.J.; Park, R.Y.; Xuan, N.T.; Kang, T.W.; Kwon, D.D.; Choi, C.; Kim, M.S.; Nam, K.I.; Ahn, K.Y.; et al. HOXB13 promotes androgen independent growth of LNCaP prostate cancer cells by the activation of E2F signaling. Mol. Cancer 2010, 9, 124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Dong, Y.; Cai, Y.; Liu, B.; Jiao, X.; Li, Z.T.; Guo, D.Y.; Li, X.W.; Wang, Y.J.; Yang, D.K. HOXA13 is associated with unfavorable survival and acts as a novel oncogene in prostate carcinoma. Future Oncol. 2017, 13, 1505–1516. [Google Scholar] [CrossRef] [PubMed]
  33. Li, L.; Chang, W.; Yang, G.; Ren, C.; Park, S.; Karantanos, T.; Karanika, S.; Wang, J.; Yin, J.; Shah, P.K.; et al. Targeting poly(ADP-ribose) polymerase and the c-Myb-regulated DNA damage response pathway in castration-resistant prostate cancer. Sci. Signal. 2015, 7, ra47. [Google Scholar] [CrossRef] [PubMed]
  34. Alghadi, A.Y.; Khalil, A.M.; Alazab, R.S.; Aldaoud, N.H.; Zyoud, A.M.-S. Unique Expression of the XBP1 Gene Correlates with with Human Prostate Cancer. J. Investig. Genom. 2017, 4, 00091. [Google Scholar] [CrossRef]
  35. Olokpa, E.; Moss, P.E.; Stewart, L.V. Crosstalk between the Androgen Receptor and PPAR Gamma Signaling Pathways in the Prostate. PPAR Res. 2017, 2017, 9456020. [Google Scholar] [CrossRef] [PubMed]
  36. Mendiratta, P.; Mostaghel, E.; Guinney, J.; Tewari, A.K.; Porrello, A.; Barry, W.T.; Nelson, P.S.; Febbo, P.G. Genomic strategy for targeting therapy in castration-resistant prostate cancer. J. Clin. Oncol. 2009, 27, 2022–2029. [Google Scholar] [CrossRef] [PubMed]
  37. Tien, A.H.; Saladar, M.D. Androgen-Responsive Gene Expression in Prostate Cancer Progression, 1st ed.; Springer: New York, NY, USA, 2013; pp. 135–153. ISBN 978-1-4614-6182-1. [Google Scholar]
  38. Beitel, L. Androgen Receptor Interacting Proteins and Coregulators Table. 2010. Available online: http://androgendb.mcgill.ca/ARinteract.pdf (accessed on 5 December 2017).
  39. Larkin, S.E.; Zeidan, B.; Taylor, M.G.; Bickers, B.; Al-Ruwaili, J.; Aukim-Hastie, C.; Townsend, P.A. Proteomics in prostate cancer biomarker discovery. Expert Rev. Proteom. 2018, 7, 93–102. [Google Scholar] [CrossRef] [PubMed]
  40. Prensner, J.R.; Rubin, M.A.; Wei, J.T.; Chinnaiyan, A.M. Beyond PSA: The next generation of prostate cancer biomarkers. Sci. Transl. Med. 2012, 4, 127rv3. [Google Scholar] [CrossRef] [PubMed]
  41. Hatakeyama, S.; Yoneyama, T.; Tobisawa, Y.; Ohyama, C. Recent progress and perspectives on prostate cancer biomarkers. Int. J. Clin. Oncol. 2017, 22, 214–221. [Google Scholar] [CrossRef] [PubMed]
  42. Intasqui, P.; Bertolla, R.P.; Sadi, M.V. Prostate cancer proteomics: Clinically useful protein biomarkers and future perspectives. Expert Rev. Proteom. 2017, 15, 65–79. [Google Scholar] [CrossRef] [PubMed]
  43. Szabo, Z.; Hamalainen, J.; Loikkanen, I.; Moilanen, A.M.; Hirvikoski, P.; Vaisanen, T.; Paavonen, T.K.; Vaarala, M.H. Sorbitol dehydrogenase expression is regulated by androgens in the human prostate. Oncol. Rep. 2010, 23, 1233–1239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Qian, X.; Li, C.; Pang, B.; Xue, M.; Wang, J.; Zhou, J. Spondin-2 (SPON2), a more prostate-cancer-specific diagnostic biomarker. PLoS ONE 2012, 7, e37225. [Google Scholar] [CrossRef] [PubMed]
  45. Zhu, B.P.; Guo, Z.Q.; Lin, L.; Liu, Q. Serum BSP, PSADT, and Spondin-2 levels in prostate cancer and the diagnostic significance of their ROC curves in bone metastasis. Eur. Rev. Med. Pharmacol. Sci. 2017, 21, 61–67. [Google Scholar] [PubMed]
  46. Vainio, P.; Wolf, M.; Edgren, H.; He, T.; Kohonen, P.; Mpindi, J.P.; Smit, F.; Verhaegh, G.; Schalken, J.; Perala, M.; et al. Integrative genomic, transcriptomic, and RNAi analysis indicates a potential oncogenic role for FAM110B in castration-resistant prostate cancer. Prostate 2011, 72, 789–802. [Google Scholar] [CrossRef] [PubMed]
  47. Honda, K. The biological role of actinin-4 (ACTN4) in malignant phenotypes of cancer. Cell Biosci. 2015, 5, 41. [Google Scholar] [CrossRef] [PubMed]
  48. Jasavala, R.; Martinez, H.; Thumar, J.; Andaya, A.; Gingras, A.C.; Eng, J.K.; Aebersold, R.; Han, D.K.; Wright, M.E. Identification of putative androgen receptor interaction protein modules: Cytoskeleton and endosomes modulate androgen receptor signaling in prostate cancer cells. Mol. Cell Proteom. 2007, 6, 252–271. [Google Scholar] [CrossRef] [PubMed]
  49. Cordonnier, T.; Bishop, J.L.; Shiota, M.; Nip, K.M.; Thaper, D.; Vahid, S.; Heroux, D.; Gleave, M.; Zoubeidi, A. Hsp27 regulates EGF/beta-catenin mediated epithelial to mesenchymal transition in prostate cancer. Int. J. Cancer 2015, 136, E496–E507. [Google Scholar] [CrossRef] [PubMed]
  50. Zoubeidi, A.; Zardan, A.; Beraldi, E.; Fazli, L.; Sowery, R.; Rennie, P.; Nelson, C.; Gleave, M. Cooperative interactions between androgen receptor (AR) and heat-shock protein 27 facilitate AR transcriptional activity. Cancer Res. 2007, 67, 10455–10465. [Google Scholar] [CrossRef] [PubMed]
  51. Rocchi, P.; Beraldi, E.; Ettinger, S.; Fazli, L.; Vessella, R.L.; Nelson, C.; Gleave, M. Increased Hsp27 after androgen ablation facilitates androgen-independent progression in prostate cancer via signal transducers and activators of transcription 3-mediated suppression of apoptosis. Cancer Res. 2005, 65, 11083–11093. [Google Scholar] [CrossRef] [PubMed]
  52. Stockley, J.; Markert, E.; Zhou, Y.; Robson, C.N.; Elliott, D.J.; Lindberg, J.; Leung, H.Y.; Rajan, P. The RNA-binding protein Sam68 regulates expression and transcription function of the androgen receptor splice variant AR-V7. Sci. Rep. 2015, 5, 13426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Frisone, P.; Pradella, D.; Di Matteo, A.; Belloni, E.; Ghigna, C.; Paronetto, M.P. SAM68: Signal Transduction and RNA Metabolism in Human Cancer. Biomed Res. Int. 2015, 2015, 528954. [Google Scholar] [CrossRef] [PubMed]
  54. Busa, R.; Paronetto, M.P.; Farini, D.; Pierantozzi, E.; Botti, F.; Angelini, D.F.; Attisani, F.; Vespasiani, G.; Sette, C. The RNA-binding protein Sam68 contributes to proliferation and survival of human prostate cancer cells. Oncogene 2007, 26, 4372–4382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zhou, L.; Xia, D.; Zhu, J.; Chen, Y.; Chen, G.; Mo, R.; Zeng, Y.; Dai, Q.; He, H.; Liang, Y.; et al. Enhanced expression of IMPDH2 promotes metastasis and advanced tumor progression in patients with prostate cancer. Clin. Transl. Oncol. 2014, 16, 906–913. [Google Scholar] [CrossRef] [PubMed]
  56. Albitar, M.; Ma, W.; Lund, L.; Albitar, F.; Diep, K.; Fritsche, H.A.; Shore, N. Predicting Prostate Biopsy Results Using a Panel of Plasma and Urine Biomarkers Combined in a Scoring System. J. Cancer 2015, 7, 297–303. [Google Scholar] [CrossRef] [PubMed]
  57. Brizzolara, A.; Benelli, R.; Vene, R.; Barboro, P.; Poggi, A.; Tosetti, F.; Ferrari, N. The ErbB family and androgen receptor signaling are targets of Celecoxib in prostate cancer. Cancer Lett. 2017, 400, 9–17. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, L.G.; Johnson, E.M.; Kinoshita, Y.; Babb, J.S.; Buckley, M.T.; Liebes, L.F.; Melamed, J.; Liu, X.M.; Kurek, R.; Ossowski, L.; et al. Androgen receptor overexpression in prostate cancer linked to Pur alpha loss from a novel repressor complex. Cancer Res. 2008, 68, 2678–2688. [Google Scholar] [CrossRef] [PubMed]
  59. Hang, S.P.; Tang, Y.H.; Smith, R. Functional diversity of hnRNPs: Past, present and perspectives. Biochem. J. 2010, 430, 379–392. [Google Scholar] [CrossRef]
  60. Rahman-Roblick, R.; Roblick, U.J.; Hellman, U.; Conrotto, P.; Liu, T.; Becker, S.; Hirschberg, D.; Jörnvall, H.; Auer, G.; Wiman, K.G. p53 targets identified by protein expression profiling. Proc. Natl. Acad. Sci. USA 2007, 104, 5401–5406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Chan, J.Y.H.; Huang, S.M.; Liu, S.T.; Huang, C.H. The transactivation domain of heterogeneous nuclear ribonucleoprotein K overlaps its nuclear shuttling domain. Int. J. Biochem. Cell Biol. 2008, 40, 2078–2089. [Google Scholar] [CrossRef] [PubMed]
  62. Vlachostergios, P.J.; Puca, L.; Beltran, H. Emerging Variants of Castration-Resistant Prostate Cancer. Curr. Oncol. Rep. 2017, 19, 32. [Google Scholar] [CrossRef] [PubMed]
  63. Lamb, A.D.; Massie, C.E.; Neal, D.E. The transcriptional programme of the androgen receptor (AR) in prostate cancer. BJU Int. 2014, 113, 358–366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Holzbeierlein, J.; Lal, P.; LaTulippe, E.; Smith, A.; Satagopan, J.; Zhang, L.; Ryan, C.; Smith, S.; Scher, H.; Scardino, P.; et al. Gene expression analysis of human prostate carcinoma during hormonal therapy identifies androgen-responsive genes and mechanisms of therapy resistance. Am. J. Pathol. 2004, 164, 217–227. [Google Scholar] [CrossRef]
  65. Coutinho, I.; Day, T.K.; Tilley, W.D.; Selth, L.A. Androgen receptor signaling in castration-resistant prostate cancer: A lesson in persistence. Endocr. Relat. Cancer 2016, 23, T179–T197. [Google Scholar] [CrossRef] [PubMed]
  66. Shiota, M.; Yokomizo, A.; Fujimoto, N.; Naito, S. Androgen receptor cofactors in prostate cancer: Potential therapeutic targets of castration-resistant prostate cancer. Curr. Cancer Drug Targets 2011, 11, 870–881. [Google Scholar] [CrossRef] [PubMed]
  67. Culig, Z. Androgen Receptor Coactivators in Regulation of Growth and Differentiation in Prostate Cancer. J. Cell. Physiol. 2015, 231, 270–274. [Google Scholar] [CrossRef] [PubMed]
  68. Obinata, D.; Takayama, K.; Takahashi, S.; Inoue, S. Crosstalk of the Androgen Receptor with Transcriptional Collaborators: Potential Therapeutic Targets for Castration-Resistant Prostate Cancer. Cancers 2017, 9, 22. [Google Scholar] [CrossRef] [PubMed]
  69. Moss, P.E.; Lyles, B.E.; Stewart, L.V. The PPARgamma ligand ciglitazone regulates androgen receptor activation differently in androgen-dependent versus androgen-independent human prostate cancer cells. Exp. Cell Res. 2010, 316, 3478–3488. [Google Scholar] [CrossRef] [PubMed]
  70. Yang, J.P.; Reddy, T.R.; Truong, K.T.; Suhasini, M.; Wong-Staal, F. Functional interaction of Sam68 and heterogeneous nuclear ribonucleoprotein K. Oncogene 2002, 21, 7187–7194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367–1372. [Google Scholar] [CrossRef] [PubMed]
  72. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar] [CrossRef] [PubMed]
  73. Cline, M.S.; Smoot, M.; Cerami, E.; Kuchinsky, A.; Landys, N.; Workman, C.; Christmas, R.; Avila-Campilo, I.; Creech, M.; Gross, B.; et al. Integration of biological networks and gene expression data using Cytoscape. Nat. Protoc. 2007, 2, 2366–2382. [Google Scholar] [CrossRef] [PubMed]
  74. Lopez-Romero, P. Pre-processing and differential expression analysis of Agilent microRNA arrays using the AgiMicroRna Bioconductor library. BMC Genom. 2011, 12, 64. [Google Scholar] [CrossRef] [PubMed]
  75. Davis, S.; Meltzer, P.S. GEOquery: A bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics 2007, 23, 1846–1847. [Google Scholar] [CrossRef] [PubMed]
  76. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 01920 g001 550
Figure 1. Effects of heterogeneous nuclear ribonucleoprotein K (hnRNP K) on the regulation of androgen receptor (AR) expression and AR transcriptional activity in LNCaP, PDB, and MDB cells. (a,d) Representative western blotting (WB) analysis carried out using antibodies against hnRNP K, AR, and prostate-specific antigen (PSA) and (b,e) quantitative analysis of total proteins extract from (a) hnRNP K-silenced (siHK) and -non-silenced (siNT) LNCaP, PDB, and MDB cell lines or (d) from control (Ctrl) PDB and MDB cell lines grown in the appropriate medium (see Material and Methods) or in restored LNCaP growth medium without bicalutamide (-BIC) and with 5-α-dihydrotestosterone(+DHT) for 2, 3 weeks (wk) or 1, 2, 3 months (m). The ordinate represents the mean ± SE (b) of the percentage of selected protein expression in hnRNP K-silenced cell lines or (e) the relative amounts of proteins determined by quantitative analysis. (c) AR transcriptional activity determined in hnRNP K-silenced and non-silenced LNCaP by the luciferase activity assay; * p < 0.05, ** p < 0.01, and *** p < 0.001 (Student’s t-test).
Figure 1. Effects of heterogeneous nuclear ribonucleoprotein K (hnRNP K) on the regulation of androgen receptor (AR) expression and AR transcriptional activity in LNCaP, PDB, and MDB cells. (a,d) Representative western blotting (WB) analysis carried out using antibodies against hnRNP K, AR, and prostate-specific antigen (PSA) and (b,e) quantitative analysis of total proteins extract from (a) hnRNP K-silenced (siHK) and -non-silenced (siNT) LNCaP, PDB, and MDB cell lines or (d) from control (Ctrl) PDB and MDB cell lines grown in the appropriate medium (see Material and Methods) or in restored LNCaP growth medium without bicalutamide (-BIC) and with 5-α-dihydrotestosterone(+DHT) for 2, 3 weeks (wk) or 1, 2, 3 months (m). The ordinate represents the mean ± SE (b) of the percentage of selected protein expression in hnRNP K-silenced cell lines or (e) the relative amounts of proteins determined by quantitative analysis. (c) AR transcriptional activity determined in hnRNP K-silenced and non-silenced LNCaP by the luciferase activity assay; * p < 0.05, ** p < 0.01, and *** p < 0.001 (Student’s t-test).
Ijms 19 01920 g001
Ijms 19 01920 g002 550
Figure 2. HnRNP K decreased expression and altered phosphorylation in resistant cell lines PDB and MDB compared with LNCaP. (a) WB analysis of hnRNP K expression in total extracts from LNCaP, PDB, and MDB cell lines. The histogram on the right side represents the mean ± SE of relative amount of proteins determined in six WBs. (b) The quantitative analysis of hnRNP K phosphorylation was carried out using 1D Phos-tag and WB analysis of LNCaP, PDB, and MDB total extracts. The histograms on the right side represent the mean percentage ± SE of phosphorylated hnRNP K (p-hK) evaluated in four experiments. The hnRNP K isoforms with minimal (phK0), intermediate (phK1), or maximal (phK2) phosphorylation are indicated; * p < 0.05 and ** p < 0.01 (Student’s t-test).
Figure 2. HnRNP K decreased expression and altered phosphorylation in resistant cell lines PDB and MDB compared with LNCaP. (a) WB analysis of hnRNP K expression in total extracts from LNCaP, PDB, and MDB cell lines. The histogram on the right side represents the mean ± SE of relative amount of proteins determined in six WBs. (b) The quantitative analysis of hnRNP K phosphorylation was carried out using 1D Phos-tag and WB analysis of LNCaP, PDB, and MDB total extracts. The histograms on the right side represent the mean percentage ± SE of phosphorylated hnRNP K (p-hK) evaluated in four experiments. The hnRNP K isoforms with minimal (phK0), intermediate (phK1), or maximal (phK2) phosphorylation are indicated; * p < 0.05 and ** p < 0.01 (Student’s t-test).
Ijms 19 01920 g002
Ijms 19 01920 g003 550
Figure 3. The hnRNP K phosphorylation status correlates with AR activity in prostate cancer cell lines. Nuclear (a) and cytoplasmic (b) profiles of hnRNP K phosphorylated isoforms in androgen-responsive LNCaP and resistant cell lines PDB, MDB, 22Rv1, and PC3, evaluated using 2D Phos-tag and WB analysis. The membranes were probed with an anti-hnRNP K antibody. The phK0, phK1, and phK2 identified in Phos-tag 1D (Figure 2b) are indicated.
Figure 3. The hnRNP K phosphorylation status correlates with AR activity in prostate cancer cell lines. Nuclear (a) and cytoplasmic (b) profiles of hnRNP K phosphorylated isoforms in androgen-responsive LNCaP and resistant cell lines PDB, MDB, 22Rv1, and PC3, evaluated using 2D Phos-tag and WB analysis. The membranes were probed with an anti-hnRNP K antibody. The phK0, phK1, and phK2 identified in Phos-tag 1D (Figure 2b) are indicated.
Ijms 19 01920 g003
Ijms 19 01920 g004 550
Figure 4. The hnRNP K interactome depends on the cellular phenotype. Protein–protein interaction networks of co-immunoprecipitates (co-IPs) obtained from LNCaP, PDB, and MDB, highlighting GO molecular function annotation restricted to prostate tissue grouped in three clusters, were realized in Cytoscape environment. The dot color refers to the Z score of protein intensity calculated on the basis of all three experiments: upregulated (red), downregulated (blue), and medium (black) value.
Figure 4. The hnRNP K interactome depends on the cellular phenotype. Protein–protein interaction networks of co-immunoprecipitates (co-IPs) obtained from LNCaP, PDB, and MDB, highlighting GO molecular function annotation restricted to prostate tissue grouped in three clusters, were realized in Cytoscape environment. The dot color refers to the Z score of protein intensity calculated on the basis of all three experiments: upregulated (red), downregulated (blue), and medium (black) value.
Ijms 19 01920 g004
Ijms 19 01920 g005 550
Figure 5. Transcriptional regulatory subnetworks of transcription factors (hK-TFs)potentially regulating the expression of differentially interacting hnRNP K proteins (DIhKPs) and selected DIhKPs, predicted and enriched by DAVID UCSC_TFBS option. The comparison of the subnetworks reveals dynamic remodeling of the hnRNP K interactome in LNCaP, PDB, and MDB. The six enriched hK-TFs are represented in colored circles, and DIhKPs in white circles, the edges are colored according to the hK-TFs. The subnetworks were grouped by the specific comparison of proteins abundance: PDBvsLNCaP (a), MDBvsLNCaP (b), MDB and PDBvsLNCaP (c), and MDBvsPDB (d). The red arrows indicate upregulated DIhKP expression; in all other cases, DIhKPs expression is downregulated.
Figure 5. Transcriptional regulatory subnetworks of transcription factors (hK-TFs)potentially regulating the expression of differentially interacting hnRNP K proteins (DIhKPs) and selected DIhKPs, predicted and enriched by DAVID UCSC_TFBS option. The comparison of the subnetworks reveals dynamic remodeling of the hnRNP K interactome in LNCaP, PDB, and MDB. The six enriched hK-TFs are represented in colored circles, and DIhKPs in white circles, the edges are colored according to the hK-TFs. The subnetworks were grouped by the specific comparison of proteins abundance: PDBvsLNCaP (a), MDBvsLNCaP (b), MDB and PDBvsLNCaP (c), and MDBvsPDB (d). The red arrows indicate upregulated DIhKP expression; in all other cases, DIhKPs expression is downregulated.
Ijms 19 01920 g005
Ijms 19 01920 g006 550
Figure 6. Schematic representation of the cross-talk network between hnRNP K, AR, the three AR-regulated hK-TFs, and the eight proteins selected for the hnRNP K–AR-related signature. The eight DIhKPs (heavy and light gray circles) directly or indirectly interact with hnRNP K (blue dotted lines). The hK-TFs (black circles) potentially control hnRNP K, and the eight DIhKPs (black dotted lines). AR coregulators (heavy gray circles) and androgen-responsive gene (ARGs) (highlighted with red *) are associated with AR (red solid lines).
Figure 6. Schematic representation of the cross-talk network between hnRNP K, AR, the three AR-regulated hK-TFs, and the eight proteins selected for the hnRNP K–AR-related signature. The eight DIhKPs (heavy and light gray circles) directly or indirectly interact with hnRNP K (blue dotted lines). The hK-TFs (black circles) potentially control hnRNP K, and the eight DIhKPs (black dotted lines). AR coregulators (heavy gray circles) and androgen-responsive gene (ARGs) (highlighted with red *) are associated with AR (red solid lines).
Ijms 19 01920 g006
Table
Table 1. Quantitative analysis of the spots reported in Figure 3. The assignment of the alternatively spliced isoforms 1 and 2, non-phosphorylated, and S116-, S284-, S353-phosphorylated forms was carried out according to the schematic representation showed in Figure S1.
Table 1. Quantitative analysis of the spots reported in Figure 3. The assignment of the alternatively spliced isoforms 1 and 2, non-phosphorylated, and S116-, S284-, S353-phosphorylated forms was carried out according to the schematic representation showed in Figure S1.
Cell LineLNCaP(AR-FL 1)PDB(AR-FL)MDB(AR-FL)22Rv1 (ARv7 2)PC3(AR-Null 3)
AR StatusActiveHypersensitiveInactiveAndrogen IndependentNot Detected
Nucleus
Alternatively spliced
isoform 1 (%)50.166.762.576.666.8
isoform 2 (%)49.933.337.523.433.2
Phosphorylated forms
pS116 (%)17.119.515.01.70.7
pS284 (%)37.645.823.236.122.8
pS353 (%)21.68.416.21.810.6
Non-phosphorylated forms (%)23.626.245.560.365.8
Cytoplasm
Alternatively spliced
isoform 1 (%)38.745.166.163.776.2
isoform 2 (%)61.354.933.936.323.8
Phosphorylated forms
pS116 (%)0.08.29.913.36.8
pS284 (%)0.00.042.041.30.0
pS353 (%)24.935.721.517.823.8
Non-phosphorylated forms (%)75.156.126.627.769.4
1 Cell line expressing AR full length; 2 Cell line expressing AR splicing isoform ARv7; 3 Cell line not expressing AR.
Table
Table 2. Selected DIhKPs associated with hK-TF, AR, PCa and differential expression in human prostate tissues, delineating a hnRNP K–AR-related signature.
Table 2. Selected DIhKPs associated with hK-TF, AR, PCa and differential expression in human prostate tissues, delineating a hnRNP K–AR-related signature.
Gene symbol/Protein NameDifferential hnRNP K InteractionAssociated hK-TFAssociation with ARHuman Prostate Tissues Differential ExpressionAssociation with PCa (Literature Evidences)
Proteomic StudiesTranscriptomic Studies
KLK3/PSAMDB&PDB vs. LNCaP ↓ *PAX5AR controlledPCa vs. BPH ↓PCa & mCRPC vs. NT ↑mCRPC vs. PCa ↓mCRPC non-AR driven vs. AR driven ↓Biomarker [40,41,42]
SORD/Sorbitol dehydrogenaseMDB&PDB vs. LNCaP ↓PPARGAR controlled PCa vs. NT ↑mCRPC vs. PCa ↓mCRPC non-AR driven vs. AR driven ↓Increased expression in high Gleason PCa; reduced expression after castration [43]
SPON2/Spondin2MDB&PDB vs. LNCaP ↓HOX13AR controlled mCRPC non-AR driven vs. AR driven ↓Biomarker [44,45]
ATP1B1/Sodium-potassium-transporting ATPase subunit beta1PDB vs. LNCaP ↑ **PPARG, XBP1, PAX5, HOX13, CMYBAR controlled Frequently overexpressed and amplified in CRPC [46]
ACTN4/Actinin4PDB vs. LNCaP ↓PPARG, XBP1, PAX5, HOX13, ZIC2, CMYBAR coactivator, direct interaction PCa vs. NT ↓ Reduced expression in high-grade PCa nuclei [42,47,48]
HSPB1/HSP27PDB vs. LNCaP ↓XBP1, PAX5AR coactivator, direct interactionPCa vs. NT & BPH ↑PCa vs. NT ↑­ Role in PCa epithelial to mesenchymal transition and resistance [49,50,51]
KHDRBS1/Sam68MDB vs. LNCaP ↓PPARG, PAX5, ZIC2, CMYBAR coactivator/corepressor Role in regulating AR splice variant, cell proliferation, and survival to chemotherapeutic agents [52,53,54]
IMPDH2/Inosine-5-monophosphate dehydrogenase 2MDB vs. LNCaP ↓PPARG, XBP1, PAX5, HOX13, CMYB PCa vs. NT ↑PCa & mCRPC vs. NT ↑ PCa vs. mCRPC non-AR driven ↑Role in PCa metastasis and progression [55,56]
↓* downregulated; ↑ ** upregulated.

Share and Cite

MDPI and ACS Style

Capaia, M.; Granata, I.; Guarracino, M.; Petretto, A.; Inglese, E.; Cattrini, C.; Ferrari, N.; Boccardo, F.; Barboro, P. A hnRNP K–AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer. Int. J. Mol. Sci. 2018, 19, 1920. https://doi.org/10.3390/ijms19071920

AMA Style

Capaia M, Granata I, Guarracino M, Petretto A, Inglese E, Cattrini C, Ferrari N, Boccardo F, Barboro P. A hnRNP K–AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer. International Journal of Molecular Sciences. 2018; 19(7):1920. https://doi.org/10.3390/ijms19071920

Chicago/Turabian Style

Capaia, Matteo, Ilaria Granata, Mario Guarracino, Andrea Petretto, Elvira Inglese, Carlo Cattrini, Nicoletta Ferrari, Francesco Boccardo, and Paola Barboro. 2018. "A hnRNP K–AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer" International Journal of Molecular Sciences 19, no. 7: 1920. https://doi.org/10.3390/ijms19071920

APA Style

Capaia, M., Granata, I., Guarracino, M., Petretto, A., Inglese, E., Cattrini, C., Ferrari, N., Boccardo, F., & Barboro, P. (2018). A hnRNP K–AR-Related Signature Reflects Progression toward Castration-Resistant Prostate Cancer. International Journal of Molecular Sciences, 19(7), 1920. https://doi.org/10.3390/ijms19071920

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop