Phosphoproteome and drug-response effects mediated by the three protein phosphatase 2A inhibitor proteins CIP2A, SET, and PME-1

Protein phosphatase 2A (PP2A) critically regulates cell signaling and is a human tumor suppressor. PP2A complexes are modulated by proteins such as cancerous inhibitor of protein phosphatase 2A (CIP2A), protein phosphatase methylesterase 1 (PME-1), and SET nuclear proto-oncogene (SET) that often are deregulated in cancers. However, how they impact cellular phosphorylation and how redundant they are in cellular regulation is poorly understood. Here, we conducted a systematic phosphoproteomics screen for phosphotargets modulated by siRNA-mediated depletion of CIP2A, PME-1, and SET (to reactivate PP2A) or the scaffolding A-subunit of PP2A (PPP2R1A) (to inhibit PP2A) in HeLa cells. We identified PP2A-modulated targets in diverse cellular pathways, including kinase signaling, cytoskeleton, RNA splicing, DNA repair, and nuclear lamina. The results indicate nonredundancy among CIP2A, PME-1, and SET in phosphotarget regulation. Notably, PP2A inhibition or reactivation affected largely distinct phosphopeptides, introducing a concept of nonoverlapping phosphatase inhibition- and activation-responsive sites (PIRS and PARS, respectively). This phenomenon is explained by the PPP2R1A inhibition impacting primarily dephosphorylated threonines, whereas PP2A reactivation results in dephosphorylation of clustered and acidophilic sites. Using comprehensive drug-sensitivity screening in PP2A-modulated cells to evaluate the functional impact of PP2A across diverse cellular pathways targeted by these drugs, we found that consistent with global phosphoproteome effects, PP2A modulations broadly affect responses to more than 200 drugs inhibiting a broad spectrum of cancer-relevant targets. These findings advance our understanding of the phosphoproteins, pharmacological responses, and cellular processes regulated by PP2A modulation and may enable the development of combination therapies.

motes in vivo tumorigenesis (9 -13), and recent studies implicate PP2A as a critical determinant to kinase inhibitor responses (8, 14 -16). Furthermore, orally bioavailable and nontoxic PP2A-reactivating small molecule compounds show robust antitumor effects both as monotherapies and in combination with kinase inhibitors (16 -19). However, despite its critical role in cancer, and in many other diseases (20), the effects of endogenous PP2A-modulating mechanisms on phospho-signaling have not yet been systematically analyzed.
PP2A activity is affected by diverse mechanisms, including mutations (21), and interactions with endogenous proteins that modulate PP2A function. Represented in this study by CIP2A, PME-1, and SET (Fig. 1A), these PP2A inhibitor proteins (PAIPs hereafter) (4) have been found to suppress PP2A function toward its phosphosubstrates (3,4,22,23). Despite their shared function as PAIPs, CIP2A, PME-1, and SET do not share structural features (24 -26), and the mechanisms by which they inhibit PP2A activity toward their selected phosphoprotein targets are vastly different (23, 24, 26 -28). Each of them may also impact PP2A activity by different mechanisms, as exemplified by PME-1 that impacts the PP2A C-subunit both by its methylesterase activity and by direct binding to the catalytic center of the PP2AC (24). However, despite their differential mode of PP2A regulation, the cancer-relevant phenotypes resulting from CIP2A, PME-1, or SET modulation can all be rescued by concomitant PP2A inhibition (5,11,29,30). Furthermore, their overexpression can substitute for viral small-t as a PP2A inhibitor in human cell transformation models (10,31,32). Inhibition of these three PAIPs suppresses malignant cell growth and tumorigenesis in vivo (5,10,11,30). Moreover, PP2A inhibition by CIP2A overexpression in normal mouse brains induced the symptoms and pathogenic changes resembling Alzheimer's disease (20). Based on these data, modulation of PP2A through interactions with CIP2A, PME-1, or SET has recently emerged as an attractive novel therapeutic approach for cancer and Alzheimer's disease (4,20,22,23,33,34). However, the processes regulated by and the functional redundancy of these PP2Amodulating proteins remain poorly understood.
To address these unresolved questions, we characterized PP2A modulated phosphotargets in HeLa cells depleted of CIP2A, PME-1, and SET as well as the scaffolding A-subunit of PP2A (PPP2R1A). We identified targets in diverse cellular processes and functions, and we showed significant differences between PAIPs in their functional outcomes. To validate the global functional relevance of these modulations, we performed a drug sensitivity screen and found a broad impact on cancer cells across multiple drug target families. This resource provides insights into phosphatase-mediated signal transduction in human cells and helps advance development of rational mono-and combination therapies.

Identification and quantification of the PP2A-modulated targets via LC-MS/MS
PPP2R1A is the predominant PP2A scaffold A-subunit that is required for functional PP2A complex formation (Fig. 1A). Unlike depletion of the catalytic PP2A subunit PPP2CA, siRNA of PPP2R1A did not cause cell lethality and, on the other hand, did not affect the PPP2CA expression (Fig. S1A). However, and as expected from the literature (35), PPP2R1A depletion resulted in destabilization of PPP2R5A (B56␣) and PPP2R2A (B55␣) B-subunits (Fig. S1B). These results confirm that PPP2R1A depletion results in comprehensive PP2A complex inhibition. To map the phosphoproteome changes regulated by the three most prevalent PAIPs, we post-transcriptionally inhibited CIP2A, PME-1, and SET (Fig. 1A). As confirmed by Western blot analysis, none of the targeted PAIPs interfered with each other's expression at the protein level ( Fig. 1B) nor had notable effects on B-subunit expression (Fig. S1B). For MS analysis, each gene was targeted by three to four independent siRNA sequences and analyzed by MS in triplicate. Changes in peptide phosphorylation were analyzed 72 h after siRNA depletion via label-free phosphoproteomics pipeline (Fig. S1C) (36). The performance of the pipeline was previously validated as we observed high reproducibility between the global phosphoproteomics data and a phospho-specific antibody survey (36). We confirmed that global phosphorylation changes were largely replicated by different siRNAs targeting the same gene (Fig.  S1D). To focus on the most reliable and robust phosphosites, we statistically selected a subset of peptides that were consistently phosphoregulated across our replicates (limma q Ͻ0.05).
To minimize the likelihood that the phosphoproteome was mapped in cells that were affected by transfection procedures, we selected a late time point of 72 h post-transfection for our proteomics assays. To confirm that our observations were not confounded by mitotic arrest events (37,38), we determined mitotic index values. Importantly, more than 90% of cells in each condition remained nonmitotic at the same 72-h time point of the phosphoproteome analysis (Fig. S1E). Therefore, our data represent phosphoproteomes of stably PP2A-modulated interphase cells that are not cofounded by mitotic arrest. Of note, our dataset represents a single time point that does not naturally allow us to distinguish direct and indirect phosphotargets. Such single time point data are fully relevant, especially for understanding the cells in which PAIPs are differentially expressed, as is often the case between normal and cancerous tissues (4). Nevertheless, in addition to long-term siRNA treatments, we validated more immediate regulation of some of the significantly regulated phosphopeptides by pharmacological PP2A inhibition and activation. By using okadaic acid as a phosphatase inhibitor, and FTY-720 or DT-061 as PP2A activators, we were able to see consistent phosphoregulation of the selected targets as compared with siRNA data and after 12-24 h of treatment (Fig. S1, F and G).

