Polo-like kinase 4 maintains centriolar satellite integrity by phosphorylation of centrosomal protein 131 (CEP131)

The centrosome, consisting of two centrioles surrounded by a dense network of proteins, is the microtubule-organizing center of animal cells. Polo-like kinase 4 (PLK4) is a Ser/Thr protein kinase and the master regulator of centriole duplication, but it may play additional roles in centrosome function. To identify additional proteins regulated by PLK4, we generated an RPE-1 human cell line with a genetically engineered “analog-sensitive” PLK4AS, which genetically encodes chemical sensitivity to competitive inhibition via a bulky ATP analog. We used this transgenic line in an unbiased multiplex phosphoproteomic screen. Several hits were identified and validated as direct PLK4 substrates by in vitro kinase assays. Among them, we confirmed Ser-78 in centrosomal protein 131 (CEP131, also known as AZI1) as a direct substrate of PLK4. Using immunofluorescence microscopy, we observed that although PLK4-mediated phosphorylation of Ser-78 is dispensable for CEP131 localization, ciliogenesis, and centriole duplication, it is essential for maintaining the integrity of centriolar satellites. We also found that PLK4 inhibition or use of a nonphosphorylatable CEP131 variant results in dispersed centriolar satellites. Moreover, replacement of endogenous WT CEP131 with an S78D phosphomimetic variant promoted aggregation of centriolar satellites. We conclude that PLK4 phosphorylates CEP131 at Ser-78 to maintain centriolar satellite integrity.

duplication. Additionally, phosphorylation of CPAP allows for procentriole assembly and centriole elongation, and phosphorylation of CEP135 allows for recruitment of asterless/CEP152 to centrioles (32). Other known substrates of PLK4 mediate centriolar satellite integrity (PCM1 (19)), trophoblast differentiation (HAND1 (33)), and cell motility (APR2 (34)). Additionally, CEP152 (28,35), CDC25C (36), CHK2 (37), and ECT2 (38) have been identified as PLK4 substrates, although the functions of these phosphorylation events remain unclear. Three previous studies have reported MS-based screens to identify PLK4 interactors. One study purified Plx4 from Xenopus egg extracts by affinity chromatography and analyzed its associated binding partners by MS (35), the second study performed PLK4 BioID to identify proximity interactions (39), and the third study immunoprecipitated tagged PLK4 from HeLa cells and examined binding partners by MS (28). To our knowledge, analysis of PLK4-regulated phosphorylation sites in human cells has not been reported. Here, we employed a chemical genetic system to perform an unbiased MS-based analysis of PLK4 phosphorylation targets in nontransformed human cells to discover novel substrates of PLK4. We identify CEP131 as a novel substrate of PLK4 and find that PLK4 phosphorylation of CEP131 is an important contributor to maintaining centriolar satellite integrity.

Generation of PLK4 analog-sensitive cell line
PLK4 is the master regulator of centriole duplication (8,9) and plays additional roles, but few substrates have been identified. To identify new substrates, we engineered a PLK4 conditional knockout cell line and a PLK4 analog-sensitive (AS) cell line. First, recombinant adeno-associated virus was used to insert loxP sites around exons 3 and 4 of one endogenous PLK4 locus. The other genomic locus of PLK4 was deleted, so the genotype of this cell line is PLK4 WT/⌬ . Two clones were identified and validated by PCR (Fig. S1, A and B) and by examining the disappearance of centrosome components by immunofluorescence after treatment with adenoviral Cre (AdCre; Fig. S1, C and D). The efficiency of centrosome loss is similar to a previous report of 39% of cells losing centrosomes 2 days after treatment with PLK4 siRNA (9). As expected in p53 WT cells (40,41), loss of PLK4 significantly reduced colony formation (Fig. S1E).
Next, we used recombinant adeno-associated virus to engineer an AS cell line in immortalized, nontransformed RPE-1 cells. We targeted exons 3 and 4 to generate a point mutation at the conserved gatekeeper residue (L89G), enlarging the ATPbinding pocket for selective inhibition by the bulky ATP analog 3-MB-PP1 (42). To ensure that only PLK4 AS was expressed, the floxed exons 3 and 4 were deleted from the PLK4 WT locus by AdCre treatment, and clones were screened by PCR (Fig. S1, F and G). The genotype of this cell line is PLK4 AS/⌬ , which is henceforth referred to as "AS" throughout.
PLK4 inhibition is expected to prevent new centriole assembly without disassembling preexisting centrioles. Consistent with this, inhibition of AS with 3-MB-PP1 for 5 days yields about 65% of cells with fewer than two centrioles (Fig. 1, A and  B), consistent with previous reports (40,41).
Coincident with centriole loss, there was a proliferation defect in the PLK4 AS cells after 5 days of incubation in 10 M 3-MB-PP1 (Fig. 1, C-E). Importantly, the effect was specific to PLK4 AS cells, as PLK4 WT cells had normal cloning efficiency and proliferated well in the presence of 3-MB-PP1 (Fig. 1, C-E). Further characterization of the proliferation defect revealed cell cycle arrest, but not cell death, as evidenced by negative trypan blue and annexin V staining (Fig. 1F) and positive ␤-gal staining (Fig. 1G). These findings are consistent with previous reports indicating that cells lacking centrioles undergo cell cycle arrest but not apoptosis (40,41). Thus, this PLK4 AS cell line allows for specific inhibition by 3-MB-PP1, resulting in loss of centrioles and cell cycle arrest.

