(species) or
resource
All mice were housed with food and water ad libidum in transparent plastic cages (267 mm long ×207 mm wide ×140 mm high; Techniplast Makrolon type 2 1264C001) with a stainless-steel wire lid (Techniplast 1264C116), kept in light- and soundproof ventilated chambers. All mice were entrained to a 12:12 hr LD cycle, and the time of day was expressed as Zeitgeber time (ZT; ZT0 lights on, ZT12 lights off). Two- to four-month-old males were used for the experiments. Housing as well as experimental procedures were performed in accordance with the guidelines of the Schweizer Tierschutzgesetz and the declaration of Helsinki. The state veterinarian of the Canton of Fribourg approved the protocol. The floxed Per2 mice ( Chavan et al., 2016 ) are available at the European Mouse Mutant Archive (EMMA) strain ID EM:10599, B6;129P2-Per2 tm1Ual /Biat.
The SDL screen was essentially performed as described earlier ( Measday et al., 2005 ; Tong, 2001 ). Briefly, the bait strain Y2454 (MATα mfa1Δ::MFA1pr-HIS3, can1Δ, his3Δ1, leu2Δ0, ura3∆0, lys2Δ0 ) carrying the plasmid YCplF2- mPer2 (that drives expression of PER2 from the galactose-inducible GAL1 promoter) was inoculated into 50 mL glucose-containing synthetic dropout medium lacking leucine (SD-Leu) and grown at 30°C overnight with shaking. Cells were then centrifuged, resuspended in 20 mL of the supernatant, poured into a sterile rectangular petri dish, spotted in a 96-well format on rectangular SD-Leu plates (coined ‘bait plates’ hereafter) using a Biomek 2000 robot (Beckman Coulter, USA), and then grown at 30°C for 3 days. In parallel, the gene deletion array in the strain BY4741 (MATa his3∆1, leu2∆0, met15∆0, ura3∆0 ) was spotted from the storage plates onto fresh G418-containing YPD plates (96-well format) and also grown at 30°C for 3 days. For the mating procedure (overnight at 30°C), colonies from bait plates were (robotically) spotted onto plates containing YPD (plus adenine) and the colonies from the gene deletion array plates were (each separately and in duplicate) spotted on top of them. The next day, the colonies were transferred to G418-containing SD plates lacking lysine, methionine, and leucine (SD-Lys/Met/Leu/+G418) to select for diploids that harbour the YCplF2- mPer2 plasmid. Diploids were then spotted onto plates containing sporulation medium (10 g L −1 potassium acetate, 1 g L −1 yeast extract, 0.1 g L −1 glucose, 2% w/v agar, supplemented with uracil, histidine, and G418) and incubated at 24°C. After 9 days, tetrads were observed and the colonies were transferred to canavanine-containing SD plates lacking arginine and histidine (SD-Arg/His/+canavanine) to select for MATa haploids. Following growth at 30°C for three days, a second haploid selection was carried out by spotting the colonies on SD-Arg/His/Leu/+canavanine plates (to select for MATa haploids containing the YCplF2- mPer2 plasmid). Following growth at 30°C for 2 days, a third haploid selection was carried out by spotting cells on SD-Arg/His/Leu/+canavanine/+G418 plates (to select for MATa haploids containing the YCplF2- mPer2 plasmid as well as the respective gene deletions of the yeast knockout collection). Following incubation at 30°C for 5 days, colonies were then spotted in parallel onto SD-Arg/His/Leu/+G418 plates and on SD-Raf/Gal-Arg/His/Leu/+G418 plates (containing 1% raffinose and 2% galactose as carbon sources) to induce expression of PER2. Both types of plates were incubated at 30°C for 4 days and photographed every day. Strains that grew significantly less well on SD-Raf/Gal-Arg/His/Leu/+G418 than on SD-Arg/His/Leu/+G418 included eap1∆ , gnd1∆ , and pho85∆ . In control experiments, the respective original yeast knockout collection mutants were transformed in parallel with the YCplF2- mPer2 or the empty YCplF2 plasmid ( Foreman and Davis, 1994 ), selected on SD-Leu plates, grown overnight in liquid SD-Leu, spotted (10-fold serial dilutions) on SD-Raf/Gal-Leu plates, and grown for 3 days at 30° ( Figure 1A ). Please note that all media containing G418 were made with glutamate (1 g L −1 ) instead of ammonium sulfate as nitrogen source, as recommended in Tong (2001) .
Adeno Associate Viruses (AAVs) were produced in the Viral Vector Facility (ETH Zurich). Plasmids used for the production are available on the VVF web site. Two constructs were produced. ssAAV-9/2-hSyn1-chI[mouse(shCdk5)]-EGFP-WPRE-SV40p(A) carried the shRNA against Cdk5 (shD, see Figure 2—figure supplement 1 and Supplementary file 2 ) which knocked down only neuronal Cdk5 . ssAAV-9/2-hSyn1-chI[1x(shNS)]-EGFP-WPRE-SV40p(A) was the scrambled control.
Stereotaxic injections were performed on 8-week-old mice under isofluorene anaesthesia using a stereotaxic apparatus (Stoelting). The brain was exposed by craniotomy and the Bregma was used as reference point for all coordinates. AAVs were injected bilaterally into the SCN (Bregma: anterior-posterior (AP) − 0.40 mm; medial-lateral (ML) ±0.00 mm; dorsal-ventral (DV) – 5.5 mm, angle + /- 3°) using a hydraulic manipulator (Narishige: MO-10 one-axis oil hydraulic micromanipulator, http://products.narishige-group.com/group1/MO-10/electro/english.html ) at a rate of 40 nL/min through a pulled glass pipette (Drummond, 10 µl glass micropipet; Cat number: 5-000-1001-X10). The pipette was first raised 0.1 mm to allow spread of the AAVs, and later withdrawn 5 min after the end of the injection. After surgery, mice were allowed to recover for 2 weeks and entrained to LD 12:12 prior to behavior and molecular investigations.
Analysis of locomotor activity parameters was done by monitoring wheel-running activity, as described in Jud et al. (2005) , and calculated using the ClockLab software (Actimetrics). Briefly, for the analysis of free-running rhythms, animals were entrained to LD 12:12 and subsequently released into constant darkness (DD). Internal period length (τ) was determined from a regression line drawn through the activity onsets of ten days of stable rhythmicity under constant conditions. Total and daytime activity, as well as activity distribution profiles, was calculated using the respective inbuilt functions of the ClockLab software (Acquisition Version 3.208, Analysis version 6.0.36). Numbers of animals used in the behavioral studies are indicated in the corresponding figure legends.
Animals used for the immunohistochemistry were killed at appropriate ZTs indicated in the corresponding figure legends. Brains were perfused with 0.9% NaCl and 4% PFA. Perfused brains were cryoprotected by 30% sucrose solution and sectioned (40 µm, coronal) using a cryostat. Sections chosen for staining were placed in 24-well plates (two sections per well), washed three times with 1x TBS (0.1 M Tris/0.15 M NaCl) and 2x SSC (0.3 M NaCl/0.03 M tri-Na-citrate pH 7). Antigen retrieval was performed with 50% formamide/2x SSC by heating to 65°C for 50 min. Then, sections were washed twice in 2x SSC and three times in 1x TBS pH 7.5, before blocking them for 1.5 hr in 10% fetal bovine serum (Gibco)/0.1% Triton X-100/1x TBS at RT. After the blocking, the primary antibodies, rabbit anti-PER2-1 1:200 (Alpha Diagnostic, Lot numb. 869900A1.2-L), mouse anti-Cdk5 clone 2H6 1:20 (Origene, Lot numb. A001), and rabbit anti-GFP 1:500 (abcam ab6556) diluted in 1% FBS/0.1% Triton X-100/1x TBS, were added to the sections and incubated overnight at 4°C. The next day, sections were washed with 1x TBS and incubated with the appropriate fluorescent secondary antibodies diluted 1:500 in 1% FBS/0.1% Triton X-100/1x TBS for 3 hr at RT. (Alexa Fluor 488-AffiniPure Donkey Anti-Rabbit IgG (H+L) no. 711–545–152, Lot: 132876, Alexa Fluor647-AffiniPure Donkey Anti-Mouse IgG (H+L) no. 715–605–150, Lot: 131725, Alexa Fluor647-AffiniPure Donkey Anti-Rabbit IgG (H+L) no. 711–602–152, Lot: 136317 and all from Jackson Immuno Research). Tissue sections were stained with DAPI (1:5000 in PBS; Roche) for 15 min. Finally, the tissue sections were washed again twice in 1x TBS and mounted on glass microscope slides. Fluorescent images were taken by using a confocal microscope (Leica TCS SP5), and images were taken with a magnification of 40x or 63x. Images were processed with the Leica Application Suite Advanced Fluorescence 2.7.3.9723 according to the study by Schnell et al. (2014) .
Immunostained sections were quantified using ImageJ version 1.49. Background was subtracted and the detected signal was divided by the area of measurement. An average value obtained from three independent areas for every section was used. The signal coming from sections obtained from silenced mice was quantified as relative amount to the scramble, which was set to 1. Statistical analysis was performed on three animals per treatment.
NIH3T3 mouse fibroblast cells (ATCCRCRL-1658) were maintained in Dulbecco's modified Eagle's medium (DMEM), containing 10% fetal calf serum (FCS) and 100 U/mL penicillin-streptomycin at 37°C in a humidified atmosphere containing 5% CO2. Cdk5 KO cells were produced applying the CRISPR/Cas9 technique according to the manufacturer’s protocol of the company (Origene, SKU # KN303042).
The following plasmids used were previously described: pSCT-1, pSCT-1mPer2, pSCT-1 mPer-V5, pSCT1 ΔPasA mPer2 -V5, pSCT1 ΔPasB mPer2 -V5 ( Langmesser et al., 2008 ) ( Schmutz et al., 2010 ). pSCT-1 Cdk5-HA, pet-15b Cdk5-HIS, Gex-4T Per2 1–576, pGex-4T Per2 577–1256 were produced for this paper. The full-length cDNA (or partial fragments) encoding PER2 and the full-length Cdk5 were previously sub-cloned in the TOPO vector according to the manufacturer’s protocol (Catalog numbers pCR2.1-TOPO vector: K4500-01). They were subsequently transferred into the plasmid pSCT-1 using appropriate restriction sites. pGex-4T Per2 1–576 S394G, S394D, pSCT-1 mPer2 S394G were obtained using site-specific mutagenesis according to the manufacturer’s protocol on the requested codon carrying the interested amino acid of interest (Agilent Catalog # 200518). For accession numbers, vectors, mutations, and primers sources, see Supplementary file 2 .
NIH 3T3 cells were transfected in 10 cm dishes at about 70% of their total confluency using linear polyethylenimine (LINPEI25; Polysciences Europe). The amounts of expression vectors were adjusted to yield comparable levels of expressed protein. Thirty hours after transfection, protein extracts were prepared. With regard to immunoprecipitation, each antibody mentioned in the paper was used in the conditions indicated by the respective manufacturer. The next day, samples were captured with 50 µL at 50% (w/v) of protein-A agarose beads (Roche) at 50% (w/v) and the reaction was kept at 4°C for 3 hr on a rotary shaker. Prior to use, beads were washed three times with the appropriate protein buffer and resuspended in the same buffer (50% w/v). The beads were collected by centrifugation and washed three times with NP-40 buffer (100 mM Tris-HCl pH7.5, 150 mM NaCl, 2 mM EDTA, 0.1% NP-40). After the final wash, beads were resuspendend in 2% SDS, 10% glycerol, 63 mM Trish-HCL pH 6.8 and proteins were eluted for 15 min at RT. Laemmli buffer was finally added, samples were boiled for 5 min at 95° C and finally loaded onto 10% SDS-PAGE gels ( Laemmli, 1970 ).
Medium was aspirated from cell plates, which were washed two times with 1x PBS (137 mM NaCl, 7.97 mM Na 2 HPO 4 × 12 H 2 O, 2.68 mM KCl, 1.47 mM KH 2 PO 4 ). Then PBS was added again and plates were kept 5 min at 37°C. NHI3T3 or HEK cells were detached and collected in tubes and washed two times with 1x PBS. After the last washing, pellets were frozen in liquid N 2 , resuspended in Ripa buffer (50 mM Tris-HCl pH7.4, 1% NP-40, 0.5% Na-deoxycholate, 0.1% SDS, 150 mM NaCl, 2 mM EDTA, 50 mM NaF) with freshly added protease or phosphatase inhibitors, and homogenized by using a pellet pestle. After that samples were centrifuged for 15 min at 16,100 g at 4°C. The supernatant was collected in new tubes and pellet discarded.
Total brain or isolated SCN tissue was frozen in liquid N 2 , and resuspended in lysis buffer (50 mM Tris-HCl pH 7.4, 150 mM NaCl, 0.25% SDS, 0.25% sodium deoxycholate, 1 mM EDTA) and homogenized by using a pellet pestle. Subsequently, samples were kept on ice for 30 min and vortexed five times for 30 s each time. The samples were centrifuged for 20 min at 12,000 rpm at 4°C. The supernatant was collected in new tubes and the pellet discarded.
Tissues or cells were resuspended in 100 mM Tris-HCl pH 8.8/10 mM DTT and homogenized with a disposable pellet pestle. After 10 min incubation on ice, the samples were centrifuged at 2500 g for 2 min at 4°C and the supernatant discarded. After adding 90 μL of completed cytoplasmic lysis buffer (10 mM EDTA, 1 mM EGTA, 10 mM Hepes pH 6.8, 0.2% Triton X-100, protease inhibitor cocktail (Roche), NaF, PMSF, ß-glycerophosphate), the pellet was resuspended by vortexing, followed by centrifugation at 5200 rpm for 2 min at 4°C. The supernatant transferred into a fresh 1.5 mL tube was the CYTOPLASMIC EXTRACT. The pellet was washed three times with cytoplasmic lysis buffer and resuspended in 45 μL 1x NDB (20% glycerol, 20 mM Hepes pH 7.6, 0.2 mM EDTA, 2 mM DTT) containing 2x proteinase and phosphatase inhibitors followed by adding 1 vol of 2x NUN (2 M Urea, 600 mM NaCl, 2% NP-40, 50 mM Hepes pH 7.6). After vortexing the samples were incubated 30 min on ice, centrifuged 30 min at 13,000 rpm at 4°C and the supernatant that was transferred into a fresh tube was the NUCLEAR EXTRACT.
A protein amount corresponding to between 400 and 800 µg of total extract was diluted with the appropriate protein lysis buffer in a final volume of 250 µL and immunoprecipitated using the indicated antibody (ratio 1:50) and the reaction was kept at 4°C overnight on a rotary shaker. The day after, samples were captured with 50 µL of 50% (w/v) protein-A agarose beads (Roche) and the reaction was kept at 4°C for 3 hr on a rotary shaker. Prior to use, beads were washed three times with the appropriate protein buffer and resuspended in the same buffer (50% w/v). The beads were collected by centrifugation and washed three times with NP-40 buffer (100 mM Tris-HCl pH7.5, 150 mM NaCl, 2 mM EDTA, 0.1% NP-40). After the final wash, beads were resuspendend in 2% SDS 10%, glycerol, 63 mM Trish-HCL pH 6.8 and proteins were eluted for 15 min at RT. Laemmli buffer was finally added, samples were boiled 5 min at 95° C and loaded onto 10% SDS-PAGE gels.
GST-fused recombinant Per2 proteins were expressed in the E. coli Rosetta strain [plasmids: GST-Per2 N-M (1-576), GST-Per2 M-C (577-1256)]. Proteins were induced with 1 mM IPTG (Sigma-Aldrich) for 3 hr at 30°C. Subsequently, proteins were extracted in an appropriate GST lysis buffer (50 mM Tris-Cl pH 7.5, 150 mM NaCl, 5% glycerol) adjusted to 0.1% Triton X-100 and purified to homogeneity with glutathione-agarose beads for 2 hr at 4°C. The beads were then incubated overnight at 4°C and washed with GST lysis buffer adjusted with 1 mM DTT. Subsequently, elution with 10 mM reduced glutathione took place for 15 min at room temperature. Elution was stopped by adding Laemmli buffer and samples were loaded onto the gel after 5 min at 95°C and WB was performed using anti-GST (Sigma no. 06–332) and anti-HA antibodies (Roche no. 11867423001) for immunoblotting.
The CRISPR/Cas9 Cdk5 cell line was produced starting from NIH3T3 cells using a kit provided by Origene ( www.origene.com ). The knock-out cell line was produced according to the manufacturer’s protocol. Briefly, cells at 80% of confluency were co-transfected with a donor vector containing the homologous arms and functional cassette, and the guide vector containing the sequence that targets the Cdk5 gene. In parallel, a scrambled negative guide was also co-transfected with a donor vector. 48 hr after transfection the cells were split 1:10 and grown for 3 days. Cells were split another seven times (this time is necessary to eliminate the episomal form of donor vector, in order to have only integrated forms). Then, single colonies were produced and clones were analyzed by PCR in order to find those clones that did not express Cdk5 RNA. Positive clones were re-amplified.
PCR primers for genomic Cdk5:
FW: 5’- tgtgagtaccacctcctctgcaa -3’
RW: 5’- ttaaacaggccaggcccc -3’
About 24 hr after seeding cells, different shRNA Cdk5 plasmids (Origene TL515615 A/B/C/D Cdk5 shRNA) were transfected to knock down Cdk5 according to the manufacturer’s instructions. The knock-down efficiency was assessed at 48 hr post transduction by western blotting. Scrambled shRNA plasmid (Origene TR30021) was used as a negative control.
NIH3T3 cells were treated with 100 µM cycloheximide 48 hr after transfection with the indicated vectors, and cells were harvested 0, 3, and 6 hr after treatment.
About 48 hr after transfection with either scrambled or shCdk5, cells where Cdk5 was silenced were treated for 12 hr with either DMSO (vehicle) or epoxomicin (Sigma-Aldrich) at a final concentration of 0.2 µM. Samples were collected, and proteins extracted followed by western blotting.
Recombinant GST-fused PER2 protein fragments were expressed and purified from the BL21 Rosetta strain of E. coli according to the manufacturer’s protocol described before, using glutathione-sepharose affinity chromatography (GE Healthcare). Each purified protein (1 µg) was incubated in the presence or absence of recombinant Cdk5/p35 (the purified recombinant N-terminal His6-tagged human Cdk5 and N-terminal GST-tagged human p25 were purchased from Millipore). Reactions were carried out in a reaction buffer (30 mM Hepes, pH 7.2, 10 mM MgCl2, and 1 mM DTT) containing [γ- 32 P] ATP (10 Ci) at room temperature for 1 hr and then terminated by adding SDS sample buffer and boiling for 10 min. Samples were subjected to SDS-PAGE, stained by Coomassie Brilliant Blue, and dried, and then phosphorylated proteins were detected by autoradiography.
CDK5 was immunoprecipitated from SCN samples at different ZTs (circa 600 µg of protein extract) ( Figure 8 ). After immunoprecipitation at 4°C overnight with 2x Protein A agarose (Sigma-Aldrich), samples were diluted in washing buffer and split in two halves. One half of the IP was used for an in vitro kinase assay. Briefly, 1 µg of histone H1 (Sigma-Aldrich) was added to the immunoprecipitated CDK5 and assays were carried out in reaction buffer (30 mM Hepes, pH 7.2, 10 mM MgCl 2 , and 1 mM DTT) containing [γ- 32 P] ATP (10 Ci) at room temperature for 1 hr and then terminated by adding SDS sample buffer and boiling for 5 min. Samples were subjected to 15% SDS-PAGE, stained by Coomassie Brilliant Blue, and dried, and then phosphorylated histone H1 was detected by autoradiography. The other half of the IP was used for Western blotting to determine the total amount of CDK5 immunoprecipitated from the SCN samples. To quantify the kinase activity at each time point, the following formula was used: ([ 32 P] H1/total H1 for each reaction)/CDK5 IP protein.
Workflow of the in vitro kinase assay performed using immunoprecipitated CDK5 from SCN protein extracts is schematized here. Seven mice were sacrificed, SCN tissues were isolated and pooled together every 4 hr starting from ZT 0 (lights on) until ZT20 (ZT12 lights off). Total protein was obtained from each pool of tissues, the quality of the extracts was checked by WB, and subsequently CDK5 was immunoprecipitated at each time point. Agarose beads detained the immunoprecipitation and one half of the precipitate was used for an in vitro kinase assay using as substrate commercial histone H1 as substrate. The other half was analyzed by WB in order to quantify the amount of protein immunoprecipitated, which was used for the kinase assay. Kinase activity around the clock was quantified using the following formula: ( 32 P-H1/total H1)/amount of immunoprecipitated CDK5.
Filter-aided in vitro kinase assays and mass spectrometry analyses were performed essentially as described ( Hatakeyama et al., 2019 ). Briefly, recombinant Cdk5/p35 (Millipore) was incubated with the GST-fused PER2 protein fragment. On 10 kDa MW-cutoff filters (PALL) samples were incubated in kinase buffer containing 50 mM Hepes, pH 7.4, 150 mM NaCl, 0.625 mM DTT, Phostop tablets (Roche), 6.25 mM MgCl 2 , and 1.8 mM ATP at 30°C for 1 hr. Samples without ATP were used as negative control. Assays were quenched by 8 M urea and 1 mM DTT. Protein digestion for MS analysis was performed overnight ( Wiśniewski et al., 2009 ). Phosphopeptides were enriched by metal oxide affinity enrichment using titanium dioxide (GL Sciences Inc, Tokyo, Japan) ( Zarei et al., 2016 ).
LC-MS/MS measurements were performed on a QExactive Plus mass spectrometer coupled to an EasyLC 1000 nanoflow-HPLC. Peptides were separated on fused silica HPLC-column tip (I.D. 75 µm, New Objective, self-packed with ReproSil-Pur 120 C18-AQ, 1.9 µm [Dr. Maisch, Ammerbuch, Germany] to a length of 20 cm) using a gradient of A (0.1% formic acid in water) and B (0.1% formic acid in 80% acetonitrile in water): loading of sample with 0% B with a flow rate of 600 nL min-1; separation ramp from 5–30% B within 85 min with a flow rate of 250 nL min-1. NanoESI spray voltage was set to 2.3 kV and ion-transfer tube temperature to 250°C; no sheath and auxiliary gas was used. Mass spectrometers were operated in the data-dependent mode; after each MS scan (mass range m/z = 370–1750; resolution: 70,000) a maximum of 10 MS/MS scans were performed using a normalized collision energy of 25%, a target value of 1000 and a resolution of 17,500. The MS raw files were analyzed using MaxQuant Software version 1.4.1.2 ( Cox and Mann, 2008 ) for peak detection, quantification and peptide identification using a full-length Uniprot Mouse database (April, 2016) and common contaminants such as keratins and enzymes used for digestion as reference. Carbamidomethylcysteine was set as fixed modification and protein amino-terminal acetylation, serine-, threonine- and tyrosine-phosphorylation, and oxidation of methionine were set as variable modifications. The MS/MS tolerance was set to 20 ppm and three missed cleavages were allowed using trypsin/P as enzyme specificity. Peptide, site and protein FDR based on a forwards-reverse database were set to 0.01, minimum peptide length was set to 7, and minimum number of peptides for identification of proteins was set to one, which must be unique. The ‘match-between-run’ option was used with a time window of 1 min. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD012068 (project name: Cyclin dependent kinase 5 (CDK5) regulates the circadian clock; project accession: PXD012068).
We raised in mouse a specific monoclonal antibody recognizing the phosphorylated form of serine 394 of PER2 in collaboration with GenScript Company. The sequence used for the immunogen preparation was: FDY {pSer} PIRFRTRNGEC. 3 Balb/c mice and 3 C57 mice were immunized with conventional strategies and antisera obtained from those animals were used for the first control experiment performed by in vitro kinase assay ( Figure 5—figure supplement 3 ). The positive antiserum was used for the cell fusions. Subsequently, a screening with 16 96-well plates (from 2 × 10E4 clones) was performed by indirect ELISA, primary screening with phospho-peptide, then counter-screening with non-phospho-peptide. The obtained supernatants were tested by in vitro kinase assay in order to screen which one was better recognized the phospho-form of PER2 S394 ( Figure 5—figure supplement 4 ). Finally, five selected positive primary clones selected were subcloned by limiting dilution and tested as final antibody ( Figure 5—figure supplement 5 ).
Statistical analysis of all experiments was performed using GraphPad Prism6 software. Depending on the type of data, either an unpaired t-test, or one- or two-way ANOVA with Bonferroni or Tukey’s post-hoc test was performed. Values considered significantly different are highlighted. [p<0.05 (*), p<0.01 (**), or p<0.001 (***)].
Data supporting the findings of this work are available within the paper and its supplementary files, and on Dryad ( https://doi.org/10.5061/dryad.4067r78 ). Non-commercial biological materials are provided upon request to the corresponding author. Proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD012068. The Per2Brdm1 mutant mouse strain is available at the Jackson Laboratory Stock No: 003819 (B6.Cg-Per2 tm1Brd Tyr c-Brd). The floxed Per2 mice are available at the European Mouse Mutant Archive (EMMA) strain ID EM:10599, B6;129P2-Per2tm1Ual/Biat.
Contribution, competing interests.
For correspondence, fondazione cenci bolognetti, instituto pasteur, schweizerischer nationalfonds zur förderung der wissenschaftlichen forschung (31003a_166682), schweizerischer nationalfonds zur förderung der wissenschaftlichen forschung (310030_166474/1), schweizerischer nationalfonds zur förderung der wissenschaftlichen forschung (316030_177088).
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Stéphanie Aebischer, Antoinette Hayoz, Cressida Harvey, Naila Ben Fredi, Jean-Charles Paterna (Viral Vector Facility, University of Zürich) and the Bioimage platform (University of Fribourg) for technical support. Funding from the Swiss National Science Foundation (31003A_166682, 310030_166474/1 and 316030_177088) is acknowledged. AB was supported by a fellowship from the Fondazione Cenci Bolognetti, Instituto Pasteur.
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Swiss Legislation and the declaration of Helsinki. The protocols were approved by the state veterinarian of the State of Fribourg (Permit Number: 2015-33).
© 2019, Brenna et al.
This article is distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use and redistribution provided that the original author and source are credited.
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Downloads (link to download the article as pdf).
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools), categories and tags.
Syntaxin-6 delays prion protein fibril formation and prolongs the presence of toxic aggregation intermediates.
Prions replicate via the autocatalytic conversion of cellular prion protein (PrP C ) into fibrillar assemblies of misfolded PrP. While this process has been extensively studied in vivo and in vitro, non-physiological reaction conditions of fibril formation in vitro have precluded the identification and mechanistic analysis of cellular proteins, which may alter PrP self-assembly and prion replication. Here, we have developed a fibril formation assay for recombinant murine and human PrP (23-231) under near-native conditions (NAA) to study the effect of cellular proteins, which may be risk factors or potential therapeutic targets in prion disease. Genetic screening suggests that variants that increase syntaxin-6 expression in the brain (gene: STX6) are risk factors for sporadic Creutzfeldt–Jakob disease. Analysis of the protein in NAA revealed, counterintuitively, that syntaxin-6 is a potent inhibitor of PrP fibril formation. It significantly delayed the lag phase of fibril formation at highly sub-stoichiometric molar ratios. However, when assessing toxicity of different aggregation time points to primary neurons, syntaxin-6 prolonged the presence of neurotoxic PrP species. Electron microscopy and super-resolution fluorescence microscopy revealed that, instead of highly ordered fibrils, in the presence of syntaxin-6 PrP formed less-ordered aggregates containing syntaxin-6. These data strongly suggest that the protein can directly alter the initial phase of PrP self-assembly and, uniquely, can act as an ‘anti-chaperone’, which promotes toxic aggregation intermediates by inhibiting fibril formation.
Proteasomes are essential molecular machines responsible for the degradation of proteins in eukaryotic cells. Altered proteasome activity has been linked to neurodegeneration, auto-immune disorders and cancer. Despite the relevance for human disease and drug development, no method currently exists to monitor proteasome composition and interactions in vivo in animal models. To fill this gap, we developed a strategy based on tagging of proteasomes with promiscuous biotin ligases and generated a new mouse model enabling the quantification of proteasome interactions by mass spectrometry. We show that biotin ligases can be incorporated in fully assembled proteasomes without negative impact on their activity. We demonstrate the utility of our method by identifying novel proteasome-interacting proteins, charting interactomes across mouse organs, and showing that proximity-labeling enables the identification of both endogenous and small-molecule-induced proteasome substrates.
Several metabolites have been shown to have independent and at times unexpected biological effects outside of their metabolic pathways. These include succinate, lactate, fumarate, and 2-hydroxyglutarate. 2-Hydroxybutyrate (2HB) is a byproduct of endogenous cysteine synthesis, produced during periods of cellular stress. 2HB rises acutely after exercise; it also rises during infection and is also chronically increased in a number of metabolic disorders. We show here that 2HB inhibits branched-chain aminotransferase enzymes, which in turn triggers a SIRT4-dependent shift in the compartmental abundance of protein ADP-ribosylation. The 2HB-induced decrease in nuclear protein ADP-ribosylation leads to a C/EBPβ-mediated transcriptional response in the branched-chain amino acid degradation pathway. This response to 2HB exposure leads to an improved oxidative capacity in vitro. We found that repeated injection with 2HB can replicate the improvement to oxidative capacity that occurs following exercise training. Together, we show that 2-HB regulates fundamental aspects of skeletal muscle metabolism.
An official website of the United States government
The .gov means it’s official. Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.
The site is secure. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.
Preview improvements coming to the PMC website in October 2024. Learn More or Try it out now .
Alison shupp.
1 Departments of Cancer Biology, Medical Oncology, Sidney Kimmel Cancer Center, Thomas Jefferson University, Philadelphia, PA, USA
2 Pennsylvania Cancer and Regenerative Medicine Research Center, Baruch S Blumberg Institute, Doylestown, PA, USA
Cyclin dependent kinases are proline-directed serine/threonine protein kinases that are traditionally activated upon association with a regulatory subunit. For most CDKs, activation by a cyclin occurs through association and phosphorylation of the CDK’s T-loop. CDK5 is unusual because it is not typically activated upon binding with a cyclin and does not require T-loop phosphorylation for activation, even though it has high amino acid sequence homology with other CDKs. While it was previously thought that CDK5 only interacted with p35 or p39 and their cleaved counterparts, Recent evidence suggests that CDK5 can interact with certain cylins, amongst other proteins, which modulate CDK5 activity levels. This review discusses recent findings of molecular interactions that regulate CDK5 activity and CDK5 associated pathways that are implicated in various diseases. Also covered herein is the growing body of evidence for CDK5 in contributing to the onset and progression of tumorigenesis.
Cyclin dependent kinases are proline-directed serine/threonine protein kinases that are traditionally activated upon association with a regulatory subunit. CDKs are a part of a kinase family that has been conserved throughout evolution and can be found in species from Saccharomyces cerevisia to humans. In humans there are 13 different CDKs (CDK1 - CDK13) that are highly expressed in mitotic cells [ 1 ]. For most CDKs, activation by a cyclin occurs through association and phosphorylation of the CDK's T-loop. Despite having high amino acid sequence homology with other CDKs, CDK5 is unusual because it is not typically activated upon binding with a cyclin and does not require T-loop phosphorylation for activation. Additionally, CDK5 has functions in both terminally differentiated and proliferating cells [ 2 ]. CDK5 was first identified in 1992 by multiple groups and was given a different name by each, including tau kinase II [ 3 ], neuronal Cdc2 like kinase [ 4 ], brain proline-directed kinase [ 5 ], PSSALRE [ 6 ], and CDK5 [ 7 ]. An isoform of CDK5, termed either CDK5-SV or CDK5-V1, was recently discovered [ 8 , 9 ]. One study reported that this splice variant lacks 32 amino acids encoded by exon 7 [ 8 ], while another study stated the missing 32 amino acids are encoded by exon 6 [ 9 ]. Although these two groups reported conflicting data, it has been suggested that the identified isoforms are in fact the same protein and the variances in their data are due to different methodologies [ 10 ].
CDK5 can be mapped to chromosome 7q36 and its expression is upregulated by the transcription factors Fos and CREB through the MEK/ERK pathway and by δFosB [ 11 , 12 ]. CDK5 plays a vital role in the central nervous system but has functions in other cell types. Outside of the nervous system, active CDK5 has been found in pancreatic β cells [ 13 ], corneal epithelial cells [ 14 ] and monocytes [ 15 ] amongst various other cell types [ 10 , 16 ]. In the nervous system, CDK5 is involved in neuron migration, neurite outgrowth and support, and synaptogenesis. CDK5's function in cells other than neurons includes the induction of cell motility, apoptosis, and cell cycle progression as well as functions involved with the immune system, lymphatic system, vascularization, and insulin secretion. A summary of CDK5 functions as discussed herein can be found in Table TableI. I . CDK5 has recently been implicated in diseases, including the development and progression of cancer and neurodegenerative diseases. For this reason, the regulation of CDK5 activity is now emerging as a candidate therapeutic target.
Biological system/process | CDK5 function | Mechanism |
---|---|---|
Central nervous system | Support growth cones | CDK5 phosphorylates CRMP2A at Ser27 during semaphorin3A stimulation. CDK5 also phosphorylates neurofilament heavy chain to promote neurofilament assembly [ – ] |
Growth cone collapse | CDK5 associates with alpha2-chimerin and phosphorylates CRMP2 at Ser522. CRMP2 further phosphorylated and inactivated by GSK3beta [ ] | |
Immune system | Increased IFNγ-induced PD-L1 expression | CDK5 expression decreases the expression of PD-L1 transcriptional repressors (IRF2 and IRF2BP) [ ] |
Insulin secretion | Reduction of insulin secretion | CDK5 phosphorylates L-VDCC and prevents exocytosis of insulin [ ] |
Vascular | Promotes angiogenesis | CDK5 expression increases abundance of HIF-1α [ ] |
Lymphatic | Lymphatic valve formation | CDK5 phosphorylates Foxc2, which regulates the expression of connexin 37 [ ] |
Cell Cycle | Increased expression of cyclins and other CDK's | Rb is a downstream target of CDK5's activity [ ] |
Reduction of CDK5 activity | Cyclin D1 and cyclin E can bind CDK5 to prevent CDK5's activation [ , ] | |
Cancer Progression | Cell proliferation | Reduction of p25 expression or CDK5 expression can prevent proliferation [ ] |
Cell migration/metastasis | CDK5 activity leads to caldesmon phosphorylation and actin polymerization. CDK5 enhances pro-migratory P13K/AKT signaling [ , ] |
Unlike other CDKs, CDK5 is not primarily activated by cyclins. Instead it is through specific binding with the proteins p35 or p39, or their respective cleaved counterparts p25 and p29, that CDK5 becomes active [ 1 , 17 , 18 ]. It was found that p35 knockout mice have defective cortical lamination and adults suffered from sporadic lethality and seizures [ 19 ], which is a less severe phenotype than that exhibited by Cdk5 knockout mice [ 20 ]. p39 -/- mice did not display any obvious abnormalities, however p35/p39 compound knockout mice displayed a phenotype identical to that of the Cdk5 -/- mice [ 21 ], suggesting that while p39 may not play a pivotal role in Cdk5 activation, it becomes necessary for nervous system development in the absence of p35.
p35 has a myristolation sequence that localizes it to phospholipid membranes [ 22 ]. Active CDK5 can phosphorylate p35 at Ser8 and Thr138. In the brain, phosphorylation of S8 is constant throughout development, but phosphorylation of T138 is found more abundantly in fetal brain tissue [ 23 ]. The phosphorylation at S8 leads to a more diffuse localization throughout the cytoplasm. This could be due to increased p35 mobility on membranes due to an altered interaction between the protein and phospholipids that constitute cell membranes [ 24 ]. p35 phosphorylation at T138 prevents its cleavage to p25 by calpain [ 23 ]. Because CDK5 has various regulatory functions in neuron development and migration, it is likely that the phosphorylation of p35 at T138 protects against aberrant CDK5 activation through formation of p25 in the fetal stage of brain development when CDK5 activity is also high [ 1 ]. Additionally, in vitro , under conditions of oxidative stress, p35 has been found conjugated to SUMO2 at Lys246 and Lys290, which led to increased p35/CDK5 activity [ 25 ].
As previously mentioned, the CDK5 activator p25 is formed through cleavage of p35 by calpain. This produces both the p25 product as well as a p10 product. Cleavage of p35 occurs under stress conditions such as amyloidβ presence, excitotoxicity, or oxidative stress [ 22 , 26 ]. This cleavage allows p25 to localize to nuclear and perinuclear regions by removing the p10 myristolation sequence [ 22 ]. Compared with p35, p25 has a longer half-life, and therefore prolongs the activation period of CDK5, leading to increased phosphorylation of CDK5's target proteins [ 22 , 27 ].
The functions of CDK5 activators p39 and p29 largely overlap with those of p35 and p25, respectively, however their expression throughout brain regions vary. p39 and p29 are mainly expressed in postnatal cerebral cortex and the hindbrain while p35 and p25 are largely expressed in the cerebral cortex of developing brains [ 27 ]. The localization of p39 to membranes is similar to that of p35 due to its conserved myristolation sequence [ 22 ]. Likewise, p39 also shows a more diffuse localization upon phosphorylation of Ser8 by CDK5 [ 24 ]. p39 can be phosphorylated by CDK5 at Ser173, a site equivalent to T138 in p35, and Thr84, however the effect of these phosphorylations on controlling protein stability have not yet been explored [ 1 , 24 ].
In addition to p35 and p39, cyclin I has also been shown to activate CDK5. Cyclin I-CDK5 binding targets CDK5 to the nucleus [ 28 ] and increases levels of anti-apoptotic proteins Bcl2 and Bcl2l1 via the MEK/ERK pathway [ 29 ] . This upregulation of Bcl2 and Bcl211 is observed only through cyclin I activation of CDK5, not activation via p35 [ 29 , 30 ]. CDK5 has been found to bind cyclin D1 and cyclin D3 in human fibroblasts, however this interaction had no influence on the activation and kinase activity of CDK5 [ 7 , 31 ].
While CDK5 is only activated by p35/p25, p39/p29, or cyclin I, the activity of CDK5 can be modulated by a variety of other proteins, as depicted in Figure Figure1. 1 . For instance, cyclin D1 can attenuate CDK5 kinase activity by competing with p35 for binding with CDK5, thereby forming an inactive complex of cyclin D1 and CDK5 (Fig. (Fig.1). 1 ). CDK5 and cyclin D1 can be found in the rat cerebellum during the first 24 days of postnatal development, albeit at varying abundances. CDK5 abundance increased while cyclin D1 decreased from day 9 on to adulthood [ 32 ]. In post-mitotic neurons, cyclin D1/CDK5 association was found to lead to cell cycle related neuronal apoptosis through sustained MEK/ERK signaling [ 33 ].
Cyclin E can directly interact with Cdk5 to reduce its activity. Cyclin E was found to sequester mouse Cdk5 away from other protein activators along with p27 KIP1 . The formation of this complex, and consequent attenuation of Cdk5 activity was found to promote synaptic plasticity, memory formation, and dendritic growth, as cyclin E -/- mice, that had increased Cdk5 activity, were deficient in these processes [ 34 ]. While this result may seem counterintuitive due to active CDK5's function with supporting neurite outgrowths, this observation could be explained by an overabundance of active CDK5 detrimentally effecting neurite outgrowth and subsequently synaptic plasticity. This theory would be consistent with findings that CDK5 expression levels are increased in certain neurodegenerative diseases, and that it is the aberrant CDK5 activity that leads to neurite collapse and death [ 35 – 37 ].
Glutathione-S-transferase (GSTP1) is another regulator of CDK5 activity that functions by competing with p35 for CDK5 binding. GSTP1 also reduces aberrant CDK5 activity by scavenging for molecules associated with oxidative stress and thereby decreasing the likelihood of p35/p39 cleavage to p25/p29 [ 38 ] (Figure (Figure1 1 ).
TP53 induced glycolysis regulatory phosphatase (TIGAR) has been shown to upregulate CDK5 expression levels in the presence of induced DNA damage (Figure (Figure1). 1 ). Knockdown of TIGAR led to decreased CDK5 expression, decreased phosphorylated ATM, and consequently increased levels of induced DNA damage. This suggests that DNA damage repair is mediated via TIGAR activation of the CDK5-ATM pathway [ 39 ].
Previously, CDK5 was thought to function in a cell cycle independent manner; however, recently the retinoblastoma protein (Rb) was discovered as a downstream target of CDK5. Expression of CDK5 leads to the phosphorylation of Rb, ultimately leading to the expression of cyclins and other cdks [ 40 ]. The protein kinase CK1 is phosphorylated by CDK5, and is involved in a wide array of signaling pathways including cell cycle, DNA repair, and apoptosis [ 41 ]. When CDK5 phosphorylates CK1, its kinase activity is subsequently reduced [ 42 ]. The functional affect of CDK5-mediated phosphorylation of CK1 on cell cycle, DNA repair, or apoptosis has yet to be explored.
In pancreatic β cells, CDK5 activity reduces insulin secretion in response to glucose abundance (Figure (Figure2). 2 ). This was demonstrated using CDK5 inhibitors, as well as inhibition of CDK5's activator p35. When CDK5 is active, it phosphorylates the L-type voltage-dependent Ca +2 channel (L-VDCC) at Ser783, which prevents the association of L-VDCC with syntaxin and SNAP-25, thereby preventing exocytosis of insulin from the cell [ 13 ].
Within the immune system, CDK5 has been implicated in IFNγ-induced programmed death ligand 1 (PD-L1) upregulation, which allows certain cells to evade detection by the immune system. Decreased CDK5 expression led to increased expression of the PD-L1 transcriptional repressors IRF2 and IRF2BP and consequent decreased PD-L1 expression (Figure (Figure2) 2 ) [ 43 ]. PD-L1 is a ligand that binds with PD-1, which is found on various immune cells. The binding of PD-L1 and PD-1 decreases an immune response by inhibiting T-cell activation and cytokine production. In normal tissues this is vital for maintaining homeostasis [ 44 ]. However tumor cells can also express PD-L1, which allows them to avoid detection and elimination by T-cells [ 45 , 46 ].
CDK5 promotes the formation of lymphatic vessels. CDK5 phosphorylates Foxc2, a protein that regulates the expression of connexin 37, which is critical for lymphatic valve formation (Figure (Figure2). 2 ). Moreover, knockout of CDK5 in the endothelium leads to lymphedema formation and embryonic lethality in mice [ 47 ].
CDK5 has previously been implicated in the migration of neurons. CDK5 knockout mice have abnormal cortical lamination, and more than 60% of CDK5 -/- mice died in utero [ 20 ]. Various studies have since implicated CDK5 in cell migration as it governs cancer metastasis. In prostate cancer cells, inhibition of CDK5 by the drug roscovitine prevented cell migration. The roscovitine treated cells did not project lamellopodia, and had reduced tubulin structures compared to untreated cells. This suggests that CDK5 inhibition prevented the establishment of cell polarity required for movement [ 48 ]. Additionally, knockdown of CDK5 in melanoma cell lines decreased cell motility and cell spreading in vitro , and decreased formation of lung and liver metastases in vivo in a mouse model of human melanoma. The decrease in CDK5 expression led to decreased phosphorylation of caldesmon, which decreased its binding affinity with actin and calmodulin (Figure (Figure2) 2 ) [ 49 ]. Another mechanism by which CDK5 may promote cell migration is by enhancing pro-migratory P13K/AKT signaling. CDK5 phosphorylates the Gα –interacting vesicle associated protein (GIV), which promotes GIV interaction with Gαi, thereby enhancing P13K/AKT signaling (Figure (Figure2) 2 ) [ 50 ]. Together, these studies demonstrate the importance of CDK5 in cell motility, a naturally occurring and necessary process. However, CDK5 mediated movement could also be an underlying driver of cancer metastasis and could be targeted in treatments to halt cancer metastasis.
An important function of CDK5, especially in neurons, is the organization of the cytoskeleton and support of cellular outgrowths (Figure (Figure2). 2 ). Expression of p35 or p39 in vitro stimulates neurite outgrowths, and a dominant negative mutant of CDK5 was found to abolish the formation of these outgrowths [ 51 ]. CDK5 supports axon and neurite outgrowth is through phosphorylation of the neurofilament heavy chain, resulting in the assembly of neurofilaments [ 52 ].
CDK5 has been shown to both prevent and promote growth cone collapse under different circumstances. CDK5 phosphorylates the protein CRMP2A at Ser27, which can be stabilized by Pin1 to support the growth of growth cones in the presence of semaphorin3A stimulation [ 53 , 54 ]. Additionally, CDK5 can promote axonal growth through indirect activation of CRMP2 by phosphorylating the protein Axin. Phosphorylated Axin inhibits GSK3β activity, leading to an increase in active, unphosphorylated CRMP2 [ 55 ] (Figure (Figure2 2 ).
Conversely, CDK5 promotes the collapse of growth cones through association with CRMP2 and α2-chimerin, an adaptor protein between CRMP2 and CDK5-p35. This association of CRMP2, α2-chimerin, and CDK5-p35 promotes the phosphorylation of CRMP2 at Ser522 by CDK5. In turn, this allows for CRMP2 to associate with and be phosphorylated at T514 by GSK3β, resulting in CRMP2 inactivation, microtubule disassembly, and ultimately growth cone collapse [ 56 ]. In this manner, CDK5 activity can both prevent and promote collapse of growth cones.
CDK5 can also reduce cellular outgrowth by regulating cytoskeletal organization through phosphorylation of p35 at T138, which prevents the polymerization of microtubules. This phosphorylation at T138 is found primarily in fetal brain tissues as opposed to adult brain [ 23 ].
Due to the many roles of CDK5 in the development of the nervous system, as well as the effects of cellular stress on CDK5 activation, CDK5 has been implicated in the progression of various neurological diseases and as a potential therapeutic target in disease treatment. For instance, while CDK5 normally phosphorylates collapsin response mediator protein 2 (CRMP2) to stimulate axon growth, it was found that hyperphosphorylation of CRMP2, as well as Tau, were implicated in the generation of neurofibrillary tangles characteristic of Alzheimer's disease [ 53 ]. Cell stress, including the presence of amyloid beta, is known to aberrantly activate CDK5 due to the formation of p25, which has been shown to cause the hyperphosphorylation of Tau, leading to atypical cell cycling, synaptotoxicity, and neuronal apoptosis [ 57 ]. Additionally, increased CDK5 activity caused by the sumoylation of p35 under oxidative stress, also contributes to neurodegeneration [ 25 ].
While CDK5 overexpression and aberrant activation are associated with neurodegenerative diseases, a loss or reduction in CDK5 activity is implicated in certain intellectual disabilities and neurodevelopmental disorders. Decreased CDK5 activity has been associated with intellectual disability in NF1 microdeletion syndrome patients [ 58 ] and schizophrenia [ 59 ]. Additionally, transgenic mice with decreased Cdk5 activity exhibited spontaneous seizures [ 60 ] as well as behaviors similar to ADHD [ 61 ].
Elevated levels of CDK5 have been found in various mouse tumors and human malignant tumors [ 40 ] [ 53 , 62 – 65 ]. The mechanisms involve effects on angiogenesis, cell proliferation and the immune system. As noted above, CDK5 enhances pRb phosphorylation and thereby cell-cycle progression [ 40 ]. Furthermore, CK1 is phosphorylated by CDK5, which in turn governs cell cycle, DNA repair, and apoptosis [ 41 ]. Increased levels of CDK5 target proteins are being considered as possible biomarkers of specific cancers. For example, an increase in CRMP2 phosphorylation could be a potential biomarker for certain lung cancers, as phosphorylated CRMP2 was found in the nuclei of biopsied lung cancer cells, but not cells in the surrounding epithelium [ 53 ].
In a transgenic mouse model of sporadic medullary thyroid carcinoma (MTC), p25 overexpression led to the development of bilateral malignant thyroid tumors, and was fatal after 30 weeks. However, arresting p25 expression at 5, 11, or 16 weeks led to 100 percent survival in all mice analyzed after 30 weeks. Similar results were discovered in vitro , in which reducing p25 expression or knocking down Cdk5 expression prevented further cell proliferation. This suggests that it is the aberrant activation of Cdk5 by p25 that leads to the progression of sporadic MTC [ 40 ].
CDK5 expression in medulloblastoma allows tumor cells to evade detection by T-cells in vivo . Conversely, decreased CDK5 expression enhanced the recruitment of CD4 + T-cells to the tumor site in mice, and increased the tumor-free survival rate of the mice. CDK5 regulates the evasion of tumors from the immune system by decreasing expression of transcriptional repressors of PD-L1 expression, thus increasing the abundance of PD-L1 [ 45 ].
Inhibiting CDK5 activity in hepatocellular carcinoma (HCC) cells prevented angiogenesis in vivo by decreasing the abundance of HIF-1α. Because HCC is a highly vascularized tumor type, inhibiting CDK5 and therefore angiogenesis, could prove a promising treatment for this tumor subtype and other highly vascularized tumors [ 64 ].
Due to the biological and clinically relevant importance of CDK5's function in multiple cell types, CDK5 presents an attractive therapeutic target for treating a variety of conditions such as diabetes, cancer, and neurodegeneration. Additionally, the upregulation of CDK5 associated with various cancers and neurodegenerative diseases further implicates its role in the development and progression of disease. Recently, tamoxifen (TMX), a drug currently used in breast cancer treatment, was found to decrease CDK5 activation by competitively binding with p35 and p25, and preventing their activation of CDK5. While the TMX inhibition of CDK5 activity could contribute to the anti-tumor effects of the drug, TMX treatment was also found to decrease Tau phosphorylation, suggesting a use for tamoxifen in treating Alzheimer's disease [ 66 ]. However, because of the broad functions of CDK5 in different cell and tissue types and the pan CDK inhibitory effect on other family members, the off target affects of a CDK5 inhibitory drug may create undesirable side effects. Nonetheless, CDK inhibitors are an intriguing clinical therapy for the treatment of various cancers. A list of current cyclin-dependent kinase inhibitors, including inhibitors of CDK5, and their associated clinical trials for the treatment of cancer can be seen in Table TableII II .
Treatment | Major Targets | Disease(s) | Clinical trial identifier |
---|---|---|---|
Terameprocol | CDK1 | Phase I: Leukemia, refractory solid tumors, lymphoma, glioma | NCT00664677, NCT00664586, NCT00404248 |
PHA-793887 | CDK1, CDK2, CDK4 | Phase I: Solid tumors | NCT00996255 |
Flavopiridol | CDK1, CDK2, CDK4, CDK7, CDK9 | Phase I-II: Various cancer including leukemia, multiple myeloma, lymphoma, sarcoma, and solid tumors (alone and in combination with other cytotoxic drugs) | NCT02520011, NCT00112723, NCT00005974, NCT00098579, NCT00007917, NCT00324480 |
BAY1000394 | CDK1, CDK2, CDK4, CDK9 | Phase I: solid tumors | NCT01188252 |
Dinaciclib | CDK1, CDK2, CDK5, CDK9 | Phase I-II: Advanced malignancies and relapsed multiple myeloma (alone and in combination with other cytotoxic drugs) | NCT01783171, NCT01624441, NCT01096342, NCT02684617, NCT01434316, NCT00871663, NCT01624441 |
P276-00 | CDK1, CDK4, CDK9 | Phase I-II: Multiple myeloma, mantle cell lymphoma, head and neck cancers, cyclin D1-positive melanoma | NCT00882063, NCT00848050, NCT00824343, NCT00899054, NCT00835419 |
AT7519 | CDK2, CDK4, CDK5, CDK9 | Phase I: Advanced or metastatic solid tumors, lymphoma | NCT02503709, NCT01652144, NCT01627054 |
R-roscovitine | CDK2, CDK5 | Phase I-II: Advanced solid tumors, non-small cell lung cancer | NCT00999401, NCT00372073 |
SNS-032 | CDK2, CDK7, CDK9 | Phase I: B-lymphoid malignancies and advanced solid tumors | NCT00446342 |
P1446A-05 | CDK4 | Phase I: Advanced refractory solid tumors and hematological tumors | NCT00840190 NCT00772876 |
PD 0332991 | CDK4, CDK6 | Phase I: Advanced cancers, mantle cell lymphoma Phase II: Multiple myeloma, advanced breast cancer, non-small cell lung cancer, ovarian cancer | NCT01522989, NCT00141297, NCT02008734, NCT02101034, NCT01976169, NCT01907607, NCT01356628, NCT01291017, NCT01536743 |
LY2835219 | CDK4, CDK6 | Phase I-II: Metastatic breast cancer, non small cell lung cancer | NCT02102490, NCT02246621, NCT02441946, NCT02450539, NCT02079636, NCT02779751 NCT02152631, NCT02675231 |
This is representative rather than a comprehensive list of past and present clinical trials in the field.
One of the most well studied CDK inhibitors being used in cancer clinical trials is flavopiridol, a drug developed by Tolero pharmaceuticals under the name Alvocidib. Flavopiridol was found to competivively bind to the ATP-binding pocket of CDK1, CDK2, CDK4, and CDK9, consequently inducing apoptosis in both dividing and quiescent cells. Early clinical trials with flavopiridol as a monotherapy proved ineffective in that there was a narrow window between no clinical response and severe, lethal tumor lysis. Ongoing trials involve combination therapies with other novel chemotherapy agents to overcome the limitations of flavorpiridol [ 67 ].
Another relatively well studied CDK inhibitor, Dinaciclib, was found to be more efficacious than flavopiridol, with IC 50 values in the low nanomolar range (1-4 nM – in various models flavopiridol's IC50 values range from 50-350 nM) [ 67 , 68 ]. Dinaciclib selectively inhibits CDK1, CDK2, CDK5, and CDK9 [ 67 ]. Preclinical studies and early clinical trials demonstrated the cytotoxicity of Dinaciclib in solid tumors and chronic lymphocytic leukemia, while not affecting T-cell function or number [ 69 ].
Roscovitine, marketed under the name Seliciclib, is an inhibitor of CDK5 and CDK2. Many of the clinical trials for Seliciclib were intiated determine dose-limiting toxicities of the drug alone or in combination with other chemotherapeutics. While roscovitine is used widely experimentally to inhibit CDK5 activity, it is not being intensively examined as a clinical cancer therapeutic [ 67 ].
To potentially reduce broad undesirable off target effects of pan-CDK inhibitors, CDK5 inhibitory peptide (CIP) has been studied as a potential therapeutic for neurodegeneration. CIP specifically targets the hyperactivated state of CDK5 as mediated by p25/p29, while allowing normal activation of CDK5 by p35/p39. CDK5 inhibitory peptide (CIP) was found to inhibit the hyperactivation of CDK5 by p25 overexpression in vivo , which reduced neurodegeneration and improved cognitive function of transgenic mice, without affecting neurodevelopment [ 70 ]. In the future, CIP could possibly be adapted to treat certain cancers caused by aberrant CDK5 activation.
CONFLICTS OF INTEREST
There is no conflict of interest.
Loading metrics
Open Access
Peer-reviewed
Research Article
Contributed equally to this work with: Bitna Lim, Yurika Matsui
Roles Data curation, Formal analysis, Investigation, Visualization
Current address: CHA Future Medicine Research Institute, Seongnam, Republic of Korea
Affiliation Department of Developmental Neurobiology, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Roles Data curation, Formal analysis, Methodology, Visualization, Writing – review & editing
Roles Formal analysis
Affiliation Center for Applied Bioinformatics, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Affiliation Center for Proteomics and Metabolomics, St. Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Roles Data curation, Formal analysis
Roles Supervision
Affiliation Department of Biostatistics, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Roles Methodology
Affiliation Department of Cell & Molecular Biology, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Roles Data curation
Affiliations Department of Developmental Neurobiology, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America, Department of Structural Biology, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
Roles Formal analysis, Supervision
Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing
* E-mail: [email protected]
53BP1 is a well-established DNA damage repair factor that has recently emerged to critically regulate gene expression for tumor suppression and neural development. However, its precise function and regulatory mechanisms remain unclear. Here, we showed that phosphorylation of 53BP1 at serine 25 by ATM is required for neural progenitor cell proliferation and neuronal differentiation in cortical brain organoids. Dynamic phosphorylation of 53BP1-serine 25 controls 53BP1 target genes governing neuronal differentiation and function, cellular response to stress, and apoptosis. Mechanistically, ATM and RNF168 govern 53BP1’s binding to gene loci to directly affect gene regulation, especially at genes for neuronal differentiation and maturation. 53BP1 serine 25 phosphorylation effectively impedes its binding to bivalent or H3K27me3-occupied promoters, especially at genes regulating H3K4 methylation, neuronal functions, and cell proliferation. Beyond 53BP1, ATM-dependent phosphorylation displays wide-ranging effects, regulating factors in neuronal differentiation, cytoskeleton, p53 regulation, as well as key signaling pathways such as ATM, BDNF, and WNT during cortical organoid differentiation. Together, our data suggest that the interplay between 53BP1 and ATM orchestrates essential genetic programs for cell morphogenesis, tissue organization, and developmental pathways crucial for human cortical development.
Citation: Lim B, Matsui Y, Jung S, Djekidel MN, Qi W, Yuan Z-F, et al. (2024) Phosphorylation of the DNA damage repair factor 53BP1 by ATM kinase controls neurodevelopmental programs in cortical brain organoids. PLoS Biol 22(9): e3002760. https://doi.org/10.1371/journal.pbio.3002760
Received: August 25, 2023; Accepted: July 19, 2024; Published: September 3, 2024
Copyright: © 2024 Lim et al. This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All sequencing data are deposited in NCBI GEO database under accession number GSE231321. Codes for analyzing sequencing data are deposited in https://doi.org/10.6084/m9.figshare.7411835 . Mass spectrometry data were deposited in ProteomXchange, with project accession number PXD041699. Numerical data are in S1 Data , and uncropped Western blot images are in S1 Raw Images .
Funding: This work was supported by the American Lebanese Syrian Associated Charities ( https://www.stjude.org/ to JCP), American Cancer Society ( https://www.cancer.org/ ; 132096-RSG-18-032-01-DDC to JCP), and NIH ( https://www.nih.gov/ ; 1R01GM134358-05 to JCP). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: ATM, ataxia telangiectasia mutated; BSA, bovine serum albumin; FDR, false discovery rate; GO, Gene Ontology; GSEA, gene set enrichment analysis; hESC, human embryonic stem cell; KO, knockout; NPC, neural progenitor cell; PSM, peptide spectral match; TMT LC-MS/MS, liquid chromatography-tandem mass spectrometry; WB, western blot; WT, wild type; 53BP1, p53 binding protein 1; 53BP1-pS25, 53BP1 phosphorylated at serine 25
Transcription ensures the proper expression of genetic information for the development and function of the organism, whereas DNA repair maintains the integrity of the genetic code. These 2 processes share cross-functional factors, including CSB, TFII, and XPG, which repair DNA damage caused by torsional stress from transcription-initiating RNA polymerase II [ 1 – 3 ]. Conversely, some proteins initially believed to function exclusively in DNA repair have been found to regulate gene expression. For example, 53BP1 (p53 binding protein 1) is a key regulator of DNA repair mechanisms, promoting nonhomologous end-joining over homologous recombination [ 4 ]. During the DNA damage response, 53BP1 plays a pivotal role in p53-mediated activation of tumor suppressive genetic programs [ 5 ]. Recent research has also revealed that 53BP1 collaborates with the chromatin modifier UTX in neural progenitor cells (NPCs), promoting an open chromatin to facilitate the activation of neurogenic or corticogenic programs [ 6 ]. Intriguingly, the 53BP1–UTX interaction is observed in humans but not in mice [ 6 ]; the mechanism is not well conserved and regulates primate neurodevelopment. These discoveries highlight the importance of 53BP1 in gene regulation for tumor suppression and neural development. However, the precise mechanisms underlying 53BP1’s role in gene regulation and its upstream mechanism are yet to be fully understood.
Studies of 53BP1 have primarily focused on its role in the DNA damage response. To localize to chromatin with double-stranded breaks, 53BP1 uses its BRCT domain to bind to γH2AX, the Tudor domain to bind to H4K20 dimethylation, and its UDR segment to bind to ubiquitinated H2AK15 [ 7 – 10 ]. Additionally, the phosphorylated SQ/TQ motif of 53BP1 coordinates the docking of RIF1 or SCAI, selectively promoting nonhomologous end-joining or reducing homologous recombination [ 11 , 12 ]. These interactions are likely relevant to the gene regulatory activities of 53BP1. For example, γH2AX recruits 53BP1 and is required for resolving R-loops, DNA demethylation, transcription activation, and transcription elongation [ 13 , 14 ]. These findings suggest that the activities of 53BP1 in DNA damage response are interconnected with its gene regulatory functions.
The studies mentioned above have contributed to a model of 53BP1, where posttranslational modifications of its different residues and domains coordinate various activities. Most prominently, numerous residues of 53BP1 are phosphorylated by ATM (ataxia telangiectasia mutated) kinase [ 10 , 15 – 17 ]. ATM-mediated phosphorylation of 53BP1 or 53BP1-interacting proteins controls protein interactions, cellular localization, and DNA repair mechanisms [ 11 , 12 , 18 – 20 ]. Despite these discoveries, the impact of phosphorylation on the gene regulatory activity of 53BP1 remains unknown. Here, we report that phosphorylation of 53BP1-serine 25 by ATM is crucial for the proper expression of genetic programs during the growth and development of cortical brain organoids. ATM-dependent phosphorylation controls the chromatin binding of 53BP1 to genomic targets functioning in several key pathways, including neuronal differentiation, cytoskeleton, p53, and ATM, BNDF, and WNT signaling pathways. These results highlight the essential role of 53BP1 phosphorylation in regulating genetic programs for the differentiation of cortical brain organoids.
Although 53BP1 is required for human embryonic stem cells (hESCs) to differentiate into NPCs [ 6 ], its levels did not change during differentiation ( Fig 1A ). Therefore, its regulation is likely posttranslational during neural differentiation. Human NPCs were analyzed by RNA-seq and immunofluorescence to validate successful NPC generation ( S1A-S1E Fig ). 53BP1 is targeted by various kinases, including ATM, and we hypothesized that 53BP1 phosphorylation regulates the differentiation of hESCs into NPCs. Intriguingly, we found that the levels of 53BP1 phosphorylated at serine 25 (53BP1-pS25) were markedly increased in NPCs compared to hESCs (Figs 1A and S1F-S1H ). The levels of the DNA damage marker γH2AX were similar in NPCs and hESCs ( S1I Fig ), suggesting that the increase in 53BP1-pS25 levels during NPC differentiation is not due to increased DNA damage.
( A ) WB of the nuclear extract of hESCs and hNPCs showed marked increase of 53BP1-pS25 in hNPCs. WB analysis of IgG, ( B ) 53BP1, and ( C ) ATM co-immunoprecipitation in the nuclear extract of hESCs. ( D ) Quantification of the relative ATM protein levels (normalized to β-ACTIN) in 5 replicate WB analyses of hESCs and hNPCs. ( E ) Schematic diagram of the cortical organoid differentiation. Aggregates were formed in the induction media for 17 days, embedded in Matrigel droplets and cultured in cortical differentiation medium for 16 days, and then cultured in cortical maturation media thereafter. ( F ) WB analysis of WT and ATM-KO cortical organoids at day 35 of differentiation. Immunofluorescence of ( G ) PAX6 and CTIP2 and ( J ) KI67 in cryosections of cortical organoids at day 35 of differentiation. Bar, 100 μm. At day 35 of differentiation, the ( H ) area and ( I ) thickness of VZ-like regions were compared between groups. Data points represent single organoids. The mean ± SEM values were compared by one-way ANOVA with Dunnett’s multiple comparisons test to yield **** indicating p < 0.0001. n = 13 organoids/group. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; DMEM, Dulbecco’s Modified Eagle Medium; GMEM, Glasgow Modified Essential Medium; hESC, human embryonic stem cell; hNPC, human neural progenitor cell; IgG, immunoglobulin G; KO, knockout; KSR, Knockout Serum Replacement; VZ, ventricular zone; WB, western blot; WT, wild type; 53BP1, p53 binding protein 1; 53BP1-pS25, 53BP1 phosphorylated at serine 25.
https://doi.org/10.1371/journal.pbio.3002760.g001
The ATM kinase phosphorylates 53BP1-S25 [ 15 ], and we thus investigated whether ATM plays a role in neural differentiation. First, we found that ATM co-immunoprecipitated with 53BP1, as did the positive control UTX, but not with the negative control SUZ12 (a core subunit of PRC2, which does not bind these proteins; Fig 1B ). Similarly, 53BP1 co-immunoprecipitated with ATM, but not with the negative control SUZ12 ( Fig 1C ). Like 53BP1-pS25, ATM levels were significantly increased in NPCs compared with hESCs ( Fig 1D ). ATM up-regulation in NPCs was shown by a previous DNA damage response study [ 21 ].
Next, we used the CRISPR-Cas9 system to generate 4 ATM -knockout (KO) hESC lines (Figs 2A , 2B , and S1J ). Data from RNA-seq and immunofluorescence showed that ATM -KO did not markedly alter hESC pluripotency ( S2C and S2D Fig ). All cell lines underwent karyotyping analysis and were characterized as karyotypically normal ( S1 Table ). ATM -KO lines had minor abnormalities, as expected due to the requirement of ATM for DNA damage repair. To analyze the role of ATM in human cortical development, we used an established protocol to differentiate wild type (WT) and ATM -KO hESCs into cortical organoids ( Fig 1E , Methods; [ 6 ]). We did not detect 53BP1-pS25 in ATM -KO D35 cortical organoids ( Fig 1F ) nor NPCs ( S2B Fig ), consistent with loss of ATM-mediated phosphorylation of 53BP1-S25 during neural differentiation of hESCs. ATM -KO modestly reduced γH2AX levels in NPCs ( S2E Fig ), suggesting that ATM promotes the phosphorylation of H2AX-S139 in NPCs.
Quantification of PAX6/CTIP2 ratios in ( A ) D28 and ( B ) D37 cortical organoids. ( C ) Quantification of NEUN/DAPI in D37 cortical organoids. ( D ) In D37 cortical organoids, 6 organoids were surveyed to count ZO-1-positive apical surfaces and proportions of PAX6-positive NPCs that are organized around the apical surfaces. ( E ) Proportions of PH3-positive cells that are adjacent to ZO-1-positive apical surfaces (“rings”). Data from 53BP1-S25A and S25D were included for comparison. **, p < 0.01; ***, p < 0.001; ****, p < 0.0001; ns, not significant by Welch’s t test in (A-C) and two-way ANOVA test in (D). From GSEA, functional terms that are highly enriched in ( F ) up-regulated and ( G ) down-regulated genes in ATM -KO vs. WT NPCs. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; GSEA, gene set enrichment analysis; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; WT, wild type; 53BP1, p53 binding protein 1.
https://doi.org/10.1371/journal.pbio.3002760.g002
By day 35 (D35) of differentiation, WT cortical organoids expressed the forebrain NPC marker PAX6 in ventricular zone–like regions that were radially organized ( Fig 1G ). In contrast, ATM -KO D35 cortical organoids displayed disorganized and smaller ventricular zone–like regions ( Fig 1G-1I ). We quantitatively compared NPC proliferation, neuronal differentiation, cell death, and cell organization in ATM -KO versus WT cortical organoids. Examination of endogenous DNA damage, by γH2AX immunofluorescence, did not reveal marked difference ( S3A Fig ), confirming our western blot (WB) results in S1I Fig . Although quantification of cell death marker cleaved-caspase 3 by FACS revealed a modest increase of cell death in D21 ATM-KO cortical organoids, FACS and immunofluorescence quantification showed that D28 and D35 ATM -KO and WT cortical organoids are similar in cell death frequencies ( S3B-S3E Fig ). Cell proliferation frequencies did not significantly differ between D28 and D35 ATM -KO and WT cortical organoids ( S4 Fig ). We next quantified immature neuronal marker NEUN and PAX6/CTIP2 ratios. Despite lower levels of immature neuronal differentiation, ATM-KO exhibited higher neuronal maturation (Figs 2A-2C and S5A ). These data suggest that ATM -KO fastens the phase of immature neuronal differentiation, leading to enhanced neuronal maturation. Finally, we quantified ZO-1-positive ventricular surfaces and the organization of PH3-positive and PAX6-positive cells around ventricular surfaces. The ATM -KO ventricular surfaces were similar to WT at D28 ( S5B-S5D Fig ), but the number was much reduced by D37 ( Fig 2D ). NPC organization around the ventricular surfaces were similarly organized in D37 (Figs 2D and S5E ); however, fewer ATM -KO proliferative cells were adjacent to ventricular surfaces ( Fig 2E ). These data suggest that ATM -KO enhances neuronal maturation and cellular disorganization in developing cortical organoids. By D55, ATM -KO cortical organoids had similar size distribution as the control ( S5F and S5G Fig ). Thus, ATM controls neuronal differentiation and cellular organization to form ventricular zone–like regions in cortical organoids.
To investigate the molecular basis of the cellular defects we observed in ATM -KO, we performed RNA-seq to compare ATM -KO to WT NPCs and D35 cortical organoids derived from WT and ATM -KO hESCs. The expression of forebrain markers was similar between WT and ATM -KO cortical organoids (and low expression of midbrain and hindbrain markers; S2 Table ), suggesting that the ATM -KO cortical organoids specified to the forebrain lineage. A false discovery rate (FDR) <0.05 was used to identify differentially expressed genes. Gene set enrichment analysis (GSEA) showed that up-regulated genes in ATM -KO NPCs were enriched in forebrain development, axis specification, and metabolic pathways (Figs 2F and S6A ), whereas down-regulated genes were enriched in neuronal differentiation, epithelial mesenchymal transition, and tube morphogenesis ( Fig 2G ). Comparison of transcriptomes profiles of D35 cortical organoids from 8 ATM -KO versus 6 WT datasets yielded similar GSEA terms ( Fig 3A and 3B ). These data suggest that ATM regulates genetic programs related to forebrain development, metabolism, and neuronal differentiation in NPCs and cortical organoids. Dysregulated genetic programs likely contributed to enhanced neuronal maturation and cellular disorganization in ATM -KO cortical organoids.
From GSEA, functional terms that are highly enriched in ( A ) up-regulated and ( B ) down-regulated genes in ATM -KO D35 cortical organoids. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. ( C ) Schematic diagram outlining TMT LC-MS/MS profiling of total proteomics and phosphoproteomics of D35 WT and ATM -KO cortical organoids. TMT signals from total proteomics were used to normalize those of phosphopeptides. ( D ) Using FC>1.5 and FDR<0.05, 198 phosphoproteins were found to be lower in 2 ATM -KO versus WT. ( E ) Normalized levels of phosphoproteins that have ATM-dependent phosphorylation in D35 cortical organoids. 53BP1 and EIF4EBP1 were known substrates of ATM. The error bars depict the mean and standard error of the mean values, which were calculated based on the normalized levels of each phosphopeptide in the protein. ( F ) Enrichment of proteins with ATM-dependent phosphorylation in specific functional categories. ( G ) Heatmap showing altered activities of kinases between D35 ATM -KO2 and WT cortical organoids. Relative changes in kinase activity are shown as row Z‐scores. Kinase activity was inferred by IKAP [ 25 ] based on normalized substrate phosphorylation levels from phosphor-proteome. The normalization was performed by dividing phosphor-peptide abundance of each protein by corresponding protein abundance [ 57 ]. Circos plots showing kinases with inferred ( H ) higher and ( I ) lower activities in D35 ATM -KO versus WT cortical organoids and their corresponding enriched pathways. Underlying numerical values for figures are found in S1_Data.xlsx. ATM, ataxia telangiectasia mutated; FC, fold-change; FDR, false discovery rate; GSEA, gene set enrichment analysis; KO, knockout; NES, normalized enrichment score; TMT LC-MS/MS, liquid chromatography-tandem mass spectrometry; WT, wild type; 53BP1, p53 binding protein 1.
https://doi.org/10.1371/journal.pbio.3002760.g003
As ATM kinase is crucial for many cell and developmental processes, we aimed to analyze its effect on the proteome and phosphoproteome of differentiating cortical organoids. First, we used multiplexed tandem mass tag-based quantification and 2D liquid chromatography-tandem mass spectrometry (TMT LC-MS/MS) to profile the proteome of WT and ATM -KO D35 cortical organoids ( S3B Fig , Methods). We quantified 10,895 proteins between 4 WT, 4 ATM -KO2, 3 ATM -KO3, and 3 ATM -KO14 D35 cortical organoid samples by using the criteria of fold change >1.5 and FDR <0.05 ( Fig 3C ). Consistency between replicate datasets is supported by principal component analysis ( S3B Fig ). GSEA showed that compared to WT, up-regulated proteins in ATM -KO were enriched in terms related to neurotransmission, neuron spine, dendrite, synapse, and axon ( S3C Fig ), whereas down-regulated proteins were enriched in BMP/TGFβ and WNT signaling, epithelial morphogenesis, and stem cell differentiation ( S3D Fig ). These data suggest that ATM controls posttranscriptional and translational gene regulation to suppress neuronal function and promote stem cell differentiation, epithelial morphogenesis, and TGFβ and WNT signaling pathways in D35 cortical organoids.
We have observed distinct patterns in the transcriptomics and proteomics data in ATM -KO versus WT cortical organoids. Interestingly, while transcriptomic programs related to neuronal differentiation were down-regulated ( Fig 3B ), proteomic programs related to neuronal function were up-regulated ( S3C Fig ) in ATM -KO versus WT cortical organoids. These differential patterns in transcriptomics and proteomics are likely a consequence of the regulatory role of ATM in multiple cellular processes. It is possible that the higher protein expression related to neuronal functions in ATM -KO lead to down-regulation of transcriptional expression of neuronal differentiation programs. This would suggest that the dysregulated transcriptomic and proteomic programs in ATM-KO cortical organoids are interconnected and result from the complex interplay of ATM’s regulation of various cellular pathways. These findings shed light on the intricate role of ATM in coordinating gene expression and protein levels, influencing neuronal differentiation and function in cortical organoids.
To investigate how ATM exerts its modulatory control during cortical organoid formation, we performed phosphoproteomics analysis of WT and ATM -KO cortical organoids. Using TMT LC-MS/MS, we quantified 22,646 phosphopeptides and normalized their abundance based on the protein abundance measured in the total proteomics analysis. A comparison between WT and ATM -KO lines revealed that 198 proteins had consistently lower levels of phosphorylation in at least 2 of the 3 ATM -KO lines (log2(fold change >1.5) and FDR <0.05; Fig 3D and S3 Table ). Among these proteins, 53BP1 and EIF4EBP1 were known substrates of ATM ( Fig 3E ) [ 15 , 16 , 22 , 23 ], validating the approach to identify putative ATM substrates in cortical organoids. However, it is essential to note that this approach does not distinguish between direct and indirect effects, and, therefore, some of the identified proteins could be phosphorylated by protein kinases that require ATM for their activity. Notably, many ATM-dependent phosphorylated proteins were found to be key neurodevelopmental regulators ( Fig 3E ) and enriched in functions related to neurodevelopment, neurogenesis, cell morphogenesis, and cytoskeleton ( Fig 3F ). These findings suggest that ATM plays a critical role in regulating the phosphorylation of proteins involved in essential processes for neurodevelopment and neuronal function in cortical organoids.
We further explored the effects of ATM by identifying protein kinases that had ATM-dependent phosphorylation. We used the IKAP machine learning algorithm [ 24 ] to analyze substrates (inferred from literature curation) and deduce the activities of those kinases. For example, in ATM -KO compared to WT, we found reduced phosphorylation of proteins related to MAPK9 activities, such as DCX, MAPT, and NFATC4 ( S7A Fig and S4 Table ) [ 24 ]. On the other hand, we found higher phosphorylation of proteins related to CDK5 activities, including ADD2, ADD3, DCX, DNM1L, DPYSL3, MAPT, and SRC ( S7B Fig and S4 Table ) [ 24 ]. In ATM -KO, we inferred lower activities in MAPK9, CDK2, CHEK1, ATR, CSNK1A1, MTOR, CAMK2A, and PRKACA (Figs 3G and S7C ), with enriched functions in ATM signaling, BNDF signaling, and axon guidance (Figs 3G , 3H , and S7C ). On the other hand, we inferred higher activities in GSK3B, MAPK3, PAK1, CSNK2A1, CDK5, CDK1, and PRKDC (Figs 3G , 3I , and S7C ), with enriched function in ATM signaling, WNT signaling, G2/M checkpoint, and p53 regulation in ATM -KO ( Fig 3H and 3I ). ATM KO leads to both lower and higher activities of kinases in ATM signaling. Additionally, some of the altered kinase activities could be secondary to ATM -KO, as CHEK1, ATR, and PRKDC were known substrates of ATM [ 23 , 25 ]. Overall, these data suggest that the activities of kinases related to ATM signaling, BNDF signaling, WNT signaling, G2/M checkpoint, and p53 regulation became dysregulated in ATM -KO D35 cortical organoids.
We thus conclude that ATM plays a crucial role in controlling key neurodevelopmental regulators. The dysregulated phosphorylation and activities of these regulators disrupt the normal transcriptomic program responsible for neuronal differentiation, leading to higher proteomic programs associated with neuronal function. As a consequence, the dysregulated programs in ATM -KO cortical organoids are likely responsible for the observed defects in neurogenesis and morphogenesis (formation of ventricular zone–like regions). These findings provide valuable insights into the role of ATM in neurodevelopment and shed light on potential molecular mechanisms underlying neurological disorders associated with ATM dysfunction.
We next examined ATM-dependent phosphorylation of 53BP1-S25. To specifically investigate the functional significance of 53BP1-pS25, we used the CRISPR-Cas9 system to mutate the endogenous 53BP1 serine 25 to alanine (S25A) or aspartic acid (S25D) (Figs 4A and S7D , Methods). The alanine substitution precludes phosphorylation, whereas aspartic acid is chemically similar to phosphoserine [ 26 ]. We generated 4 53BP1-S25A hESC lines (34–3, 34–4, 79–1, and 79–3) and 4 53BP1-S25D hESC lines (14–3, 14–15, 14–19, and 17). The total levels of 53BP1 were similar in WT, 53BP1-S25A, and 53BP1-S25D NPCs, and we did not detect pS25 in 53BP1-S25A NPCs, as expected ( S1D and S7E Figs). Control, 53BP1-S25A, and 53BP1-S25D hESC lines displayed similar transcriptomic profiles and pluripotency marker expression ( S2C , S7F , and S8A Figs), suggesting that 53BP1-S25A and 53BP1-S25D do not affect hESC self-renewal.
ATM is required for the phosphorylation of many neurodevelopmental regulators ( Fig 3F ). As the role of 53BP1-S25 beyond DNA damage repair is not known, we seek to analyze its role in human cortical development. We differentiated control WT, 53BP1-S25A, and 53BP1-S25D hESCs into cortical organoids ( Fig 1E , Methods). The 53BP1-S25A and 53BP1-S25D D35 cortical organoids displayed smaller sizes compared to WT controls ( Fig 4B and 4C ), suggesting that phosphorylation at S25 is essential for cortical organoid growth and development. Analysis of the ventricular zone–like regions showed 53BP1-S25A and 53BP1-S25D are significantly smaller than those in WT ( Fig 4D-4F ). Fewer cells were positive for KI67 (proliferation marker) or phosphorylated-serine 10 histone H3 (mitotic chromatin marker) ( S8B–S8D Fig ). Examination of endogenous DNA damage and cell death, assessed by γH2AX and cleaved-caspase 3, respectively, did not reveal significant differences between 53BP1-S25A, S25D, and WT ( S9 Fig ). To explore the developmental timing of the cellular phenotypes, we quantified KI67, NPC marker PAX6, and neuronal marker CTIP2 in D14, D21, D28, and D35 cortical organoids. At D28, 53BP1-S25A and S25D cortical organoids had significantly lower cell proliferation and higher neuronal differentiation ( Fig 5 ). Quantification of the tight junction protein ZO-1 showed significantly fewer ZO-1-positive ventricular surfaces in D28 53BP1-S25A and S25D cortical organoids compared to WT ( S10 Fig ). For the ventricles that did form in D28 53BP1-S25A and S25D, their surface areas did not significantly differ from those of WT ( S10B Fig ). These data suggest that lower ventricle formation, lower cell proliferation, and higher neuronal differentiation contributed to the depletion of progenitor pools and smaller cortical organoids in 53BP1-S25A and S25D.
( A ) In the endogenous 53BP1 locus, the codon TCT encoding serine-25 in was mutated to GCT and GAT encoding alanine and glutamate, respectively. ( B ) Bright-field images of cortical organoids formed by 4 53BP1-S25A lines, 4 53BP1-S25D lines, and 2 WT control at day 35 of differentiation. Bar, 1.5 mm. At day 35 of differentiation, the ( C ) organoid size and ( F ) area of ventricular zone–like region were compared between groups. Data points represent single organoids. The mean ± SEM values were compared by one-way ANOVA with Dunnett’s multiple comparisons test to yield ****, ***, **, *, and ns indicating p < 0.0001, 0.001, 0.01, 0.05, and not significant, respectively. n = 39–47 organoids/group for ( C ) and 15–33 organoids/group for ( F ). ( D ) Immunofluorescence of PAX6 and CTIP2 in cryosections of cortical organoids at day 35 of differentiation. Bar, 100 μm. ( E ) Illustration of ventricular zone–like areas in cortical organoids. Underlying numerical values for figures are found in S1 Data . WT, wild type; 53BP1, p53 binding protein 1; 53BP1-pS25, 53BP1 phosphorylated at serine 25.
https://doi.org/10.1371/journal.pbio.3002760.g004
FACS quantified ratios of ( A ) KI67, ( B ) PAX6, and ( C ) CTIP2 to total cells in D28 cortical organoids. ( D ) Immunofluorescence of PAX6 and CTIP2 in D28 cortical organoids. Bar, 100 μm. Quantification of immunofluorescence signals of ( E ) PAX6/DAPI, ( F ) CTIP2/DAPI, and ( G ) PAX6/CTIP2 in D28 cortical organoids. Each data point represents quantification of cells in 1 cortical organoid. Quantification of ( H ) KI67/PAX6 and ( I ) PAX6/CTIP2 ratios in immunofluorescence of D35 cortical organoids. Each data point represents quantification of cells in 1 cortical organoid. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001; ns, not significant by two-way ANOVA test. Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.g005
At D55, the 53BP1-S25A and S25D cortical organoids remained significantly smaller than WT ( S12 Fig and S5 and S6 Tables). These data suggest that the cell biological effects of the S25A and S25D mutations were similar, despite the aspartic acid mutation (S25D) being chemically similar to phosphoserine, which is the phosphorylated form of S25. The S25D mutation may act as an inhibitory mimic of phosphorylation, akin to the S25A mutation. Consequently, the absence of S25 phosphorylation impacts NPC proliferation and overall cortical organoid growth.
Using RNA-seq, we examined the transcriptomes of WT (6 samples), 53BP1-S25A (8 samples), and 53BP1-S25D (8 samples) D35 cortical organoids. We analyzed expressed genes with counts per million values >1 and observed few differences in gene expression between 53BP1-S25A and 53BP1-S25D D35 cortical organoids, using FDR <0.05 ( Fig 6A ). When comparing the transcriptomes of 53BP1-S25A and 53BP1-S25D organoids to WT, there were high concordant changes in gene expression, with over 87% of differentially expressed genes in 53BP1-S25A also being altered in 53BP1-S25D (Figs 6B and S12D ). However, 53BP1-S25D disrupted the expression of 2- to 3-fold more genes than 53BP1-S25A, suggesting a gain-of-function effect for the 53BP1-S25D mutation. To explore this further, we performed GSEA and found that the top terms enriched in the up-regulated genes of 53BP1-S25A and S25D organoids were highly overlapping ( Fig 6C ). In contrast, there was low overlap of the top terms in the down-regulated genes in 53BP1-S25A versus WT and 53BP1-S25D versus WT ( S12E Fig ). Both mutations led to the up-regulation of genes related to synapse, axon, and neurotransmitter functions, suggesting a shared effect on enhancing neuronal function ( Fig 6C ). The S25D mutation specifically up-regulated more genes involved in neuronal function compared to S25A, indicating a stronger impact on this aspect of gene regulation ( Fig 6C and 6D ). These findings highlight the significance of the S25 phosphorylation site in 53BP1 for the regulation of genes involved in neuronal function and support that the S25D mutation results in a gain-of-function effect, leading to more pronounced changes in gene expression related to neuronal processes.
( A ) Number of differentially expressed genes identified by pairwise comparisons at FDR <0.05. At day 35 of differentiation, 53BP1-S25A and S25D cortical organoids are molecularly similar. ( B ) Differentially expressed genes in 53BP1-S25D versus WT overlap 87% (764/875) and 91% (361/396) of those in 53BP1-S25A versus WT. ( C ) Extensive overlap of up-regulated GSEA terms between 53BP1-S25A versus WT and 53BP1-S25D versus WT. Most terms relate to axon, synapse, and neurotransmitter. ( D ) Of 53BP1 target genes up-regulated by S25A and S25D, 212 genes require WT 53BP1 for expression in cortical organoids. ( E ) The 212 genes are enriched in functions related to transcriptional regulation, neuron projection, axonogenesis, synapse, neurotransmitter synthesis and transport, and membrane depolarization. ( F ) Venn diagrams depict high overlaps between down-regulated genes in all 3 groups of mutant versus WT pairwise comparisons. ( G ) GSEA graphs showed that down-regulated genes in 53BP1-S25A or S25D vs. WT had significant enrichment in down-regulated genes of ATM-KO vs. WT cortical organoids. P values were calculated by the hypergeometric test, assuming normal data distribution. ( H ) GSEA terms of the 115 genes that were down-regulated in all 3 groups (versus WT) revealed the genetic programs copromoted by ATM and 53BP1-pS25. ATM, ataxia telangiectasia mutated; FDR, false discovery rate; GSEA, gene set enrichment analysis; KO, knockout; WT, wild type; 53BP1, p53 binding protein 1.
https://doi.org/10.1371/journal.pbio.3002760.g006
It remained unclear whether the higher expression of neuronal genetic programs in the 53BP1 mutants occurred in NPCs or neurons. Therefore, we compared the transcriptomes of 53BP1-S25A and S25D to WT NPCs, which had similar expression of NPC markers PAX6 and NES ( S12A Fig ). Up-regulated genetic programs in 53BP1-S25A and S25D NPCs shared categories such as translation control and ribosome ( S12B and S12C Fig ), whereas 53BP1-S25A NPCs also up-regulated cell cycle control and chromosome segregation ( S12C Fig ). Surprisingly, down-regulated genetic programs in 53BP1-S25A and S25D NPCs were highly enriched in neuronal differentiation ( S12D and S12E Fig ). The down-regulated genetic programs in NPCs are similar to neuronal programs that became up-regulated in 53BP1-S25A and S25D versus WT cortical organoids. These data suggest that 53BP1-S25 phosphorylation promotes the appropriate expression of neurogenic programs in NPCs and modulates the expression of the same programs in differentiating neurons in cortical organoids.
To dig deeper into analyses, we compared our data with previously published transcriptomic data that compared 53BP1 -KO and WT cortical organoids, which support a requirement of 53BP1 for activating neurogenic genes [ 6 ]. We observed that gene categories up-regulated by 53BP1-S25A and S25D were similar to those down-regulated in 53BP1 -KO cortical organoids. This was a significant overlap of 212 genes up-regulated by 53BP1-S25A and 53BP1-S25D with 53BP1-bound target genes that were down-regulated in 53BP1 -KO versus WT ( p = 0 by empirical estimation; Fig 6D ). The 212 genes were enriched in functions related to regulation of transcription, neurogenesis, neuronal projection, axonogenesis, synapse organization, and membrane depolarization ( Fig 6E ). This suggests that the expression of these genes is dependent on and up-regulated by 53BP1 phosphorylated at S25 in cortical organoids.
We next examined how transcriptomic changes in 53BP1-S25A and S25D compared to those in ATM -KO cortical organoids. We observed little overlap between the down-regulated genes in ATM -KO and the up-regulated genes in 53BP1-S25A and 53BP1-S25D. In contrast, we observed a greater overlap in concordant gene expression changes in ATM -KO, 53BP1-S25A, and 53BP1-S25D versus WT (Figs 6F , S12F , and S12G ). GSEA showed a significant enrichment of concordantly differentially expressed genes among ATM -KO, 53BP1-S25A, and 53BP1-S25D versus WT (Figs 6G and S13A-S13C ). Notably, all 3 mutant types shared down-regulated genes that were enriched functions related to TNFα signaling via NFκB, p53 pathway, IRE1-mediated unfolded protein response, FGFR signaling, TGFβ signaling, apoptosis, regulation of cell proliferation, and epithelial mesenchymal transition ( Fig 6G ). These data suggest that both ATM and 53BP1-pS25 promote the expression of these genes. From these findings, we can infer that ATM likely promotes the expression of these genes via phosphorylating 53BP1 at S25 in D35 cortical organoids. This suggests that ATM and 53BP1 may function together in a coordinated manner to regulate the expression of genes involved in critical signaling pathways and cellular processes during cortical development.
To obtain further mechanistic insights into the role of 53BP1 in controlling gene expression, we reanalyzed 53BP1 ChIP-seq data (using 2 separate anti-53BP1 antibodies) in WT NPCs [ 6 ]. Using SICER [ 27 ] and MACS2 [ 28 ] with a criterion of FDR <0.05, we identified 37,519 targets bound by 53BP1. About 41% of these 53BP1 targets localize to promoter regions, suggesting a transcriptional regulatory role of 53BP1 ( S14D Fig ). Remarkably, more than 82% of the differentially expressed genes in 53BP1-S25A and 53BP1-S25D D35 cortical organoids were found to be targets bound by 53BP1 (Figs 7A , S13E , and S13F ). 53BP1 target genes with increased transcript levels in the mutant organoids were highly enriched in neuronal development, axonogenesis, neuron projection, synapse organization, and neurotransmitter transport, transmission, and signaling ( Fig 7B ). On the other hand, 53BP1 targets with reduced transcript levels in the mutant organoids were enriched in IRE1-mediated unfolded protein response, cellular response to stress, iron import, and apoptosis regulation ( S13G Fig ). Of note, genes involved in IRE1-mediated unfolded protein response and apoptosis regulation showed reduced expression upon loss of ATM or mutation of 53BP1-S25 and were identified as direct targets of 53BP1 in NPCs. This suggests that ATM-mediated phosphorylation of 53BP1-S25 directly promotes the expression of these genes to maintain NPCs during formation of cortical organoids.
( A ) More than 82% of differentially expressed genes in 53BP1-S25A or S25D versus WT are chromatin targets bound by 53BP1 in WT NPCs. ( B ) 53BP1-S25A and S25D up-regulate 53BP1 targets that are involved in neuron development and projection, axonogenesis, synapse, and neurotransmitter synthesis and transport. ( C ) Heatmaps aligning peaks with 53BP1-pS25 CUT&RUN and 53BP1 ChIP-seq signals in WT NPCs. Input track was included as a negative control. n = numbers of peaks with differential and overlapped bindings. Criteria of FC>2 and p < 0.05 were used for comparison. ( D ) GSEA graph of 53BP1-pS25 CUT&RUN signals in genes that were lower in ESCs vs. NPCs, which were up-regulated in NPCs. P values were calculated by the hypergeometric test, assuming normal data distribution. Heatmaps aligning peaks with significantly different 53BP1 ChIP-seq signals in ( E ) 53BP1-S25A vs. WT and ( F ) 53BP1-S25D vs. WT, using the criterion of FC>2 and p < 0.05. Control peaks are those, after voom normalization, showed the least changes and served as semi-independent validation of differential ChIP-seq analysis. Bubble graphs present top enriched categories of genes that had significantly higher 53BP1 ChIP-seq in ( G ) 53BP1-S25A vs. WT and ( H ) 53BP1-S25D vs. WT. ESC, embryonic stem cell; FC, fold-change; GSEA, gene set enrichment analysis; NPC, neural progenitor cell; WT, wild type; 53BP1, p53 binding protein 1; 53BP1-pS25, 53BP1 phosphorylated at serine 25.
https://doi.org/10.1371/journal.pbio.3002760.g007
We wanted to test whether ATM alters 53BP1 binding, considering ATM is required for 53BP1-pS25 in D35 cortical organoids ( Fig 1F ) and NPCs ( S2B Fig ). A comparison of 53BP1 ChIP-seq in WT and ATM -KO NPCs showed that ATM -KO altered 53BP1 binding to chromatin ( S14A-S14C Fig ). ATM -KO reduced 53BP1 binding at specific sites, with 96.3% of these sites being promoters ( S14C Fig ). To explore the impact of 53BP1-pS25, we performed CUT&RUN in 2 separate WT NPC lines. Our analysis revealed that 67.1% of 53BP1-pS25 targets localize to promoter regions, suggesting a transcriptional regulatory role ( S14D Fig ). Under the criteria of fold-change >2 and p < 0.05, 58.6% (3,390/5,789) of 53BP1-pS25 targets overlapped with 53BP1 targets ( Fig 7C ); the nonoverlapped sites may be attributed to differences in ChIP versus CUT&RUN procedures and the accessibility of 53BP1 versus 53BP1-pS25 antibodies. 53BP1-pS25 targets were significantly enriched in 414 up-regulated genes in NPCs versus ESCs ( Fig 7D ), suggesting a role of 53BP1-pS25 in promoting their expression in NPCs. Genes having overlapped 53BP1 ChIP-seq and 53BP1-pS25 CUT&RUN signals were enriched in chromatin remodeling, DNA metabolism, RNA splicing, translation, transcription, cell cycle, and neuron development ( S14E Fig ). These data suggest that ATM can alter 53BP1 binding and that 53BP1-pS25 is enriched in the promoters of genes regulating cellular processes and neurodevelopment.
To investigate the impact of 53BP1-S25 on the genomic distribution of 53BP1, we performed ChIP-seq in 53BP1-WT, S25A, and S25D NPCs. Two independent NPC lines were used for each group, and the ChIP-seq data were subjected to principal component analysis, which showed high consistency between the replicate dataset ( S8A Fig ). We used SICER [ 27 ] and MACS2 [ 28 ] with the criteria of fold-change >2 and p < 0.05 to perform pairwise comparisons of the merged datasets from 53BP1-WT, S25A, and S25D ChIP-seq experiments. The pairwise comparisons identified thousands of 53BP1-bound regions that were significantly different between 53BP1-WT, S25A, and S25D. Notably, the regions that significantly gained binding in 53BP1-S25A or S25D versus WT were highly enriched at promoters (within 2 kb of transcription start sites), constituting 82% and 71.1%, respectively ( S8B and S8C Fig ). In contrast, the regions that significantly lost binding in 53BP1-S25A or S25D versus WT were not as enriched at promoters, constituting 32.6% and 33%, respectively ( S8B and S8C Fig ). We generated heatmaps to visualize the genomic regions with significantly different 53BP1 binding intensity (compared against control regions). The heatmaps confirmed consistent changes in 53BP1 binding patterns between 53BP1-S25A and S25D versus WT, and between 53BP1-S25A versus S25D (Figs 5C , 5D , and S8D ). These data support that 53BP1-S25 and its phosphorylation control the genomic distribution of 53BP1 on chromatin.
We next set out to examine the correlation between changes in 53BP1 distribution on chromatin and changes in gene expression in 53BP1-WT, S25A, and S25D cortical organoids. We performed GSEA and made some notable observations. Firstly, regions that gained 53BP1 binding in 53BP1-S25A or S25D cortical organoids, as compared to WT, were enriched with up-regulated genes ( Fig 5E and 5F ). Similarly, regions that had lower 53BP1 binding were enriched with down-regulated genes in 53BP1-S25A or S25D versus WT ( S8E and S8F Fig ). These results suggest that the 53BP1-S25A or S25D mutation directly influences 53BP1 binding and gene expression and subsequently regulates gene expression, particularly at promoters where higher 53BP1 binding leads to higher gene expression.
Interestingly, genes that lost 53BP1-S25A or S25D protein binding had minimal overlap in GSEA terms, except promoters occupied with H3K4me3 and regulation of epithelial-mesenchymal transition ( S8E and S8F Fig ). In contrast, the genes that gained 53BP1-S25A or S25D protein binding were enriched with promoters marked by bivalent histone marks (H3K4me3 and H3K27me3) or occupied by H3K27me3 alone [ 29 ] ( Fig 5E and 5F ), suggesting that 53BP1-S25D or S25A proteins preferentially bind to these promoters and subsequently up-regulate gene expression. Moreover, genes that gained 53BP1-S25A or S25D binding shared common functions related to sodium ion transmembrane transporter, DNA replication, positive regulation of cell division, and regulation of histone H3K4 methylation ( Fig 5E and 5F ). This suggests that despite the 1,187 regions showing different 53BP1 bindings between 53BP1-S25A and S25D ( S8D Fig ), both mutations impact genes involved in neuronal functions and cell proliferation. Altogether, these findings show that 53BP1-S25A and S25D mutations have a direct impact on 53BP1 binding to chromatin and subsequently affecting gene regulation. We propose that 53BP1-pS25 likely inhibits 53BP1 binding to promoters associated with bivalent and H3K27me3-occupied promoters. This inhibition may lead to the reduced expression of genes involved in the regulation of H3K4me3, neuronal functions, and cell proliferation.
We next tried to identify a regulation of ATM, whose protein levels increased in NPCs ( Fig 1D ). This led us to test whether and how inhibitors of TGFb, WNT, and HH signaling control protein levels of ATM, 53BP1, and pS25-53BP1 by removing one inhibitor at a time from the cortical organoid differentiation media ( S16A Fig ). As we could not successfully identify physical presence of ATM at promoters, WB analysis is most apt to study ATM level and activity. By day 4 of neural differentiation, although 53BP1 protein levels were reduced by the withdrawal of SB431542 (TGFβ inhibitor) or IWR1-endo (WNT inhibitor), pS25-53BP1 was not altered ( S16B and S16C Fig ). By day 10 of neural differentiation, the withdrawal of cyclopamine (HH inhibitor) reduced pS25-53BP1 level (but not ATM or 53BP1 proteins; S16D and S16E Fig ). These signaling pathways may affect pS25-53BP1 or ATM activities during neural differentiation.
Next, we tested whether another DNA damage response factor, apart from ATM, influences 53BP1-pS25. RNF168 plays a central role in the γH2AX-MDC1-RNF8-RNF168-H2AK15ub axis, which governs the binding of 53BP1 to chromatin with DNA damage [ 30 ]. We generated RNF168 -KO hESC clone 44, which maintained pluripotency and genome integrity ( S17A-S17D Fig and S1 Table ). RNF168 -KO hESCs were differentiated to NPCs, which expressed NPC markers similar to WT NPCs ( S17E Fig ). RNA-seq analysis comparing 2 datasets each from RNF168 -KO44 and WT NPCs revealed that up-regulated genes were enriched in neuronal differentiation, translation and ribosome, and cell cycle transition ( S17F Fig ), while down-regulated genes were enriched in cilium movement, H3K27me3 targets, H3K4me3 targets, astrocyte markers, signaling pathways, and positive regulation of NPC proliferation ( S17G Fig ). We performed 53BP1-pS25 CUT&RUN and showed that RNF168 -KO disrupted 53BP1-pS25 localization on chromatin ( S17H Fig ). RNF168 -KO increased 53BP1-pS25 levels at genes enriched in neuronal differentiation, cell morphogenesis, and stem cell maintenance, whereas RNF168 -KO decreased 53BP1-pS25 levels at genes enriched in cell cycle transition, signaling receptor regulation, anterior-posterior patterning, and transcription activator ( S17I Fig ). The altered 53BP1-pS25 localization correlated with differential gene expression in RNF168-KO versus WT NPCs ( S17J Fig ). Altogether, these data suggest that DNA damage signaling regulates 53BP1 binding to chromatin, affecting genetic programs related to signaling pathways, protein translation, and NPC proliferation and differentiation.”
In our study, we made significant discoveries regarding the role of ATM and 53BP1-pS25 in controlling gene expression during the differentiation of hESCs into cortical organoids. We revealed that ATM exerts a strong influence over various aspects of gene regulation, including transcriptional, posttranscriptional, and translational control. While our in vitro model may not fully recapitulate neurodevelopment in vivo, it provides valuable insights into corticogenesis. We have shown that neural differentiation promotes ATM protein levels, and ATM-dependent phosphorylation predominantly impacts factors involved in neurogenesis, neuronal differentiation, cell morphogenesis, and microtubule cytoskeleton. Dysregulation of these processes led to the cellular defects in ATM -KO cortical organoids. We showed that key signaling pathways may affect ATM during neural induction. The activity of ATM can be regulated by DNA damage response, reactive oxygen species, hypothxia, hypothermia, and phosphatase WIP1 [ 31 – 33 ]. The exact clarification of mechanisms promoting ATM activities, especially in directing its kinase activity at specific promoters, is beyond the scope of this study. Additionally, we have identified kinases involved in ATM, BDNF, and WNT signaling, G2/M checkpoint, and p53 regulation as being influenced by ATM-dependent phosphorylation during cortical organoid differentiation. These molecular pathways may function in diseases associated with ATM, including ataxia telangiectasia [ 34 – 36 ].
We recognized the diverse effects of ATM and decided to focus our studies on 53BP1-pS25, a phosphorylation event dependent on ATM. We found that 53BP1-pS25 regulates genetic programs including signaling pathways, p53 regulation, apoptosis, and cell proliferation. To understand the mechanisms underlying 53BP1’s involvement in gene regulation, we built a model that incorporates current knowledge about 53BP1 functions in the DNA damage response. We propose that ATM phosphorylates H2AX at transcription start sites [ 13 , 14 ], facilitating the recruitment of 53BP1 and subsequent phosphorylation of 53BP1-S25. RNF168, key to DNA damage response signaling [ 30 ], also regulates 53BP1-pS25 on chromatin and genetic programs crucial to neural differentiation. Phosphorylation of 53BP1-S25 inhibits the recruitment of 53BP1 to bivalent or H3K27me3-occupied promoters for suppressing the expression of genes involved in the regulation of H3K4me3, neuronal functions, and cell proliferation. The fidelity of gene expression in cortical brain organoids requires dynamic changes in the phosphorylation of 53BP1-S25. This process is likely to involve the interactions of 53BP1 with other proteins, including RIF1, SCAI, and UTX [ 6 , 11 , 12 ]. These interactors have known roles in chromatin alterations and gene regulation. Notably, UTX is an H3K27me3 demethylase that can modify bivalent or H3K27me3-occupied promoters and has been shown to partner with 53BP1 to promote neurogenesis in humans but not in mice [ 6 ]. Given our findings, we propose that 53BP1-pS25 may influence the activities of 53BP1–UTX at bivalent or H3K27me3-occupied promoters, thus modulating gene expression and contributing to the timing of neuronal differentiation.
Our studies have uncovered the remarkable role of ATM–53BP1 in regulating neurodevelopmental programs. Its impact is multifaced. Firstly, ATM–53BP1 plays a crucial role in maintaining NPCs and controlling the size of cortical organoids. Secondly, ATM–53BP1 is involved in driving and modulating programs related to synapse formation, axon development, and neurotransmitter regulation, processes fundamental for establishing neuronal networks and communication within the brain. Thirdly, our findings reveal a temporal component in the regulation of neurodevelopmental programs by ATM–53BP1. As cortical organoids progress in differentiation, there is a temporal regulation of neuronal differentiation and function. This switch involves ATM and the 53BP1-pS25 dynamics to specifically control genes associated with synapse, axon, and neurotransmitter, which are crucial to cognition. In the future, elucidation of this mechanism will provide valuable insights into the molecular control of corticogenesis. Beyond 53BP1, ATM-dependent phosphorylation likely controls many other key neurodevelopmental regulators. Future studies of how ATM selects substrates to exert its multiple influences will significantly advance our understanding of the epigenetic programming underlying human neurodevelopment.
PBS: 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4 (pH 7.4). PBST: PBS with 0.1% Triton X-100. HEPM: 25 mM HEPES (pH 6.9), 10 mM EGTA, 60 mM PIPES, 2 mM MgCl2. Immunofluorescence blocking solution: 1/3 Blocker Casein (Thermo Fisher Scientific), 2/3 HEPM with 0.05% TX-100. Buffer A: 10 mM HEPES (pH 7.9), 10 mM KCl, 1.5 mM MgCl2, 0.34 M sucrose, 10% glycerol. Buffer B: 3 mM EDTA, 0.2 mM EGTA. Buffer D: 400 mM KCl, 20 mM HEPES, 0.2 mM EDTA, 20% glycerol. ChIP lysis buffer 3: 10 mM Tris-HCl (pH 8.0), 100 mM NaCl, 1 mM EDTA, 0.5 mM EGTA, 0.1% sodium deoxycholate, 0.5% N-Lauroylsarcosine. ChIP wash buffer: 50 mM HEPES (pH 7.5), 500 mM LiCl, 1 mM EDTA, 1% NP-40, 0.7% Na-deoxycholate. ChIP elution buffer: 50 mM Tris-HCl (pH 8.0), 10 mM EDTA, 1% SDS. CUT&RUN Wash buffer: 20 mM HEPES (pH 7.5), 150 mM NaCl, 0.5 mM spermidine, protease inhibitor cocktail (Sigma-Aldrich 11873580001). CUT&RUN Binding buffer: 20 mM HEPES-KOH (pH 7.9), 10 mM KCl, 1 mM CaCl2, 1 mM MnCl2. CUT&RUN Digitonin buffer: CUT&RUN Wash buffer with 0.01% digitonin. CUT&RUN Antibody buffer: CUT&RUN Digitonin buffer with 2 mM EDTA. CUT&RUN 2X Stop buffer: 340 mM NaCl, 20 mM EDTA, 4 mM EGTA. CUT&RUN Stop buffer: Into 1 mL of 2X Stop buffer stock, add 5 μL of 10 mg/mL RNase A and 3.3 μL of 15 mg/mL. GlycoBlue Coprecipitant (Thermo Fisher AM9516)
S7 Table lists all antibodies and conditions used in this study.
H9/WA09 (WiCell) hESCs were grown on Matrigel with reduced growth factors (Thermo Fisher Scientific, #35423) in mTeSR1 medium (STEMCELL Technologies, #85850) at 37°C and 5% CO 2 . The 53BP1 knock-in cell lines (53BP1 S25A 34–3, 34–4, 79–1, 79–3 and S25D 14–3, 14–15, 14–19, 17) and ATM KO cell lines (ATM-KO2, 3, 14, and 43) were generated using CRISPR/Cas9 gene-editing technology. Genome editing reagents were designed and validated in the Center for Advanced Genome Engineering at St. Jude Children’s Research Hospital. Briefly, a chemically modified sgRNA (Synthego) was precomplexed with SpCas9 protein (St. Jude Protein Production Core) and cotransfected with an ssODN donor template containing the desired modification into H9/WA09 cells via nucleofection (Amaxa P3 primary cell 4D nucleofector X kit L, Lonza) using the manufacturer’s recommended protocol. Transfected cells were sorted (BD FACSAria Fusion) onto Matrigel and allowed to grown single-cell clones. Clones were identified via targeted mi-seq using a 2-step PCR library setup as previously described [ 37 ]. Samples were demultiplexed using the index sequences, fastq files were generated, and NGS analysis was performed using CRIS.py [ 38 ]. S8 Table lists genome editing reagents and associated primers.
ESCs were seeded onto AggreWell800 plates (STEMCELL Technologies, #34811) and fed with neural induction medium (STEMCELL Technologies, #05835) to form embryoid bodies. On day 5, embryoid bodies were replated onto Matrigel-treated 6-well plates in the same media. On day 17, cells were harvested as NPCs.
ESCs and NPCs were incubated in Buffer A + PI + DTT for 5 min on ice. After centrifugation at 1,750 g for 2 min at 4°C, the nuclei pellet was washed in Buffer A and subsequently incubated for approximately 25 min in Buffer D + PI + DTT at 4°C with rotation to obtain the nuclear fraction. Nuclear extracts were separated by SDS–PAGE and transferred onto a nitrocellulose membrane (Bio-Rad). Membranes were blocked with 3% bovine serum albumin (BSA) in HEPM, incubated in primary antibodies (HEPM containing 1% BSA and 0.1% Triton X-100) overnight at 4°C, washed in PBS-T, incubated in IRDye-conjugated secondary antibodies (LI-COR), and imaged on an Odyssey Fc imaging system (LI-COR). Signals were quantitated with the Image Studio software (version 1.0.14; LI-COR).
Antibody was bound to protein A and protein G Dynabeads (Thermo Fisher 10002D and 10004D) for 2 h at room temperature. Nuclear extract was incubated with the Dynabeads-antibody complex for 5 h at 4°C, washed with PBST, and eluted with 0.1 M glycine (pH 2.3). Eluates were neutralized with 1/10 volume of 1.5 M Tris buffer (pH 8.8).
Cortical organoids were generated based on previously published methods with minor modifications [ 39 , 40 ]. In brief, hESC lines were expanded and dissociated to single cells using Accutase, seeded onto low-attachment V-bottom 96-well plates (Costar, #7007) at a density of 9,000 cells per well to aggregate into embryoid bodies. The embryoid bodies formation medium (DMEM/F-12 with 20% KO serum replacement, 3% ESC-quality FBS, 2 mM GlutaMAX, 0.1 mM nonessential amino acids) was supplemented with dorsomorphin (2 μM), WNT inhibitor (IWR1, 3 μM), TGF-β inhibitor (SB431542, 5 μM), and Rho kinase inhibitor (Y-27623, 20 μM). Starting from day 4, embryoid bodies were fed with cortical differentiation medium (Glasgow-MEM, 20% KSR, 0.1 mM NEAA, 1 mM sodium pyruvate, 0.1 mM β-ME, and 1% anti-anti), supplemented with WNT inhibitor (IWR1, 3 μM), TGF-β inhibitor (SB431542, 5 μM), cyclopamine (2.5 μM) and Rho kinase inhibitor (Y-27623, 20 μM). On day 17, embryoid bodies were embedded in Matrigel droplets and transferred onto low-attachment 6-wells and cultured in suspension using DMEM/F-12 supplemented with 1% N2 supplement, 1% lipid concentrate, 2% B27 supplement without vitamin A, and 1% anti-anti under 40% O 2 /5% CO 2 conditions on shaker. Starting from day 30, medium was changed to 50% DMED/F-12, 50% neurobasal media, 0.5% N2 supplement, 1% GlutaMax, 0.05 mM NEAA, 0.025% human insulin, 0.1 mM β-ME, and 1% anti-anti, supplemented with 2% B27.
Cells and cryosectioned organoids were blocked with IF blocking solution for 2 h at room temperature and primary antibodies (diluted in blocking buffer) added and incubated O/N at 4°C. After 3 washes in PBS-T, fluorescent dye-conjugated secondary antibodies (1:500, Alexa Fluor-CONJUGATED antibodies, Thermo Fisher Scientific) were added and incubated for 3 h at room temperature. Secondary was washed with PBS-T 3 times, and samples were washed and coverslips mounted with Prolong Glass Mounting Reagent (Thermo Fisher Scientific), which contains DAPI. Images were acquired with Zeiss LSM780.
At days 35 and 55, bright-field images of organoids were captured with Axiocam 208 (Zeiss). Area of organoids, area of ventricular zone–like regions, and marker-positive cells were quantified by using the software FIJI: Signals-positive cells were identified based on signal and width thresholds. For ventricular zone–like region quantification, inner and outer edges of the regions in the image were manually traced, based on CTIP2-positive cells encircling the outer edges. FIJI was used to quantify area, perimeter, major and minor axes of the inner and outer traces. Mean perimeter and the difference between the major axes of the inner and outer traces were used to estimate the thickness of the structure. Mean Perimeter = (outer perimeter + inner perimeter) / 2. MajorAxisDiff = (outer major axis − inner major axis) / 2. MinorAxisDiff = (outer minor axis − inner minor axis) / 2. To quantify ZO-1-positive ventricular surfaces, ZO-1 signals were normalized by the Integral Image Filters plugin, and surface areas were manually traced for quantification. The VZ/SVZ structure was considered organized if PAX6-positive nuclei were densely packed with radial organization around ZO-1-positive ventricular surfaces. Ilastik [ 41 ] was used to quantify nuclear areas positive for different markers, using segmentation via a machine learning-based package and area quantification of segmented areas. Marker ratios were then calculated based on quantified areas.
Nine to 12 organoids of each line were dissociated using the papain dissociation system (Worthington LK003153). Dissociated cells were fixed in 4% formaldehyde solution at 4°C overnight and washed once in 1X PBS. Then, cells were permeabilized in 1X PBST for 2 h at room temperature on an orbital shaker. Cells were blocked in IF blocking buffer (1/3 Blocker Casein (Thermo Fisher 37528), 2/3 HEPM with 0.05% Triton X-100) for 2 h at room temperature on a shaker. Primary antibodies in IF blocking buffer were mixed with cells at 4°C overnight followed by washing twice with 1X PBST. Secondary antibodies in IF blocking buffer were mixed with cells for 2 h at room temperature on a shaker. After washing cells once, a conjugated antibody was added and incubated for 2 h at room temperature on a shaker. Cells were washed one last time before resuspended in 1X PBS for FACS. FACSymphony A1 sorter was used for analysis. All the centrifugation steps were done at 500 × g for 4 min at room temperature. All washes were performed by incubating the cells with 1X PBS (after fixation) or PBST (after antibody staining) for 5 min at room temperature on a shaker. Primary antibodies used are Ki67 (Cell Signaling 9129), PAX6 (DSHB supernatant 1mL), CTIP2 (Abcam 18465), and cleaved Caspase3-AF405-conjugated (R&D Systems IC835V).
Total RNA was extracted with TRIzol reagent (Invitrogen, #15596026) and Direct-zol RNA Microprep (Zymo Research, # R2062) by following manufacturer’s instructions. DNA digestion with DNase I was performed during RNA extraction. Paired-end 100-cycle sequencing was performed on NovaSeq6000 sequencer by following the manufacturer’s instructions (Illumina). Raw reads were first trimmed using TrimGalore (version 0.6.3) available at: https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/ , with parameters ‘—paired—retain_unpaired’. Filtered reads were then mapped to the Homo sapiens reference genome (GRCh38 + Gencode-v31) using STAR (version 2.7.9a) [ 42 ]. Gene-level read quantification was done using RSEM (version 1.3.1) [ 43 ]. To identify the differentially expressed genes between control and experimental samples, the variation in the library size between samples was first normalized by trimmed mean of M values (TMM) and genes with CPM < 1 in all samples were eliminated. Then, the normalized data were applied to linear modeling with the voom from the limma R package [ 44 ]. GSEA was performed against using the MSigDB database (version 7.1), and differentially expressed genes were ranked based on log 2 (FC) [ 45 , 46 ].
Organoids were harvested on day 35, and the Matrigel droplets were eliminated by multiple ice-cold PBS washes. The organoid pellet was extracted in the lysis buffer (50 mM HEPES (pH 8.5), 8 M urea, and 0.5% sodium deoxycholate, 100 μl buffer per 10 mg tissue) with 1x PhosSTOP phosphatase inhibitor cocktail (Sigma-Aldrich). Protein concentration was estimated by a Coomassie stained short gel with BSA as a standard. About 600 μg each of protein samples was digested with LysC (Wako) at an enzyme-to-substrate ratio of 1:100 (w/w) for 2 h at room temperature in the presence of 1 mM DTT. The samples were then diluted to a final 2 M urea concentration with 50 mM HEPES (pH 8.5) and digested with Trypsin (Promega) at an enzyme-to-substrate ratio of 1:50 (w/w) for 3 h. The peptides were reduced by adding 1 mM DTT for 30 min at room temperature followed by alkylation with 10 mM iodoacetamide for 30 min in the dark at room temperature. The unreacted iodoacetamide was quenched with 30 mM DTT for 30 min. Finally, the digestion was terminated and acidified by adding trifluoroacetic acid to 1%, peptides desalted using Sep-Pak C18 cartridge (Waters), and dried by speed vac. The purified peptides were resuspended in 50 mM HEPES (pH 8.5) and labeled with 16-plex Tandem Mass Tag (TMTpro) reagents (Thermo Scientific) following the manufacturer’s recommendation. The TMT labeled samples were mixed equally, desalted using Sep-Pak C18 cartridge (Waters), and dried by speed vac.
The dried TMT mix was resuspended and fractionated on an offline HPLC (Agilent 1220) using basic pH reverse phase liquid chromatography (pH 8.0, XBridge C18 column, 4.6 mm × 25 cm, 3.5 μm particle size, Waters). A total of 160 one-minute fractions were collected and concatenated to 80 fractions. For whole proteome analysis, 10% of these 80 fractions was used. The remaining 90% of the 80 fractions were concatenated to 20 fractions for phophopeptide enrichment. Phosphopeptide enrichment was performed according to a previously published protocol [ 47 ]. The phosphopeptide enrichment eluents and the total proteome fractions were dried and resuspended in 5% formic acid and analyzed by acidic pH reverse phase LC-MS/MS analysis. The peptide samples were loaded on a nanoscale capillary reverse phase C18 column (New objective, 75 μm ID × approximately 15 cm, 1.9 μm C18 resin from Dr. Maisch GmbH) by a HPLC system (Thermo Ultimate 3000) and eluted by either a 125-min gradient (phosphofractions) or 110-min gradient for total proteome fractions. The eluted peptides were ionized by electrospray ionization and detected by an inline Orbitrap Fusion mass spectrometer (Thermo Scientific). For total proteome fractions, the mass spectrometer is operated in data-dependent mode with a survey scan in Orbitrap (60,000 resolution, 2 × 10 5 AGC target and 50 ms maximal ion time) and MS/MS high-resolution scans (60,000 resolution, 1 × 10 5 AGC target, 150 ms maximal ion time, 36.5 HCD normalized collision energy, 1 m/z isolation window, and 15-s dynamic exclusion). For phosphoproteome fractions, the mass spectrometer is operated in data-dependent mode with a survey scan in Orbitrap (60,000 resolution, 3 × 10 5 AGC target and 50 ms maximal ion time) and MS/MS high-resolution scans (60,000 resolution, 1 × 10 5 AGC target, 150 ms maximal ion time, 36.5 HCD normalized collision energy, 1 m/z isolation window, and 10-s dynamic exclusion).
The MS/MS raw data were processed by a tag-based hybrid search engine JUMP [ 48 ]. The data were searched against the UniProt human database (168,305 protein entries; downloaded in April 2020) concatenated with a reversed decoy database for evaluating FDR. Searches were performed using a 15-ppm mass tolerance for fragment ions, fully tryptic restriction with 2 maximal missed cleavages, 3 maximal modification sites, and the assignment of b and y ions. TMT tags on Lysine residues and N-termini (+304.2071453 Da) were used for static modifications and Met oxidation (+15.99492 Da) was considered as a dynamic modification. Phosphorylation (+79.96633 Da) was considered as a dynamic modification for STY residues. Putative peptide spectral matches (PSMs) were filtered by mass accuracy and then grouped by precursor ion charge state and filtered by JUMP-based matching scores (Jscore and ΔJn) to reduce FDR below 1% for proteins during the whole proteome analysis or 1% for phosphopeptides during the phosphoproteome analysis. Phosphosites were further evaluated by JUMPl program using the concept of the phosphoRS algorithm [ 49 ] to calculate phosphosite localization scores (Lscore, 0% to 100%) for each PSM.
TMT reporter ion intensities of each PSM were extracted and corrected based on isotopic distribution of each labeling reagent. Those PSMs with very low intensities (e.g., minimum intensity of 1,000 and median intensity of 5,000) were excluded for quantification. Sample loading bias was mitigated by normalization with the trimmed median intensity of all PSMs. Protein or phosphopeptide relative intensities were calculated by dividing the intensity of each channel by the mean intensity. Protein or phosphopeptide absolute intensities were computed by multiplying the relative intensities by the grand-mean of 3 most highly abundant PSMs.
Differentially expressed proteins between the 2 strains and 2 different doses were identified by the limma R package [ 50 ]. The Benjamini–Hochberg method was used to control multiple-testing correction, and proteins with an adjusted p -value of <0.05 and log2 fold change of >1.5 were defined as differentially expressed.
Pathway enrichment analysis was carried out to infer functional groups of proteins that were enriched in a given dataset. The 4 common pathway databases were used, including Gene Ontology (GO), KEGG, Hallmark, and Reactome. The analysis was performed using Fisher’s exact test with the Benjamini–Hochberg correction for multiple testing. A cutoff of adjusted p -value < 0.2 was used to identify significantly enriched pathways.
Kinase activity was inferred based on known substrates in the PhosphoSitePlus database [ 51 ] using the IKAP algorithm [ 24 ]. The phosphoproteome data were normalized against the whole proteome. We performed 100 times of calculations to overcome the potential problem of local optimization.
Cells were harvested in PBS. Cytoplasmic fractions were extracted using buffer A with 1× protease inhibitors and 1 mM DTT. Nuclear pellets were cross-linked by 1.1% formaldehyde in buffer B with 1× protease inhibitors and 1 mM DTT; washed; and lysed in lysis buffer 3 with 1× protease inhibitors, 1 mM DTT, and 1 mM PMSF. The fixed and lysed nuclear extract was sonicated with Bioruptor Pico (Diagenode) 10 times for 15 s each, with 45-s intervals. Chromatin was added to Dynabeads (Life Technologies) prebound with 4 μg of antibodies for overnight incubation. After incubation, beads were washed and immunoprecipitates were eluted. DNA from eluates was recovered by the GeneJET FFPE DNA purification kit (Thermo Fisher Scientific, #K0882). DNA libraries were generated using the NEBNext Ultra DNA Library Prep kit (NEB, #E7370S) and sequenced at the St. Jude Hartwell Center.
Approximately 5 × 10 5 live cells were mixed with 5 × 10 4 Drosophila S2 cells per reaction. For CUT&RUN, we followed EpiCypher CUTANA protocol. In brief, we first isolated nuclei by incubating cells on ice for 5 min in Buffer A with protease inhibitor and 0.1% Triton X-100. After centrifugation at 1,750 × g for 2 min at 4°C, nuclei were resuspended in Wash buffer. Bio-Mag Plus Concanavalin-A (Con A) coated beads (Bangs Laboratories BP531) activated in Binding buffer were then added to the nuclei and rotated for 10 min at room temperature. About 1 μg primary antibody with 0.25 μg Spike-in antibody (Active Motif 61686) diluted in Antibody buffer was added to the bead-nuclei mixture and incubated for 2 h at room temperature. Beads were washed twice with Digitonin buffer and incubated with pAG-MNase for 10 min at room temperature. Beads were then washed twice with Digitonin buffer, incubated with 2 mM CaCl 2 for 2 h at 4°C, and quenched by adding Stop buffer. DNA was released from the beads by incubating them for 10 min at 37°C and purified by CUTANA DNA purification kit (EpiCypher SKU:14–0050). Libraries were constructed using xGen ssDNA and Low-Input DNA Prep by following the manufacturer’s instructions (IDT 10009817) and sequenced at the St. Jude Hartwell Center.
Approximately 50 bp single-end reads were obtained and aligned to human genome hg38 by BWA (version 0.7.170.7.12, default parameter). Duplicated reads were marked by the bamsormadup from the biobambam tool (version 2.0.87) available at https://www.sanger.ac.uk/tool/biobambam/ . Uniquely mapped reads were kept by samtools (parameter “-q 1 -F 1804,” version 1.14). Fragments <2,000 bp were kept for peak calling, and bigwig files were generated for visualization. SICER [ 27 ] and macs2 [ 28 ] were both used for peak calling to identify both the narrow and broad peak correctly. With SICER, we assigned peaks that were at the top 1 percentile as the high-confidence peaks and the top 5 percentile as the low-confidence peaks. Two sets of peaks were generated: Strong peaks called with parameter “FDR < 0.05” by at least 1 method (macs2 or SICER) and weak peaks called with parameter “FDR < 0.5” by at least 1 method (macs2 or SICER). Peaks were considered reproducible if they were supported by 1 strong peaks and at least 1 weak peak in other replicates. For downstream analyses, heatmaps were generated by deepTools [ 52 ], and gene ontology was performed with Enrichr [ 53 , 54 ] and GSEA, in addition to custom R scripts. For differential peak analysis, peaks from 2 replicates were merged and counted for number of overlapping extended reads for each sample (bedtools v2.24.0) [ 55 ]. Then, we detected the differential peaks by the empirical Bayes method (eBayes function from the limma R package) [ 44 ]. For downstream analyses, heatmaps were generated by deepTools (v3.5.0) [ 56 ]. Peaks were annotated based on Gencode following this priority: “Promoter.Up”: if they fall within TSS– 2 kb, “Promoter.Down”: if they fall within TSS– 2 kb, “Exonic” or “intronic”: if they fall within an exon or intron of any isoform, “TES peaks”: if they fall within TES ± 2 kb, “distal5” or “distal3” if they are with 50 kb upstream of TSS or 50 kb downstream of TES, respectively, and they are classified as “intergenic” if they do not fit in any of the previous categories.
S1 data. numerical data used to generate summary data in this study..
https://doi.org/10.1371/journal.pbio.3002760.s001
https://doi.org/10.1371/journal.pbio.3002760.s002
( A ) Schematic diagram of neural differentiation of hESCs: neural induction, differentiation, and maturation media to form EBs, rosettes, NPCs, and neurons. ( B ) Principal component analysis of WT ESCs, NPCs, day 10 (D10) cortical organoids, and D17 cortical organoids. GSEA terms that are highly enriched in significantly ( C ) down-regulated and ( D ) up-regulated genes in WT NPCs compared to ESCs. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. ( E ) Immunofluorescence of NPC markers PAX6 and NESTIN. Bar, 50 μm. ( F ) Quantification of 53BP1-pS25-positive hESCs or hNPCs. Data are presented as the mean ± SEM, with p < 0.0001. ( G ) WB analysis of control cells and 53BP1-KO clones 415, 416, and 209, which are clones KO1, KO2, and KO3 in Yang and colleagues’ study [ 6 ]. (H) WB analysis of control and 53BP1-S25A hNPCs. The S25A mutation prohibits phosphorylation. ( I ) WB analysis of hESCs and hNPCs and quantification. ( J ) Schematic diagram of genome editing in hESCs. Guide RNA 6 were complexed with Cas9 proteins and used along single-stranded nucleotide donors to transfect hESCs. Individual clones from transfection were cultured, sequenced by mi-seq across the targeted 53BP1 locus, and established as >99% pure clonal lines. Diagram was generated using open-sourced images available at biorender.com . Underlying numerical values for figures are found in S1 Data . EB, embyoid body; ESC, embryonic stem cell; GSEA, gene set enrichment analysis; hESC, human embryonic stem cell; hNPC, human neural progenitor cell; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; WB, western blot; WT, wild type; 53BP1-pS25, 53BP1 phosphorylated at serine 25.
https://doi.org/10.1371/journal.pbio.3002760.s003
( A ) Alignment of WT and ATM -KO mutation sequences on 2 alleles (al) in the ATM locus. Red indicates the gRNA sequence. ( B ) WB analysis of WT and 4 ATM -KO hNPCs. ( C ) Principal component analysis showed the intermixing and similar RNA-seq profiles from hESCs of 7 WT, 4 53BP1-S25A, 4 53BP1-S25D, 4 ATM-KO, and 4 53BP1-KO lines. ( D ) Immunofluorescence showed similar expression of OCT4 and SSEA4 proteins in control and ATM -KO hESCs. Bar, 100 μm. ( E ) WB analysis of WT and 2 ATM -KO hNPCs. Quantification suggests reduction of γH2AX in ATM -KO hNPCs. Welch’s t test was used to perform pairwise comparisons of WT and ATM -KO. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; hESC, human embryonic stem cell; hNPC, human neural progenitor cell; KO, knockout; WB, western blot; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s004
( A ) Immunofluorescence showed D35 ATM -KO and WT cortical organoids had similar γH2AX foci. Bar, 100 μm. FACS analysis of CC3 in ( B ) D21 and ( C ) D28 cortical organoids. Two biological replicates were done, and each data point was based on 3 technical replicate analyses of 10–12 cortical organoids. ( D , E ) Immunofluorescence and quantification of CC3 in D28 cortical organoids. Bar, 100 μm. Graphs are presented in ratios (out of 1), with **, p < 0.01; ****, p < 0.0001; ns, not significant by two-way ANOVA test. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; CC3, cleaved-caspase 3; KO, knockout; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s005
FACS analysis of PAX6 and KI67 in ( A , B ) D28 and ( C ) D35 cortical organoids. Each data point was based on the 3 technical replicate analyses of 10 cortical organoids. ( D ) Quantification of KI67/PAX6 ratios in immunofluorescence of D35 cortical organoids. Each data point represents quantification of cells in 1 cortical organoid. ( E , F ) Immunofluorescence and quantification of H3-pS10 (PH3) in D28 cortical organoids. Bar, 100 μm. ( F - H ) Immunofluorescence and quantification of PH3 and KI67 in D35 cortical organoids. Bar, 100 μm. ***, p < 0.001; ns, not significant by two-way ANOVA test. Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.s006
Immunofluorescence of ( A ) NEUN and ( E ) ZO-1 and PAX6 in D37 cortical organoids. Bar, 100 μm. ( B ) Immunofluorescence of ZO-1 in D28 cortical organoids. Bar, 100 μm. Quantification of the ( C ) number and ( D ) surface area of ZO-1-positive ventricles in D28 cortical organoids. *, p < 0.05; ***, p <0 .001; ns, not significant by two-way ANOVA test. ( F ) Bright-field images of cortical organoids formed by ATM -KO2, 3, 14, 43, and WT control at day 55 of differentiation. Bar, 1.5 mm. ( G ) The size of cortical organoids was compared between groups by one-way ANOVA with Dunnett’s multiple comparisons test, with ns, not significant and ***, p < 0.001. n = 13 organoids/group. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; NPC, neural progenitor cell; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s007
( A ) Immunofluorescence of PAX6 and NES in NPCs. Bar, 50 μm. ( B ) Principal component analysis of proteomics data of D35 WT and ATM -KO cortical organoids. GSEA terms that are highly enriched in significantly ( C ) higher and ( D ) lower total proteins in D35 ATM -KO versus WT cortical organoids. ( E ) GSEA terms that are highly enriched in significantly higher phosphoproteins, which were normalized to total proteomics, in D35 ATM -KO versus WT cortical organoids. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; GSEA, gene set enrichment analysis; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s008
Heatmaps showing relative phosphorylation levels of ( A ) 7 MAPK9 substrates that are significantly lower and ( B ) 7 CDK5 substrates that are significantly higher in D35 ATM -KO versus WT cortical organoids. ( C ) Heatmaps showing activity of selected protein kinases between ATM-KO3, ATM-KO4, and WT cell lines. ( D ) Alignment of WT and 53BP1-S25A and S25D mutation sequences on 2 alleles (al). Red indicates the gRNA sequence. Underline indicates codon encoding the WT serine 25, mutant alanine, or mutant aspartic acid. ( E ) WB analysis of control and 53BP1-S25D hNPCs, which have comparable levels of 53BP1 protein. ( F ) Transcripts per million values of 10 pluripotency genes were used for comparison to show that control, 53BP1-S25A, and 53BP1-S25D hESCs did not differ in pluripotency. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; hESC, human embryonic stem cell; hNPC, human neural progenitor cell; KO, knockout; WB, western blot; WT, wild type; 53BP1, p53 binding protein 1.
https://doi.org/10.1371/journal.pbio.3002760.s009
( A ) Immunofluorescence showed similar expression of OCT4 and SSEA4 proteins in WT, 53BP1-S25A, and 53BP1-S25D hESCs. Bar, 100 μm. Immunofluorescence of ( B ) KI67 and ( D ) PH3 in cryosections of cortical organoids at day 35 of differentiation. Bar, 100 μm. ( C ) Quantification of KI67-positive cells in D35 cortical organoids. Data points represent single organoids. The mean ± SEM values were compared by one-way ANOVA with Dunnett’s multiple comparisons test to yield ****, ***, and ** indicating p < 0.0001, 0.001, and 0.01, respectively. n = 3 organoids/group. Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.s010
Immunofluorescence of ( A ) γH2AX in D35 cortical organoids and ( D ) CC3 in D28 cortical organoids. Bar, 100 μm. CC3 quantification by FACS of ( B ) D21 and ( C ) D28 cortical organoids. For each datapoint, 10–12 organoids from each line were analyzed via 3 technical replicates, and data from 4 mutant lines were consolidated to achieve rigorous comparisons. **, p < 0.01 and ns, not significant by two-way ANOVA test. ( E ) CC3 quantification of immunofluorescence images of D28 cortical organoids. For each line, 4–6 images and >10,000 cells were analyzed. *, p < 0.05; **, p < 0.01; ns, not significant by two-way ANOVA test. Graphs in ( B , C , E ) are presented in ratios (out of 1). Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.s011
Quantification of the ( A ) number and ( B ) surface area of ZO-1-positive ventricles. ( C ) Immunofluorescence of ZO-1 in D28 cortical organoids. Bar, 100 μm. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ns, not significant by two-way ANOVA test. Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.s012
( A ) Bright-field images of cortical organoids formed by cell lines 53BP1-S25A 34–3, 34–4, 79–1, 79–3 and S25D 14–3, 14–15, 14–19, 17, and 2 WT control at day 55 of differentiation. Bar, 1.5 mm. Blue transparent structures around organoids are Matrigel embedment. ( B ) At day 55 of differentiation, the size of cortical organoids was compared between groups. ( C ) The growth (comparing organoids at days 35 and 55) of cortical organoids were compared between groups. Data points represent single organoids. The mean ± SEM values were compared by one-way ANOVA with Dunnett’s multiple comparisons test to yield **** and ** indicating p < 0.0001 and 0.01, respectively. n = 15–36 organoids/group. ( D ) Two genes overlapped between up-regulated genes in 53BP1-S25A versus WT and down-regulated genes in 53BP1-S25D versus WT cortical organoids. No gene overlapped between down-regulated genes in 53BP1-S25A versus WT and up-regulated genes in 53BP1-S25D versus WT cortical organoids. ( E ) Down-regulated GSEA terms between 53BP1-S25A versus WT and 53BP1-S25D versus WT were not highly overlapped. Ten GSEA terms were specific to 53BP1-S25D versus WT. Underlying numerical values for figures are found in S1 Data .
https://doi.org/10.1371/journal.pbio.3002760.s013
( A ) Immunofluorescence of NPC markers PAX6 and NESTIN. Bar, 50 μm. GSEA identified top enrichment of differentially expressed genes in ( B , D ) 53BP1-S25A or ( C , E ) S25D versus WT NPCs. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. Venn diagrams depict overlaps between down-regulated genes in ATM -KO with 53BP1- ( F ) S25A or ( G ) S25D cortical organoids. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; GSEA, gene set enrichment analysis; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s014
( A ) GSEA graphs showed that up-regulated genes in 53BP1-S25A or S25D vs. WT had significant enrichment in down-regulated genes of ATM-KO vs. WT cortical organoids. P values were calculated by the hypergeometric test, assuming normal data distribution. ( B ) Concordantly differential expression of genes in 53BP1-S25D vs. WT were enriched in those in 53BP1-S25A vs. WT. ( C ) Concordantly differential expression of genes in 53BP1-S25A vs. WT were enriched in those in 53BP1-S25D vs. WT. For ( A - C ), P values were calculated by the hypergeometric test, assuming normal data distribution. ( D ) Proportions of 53BP1 binding to genomic features. 53BP1 ChIP-seq tracks at loci of representative ( E ) up-regulated and ( F ) down-regulated genes in 53BP1-S25A and S25D versus WT D35 cortical organoids. ( G ) S25A and S25D down-regulate 53BP1 targets that are enriched in IRE1-mediated unfolded protein response, regulation of cellular response to stress, iron import into cells, and regulation of apoptosis. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; GSEA, gene set enrichment analysis; KO, knockout; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s015
( A ) MA plot displays 53BP1 ChIP-seq signals at genomic sites that are significantly different in ATM -KO vs. WT NPCs. Proportions of genomic features and gene ontology of genes with ( B ) higher or ( C ) lower 53BP1 binding in ATM -KO vs. WT NPCs. ( D ) Proportions of 53BP1-pS25 binding to genomic features. ( E ) GSEA identified top enrichment of genes occupied by 53BP1-pS25 in WT NPCs. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; GSEA, gene set enrichment analysis; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s016
( A ) Principal component analysis of top 3,000 most variable peaks in 53BP1 ChIP-seq of 53BP1-WT, S25A, and S25D NPCs. Two independent cell lines for each group were used for ChIP-seq. Proportions of genomic features in regions with significantly different 53BP1 ChIP-seq in ( B ) 53BP1-S25A vs. WT and ( C ) 53BP1-S25D vs. WT, using the criterion of FC>2 and p < 0.05. ( D ) Heatmaps aligning peaks with significantly different 53BP1 ChIP-seq in 53BP1-S25A vs. S25D. Control regions are those, after voom normalization, showed the least changes and served as semi-independent validation of differential ChIP-seq analysis. Bubble graphs present top enriched categories of genes that had significantly lower 53BP1 ChIP-seq in ( E ) 53BP1-S25A vs. WT and ( F ) 53BP1-S25D vs. WT. Underlying numerical values for figures are found in S1 Data . FC, fold-change; NPC, neural progenitor cell; WT, wild type.
https://doi.org/10.1371/journal.pbio.3002760.s017
( A ) Schematic diagram of neural specification of hESCs with HH (SB421542), TGFβ (dorsomorphin), and WNT (IWR1e and cyclopamine) signaling inhibitors. Nuclear extract was harvested on day 4 and day 10. ( B ) WB analysis of day 4 samples. ( C ) Quantification of day 4 WB. Data are presented as the mean ± SEM, and Student t test was performed for pairwise comparisons. n.s., *, and ** indicate not significant, p < 0.05, and p < 0.01, respectively. ( D ) WB analysis of day 10 samples. ( E ) Quantification of day 10 WB. Data are presented as the mean ± SEM, and Student t test was performed for pairwise comparisons. n.s., *, and ** indicate not significant, p < 0.05, and p < 0.01, respectively. Underlying numerical values for figures are found in S1 Data . ATM, ataxia telangiectasia mutated; hESC, human embryonic stem cell; WB, western blot.
https://doi.org/10.1371/journal.pbio.3002760.s018
( A ) Alignment of WT and RNF168 -KO mutation sequences in the RNF168 locus. Red indicates the gRNA sequences. ( B ) WB analysis of WT and RNF168 -KO hESCs. ( C ) RT-qPCR analysis showing that pluripotent genes in RNF168 -KO were expressed higher or the same as those in WT. RNF168 -KO did not reduce pluripotent gene expression. *, p < 0.05; ns, not significant by two-way ANOVA text. ( D ) Immunofluorescence showed similar expression of OCT4 and SSEA4 proteins in WT and RNF168 -KO hESCs. Bar, 100 μm. ( E ) Immunofluorescence of showed similar expression of PAX6 and NES in NPCs. WT and RNF168 -KO NPCs. Bar, 50 μm. Functional terms that are highly enriched in ( F ) up-regulated and ( G ) down-regulated genes in RNF168 -KO D35 cortical organoids. % Match, % of genes in the enriched term that overlap the differentially expressed genes or proteins. ( H ) Heatmaps aligning peaks with 53BP1-pS25 CUT&RUN signals that were gained, the same, or lost in RNF168 -KO vs. WT NPCs, using the criterion of FC>2 and p < 0.05. n = numbers of peaks. Regions with the same signals, are n = 899, which showed the least changes after voom normalization and served as semi-independent validation of differential ChIP-seq analysis. ( I ) Functional terms of 53BP1-pS25-bound genes in WT NPCs. % Match, % of genes in the enriched term that overlap the differentially bound genes. ( J ) Number of differentially expressed genes identified by comparison of RNF168 -KO vs. WT NPCs at p < 0.05. Of these genes, we list the numbers of 53BP1-pS25-bound targets and targets with higher or lower 53BP1-pS25 CUT&RUN signals in RNF168 -KO NPCs. Underlying numerical values for figures are found in S1 Data . FC, fold-change; hESC, human embryonic stem cell; KO, knockout; NES, normalized enrichment score; NPC, neural progenitor cell; RT-qPCR, quantitative reverse transcription PCR; WB, western blot; WT, wild type; 53BP1-pS25, 53BP1 phosphorylated at serine 25.
https://doi.org/10.1371/journal.pbio.3002760.s019
Typically normal karyotypes and 3 abnormalities are shown.
https://doi.org/10.1371/journal.pbio.3002760.s020
Data suggest that D35 ATM -KO cortical organoids specified to the forebrain lineage.
https://doi.org/10.1371/journal.pbio.3002760.s021
https://doi.org/10.1371/journal.pbio.3002760.s022
https://doi.org/10.1371/journal.pbio.3002760.s023
Data from WT and 53BP1 mutants are compared pairwise by using the two-sample t test. The sizes of organoids are significantly different between each comparison pair (all p < 0.05).
https://doi.org/10.1371/journal.pbio.3002760.s024
This table lists the calculation for different combinations of data and the descriptive statistics.
https://doi.org/10.1371/journal.pbio.3002760.s025
https://doi.org/10.1371/journal.pbio.3002760.s026
https://doi.org/10.1371/journal.pbio.3002760.s027
The authors thank A. Andersen and I. Chen for discussions and editing the manuscript; A. N. Kettenbach for advice; J. Houston and K. Lowe for FACS; P. Sinojia and E. Rivera-Peraza for preliminary experiments and data analyses. Sequencing was performed at the Harwell Center for Biotechnology, images were acquired at the Cell & Tissue Imaging Center, and karyotyping was analyzed by J. Wilbourne and V. Valentine at the Cytogenetics Core; all are supported by SJCRH and NCI P30 (CA021765).
Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Nature ( 2024 ) Cite this article
48 Altmetric
Metrics details
CDK1 has been known to be the sole cyclin-dependent kinase (CDK) partner of cyclin B1 to drive mitotic progression 1 . Here we demonstrate that CDK5 is active during mitosis and is necessary for maintaining mitotic fidelity. CDK5 is an atypical CDK owing to its high expression in post-mitotic neurons and activation by non-cyclin proteins p35 and p39 2 . Here, using independent chemical genetic approaches, we specifically abrogated CDK5 activity during mitosis, and observed mitotic defects, nuclear atypia and substantial alterations in the mitotic phosphoproteome. Notably, cyclin B1 is a mitotic co-factor of CDK5. Computational modelling, comparison with experimentally derived structures of CDK–cyclin complexes and validation with mutational analysis indicate that CDK5–cyclin B1 can form a functional complex. Disruption of the CDK5–cyclin B1 complex phenocopies CDK5 abrogation in mitosis. Together, our results demonstrate that cyclin B1 partners with both CDK5 and CDK1, and CDK5–cyclin B1 functions as a canonical CDK–cyclin complex to ensure mitotic fidelity.
This is a preview of subscription content, access via your institution
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
24,99 € / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
185,98 € per year
only 3,65 € per issue
Buy this article
Prices may be subject to local taxes which are calculated during checkout
Data availability.
All data supporting the findings of this study are available in the Article and its Supplementary Information . The LC–MS/MS proteomics data have been deposited to the ProteomeXchange Consortium 60 via the PRIDE 61 partner repository under dataset identifier PXD038386 . Correspondence regarding experiments and requests for materials should be addressed to the corresponding authors.
Wieser, S. & Pines, J. The biochemistry of mitosis. Cold Spring Harb. Perspect. Biol. 7 , a015776 (2015).
Article PubMed PubMed Central Google Scholar
Dhavan, R. & Tsai, L. H. A decade of CDK5. Nat. Rev. Mol. Cell Biol. 2 , 749–759 (2001).
Article CAS PubMed Google Scholar
Malumbres, M. Cyclin-dependent kinases. Genome Biol. 15 , 122 (2014).
Coverley, D., Laman, H. & Laskey, R. A. Distinct roles for cyclins E and A during DNA replication complex assembly and activation. Nat. Cell Biol. 4 , 523–528 (2002).
Desai, D., Wessling, H. C., Fisher, R. P. & Morgan, D. O. Effects of phosphorylation by CAK on cyclin binding by CDC2 and CDK2. Mol. Cell. Biol. 15 , 345–350 (1995).
Article CAS PubMed PubMed Central Google Scholar
Brown, N. R. et al. CDK1 structures reveal conserved and unique features of the essential cell cycle CDK. Nat. Commun. 6 , 6769 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Strauss, B. et al. Cyclin B1 is essential for mitosis in mouse embryos, and its nuclear export sets the time for mitosis. J. Cell Biol. 217 , 179–193 (2018).
Gavet, O. & Pines, J. Activation of cyclin B1-Cdk1 synchronizes events in the nucleus and the cytoplasm at mitosis. J. Cell Biol. 189 , 247–259 (2010).
Barbiero, M. et al. Cell cycle-dependent binding between cyclin B1 and Cdk1 revealed by time-resolved fluorescence correlation spectroscopy. Open Biol. 12 , 220057 (2022).
Pines, J. & Hunter, T. Isolation of a human cyclin cDNA: evidence for cyclin mRNA and protein regulation in the cell cycle and for interaction with p34cdc2. Cell 58 , 833–846 (1989).
Clute, P. & Pines, J. Temporal and spatial control of cyclin B1 destruction in metaphase. Nat. Cell Biol. 1 , 82–87 (1999).
Potapova, T. A. et al. The reversibility of mitotic exit in vertebrate cells. Nature 440 , 954–958 (2006).
Basu, S., Greenwood, J., Jones, A. W. & Nurse, P. Core control principles of the eukaryotic cell cycle. Nature 607 , 381–386 (2022).
Santamaria, D. et al. Cdk1 is sufficient to drive the mammalian cell cycle. Nature 448 , 811–815 (2007).
Article ADS CAS PubMed Google Scholar
Zheng, X. F. et al. A mitotic CDK5-PP4 phospho-signaling cascade primes 53BP1 for DNA repair in G1. Nat. Commun. 10 , 4252 (2019).
Article ADS PubMed PubMed Central Google Scholar
Fagerberg, L. et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol. Cell. Proteom. 13 , 397–406 (2014).
Article CAS Google Scholar
Pozo, K. & Bibb, J. A. The emerging role of Cdk5 in cancer. Trends Cancer 2 , 606–618 (2016).
Sharma, S. & Sicinski, P. A kinase of many talents: non-neuronal functions of CDK5 in development and disease. Open Biol. 10 , 190287 (2020).
Sun, K. H. et al. Novel genetic tools reveal Cdk5’s major role in Golgi fragmentation in Alzheimer’s disease. Mol. Biol. Cell 19 , 3052–3069 (2008).
Sharma, S. et al. Targeting the cyclin-dependent kinase 5 in metastatic melanoma. Proc. Natl Acad. Sci. USA 117 , 8001–8012 (2020).
Nabet, B. et al. The dTAG system for immediate and target-specific protein degradation. Nat. Chem. Biol. 14 , 431–441 (2018).
Simpson, L. M. et al. Target protein localization and its impact on PROTAC-mediated degradation. Cell Chem. Biol. 29 , 1482–1504 e1487 (2022).
Vassilev, L. T. et al. Selective small-molecule inhibitor reveals critical mitotic functions of human CDK1. Proc. Natl Acad. Sci. USA 103 , 10660–10665 (2006).
Janssen, A. F. J., Breusegem, S. Y. & Larrieu, D. Current methods and pipelines for image-based quantitation of nuclear shape and nuclear envelope abnormalities. Cells 11 , 347 (2022).
Thompson, S. L. & Compton, D. A. Chromosome missegregation in human cells arises through specific types of kinetochore-microtubule attachment errors. Proc. Natl Acad. Sci. USA 108 , 17974–17978 (2011).
Kline-Smith, S. L. & Walczak, C. E. Mitotic spindle assembly and chromosome segregation: refocusing on microtubule dynamics. Mol. Cell 15 , 317–327 (2004).
Prosser, S. L. & Pelletier, L. Mitotic spindle assembly in animal cells: a fine balancing act. Nat. Rev. Mol. Cell Biol. 18 , 187–201 (2017).
Zeng, X. et al. Pharmacologic inhibition of the anaphase-promoting complex induces a spindle checkpoint-dependent mitotic arrest in the absence of spindle damage. Cancer Cell 18 , 382–395 (2010).
Warren, J. D., Orr, B. & Compton, D. A. A comparative analysis of methods to measure kinetochore-microtubule attachment stability. Methods Cell. Biol. 158 , 91–116 (2020).
Gregan, J., Polakova, S., Zhang, L., Tolic-Norrelykke, I. M. & Cimini, D. Merotelic kinetochore attachment: causes and effects. Trends Cell Biol 21 , 374–381 (2011).
Etemad, B., Kuijt, T. E. & Kops, G. J. Kinetochore-microtubule attachment is sufficient to satisfy the human spindle assembly checkpoint. Nat. Commun. 6 , 8987 (2015).
Tauchman, E. C., Boehm, F. J. & DeLuca, J. G. Stable kinetochore-microtubule attachment is sufficient to silence the spindle assembly checkpoint in human cells. Nat. Commun. 6 , 10036 (2015).
Mitchison, T. & Kirschner, M. Microtubule assembly nucleated by isolated centrosomes. Nature 312 , 232–237 (1984).
Fourest-Lieuvin, A. et al. Microtubule regulation in mitosis: tubulin phosphorylation by the cyclin-dependent kinase Cdk1. Mol. Biol. Cell 17 , 1041–1050 (2006).
Ubersax, J. A. et al. Targets of the cyclin-dependent kinase Cdk1. Nature 425 , 859–864 (2003).
Yang, C. H., Lambie, E. J. & Snyder, M. NuMA: an unusually long coiled-coil related protein in the mammalian nucleus. J. Cell Biol. 116 , 1303–1317 (1992).
Yang, C. H. & Snyder, M. The nuclear-mitotic apparatus protein is important in the establishment and maintenance of the bipolar mitotic spindle apparatus. Mol. Biol. Cell 3 , 1259–1267 (1992).
Kotak, S., Busso, C. & Gonczy, P. NuMA phosphorylation by CDK1 couples mitotic progression with cortical dynein function. EMBO J. 32 , 2517–2529 (2013).
Kitagawa, M. et al. Cdk1 coordinates timely activation of MKlp2 kinesin with relocation of the chromosome passenger complex for cytokinesis. Cell Rep. 7 , 166–179 (2014).
Schrock, M. S. et al. MKLP2 functions in early mitosis to ensure proper chromosome congression. J. Cell Sci. 135 , jcs259560 (2022).
Sun, M. et al. NuMA regulates mitotic spindle assembly, structural dynamics and function via phase separation. Nat. Commun. 12 , 7157 (2021).
Chen, Q., Zhang, X., Jiang, Q., Clarke, P. R. & Zhang, C. Cyclin B1 is localized to unattached kinetochores and contributes to efficient microtubule attachment and proper chromosome alignment during mitosis. Cell Res. 18 , 268–280 (2008).
Kabeche, L. & Compton, D. A. Cyclin A regulates kinetochore microtubules to promote faithful chromosome segregation. Nature 502 , 110–113 (2013).
Hegarat, N. et al. Cyclin A triggers mitosis either via the Greatwall kinase pathway or cyclin B. EMBO J. 39 , e104419 (2020).
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596 , 583–589 (2021).
Wood, D. J. & Endicott, J. A. Structural insights into the functional diversity of the CDK-cyclin family. Open Biol. 8 , 180112 (2018).
Brown, N. R., Noble, M. E., Endicott, J. A. & Johnson, L. N. The structural basis for specificity of substrate and recruitment peptides for cyclin-dependent kinases. Nat. Cell Biol. 1 , 438–443 (1999).
Tarricone, C. et al. Structure and regulation of the CDK5-p25 nck5a complex. Mol. Cell 8 , 657–669 (2001).
Poon, R. Y., Lew, J. & Hunter, T. Identification of functional domains in the neuronal Cdk5 activator protein. J. Biol. Chem. 272 , 5703–5708 (1997).
Oppermann, F. S. et al. Large-scale proteomics analysis of the human kinome. Mol. Cell. Proteom. 8 , 1751–1764 (2009).
van den Heuvel, S. & Harlow, E. Distinct roles for cyclin-dependent kinases in cell cycle control. Science 262 , 2050–2054 (1993).
Article ADS PubMed Google Scholar
Nakatani, Y. & Ogryzko, V. Immunoaffinity purification of mammalian protein complexes. Methods Enzymol. 370 , 430–444 (2003).
Tyanova, S., Temu, T. & Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 11 , 2301–2319 (2016).
Tyanova, S. et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 13 , 731–740 (2016).
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43 , e47 (2015).
R Core Team. R: a language and environment for statistical computing (2021).
Wickham, H. ggplot2: elegant graphics for data analysis (2016).
Slowikowski, K. ggrepel: automatically position non-overlapping text labels with “ggplot2” (2018).
Wu, T. et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation 2 , 100141 (2021).
CAS PubMed PubMed Central Google Scholar
Deutsch, E. W. et al. The ProteomeXchange consortium in 2020: enabling ‘big data’ approaches in proteomics. Nucleic Acids Res. 48 , D1145–D1152 (2020).
CAS PubMed Google Scholar
Perez-Riverol, Y. et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 47 , D442–D450 (2019).
Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26 , 139–140 (2010).
Nagahara, H. et al. Transduction of full-length TAT fusion proteins into mammalian cells: TAT-p27Kip1 induces cell migration. Nat. Med. 4 , 1449–1452 (1998).
Mirdita, M. et al. ColabFold: making protein folding accessible to all. Nat. Methods 19 , 679–682 (2022).
Lu, C. et al. OPLS4: improving force field accuracy on challenging regimes of chemical space. J. Chem. Theory Comput. 17 , 4291–4300 (2021).
Obenauer, J. C., Cantley, L. C. & Yaffe, M. B. Scansite 2.0: proteome-wide prediction of cell signaling interactions using short sequence motifs. Nucleic Acids Res. 31 , 3635–3641 (2003).
Download references
We thank D. Pellman for comments on the manuscript; W. Michowski, S. Sharma, P. Sicinski, B. Nabet and N. Gray for the reagents; J. A. Tainer for providing access to software used for structural analysis; and S. Gerber for sharing unpublished results. D.C. is supported by grants R01 CA208244 and R01 CA264900, DOD Ovarian Cancer Award W81XWH-15-0564/OC140632, Tina’s Wish Foundation, Detect Me If You Can, a V Foundation Award, a Gray Foundation grant and the Claudia Adams Barr Program in Innovative Basic Cancer Research. A. Spektor would like to acknowledge support from K08 CA208008, the Burroughs Wellcome Fund Career Award for Medical Scientists, Saverin Breast Cancer Research Fund and the Claudia Adams Barr Program in Innovative Basic Cancer Research. X.-F.Z. was an American Cancer Society Fellow and is supported by the Breast and Gynecologic Cancer Innovation Award from Susan F. Smith Center for Women’s Cancers at Dana-Farber Cancer Institute. A. Syed is supported by the Claudia Adams Barr Program in Innovative Basic Cancer Research. B.T. was supported by the Polish National Agency for Academic Exchange (grant PPN/WAL/2019/1/00018) and by the Foundation for Polish Science (START Program). A.D.D is supported by NIH grant R01 HL52725. A.G.P. by National Cancer Institute grants U01CA214114 and U01CA271407, as well as a donation from the Aven Foundation; J.R.W. by National Cancer Institute grant R50CA211499; and K.S. by NIH awards 1R01-CA237660 and 1RF1NS124779.
Bartłomiej Tomasik
Present address: Department of Oncology and Radiotherapy, Medical University of Gdańsk, Faculty of Medicine, Gdańsk, Poland
These authors contributed equally: Xiao-Feng Zheng, Aniruddha Sarkar
Division of Radiation and Genome Stability, Department of Radiation Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
Xiao-Feng Zheng, Aniruddha Sarkar, Aleem Syed, Huy Nguyen, Bartłomiej Tomasik, Kaimeng Huang, Feng Li, Alan D. D’Andrea, Alexander Spektor & Dipanjan Chowdhury
Department of Chemistry and Purdue University Center for Cancer Research, Purdue University, West Lafayette, IN, USA
Humphrey Lotana & Kavita Shah
Translational Science and Therapeutics Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
Richard G. Ivey, Jacob J. Kennedy, Jeffrey R. Whiteaker & Amanda G. Paulovich
Department of Biostatistics and Translational Medicine, Medical University of Łódź, Łódź, Poland
Broad Institute of Harvard and MIT, Cambridge, MA, USA
Kaimeng Huang, Alan D. D’Andrea, Alexander Spektor & Dipanjan Chowdhury
Department of Biological Chemistry & Molecular Pharmacology, Harvard Medical School, Boston, MA, USA
Dipanjan Chowdhury
You can also search for this author in PubMed Google Scholar
X.-F.Z., A. Sarkar., A. Spektor. and D.C. conceived the project and designed the experiments. X.-F.Z. and A. Sarkar performed the majority of experiments and associated analyses except as listed below. H.L. expressed relevant proteins and conducted the kinase activity assays for CDK5–cyclin B1, CDK5–p35 and CDK5(S46) variant complexes under the guidance of K.S.; A. Syed performed structural modelling and analysis. R.G.I., J.J.K. and J.R.W. performed MS and analysis. B.T. and H.N. performed MS data analyses. K.H. provided guidance to screen CDK5(as) knocked-in clones and performed sequence analysis to confirm CDK5(as) knock-in. F.L. and A.D.D. provided reagents and discussion on CDK5 substrates analyses. X.-F.Z., A. Sarkar, A. Spektor and D.C. wrote the manuscript with inputs and edits from all of the other authors.
Correspondence to Alexander Spektor or Dipanjan Chowdhury .
Competing interests.
A.D.D. reports consulting for AstraZeneca, Bayer AG, Blacksmith/Lightstone Ventures, Bristol Myers Squibb, Cyteir Therapeutics, EMD Serono, Impact Therapeutics, PrimeFour Therapeutics, Pfizer, Tango Therapeutics and Zentalis Pharmaceuticals/Zeno Management; is an advisory board member for Cyteir and Impact Therapeutics; a stockholder in Cedilla Therapeutics, Cyteir, Impact Therapeutics and PrimeFour Therapeutics; and reports receiving commercial research grants from Bristol Myers Squibb, EMD Serono, Moderna and Tango Therapeutics. The other authors declare no competing interests.
Peer review information.
Nature thanks Yibing Shan and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data fig. 1 inhibition of cdk5 in analogue-sensitive (cdk5- as ) system..
a , Schematics depicting specific inhibition of the CDK5 analogue-sensitive ( as ) variant. Canonical ATP-analogue inhibitor (In, yellow) targets endogenous CDK5 (dark green) at its ATP-binding catalytic site nonspecifically since multiple kinases share structurally similar catalytic sites (left panel). The analogue-sensitive ( as , light green) phenylalanine-to-glycine (F80G) mutation confers a structural change adjacent to the catalytic site of CDK5 that does not impact its catalysis but accommodates the specific binding of a non-hydrolysable bulky orthogonal inhibitor 1NM-PP1(In*, orange). Introduction of 1NM-PP1 thus selectively inhibits CDK5- as variant (right panel). b , Immunoblots showing two clones (Cl 23 and Cl 50) of RPE-1 cells expressing FLAG-HA-CDK5- as in place of endogenous CDK5. Representative results are shown from three independent repeats. c , Proliferation curve of parental RPE-1 and RPE-1 CDK5- as cells. Data represent mean ± s.d. from three independent repeats. p -value was determined by Mann Whitney U test. d , Immunoblots showing immunoprecipitated CDK1-cyclin B1 complex or CDK5- as -cyclin B1 complex by the indicated antibody-coupled agarose, from nocodazole arrested RPE-1 CDK5- as cells with treated with or without 1NM-PP1 for inhibition of CDK5- as , from three independent replicate experiments. e , In-vitro kinase activity quantification of immunoprecipitated complex shown in d . Data represent mean ± s.d. from three independent experiments. p -values were determined by unpaired, two-tailed student’s t-test. f , Immunoblots of RPE-1 CDK5- as cells treated with either DMSO or 1NM-PP1 for 2 h prior to and upon release from RO-3306 and collected at 60 min following release. Cells were lysed and blotted with anti-bodies against indicated proteins (upper panel). Quantification of the relative intensity of PP4R3β phosphorylation at S840 in 1NM-PP1-treated CDK5- as cells compared to DMSO-treatment (lower panel). g , Experimental scheme for specific and temporal abrogation of CDK5 in RPE-1 CDK5- as cells. Data represent mean ± S.D from quadruplicate repeats. p -value was determined by one sample t and Wilcoxon test. h , Hoechst staining showing primary nuclei and micronuclei of RPE-1 CDK5- as with indicated treatment; scale bar is as indicated (left panel). Right, quantification of the percentage of cells with micronuclei after treatment. Data represent mean ± s.d. of three independent experiments from n = 2174 DMSO, n = 1788 1NM-PP1 where n is the number of cells. p- values were determined by unpaired, two-tailed student’s t-test. Scale bar is as indicated. Uncropped gel images are provided in Supplementary Fig. 1 .
a , Schematic depicting the dTAG-13-inducible protein degradation system. Compound dTAG-13 links protein fused with FKBP12 F36V domain (dTAG) to CRBN-DDB1-CUL4A E3 ligase complex, leading to CRBN-mediated degradation. b , Immunoblots showing two clones of RPE-1 cells that express dTAG -HA-CDK5 in place of endogenous CDK5 (Cl N1 and Cl N4). Representative results are shown from three independent repeats. c , Proliferation curve of parental RPE-1 and RPE-1 CDK5-dTAG. Data represent mean ± s.d. of three independent repeats. p -value was determined by Mann Whitney U test. d and e , Representative images of RPE-1 CDK5- dTAG clone 1 (N1) ( d ) and RPE-1 CDK5- dTAG clone 4 (N4) ( e ) treated with DMSO or dTAG-13 for 2 h prior to and upon release from G2/M arrest and fixed at 120 min after release (top panel); quantification of CDK5 total intensity per cell (lower panels). Data represent mean ± s.d. of at least two independent experiments from n = 100 cells each condition. p- values were determined by unpaired, two-tailed student’s t-test. f , Immunoblots showing level of indicated proteins in RPE-1 CDK5- dTAG cells. Cells were treated with either DMSO or dTAG-13 for 2 h prior to and upon release from RO-3306 and lysed at 60 min following release (upper panel). Quantification of the relative intensity of PP4R3β phosphorylation at S840 in dTAG13-treated CDK5- dTAG cells compared to DMSO-treatment (lower panel). Data represent mean ± s.d. of four independent experiments. p -value was determined by one sample t and Wilcoxon test. g , Experimental scheme for specific and temporal abrogation of CDK5 in RPE-1 CDK5- dTAG cells. h , Hoechst staining showing primary nuclei and micronuclei of RPE-1 CDK5- dTAG with indicated treatment; scale bar is as indicated (left panel). Right, quantification of the percentage of cells with micronuclei after treatment. Data represent mean ± s.d. of three independent experiments from n = 2094 DMSO and n = 2095 dTAG-13, where n is the number of cells. p- values were determined by unpaired, two-tailed student’s t-test. Scale bar is as indicated. Uncropped gel images are provided in Supplementary Fig. 1 .
a and b , Live-cell imaging snapshots of RPE-1 CDK5- as cells ( a ) and RPE-1 CDK5- dTAG cells ( b ) expressing mCherry-H2B and GFP-α-tubulin, abrogated of CDK5 by treatment with 1NM-PP1 or dTAG-13, respectively. Imaging commenced in prophase following release from RO-3306 into fresh media containing indicated chemicals (left); quantification of the percentage of cells with abnormal nuclear morphology (right). c and d , Representative snapshots of the final frame prior to metaphase-to-anaphase transition from a live-cell imaging experiment detailing chromosome alignment at the metaphase plate of RPE- CDK5- as (c) and RPE-1 CDK5- dTAG ( d ) expressing mCherry-H2B, and GFP-α-tubulin (left); quantification of the percentage of cells displaying abnormal chromosome alignment following indicated treatments (top right). e , Representative images showing the range of depolymerization outcomes (low polymers, high polymers and spindle-like) in DMSO- and 1NM-PP1-treated cells, as shown in Fig. 2e , from n = 50 for each condition, where n is number of metaphase cells . f , Quantifications of mitotic duration from nuclear envelope breakdown (NEBD) to anaphase onset of RPE-1 CDK5- as (left ) and RPE-1 CDK5- dTAG (right) cells, following the indicated treatments. Live-cell imaging of RPE-1 CDK5- as and RPE-1 CDK5- dTAG cells expressing mCherry-H2B and GFP-BAF commenced following release from RO-3306 arrest into fresh media containing DMSO or 1NM-PP1 or dTAG-13. g , Quantifications of the percentage of RPE-1 CDK5- as (left) and RPE-1 CDK5- dTAG (right) cells that were arrested in mitosis following the indicated treatments. Imaging commenced in prophase cells as described in a , following release from RO-3306 into fresh media in the presence or absence nocodazole as indicated. The data in a, c , and g represent mean ± s.d. of at least two independent experiments from n = 85 DMSO and n = 78 1NM-PP1 in a and c ; from n = 40 cells for each treatment condition in g . The data in b , d , and f represent mean ± s.d. of three independent experiments from n = 57 DMSO and n = 64 dTAG-13 in b and d ; from n = 78 DMSO and n = 64 1NM-PP1; n = 59 DMSO and n = 60 dTAG-13, in f , where n is the number of cells. p- values were determined by unpaired, two-tailed student’s t-test. Scale bar is as indicated.
a, b , Immunostaining of RPE-1 cells with antibodies against CDK1 and α-tubulin ( a ); and CDK5 and α-tubulin ( b ) at indicated stages of mitosis. c, d , Manders’ overlap coefficient M1 (CDK1 versus CDK5 on α-tubulin) ( c ); and M2 (α-tubulin on CDK1 versus CDK5) ( d ) at indicated phases of mitosis in cells shown in a and b . The data represent mean ± s.d. of at least two independent experiments from n = 25 cells in each mitotic stage. p- values were determined by unpaired, two-tailed student’s t-test.
a , Scheme of cell synchronization for phosphoproteomics: RPE-1 CDK5- as cells were arrested at G2/M by treatment with RO-3306 for 16 h. The cells were treated with 1NM-PP1 to initiate CDK5 inhibition. 2 h post-treatment, cells were released from G2/M arrest into fresh media with or without 1NM-PP1 to proceed through mitosis with or without continuing inhibition of CDK5. Cells were collected at 60 min post-release from RO-3306 for lysis. b , Schematic for phosphoproteomics-based identification of putative CDK5 substrates. c , Gene ontology analysis of proteins harbouring CDK5 inhibition-induced up-regulated phosphosites. d , Table indicating phospho-site of proteins that are down-regulated as result of CDK5 inhibition. e , Table indicating the likely kinases to phosphorylate the indicated phosphosites of the protein, as predicted by Scansite 4 66 . Divergent score denotes the extent by which phosphosite diverge from known kinase substrate recognition motif, hence higher divergent score indicating the corresponding kinase is less likely the kinase to phosphorylate the phosphosite.
a , Endogenous CDK5 was immunoprecipitated from RPE-1 cells collected at time points corresponding to the indicated cell cycle stage. Cell lysate input and elution of immunoprecipitation were immunoblotted by antibodies against the indicated proteins. RPE-1 cells were synchronized to G2 by RO-3306 treatment for 16 h and to prometaphase (M) by nocodazole treatment for 6 h. Asynch: Asynchronous. Uncropped gel images are provided in Supplementary Fig. 1 . b , Immunostaining of RPE-1 cells with antibodies against the indicated proteins at indicated mitotic stages (upper panels). Manders’ overlap coefficient M1 (Cyclin B1 on CDK1) and M2 (CDK1 on Cyclin B1) at indicated mitotic stages for in cells shown in b (lower panels). The data represent mean ± s.d. of at least two independent experiments from n = 25 mitotic cells in each mitotic stage. p- values were determined by unpaired, two-tailed student’s t-test. c , Table listing common proteins as putative targets of CDK5, uncovered from the phosphoproteomics anlaysis of down-regulated phosphoproteins upon CDK5 inhibition (Fig. 3 and Supplementary Table 1 ), and those of cyclin B1, uncovered from phosphoproteomics analysis of down-regulated phospho-proteins upon cyclin B1 degradation (Fig. 6 and Table EV2 in Hegarat et al. EMBO J. 2020). Proteins relevant to mitotic functions are highlighted in red.
a , Predicted alignment error (PAE) plots of the top five AlphaFold2 (AF2)-predicted models of CDK5-cyclin B1 (top row) and CDK1-cyclin B1 (bottom row) complexes, ranked by interface-predicted template (iPTM) scores. b , AlphaFold2-Multimer-predicted structure of the CDK5-cyclin B1 complex. c , Structural comparison of CDK-cyclin complexes. Left most panel: Structural-overlay of AF2 model of CDK5-cyclin B1 and crystal structure of phospho-CDK2-cyclin A3-substrate complex (PDB ID: 1QMZ ). The zoomed-in view of the activation loops of CDK5 and CDK2 is shown in the inset. V163 (in CDK5), V164 (in CDK2) and Proline at +1 position in the substrates are indicated with arrows. Middle panel: Structural-overlay of AF2 model of CDK5-cyclin B1 and crystal structure of CDK1-cyclin B1-Cks2 complex (PDB ID: 4YC3 ). The zoomed-in view of the activation loops of CDK5 and CDK1 is shown in the inset. Cks2 has been removed from the structure for clarity. Right most panel: structural-overlay of AF2 models of CDK5-cyclin B1 and CDK1-cyclin B1 complex. The zoomed view of the activation loops of CDK5 and CDK1 is shown in the inset. d , Secondary structure elements of CDK5, cyclin B1 and p25. The protein sequences, labelled based on the structural models, are generated by PSPript for CDK5 (AF2 model) ( i ), cyclin B1 (AF2 model) ( ii ) and p25 (PDB ID: 3O0G ) ( iii ). Structural elements ( α , β , η ) are defined by default settings in the program. Key loops highlighted in Fig. 4d are mapped onto the corresponding sequence.
a , Structure of the CDK5-p25 complex (PDB ID: 1h41 ). CDK5 (blue) interacts with p25 (yellow). Serine 159 (S159, magenta) is in the T-loop. b , Sequence alignment of CDK5 and CDK1 shows that S159 in CDK5 is the analogous phosphosite as that of T161 in CDK1 for T-loop activation. Sequence alignment was performed by CLC Sequence Viewer ( https://www.qiagenbioinformatics.com/products/clc-sequence-viewer/ ). c , Immunoblots of indicated proteins in nocodazole-arrested mitotic (M) and asynchronous (Asy) HeLa cell lysate. d , Myc-His-tagged CDK5 S159 variants expressed in RPE-1 CDK5- as cells were immunoprecipitated from nocodazole-arrested mitotic lysate by Myc-agarose. Input from cell lysate and elution from immunoprecipitation were immunoblotted with antibodies against indicated protein. EV= empty vector. In vitro kinase activity assay of the indicated immunoprecipitated complex shown on the right panel. Data represent mean ± s.d. of four independent experiments. p -values were determined by unpaired two-tailed student’s t-test. e , Immunoblots showing RPE-1 FLAG-CDK5- as cells stably expressing Myc-His-tagged CDK5 WT and S159A, which were used in live-cell imaging and immunofluorescence experiments to characterize chromosome alignment and spindle architecture during mitosis, following inhibition of CDK5- as by 1NM-PP1, such that only the Myc-His-tagged CDK5 WT and S159A are not inhibited. Representative results are shown from three independent repeats. f , Hoechst staining showing nuclear morphology of RPE-1 CDK5- as cells expressing indicated CDK5 S159 variants following treatment with either DMSO or 1NMP-PP1 and fixation at 120 min post-release from RO-3306-induced arrest (upper panel); quantification of nuclear circularity and solidity (lower panels) g , Snapshots of live-cell imaging RPE-1 CDK5- as cells expressing indicated CDK5 S159 variant, mCherry-H2B, and GFP-α-tubulin, after release from RO-3306-induced arrest at G2/M, treated with 1NM-PP1 2 h prior to and upon after release from G2/M arrest (upper panel); quantification of cells displaying abnormal chromosome alignment in (lower panel). Representative images are shown from two independent experiments, n = 30 cells each cell line. h , Representative images of RPE-1 CDK5- as cells expressing indicated CDK5 S159 variants in metaphase, treated with DMSO or 1NM-PP1 for 2 h prior to and upon release from RO-3306-induced arrest, and then released into media containing 20 µM proTAME for 2 h, fixed and stained with tubulin and DAPI (upper panel); metaphase plate width and spindle length measurements for these representative cells were shown in the table on right; quantification of metaphase plate width and spindle length following the indicated treatments (lower panel). Data in f and h represent mean ± s.d. of at least two independent experiments from n = 486 WT, n = 561 S159A, and n = 401 EV, where n is the number of cells in f ; from n = 65 WT, n = 64 S159A, and n = 67 EV, where n is the number of cells in h . Scale bar is as indicated. Uncropped gel images are provided in Supplementary Fig. 1 .
a , Structure of the CDK5-p25 complex (PDB ID: 1h41 ). CDK5 (blue) interacts with p25 (yellow) at the PSSALRE helix (green). Serine 46 (S46, red) is in the PSSALRE helix. Serine 159 (S159, magenta) is in the T-loop. b , Sequence alignment of CDK5 and CDK1 shows that S46 is conserved in CDK1 and CDK5. Sequence alignment was performed by CLC Sequence Viewer ( https://www.qiagenbioinformatics.com/products/clc-sequence-viewer/ ). c , Immunoblots of CDK5 immunoprecipitation from lysate of E. coli BL21 (DE3) expressing His-tagged human CDK5 WT or CDK5 S46D, mixed with lysate of E. coli BL21 (DE3) expressing His-tagged human cyclin B1. Immunoprecipitated CDK5 alone or in the indicated complex were used in kinase activity assay, shown in Fig. 5b . Representative results are shown from three independent repeats. d , Immunoblots showing RPE-1 FLAG-CDK5- as cells stably expressing Myc-His-tagged CDK5 S46 phospho-variants, which were used in live-cell imaging and immunofluorescence experiments to characterize chromosome alignment and spindle architecture during mitosis, following inhibition of CDK5- as by 1NM-PP1, such that only the Myc-His-tagged CDK5 S46 phospho-variants are not inhibited. Representative results are shown from three independent repeats. e , Immunostaining of RPE-1 CDK5- as cells expressing Myc-His-tagged CDK5 WT or S46D with anti-PP4R3β S840 (pS840) antibody following indicated treatment (DMSO vs 1NM-PP1). Scale bar is as indicated (left). Normalized intensity level of PP4R3β S840 phosphorylation (right). Data represent mean ± s.d. of at least two independent experiments from n = 40 WT and n = 55 S46D, where n is the number of metaphase cells. p- values were determined by unpaired two-tailed student’s t-test. f , Immunoblots showing level of indicated proteins in RPE-1 CDK5- as cells expressing Myc-His-tagged CDK5 WT or S46D. Cells were treated with either DMSO or 1NM-PP1 for 2 h prior to and upon release from RO-3306 and collected and lysed at 60 min following release (left). Quantification of the intensity of PP4R3β phosphorylation at S840 (right). Data represent mean ± s.d. of four independent experiments. p -values were determined by two-tailed one sample t and Wilcoxon test. g , Representative snapshots of live-cell imaging of RPE-1 CDK5- as cells harbouring indicated CDK5 S46 variants expressing mCherry-H2B and GFP-α-tubulin, treated with 1NM-PP1, as shown in Fig. 5d , from n = 35 cells. Imaging commenced in prophase following release from RO-3306 into fresh media containing indicated chemicals. Uncropped gel images are provided in Supplementary Fig. 1 .
Immunostaining of RPE-1 CDK5- as cells stably expressing Myc-His CDK5-WT ( a ), S46A ( b ), and S46D ( c ) with antibodies against indicated protein in prophase, prometaphase, and metaphase. Data represent at least two independent experiments from n = 25 cells of each condition in each mitotic stage.
a , Chromatogram showing RPE-1 that harbours the homozygous CDK5- as mutation F80G introduced by CRISPR-mediated knock-in (lower panel), replacing endogenous WT CDK5 (upper panel). b , Immunoblots showing level of CDK5 expressed in parental RPE-1 and RPE-1 that harbours CDK5- as F80G mutation in place of endogenous CDK5. c , Representative images of CDK5- as knocked-in RPE-1 cells exhibiting lagging chromosomes following indicated treatments. d , Quantification of percentage of cells exhibiting lagging chromosomes following indicated treatments shown in (c). Data represent mean ± s.d. of three independent experiments from n = 252 DMSO, n = 220 1NM-PP1, where n is the number of cells. p -value was determined by two-tailed Mann Whitney U test.
a , CDK5 RNAseq expression in tumours (left) with matched normal tissues (right). The data are analysed using 22 TCGA projects. Note that CDK5 expression is higher in many cancers compared to the matched normal tissues. BLCA, urothelial bladder carcinoma; BRCA, breast invasive carcinoma; CESC cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; ESCA, esophageal carcinoma; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; STAD, stomach adenocarcinoma; THCA, thyroid carcinoma; THYM, thymoma; and UCEC, uterine corpus endometrial carcinoma. p -value was determined by two-sided Student’s t-test. ****: p <= 0.0001; ***: p <= 0.001; **: p <= 0.01; *: p <= 0.05; ns: not significant, p > 0.05. b , Scatter plots showing cells of indicated cancer types that are more dependent on CDK5 and less dependent on CDK1. Each dot represents a cancer cell line. The RNAi dependency data (in DEMETER2) for CDK5 and CDK1 were obtained from the Dependency Map ( depmap.org ). The slope line represents a simple linear regression analysis for the indicated cancer type. The four indicated cancer types (Head/Neck, Ovary, CNS/Brain, and Bowel) showed a trend of more negative CDK5 RNAi effect scores (indicative of more dependency) with increasing CDK1 RNAi effect scores (indicative of less dependency). The p value represents the significance of the correlation computed from a simple linear regression analysis of the data. Red circle highlights the subset of the cells that are relatively less dependent on CDK1 but more dependent on CDK5. c , Scatter plots showing bowel cancer cells that expresses CDK5 while being less dependent on CDK1. Each dot represents a cancer cell line. The data on gene effect of CDK1 CRISPR and CDK5 mRNA level were obtained from the Dependency Map ( depmap.org ). The slope line represents a simple linear regression analysis. Red circle highlights the subset of cells that are relatively less dependent on CDK1 but expresses higher level of CDK5. For b and c , solid line represents the best-fit line from simple linear regression using GraphPad Prism. Dashed lines represent 95% confidence bands of the best-fit line. p -value is determined by the F test testing the null hypothesis that the slope is zero. d , Scatter plots showing rapidly dividing cells of indicated cancer types that are more dependent on CDK5 and less dependent on CDK1. Each dots represents a cancer cell line. The doubling time data on the x-axis were obtained from the Cell Model Passports ( cellmodelpassports.sanger.ac.uk ). The RNAi dependency data (in DEMETER2) for CDK5, or CDK1, on the y-axis were obtained from the Dependency Map ( depmap.org ). Only cell lines with doubling time of less than 72 h are displayed and included for analysis. Each slope line represents a simple linear regression analysis for each cancer type. The indicated three cancer types were analysed and displayed because they showed a trend of faster proliferation rate (lower doubling time) with more negative CDK5 RNAi effect (more dependency) but increasing CDK1 RNAi effect (less dependency) scores. The p value represents the significance of the association of the three cancer types combined, computed from a multiple linear regression analysis of the combined data, using cancer type as a covariate. Red circle depicts subset of fast dividing cells that are relatively more dependent on CDK5 (left) and less dependent on CDK1 (right). Solid lines represent the best-fit lines from individual simple linear regressions using GraphPad Prism. p -value is for the test with the null hypothesis that the effect of the doubling time is zero from the multiple linear regression RNAi ~ Intercept + Doubling Time (hours) + Lineage.
Supplementary figure 1.
Full scanned images of all western blots.
Peer review file, supplementary table 1.
Phosphosite changes in 1NM-PP1-treated cells versus DMSO-treated controls as measured by LC–MS/MS.
Global protein changes in 1NM-PP1-treated cells versus DMSO-treated controls as measured by LC–MS/MS.
RPE-1 CDK5(as) cell after DMSO treatment, ×100 imaging.
RPE-1 CDK5(as) cell after 1NM-PP1 treatment (example 1), ×100 imaging.
RPE-1 CDK5(as) cell after 1NM-PP1 treatment (example 2), ×100 imaging.
RPE-1 CDK5(dTAG) cell after DMSO treatment, ×100 imaging.
RPE-1 CDK5(dTAG) cell after dTAG-13 treatment (example 1), ×100 imaging.
RPE-1 CDK5(dTAG) cell after dTAG-13 treatment (example 2) ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(WT) after 1NM-PP1 treatment, ×20 imaging.
RPE-1 CDK5(as) cells expressing MYC-EV after 1NM-PP1 treatment, ×20 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S159A) after 1NM-PP1 treatment (example 1), ×20 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S159A) after 1NM-PP1 treatment (example 2), ×20 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(WT) after 1NM-PP1 treatment, ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S46A) after 1NM-PP1 treatment (example 1), ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S46A) after 1NM-PP1 treatment (example 2), ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S46D) after 1NM-PP1 treatment (example 1), ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-CDK5(S46D) after 1NM-PP1 treatment (example 2), ×100 imaging.
RPE-1 CDK5(as) cells expressing MYC-EV after 1NM-PP1 treatment,×100 imaging.
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Reprints and permissions
Cite this article.
Zheng, XF., Sarkar, A., Lotana, H. et al. CDK5–cyclin B1 regulates mitotic fidelity. Nature (2024). https://doi.org/10.1038/s41586-024-07888-x
Download citation
Received : 24 March 2023
Accepted : 30 July 2024
Published : 04 September 2024
DOI : https://doi.org/10.1038/s41586-024-07888-x
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.
IMAGES
VIDEO
COMMENTS
Previous experiments indicate that CDK5 is active only when attached to a protein called p35. Which of the following best predicts how p35 might play a role in regulating neuron function? (A) Elevated intracellular levels of p35 result in increased synaptic activity. (B) Degradation of p35 results in increased synaptic activity. ...
Discovery of Cdk5. Cyclin-dependent kinase 5 (Cdk5) is a proline-directed serine/threonine protein kinase. Thirty years ago, Cdk5 was first discovered on the basis of its close sequence homology to the human cell division cycle protein 2 (Cdc2, also known as Cdk1), a regulator of cell cycle progression [1-3].The kinase activity of Cdc2 is detected in proliferating cells and activates prior ...
Cdk5 is a protein kinase that regulates various cellular processes in post-mitotic neurons and non-neuronal cells. Learn about its discovery, regulation, functions, and involvement in neurodegenerative diseases in this review article.
p35. p35 was the first Cdk5 activator to be identified (Fig. 1).It forms a tertiary structure that is similar to the cyclin-box fold domain that is required for Cdk activation, although the primary sequence of p35 is distinct from that of cyclins (Tang et al., 1997).In fact, the Cdk5-p35 complex has been shown to adopt a ternary structure that is similar to the complex between Cdk2 and ...
Gene-targeting experiments demonstrate an essential role of CDK5 in the cytoarchitecture of the CNS. ... These observations indicate that CDK5 is indispensable in neuregulin-mediated development ...
Abstract. The atypical cyclin-dependent kinase 5 (CDK5) is considered as a neuron-specific kinase that plays important roles in many cellular functions including cell motility and survival. The activation of CDK5 is dependent on interaction with its activator p35, p39, or p25. These activators share a CDK5-binding domain and form a tertiary ...
The atypical cyclin-dependent kinase 5 (CDK5) is considered as a neuron-specific kinase that plays important roles in many cellular functions including cell motility and survival. The activation of CDK5 is dependent on interaction with its activator p35, p39, or p25. These activators share a CDK5-binding domain and form a tertiary structure ...
Since Cdk5 activity is tightly regulated, a method for measuring its kinase activity is needed to fully understand the precise role of Cdk5 in developmental and disease processes. This article includes methods for detecting Cdk5 kinase activity in cultured cells or tissues, identifying new substrates, and screening for new kinase inhibitors.
Our previous studies indicated that Cdk5 activity is ... These data indicate that p25/Cdk5 may be more potent than p35 ... We also found in COS-7 cell transfection experiments that many cells died ...
These data indicate that NR2B, Cdk5, p35 and calpain are in a postsynaptic complex and Cdk5 mediates the interactions between NR2B and p35. NR2B-PSD-95 interaction ratios were unaffected by Cdk5 ...
p35 is a protein that binds and activates cyclin-dependent kinase 5 (Cdk5), which is essential for neuronal development and function. The authors show that p35 is a short-lived protein that is stabilized by Cdk5 activation and phosphorylation, and that it is degraded by the ubiquitin-proteasome pathway.
These data indicate that CDK5 is a substrate of presynaptic homeostatic plasticity and that changes in recycling versus resting pools driven by silencing can be accounted for by changes in presynaptic CDK5 levels. ... independent of CDK5 itself. These previous experiments however had not examined the impact of these inhibitors in the absence of ...
These results indicate that Cdk5 is localized in both the cytoplasm and nucleus of neurons. ... (mean ± SEM from three independent experiments; ***p < 0.005, cdk5 −/ ... Because our previous mass spectrometry screening identified MeCP2 as a potential substrate of Cdk5 and because the activity of nuclear Cdk5 was elevated by membrane ...
These data indicate that CDK5 is a substrate of presynaptic homeostatic plasticity and that changes in recycling versus resting pools driven by silencing can be accounted for by changes in presynaptic CDK5 levels. ... independent of CDK5 itself. These previous experiments however have not examined the impact of these inhibitors in the absence ...
The previous findings show a critical role of Cdk5 catalytic activity in 5-HT 6 R-elicited neurite extension and that the growth-promoting effect of Cdk5 depends on 5HT 6 R phosphorylation at ...
The signal for CDK5/p35 alone indicates CDK5 auto-phosphorylation seen in all lanes when CDK5 is present. (B) Annotated mass spectrum of the tryptic peptide PER2 ... a proteasome inhibitor, or with the solvent DMSO. In line with our previous experiments, shCdk5 treatment efficiently knocked down CDK5 and reduced PER2 levels compared with ...
Cdk5 assay was performed using Dynabeads™ with anti-Cdk5 antibody (~2 μg) from different companies. Cdk5 was immunoprecipitated from wild-type mouse brains. Lane 1 is C-8 (SantaCruz Biotechnology #sc-173), which has now been discontinued. The commercially available antibodies in Lanes 2 and 3 were not efficient in pulling down Cdk5.
A set of flashcards for AP Biology students to review cell regulation, neurotransmitter release, and cell cycle topics. The answer to the query question is option D, which explains how CDK5 inhibits neurotransmitter release.
Previous studies of Cdk5 in the mitochondria have mainly focused on neuronal cells ... Values are means ± SEM from three independent experiments. * indicates statistically significant difference ...
This web page contains multiple-choice questions (MCQ) on various topics in biology, such as neurotransmitter release, chemosynthesis, and evolution. The answer to the query is option D, which states that inhibition of CDK5 activity in neurons increases the movement of synaptic vesicles to the plasma membrane in response to a specific stimulus.
An important function of CDK5, especially in neurons, is the organization of the cytoskeleton and support of cellular outgrowths (Figure 2). Expression of p35 or p39 in vitro stimulates neurite outgrowths, and a dominant negative mutant of CDK5 was found to abolish the formation of these outgrowths [51].
ATM up-regulation in NPCs was shown by a previous DNA damage response study . Next, we used the CRISPR-Cas9 ... we found higher phosphorylation of proteins related to CDK5 activities, including ADD2, ADD3, DCX ... -change >2 and p < 0.05 to perform pairwise comparisons of the merged datasets from 53BP1-WT, S25A, and S25D ChIP-seq experiments ...
This hypothesis is consistent with previous data, ... recombinant Grb2-SH3 domain and then phosphorylated by Cdk5. Experiments in c and d were ... These experiments also indicate that dynamin ...
The experiments using the CDK5 dominant negative mutant demonstrated that CDK5 affected cytoskeletal protein F-actin ... Previous studies indicate that CDK5 regulates several processes in ...
d, Immunoblots showing RPE-1 FLAG-CDK5-as cells stably expressing Myc-His-tagged CDK5 S46 phospho-variants, which were used in live-cell imaging and immunofluorescence experiments to characterize ...