Global phosphorylation changes by PP2A modulations
In support of the global relevance in cellular Ser/Thr phosphorylation regulation, 43% of the phosphopeptides were significantly regulated by either inhibition of PAIPs or PPP2R1A (q Ͻ0.05), and 57% of the phosphoproteins had at least one differentially-regulated phosphopeptide. The catalogue of differentially-regulated phosphopeptides (q Ͻ0.05) is included as Table S1. Depletion of the PPP2R1A resulted in a robust increase in the average phosphorylation level across the HeLa cell phosphoproteome (Fig. 1C and Fig. S1, D and H). In con-Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes trast to PPP2R1A, cells that were depleted of any of the PAIPs displayed global dephosphorylation ( Fig. 1C and Fig. S1, D and H). These results indicate that the dominant function of PAIPs is indeed the inhibition of dephosphorylation activity, thus validating the previous findings that the functional effects of these proteins can be rescued by concomitant PP2A inhibition (5,11,29,30) and that they can substitute for small-t as a PP2A inhibitor (10,31,32).
Among PAIPs, SET was the most robust inhibitor of protein dephosphorylation that regulated 30.5% of the phosphopeptides (Fig. S1H). In accord with a common regulatory mechanism with SET (i.e. PP2A inhibition), the majority of the CIP2A and PME-1 target peptides overlapped with the target peptides of SET (Fig. S1I). Nevertheless, CIP2A regulated 70 unique targets not significantly regulated by PME-1, and PME-1 regulated 115 targets that were not regulated by CIP2A (Fig. S1I). These nonredundant targets could be apparent due to the differential selectivity of CIP2A and PME-1 for different PP2A trimer complexes (23, 24, 26 -28). The larger number of SET targets compared with CIP2A and PME-1 (Fig. S1I) may be at least partly due to the direct inhibition of the catalytic activity of PP2Ac by SET (27,28) and the limited selectivity toward specific B-subunit-containing PP2A trimers, as is the case with PME-1 or CIP2A (10,26,39). Importantly, as only few phosphopeptides that were consistently and significantly regulated by PAIPs exhibited increased phosphorylation (0.5-2.7%) (Fig. S1H), these data indicate that PAIP inhibition cannot be compensated by increased kinase activity toward another set of phosphotargets, at least at the 72-h time point.

PP2A inhibition-and activation-responsive phosphotargets
Accumulating data show that specific phosphorylation sites tend to be predominantly phosphorylated, or dephosphorylated, depending on the amino acid context of the phosphorylated   Figure 1. PPP2R1A or PAIP inhibition regulates largely nonoverlapping phosphosites. A, schematic presentation of PP2A complex and its activity regulation by subunit composition and by endogenous inhibitor proteins (PAIPs). Whereas B-subunits regulate substrate specificity of the trimeric PP2A complex, each PAIP (CIP2A, PME-1, and SET) regulates PP2A complex activity by a distinct mechanism. B, Western blot analysis of selectivity of siRNAs for the studied PP2A proteins in HeLa cells 72 h after siRNA transfection. C, PPPP2R1A or PAIP inhibition causes global shift in HeLa cell phosphoproteome. Shown are fold change distributions of representative replicates (top) and mean log 2-fold changes of all replicates (bottom). Fig. S1H shows fractions of differentially regulated phosphosites for each condition. D, Venn diagram shows the overlap between PPP2R1A and PAIP-regulated peptides at q Ͻ0.05. Total number of significantly-regulated peptides are shown in parentheses. E, analysis of targets of different kinases by phosphomotif antibodies arranged from the most PPP2R1A-responsive phosphoproteins (PIRS) to the most SET-responsive phosphoproteins (PARS) in descending order. F, unsupervised soft clustering analysis of the high-confidence phosphoproteome data into six clusters. Cluster membership percentage is indicated by a color scale, and representative peptides are listed inside the plots. Cluster centers (indicated by black dots) exhibit mainly either up-regulation (cluster 2) or down-regulation (clusters 3-6). Cluster 1 instead shows the pattern of phosphoregulation for peptides that respond to both PP2A inhibition and activation. Cluster membership of each clustered phosphopeptide is shown in Table S1.

Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes
site (40 -44). Furthermore, previous global phosphoproteome analyses demonstrated that a majority of phosphosites had either high or low stoichiometry, although intermediate stoichiometry was rare (1,42). Therefore, we hypothesized that PP2A inhibition would primarily affect those sites that are not constitutively occupied by a phosphate (phosphatase inhibition-responsive site (PIRS)), whereas inhibition of PAIPs might preferentially regulate high-stoichiometry phosphosites (phosphatase activation-responsive site (PARS)).
In support of this hypothesis, the majority of the same peptides that were regulated by PPP2R1A inhibition were not regulated by depletion of PAIPs and vice versa (211/320 PPP2R1A targets (66%) and 568/677 PAIP targets (84%)) (Fig. 1D). These findings were methodologically and independently validated by phospho-motif antibodies showing that PPP2R1A depletion increased the phosphorylation of a different set of phosphomotifs compared with those that were dephosphorylated by increased PP2A activity in SET-depleted cells (Fig. 1E).
To analyze the characteristics of PIRS and PARS phosphosites, we performed an unsupervised soft clustering analysis (45). Consistent with the observation of distinct targets of PPP2R1A and PAIPs (Fig. 1D), the peptides in five of the six unsupervised clusters were regulated exclusively toward either increased phosphorylation (by PPP2R1A inhibition in cluster 2) or dephosphorylation (by inhibition of PAIPs in clusters 3-6) (Fig. 1F). As an exception, the peptides in cluster 1, consisting of 8.5% of all clustered phosphopeptides, followed the pattern of being responsive to both PP2A inhibition and activation (Fig.  1F). Furthermore, we observed that the target peptides in clusters 5 and 6 were most potently impacted by SET, followed by PME-1, and then CIP2A (Fig. 1F), supporting the partially dif-ferential outcomes via inhibition of each PAIP (Fig. S1, H and I).
Notably, the clustering analysis remained robust when normalizing the phosphopeptides in each cluster relative to total protein abundance of the nonenriched sample for proteins in which we could identify Ն2 unenriched peptides (Fig. S2). Thus, the phosphosite clustering represents true phosphorylation stoichiometry changes rather than global changes in protein abundance. Finally, to ensure that these findings were not limited to the current model of HeLa cells, we performed soft clustering of PP2A-modulated phosphoproteome from A549 lung cancer cell line from a previous study (16). Similar to HeLa cells, the clusters were characterized by either up-regulation or down-regulation but not both in the same cluster. Consistently, most prominent up-regulations were observed in PPP2R1Adepleted cells and down-regulations were observed in PAIPdepleted cells (Fig. S3).