CEP131 is a substrate of PLK4
To identify novel substrates of PLK4 using this system, we employed quantitative high-resolution MS for phosphoproteome analysis of PLK4 AS cells arrested at the G 1 /S boundary ( Fig. 2A) with or without 3-MB-PP1 challenge. Isobaric labeling with 6-plex tandem mass tags (TMTs) allowed us to run three replicates each of AS cells with and without 3-MB-PP1 in a single MS experiment (43)(44)(45)(46). The G 1 /S boundary was selected for this analysis ( Fig. 2A) because centriole duplication occurs at this phase of the cell cycle (14,47). Using this approach, we identified 11,501 phosphopeptides (Table S1), corresponding to 3711 proteins, of which 15 were significantly down-regulated (p value Ͻ0.05 and -fold change Ͼ1.5) by inhibition of PLK4. The most strongly regulated phosphopeptides were found in RUNX1, PTPN12, IL6ST, TRIM3, and SCRIB. Given that none of these proteins are known to localize to the centrosome or play any role in centrosome or microtubule biology, these phosphorylation events are likely to be indirectly regulated by PLK4. To increase the likelihood of identifying direct substrates, we therefore focused on proteins known to localize to the centrosome (Fig. 2B) and reduced our -fold change and statistical significance thresholds, as PLK4 is a lowabundance kinase with rare phosphorylation events that may not reach robust -fold change or statistical significance thresholds. We considered proteins to be centrosomal if their Uniprot entry mentioned localization to the centrosome, centriole, or minus ends of microtubules or if the protein was identified in a phosphoproteomic screen of purified centrosomes (48) (Fig.  2C). Many such phosphopeptides had only modest changes in abundance, possibly due to the moderating effect of feedback regulation of PLK4 catalytic activity (42). Two centrosome proteins identified in our screen, CEP152 and PCM1, were previously identified as PLK4 substrates (Table 1), providing credence that this screening approach identified PLK4 substrates. From the remaining list, we selected the proteins known to be associated with centriole duplication: CEP131 (7,49), CEP215/ CDK5RAP2 (50), and CEP350 (51). None of these proteins have been previously identified as PLK4 substrates. Among the potential phosphorylation sites, we identified residues that are evolutionarily conserved: Ser-78 and Ser-89 on CEP131, Ser-1238 on CEP215, and Ser-1648 on CEP350 (Fig. S2). To evaluate these as potential phosphorylation sites, direct phosphory-lation was tested with GST-tagged 13-mer peptide fragments with the serine of interest located centrally (sequences listed in Table S2). STIL (Ser-1108 Ϯ 6 amino acids) was used as a pos-itive control (29). Of the potential fragments identified, only CEP131 Ser-78 produced a strong signal when incubated with WT PLK4 (Fig. 2D). We also incubated PLK4 with a larger . C, PLK4 AS/⌬ and PLK4 WT/⌬ were subcloned by limiting dilution, and cloning efficiency was assessed using a Poisson correction. D, proliferation was assessed by culturing cells in the indicated concentration of 3-MB-PP1 for 7 days and then staining with crystal violet. E, 1000 cells were plated per well in a 6-well plate (day Ϫ2), and 3-MB-PP1 was added on day 0. Cells were counted every 2 days. Each point represents the average of three technical replicates Ϯ S.E. F, cells were cultured for 7 days with DMSO or 3-MB-PP1, and apoptosis was assessed by both trypan blue exclusion assay and staining with propidium iodide and anti-annexin V followed by flow cytometric detection. The percentage of cells staining negative for both propidium iodide and annexin V is plotted. G, senescence was assessed by ␤-gal staining after 14 days of culture with 3-MB-PP1. The micrographs are representative images of the staining, and the bar graphs quantify the percentage of senescent cells. All experiments in this figure utilized RPE cells. *, p Ͻ 0.05.

PLK4 phosphorylates CEP131
N-terminal fragment of CEP131 (amino acids 1-250) and observe that PLK4 is capable of phosphorylating this fragment (Fig. S3). Importantly, phosphorylation is abolished when the PLK4 inhibitor, centrinone B (40), is added to the reaction (Fig.  2D). Moreover, a CEP131 fragment containing the S78A mutation is not phosphorylated (Fig. 2E). We conclude that CEP131 Ser-78 is directly phosphorylated by PLK4 in vitro. To assess which domain of CEP131 interacts with PLK4, we performed co-immunoprecipitation experiments. Full-length PLK4 was co-expressed with different domains of CEP131. We find that PLK4 interacts with the N terminus of CEP131 (Fig. S5).

