- 1Center for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
- 2Basic Sciences Department, Faculty of Medicine and Health Sciences, Universitat Internacional de Catalunya, Barcelona, Spain
- 3Universitat Pompeu Fabra (UPF), Barcelona, Spain
Background: Breast cancer cells enter into the cell cycle following progestin exposure by the activation of signalling cascades involving a plethora of enzymes, transcription factors and co-factors that transmit the external signal from the cell membrane to chromatin, ultimately leading to a change of the gene expression program. Although many of the events within the signalling network have been described in isolation, how they globally team up to generate the final cell response is unclear.
Methods: In this study we used antibody microarrays and phosphoproteomics to reveal a dynamic global signalling map that reveals new key regulated proteins and phosphor-sites and links between previously known and novel pathways. T47D breast cancer cells were used, and phospho-sites and pathways highlighted were validated using specific antibodies and phenotypic assays. Bioinformatic analysis revealed an enrichment in novel signalling pathways, a coordinated response between cellular compartments and protein complexes.
Results: Detailed analysis of the data revealed intriguing changes in protein complexes involved in nuclear structure, epithelial to mesenchyme transition (EMT), cell adhesion, as well as transcription factors previously not associated with breast cancer cell proliferation. Pathway analysis confirmed the key role of the MAPK signalling cascade following progesterone and additional hormone regulated phospho-sites were identified. Full network analysis shows the activation of new signalling pathways previously not associated with progesterone signalling in T47D breast cancer cells such as ERBB and TRK. As different post-translational modifications can mediate complex crosstalk mechanisms and massive PARylation is also rapidly induced by progestins, we provide details of important chromatin regulatory complexes containing both phosphorylated and PARylated proteins.
Conclusions: This study contributes an important resource for the scientific community, as it identifies novel players and connections meaningful for breast cancer cell biology and potentially relevant for cancer management.
Background
Female steroid hormones, oestrogen and progesterone play a key role not only in the normal development of target tissues during puberty, pregnancy and menopause but also in breast and endometrium cancer cell proliferation. Breast cancer cells respond to progestin exposure with two intermingled pathways that culminate in extensive gene expression changes and entry in the cell cycle. The classical view is that the hormone diffuses through the cell membrane and binds to intracellular progesterone receptors (PR), which are maintained in an inactive state by a chaperone complex, including Heat Shock Proteins 70 and 90 (HSP70/90). Upon hormone binding, PR weakens its interaction with the chaperones, dimerizes and moves to chromatin where it eventually binds to palindromic DNA sequences called progesterone responsive elements (PREs) (1). Once bound to chromatin, PR recruits various co-regulators and chromatin remodellers that modulate access for the transcription machinery including RNA polymerase II (2).
This simplified model was further complicated later by the finding that a tiny fraction of PR (3-5%) is attached to the cell membrane via palmitoylation at C820 (3, 4), forming a complex with estrogen receptor alpha (ERa) (5). Upon binding progestins, the membrane anchored PR activates SRC, either directly (6) or via ERa, initiating a kinase signalling pathway that ends in activation of the extracellular signal-regulated kinase (ERK) (5). ERK1 phosphorylates intracellular PR at S294 favoring its dissociation from the chaperone complex. In the cell nucleus ERK1 activates MSK1 (Mitogen-and stress activated protein kinase 1), resulting in the formation of a ternary complex of PR-ERK1-MSK1, which is the active form of PR able to regulate chromatin structure and gene expression. ERK also activates cyclin dependent kinase 2 (CDK2) that in turns activates ARTD1 (ADP-ribose transferase 1) by phosphorylating two serines in the NAD+ binding pocket (7). Phosphorylation contributes to dissociation of histone H1 and H2A/H2B dimers (7, 8) and to local chromatin opening by further recruitment of transcription factors, co-regulators, histone modifiers (PCA, P3000) and ATP-dependent chromatin remodellers (NURF and BAF), ultimately leading to the activation of gene expression changes (9–11). Moreover, there is also evidence for the activation by progesterone of other signaling pathways induced by progesterone, such as AKT (12), cAMP (13, 14), GSK3 (15) and STAT (16). However, many of these studies used different types of cells derived from endometrial or ovarian tissues (17, 18).
In addition to the key role of the kinase cascades, we have identified a pivotal role for another post-translational modification in progestin induced gene regulation; namely Poly-ADP-ribosylation (PARylation). As discussed above the PAR polymerase PARP1, also known as ADP-ribosyltransferase 1 (ARTD1), is activated within the initial minutes following hormone exposure via phosphorylation by CDK2 (7), giving rise to a large increase in PARylation within the cell nucleus (7). PARylation of ARTD1 itself and of chromatin proteins is essential for the initial dissociation of histone H1 (19, 20). We also found that degradation of PAR to ADP-Ribose by PAR glycohydrolase (PARG) is required for complete chromatin remodelling and activation of the gene expression network (21). Mass spec analysis of the proteins interacting with PAR in T47D cells exposed to progestins revealed structural proteins, DNA damage response proteins and chromatin modifying enzymes (21). One key enzyme identified in this study was NUDT5 or NUDIX5 (Nudix hydrolase 5), which hydrolyses ADPR to AMP and ribose-5-phosphate. Subsequently, we found that upon dephosphorylation at T45, NUDT5 can use ADPR and diphosphate for the synthesis of ATP (21). In this way, part of the ATP consumed during the synthesis of NAD+ and stored in PAR is recovered and used for chromatin remodeling and changes in gene expression. The synthesis of nuclear ATP is transient, peaking at 40 minutes after hormone exposure and returning to basal levels after ~60 minutes. However, although we know that nuclear ATP synthesis is essential for the initial chromatin remodelling, the role of the nuclear ATP at later time points is unclear. We can envision several hypotheses; as a direct, local source of ATP for the massive amount of ATP-dependent chromatin remodelling and 3D conformational changes induced by hormone (9, 22) or to facilitate phase separation of chromatin fiber (23). In any case, we know that nuclear ATP synthesis by NUDT5 is essential for the generation and maintenance of the cancer stem cell population (24).
Over the past years, there has been a large number of studies investigating the role and mechanism of action of one or more of the pathway components in response to progesterone exposure, revealing a dynamic crosstalk between canonical pathways (25–28). However, these studies focus on one or very few components at a time, and do not explain how these various pathways interact and coordinate the cell response to hormones. The work described here aims to provide a more comprehensive map of progesterone signalling in breast cancer using luminal A T47D breast cancer cells as a model, combining antibody arrays technology, shotgun proteomics, and previously published PARylation datasets to develop for the first-time a global map of the dynamic signalling events induced by progestins in T47D breast cancer cells.
Methods
Cell culture
The hormone receptor positive breast cancer cell line T47DM (CLS Cat# 300353/p525_T-47D, RRID : CVCL_0553) was used in all experiments unless otherwise stated. T47DM cells were routinely grown in RPMI (Supplemented with 10% foetal bovine serum (FBS), penicillin/streptomycin (pen/strep), L-glutamine (L-glut) as previously described (21). For hormone induction experiments, cells were seeded at a concentration of 5 x 106 per 150mm cell culture dish in RPM1 white (15% charcoal stripped FBS, Pen/strep, L-glut) for 48 hours. 16 hours prior to hormone induction for time indicated (10nM R5020), the medium was replaced with RPM1 white (0% FBS, Pen/strep, L-glut). Samples were harvested at the time points indicated.
BCA assay
The total protein content of the samples was calculated prior to antibody array, mass spec or western blotting analysis using BCA assay (Thermo Fisher, catalogue number 23227) according to manufacturer’s instructions.
Protein visualization
Changes in phosphorylation sites within individual proteins identified was confirmed by western blotting as previously described (19) using specific antibodies; Progesterone receptor (PGR) phospho-S162 (Abcam Cat# ab58564, RRID : AB_883089), and as a loading control, total PGR (Santa Cruz Biotechnology Cat# sc-7208, RRID : AB_2164331), total CDK2 (Santa Cruz Biotechnology Cat# sc-6248, RRID : AB_627238) or CDK2 phospho-T160 (Abcam Cat# ab47330, RRID : AB_869087).
Antibody microarray
Phosphorylation antibody array analysis was carried out by Kinexus™ using Kinexus ™ Antibody Microarray (KAM) technology. For each time point 3 biological replicates were prepared independently. For each replicate, 50ug of protein lysate was prepared and samples prepared by Kinexus™ in house (Kinexus Bioinformatics Corporation, RRID : SCR_012553). Signal quantification was performed using ImaGene 8.0 (ImaGene, RRID : SCR_002178) from BioDiscovery (BioDiscovery, RRID : SCR_004557). Background corrected raw intensity data was logarithmically transformed with base 2 and Z scores calculated (29). Any poor-quality spots based on morphology and/or background, were flagged as unreliable and removed from any subsequent analysis.
Mass spec sample preparation
Samples (triplicates from independent experiments, per time point) were precipitated with 6 volumes of cold acetone and the pellet was dissolved in 6M Urea/200mM ammonium bicarbonate at a concentration of 1 µg/µL. 260 µg of each sample was reduced with dithiothreitol (DTT, 10mM, 37°C, 60 min) and alkylated with iodoacetamide (IAM, 20mM, 25°, 30 min). The resulting protein extract was first diluted to 2M urea with 200 mM ammonium bicarbonate for digestion with endoproteinase LysC (1:10 w:w, 37°C, o/n, Wako, cat # 129-02541), and then diluted 2-fold with 200 mM ammonium bicarbonate for trypsin digestion (1:10 w:w, 37°C, 8h, Promega cat # V5113). Digested peptides were subjected to phosphopeptide enrichment using the “TiO2 Phosphopeptide Enrichment and Clean-up Kit (Pierce)”
45% of each enriched sample was analyzed using an Orbitrap Fusion Lumos with an EASY-Spray nanosource coupled to a nano-UPLC system (EASY-nanoLC 1000 liquid chromatograph) equipped with a 50-cm C18 column (EASY-Spray; 75 µm id, PepMap RSLC C18, 2 µm particles, 45 °C). Chromatographic gradients started at 5% buffer B with a flow rate of 300 nL/min and gradually increased to 22% buffer B in 105 min and to 32% in 10 minutes. After each analysis, the column was washed for 10 min with 95% buffer B (Buffer A: 0.1% formic acid in water. Buffer B: 0.1% formic acid in acetonitrile).
The mass spectrometer was operated in data-dependent acquisition mode, with full MS scans over a mass range of m/z 350–1500 with detection in the Orbitrap (120K resolution) and with auto gain control (AGC) set to 100,000. In each cycle of data-dependent acquisition analysis, following each survey scan, the most intense ions above a threshold ion count of 10,000 were selected for fragmentation with HCD at normalized collision energy of 28%. The number of selected precursor ions for fragmentation was determined by the “Top Speed” acquisition algorithm (maximum cycle time of 3 seconds), and a dynamic exclusion of 60 s was set. Fragment ion spectra were acquired in the ion trap with an AGC of 4000 and a maximum injection time of 300 ms.
Raw files, quantified peptides, proteins and phosphosites (5% FDR) are deposited and publicly available in the UCSD Data depository, MassIVE (Project ID: MSV000089820, ftp://MSV000089820@massive.ucsd.edu).
Phosphopeptide analysis
Acquired data were analyzed using the Proteome Discoverer software suite (v2.0, Thermo Fisher Scientific), and the Mascot search engine (v2.5.1, Matrix Science) was used for peptide identification. Data were searched against a Human protein database derived from SwissProt plus the most common contaminants. A precursor ion mass tolerance of 7 ppm at the MS1 level was used, and up to three missed cleavages for trypsin were allowed. The fragment ion mass tolerance was set to 0.5Da. Oxidation of Methionine, N-terminal protein acetylation and phosphorylation in Serine, Threonine and Tyrosine were defined as variable modification and carbamidomethylation of Cysteines was set as fixed modification. The identified peptides were filtered by 5%FDR. Peptide areas were obtained with the area under from extracted ion chromatogram using the “Precursor Ions Area Detector” node from Proteome Discoverer.
Bioinformatic procedures
Gene Ontology (GO) GO-Biological process (GO-BP), GO-Molecular Process (GO-MF), KEGG (KEGG, RRID : SCR_012773) and Biocarta (BioCarta Pathways, RRID : SCR_006917) pathway analysis of networks was carried out using GeneMania application (GeneMANIA, RRID : SCR_005709) (30) within Cytoscape (Cytoscape, RRID : SCR_003032). Clustering analysis, similarity analysis was carried out using GeneE. Network analysis was performed using Cytoscape v3.5 (31). The initial prior knowledge network (PKN) was generated based on known protein-protein interactions only validated experimentally. Significantly enriched pathways were analyzed within the network using CytoKEGG application within Cytoscape. The parent network and each of the individual pathway networks are available for visualisation and further analysis (cytoscape session link Supplementary Materials). Comprehensive resource of mammalian protein complexes (Corum analysis) was carried out using online tool http://mips.helmholtz-muenchen.de/corum/CORUM, RRID : SCR_002254 (Giurgiu et al., 2018). Functional classification GO biological process (BP), molecular function (MF) and cellular component (CC) were carried out using molecular signatures database (MSigD) within Gene Set Enrichment (GSEA) tool (Gene Set Enrichment Analysis, RRID : SCR_003199) and terms with a p value of less than 0.001 were considered significantly enriched (32, 33).
Kaplan Meyer and protein expression in clinical samples
Analysis of the overall survival of breast cancer patients using a Kaplan-Meier plot were carried out using KMPlotter (34), https://kmplot.com/analysis/to assess the correlation between the expression of the gene (mRNA) and survival in 4929 breast cancer patient samples. Sources for the databases include GEO, EGA, and TCGA. All patients’ samples were included in the analysis shown, i.e ER, PR status, subtype, lymph node status and grade. Analysis of protein expression levels in breast tumor versus normal samples were representative of those within the Human Protein Atlas database (35) http://www.proteinatlas.org.
Results
Prior knowledge network
Before starting to add quantitative dynamic data to the already existing knowledge of progesterone signalling events in T47D cells, we have generated a “prior knowledge network (PKN)”, based on the published literature (Figures S1A, B, and Supplementary Material; Cytoscape Session 1). Each protein-protein interaction is characterized based on type (interaction, phosphorylation or dissociation) and is displayed as a unique edge. The corresponding literature is given in Figure S1B and within the Cytoscape session. The PKN already shows the key role played by kinases in response to progestins. First, progestins via ERa activate SRC1 that phosphorylates MAPKK1, that activates ERK1, that phosphorylates PGR, resulting in dissociation from the HSP90A and B proteins (36, 37). Activated ERK1 also phosphorylates ERa at S118 (38). ERK1 in association with hormone receptors translocates to the cell nucleus where it phosphorylates MSK1 (39), leading to the formation of an active complex PR-ERK-MSK1 that interacts with chromatin containing accessible PRE. Activated PR also interacts with PLK1 that activates MLL2 (40), with CDK2 that phosphorylates and activates ARTD1 (7), and with JAK2 that activates STAT5 (16). Simultaneously, membrane activated SRC1, also activates RAS and EGFR (41), which feeds back activating the MAPK cascade. Membrane associated ERa also activates PI3K and cAMP, which upon binding with AKT and PKA respectively lead to the activation (via interaction and direct phosphorylation) of GSK3, mTOR (42, 43) and the arginine methyltransferases CARM1 and PRMT1 within the nucleus (44–46). This brief description of the PKN shows that it already encompasses a great degree of complexity and complementary connections that need additional data to be resolved.
Microarrays of antibodies to phosphorylated sites in proteins
Our plan was to combine antibody microarray technology and shotgun phosphoproteomics in T47DM breast cancer cells exposed to 10 nM R5020 for different lengths of time, as previously described (7). For each experimental approach and exposure time total protein extracts were harvested in triplicate. For the antibody arrays, data was collected, filtered for quality control and summarized as log2 ratio over time zero (as described in materials and methods, Figure S2A). This dataset provided 246 unique phosphorylation sites corresponding to 155 proteins (Figure 1A; Supplementary Table S1). The majority of proteins contain 1 phosphorylation site, although for several proteins (Tau, RB1, MAP2K1, PTK2 and the protein kinase RPS6KA1) 7 or more significantly regulated phosphorylation sites were identified (Figure S2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g001.jpg)
Figure 1 Targeted Antibody Array Phosphorylation data following hormone. (A) Schematic overview of experimental procedure. Synchronized T47D cells were exposed to hormone for the length of time indicated. Triplicate samples were harvested and phosphorylated proteins were identified using an antibody microarray (see materials and methods). Data was log2 normalized resulting in a total of 246 significant phosphosites from 155 unique proteins. (B) Number of regulated sites per time point, log2FC (Fold change) >0.6<-0.6 versus time 0. (C) Breakdown of up (>0.6log2FC) and down (<-0.6 log2 FC) per time point versus T0. (D) Functional classification of the proteins identified as significantly phosphorylated across all time points, individual time point functional analysis (per time point see Figure S2C). (E) Venn diagram showing the overlap of significantly regulated phosphorylation sites across all time points. (F) Log2 FC following hormone of phosphorylated Mnk1 (T197/202) and EIF4E S209. (G) Expression of EIF4E in normal versus breast tumor samples from breast cancer patients (Protein Atlas, see materials and methods). (H) Kaplan Meyer overall survival stratifying patients based on the expression level of EIF4E in breast cancer data set (p=5.3e-10). (I) mRNA expression level of EIF4E in T47D cells treated with hormone. Dynamics of (J) p27/KIP T187 (K) CDK2 T160, ERK Y202/204. (L) Expression level of total PR and CDK2 and phospho-CDK2 T160 in T47D breast cancer cells exposed to hormone in the presence or absence of ERK inhibitor (ERKii) as determined by western blotting using specific antibodies. (M) CDK2 T14/Y15 phosphorylation in response to hormone as determined by antibody array. Dynamics of (N) Cdc25C S216 and (O) MAPKAPK2 T222 phosphorylation in response to hormone as determined by antibody array. (P). Model for CDK2/ERK dynamic activation and deactivation in response to hormone based on the data presented in (J–O).
Analysis of the number of phosphorylation sites clearly showed a rapid activation already 1-minute following hormone (68 significant phosphorylation events Figure 1B). Signalling persists throughout the time course showing two peaks at 30- and 360-minutes following hormone (Figure 1B). Phosphorylation sites were characterized as up or down-regulated, using a threshold for the log2 fold change with respect to time 0 of -0.6< or >0.6 respectively (Figure 1C). We see a trend for early phosphorylation sites to be dynamically increased compared to time zero, in contrast to later time points where protein phosphorylation sites as a whole decrease compared to time zero (Figure 1C). The majority of phosphorylation sites identified belong to protein kinases (45%), co-factors (11%), transcription factors (18%) and structural proteins (13%). (Figure 1D). Combining the identified phosphorylation sites over the time course reveals that the majority of sites are regulated at more than one time point (Figure 1E), however the protein function enrichment does not alter significantly over time, with kinases and transcription factors being the main protein groups where the phosphorylation sites are observed (Figure S2C).
Pathway and gene ontology (GO) for biological function (BP) molecular function (MF) analysis revealed a significant increase in Cancer pathways (Figure S2D), signal transduction, biopolymer metabolic process and kinase activity (Figures 2SE and F). Within this dataset we observed a strongly upregulated phosphorylation of the MAPK Signal-Integrating Kinase 1, MNK1 at T250/T255 (Figure 1F) in T47D in response to progesterone stimulation. Phosphorylation of MNK1 at T250/T255 by ERK induces the activity of MNK1 (47). Once activated, MNK1 phosphorylates its targets, including the proto-oncogene Eukaryotic Translation Initiation Factor 4E (EIF4E), for which we also observed a modest phosphorylation which follows a similar pattern to MNK1 (Figure 1F). Activation of MNK1 has been shown to promote cell proliferation thus MNK1 inhibitors appear as an exciting opportunity for cancer therapy. MNK1 signalling play a key role in invasive breast cancer growth (48), MNK1 inhibitors have been shown to block breast cancer proliferation in multiple cell lines (49), and its downstream target EIF4E is overexpressed in tumor versus normal samples from breast cancer patients (Figure 1G) and associated with a poor overall survival (Figure 1H). Our results are the first indication that MNK1 activation may be relevant for progesterone induced breast cancer cell proliferation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g002.jpg)
Figure 2 Phosphosite enriched Shotgun Proteomics following hormone. (A) T47D cells were treated with hormone at the times indicated. Biological triplicates were enriched for phosphopeptides using TiO2 followed by LC-MS-MS peptide identification. Data was log2 normalised resulting in a total of 310 phosphosites from 264 unique proteins. (B) Volcano plots showing phosphopeptide log2FC versus p-value for each of the time points following hormone. (C) Number of significant phosphosites identified per time point. (D) Analysis of the proportion of threonine, tyrosine and serine phosphorylated residues identified. (E) Breakdown of up (>0.6log2FC) and down (<-0.6 log2 FC) per time point versus T0. (F) Venn diagram showing the overlap of significantly regulated phosphorylation sites over time. (G) Phosphorylation of progesterone receptor PR (S162), following progesterone validated by western blotting in presence or absence of CDK2 inhibitor (H). Total PR levels are shown as a loading control.
CDK2 plays an important role in progesterone signalling, activating ARTD1, and phosphorylating histone H1 (7). CDK2 activity is controlled by the formation of an active complex with the cyclin partner; either Cyclin E or A. In addition to binding the cyclin partner, CDKs are also controlled via interactions with Kinase Inhibitory Proteins (KIPs). p27/KIP is rapidly dephosphorylated at T187 in response to hormone, dropping sharply at 1 minute after hormone exposure (Figure 1J), when CDK2 is phosphorylated and activated. Phosphorylation of p27 at T187 results in its ubiquitination and degradation and inhibits the interaction with CDK2 (50), which would result in the release of CDK2 from the inhibitory protein complex resulting in the activation of ARTD1 and subsequent nuclear effects. In addition, we observe the coordinated activation of the upstream kinase of CDK2 at T160; ERK at Y202/204 (Figure 1K) and could validate the phosphorylation of CDK2 T160 and its dependence on ERK activity by western blotting in the presence of ERK inhibition (Figure 1L). CDK2 at T160 is the active phosphorylation site of CDK2 peaking at 1-minute following hormone exposure (Figure 1K) in contrast to the inactive phosphorylation site of CDK2 (T14/Y15) which peaks at 60 minutes following hormone exposure to silence the kinase (Figure 1M). The activity of the phosphatase CDC25C is key for the removal of the inhibitory T14/Y15 phosphorylation sites of CDK2. The phosphatase itself is inactivated by phosphorylation at S216. We observe a peak in CDC25C phosphorylation prior to and following 60 minutes of hormone exposure, which would permit the phosphorylation of the inhibitory phosphorylation site in CDK2 (Figures 1M, N). Going one step further within this pathway; MAPKAPK2, the kinase which phosphorylates CDC25C at S216 is activated following the same time dynamic as its target (Figure 1O). Although the importance of CDK2 in progestin induced cell proliferation has been studied (7, 51) the complex mechanism of CDK2 activation; phosphorylation of active/inactive marks, activation and regulation of upstream phosphatases and kinases was not clear until now (Figure 1P). These examples of the dynamic phosphorylation of MNK1 and CDK2 highlight the insight that can be gained by this type of global signalling datasets.
Shotgun phosphoproteomics
To complement the microarray dataset, we performed shotgun phosphoproteomic analysis using mass spec. Phospho-peptides from T47DM cells exposed to 10 nM R5020 for the same duration as in the array experiments, were enriched using TiO2 and phosphorylated peptides identified by LC-MS-MS (Figure 2A; Figure S3A; Supplementary Table S1, S2). We identified 7482 phosphopeptides with a FDR 5% (Table S2) and dynamic, significant changes in 310 unique phosphorylation sites within 264 unique proteins (Figures 2B, C). The majority of proteins exhibited regulation of a single phosphosite, except for the serine/arginine repetitive matrix protein, SRRM1, involved in mRNA processing and the TP53 enhancing protein TP53BP1, that exhibited 8 and 10 regulated phosphorylation sites respectively (Figure S3B). Most phosphorylation sites identified were phosphor-serine consistent with the biological ratio of residue specific phosphorylation (Figure 2D). Over the time course, changes at each time point were identified as either up (log2FC>0.6) or down (log2FC<-0.6) regulated (Figure 2E). Up-regulated sites prevailed at early time points and many of these phosphorylation sites were significantly regulated at more than one time point (Figure 2F). Pathway and GO-BP (Biological Process) and MF (Molecular Function) enrichment analysis was consistent with the antibody array enrichment and revealed an increase in pathways in cancer, biopolymer metabolic process and kinase activity (Figures S3C–E).
PR S294 is rapidly phosphorylated in response to hormone resulting in its activation and dissociation from chaperone complexes and increased protein turnover (52). In recent years it has been shown that clinical samples assigned as “PR low” actually have elevated levels of phosphorylated PR S294 and that this phosphorylation is associated with a genetic signature linked to cancer stem cell growth and increased recurrence which may have implications for the treatment of PR low patients with anti-progestins (53). Phosphorylation of PR S162 in the hinge region showed a strong hormone induced increase by mass spec (Figures 2G, H). Phosphorylation within this region of PGR has been previously reported to be mediated by CDK2 (54), which we were able to confirm as the specific phosphorylation of S162 PR in response to progesterone was strongly decreased in the presence of CDK2 inhibition (Figure 2H).
Bioinformatic functional analysis
In order to investigate the dynamics of progestin signalling over time and with the aim of avoiding inherent biases generated from either technical approach, we combined the significantly regulated phosphorylation sites from both datasets (Figures 1, 2) resulting in a list of 420 unique phosphorylation sites within 390 proteins (Figure S4A). PCA analysis of the samples reveals a clear separation of the phosphorylation data at 6 hours following the hormone, given the majority of phosphorylation sites are rapid, this separation of the latest time point may reveal changes in protein abundance at this time point. The majority of these proteins showed the regulation of a single phosphorylation event with the exception of several highlighted proteins, including FAK, MAPT and EGFR (Figure S4C). As in the individual analysis, phosphorylation sites were significantly regulated over several time points (Figure S4D) and showed a switch from up-regulated sites early after hormone exposure to down-regulated sites at later time points (Figure S4D). KEGG pathway analysis shows a significant enrichment in MAPK, PI3K-AKT, neurotrophin (TRK) and ERBB signalling pathways (Figure 3A; Supplementary Table S2), in addition to pathways key in the progression of cancer, specifically cancer stem cells, such as focal adhesion (Figure 3A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g003.jpg)
Figure 3 Combining Target Antibody Arrays and Shotgun Phosphoproteomic datasets following hormone. (A) KEGG pathway enrichment analysis of proteins identified as regulated by phosphorylation in response to hormone. (B) Cellular component analysis of phosphosites enriched per time point. Showing the hormone induced phosphorylation of the nucleoplasm and cytosol across all time points (group I) the activation of membrane raft proteins enriched at 1 minute (Group II) and phosphorylation of mitochondrial proteins enriched at 15 minutes (Group III), activation of nuclear structures; PML bodies and the nuclear matrix at 60 minutes (Group IV) and the activation of the cell-cell junctions and microtubules at 360 minutes (Group V). (C) K mean clustering of all significantly regulated phosphorylation sites over time reveals 6 distinct clusters. (D) Similarity matrix of clusters 1-6 reveals similar dynamics for clusters 1 and 4 and an opposing similarity in phosphorylation dynamics for clusters 3 and 6. Red indicates highly similar, well correlated, blue inversely correlated patterns of regulation. (E) Word cloud showing the enrichment of GO-biological processes associated with proteins identified in similar clusters 1 and 4 “Early risers” which are regulated rapidly after hormone. (F) Graph showing the opposing phosphorylation dynamic of proteins within clusters 3 and 6. (G) Venn diagram showing the overlap of significantly identified Corum protein complexes identified in clusters 3 and 6. (H) Phosphorylation dynamic in response to hormone of Myc S373, and T58 and BAF53 S233. (I) Expression level of BAF53 in tumor versus normal tissue within the TGCA dataset. (J) Phosphorylation of AMPL T183 decreases rapidly in response to hormone. (K) Model showing the key role of AMPK dephosphorylation in response to hormone in breast cancer cells, AMPK dephosphorylation is required in order for subsequent signalling cascades including NFkB, insulin, Hippo, JAK/STAT and mTOR to continue and the phosphorylation of PR S294 to take place.
GO cellular component analysis reveals a dynamic pattern of specific cellular compartments over time (Figure 3B; Supplementary Table S3). As expected, over the whole-time course, proteins are mainly found within the cytosol and nucleoplasm. However, prior to hormone exposure, phosphorylated proteins are enriched in RNA transcription repression complex and nuclear chromatin. The addition of hormone rapidly induces the phosphorylation of the membrane rafts, components of focal adhesion and protein kinases consistent with published works whereby signalling initiates from the plasma membrane. This transient phosphorylation of the membrane rafts diminish after 1 minute and is followed by the phosphorylation of transcription factors and proteins within the cytoskeleton (Figure 3B). Interestingly, in line with our findings showing the generation of nuclear ATP synthesis independent of mitochondrial supplementation at 30 minutes after hormone we observe an enrichment in phosphorylated proteins located within the mitochondrial membrane at 15 minutes (Figure 3B). The dynamic regulation of these proteins; CYB5B (cytochrome b5), the transcriptional activator ATF2 (Cyclic AMP-dependent transcription factor ATF-2), RPS6KB1 (Ribosomal protein S6 kinase beta-1) and PI4KB (Phosphatidylinositol 4-kinase beta) (Figure S4F) may suggest an as yet undiscovered crosstalk between the nuclear and mitochondrial ATP synthesis pathways. Mitochondrial PR (PR-M) is a truncated isoform of the nuclear progesterone receptors PRB and PRA, which lacks the N-terminal DNA binding domain present in PRA and PRB but does contain the hinge region responsible for dimerization and the ligand binding domain (55). PR-M has been shown to increase cellular respiration hence cell energy levels in response to ligand in various physiological situations and animal models (56). Therefore, the coordinated phosphorylation of proteins within the mitochondria in response to ligand (Figure S4F) may provide an interesting insight into a possible crosstalk between mitochondrial PR-M and the nuclear receptors PRA and PRB.
At 60 minutes following hormone exposure the main localization of phosphorylation changes and shifts again to nuclear matrix proteins and proteins found within distinct regions of the nucleus, such as PML bodies (Figure 3B group IV), which may be involved in the reorganization of chromatin in response to progestins (22). At 6 hours following hormone exposure cells enter the early stages of the cell cycle and cellular movement is increased. This is also evident by the enrichment of phosphorylation sites in proteins within cell-cell junctions, the cytoskeleton and microtubules (Figure 3B group V).
Protein class analysis
The majority of identified phosphorylated proteins (60%) were assigned to one class, however due to the promiscuous nature of enzymes nearly 40% were assigned to more than one class (Figure S4G). Taking only the parent class into account, we observed 5 distinct functions; 1) nucleic acid binding, 2) enzymes, 3) structural proteins, 4) protein modulators, and 5) proteins involved in signalling, membrane and cell-cell contacts (Figure S4H). Each function class consists of sub-groups (Figures S5A–F). The Nucleic Acid binding class includes DNA binding proteins, helicases, nucleases and RNA binding protein subgroups (Figure S5A). The enzyme class is dominated by kinases but also includes histone modifying enzymes, hydrolases, ligases and oxidoreductases (Figure S5D). The structural class is dominated by cytoskeleton proteins (Figure S5F). The protein modulator class includes chaperones and various kinases and G proteins regulators (Figure S5C). The cell signalling and the membrane/cell-cell contact classes are more complex and include many specialized proteins such as signalling molecules, receptors and transporters (Figures S5B, E).
Gene Ontology of Biological Processes (GO-BP) and Molecular Function (GO-MF) showed an enrichment in signal transduction and general biological processes across the entire time course (Figures S6A, B; Supplementary Tables S4, S5). However, several interesting dynamic functions were identified. For instance, transcription cofactors, transcriptional repressors and transcription factor binding were already enriched 1 minute after hormone exposure (Figure S6B) consistent with our previous observations of rapid transcription factor recruitment following hormone exposure (9, 19). We observed enrichment in ATP binding and Adenyl-ribonucleotide binding after 5 and 60 minutes of hormone exposure (Figure S6B), which may represent regulation of the two cycles of ATP dependent chromatin modifiers in response to progesterone (21, 57). KEGG pathway analysis reveals a significant enrichment in signalling and in many cancer pathways, including Prostate, Glioma, CML, lung, AML, endometrial and pancreatic cancer, as well as focal adhesion and tight junctions (Figure S6C). Annotated signalling cascades were significantly enriched at all time points in response to progestin, including MAPK, neurotrophin (TRK), ERBB, FC-receptor and insulin signalling (Figure S6C).
Specific pathways: roles of AMPK, insulin TNFa, and PIK3
K Means clustering analysis revealed six patterns of regulation over the time course (Figure 3C). Similarity analysis of all phosphorylation sites within all clusters shows several interesting dynamics. First, “Early-risers” cluster 1 and 4, are positively correlated on the similarity matrix and show their initial increase in phosphorylation early at 1 and 5 minutes, respectively (Figure 3D). GO-BP analysis of the proteins contained within the “early riser” clusters shows an enrichment in signal regulation, and signalling cascades including Hippo, NFk-B and MAPK pathways (Figure 3E; Supplementary Table S9). Second, clusters 3 and 6 show an opposing nature (negative correlation Figure 3D). This antagonistic behavior of the two clusters is clearly shown averaging the signal of all phosphorylation within each cluster (Figure 3F). Corum (comprehensive resource of mammalian protein complexes) analysis of the significantly enriched protein complexes contained within clusters 3 and 6 (Supplementary Table S6) showed that most protein complexes were enriched in one cluster or the other (Figure 3G), likely representing crosstalk. Ten protein complexes were found to be enriched in both cluster 3 and 6, having phosphorylation sites within the same protein complex regulated in an opposite manner (Figure 3G).
One such complex was the cMyc-ATPase-Helicase complex, which contains 5 proteins; cMyc, the chromatin remodelling component BAF53, the ATP-dependent helicases RUVBL1 and RUVBL2 (also known as TIP48 and 49) and the histone acetyltransferase, TRRAP. This complex is involved in chromatin organization, histone acetylation and transcriptional regulation (58). Analysis of the phosphorylation sites showed that two sites (S373, T58) within Myc were increased early after hormone exposure, and decreased after 60 minutes, whereas one site of BAF53 (S233) shows the opposite dynamic (Figure 3H). Database analysis also reveals a strong overexpression of BAF53 in tumor versus normal samples in multiple cancer types (Figure 3I). Myc has an important role in breast cancer growth via the activation of AMPK (59).
The AMP-activated protein kinase (AMPK), exhibited a decrease in T183 phosphorylation in response to hormone. This site is phosphorylated by CAMKK1 or 2 (60). AMPK is a master sensor, and its activation inhibits several kinase pathways including mTOR, NfkB, JAK/STAT, insulin and Hippo (61–64). In addition, active AMPK inhibits the phosphorylation of PR S294, PR recruitment to chromatin and the activation of progesterone regulated genes (65). Activation of the kinase, specifically requires the phosphorylation of AMPK at T183 by CAMKKs, and de-phosphorylation of this site has been shown to be induced by estrogens and androgens in adipocytes (66, 67). Previous published results and the data presented here suggests a model where AMPK must be silenced in order for regulatory pathways described and PR itself to be active, which is what we observe within 1 minute of progesterone stimulation (Figures 3J, K).
Further in-depth analysis of the complexes which are regulated by phosphorylation in response to progestin revealed a full list of complexes with at least 2 proteins phosphorylated in response to hormone. The Sam68-p85 P13K-IRS-1-IR signalling complex, which encompasses the insulin receptor (INSR), the insulin receptor substrate 1 (IRS1), the KH domain containing transduction-associated protein 1 (Sam68) and the phosphatidylinositol 3-kinase regulatory subunit alpha (GRB1). This protein complex is involved in insulin signalling and has been proposed to provide a link between the PI3K pathway and other signalling cascades of insulin or p21/RAS (68). We observed a dynamic phosphorylation of several sites within the complex (Figure 4A), including 4 distinct phosphorylation events within IRS1, two of which peak at 1 minute (S312, S639) and two sites where the peak in phosphorylation is observed at 60 minutes (Y1179, and Y612). S312 has been shown to be directly phosphorylated by c-Jun N-terminal kinase (JNK1) in breast cancer signalling (69) and this phosphorylation inhibits its interaction with IKKA. S639 is phosphorylated by mTOR has been linked to PI3K/Akt/mTOR signalling in breast cancer (70, 71) and effects the intracellular localization of IRS1 (72). Y1179 has been reported to be phosphorylated by IGF1R or INSR itself (73) and Y612 phosphorylation activates the interaction with PIK3R1 (74).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g004.jpg)
Figure 4 Protein Complex analysis and Overlap of PARylation and Phosphorylation in response to Progesterone. Heatmaps showing the phosphorylation of proteins within the Sam68-p85 IRS (A) and NF-kappa B (B) signalling complexes in response to hormone over time. (C) Heatmap showing the phosphorylation of proteins of the p130 Cas-ER-Src-PI3K complex in response to hormone over time, the coordinated phosphorylation of each phosphosite individually is represented as a line graph (right panel). (D) Venn diagram showing the overlap of proteins which contain either a phosphorylation site (379), PARylation site (1187) or both PTMs within the same protein after hormone exposure in T47D cells (52). (E) Word cloud representation showing the GO-cellular component enrichment analysis of the 52 proteins identified as phosphorylated and PARylated in response to hormone (D). (F) Heatmap showing the phosphorylation of components of the KRS1-RAF1-MEK signalling complex in response to hormone over time, all proteins shown are phosphorylated and PARylated and the dynamics of individual sites is shown on the right panel. (G) Heatmap showing the phosphorylation proteins of the Emerin complex in response to hormone over time. (H) Heatmap showing the phosphorylation of proteins (5/15) within the Kinase Maturation complex in response to hormone over time, phosphorylation dynamics of individual sites is shown (lower panel). (I) Venn diagram showing the phosphorylation and or PARylation of proteins contained within the kinase maturation complex; 14/15 protein components of the complex contain at least one of the PTMs. (J) Schematic representation of the complex components, PARylated proteins are indicated by blue circle, Phosphorylated by red (right panel). *Indicates proteins which are also PARylated.
The phosphorylation of several components of the TNFa/NFkB signalling complex were also identified (Figure 4B). This complex is involved in I-kB kinase/NF-kB signalling in tumor progression. Indeed, complex components IKKα, RelB and p52 are associated with decreased cancer-specific survival in ERa-positive breast cancer (75). This may be linked to the cancer stem cell niche, which we showed recently was present in T47D cells grown in 3D cultures (24). NFkB regulates self-renewal in breast cancer stem cell (BCSC) models and deletion of IKKα in mammary-gland epithelial cells affects progestin-driven breast cancer (76, 77). Indeed, the upstream activator RANK ligand (RANKL) and hence the RANK pathway promotes mammary tumor formation, (78), (76). Another example is the P130Cas-ER-cSrc-PIK3 kinase complex (Figure 4C). Which has been shown to induce transcriptional changes in response to oestrogen and mammary proliferation in breast cancer. The authors showed that estradiol triggers the association of ERa, c-Src, the p85 subunit of PI 3-kinase (PI3K) and p130Cas in a macromolecular complex and activates the c-Src kinase leading to p130Cas-dependent Erk1/2 phosphorylation (79, 80). Given the similarity of the phosphorylation dynamics, peaking early at 5 and 15 minutes across Src, PIK3 and ESR1 within the complex induced by progestin (Figure 4C right panel) this may (similarly to the induction by oestrogen shown by others) present a novel ERa-ERK-cSrc activation mechanism in response to progestin in breast cancer cells.
Crosstalk between progestin induced phosphorylation and PARylation
Recently, it has been shown that the majority of PARylation events on eukaryotic nuclear proteins take place on serine residues rather than acidic residues as previously accepted (81–86). Given the enrichment of serine in the phosphorylation dataset (Figure 2D) and the importance of both PTMs in progestin gene regulation (21), we investigated the overlap between PARylation sites and phosphorylation sites within protein complexes. We identified 52 proteins, which were both phosphorylated and PARylated in response to progestins in T47D breast cancer cells (Figure 4D). Cellular component analysis of this set of 52 proteins indicates a significant enrichment in nuclear, cytoskeleton and chromatin associated proteins (Figure 4E; Supplementary Table S7), in line with the well described role of PAR in the nucleus, nuclear organization, chromatin organization, relaxation and transcriptional regulation (87–90). This finding may indicate a crosstalk between PARylation and phosphorylation with regards to nuclear structure and chromatin organization. Analysis of complexes significantly enriched within this group of proteins revealed 12 protein complexes (Supplementary Table S8), which contained proteins both PARylated and phosphorylated. One of them is the KSR1-RAF1-MEK complex composed of MEK1 and 2, both PARylated and phosphorylated, and RAF1 which is phosphorylated (Figure 4F). This complex is involved in the MAKPKKK cascade, and in response to EGF it activates BRAF mediated phosphorylation of MEK1, at 3 sites, and MEK2, which activate MAPK1 and 3 (91). In our dataset we observe a clear change in phosphorylation of all members of the complex in response to progestin (Figure 4F right panel).
We also observed the phosphorylation and PARylation of the Emerin complex (Figure 4G). This complex is involved in DNA replication, transcription and structural integrity of the nucleus, specifically of the inner nuclear membrane (92). Depletion of Emerin results in changes in the organization and dynamics of the nucleus, increased chromatin mobility and a mis-localization of chromosome territories (93). Within this complex we find proteins phosphorylated, PARylated, or phosphorylated and PARylated (Figure 4G). Given the role of PAR in the structure of the nucleus, this complex may present an interesting example for studying the PARylation, phosphorylation crosstalk.
As discussed, prior to hormone exposure PR is present in an inactive complex with the HSP70 and 90 proteins as part of the Kinase Maturation Complex. We know that progestins promote the phosphorylation and dimerization of the receptor and we found that phosphorylation of the HSP90 and 70, along with other members of the complex, is initiated within 1 minute of hormone exposure (Figure 4H), again showing a rapid and concerted phosphorylation of several members of the complex (Figure 4H). In addition, the HSPs are also PARylated as compared to other components where only phosphorylation (MARK2, MAP2K5) or PARylation (14-3-3 components) are present (Figures 4I, J). Further investigation regarding the crosstalk between PARylation and phosphorylation within protein complexes will be the focus of future studies.
Progesterone signalling network generation
In order to understand the crosstalk between the signalling pathways activated by Progesterone in breast cancer cells, a Protein-Protein Interaction (PPI) network was generated using all identified phosphorylation sites (Supplementary File: Network session 2), based on known PPI (evidence based). The resulting network consists of 427 nodes (proteins) and 4309 unique interactions (edges) (Figure S7A). Pathway analysis was carried out on this network, using Genemania™, and 23 statistically significant (p<0.01) pathways were identified (Supplementary Tables S9 and S10). The proteins and interactions (nodes and edges) associated with each pathway were selected and new networks generated (Supplementary File: Network session 2). Several of which are shown in Figures S7B–G and discussed briefly below.
One such pathway, the Fc receptor signalling pathway was enriched (Figure S6C) and the PPI network is shown in Figure S7B. Fc receptors are cell surface proteins that recognize the FC fragment of antibodies, mainly on immune cells. However, recent studies have shown that different subsets of Fc receptors may play a role in tumor cells (94). In particular, it was shown that T47D cells express the FcyRI (CD64). These FC-receptor expressing breast cancer cells can activate the tyrosine kinase signal transduction pathway. Indeed, T47D cells treated with selective tyrosine kinase inhibitors do not proliferate in a FC receptor- tyrosine kinase signalling dependent manner (94).
As mentioned before (Figure S6C), another pathway identified as activated in response to progestin is the ERBB-EGF network (Figures S7D, H). ERBB2 (HER2) is overexpressed in 15-20% of breast cancer in response to EGF activation, and plays a major role in EMT (95). PR interacts with ERBBs and induces the translocation of ERBB2-PR-STAT3 complex to the nucleus. ERBB2 acts as a co-activator of STAT3 and drives the activation of progestin regulated genes, especially genes such as Cyclin D1 that do not contain HREs (96, 97). Blocking PR signalling in PR-ERBB2 positive breast cancer patients has been suggested as a treatment (98). The ERBB pathway may represent a new mechanism for further study to understand the activation of these “non-classical” PR dependent genes in response to progestin. In addition, to the role of ERBB2, the role of ER activation in response to progesterone in breast cancer cells is also critical, as shown in Figure 4C we observe the coordinated activation of the ESR1-Src-PIK3 complex peaking at 15 minutes following hormone exposure. The phosphorylation site of ESR1 which increases is S104. ER S104 phosphorylation is essential for ER activity (99) and it has been suggested that hyperphosphorylation of ER at these sites may contribute to resistance to tamoxifen in hormone receptor positive breast cancer (100–102). ER S104 phosphorylation by ERK has been shown previously in response to oestrogen and EGF but not progesterone exposure. In addition, ER S104 has been implicated in mTOR signalling (103). Given the phosphorylation of ER, the dynamic activation of the ER membrane complex and the role of mTOR in AMPK and insulin signalling described earlier (Figures 3J and 4C) this phosphorylation site may present a key step in the cellular response to progesterone in breast cancer.
A pathway exhibiting strong activation by progestins is the Insulin signalling (Figure S6C; Figure S7F). Insulin-like growth factors (IGFs) and progestins both play a major role in normal mammary gland development and R5020 has been shown to induce the expression of insulin receptor substrate-2 in MCF7 cells (104, 105). Moreover, IGF signalling via IRS2 is known to be essential for breast cancer cell migration. It has also been shown that R5020 pretreatment followed by IGF stimulation increases binding of IRS to PI3K-p85 regulatory complex, which in turn activates ERK and AKT signalling (106). Interestingly, not only do we observe the activation of the insulin pathway in network analysis (Figure S7F, but also the coordinated phosphorylation of all members of the IRS-PIK3 complex (Figure 4A), indicating that indeed progesterone stimulation of T47D cells activates not only the insulin pathway in general terms but also triggers the coordinated regulation of complexes within it.
Discussion
The data presented here provides a source of knowledge for the scientific community with regards to progesterone induced gene expression, and the signalling pathways involved. We have shown the rapid induction of phosphorylation using two distinct technologies (Figures 1, 2). Pathway analysis showed a strong enrichment in pathways associated with cancer, known and novel Pg-dependent signalling events (Figure 3A). But also identified signalling pathways not previously known to mediate progesterone action in breast cancer cells, such as MNK1/EIF4E pathway and the connection between CDK2, Cdc25 and the MAPK pathway.
Cellular component analysis confirmed our expectations and the statistically significant activation of the cell membrane within 1 minute of hormone exposure (Figure 3B), but also revealed a consistent (over all members) phosphorylation peak within proteins associated within the mitochondria at 15 minutes after hormone exposure (Figures 3B and S4E). Mitochondrial activation in response to progesterone in breast cancer cells has not been extensively studied yet. Indeed, ATP synthesis 45-60 minutes after hormone stimulation is independent of mitochondrial involvement (21). However, there are some interesting findings in the literature. Following the observation that the PR negative cell line MCF10A exhibits a progestin-induced cell proliferation (107). Behera and colleagues showed that MCF10A responded to R5020 with an increase in mitochondrial activation (108). Given the absence of the nuclear PR in these cells they hypothesized that the activation of progestin-induced cell growth was due to non-genomic metabolic effects, mediated by a yet undiscovered receptor. We propose that the observed mitochondrial activation in T47D in response to progestin (Figures 3B and S4E) suggests the existence of a third and interconnected hormonal signal transduction pathway via the mitochondria (109, 110).
Dynamic phosphorylation analysis over time reveals distinct groups or clusters of phosphorylation events which follow a similar time response; such as early risers, sustained or late (Figure 3C). Similarity analysis of these dynamic phosphorylation sites reveals some interesting crosstalk between protein complexes not previously identified as players in progesterone signalling (Supplementary Table S6), and complexes where a mobilization of phosphorylation (showing similar dynamics) was observed within the whole macromolecular complex; such as PIK3, NFkB (Figures 4B, C).
Overlap of phosphorylation sites with existing PARylation, revealed 52 proteins for which both phosphorylation and PARylation was found. The data also clearly shows an enrichment in protein complexes that play a role in the structural organization of the nucleus (Figures 4D, E), specifically the Emerin complex and Lamin (Figure 4G). The key location of these complexes at the nuclear membrane suggests that perhaps these two PTMs may affect and play a role in the dynamic structure of the nucleus. This could be explored in the future by global chromatin proximity Hi-C experiments. The complexes identified in this study and the dual post translational modification of proteins with known important roles within the cell may provide exciting opportunities for future studies which aim to understand the crosstalk between Serine PARylation and phosphorylation in the context of nuclear architecture, signalling and breast cancer progression (Supplementary Table S7).
In addition to pathway analysis at the single network level (Figure S7), wherein we identified pathways such as insulin, Fc-receptor and ERBB signalling, it is also clearly important to consider the connection between pathways and networks as a whole. One such example is the connections between the phosphorylation events within the cytoskeleton, membrane raft and proteins associated with cell adhesion. We observe phosphorylation events in multiple proteins within both cell adhesion and the membrane raft, forming tight strongly connected PPI networks (Figure 5A). A network merge of these two pathways reveals 4 key proteins which are present in both (JAK2, SRC1, LYN and KDR), indicating a strongly connected network (Figure 5B), which in addition to common members exhibits a large first neighbor selection between the two initial pathways (selection of only direct PPI) (Figure 5B right panel) with a similar phospho kinetic pattern (Figure 5C). Incorporation of the significantly enriched phosphorylated proteins within the cytoskeleton (Figure 3B; Figure S6D; Supplementary Table S3) into the merged network (Figure 5B) results in a larger global connected network (Figure 5D) which supports the activation within the membrane proteins after 1-minute following hormone exposure that triggers the subsequent cascades of phosphorylation in the cytoskeleton (Figure 5E). These findings clearly show the importance of studying the pathways not in isolation, but rather in connection with each other.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g005.jpg)
Figure 5 Combining PPI networks from distinct cellular compartments reveals a coordinated crosstalk. (A) PPI network showing the significantly regulated phosphorylated proteins located in cell adhesion (blue) and the membrane raft (red) identified in response to hormone. (B) Merge of the cell adhesion network (A) blue and membrane raft network (A) red. The two networks connect based on known PPI however no protein was identified as annotated in both sets. This integration of the two networks is highlighted (right panel) where proteins from each network were selected based on having a first neighbor with a protein of the other network. (C) Violin plot showing the average phosphorylation of proteins over time in response to hormone within the membrane raft or cell adhesion networks. Data is normalised to time 0 = 1. (D) Merge of Cell Adhesion-Membrane (B) and the cytoskeleton networks. The two networks are merged based on known PPI. Proteins annotated in more than one function are colored based on the Venn diagram (i.e. cytoskeleton and cell Adhesion; light green, membrane raft and cytoskeleton; brown). (E) Heatmap showing the average phosphorylation of all proteins within each network in response to hormone over time, showing the activation of the membrane raft first at 1 minute followed by the cytoskeleton and cell adhesion.
Other examples of connectivity are observed between the ERK subgroup and the MAPK cascades and the FC Receptor and TRK signalling. ERK and MAPK form a strong network (Figures 6A, B). As discussed, earlier Fc-receptor signalling shows a strong activation (Figure S7B). The tropomyosin receptor tyrosine kinases (TRKs) are primarily known for their roles in neuronal differentiation and survival. However, increasing evidence shows that TRK receptors can be found in a host of mammalian cell types to drive several cellular responses (111, 112). Aberrations in TRK signalling, which can occur through events such as protein overexpression, alternative splicing, or gene amplification, can lead to disease such as cancer (113–115). The receptor tyrosine kinase NTRK2, activates GRB2-Ras-MAPK cascade in neurons and increases secretion by epithelial cells in culture in response to oestrogen or progestin treatment and NTRK2 was identified as differentially expressed between stromal and epithelial breast cells which may have implications in invasion and metastasis (116). Merging of Fc-receptor and TRK signalling pathways (Figure 6C) shows a strong protein overlap and a dense connected network with FOXO1 at the center (Figure 6D). FOXO1 is phosphorylated after 1 to 5 minutes of progestin exposure (Figure 6E). Phosphorylation of FOXO1 by PKB/Akt has been shown to be important for the binding to 14-3-3 proteins on chromatin (117–120). The role of these two pathways in progesterone induced gene regulation has not been shown previously. FOXO factors have a key role to play in tumor resistance to therapy and patient outcome (121). Interestingly, from a clinical perspective, stratifying patients based on either the expression levels of NTRK2 (TRKB) or FOXO1 is predictive of a good prognosis (overall survival) in breast cancer, similar to prognosis based on PR expression (Figures 6F, G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/888802/fendo-13-888802-HTML/image_m/fendo-13-888802-g006.jpg)
Figure 6 Network Integration of signalling networks identified in response to hormone. (A) PPI network showing the phosphorylated proteins present within the ERK signalling cascade (green) and the MAPK cascade (blue) identified in response to hormone. (B) Merge of ERK-MAPK networks (A). The two networks are merged based on known PPI. Proteins annotated in both pathways are colored based on the Venn diagram (fuchsia). (C) PPI network showing the phosphorylated proteins present within the FC-receptor (yellow) and TRK-neurotrophin (red) signalling pathways (yellow) identified in response to hormone (left and middle panel). (D) Merge of FC-receptor and TRK neurotrophin networks. The two networks are merged based on known PPI. Proteins annotated in both pathways are colored based on the Venn diagram (orange). (E) Rapid and coordinated phosphorylation of FOXO1 S256 and FOXO S319 in response to hormone. Kaplan Meyer overall survival of patients stratified based on the expression of NTRK2 (F) and FOXO1 (G) in breast cancer patients (p=9.4E-10 and 6.8e-12 respectively). All networks, PPIs and integrated cascades are supplied in Cytoscape session 2.
The examples described here in addition to other examples contained within the data for future discovery show the importance of network connectivity in trying to understand not only individual proteins or pathways but the significant overlap between the pathways within the signalling network activated by progesterone in breast cancer. Further analysis of the detected connections and identification of the key regulators may provide a source of targets for drug discovery aiming at the treatment of hormone receptor positive breast cancer patients.
Conclusions
The data presented here, reveals a high level of complexity in progesterone signalling in T47D breast cancer cells, shedding new light on known proteins and signalling pathways. Functional analysis reveals the activation of known pathways such as MAPK cascade but also the activation of signalling cascades not previously associated with progesterone signalling such as TRK, TNFa and ERBB. Our analysis indicates that there is a full cellular coordinated response, with proteins activated in different cellular compartments at different times following hormone exposure, in addition to the activation of whole protein complexes previously not associated with progesterone signalling. One limitation of our study is the use of a single breast cancer cell line, future work will aim to address this limitation with the use of additional cell lines. However, we believe that this signalling network and the phosphosites identified represent a rich resource for the breast cancer research community, opening up new lines of research and ideas for possible drug discovery projects for the benefit of breast cancer patients.
Data availability statement
‘The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Authors contributions
Experimental design, RW and MB. Bioinformatic analysis, JO, JC-C, and RW. Manuscript writing and editing, RW and MB. Experiments, RW. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by European Research Council (Project “4D Genome” 609989), the Ministerio de Economía y Competitividad (Project G62426937) and the Generalitat de Catalunya (Project AGAUR SGR 575 and AGAUR 2019PROD00115/IU68-016733), European Research Council -Proof Of Concept (Project “Impacct” 825176).
Acknowledgments
We acknowledge the support of all members of the Chromatin, Gene Regulation laboratory and members of the Gene Regulation Cancer and Stem Cells department at Centre for Genomic Regulation (CRG, Barcelona Spain) and the CRG/UPF proteomics unit.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.888802/full#supplementary-material
Supplementary Figure 1 | Prior Knowledge Network (PKN) Progesterone Signalling. (A) Edge directed PKN network was manually curated from the literature. Annotated phosphorylation events, interactions, dissociations and cellular compartment are indicated. Network is available as a cytoscape network session (cys) or interaction (sys) file containing references for all edges present as shown in (B) (See Supplementary File Network 1).
Supplementary Figure 2 | Antibody Array controls and data analysis. (A) Schematic indicating the experimental procedure, quality control checks, and filtering applied to the antibody array experiments. (B) Number of phosphosites identified per protein. Tau, PTK2, RPS6KA1 and RB1 are highlighted as they have multiple sites identified. (C) Functional classification of the proteins identified as significantly phosphorylated at each time point (D) KEGG pathway analysis, showing significant pathways (-log10 p-value) at each time point. (E) Heatmap representation of GO biological process data, showing significant (-log10 p-value) processes at each time point. (F) Heatmap representation of GO molecular function data, showing significant (-log10 p-value) functions at each time point.
Supplementary Figure 3 | Phosphoproteomic data acquisition and controls. (A) Correlation of triplicate samples from each of the time points. (B) Number of phosphosites identified per protein, the proteins showing multiple sites per protein are highlighted. (C) KEGG pathway analysis, showing significant pathways (-log10 p-value) at each time point. (D) Heatmap representation of GO-biological process data, showing significant (-log10 p-value) processes at each time point. (E) Heatmap representation of GO-molecular function data, showing significant (-log10 p-value) functions per time point.
Supplementary Figure 4 | Combining Antibody Array and Phosphoproteomic LC-MS-MS datasets. (A) Schematic representation showing the methodology and overlap combining antibody array and LC-MS-MS datasets. (B) PCA analysis of phosphorylation datasets. (C) Number of phosphosites identified per protein, the names of proteins showing multiple sites per protein are highlighted. (D) Venn diagram showing the overlap of phosphosites per time point. (E) Up and down regulated phosphorylation sites identified per time point. (F) Phosphorylation levels of the proteins identified as significantly regulated after hormone located within the mitochondria. (G) Analysis of the number of functions to which each unique protein was assigned (H) Venn diagram showing the overlap of protein functional class; Enzymes, Structural protein, Membrane-cell-cell contact, protein modulators and proteins with nucleic acid binding capacities.
Supplementary Figure 5 | Functional Analysis of Proteins Identified. All proteins were assigned one or more function based on the GSEA database. Both Parent (outside/title), and children (within) are shown for each class and the proteins identified within that sub-group are shown. Nucleic acid binding (A), Membrane/Cell-cell contact (B), Protein Modulators (C), Enzymes (D), Cell signalling (E), and Structural proteins (F).
Supplementary Figure 6 | Gene Ontology and Pathway analysis of combined dataset. (A) Heatmap representation of GO biological process data, showing significant (-log10 p-value) biological processes enrichment based on the protein phosphorylation at each time point. (B) Heatmap representation of GO molecular function enrichment, showing significant (-log10 p-value) functions at each time point following hormone. (C) KEGG pathway analysis, showing significant pathways (-log10 p-value) enriched at each time point following hormone exposure. (D) Protein protein interaction (PPI) network generated using proteins identified as phosphorylated following hormone and were assigned as cytoskeleton located.
Supplementary Figure 7 | Pathway Network Generation in Breast Cancer cells in response to Hormone. (A) Protein protein interaction (PPI) network was generated using a full phosphorylation dataset encompassing 321 proteins (Supplementary Material Network session 2) in Cytoscape using Genemania™ only considering protein-protein interactions with experimental evidence (Supp. Materials and methods), each node represents and individual protein and interactions are represented by edges. Functional analysis was carried out to identify key pathways enriched within the full network (Full list Supplementary Table 10). Individual networks were generated from each function individually and are available within additional Network session 2. Graphs of several pathways determined to be enriched within the dataset are shown (B) Fc receptor, (C) MAPK, (D) EGF (E) ERK (F) Insulin, (G) TRK signalling, (H) ERBB.
Supplementary Table 1 | Uniprot IDs of phosphorylated proteins identified in response to hormone. Time after hormone (minutes), data is normalised 0-1 row maximum and minimum.
Supplementary Table 2 | KEGG pathway enrichment; the pathway term, p value and the proteins associated with the pathway are shown.
Supplementary Table 3 | Cellular component enrichment analysis of phosphorylated proteins. The time after hormone in which they peak, the adjusted p value, and proteins associated with each specific cellular component are given.
Supplementary Table 4 | Gene Ontology Biological Process enrichment analysis of phosphorylated proteins. The cluster in which the term is enriched, the adjusted p value, and proteins associated with each specific biological process are given.
Supplementary Table 5 | Gene Ontology Molecular Function enrichment analysis of phosphorylated proteins. The cluster in which the term is enriched, the adjusted p value, and proteins associated with each specific molecular function are given.
Supplementary Table 6 | Corum enrichment analysis of phosphorylated proteins. The p-value, and proteins associated with each complex are given.
Supplementary Table 7 | Cellular component enrichment analysis of phosphorylated and PARylated proteins. The adjusted p value, and proteins associated with each specific cellular component are given.
Supplementary Table 8 | Corum enrichment analysis of phosphorylated and PARylated proteins. The p-value, and proteins associated with each complex are given, phosphorylated proteins are highlighted in yellow.
Supplementary Table 9 | Genemania analysis of phosphorylated proteins, all protein IDs are listed along with the GO: IDs for which they are associated.
Supplementary Table 10 | Genemania analysis of phosphorylated proteins, the pathways enriched in Network 2 are shown. The q value and the number of occurrences in the network versus the occurrences in the Network are shown.
Network Session 1 | Edge directed PKN network was manually curated from the literature. Annotated phosphorylation events, interactions, dissociations and cellular compartment are indicated.
Network Session 2 | Protein interaction (PPI) network was generated using a full phosphorylation dataset encompassing 321 proteins in Cytoscape using Genemania™ only considering protein-protein interactions with experimental evidence each node represents and individual protein and interactions are represented by edges. Functional analysis was carried out to identify key pathways enriched within the full network. Individual networks were generated from each function individually and are available as unique networks within the Network session.
Abbreviations
AKT, AKT serine/threonine kinase 1; ARTD1, ADP-ribosyltransferase diphtheria toxin-like 1; ATP, Adenosine triphosphate; BAF, BRG1-associated factor 53A; BCA, Bicinchoninic acid; CAMKK1, Calcium/calmodulin-dependent protein kinase kinase 1; CDK1/2, Cyclin dependent kinase; EMT, Epithelial to mesenchyme transition; ER, Oestrogen/Estrogen receptor; ERBB, Receptor tyrosine-protein kinase erbB-2; ERK1, Extracellular signal-regulated kinase 1; KH, K homology; MNK1, MAP kinase signal-integrated kinase 1; MSK, Mitogen and stress activated protein kinase 1; NUDT5, Nucleoside diphosphate-linked moiety X motif 5; NURF, Nucleosome-remodelling factor subunit; PAR, Poly-ADP-ribose; PARG, Poly-ADP-ribose glycohydrolase; PARP1, Poly-ADP-ribose polymerase; PCA, Principal component analysis; PGR, Progesterone Receptor; PKN, Prior knowledge network; PML, Prommyelocytic leukemia; PPI, Protein protein interaction; PRE, Progesterone responsive element.
References
1. Pina B, Bruggemeier U, Beato M. Nucleosome positioning modulates accessibility of regulatory proteins to the mouse mammary tumor virus promoter. Cell (1990) 60:719–31. doi: 10.1016/0092-8674(90)90087-U
2. Beato M, Herrlich P, Schutz G. Steroid hormone receptors: many actors in search of a plot. Cell (1995) 83:851–7. doi: 10.1016/0092-8674(95)90201-5
3. Migliaccio A, Piccolo D, Castoria G, Di Domenico M, Bilancio A, Lombardi M, et al. Activation of the Src/p21ras/Erk pathway by progesterone receptor via cross-talk with estrogen receptor. EMBO J (1998) 17(7):2008–18. doi: 10.1093/emboj/17.7.2008
4. Pedram A, Razandi M, Sainson RC, Kim JK, Hughes CC, Levin ER. A conserved mechanism for steroid receptor translocation to the plasma membrane. J Biol Chem (2007) 282:22278–88. doi: 10.1074/jbc.M611877200
5. Ballare C, Uhrig M, Bechtold T, Sancho E, Di Domenico M, Migliaccio A, et al. Two domains of the progesterone receptor interact with the estrogen receptor and are required for progesterone activation of the c-Src/Erk pathway in mammalian cells. Mol Cell Biol (2003) 23:1994–2008. doi: 10.1128/MCB.23.6.1994-2008.2003
6. Boonyaratanakornkit V, Scott MP, Ribon V, Sherman L, Anderson SM, Maller JL, et al. Progesterone receptor contains a proline-rich motif that directly interacts with SH3 domains and activates c-src family tyrosine kinases. Mol Cell (2001) 8:269–80. doi: 10.1016/S1097-2765(01)00304-5
7. Wright RH, Castellano G, Bonet J, Le Dily F, Font-Mateu J, Ballare C, et al. CDK2-dependent activation of PARP-1 is required for hormonal gene regulation in breast cancer cells. Genes Dev (2012) 26:1972–83. doi: 10.1101/gad.193193.112
8. Vicent GP, Ballare C, Nacht AS, Clausell J, Subtil-Rodriguez A, Quiles I, et al. Induction of progesterone target genes requires activation of erk and msk kinases and phosphorylation of histone H3. Mol Cell (2006) 24:367–81. doi: 10.1016/j.molcel.2006.10.011
9. Vicent GP, Nacht AS, Font-Mateu J, Castellano G, Gaveglia L, Ballare C, et al. Four enzymes cooperate to displace histone H1 during the first minute of hormonal gene activation. Genes Dev (2011) 25:845–62. doi: 10.1101/gad.621811
10. Vicent GP, Zaurin R, Ballare C, Nacht AS, Beato M. Erk signaling and chromatin remodeling in MMTV promoter activation by progestins. Nucl Recept Signal (2009) 7:e008. doi: 10.1621/nrs.07008
11. Vicent GP, Zaurin R, Nacht AS, Font-Mateu J, Le Dily F, Beato M. Nuclear factor 1 synergizes with progesterone receptor on the mouse mammary tumor virus promoter wrapped around a histone H3/H4 tetramer by facilitating access to the central hormone-responsive elements. J Biol Chem (2010) 285:2622–31. doi: 10.1074/jbc.M109.060848
12. Fu XD, Goglia L, Sanchez AM, Flamini M, Giretti MS, Tosi V, et al. Progesterone receptor enhances breast cancer cell motility and invasion via extranuclear activation of focal adhesion kinase. Endocr Relat Cancer (2010) 17:431–43. doi: 10.1677/ERC-09-0258
13. Garg D, Ng SSM, Baig KM, Driggers P, Segars J. Progesterone-mediated non-classical signaling. Trends Endocrinol Metab (2017) 28:656–68. doi: 10.1016/j.tem.2017.05.006
14. Takahashi A, Kato K, Kuboyama A, Inoue T, Tanaka Y, Kuhara A, et al. Induction of senescence by progesterone receptor-b activation in response to cAMP in ovarian cancer cells. Gynecol Oncol (2009) 113:270–6. doi: 10.1016/j.ygyno.2008.12.032
15. Rider V, Isuzugawa K, Twarog M, Jones S, Cameron B, Imakawa K, et al. Progesterone initiates wnt-beta-catenin signaling but estradiol is required for nuclear activation and synchronous proliferation of rat uterine stromal cells. J Endocrinol (2006) 191:537–48. doi: 10.1677/joe.1.07030
16. Hagan CR, Knutson TP, Lange CA. A common docking domain in progesterone receptor-b links dusp6 and ck2 signaling to proliferative transcriptional programs in breast cancer cells. Nucleic Acids Res (2013) 41:8926–42. doi: 10.1093/nar/gkt706
17. Lee II, Kim JJ. Influence of AKT on progesterone action in endometrial diseases. Biol Reprod (2014) 91:63. doi: 10.1095/biolreprod.114.119255
18. Wang Y, Feng H, Bi C, Zhu L, Pollard JW, Chen B. GSK-3beta mediates in the progesterone inhibition of estrogen induced cyclin D2 nuclear localization and cell proliferation in cyclin D1-/- mouse uterine epithelium. FEBS Lett (2007) 581:3069–75. doi: 10.1016/j.febslet.2007.05.072
19. Nacht AS, Pohl A, Zaurin R, Soronellas D, Quilez J, Sharma P, et al. Hormone-induced repression of genes requires BRG1-mediated H1.2 deposition at target promoters. EMBO J (2016) 35:1822–43. doi: 10.15252/embj.201593260
20. Vicent GP, Wright RH, Beato M. Linker histones in hormonal gene regulation. Biochim Biophys Acta (2016) 1859:520–5. doi: 10.1016/j.bbagrm.2015.10.016
21. Wright RH, Lioutas A, Le Dily F, Soronellas D, Pohl A, Bonet J, et al. ADP-ribose-derived nuclear ATP synthesis by NUDIX5 is required for chromatin remodeling. Science (2016) 352:1221–5. doi: 10.1126/science.aad9335
22. Le Dily F, Vidal E, Cuartero Y, Quilez J, Nacht AS, Vicent GP, et al. Hormone-control regions mediate steroid receptor-dependent genome organization. Genome Res (2019) 29:29–39. doi: 10.1101/gr.243824.118
23. Wright RHG, Le Dily F, Beato M. ATP, Mg(2+), nuclear phase separation, and genome accessibility. Trends Biochem Sci (2019) 44:565–74. doi: 10.1016/j.tibs.2019.03.001
24. Pickup KE, Pardow F, Carbonell-Caballero J, Lioutas A, Villanueva-Canas JL, Wright RHG, et al. Expression of oncogenic drivers in 3d cell culture depends on nuclear atp synthesis by nudt5. Cancers (Basel) (2019) 11(9):1337. doi: 10.3390/cancers11091337
25. Boonyaratanakornkit V, Bi Y, Rudd M, Edwards DP. The role and mechanism of progesterone receptor activation of extra-nuclear signaling pathways in regulating gene transcription and cell cycle progression. Steroids (2008) 73:922–8. doi: 10.1016/j.steroids.2008.01.010
26. Faivre E, Skildum A, Pierson-Mullany L, Lange CA. Integration of progesterone receptor mediated rapid signaling and nuclear actions in breast cancer cell models: role of mitogen-activated protein kinases and cell cycle regulators. Steroids (2005) 70:418–26. doi: 10.1016/j.steroids.2005.02.012
27. Qiu M, Olsen A, Faivre E, Horwitz KB, Lange CA. Mitogen-activated protein kinase regulates nuclear association of human progesterone receptors. Mol Endocrinol (2003) 17:628–42. doi: 10.1210/me.2002-0378
28. Skildum A, Faivre E, Lange CA. Progesterone receptors induce cell cycle progression via activation of mitogen-activated protein kinases. Mol Endocrinol (2005) 19:327–39. doi: 10.1210/me.2004-0306
29. Cheadle C, Cho-Chung YS, Becker KG, Vawter MP. Application of z-score transformation to affymetrix data. Appl Bioinf (2003) 2:209–17.
30. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res (2010) 38:W214–220. doi: 10.1093/nar/gkq537
31. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13:2498–504. doi: 10.1101/gr.1239303
32. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics (2011) 27:1739–40. doi: 10.1093/bioinformatics/btr260
33. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
34. Lanczky A, Gyorffy B. Web-based survival analysis tool tailored for medical research (KMplot): Development and implementation. J Med Internet Res (2021) 23(7):e27633. doi: 10.2196/27633
35. Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. tissue-based map of the human proteome. Science (2015) 347:1260419. doi: 10.1126/science.1260419
36. Haverinen M, Passinen S, Syvala H, Pasanen S, Manninen T, Tuohimaa P, et al. Heat shock protein 90 and the nuclear transport of progesterone receptor. Cell Stress Chaperones (2001) 6:256–62. doi: 10.1379/1466-1268(2001)006<0256:HSPATN>2.0.CO;2
37. Smith DF. Dynamics of heat shock protein 90-progesterone receptor binding and the disactivation loop model for steroid receptor complexes. Mol Endocrinol (1993) 7:1418–29. doi: 10.1210/mend.7.11.7906860
38. Kato S, Endoh H, Masuhiro Y, Kitamoto T, Uchiyama S, Sasaki H, et al. Activation of the estrogen receptor through phosphorylation by mitogen-activated protein kinase. Science (1995) 270:1491–4. doi: 10.1126/science.270.5241.1491
39. Reyes D, Ballaré C, Castellano G, Soronellas D, Bagó JR. Activation of mitogen- and stress-activated kinase 1 is required for proliferation of breast cancer cells in response to estrogens or progestins. Oncogene (2014) 33(12):1570–80. doi: 10.1038/onc.2013.95
40. Wierer M, Verde G, Pisano P, Molina H, Font-Mateu J, Di Croce L, et al. PLK1 signaling in breast cancer cells cooperates with estrogen receptor-dependent gene transcription. Cell Rep (2013) 3:2021–32. doi: 10.1016/j.celrep.2013.05.024
41. Boonyaratanakornkit V, McGowan E, Sherman L, Mancini MA, Cheskis BJ, Edwards DP. The role of extranuclear signaling actions of progesterone receptor in mediating progesterone regulation of gene expression and the cell cycle. Mol Endocrinol (2007) 21:359–75. doi: 10.1210/me.2006-0337
42. Ciruelos Gil EM. Targeting the PI3K/AKT/mTOR pathway in estrogen receptor-positive breast cancer. Cancer Treat Rev (2014) 40:862–71. doi: 10.1016/j.ctrv.2014.03.004
43. Ortega MA, Fraile-Martinez O, Asunsolo A, Bujan J, Garcia-Honduvilla N, Coca S. Signal transduction pathways in breast cancer: The important role of PI3K/Akt/mTOR. J Oncol (2020) 2020:9258396. doi: 10.1155/2020/9258396
44. Lange CA. Integration of progesterone receptor action with rapid signaling events in breast cancer models. J Steroid Biochem Mol Biol (2008) 108:203–12. doi: 10.1016/j.jsbmb.2007.09.019
45. Li X, Wong J, Tsai SY, Tsai MJ, O’Malley BW. Progesterone and glucocorticoid receptors recruit distinct coactivator complexes and promote distinct patterns of local chromatin modification. Mol Cell Biol (2003) 23:3763–73. doi: 10.1128/MCB.23.11.3763-3773.2003
46. Malbeteau L, Poulard C, Languilaire C, Mikaelian I, Flamant F, Le Romancer M, et al. PRMT1 is critical for the transcriptional activity and the stability of the progesterone receptor. iScience (2020) 23:101236. doi: 10.1016/j.isci.2020.101236
47. Dolniak B, Katsoulidis E, Carayol N, Altman JK, Redig AJ, Tallman MS, et al. Regulation of arsenic trioxide-induced cellular responses by Mnk1 and Mnk2. J Biol Chem (2008) 283:12034–42. doi: 10.1074/jbc.M708816200
48. Guo Q, Li VZ, Nichol JN, Huang F, Yang W, Preston SEJ, et al. MNK1/NODAL signaling promotes invasive progression of breast ductal carcinoma in situ. Cancer Res (2019) 79:1646–57. doi: 10.1158/0008-5472.CAN-18-1602
49. Wheater MJ, Johnson PW, Blaydes JP. The role of MNK proteins and eIF4E phosphorylation in breast cancer cell proliferation and survival. Cancer Biol Ther (2010) 10:728–35. doi: 10.4161/cbt.10.7.12965
50. Grimmler M, Wang Y, Mund T, Cilensek Z, Keidel EM, Waddell MB, et al. Cdk-inhibitory activity and stability of p27Kip1 are directly regulated by oncogenic tyrosine kinases. Cell (2007) 128:269–80. doi: 10.1016/j.cell.2006.11.047
51. Trevino LS, Bolt MJ, Grimm SL, Edwards DP, Mancini MA, Weigel NL. Differential regulation of progesterone receptor-mediated transcription by CDK2 and DNA-PK. Mol Endocrinol (2016) 30:158–72. doi: 10.1210/me.2015-1144
52. Lange CA, Shen T, Horwitz KB. Phosphorylation of human progesterone receptors at serine-294 by mitogen-activated protein kinase signals their degradation by the 26S proteasome. Proc Natl Acad Sci U.S.A. (2000) 97:1032–7. doi: 10.1073/pnas.97.3.1032
53. Knutson TP, Truong TH, Ma S, Brady NJ, Sullivan ME, Raj G, et al. Posttranslationally modified progesterone receptors direct ligand-specific expression of breast cancer stem cell-associated gene programs. J Hematol Oncol (2017) 10:89. doi: 10.1186/s13045-017-0462-7
54. Knotts TA, Orkiszewski RS, Cook RG, Edwards DP, Weigel NL. Identification of a phosphorylation site in the hinge region of the human progesterone receptor and additional amino-terminal phosphorylation sites. J Biol Chem (2001) 276:8475–83. doi: 10.1074/jbc.M009805200
55. Price TM, Dai Q. The role of a mitochondrial progesterone receptor (PR-m) in progesterone action. Semin Reprod Med (2015) 33:185–94. doi: 10.1055/s-0035-1552583
56. Dai Q, Likes CE 3rd, Luz AL, Mao L, Yeh JS, Wei Z, et al. A mitochondrial progesterone receptor increases cardiac beta-oxidation and remodeling. J Endocr Soc (2019) 3:446–67. doi: 10.1210/js.2018-00219
57. Vicent GP, Zaurin R, Nacht AS, Li A, Font-Mateu J, Le Dily F, et al. Two chromatin remodeling activities cooperate during activation of hormone responsive promoters. PloS Genet (2009) 5:e1000567. doi: 10.1371/journal.pgen.1000567
58. Park J, Wood MA, Cole MD. BAF53 forms distinct nuclear complexes and functions as a critical c-myc-interacting nuclear cofactor for oncogenic transformation. Mol Cell Biol (2002) 22:1307–16. doi: 10.1128/MCB.22.5.1307-1316.2002
59. von Eyss B, Jaenicke LA, Kortlever RM, Royla N, Wiese KE, Letschert S, et al. A MYC-driven change in mitochondrial dynamics limits yap/taz function in mammary epithelial cells and breast cancer. Cancer Cell (2015) 28:743–57. doi: 10.1016/j.ccell.2015.10.013
60. Hurley RL, Anderson KA, Franzone JM, Kemp BE, Means AR, Witters LA. The Ca2+/calmodulin-dependent protein kinase kinases are AMP-activated protein kinase kinases. J Biol Chem (2005) 280:29060–6. doi: 10.1074/jbc.M503824200
61. Hadad SM, Fleming S, Thompson AM. Targeting AMPK: a new therapeutic opportunity in breast cancer. Crit Rev Oncol Hematol (2008) 67:1–7. doi: 10.1016/j.critrevonc.2008.01.007
62. Montero JC, Esparis-Ogando A, Re-Louhau MF, Seoane S, Abad M, Calero R, et al. Active kinase profiling, genetic and pharmacological data define mTOR as an important common target in triple-negative breast cancer. Oncogene (2014) 33:148–56. doi: 10.1038/onc.2012.572
63. Yamaguchi H, Taouk GM. A potential role of YAP/TAZ in the interplay between metastasis and metabolic alterations. Front Oncol (2020) 10:928. doi: 10.3389/fonc.2020.00928
64. Zhao H, Orhan YC, Zha X, Esencan E, Chatterton RT, Bulun SE. AMP-activated protein kinase and energy balance in breast cancer. Am J Transl Res (2017) 9:197–213.
65. Wu L, Huang XJ, Yang CH, Deng SS, Qian M, Zang Y, et al. 5’-AMP-activated protein kinase (AMPK) regulates progesterone receptor transcriptional activity in breast cancer cells. Biochem Biophys Res Commun (2011) 416:172–7. doi: 10.1016/j.bbrc.2011.11.018
66. McInnes KJ, Brown KA, Hunger NI, Simpson ER. Regulation of LKB1 expression by sex hormones in adipocytes. Int J Obes (Lond) (2012) 36:982–5. doi: 10.1038/ijo.2011.172
67. McInnes KJ, Corbould A, Simpson ER, Jones ME. Regulation of adenosine 5’,monophosphate-activated protein kinase and lipogenesis by androgens contributes to visceral obesity in an estrogen-deficient state. Endocrinology (2006) 147:5907–13. doi: 10.1210/en.2006-0879
68. Sanchez-Margalet V, Najib S. Sam68 is a docking protein linking GAP and PI3K in insulin receptor signaling. Mol Cell Endocrinol (2001) 183:113–21. doi: 10.1016/S0303-7207(01)00587-1
69. Mamay CL, Mingo-Sion AM, Wolf DM, Molina MD, Van Den Berg CL. An inhibitory function for JNK in the regulation of IGF-I signaling in breast cancer. Oncogene (2003) 22:602–14. doi: 10.1038/sj.onc.1206186
70. Eto I. Upstream molecular signaling pathways of p27(Kip1) expression: effects of 4-hydroxytamoxifen, dexamethasone, and retinoic acids. Cancer Cell Int (2010) 10:3. doi: 10.1186/1475-2867-10-3
71. Tzatsos A. Raptor binds the SAIN (Shc and IRS-1 NPXY binding) domain of insulin receptor substrate-1 (IRS-1) and regulates the phosphorylation of IRS-1 at ser-636/639 by mTOR. J Biol Chem (2009) 284:22525–34. doi: 10.1074/jbc.M109.027748
72. Hiratani K, Haruta T, Tani A, Kawahara J, Usui I, Kobayashi M. Roles of mTOR and JNK in serine phosphorylation, translocation, and degradation of IRS-1. Biochem Biophys Res Commun (2005) 335:836–42. doi: 10.1016/j.bbrc.2005.07.152
73. Xu B, Bird VG, Miller WT. Substrate specificities of the insulin and insulin-like growth factor 1 receptor tyrosine kinase catalytic domains. J Biol Chem (1995) 270:29825–30. doi: 10.1074/jbc.270.50.29825
74. Valverde AM, Arribas M, Mur C, Navarro P, Pons S, Cassard-Doulcier AM, et al. Insulin-induced up-regulated uncoupling protein-1 expression is mediated by insulin receptor substrate 1 through the phosphatidylinositol 3-kinase/Akt signaling pathway in fetal brown adipocytes. J Biol Chem (2003) 278:10221–31. doi: 10.1074/jbc.M209363200
75. Paul A, Edwards J, Pepper C, Mackay S. Inhibitory-kappaB kinase (IKK) alpha and nuclear factor-kappaB (NFkappaB)-inducing kinase (NIK) as anti-cancer drug targets. Cells (2018) 7(10):176. doi: 10.3390/cells7100176
76. Schramek D, Leibbrandt A, Sigl V, Kenner L, Pospisilik JA, Lee HJ, et al. Osteoclast differentiation factor RANKL controls development of progestin-driven mammary cancer. Nature (2010) 468:98–102. doi: 10.1038/nature09387
77. Shostak K, Chariot A. NF-kappaB, stem cells and breast cancer: the links get stronger. Breast Cancer Res (2011) 13:214. doi: 10.1186/bcr2886
78. Gonzalez-Suarez E, Jacob AP, Jones J, Miller R, Roudier-Meyer MP, Erwert R, et al. RANK ligand mediates progestin-induced mammary epithelial proliferation and carcinogenesis. Nature (2010) 468:103–7. doi: 10.1038/nature09495
79. Cabodi S, Moro L, Baj G, Smeriglio M, Di Stefano P, Gippone S, et al. p130Cas interacts with estrogen receptor alpha and modulates non-genomic estrogen signaling in breast cancer cells. J Cell Sci (2004) 117:1603–11. doi: 10.1242/jcs.01025
80. Cabodi S, Tinnirello A, Di Stefano P, Bisaro B, Ambrosino E, Castellano I, et al. p130Cas as a new regulator of mammary epithelial cell proliferation, survival, and HER2-neu oncogene-dependent breast tumorigenesis. Cancer Res (2006) 66:4672–80. doi: 10.1158/0008-5472.CAN-05-2909
81. Bonfiglio JJ, Colby T, Matic I. Mass spectrometry for serine ADP-ribosylation? think o-glycosylation! Nucleic Acids Res (2017) 45:6259–64. doi: 10.1093/nar/gkx446
82. Bonfiglio JJ, Fontana P, Zhang Q, Colby T, Gibbs-Seymour I, Atanassov I, et al. Serine ADP-ribosylation depends on HPF1. Mol Cell (2017) 65:932–940.e936. doi: 10.1016/j.molcel.2017.01.003
83. Leidecker O, Bonfiglio JJ, Colby T, Zhang Q, Atanassov I, Zaja R, et al. Serine is a new target residue for endogenous ADP-ribosylation on histones. Nat Chem Biol (2016) 12:998–1000. doi: 10.1038/nchembio.2180
84. Liu Q, Florea BI, Filippov DV. ADP-ribosylation goes normal: serine as the major site of the modification. Cell Chem Biol (2017) 24:431–2. doi: 10.1016/j.chembiol.2017.04.003
85. Martello R, Leutert M, Jungmichel S, Bilan V, Larsen SC, Young C, et al. Proteome-wide identification of the endogenous ADP-ribosylome of mammalian cells and tissue. Nat Commun (2016) 7:12917. doi: 10.1038/ncomms12917
86. Messner S, Altmeyer M, Zhao H, Pozivil A, Roschitzki B, Gehrig P, et al. PARP1 ADP-ribosylates lysine residues of the core histone tails. Nucleic Acids Res (2010) 38:6350–62. doi: 10.1093/nar/gkq463
87. Hassa PO, Haenni SS, Elser M, Hottiger MO. Nuclear ADP-ribosylation reactions in mammalian cells: where are we today and where are we going? Microbiol Mol Biol Rev (2006) 70:789–829. doi: 10.1128/MMBR.00040-05
88. Hoch NC, Polo LM. ADP-ribosylation: from molecular mechanisms to human disease. Genet Mol Biol (2019) 43:e20190075. doi: 10.1590/1678-4685-GMB-2019-0075
89. Leung AK. Poly(ADP-ribose): an organizer of cellular architecture. J Cell Biol (2014) 205:613–9. doi: 10.1083/jcb.201402114
90. Thomas C, Tulin AV. Poly-ADP-ribose polymerase: machinery for nuclear processes. Mol Aspects Med (2013) 34:1124–37. doi: 10.1016/j.mam.2013.04.001
91. McKay MM, Ritt DA, Morrison DK. Signaling dynamics of the KSR1 scaffold complex. Proc Natl Acad Sci U.S.A. (2009) 106:11022–7. doi: 10.1073/pnas.0901590106
92. Holaska JM, Wilson KL. An emerin “proteome”: purification of distinct emerin-containing complexes from HeLa cells suggests molecular basis for diverse roles including gene regulation, mRNA splicing, signaling, mechanosensing, and nuclear architecture. Biochemistry (2007) 46:8897–908. doi: 10.1021/bi602636m
93. Ranade D, Pradhan R, Jayakrishnan M, Hegde S, Sengupta K. Lamin A/C and emerin depletion impacts chromatin organization and dynamics in the interphase nucleus. BMC Mol Cell Biol (2019) 20:11. doi: 10.1186/s12860-019-0192-5
94. Nelson MB, Nyhus JK, Oravecz-Wilson KI, Barbera-Guillem E. Tumor cells express FcgammaRI which contributes to tumor cell growth and a metastatic phenotype. Neoplasia (2001) 3:115–24. doi: 10.1038/sj.neo.7900140
95. Elizalde PV, Cordo Russo RI, Chervo MF, Schillaci R. ErbB-2 nuclear function in breast cancer growth, metastasis and resistance to therapy. Endocr Relat Cancer (2016) 23:T243–57. doi: 10.1530/ERC-16-0360
96. Beguelin W, Diaz Flaque MC, Proietti CJ, Cayrol F, Rivas MA, Tkach M, et al. Progesterone receptor induces ErbB-2 nuclear translocation to promote breast cancer growth via a novel transcriptional effect: ErbB-2 function as a coactivator of Stat3. Mol Cell Biol (2010) 30:5456–72. doi: 10.1128/MCB.00012-10
97. Hsu JL, Hung MC. The role of HER2, EGFR, and other receptor tyrosine kinases in breast cancer. Cancer Metastasis Rev (2016) 35:575–88. doi: 10.1007/s10555-016-9649-6
98. Proietti CJ, Rosemblit C, Beguelin W, Rivas MA, Diaz Flaque MC, Charreau EH, et al. Activation of Stat3 by heregulin/ErbB-2 through the co-option of progesterone receptor signaling drives breast cancer growth. Mol Cell Biol (2009) 29:1249–65. doi: 10.1128/MCB.00853-08
99. Thomas RS, Sarwar N, Phoenix F, Coombes RC, Ali S. Phosphorylation at serines 104 and 106 by Erk1/2 MAPK is important for estrogen receptor-alpha activity. J Mol Endocrinol (2008) 43(4):173–84. doi: 10.1038/onc.2013.95
100. de Leeuw R, Neefjes J, Michalides R. A role for estrogen receptor phosphorylation in the resistance to tamoxifen. Int J Breast Cancer (2008) 2011:232435. doi: 10.4061/2011/232435
101. Jeffreys SA, Powter B, Balakrishnar B, Mok K, Soon P, Franken A, Neubauer H, et al. Endocrine resistance in breast cancer: The role of estrogen receptor stability. Cells (2011) 9(9):2077. doi: 10.3390/cells9092077
102. Skliris GP, Nugent ZJ, Rowan BG, Penner CR, Watson PH, Murphy LC. A phosphorylation code for oestrogen receptor-alpha predicts clinical outcome to endocrine therapy in breast cancer. Endocr Relat Cancer (2011) 17(3):589–97. doi: 10.1677/ERC-10-0030
103. Alayev A, Salamon RS, Berger SM, Schwartz NS, Cuesta R, Snyder RB, et al. mTORC1 directly phosphorylates and activates ERα upon estrogen stimulation. Oncogene (2016) 7(35):3535–4397. doi: 10.1038/onc.2015.414
104. Cui X, Lazard Z, Zhang P, Hopp TA, Lee AV. Progesterone crosstalks with insulin-like growth factor signaling in breast cancer cells via induction of insulin receptor substrate-2. Oncogene (2003) 22:6937–41. doi: 10.1038/sj.onc.1206803
105. Cui X, Zhang P, Deng W, Oesterreich S, Lu Y, Mills GB, et al. Insulin-like growth factor-I inhibits progesterone receptor expression in breast cancer cells via the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin pathway: progesterone receptor as a potential indicator of growth factor activity in breast cancer. Mol Endocrinol (2003) 17:575–88. doi: 10.1210/me.2002-0318
106. Ibrahim YH, Byron SA, Cui X, Lee AV, Yee D. Progesterone receptor-b regulation of insulin-like growth factor-stimulated cell migration in breast cancer cells via insulin receptor substrate-2. Mol Cancer Res (2008) 6:1491–8. doi: 10.1158/1541-7786.MCR-07-2173
107. Kramer EA, Seeger H, Kramer B, Wallwiener D, Mueck AO. Characterization of the stimulatory effect of medroxyprogesterone acetate and chlormadinone acetate on growth factor treated normal human breast epithelial cells. J Steroid Biochem Mol Biol (2006) 98:174–8. doi: 10.1016/j.jsbmb.2005.11.002
108. Behera MA, Dai Q, Garde R, Saner C, Jungheim E, Price TM. Progesterone stimulates mitochondrial activity with subsequent inhibition of apoptosis in MCF-10A benign breast epithelial cells. Am J Physiol Endocrinol Metab (2009) 297:E1089–1096. doi: 10.1152/ajpendo.00209.2009
109. Demonacos CV, Karayanni N, Hatzoglou E, Tsiriyiotis C, Spandidos DA, Sekeris CE. Mitochondrial genes as sites of primary action of steroid hormones. Steroids (1996) 61:226–32. doi: 10.1016/0039-128X(96)00019-0
110. Hatzoglou E, Sekeris CE. The detection of nucleotide sequences with strong similarity to hormone responsive elements in the genome of eubacteria and archaebacteria and their possible relation to similar sequences present in the mitochondrial genome. J Theor Biol (1997) 184:339–44. doi: 10.1006/jtbi.1996.0285
111. Huang EJ, Reichardt LF. Trk receptors: roles in neuronal signal transduction. Annu Rev Biochem (2003) 72:609–42. doi: 10.1146/annurev.biochem.72.121801.161629
112. Reichardt LF. Neurotrophin-regulated signalling pathways. Philos Trans R Soc Lond B Biol Sci (2006) 361:1545–64. doi: 10.1098/rstb.2006.1894
113. Jin W. Roles of trkc signaling in the regulation of tumorigenicity and metastasis of cancer. Cancers (Basel) (2020) 12(1):147. doi: 10.3390/cancers12010147
114. Meng L, Liu B, Ji R, Jiang X, Yan X, Xin Y. Targeting the BDNF/TrkB pathway for the treatment of tumors. Oncol Lett (2019) 17:2031–9. doi: 10.3892/ol.2018.9854
115. Regua AT, Doheny D, Arrigo A, Lo HW. Trk receptor tyrosine kinases in metastasis and cancer therapy. Discovery Med (2019) 28:195–203.
116. Wang R, Fu L, Li J, Zhao D, Zhao Y, Yin L. Microarray analysis for differentially expressed genes between stromal and epithelial cells in development and metastasis of invasive breast cancer. J Comput Biol (2020) 27(12):1631–43. doi: 10.1089/cmb.2019.0154
117. Dobson M, Ramakrishnan G, Ma S, Kaplun L, Balan V, Fridman R, et al. Bimodal regulation of FoxO3 by AKT and 14-3-3. Biochim Biophys Acta (2011) 1813:1453–64. doi: 10.1016/j.bbamcr.2011.05.001
118. Pennington KL, Chan TY, Torres MP, Andersen JL. The dynamic and stress-adaptive signaling hub of 14-3-3: emerging mechanisms of regulation and context-dependent protein-protein interactions. Oncogene (2018) 37:5587–604. doi: 10.1038/s41388-018-0348-3
119. Tzivion G, Dobson M, Ramakrishnan G. FoxO transcription factors; regulation by AKT and 14-3-3 proteins. Biochim Biophys Acta (2011) 1813:1938–45. doi: 10.1016/j.bbamcr.2011.06.002
120. Tzivion G, Hay N. PI3K-AKT-FoxO axis in cancer and aging. Biochim Biophys Acta (2011) 1813:1925. doi: 10.1016/j.bbamcr.2011.08.014
Keywords: progesterone, breast cancer, phosphoproteome, signalling, MAPK/ERK signalling, chromatin, PARylation, cell proliferation
Citation: Wright RHG, Vastolo V, Oliete JQ, Carbonell-Caballero J and Beato M (2022) Global signalling network analysis of luminal T47D breast cancer cells in response to progesterone. Front. Endocrinol. 13:888802. doi: 10.3389/fendo.2022.888802
Received: 03 March 2022; Accepted: 11 July 2022;
Published: 11 August 2022.
Edited by:
Norbert Nass, Städtische Klinikum Dessau, GermanyReviewed by:
Giovanni Nassa, University of Salerno, ItalyEdwin Lasonder, Northumbria University, United Kingdom
Copyright © 2022 Wright, Vastolo, Oliete, Carbonell-Caballero and Beato. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roni H. G. Wright, rhgwright@uic.es; Miguel Beato, miguel.beato@crg.eu