Phosphorylation stoichiometry and amino acid context of PIRS and PARS phosphopeptides
Next, we further assessed phosphorylation stoichiometry and amino acid context of PP2A-modulated phosphopeptides in each cluster. To determine baseline phosphorylation stoichiometry, we assessed cells with no modulation of PP2A activity (i.e. cells transfected with scrambled siRNAs). Supportive of our conclusions, peptides belonging to clusters 1 and 2 (Fig. 1F) had a significantly lower baseline phosphopeptide/total protein abundance than peptides in clusters 3-6 ( Fig. 2A). The broad classification of most of the PP2A-modulated peptides as PIRS or PARS (except for those in cluster 1) was also supported by analysis of peptide sequences for determinants known to impact phosphate occupancy (40 -44). Previous studies indi- Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes cate that a higher PP2A/kinase activity ratio is associated with threonines rather than serines (40,43,44), whereas acidophilic kinase target sites exhibit a higher phosphate occupancy on average (42). To analyze our data according to these parameters, clusters 1 and 2 (cluster 1 & 2 hereafter) and clusters 5 and 6 (cluster 5 & 6 hereafter) were combined to yield four nearly equally-sized clusters together with clusters 3 and 4. Using these four clusters, we found that peptides from cluster 1 & 2, having lower phosphorylation stoichiometry (PIRS) ( Fig. 2A), were indeed enriched in phosphothreonines (Fig. 2B), and somewhat on basophilic kinase targets (Fig. S4). However, acidophilic kinase targets were enriched in PARS peptides ( Fig. 2C and Fig. S4) that had higher phosphate stoichiometry in unperturbed conditions ( Fig. 2A). Finally, we observed that peptides having three or four phosphates were enriched among PARS clusters, whereas peptides with only one phosphosite were associated with PIRS clusters, and with cluster 3 exhibiting less prominent PAIP effects. Collectively, these results (Figs. 1 and 2) are consistent with the previously identified determinants of phosphorylation stoichiometry (1, 40 -44). In the context of PP2A-mediated phosphoregulation, these data show that a significant fraction of serine/threonine phosphosites are primarily sensitive either to PP2A inhibition or activation, but not to both. Furthermore, we reveal that phosphosite amino acid sequence can be predictive of substrates regulated by different PAIPs. As an example, the SET-dominated cluster 5 & 6 ( Fig. 1F) were most strongly enriched for acidophilic kinase targets ( Fig. 2C and Fig. S4) and for peptides with several clustered phosphosites (Fig. 2D).

PP2A-mediated control of cellular processes
To investigate the biological processes that are affected in cells with long-term PP2A modulation, we analyzed the enrichment of pathways and biological processes using the g:Profiler software (46). PP2A-modulated targets were found in a large number of cellular processes and functions, including regulation of translation, RNA metabolism, DNA repair, transcription, heat response, apoptosis, and host-virus interactions (Fig.  3A). Interestingly, most of phosphopeptides from the similarlybehaving cluster (Fig. 1F) are associated with the same cellular processes, indicating that proteins involved in the same biological process respond similarly to PP2A modulations (Fig. 3A). To exclude the possibility of technical variation affecting the clustering analysis, we examined whether selective association of cellular processes with either PP2A or PAIP inhibition was also observed by using significantly-regulated individual peptides. Consistent with the analysis with clustered phosphopeptides, most of the cellular processes were associated selectively either by PPP2R1A or PAIP inhibition, but not by both (Fig. S5). The data also identified processes selectively regulated by SET, PME-1, and CIP2A ( Fig. 3B and Fig. S5), providing clear indications for nonredundant cellular functions for PAIPs.
Interestingly, in contrast to most other cellular functions that responded primarily to either PPP2R1A or PAIP inhibition ( Fig.  3A and Fig. S5), mRNA processing (including e.g. Regulation of RNA stability and RNA processing and splicing) was associated with both inhibition and activation of PP2A ( Fig. 3A and Fig.  S5). Nevertheless, when assessed at the level of individual pro-teins, PP2A inhibition almost exclusively increased phosphorylation of cancer-related splicing factors (e.g. HNRNPA1, SRSF3, and RBM4) (47), whereas a distinct set of splicing factors was found to be dephosphorylated upon depletion of the PAIPs (Fig. 4A). RBM25 was the only target in which the same peptide (Ser-677/Ser-683) was subject to regulation by both PPP2R1A and PAIPs. Therefore, the division of PP2A-modulated target proteins to either as PIRS or PARS, holds true also for splicing factors, albeit as a process RNA splicing is most likely responsive to both PP2A inhibition and activation through these nonoverlapping targets.
As RNA-binding proteins were profoundly represented in PP2A-modulated targets, and RNA processing has an important role in cancer (48), we looked for an example where PP2Amediated phosphorylation changes might impact the cancerrelevant function of an RNA-binding protein. Nucleophosmin (NPM) is a nucleolar RNA-binding protein implicated in cancer (49). In our data, serines 254 and 260 in the C-terminal domain of NPM were significantly regulated by PP2A inhibition ( Fig. 4B). Intriguingly, upon transient transfection WT and double-alanine mutants of NPM-GFP displayed the prototypical nucleolar localization pattern in MDA-MB-231 cells; phosphomimic aspartate mutations of these PP2A-regulated sites induced similar translocation of NPM to nucleoplasm as is seen with the oncogenic mutations of the NPM C-terminal domain (Fig. 4, B and C) (49). As a control, all exogenous NPM proteins were tested to be expressed at comparable protein levels (Fig. S6A). Thus, in cancer cells with inhibited PP2A, but are genetically WT for NPM, PP2A inhibition-induced hyperphosphorylation of NPM could result in similar oncogenic consequences as with NPM C-terminal mutations (49).
Enrichment of canonical signaling pathways and kinase targets between different clusters was additionally evaluated by Ingenuity Pathway Analysis software and the NetworKIN kinase target prediction tool (50). The target proteins in each pathway predicted to be regulated based on NetworKIN are listed in Table S2. Interestingly, the PIRS and PARS clusters associated with distinct canonical signaling pathways. Among the kinase signaling pathways, PTEN signaling, AKT targets, and PKC targets were most clearly associated with PIRS clusters, whereas ERK signaling associated predominantly with PARS cluster 3 (Fig. 4, D and E). ERK pathway inhibition by PAIP inhibition is consistent with recent cancer data using either PME-1 siRNA (51) or PP2A-reactivating compounds (17). In contrast to PIRS clusters enriched in most kinase-signaling pathways, PARS clusters were associated with processes related to DNA repair. Canonical DNA repair-associated pathways and DNA-PK targets were enriched in clusters 4 -6 ( Fig. 4, D and E). DNA-PK was also among the most enriched kinases predicted to phosphorylate peptides with multiple phosphorylation sites (Fig. S6B). The analysis indicates that DNA damage-response proteins are dephosphorylated upon PAIP depletion; however, their phosphorylation cannot be increased from basal line levels by PPP2R1A inhibition. This finding is consistent with recent studies characterizing DNA damage response as a cellular process regulated by dephosphorylation of constitutively-phosphorylated sites (52). Furthermore, a Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes recent study demonstrated that PP2A activity counteracts DNA damage response in yeast (53).
Collectively, PP2A modulations were found to impact the proteome in diverse cellular pathways and processes. Consistent with differential regulation of individual target peptides, different cellular functions and pathways can also be divided based on whether they primarily respond to PPP2R1A or PAIP inhibition (see summary scheme in Fig. 7).