CEP131 Ser-78 phosphorylation is required for centriolar satellite integrity but is dispensable for centriole duplication and ciliogenesis
CEP131 is a component of the centrosome and centriolar satellites and was previously identified as a PLK4 interactor in BioID experiments (39). Furthermore, CEP131 and PLK4 colocalize (Fig. S4). CEP131 is known to have three major functions in the cell: centriole duplication, ciliogenesis, and centriolar satellite formation. Moreover, PLK4 activity is implicated in each of these activities (8,9,19). Therefore, we tested whether PLK4 phosphorylation of CEP131 is required for these functions. To do this, we inhibited PLK4 in RPE-1 AS cells with 3-MB-PP1 and in HeLa cells using centrinone B. As expected, inhibition of PLK4 blocked centriole duplication (Fig. S6, A and To test whether CEP131 Ser-78 phosphorylation is important for any of these three functions, we expressed siRNA-resistant, enhanced GFP-tagged, nonphosphorylatable (S78A) or phosphomimetic (S78D) CEP131 mutants in HeLa cells. These mutants properly localize to centrioles and centriolar satellites (Fig. 3A), suggesting that this phosphorylation event is not required for CEP131 localization. Next, we depleted endogenous CEP131 with siRNA (Fig. 3B). Knockdown-addback efficiency was also assessed by quantitative RT-PCR (qRT-PCR) (Fig. 3C).
We next analyzed the ability of these cells to duplicate centrioles and nucleate primary cilia. No defects in centriole numbers were observed among the conditions (Fig. 3, D and E). With regard to primary cilia formation, CEP131-depleted cells exhibited a reduction in primary cilia formation, consistent with previous reports (52,53). However, defects in primary cilia formation did not occur in either the nonphosphorylatable (S78A) or phosphomimetic (S78D) mutant cell lines (Fig. 3, F and G). Combined, these data suggest that PLK4 phosphorylation of CEP131 Ser-78 is dispensable for centriole duplication and primary cilia formation.
PCM1 comprises a structural platform for centriolar satellites (54), and PLK4 phosphorylation of PCM1 is required for A, cell cycle profile of cells arrested with aphidicolin to enrich RPE PLK4 AS cells in the G 1 /S transition, which is when centriole duplication occurs and therefore when PLK4 is active. These cells were utilized for MS. B, volcano plot of Ϫlog 10 (p value) versus log 2 (-fold change) for the centrosome components (red dots) that were identified in the MS screen. C, diagram of the centrosome showing where each of these proteins is thought to localize; derived from Ref. 82. D, GST-tagged peptides and PLK4 kinase domain were purified from Escherichia coli and utilized for in vitro kinase assays. GST serves as a negative control, and STIL serves as a positive control. Centrinone B is a PLK4 inhibitor. E, an S78A mutation was created in the CEP131 peptide and purified for in vitro kinase assay to confirm that PLK4 phosphorylates CEP131 Ser-78.

PLK4 phosphorylates CEP131
proper organization of centriolar satellites (19). To test whether CEP131 Ser-78 phosphorylation is also important for centriolar satellite integrity, we quantified PCM1 as either normal, dispersed, or aggregated, as has been reported previously (19,55). Nocodazole served as positive control for dispersion (56). We find that depletion of CEP131 increases the dispersion of centriolar satellites (Fig. 3, H and I). Additionally, addback of WT and S78D CEP131 rescues this phenotype, but addback of the S78A nonphosphorylatable mutant does not rescue this phenotype. Furthermore, we devised a more objective method of quantifying centriolar satellite organization by computing the percentage of PCM1 intensity found within 2.5 m of the center of the centrosome (Fig. 3J). Similarly, we find that addback of the S78A nonphosphorylatable mutant does not rescue centriolar satellite dispersion (Fig. 3K). Additionally, addback of the S78D phosphomimetic mutant causes increased aggregation of centriolar satellites. These data support a model where phosphorylation of Ser-78 by PLK4 helps restrain dispersion of centriolar satellites. CEP131 Ser-78 was previously reported to be phosphorylated in response to UV irradiation by mitogen-activated protein kinase-activated protein kinase 2 (known as MK2 or MAP-KAP2) (57). In the presence of UV and other cellular stresses, centriolar satellites disperse in a manner dependent on p38mediated MK2 phosphorylation of CEP131 at both Ser-47 and Ser-78, thereby allowing 14-3-3 proteins to bind and sequester CEP131 to the cytoplasm. Given our findings, we tested whether PLK4 is important for this stress-induced remodeling. Indeed, p38 is required for this remodeling in response to UV, but PLK4 activity is dispensable (Fig. 4). This suggests that PLK4 may also regulate centriolar satellite integrity via CEP131 phosphorylation in a manner similar to MK2, albeit under different conditions. CEP131 is known to enhance CEP152 recruitment to the centrosome. CEP152 then recruits WDR62 and CEP63 and promotes the centrosomal localization of CDK2 (7). We therefore tested whether CEP131 Ser-78 phosphorylation alters the recruitment of CEP152 to the centrosome by quantitative im-

PLK4 phosphorylates CEP131
munofluorescence in our knockdown-depletion experiment. However, we did not observe any differences in CEP152 localization to the centrosome (Fig. 5).

CEP131 knockout cell lines exhibit reduced centriolar satellite organization and increased sensitivity to nutrient stresses
To assess the effects of chronic loss of endogenous CEP131, we engineered CEP131 knockout cell lines using CRISPR/Cas9 editing; this is possible as it is not essential (58). As expected, we recovered monoclonal CEP131 knockout cell lines, and confirmed these by CEP131 immunoblotting (Fig. 6A) and immunofluorescence (Fig. 6B). We utilized two different knockout clones for our experiments, denoted KO4E and KO9. As expected, knockout cells had dispersed centriolar satellites, as visualized by PCM1 (Fig. 6, C and D), a decrease in CEP152 expression at the centrosome (Fig. 5B), and an increase in mitotic errors and multinuclei (Fig. S7), consistent with previous reports (7,49). However, no differences were seen in centriole number (Fig. 6E) or formation of primary cilia (Fig. 6, F and G).
Next, we assessed proliferation in CEP131 knockout cell lines and found slightly reduced proliferation (Fig. 6H), as shown previously by CEP131 knockdown (49). Centriolar satellites are important for response to cellular stresses (54,57,59). Therefore, we tested whether our knockout cell lines were more sensitive to nutrient stresses. Indeed, both knockout cell lines are more sensitive to both serum-free medium (Fig. 6I) and 1% serum medium (Fig. 6J), as the CEP131 knockout cell lines demonstrated a greater decrease in cell proliferation in response to these treatments.
To identify the role of phosphorylation at Ser-78, we expressed GFP-tagged WT, S78A, and S78D CEP131 constructs in the CEP131 knockout cell lines and evaluated the integrity of centriolar satellites. Only cells with visible GFP signal were analyzed. Consistent with previous experiments, addback of WT and S78D rescues the dispersed satellites, whereas addback of the S78A mutant does not (Fig. 6, K-M). Additionally, there is an increase in centriolar satellite aggregation in the S78D addback. To confirm that this effect observed using PCM1 as a centriolar satellite marker is indeed a disruption in centriolar satellite integrity versus other causes (e.g. disruption in CEP131 interaction with PCM1), we assessed centriolar satellites in these same conditions using a different marker of centriolar satellites: CEP72. Similarly, we find that CEP72 is dispersed in CEP131 knockout cells and is rescued by addback of WT or S78D transgenes but not S78A (Fig. 7).
Next, we assessed whether this aggregation defect observed in S78D mutants depends on dynein. Previous reports have demonstrated that PCM1 (5, 56, 60) and CEP131 (49) localization to centriolar satellites is dynein-dependent. We expressed WT CEP131 and CEP131 S78D phosphomimetic mutant in

PLK4 phosphorylates CEP131
CEP131 knockout cells and treated with the dynein inhibitor ciliobrevin D. Ciliobrevin D caused dispersion of PCM1 and/or accumulation of PCM1 aggregates outside of the centrosome (Fig. 8, A and B). Further, treatment with ciliobrevin D in CEP131 knockout cells expressing the S78D mutant significantly reduced the aggregation of PCM1 at the centrosome (Fig. 8, C and D). These findings suggest that aggregation of satellites by PLK4 phosphorylation of CEP131 is dynein-dependent.
To verify that these observed effects of CEP131 Ser-78 mutants depend on PLK4 activity, we overexpressed WT, S78A, or S78D constructs and GFP control in WT HeLa cells (Fig. 9A), treated with the PLK4 inhibitor centrinone B, and observed the effects on centriolar satellite organization by PCM1 immunofluorescence (Fig. 9B). Consistent with previous findings, PLK4 inhibition increased the dispersion of centriolar satellites, and the S78D phosphomimetic mutant restored clustered centriolar satellites in the presence of PLK4 inhibition (Fig. 9, C and D).

Additive effects of PLK4 phosphorylation of PCM1 and CEP131 on centriolar satellite organization
The effect of PLK4 phosphorylation of CEP131 on centriolar satellite organization is reminiscent of the previously reported effect of PLK4 phosphorylation of PCM1 (19). Therefore, we tested whether these two phosphorylation events cooperate to control centriolar satellite integrity. We overexpressed both phosphomimetic versions of FLAG-CEP131 (S78D) and GFP-PCM1 (S372D) in WT HeLa cells (Fig. 10A) and treated with centrinone B. We observed an increase in centriolar satellite aggregation in cells co-expressing both CEP131 S78D and PCM1 S372D phosphomimetic constructs compared with cells expressing just one of these phosphomimetic constructs (Fig.  10, B and C), suggesting that phosphorylations of both CEP131 and PCM1 by PLK4 are important for centriolar satellite integrity. In summary, PLK4 phosphorylates CEP131 and PCM1 to maintain the integrity of centriolar satellites in interphase (Fig.  10D).

Discussion
Herein, we have coupled a chemical genetic system of selective PLK4 inhibition and an unbiased, multiplex phosphoproteomic screen to identify novel substrates of PLK4. We identify CEP131 Ser-78 as a substrate of PLK4 and find that PLK4 phosphorylation of this site is important for maintaining the integrity of centriolar satellites in interphase. Centriolar satellites are typically dispersed during mitosis (61), so we speculate that a cellular phosphatase regulates this process by reversing phosphorylation events prior to mitotic onset. However, other regulatory mechanisms are also possible. Additionally, our phosphoproteomic screen has uncovered many other potential substrates that could be explored. This is the first study to our knowledge to characterize cellular knockout of human PLK4 protein using a Cre/lox conditional system versus prior work that focused solely on chemical inhibition or knockdown. However, the phenotype is highly similar to knockdown PLK4 with RNAi (8), auxin-inducible degradation (41), or inhibited PLK4 catalytic activity (29,40). This is most consistent with PLK4 playing a primarily catalytic role without additional structural roles that would impart additional differences in phenotype between chemical inhibition and depletion or knockout.
PLK4 regulates other centrosomal activities independent of centriole duplication, including phosphorylation of PCM1 to

PLK4 phosphorylates CEP131
maintain centriolar satellite integrity and primary cilia formation (19) and phosphorylation of CEP135 to regulate its interaction with CEP152 and influence the radial positioning of CEP152 on centrioles (32). Herein we have contributed to our understanding of one such noncanonical role of PLK4: the organization of centriolar satellites. Centriolar satellites are structures located around the centrosome and are important for proper microtubule nucleation, centrosome assembly, and primary cilia formation (5,54,56). Our data suggest that PLK4 phosphorylation of both PCM1 and CEP131 is required to maintain centriolar satellite organization and/or integrity. However, PLK4 phosphorylation of CEP131 is dispensable for

PLK4 phosphorylates CEP131
centriole duplication and primary cilia formation. Given the previously reported findings that intact, functional centriolar satellites and/or their components are important for centriole duplication (7,49) and ciliogenesis (19,62), it is likely that different methods of centriolar satellite perturbations affect these pathways differently. It will be important to determine whether

PLK4 phosphorylates CEP131
it is centriolar satellites or the individual components themselves that are important for these processes.
A number of centriolar satellite components are important for accurate mitosis. For example, we and others (49) have demonstrated that depletion of CEP131 increases multipolar spindles and chromosome missegregation. Given this increase in chromosomal instability with loss of CEP131, it might be suspected that CEP131 is a tumor suppressor. However, data from the Cancer Genome Atlas indicate that genomic deletion of CEP131 or decreased CEP131 gene expression are extremely rare events in cancer and are likely not drivers of chromosomal instability in cancer. Furthermore, increased CEP131 expression has been reported in hepatocellular carcinoma and correlates with worse outcomes (63). The specific mechanism by which centriolar satellites regulate proper chromosome segregation remains to be elucidated.
As the cell enters mitosis, centriolar satellites disperse, and upon completion of mitosis and entry into G 1 , the satellites reorganize around the centrosome (61). PLK4 phosphorylation of PCM1 and CEP131 may help regulate the reorganization of centriolar satellites in G 1 in preparation for centriole duplication. Previous evidence has suggested that centriolar satellites are required for recruiting many key regulators of centriole duplication to the centrosome. Other possible mechanisms explaining centriolar satellite dispersion at the onset of mitosis include NEK2A kinase-dependent loss of C-NAP1 from the centrosome leading to reduced centriolar satellite density (64). There may be additional mechanisms of centriolar satellite dispersion in mitosis and reorganization in G 1 that remain to be elucidated, and this is an important area of future research. One additional question is why and how some centriolar satellite proteins remain at the centrosome when centriolar satellites

PLK4 phosphorylates CEP131
are dispersed in mitosis and with microtubule depolymerization (e.g. OFD1 and CEP290), whereas others (e.g. BBS4 and PCM1) do not remain at the centrosome under these conditions (6). Further work is required to elucidate this mechanism of centriolar satellite dispersion and aggregation during the cell cycle and its functional significance.
Centriolar satellites are also important for response to stresses (54), such as UV, heat stock, proteotoxic reagents, and nutrient deprivation. In response to serum starvation, MIB1 E3 ubiquitin ligase activity is reduced, which decreases PCM1 and CEP131 ubiquitination (59). Furthermore, loss of centriolar satellites by depletion of PCM1 sensitizes glioblastoma cells to temozolomide (65), a DNA-alkylating agent commonly used to treat glioblastoma. Further work will be needed to further elucidate these molecular mechanisms regarding centriolar satellites and response to cellular stresses. For example, because PLK4 inhibition causes centriolar satellite dispersion, does PLK4 inhibition also sensitize glioblastoma and other cancer types to temozolomide and other stress-inducing agents? As PLK4 inhibition is a potential strategy to treat cancer (40, 66 -68), this is an important area of future research.

PLK4 phosphorylates CEP131
CEP131 at the Ser-47 and Ser-78 residues by the p38 effector kinase MK2 promotes centriolar satellite remodeling and dissociation of CEP131 from PCM1 in response to UV light (57). Surprisingly, this previous report found that expression of both the nonphosphorylatable and phosphomimetic CEP131 prevented dispersion in response to UV light, whereas we have found that the nonphosphorylatable mutant causes centriolar satellite dispersion, whereas the phosphomimetic mutant promotes centriolar satellite aggregation. Although the reason for the discrepancy for the function of the phosphomimetic is unclear, we speculate that the regulatory mechanisms during DNA damage stresses, such as UV, may differ from non-DNA damage cellular stresses.
CEP131 is thought to be involved in centriole duplication or, at least, controlling centriole numbers. However, one report demonstrated that CEP131 depletion increases centriole numbers (49), whereas another report demonstrated that CEP131 depletion decreases centriole numbers (7). Additional evidence from Cep131 mutant MEFs found no differences in centrin foci or centrosome amplification (percentage of cells with more than two ␥-tubulin foci) (52). Furthermore, USP9X is a deubiquitinase that was shown to stabilize CEP131. This prior report also found that USP9X is required for centriole duplication, providing some indirect evidence that CEP131 may be required for centriole duplication (69). However, herein we find that depletion of CEP131 does not significantly alter centriole number, regardless of the degree or duration of knockdown.
A limitation of our experiments is that PLK4 has low abundance with rare phosphorylation events, making it particularly challenging to detect direct substrates. As a result, we reduced our -fold change and statistical significance thresholds and focused on centrosomal proteins. Because of this, any future use of our phosphoproteomic data set to determine potential substrates of PLK4 will also require robust in vitro and in vivo validation. Furthermore, our analysis did not identify some of the previously identified PLK4 substrates, such as STIL. Despite this limitation, we detected other known substrates, including CEP152 and PCM1, as well as CEP131, which is validated as a direct substrate in biochemical assays. Another limitation of our work is that we narrowly focused our potential substrates (e.g. evolutionary conservation, role in centriole duplication), and these criteria may have been too strict to uncover novel PLK4 functions and substrates. Last, the identified phosphosite on CEP131, Ser-78, does not conform to the canonical PLK4 consensus motif (Ser/Thr followed by two hydrophobic residues) (70). However, this consensus motif does not hold in all cases. Despite this, we did validate that direct phosphorylation occurs with purified proteins.
In summary, we developed a chemical genetic method to interrogate PLK4 signaling and utilized multiplex phosphoproteomics to discover novel substrates of PLK4. We identified CEP131 as a novel substrate of PLK4 and found that PLK4 phosphorylation of CEP131 Ser-78 is important for proper organization and integrity of centriolar satellites and may help explain the changes seen in centriolar satellite organization during the cell cycle.

Generation of PLK4 conditional knockout and analogsensitive cell lines
Gene editing was performed using adeno-associated virus in RPE-1 (ATCC, Manassas, VA) an immortalized, nontransformed human epithelial cell line. Procedures for preparation of infectious adeno-associated virus particles, transduction of RPE-1 cells, and isolation of properly targeted clones were performed as described previously (71,72). Briefly, a targeting vector was made by inserting a neo cassette surrounded by FRT sites in the intron upstream of exons 3 and 4 of PLK4, all flanked

PLK4 phosphorylates CEP131
by loxP sites (Fig. S1). To generate adenovirus producing this targeting vector, HEK293 cells were transfected with pRC and pHelper. RPE-1 cells were infected with adenovirus collected from the supernatant and subcloned in the presence of G418 (0.4 g/ml) and screened/verified. This cell line was targeted in two ways: 1) to make a conditional knockout cell line and and 2) to make an analog sensitive cell line. To make the conditional knockout cell line, one copy of PLK4 was deleted by treating with Cre, and then cells were retargeted with the same initial vector (neo cassette, exons 3 and 4, all flanked by loxP sites). To make the analog-sensitive cell line, cells were transfected with pFLIPe to remove the neo cassette, subcloned, and retargeted with the initial targeting vector in which the gatekeeper residue was mutated to a glycine (L89G) to enlarge the ATP-binding pocket and genetically encode chemical sensitivity to the bulky ATP analog 3-MB-PP1 (42) (Fig. S1). The resulting ASflox/ WTflox cell lines were treated with Cre and subcloned in the presence G418. Clones were screened by PCR with primers both inside and outside the neo cassette (Fig. S1). In addition, RNA was isolated from the AS cells and converted to cDNA, and exon 3 was PCR-amplified and sequenced to confirm the presence of only the mutant (L89G).