Cytoplasm-nuclear gradient of PP2A modulation
Pathway analysis of PP2A-modulated targets revealed an interesting observation that cytoplasmic kinase signaling asso-ciated with PIRS clusters, whereas nuclear processes seemed to associate with PARS clusters (Figs. 3A and 4, D and E). PPP2R1A localizes to both the cytoplasm and nucleus; CIP2A is predominantly cytoplasmic, and both SET and PME-1 are nuclear (Fig. 5A). Thus, we postulated that differential regulation of cytoplasmic and nuclear phosphoproteome targets could be at least partly linked to differential localization of PAIPs. To examine this uncharacterized aspect of PP2A biology, we assessed the subcellular protein localization of the phosphoproteome target proteins in different clusters. Whereas the targets of the PIRS cluster appeared to be indeed The majority of enriched cellular functions was associated predominantly with only one cluster, indicating that biological functions can also be roughly divided to be responsive to either PP2A inhibition or activation. B, five most significantly-enriched biological process GO terms based on phosphoproteomes significantly (q Ͻ0.05) regulated by CIP2A, PME-1, or SET.    Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes dominantly cytoplasmic, the PARS clusters were highly enriched in nuclear proteins (Fig. 5B). Nuclear PME-1 and SET had the most prominent effects on clusters 4 -6 ( Fig. 1F), and these clusters were particularly enriched for nuclear proteins (Fig. 5B). In contrast, the predominantly cytoplasmic CIP2A had a less potent effect in clusters 4 -6. The correlation between the nuclear localization of PME-1 or SET and their nuclear target proteins was also apparent at the level of the differentially-regulated peptides (Fig. 5C).

Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes
To assess how differential subcellular distribution of processes associated with PIRS and PARS clusters shows up at the individual target proteins, we visualized PP2A-modulated targets in either the cytoplasm-nuclear border or in nuclear protein complexes. As indicated previously by g:Profiler analysis of cellular processes (Fig. 3A and Fig. S5), nuclear lamina organization was one of the rare processes that responded to both PP2A and PAIP inhibition indicating for nearly equal kinase and phosphatase activities at the nuclear-cytoplasmic interphase. Previously, PAIPs CIP2A (54) and PME-1 (55) were shown to interact with one of the main constituents of nuclear lamina LMNA, and here we further validated the physical proximity of PPP2R1A with LMNA on the nuclear envelope using proximity ligation assay (Fig. 5D, inset). Phosphorylation of LMNA at phosphosites Ser-22, Ser-628, and Ser-632 was indeed significantly affected by both PP2A inhibition and activation (Fig. 5D). In addition, we identified regulation of LMNA Ser-3 and Ser-652 only by PPP2R1A and of Ser-390 and Ser-392 only by SET (Table S1). In addition to LMNA, LMNB1 and LMNB2, also the nuclear pore proteins NUP98 and NUP210 were regulated by both PP2A activation and inhibition at the same phosphosite (Fig. 5D). Additionally, the inhibition of PAIPs resulted in dephosphorylation of the nuclear lamina anchor proteins LAP2A and LAP2B (Fig. 5D). To validate the MS data, LMNA phosphorylation was shown by Western blotting to be inhibited by either PME-1 siRNA or CRISPR/Cas9 targeting (Fig. S6, C and D).
As an example of nuclear PARS processes, we focused on "chromosome organization," which in g:Profiler analysis was strongly enriched in the SET target proteins (Fig. S5). As a remarkable example, in the p53-associated protein complex, phosphorylation of every individual protein was regulated by SET (Fig. 5E). We validated SET-mediated inhibition of dephosphorylation of one target, Ser-714 of DNMT1, using a phosphospecific antibody in A549 lung cancer cells (Fig. 5E, inset). About 50% of the proteins associated with chromosome organization were also regulated by the nuclear PME-1 and were consistent with the gradient model; both CIP2A and PPP2R1A had fewer targets in these processes (Fig. 5E).
In summary, these data reveal a cytoplasm-nuclear gradient of PP2A activity. Based on data, the highly-phosphorylated peptides in clusters 4 -6 ( Fig. 2A), which are enriched both on nuclear proteins (Fig. 5B) and CK2 targets (Fig. 4E), are most responsive to SET and PME-1 inhibition (Fig. 5C). Cytoplasmic peptides instead were more frequently found to be regulated by CIP2A (Fig. 5C) and to be associated with low phosphate occupancy ( Fig. 2A), and PKC target sequences (Fig. 4E). These data are consistent with previously shown localization of CK2 and PKC (56,57). Also consistent with the gradient model, nuclear envelope target peptides, between the cytoplasm and nucleus, were found to be subject to phosphoregulation by both PP2A inhibition and activation.

Drug responses as surrogate measure for functional impact of PP2A across diverse cellular pathways
Results above reveal a role for PP2A modulations in a wide range of cellular processes and pathways, many of them involved in cancer cell drug responses. This indicates that PP2A modulations might result in more efficient drug responses. As a systematic approach to assess how globally PP2A affects cellular pathways and processes, we used impact of PP2A modulations on cellular responses to 306 individual drugs as a surrogate functional readout. For this purpose, we used the drug sensitivity and resistance testing (DSRT) platform (58) in which combined effects of drugs and PP2A modulations on cell viability were measured by the CellTiter-Glo assay. By using the same platform, we recently demonstrated a link between PP2A activity and drug responses in KRAS-mutant lung cancer cell across 200 kinase inhibitors (16). The RNAi treatments were identical to those applied for phosphoproteome profiling, and the differential drug-sensitivity score (⌬DSS) (59) was calculated between the treated and control samples. We excluded 68/306 tested drugs that elicited no significant effect on cell viability, regardless of the PP2A modulation (Table S3). The shift in the dephosphorylation balance (Fig. S1H) was used as an indicator of dephosphorylation activity, and the drugs were ranked by uncentered Pearson's correlation between the dephosphorylation activity and the changes in ⌬DSS profile (Fig. 6A). The drugs that exhibited synergy with PAIP depletion appeared at the top of the "correlation rank" (Fig. 6A, red), whereas the drugs at the bottom exhibited synergy with PPP2R1A inhibition (blue). Enrichment scores were calculated for selected drug  (clusters 4 -6). Based on lower frequency of nuclear proteins in cluster 1 & 2, it is apparent that cytoplasmic proteins are under higher dephosphorylation activity than nuclear proteins. C, fraction of nuclear proteins in the proteins with differentially-regulated phosphopeptides for each condition. A and C, **, p Ͻ 0.01, and ***, p Ͻ 0.001, in 2 test. D, differentially-regulated phosphosites on proteins in the enriched category Nuclear envelope organization. Proximity ligation assay of PPP2R1A and Lamin A/C protein association at nuclear lamina is shown in the inset; with PPP2R1A siRNA as a PLA staining specificity control. Red dot indicates close proximity of PPP2R1A and Lamin A/C proteins. Scale bar, 10 m. E, Chromosome organization indicates an example of chromatin-associated protein complexes regulated dominantly by nuclear PP2A inhibitors SET and PME-1. Inset shows Western blot analysis of DNMT1 Ser-714 phosphorylation from A549 cells transfected with SET siRNA for 72 h. * denotes unspecific band. F, schematic presentation of the observed differential cytoplasm-nuclear association of phosphorylation patterns of PP2A targets. Based on data, highly-phosphorylated peptides in clusters 4 -6 that are most responsive to SET and PME-1 inhibition are enriched in nucleus. Highly-phosphorylated peptides associated with cluster 5 & 6 were enriched on CK2 targets. Cytoplasmic peptides instead were more frequently found to be regulated by CIP2A and to be associated with low phosphate occupancy and PKC target sequences. Together, these results indicate for a gradient where cytoplasmic PP2A is more weakly inhibited (e.g. by CIP2A) than nuclear PP2A (e.g. by SET and PME-1).

Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes
groups in the ranked list (Fig. 6A, right panel). All individual drug responses are listed in Table S3.
Notably, PP2A modulations impacted cellular viability responses to drugs targeting a wide range of cellular pathways and cancer-relevant targets (Fig. 6A). Thereby these results provide functional validation of global impact of PP2A on pathway functions. Consistently with previous results implicating PP2A inhibition as a drug resistance mechanism (5,8,14,16,30), PPP2R1A depletion increased average drug resistance across all 238 drugs (Fig. 6, A and B). On the other hand, the most potent induction of protein dephosphorylation by SET depletion (Fig. S1H) correlated with sensitization of cells to all the drugs on average (Fig. 6, A and B). Consistent with the weaker influence of CIP2A and PME-1 on phosphoproteome as compared with SET (Fig. S1H), their effects on the drug sensitivities were in many cases similar to that of SET, albeit of lower magnitudes (Fig. 6A). The most evident correlation between overall dephosphorylation activity and drug response was observed with aurora kinase (AURK) inhibitors (p ϭ 0.0057 for enrichment) and JAK2 inhibitors (p ϭ 0.0142) (Fig. 6, A, C, and  D). PPP2R1A inhibition induced resistance to these drugs, whereas increased dephosphorylation potentiated the drug effects. Consistently, small molecule PP2A activator DT-061 potentiated the effects of the AURK inhibitor barasertib in HeLa and MV4-11 AML cells (Fig. S7A). In line with the broad impact of PP2A modulations on kinase pathway activities (Fig.  4, D and E), PPP2R1A depletion induced resistance also to broad selectivity kinase inhibitor midostaurin, whereas PME-1 and SET depletion increased sensitivity to this staurosporin derivative (Fig. S7B). Dephosphorylation activity also correlated with response to four HSP90 inhibitors with a pattern that was consistent with kinase inhibitor responses (Fig. S7C). This is likely to be related to HSP90 function as a chaperone for a majority of human kinome (60). Reassuringly, PP2A modulations were shown to impact heat-response processes in the pathway analysis ( Fig. 3A and Fig. S5), providing a link between these two systematic PP2A modulation analyses. Furthermore, consistent with the prominent association of PP2A modulations with DNA and RNA synthesis pathways ( Fig. 3A and Fig.  S5), specific drug interactions were observed between PP2A modulations and several nucleoside synthesis inhibitors such as
We found that certain drug responses were modulated only by PPP2R1A inhibition, similarly to its effects on the cellular processes observed above. Induction of resistance by PPP2R1A inhibition was observed for HDAC and bromodomain and extraterminal domain (BET) bromodomain inhibitors, while no sensitization effect with inhibition of any of the PAIPs was seen (Fig. 6, E and F). Notably, targets of these drugs, chromatin remodeling factors, were strongly enriched in PP2A-modulated targets, linking findings from the drug screen to cellular processes affected by PP2A modulations at the phosphoproteome level ( Fig. 3A and Fig. S5). Moreover, the BET inhibitor results are consistent with the recently reported role of PP2A inhibition in JQ1 resistance in triple-negative breast cancer (61). Another interesting but atypical dependence was observed with mTOR that exhibited synergy with PPP2R1A inhibition, in contrast to the majority of drugs we screened (Fig.  6G). The observed sensitization to mTOR inhibitors may reflect an acquired dependence on increased mTOR activity in PPP2R1A-depleted cells, consistent with the up-regulation of PI3K/AKT/mTOR signaling in PPP2R1A-inhibited cells (16). The role of SET in mTOR inhibitor response was particularly interesting because only ATP-competitive mTOR inhibitors exhibited sensitization upon SET depletion (Fig. 6G). This sensitization to kinase inhibitor by SET inhibition can again be explained by phosphoproteomics results, as SET depletion was surprisingly found to activate PI3K/AKT/mTOR signaling (Fig.  S7E).
Altogether, these data provide a systematic functional validation for impact of PP2A modulations on a wide range of cellular pathways and processes. In addition, these findings provide a knowledge resource for the development of combination strategies with the emerging PP2A-targeted therapies aiming either to reactivate (17,18) or inhibit (62) PP2A (Fig. 7). This is particularly appealing as synergy between different drug families and pharmacological PP2A modulations has been recently established (16,63).