Cell culture and assays
All cell lines were propagated at 37°C in 5% CO 2 . RPE-1derived cell lines were grown in a 1:1 mixture of Dulbecco's modified Eagle's medium and Ham's F-12 medium supplemented with 2.5 mM L-glutamine, with 10% fetal bovine serum and 100 units/ml penicillin-streptomycin. HeLa cells (ATCC) were grown in highglucose Dulbecco's modified Eagle's medium supplemented with 2.5 mM L-glutamine, with 10% fetal bovine serum and 100 units/ml penicillin-streptomycin. For stable retroviral transduction, constructs were co-transfected with a VSV-G envelope plasmid into Phoenix cells. Fresh medium was applied at 24 h post-transfection, harvested 24 h later, clarified by centrifugation and filtration through a 0.45-m membrane to remove cell debris, and diluted 1:1 with complete medium containing 10 g/ml Polybrene. Target cells were infected at 40 -60% confluence for 24 h. Polyclonal transductants were further purified by limiting dilution to obtain individual clones.
To test for a proliferative defect, subcloning was performed in the presence of 10 M 3-MB-PP1. Poisson correction was used to determine cloning efficiency. Cellular senescence was assayed using a pH-dependent ␤-gal staining kit (Cell Signaling Technology) according to the manufacturer's instructions.
For CEP131 knockdown-addback experiments, HeLa cells were grown to 80 -90% confluence in 6-well plates. 2 g of CEP131 cDNA plasmid and 4 l of 20 M siRNA were co-transfected simultaneously using Lipofectamine 2000 according to the manufacturer's instructions. Cells were given fresh medium 4 h after the transfection. 24 h after the transfection, cells were either transferred to coverslips for fixation and immunofluorescence or harvested for Western blotting and qRT-PCR.
For ciliogenesis, 5 ϫ 10 4 cells were plated on coverslips, allowed to adhere overnight, washed three times with Hank's balanced salt solution, given serum-free medium for 48 h, and then fixed with methanol and stained for acetylated tubulin to mark primary cilia.

CRISPR/Cas9
The CRISPR/Cas9 system was used to create CEP131 knockout cell lines in HeLa cells. Four guide RNAs targeting CEP131 were cloned into lentiCRISPRv2 (73,74) from Feng Zhang (Addgene plasmid 52961). The four guides used were GCG-CTCCGGGACGCTGCCGA, CCGGGACGCTGCCGATG-GCC, TTGACTCACCAGCGGGCCCG, and GGCACTGG-AACTCCGGGCAT. Cutting efficiency of each guide was assessed by a surveyor assay (IDT, 706025) following PCR of 600-bp fragments from genomic DNA encompassing the guide RNA locus. lentiCRISPRv2 vectors were transfected into HeLa cells using Fugene (Promega). Prior to transfection, medium was replaced with serum-and antibiotic-free medium. 24 h after transfection, cells were selected with puromycin (1 g/ml) for 72 h and then allowed to recover for 24 h in fresh medium. The guide RNA yielding the greatest cutting by surveyor assay (TTGACTCACC-AGCGGGCCCG) was transfected again into HeLa cells, which were then subcloned by limiting dilution. Colonies were selected and screened by PCR of the region flanking the cut site and sequencing. Verification of knockout lines was performed with CEP131 Western blotting and immunofluorescence.

Immunofluorescence and microscopy
Immunofluorescence and imaging were carried out as described previously (75,76). Cells were seeded on glass coverslips in 24-well plates and fixed with 100% ice-cold methanol for 30 min. Fixed cells were then blocked for 30 min in 3% BSA and 0.1% Triton X-100 in PBS (PBSTx ϩ BSA). Primary antibodies were incubated in PBSTx ϩ BSA for 1 h at room temperature and washed three times in PBSTx, followed by secondary antibody incubation in PBSTx ϩ BSA for 30 min at room temperature and two washes with PBSTx. Cells were counterstained with 4Ј,6-diamidino-2-phenylindole and mounted on glass slides with Prolong Gold antifade medium (Invitrogen). Image acquisition was performed on a Nikon Eclipse Ti inverted microscope and Hamamatsu ORCA Flash 4.0 camera. Where indicated in the figure legends, optical sections were taken at 0.2-m intervals and deconvolved using Nikon Elements. Where appropriate, the observer was blinded to treatment conditions during image acquisition and analysis. Images were processed and analyzed using Nikon Elements.
Qualitative analysis of centriolar satellites was performed by categorizing cells as either normal, dispersed, or aggregated, as described previously (19). Further, quantitative analysis of centriolar satellites was performed by taking z stack images, deconvolving, and quantifying the percentage of PCM1 or CEP72