Discussion
Despite the essential importance of PP2A for growth and development, and for the pathogenesis of many diseases, our current understanding of the PP2A at the phosphoproteome level remains limited. Here, we characterized global phosphoproteome changes in PP2A-modulated cancer cells by inhibition of the structural PP2A scaffold protein PPP2R1A and three distinct endogenous PP2A-modulating proteins, the functions of which in cancer can be rescued by concomitant PP2A inhibition. Importantly, in addition to specific information related to pathways and processes, our data reveal generally relevant insights into the regulation of serine/threonine phosphorylation. Our biological insights are not only important for understanding cellular phosphosignaling, but also enable translational pharmacological studies to precisely modulate PP2A activity for single and combination therapies in cancer.
In addition to discovering a role for PP2A in a number of cancer-related processes, such as RNA processing (48,49), one of our key discoveries is the differential dominance of either PPP2R1A or PAIP inhibition in the regulation of phosphosites, proteins, and cellular processes (Fig. 7). Although in our study the phosphorylation stoichiometry analysis ( Fig. 2A and Fig. S2) was restricted only to a fraction of observed phosphopeptides, our results are in line with previous reports of global bimodal distribution between dephosphorylated and phosphorylated peptides (1,42). The observed differential amino acid context of peptides regulated either by PP2A inhibition or reactivation is also consistent with previously reported preferential amino acid contexts of phosphosites regulated by kinases and phosphatases (40,41,43,44). Appreciation of different phosphosites differentially regulated by phosphatase inhibition or activation may have in the future pharmacological relevance in predicting the outcomes of pharmacological PP2A targeting (17,18,62). Based on our results, basophilic targets would have a higher probability of responding to pharmacological PP2A inhibition, whereas PP2A reactivation therapies would primarily affect acidophilic target sites (Fig. 7). Our data also imply that the effects of pharmacological phosphatase reactivation cannot be directly deconvoluted from existing kinase inhibitor data. In a simplified model, phosphatase reactivation therapies would primarily affect sites occupied with phosphate, while kinase inhibition would not affect those sites that are already occupied with phosphate. Proof of principle for this model was recently shown in KRAS-mutant lung cancer cells in which PP2A inhibition up-regulated a different set of ERK targets than those that were inhibited by the MEK inhibitor (16).
One of the key unanswered question addressed in this study was (non)redundancy of PAIPs in regulation of cellular phosphoproteome. Previous reports (8, 10, 11, 15, 20, 22, 23, 27, 30 -32) have identified several different functions regulated by CIP2A, PME-1, and SET, but thus far there has not been a systemic analysis of the degree of redundancy in the cellular  Figure 7. Summary scheme of phosphoregulation rules by PP2A modulations. Based on the results, most phosphosites and cellular processes appear to be sensitive to either phosphatase inhibition or activation. In contrast, only few processes show characteristics of bidirectional regulation. The preferable amino acid context for PIRS and PARS sites is also depicted. PIRS and PARS processes show differential spatial distribution between the cytoplasm and nucleus. Based on these results, efficient pharmacological interference with different cellular processes requires understanding whether those processes are either under PIRS or PARS dominance. As an example, phosphatase activation would mostly affect nuclear processes, whereas cytoplasmic processes would be most amenable to modulation by phosphatase inhibition.

Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes
functions they regulate. At the global level, our data demonstrate that among these three PAIPs, SET has the most profound effects on cellular phosphoproteome (Fig. S1, H and I).
The overlapping targets regulated by all of the PAIPs can be explained by direct inhibition of catalytic activity of PP2AC by both SET (27) and PME-1 (24), whereby they would also regulate the CIP2A-targeted PP2A complexes that, by current understanding, are restricted to B56-containing PP2A complexes (26). The nonoverlapping phosphoproteomes of CIP2A and PME-1 can instead be explained by the enzymatic activity of PME-1 as PP2AC C-terminal methylesterase and the role of methylation in promoting B55 subunit binding to PP2A complex (23,39). Functionally differential roles for CIP2A and PME-1 were previously shown to be relevant for staurosporin resistance of glioblastoma cells that was driven by PME-1mediated PP2A inhibition, but was not sensitive to CIP2A modulation (30). These conclusions were confirmed by drug screening in this study that demonstrated sensitization of cells to the staurosporin derivative midostaurin by depletion of PME-1 and SET, but not CIP2A (Fig. S7B). Our data also implicate that nuclear PP2A-modulating proteins may be a far-more important function in chromatin biology than previously appreciated. This could particularly be true for SET, as it associates with histones as a member of the INHAT complex with two other proteins (ANP32A and ANP32E) that also possess PP2A inhibitory activity (28,65,66). In conjunction with the different spatial distribution of PIRS and PARS targets discovered here, and the role of PAIPs in setting this balance (Fig. 5), we envision that SET-mediated PP2A inhibition at chromatin facilitates spatially restricted inhibition of dephosphorylation of chromatin remodeling complexes. Consequently, SET modulation by emerging small molecule approaches (67) could provide novel opportunities for targeting chromatin remodeling and epigenetic complexes for cancer therapy. Mechanistically, the widespread and coordinated dephosphorylation of protein complex components by PP2A may be explained by the recently characterized focal points of PP2A dephosphorylation activity that are facilitated by docking of PP2A to specific target proteins inside multiprotein complexes (44,68,69). This study extended the recently demonstrated widespread dependence of kinase inhibitor responses on PP2A activity (16) across drug classes with diverse mechanisms of action. Besides insights to the design of combinatorial therapies with PP2Amodulating compounds (14,33), our data may also facilitate novel patient therapy stratification strategies based on PPP2R1A mutations (21) or the expression of PAIPs (4) in patient tumors. In support of that, PME1-negative glioma cells were previously found to be hypersensitive to indolocarbazole kinase inhibitors in vitro and in vivo (30), and CIP2A positivity in chronic myeloid cells is a strong predictor of first generation tyrosine kinase inhibitor resistance (8). The current data with phosphoproteome profiles of cells with long-term suppression of either PP2A or PAIPs may thus help identify novel predictive phosphoprotein biomarkers and provide mechanistic explanations for the observed drug hypersensitivities.
In summary, this study provides the largest currently available resource for understanding the target processes affected by modulation of the major serine/threonine phosphatase PP2A. This systematic analysis of the biological redundancy and downstream targets of human PAIPs (CIP2A, PME-1, and SET) also provides the unbiased view on processes that would be most affected by targeting of these oncogenic proteins. We also provide an important resource for understanding how PP2A modulations impact cellular responses across a wide range of drugs and pathway inhibitors. Therapeutic modulation of PP2A has become realistic recently (17,33,62). Therefore, we anticipate that applying our datasets in the development of novel therapeutic agents may enable strategies of precision medicine for multiple human diseases (14).

Cell culture and transfection
HeLa, MDA-MB-231, PC-3, and MV4-11 cells were obtained from ATCC and cultured in Dulbecco's modified Eagle's medium (HeLa and MDA-MB-231) or RPMI (MV4-11) containing 10% fetal bovine serum, 2 mM glutamine, 50 IU/ml penicillin, and 50 g/ml streptomycin. The cells were transfected with Oligofectamine (250 nM siRNA) according to the manufacturer's protocol (Life Technologies, Inc.). CIP2A was targeted with four siRNAs and the other genes with three siRNAs. Knockdown was validated by Western blotting and RT-quantitative PCR. siRNA sequences are presented in Table  S4. PME-1 targeting by Crispr/Cas9 in PC-3 cells has been described in Ref. 70.

Phosphoproteomics
Sample preparation, LC-MS/MS analysis, and the pairwise normalization for label-free quantitative phosphoproteomics were performed according to the previously published protocol (36). For RNAi, one million cells were plated on a 10-cm dish the day before transfection. siRNAs were transfected in a 7.5-ml volume using Oligofectamine and following the manufacturer's protocol. Each siRNA treatment was performed in triplicate for HeLa cells for a total of 9 (PPP2R1A, SET, and PME-1) or 12 (CIP2A) replicate samples per gene. Samples were collected 72 h after the transfection prior to reaching full confluency. Briefly, cells were washed once with ice-cold PBS, scraped into ice-cold PBS, and pelleted in a pre-cooled centrifuge. The pellets were snap-frozen and stored at Ϫ80°C for later processing. Following steps were performed on ice or at 4°C. Cells were lysed with buffer containing 50 mM Tris, pH 7.5, 8 M urea, 2 mM EGTA, 5 mM EDTA, 30 mM sodium fluoride, 60 mM ␤-glycerophosphate, 20 mM sodium pyrophosphate, 1 mM sodium orthovanadate, 5 M pepstatin A, and a Complete protease inhibitor mixture tablet (Roche Applied Science). Lysates were sonicated with Bioruptor (Diagenode) at high intensity for 5 min with 15-s intervals and centrifuged for 35 min at 100,000 ϫ g. Protein concentration of the supernatant was determined by measuring the absorbance at 280 nm All phosphoproteomics sample preparation and LC-MS/MS analysis were performed at Turku Proteomics Facility according to Ref. 36, with some modifications (refer to Fig. S1C). Briefly, the cell lysates (1 mg protein) were mixed with a phosphoprotein standard bovine ␣-casein (10 g, Sigma), reduced with DTT, alkylated with iodoacetamide, and then digested with trypsin (Promega). After acidification with trifluoroacetic Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes acid (TFA), the samples were stored at Ϫ20°C overnight or longer. Aliquots of the digests (10 of 1000 l) were desalted with a microcolumn packed with Empore C18 disk (3M). For LC-MS/MS analysis, the samples (nonenriched digests) were reconstituted in 50 l of 0.1% formic acid (FA), of which 6 l was transferred into an LC sample vial for a 5-l injection. The remaining digests (990 l) were desalted with an Empore C18-SD 10-mm/6-ml cartridge (3M), followed by phosphopeptide enrichment with a microcolumn packed with Sachtopore-NP TiO 2 beads (20 m, 300 Å; ZirChrom). The enriched phosphopeptides were immediately desalted with the C18 microcolumn. For LC-MS/MS analysis, the samples (TiO 2 -enriched phosphopeptides) were reconstituted in 11 l of 0.1% FA, of which 5.5 l was transferred into an LC sample vial for a 5-l injection.
LC-MS/MS analysis was performed using an EASY-nLC 1000 nanoflow LC instrument coupled to a Q Exactive quadrupole-orbitrap mass spectrometer (Thermo Fisher Scientific). carbamidomethylation of cysteine as a fixed modification; oxidation of methionine and phosphorylation of serine/threonine/ tyrosine and acetylation of the protein N terminus as variable modifications; peptide mass tolerance of 5 ppm; and MS/MS ion tolerance of 0.02 Da. For phosphorylation site localization, phospho-RS (version 3.0, the neutral loss option disabled) was enabled. In parallel, peak lists (mgf files) for the TiO 2 -enriched samples were exported.
MS/MS spectra of all possible isoforms of singly phosphorylated peptides were simulated using SimPhospho program based on MS/MS identification of the dephosphorylated peptides (71). As the dephosphorylated peptides, nine LC-MS/MS data acquired in this study and two LC-MS/MS data previously reported (71) were used. Dephosphorylated peptides were searched with Mascot, where PeptideProphet probability of Ն0.95, Mascot expectation value of Յ0.05, and ⌬ score of Ն10 were applied as a cutoff. After the phosphorylation simulation, a consensus spectral library (SimHeLa) was built using Spec-traST (version 4.0, released beta), which consists of 21,922 spectra for 4774 peptide sequences. SpectraST searching against the SimHeLa library (SimSpectraST searching (71)) was performed for the TiO 2 -enriched samples via Proteome Discoverer by submitting the mgf files. The SwissProt database without the reversed sequences was used, and peptide mass tolerance was set to 3 Da. Percolator was enabled for estimating probability. The search was repeated against the library supplemented with decoy entries, using the stand-alone SpectraST (after adjusting the mgf file format) to obtain its F-value and recalculated ⌬Dot score that are required for phosphorylation site localization.
All the Mascot and SimSpectraST results were merged on Proteome Discoverer. An xlsx file was exported after applying a Mascot expectation value of Յ0.05, Percolator PEP of Յ0.05, minimum 7 amino acid length, search engine rank 1, and protein grouping. As cutoffs for 1% false localization rate (i.e. highconfidence phosphorylation sites), phospho-RS probability of Ն99.3% was used, while a SpectraST (standalone) F-value of Ն0.49 and a recalculated ⌬Dot score of Ն0.005 were applied only when a singly-phosphorylated peptide contained Ͼ1 Ser/ Thr/Tyr in its sequence (71). Based on the localization confidence (high or low), a new score was given for each phosphopeptide spectral match. When Mascot and SimSpectraST searches disagreed on high-confidence phosphosites, those were considered as low confidence. The scoring scheme is summarized in Table S5. From the TiO 2 -enriched samples, 11,759 unique phosphopeptides were identified, at a false discovery rate of 0.11% estimated using the target-decoy strategy at a phosphopeptide spectral match level, i.e. 2ϫ decoy matches/ total matches. Label-free quantification was performed using Progenesis LC-MS (version 4.1). The TiO 2 -enriched and -nonenriched samples were processed separately. All the chromatographic data were aligned automatically and further adjusted manually. Peptide ion features were detected in the automatic mode with the highest sensitivity. The features were assigned by importing the xlsx identification file, followed by applying a 5-ppm mass tolerance filter. All the features assigned to human peptides were used for global median centering of abundance ratios.

Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes
Protein abundance was quantified based on the sum of ion abundances of peptides unique to a protein. After annotating phosphorylation sites on proteins, the abundance of phosphosites was quantified based on the sum of ion abundances of phosphopeptide variants (i.e. different charge, missed cleavage, oxidation, and/or acetylation states). When a phosphosite was assigned to both high and low localization confidence features, only the former was taken into account.
The directionality of regulation by PPP2R1A and PAIPs in the HeLa samples was robust to the selection of normalization method ( Fig. 1C and Fig. S7F); however, we observed positive skewing of the fold change distribution by PPP2R1A depletion and negative skewing by depletion of PAIPs ( Fig. 1C and Fig.  S7F). This is likely a result of up-regulation and down-regulation, respectively, of a significant fraction of phosphosites in these samples (Fig. 1C) and supports the selection of pairwise normalization as described in Ref. 36. The directionality of regulation was also consistent across different siRNAs (Fig. S1D). The HeLa cell data consist of two sets analyzed at different times, and the peptides displaying the batch effect were filtered out. Specifically, linear regression model residuals in phosphopeptide sample triplicate abundances displayed bimodality after imputing the smallest nonzero values for nondetected abundances, coupled with log2-transformation, and the phosphopeptides displaying high variance were filtered out before downstream analysis. The differential expression of the filtered dataset was then examined using the popular variance-pooling linear regression R-package limma (72), where each case triplicate was compared with the corresponding scrambled control triplicate

Statistics
Statistical analyses were performed using JMP Pro (version 12.0 for Windows software, SAS Institute Inc.), Graphpad PRISM (version 6, Graphpad Software Inc.), and R (RRID: SCR_001905). Continuous variables were compared using the t test and limma. False discovery rate of limma was controlled by Storey procedure (q value). All the statistical tests were twosided and considered significant at p Յ 0.05. When necessary for the assumption of normal distribution, log transformations were performed for continuous variables. For categorical data, frequencies were compared using the 2 test.