PLK4 phosphorylates CEP131
intensity within a 2.5-m radius of the center of the centrosome (identified by ␥-tubulin).

Western blotting and immunoprecipitation
Cells were lysed in buffer (50 mM HEPES, pH 7.5, 100 mM NaCl, 0.5% Nonidet P-40, 10% glycerol) containing phosphatase inhibitors (10 mM sodium pyrophosphate, 5 mM ␤-glycerol phosphate, 50 mM NaF, 0.3 mM Na 3 VO 4 ), 1 mM phenylmethylsulfonyl fluoride, protease inhibitor mixture (Thermo Scientific), and 1 mM DTT. Samples were sonicated and heated in SDS buffer. Proteins were separated by SDS-PAGE, transferred to polyvinylidene difluoride membrane (Millipore), and blocked for at least 30 min in 5% milk and 0.1% Tween 20 in TBS, pH 7.4 (TBST ϩ milk). Membranes were incubated 1 h at room temperature or overnight at 4°C with primary antibodies diluted in TBST ϩ 5% milk, washed three times with TBST, and incubated for 1 h at room temperature in secondary antibodies conjugated to horseradish peroxidase in TBST ϩ 5% milk. Membranes were washed and developed with luminol/peroxide (Millipore) and visualized with film.
For immunoprecipitation experiments, T75 flasks of HeLa cells were transfected with FLAG-PLK4 and GFP-tagged CEP131 fragments. Two days after transfection, cells were lysed in the aforementioned lysis buffer used for Western blotting. Extracts were incubated with anti-FLAG M2 agarose beads (Sigma). Beads were washed three times with lysis buffer, and immunoprecipitated protein complexes were analyzed by SDS-PAGE and immunoblotting with anti-FLAG and anti-GFP antibodies.

Preparation of samples for mass spectrometry
PLK4 AS cells were cultured in T175 plates. At 75% confluence, cells were synchronized at the G 1 -S phase with aphidicolin (5 M) for 24 h. To assess the cell cycle, a fraction of cells from MS sample preparations were taken for flow cytometry. Cell pellets were fixed in cold 70% EtOH for at least 24 h, washed once in PBS, and resuspended in PBS with 0.5 mg/ml RNase A and 50 mg/ml propidium iodide. Samples were incubated at 4°C overnight and analyzed on a flow cytometer (FACSCalibur, BD Biosciences). 3-MB-PP1 (10 M) was added for the last 6 h to half of the plates. Three independent replicates were utilized for each of the two conditions. Cells were collected, pelleted, snap-frozen in liquid nitrogen, and stored at Ϫ80°C until they were ready for lysis. Cell pellets were then lysed in buffer (8 M urea, 50 mM Tris, pH 8.0, 100 mM CaCl 2 ) containing dissolved protease and phosphatase inhibitor tablets (Roche Applied Science) and then sonicated for 20 min. Protein concentrations were determined using a BCA kit (Thermo Pierce). Cell lysates were then reduced by the addition of DTT (final concentration of 5 mM), incubated at 37°C for 45 min, and alkylated by adding iodoacetamide (15 mM final concentration). Next, lysates were incubated for 45 min in the dark at room temperature, and remaining iodoacetamide was quenched by bringing each lysate back to a final 5 mM DTT concentration. Lysates were diluted to a 1.5 M urea concentration using a 50 mM Tris and 100 mM CaCl 2 solution. Trypsin was added to each lysate in a 50:1 (protein/enzyme) ratio and digested overnight at ambient temperature. The six samples were desalted using 100-mg C18 Sep-Pak cartridges (Waters) and dried down using a vacuum centrifuge to obtain tryptically digested peptides. 1 mg of peptides for each of the six samples was incubated with 6-plex TMT reagents (Thermo Scientific) for 3 h at room temperature. An aliquot of each sample was mixed in a 1:1 ratio and run on an Orbitrap Elite mass spectrometer (Thermo Scientific) to ensure complete TMT peptide labeling. The six samples were mixed in a final 1:1 ratio across all six TMT channels and desalted using a 500-mg C18 Sep-Pak cartridge (Waters) to produce a single pooled sample containing chemically labeled peptides from all six samples. The pooled sample was fractionated using strong cation exchange chromatography to produce 12 total peptide fractions, which were subsequently lyophilized and desalted. The resultant 12 peptide fractions were each enriched using immobilized metal affinity chromatography nickel-nitrilotriacetic acid magnetic agarose beads (Qiagen), leading to 12 final strong cation exchange-fractionated enriched phosphopeptide and unenriched peptide fractions. Each fraction was dried down using a vacuum centrifuge and resuspended in 0.2% formic acid for MS analyses.