Soft clustering, pathway analyses, and motif enrichment
Soft clustering was performed in R using Mfuzz package (73). After obtaining similar results with a range of different cluster size and fuzzification parameters, we selected the parameters c ϭ 6 and m ϭ 2 for the analysis similarly to Olsen et al. (74). Soft clustering was also performed on a previously released dataset on A549 cells (ProteomeXchange accession number: PXD009900). Sample preparation for these data has been described previously (16), and data were analyzed similarly to HeLa cell data. Pathway enrichment analysis was performed using the g:Profiler R package (46). We selected biological processes of Gene Ontology and molecular pathways of Reactome as data sources for pathway analysis. We used the default method of multiple testing correction and filtered results by corrected significance cutoff of q Ͻ0.05. We also limited gene sets to include between 3 and 500 genes. The background list for pathway enrichment analyses included the unique list of all proteins from our dataset combined with all known human phosphoproteins. The phosphoproteins were retrieved from the PhosphoSitePlus database (75) and were filtered to include previously published results ("MS_LIT" and "LT_LIT"). Protein names in the background set were converted to gene IDs of the Ensembl database using the g:Convert feature of g:Profiler. This procedure provided a list of enriched pathways for each condition. The lists were merged into a unified pathway-condition table and visualized with the Enrichment Map app (76) in Cytoscape (77). Finally, we manually reviewed the Enrichment Map and grouped into the most representative functional themes.
The following filtering steps were performed for pathway enrichment analysis. First, a list of unique proteins was constructed by selecting the peptide with the lowest p value for every protein. Second, unique proteins were filtered using a corrected p value cutoff (q Յ0.05). Third, the final list of proteins was ranked by the significance (p value) for subsequent pathway enrichment analysis. Analysis of phosphopeptide clusters was carried out similarly. Clusters were compiled by assigning proteins to clusters based on maximum probability. For proteins with multiple identified peptides, we selected the peptide with highest cluster membership (Table S1) to represent the protein and thus obtained nonredundant lists of unique proteins. The final lists of proteins were ranked according to the cluster membership probabilities.
Enrichment of canonical pathway components in the clusters, as well as subcellular localization of the proteins, was analyzed with Ingenuity Pathway Analysis, Spring 2015 release version (Qiagen). Kinase target predictions were performed with NetworKIN and NetPhorest (78). Enriched phosphosite motifs were analyzed with motif-x (79). The GO term enrichments shown in Fig. 3B were obtained by STRING database (version 11) (80) using significantly regulated peptides (q Ͻ0.05) from Table S1.

Drug-sensitivity testing
For the high-throughput DSRT analyses, cells were transfected 3 days prior to plating them on the drug-containing 384well plates at a confluency of 1000 cells/well. The subsequent analyses were carried out as described previously (58). Analysis was repeated with three (PPP2R1A, SET, and PME-1) or four (CIP2A) different siRNAs for each gene. The ⌬DSS (59) was calculated by comparing the response of each siRNA-depleted sample to the average of three control samples. In cell-viability experiments with Barasertib (Sigma) and DT-061 (MedChem-Express, HY-112929), 5000 (HeLa) or 50,000 (MV4-11) cells per well were plated on 96-well plates and cultured for 3 days in the presence of drugs and measured with WST1 assay according to the manufacturer's recommendation (Roche Applied Science). After 3 days of culture, cells were incubated for 30 min with WST1 reagent, shaken at 250 rpm for 2 min, and followed by measurement of absorbance at 440 nm. Five independent experiments were performed for both cell lines.
In Western blot assays for PP2A target validation, indicated cell types were seeded on either 6-well plates or 10-cm dishes 24 h prior to drug treatment. Cells were treated for the indi-

Enrichment of drug groups
The numbers of up-regulated phosphopeptides for each gene were subtracted from the number of down-regulated phosphopeptides, and the difference values were used as a PP2A activity index. Uncentered Pearson's correlation coefficients were calculated between the PP2A activity index and DSS drug response values for each drug. Drugs were then ranked by correlation with PP2A activity index. Enrichment scores for selected drug groups in the ranked lists were calculated similarly to Gene Set Enrichment Analysis (64).

Cell stainings
For proximity ligation assay, cells were plated on coverslips 2 days after transfection. On the following day, subconfluent cells were fixed with 3:1 acetone methanol in Ϫ20°C for 10 min, and proximity ligation assay was performed with Duolink kit (Sigma) according to the manufacturer's protocol. Briefly, cells were blocked with blocking solution and incubated in a preheated humidity chamber for 30 min at 37°C, followed by incubation with primary antibodies (PPP2R1A, sc-15355, Santa Cruz Biotechnology; Lamin A/C, sc-6215, Santa Cruz Biotechnology) overnight at 4°C. After the last wash, slides were mounted with Mowiol mounting medium and imaged on an LSM780 confocal microscope (Carl Zeiss). To study the mitotic fraction of PP2A-modulated HeLa cells, the cells were transfected in six-well plates, and 48 h later 20,000 cells were reseeded in chambered ibiTreat -Slide coverslip, which has a hydrophilic surface for optimal cell adhesion (80826, Ibidi). 24 h later cells were fixed using 4% paraformaldehyde (prewarmed at 37°C) for 15 min at room temperature. Permeabilization was done using 1% Triton X-100 (Sigma) in PBS at room temperature for 15 min, and cells were blocked in 10% goat serum (ab7481, Abcam) in PBS for 30 min at room temperature. For visualizing serine 10-phosphorylated histones, cell were incubated overnight at 4°C in rabbit phospho-H3 (Millipore, 06-570) antibody diluted to 1:500 in 10% goat serum. Cells were washed three times in PBS and incubated in secondary antibody, Alexa Fluor 488 goat anti-rabbit IgG (A-11008, Invitrogen) diluted to 1:1000 in 10% goat serum for 1 h at room temperature, and then in 4,6-diamidino-2-phenylindole (D1306, Invitrogen) for 10 min at room temperature. Cells were imaged using confocal microscope (3i CSU-W1 spinning disk) and quantified by software FIJI. Mitotic index was calculated as a percentage of histone positive cells to the total number of cells.

Nucleophosmin mutants
The following primer pairs were used to generate NPM-1 mutants, written in the 5Ј-3Ј direction for forward and reverse primer, respectively: S254A, CAA GCA GCG ATA GAA AAA GGT GGT TCT CTT C and C TTT TTC TAT CGC TGC TTG CAT TTT TGC TTT AAT GTC; S260A, GGT GGT GCG CTT CCC AAA GTG GAA GCC AAA TTC and GGG AAG CGC ACC ACC TTT TTC TAT CGC TGC; S254D, CAA GCA GAT ATA GAA AAA GGT GGT TCT CTT C and C TTT TTC TAT ATC TGC TTG CAT TTT TGC TTT AAT GTC; S260D, GGT GGT GAT CTT CCC AAA GTG GAA GCC AAA TTC and GGG AAG ATC ACC ACC TTT TTC TAT CGC TGC. All constructs were generated by mutagenesis using PCR with Phusion Green Hot Start High-Fidelity PCR Master Mix (F-5665, Thermo Fisher Scientific) and DpnI enzyme from QuikChange Lightning site-directed mutagenesis kit (210518-5, Agilent Technologies). GFP-NPM-1 WT was from Addgene (17578). DH5␣ cells were routinely used for plasmid propagation. For DNA extraction and purification, the following commercial kits were used according to the manufacturer's instructions: NucleoSpin plasmid (740588.50) and NucleoBond Xtra Maxi Plus EF (740426.50, both from Biotop). Samples were verified by sequencing at FIMM (Helsinki, Finland). For cell transfection, the MDA-MB-231 cells were plated in a 12-well plate Systemic analysis of CIP2A, SET and PME-1 phosphoproteomes format and transfected using Lipofectamine 2000 (Thermo Fisher Scientific). Per well, 1 g of DNA and 3 l of Lipofectamine 2000 were diluted in 100 l of Opti-MEM while the cells were left in full growth medium. 24 h post-transfection, the cells were imaged using EVOS cell-imaging system, with a filter for GFP.

Data availability
Raw MS data, protein sequence database, spectral libraries, and search results have been deposited to the jPOST repository with the data set identifier JPST000692 (ProteomeXchange identifier PXD016102).