LC-MS/MS
Each sample was introduced to an Orbitrap fusion mass spectrometer (Thermo Scientific) during a 90-min nano-LC separation using a nanoAcquity UPLC (Waters). A "Top N" fusion method was used to analyze eluting peptides, using a 15,000 resolving power survey scan followed by MS/MS scans collected at 15,000 resolving power. Peptides were fragmented using higherenergy collisional dissociation (HCD) at a normalized collision energy of 35%. Phosphopeptide fractions were analyzed with 200-ms maximum injection times for MS scans and 120-ms maximum injection times for MS/MS scans, whereas unenriched fractions were analyzed with 100-ms maximum injection times for MS scans and 75-ms maximum injection times for MS/MS scans. Only peptides with charge states from ϩ2 to ϩ8 were selected for MS/MS with an exclusion duration of 30 s. The 12 phosphopeptide samples were run in duplicate.

MS data analysis
Data were searched using Proteome Discoverer version 1.4.1.14 (Thermo Fisher Scientific) with the Sequest search algorithm. Thermo RAW files were searched against a Homo sapiens target-decoy database (UniProt, downloaded November 6, 2014). Peptide and phosphopeptide data sets were searched using a 50-ppm precursor mass tolerance and 0.02-Da fragment tolerance for b-and y-type ions produced by HCD fragmentation. All fractions were searched with static carbamidomethyl of cysteine residues, static TMT 6-plex modifications of peptide N termini and lysines, dynamic methionine oxidation, and dynamic TMT 6-plex modification of tyrosine residues. Phosphopeptide fractions were searched with additional dynamic phosphorylation modifications of serine, threonine, and tyrosine residues. Resulting peptide identifications were filtered to 1% false discovery rate and exported to tabdelimited text files compatible with the COMPASS software suite (77). COMPASS calculated the 6-plex TMT protein and phosphopeptide quantitation for all of the 12 peptide fractions and 12 phosphopeptide fractions. Peptides were mapped back to their parent proteins using COMPASS, and Phospho RS (78) was used to localize phosphorylation to amino acid residues with a fragment tolerance of 0.02 Da, automatically considering neutral loss peaks for HCD and considering a maximum of 200 maximum position isoforms per phosphopeptide. The raw 6-plex reporter ion intensities of localized phosphopeptide isoforms were log 2 -transformed and normalized against inhibited PLK4 (ϩ3-MB-PP1) to obtain the relative phosphopeptide and protein quantitation for each cell line and condition.

Kinase assays
The kinase domain (amino acids 1-390) of human PLK4 was cloned into pGEX-6P-1 vector (GE Healthcare). The catalytically inactive mutant, K41M, was generated using Phusion site-directed QuikChange mutagenesis. Fifteen-amino acid fragments of the identified potential PLK4 substrates (putative phosphosite Ϯ 7 amino acids) were also cloned into pGEX-6P-1. Sequences of these peptides can be found in Table S2. Protein was expressed in Rosetta DE3 bacteria, extracted with GSH Sepharose 4B (GE Healthcare, 17-0756-01), and eluted by the addition of GSH for the production of the intact GST fusion proteins.
All kinase assay reactions were incubated at 30°C for 30 min in buffer (50 mM Tris-HCl, pH 7.5, 10 mM MgCl 2 , 0.1 mM NaF, 10 M Na 3 VO 4 ) with 1 mM DTT, 1 M cold ATP, 2 Ci of [␥-32 P]ATP, 5 g of substrate, and 100 ng of PLK4. Reactions were resolved by SDS-PAGE. Dried gels were exposed to a storage phosphor screen (GE Healthcare), and ␥-32 P incorporation was visualized by a Typhoon TRIO imager (GE Healthcare).
To confirm that CEP131 Ser-78 is a substrate of PLK4, the S78A mutant (15-amino acid fragment, as above) was generated using Phusion site-directed mutagenesis, purified, and utilized in kinase assays, as described above.

Molecular biology
CEP131 cDNA was obtained as a kind gift from Drs. Spencer Collis and Katie Myers (corresponds to Uniprot isoform 2, Q9UPN4-2). Phusion site-directed mutagenesis was used to generate S78A and S78D mutants and to make CEP131 siRNA-resistant (mutated AGGCCCTCAAGGCCAACAA to AAGCATTGAAAGCAAATA).

qRT-PCR
RNA was isolated from cells using TRI Reagent (Molecular Research Center) per the manufacturer's protocol. DNase treatment was used to remove DNA contamination (Promega). RNA was then converted to cDNA using the Applied Biosystems High Capacity cDNA Reverse Transcription Kit and the manufacturer's protocol. Primer sequences are provided in Table S3. Primers were used at a final concentration of 0.67 M each, and 25 ng of cDNA was used per reaction. A StepOne Plus (Applied Biosystems) real-time PCR thermal cycler was used for amplification with iQ SYBR Green Supermix (Bioard) per the manufacturer's protocol. Quantification of CEP131 mRNA was normalized to three housekeeping genes (RRN18S, GAPDH, and ACTB). The ⌬⌬C T method was employed to calculate the -fold change in expression (79).

Drug screening
To assess the sensitivity of CEP131 knockout cell lines to chemicals, 1000 cells/well were plated in 96-well plates and allowed to adhere for 24 h. Chemicals were added the next day. Four days later, medium was removed, and cell viability was determined using Vita-Orange cell viability reagent (Biotool) according to the manufacturer's instructions. The Vita-Orange reagent detects viable, metabolically active cells. Relative proliferation was determined by subtracting absorbance values of wells containing only Vita-Orange reagent from the experimental wells and then normalizing them to untreated control wells for each cell line.

Statistics
Statistical evaluations were performed using Prism software (GraphPad). Two-tailed t tests and one-way analysis of variance were used for comparisons. p values Ͻ 0.05 were considered significant for all tests, and statistical significance is indicated in the figures.