Mammalian circadian clock proteins form dynamic interacting microbodies distinct from phase separation

Abstract

Liquid-liquid phase separation (LLPS) underlies diverse biological processes. Because most LLPS studies were performed in vitro or in cells that overexpress protein, the physiological relevance of LLPS is unclear. PERIOD proteins are central mammalian circadian clock components and interact with other clock proteins in the core circadian negative feedback loop. Different core clock proteins were previously shown to form large complexes. Here we show that when transgene was stably expressed, PER2 formed nuclear phosphorylation-dependent LLPS condensates that recruited other clock proteins. Super-resolution microscopy of endogenous PER2, however, revealed formation of circadian-controlled, rapidly diffusing microbodies that were resistant to protein concentration changes, hexanediol treatment, and loss of phosphorylation, indicating that they are distinct from the LLPS condensates caused by overexpression. Surprisingly, only a small fraction of endogenous PER2 microbodies transiently interact with endogenous BMAL1 and CRY1, a conclusion that was confirmed in cells and in mice tissues, suggesting an enzyme-like mechanism in the circadian negative feedback process. Together, these results demonstrate that the dynamic interactions of core clock proteins is a key feature of mammalian circadian clock mechanism and the importance of examining endogenous proteins in LLPS and circadian studies.

Keywords: PER proteins; circadian clock; phase separation; phosphorylation; protein complex.

Figure 1.. LLPS behavior of PER2-EGFP when it is constitutively overexpressed in U2OS cells.

(A) Top: A diagram depicting human PER2 domains and phosphorylation sites. CKBD :Casein Kinase Binding domain. CBD: CRY binding domain. Location of the known phosphorylation sites were indicated based on a previously published data(12). Bottom: Disorder score across the human PER2 protein as predicted by IUPred. (B) Representative bioluminescence rhythms of the U2OS cells that stably express either EGFP or PER2-EGFP as well as an Per2(E2)-Luc reporter. Rhythms were monitored by LumiCycle and baseline-subtracted bioluminescence data are plotted (n=3 from three independent experiments). (C) Representative confocal images of cells stably express PER2-EGFP (high and low levels) and of control cells that express EGFP. Scale bar: 5 μm. (D) Images of a live cell that stably expresses PER2-EGFP over time. Two small condensates that fused during the time course are indicated with arrowheads. Scale bar: 5 μm. (E) Images of a live U2OS cell that stably expresses PER2-EGFP before and after photobleaching of regions marked by red squares. Scale bar: 5 μm. (F) Relative fluorescence intensity of PER2-EGFP condensates after photobleaching in the PER2-EGFP stable cells. Data are presented as mean ± SD (n=11). (G) Top: Images of live cells that express PER2-EGFP at 0 min and 10 min after treatment with 1,6-hexanediol (1.5%). Bottom: Quantification of mean fluorescence intensity of PER2-EGFP nuclei at 0 min and 10 min after treatment with 1,6-hexanediol. n= 8 cells. Data are presented as mean ± SD, unpaired two-tailed Student’s t test, **p < 0.01. (H) Representative confocal images of U2OS cells that overexpress PER2-EGFP, PER2(500–950)-EGFP, PER2(1–500)-EGFP or PER2(951–1255)-EGFP. PER2(1–500)-EGFP and PER2(951–1255)-EGFP protein are accumulated in the cytoplasm. (I) Representative confocal images cells that constitutively express PER2(L730G)-EGFP or PER2(527–818 S/T-A)-EGFP (PER2 with known and potential phosphorylation sites within aa 527–818 mutated to alanines). Scale bar: 5 μm. (J) Fluorescence immunostaining of endogenous BMAL1 in cells that stably overexpresses PER2-EGFP. Endogenous BMAL1 was stained with a BMAL1 antibody, and PER2-EGFP was monitored directly using the GFP channel. Scale bar: 5 μm.

Figure 2.. CRISPR/Cas9-mediated knock-in to tag the endogenous PER2 with EGFP.

(A) Schematic of PER2-EGFP knock-in donor plasmid and genome editing strategy. (B) Live cell confocal imaging of PER2-EGFP KI U2OS cells and U2OS cells without KI. Scale bar: 5 μm. (C) Lambda scan of PER2-EGFP KI cells and control cells. The GFP wavelength was specifically extracted and all other wavelengths were merged as background. Scale bar: 5 μm. (D) Representative bioluminescence rhythms of the Per2(E2)-Luc reporter in the PER2-EGFP KI and control U2OS cells. Baseline-subtracted bioluminescence data are plotted (n=3). (E) Western blot showing the PER2-EGFP levels in the PER2-EGFP KI cells and cells that stably overexpress PER2-EGFP. The protein sample from the cells that stably overexpress PER2-EGFP was diluted 20 fold.

Figure 3.. Super-resolution imaging of PER2-EGFP KI cells showing PER2 is concentrated in nuclear microbodies.

(A) Airyscan confocal imaging of live (left) and fixed (right) PER2-EGFP KI cells in the nucleus. Scale bar: 5 μm. (B) Airyscan imaging of fixed PER2-EGFP KI cells (left), PER2-EGFP KI cells treated with Per2-specific siRNA (si-PER2), and Per2^KO cells (Per2 gene disrupted by CRSPR/Cas9). Scale bar: 5 μm. (C) Time course 3D imaging results of nuclear GFP signals in PER2-EGFP KI cells after dexamethasone synchronization. Representative 2D images from the 3D results are shown. (D) Numbers of nuclear PER bodies at different time points after dexamethasone synchronization (n=8–17 cells). Data are presented as mean ± SD, n = （8–17） cells, unpaired two-tailed Student’s t test, **p < 0.01; ****p < 0.0001. (E) Quantification of the PER2-EGFP microbodies of time course 2D imaging results of the PER2-EGFP KI cells after dexamethasone synchronization from hr 24–72. (F) Quantification of PER2-EGFP microbodies of time course 2D imaging results of of the PER2-EGFP KI cells after temperature cycle synchronization from hr 0–48. (G) Top: Western blot of EGFP in PER2-EGFP KI cells at indicated time points after dexamethasone synchronization. Equal amount of protein was loaded in each lane. Bottom: Quantification of PER2-EGFP levels over time of the western blot result. (H) Nuclear PER body mean fluorescence intensities at indicated time points after dexamethasone synchronization (n=8–17 cells).

Figure 4.. Most PER2 bodies diffuse freely in nucleus and are formed by mechanism distinct from overexpression-caused LLPS condensates.

(A) Representative image of nuclear PER bodies in PER2-EGFP KI cell under light-sheet fluorescence microscopy at hr 10 after synchronization. Scale bar: 5 μm. (B) Representative movement tracks of PER bodies. Tracks 1–4: freely diffusing tracks. Tracks 5: immobile tracks. Track 6: a diffusing PER body became immobile. Scale bar: 5 μm. (C) The averaged mean square displacement (MSD) of PER bodies computed from data collected from 3 cells. Data are presented as means ± SEM. (D) Mean speeds of different tracks (n=3 cells). (E) Duration of immobility of PER bodies of the identified immobile PER bodies (data collected from 3 cells). (F) Representative image of nuclear PER bodies at 0 min and 10 min after treatment with 1,6-hexanediol in fixed KI cells at hr 10 after synchronization. Scale bar: 5 μm. (G) Comparison of nuclear PER2-EGFP fluorescence intensities of KI cells between 0 min and 10 min after treatment with 1,6-hexanediol at hr 10 after synchronization. Data are presented as median fluorescence intensity in each cells. n=14–18 cells. (H) Western blot of PER2-EGFP in the PER2-EGFP KI and the PER2(Δ727–731)-EGFP KI; Per1^KO cells at hr 10 after synchronization. The arrow indicates the phosphorylated PER2-EGFP species. * indicates a background protein band. (I) Representative Airyscan images of PER bodies in fixed PER2-EGFP KI and PER2(Δ727–731)-EGFP KI; Per1^KO cells at hr 10 after synchronization. Scale bar: 5 μm. (J) Nuclear PER2 body numbers in the PER2-EGFP KI, PER2-EGFP KI; Per1^KO, and PER2(Δ727–73)-EGFP KI; Per1^KO cells at 10 h and 16 h after dexamethasone synchronization.

Figure 5.. Airyscan imaging analyses of dual fluorescence reporter U2OS cells in which the endogenous PER2 and BMAL/CRY1 are tagged by EGFP and mScarlet-I, respectively.

(A) Diagrams showing the design of BMAL1-mScarlet-I and CRY1-mScarlet-I knock-in donor plasmids. (B) Time course Airyscan fluorescence imaging of PER2 and BMAL1 in fixed PER2-EGFP and BMAL1-mScarlet-I double KI cells after dexamethasone synchronization. Scale bar: 5 μm. (C) Time course Airyscan imaging of PER2 and CRY1 in fixed PER2-EGFP and CRY1-mScarlet-I double KI cells after dexamethasone synchronization. Scale bar: 5 μm. (D) Airyscan fluorescence imaging of PER2 and BMAL1 in live PER2-EGFP and BMAL1-mScarlet-I double KI cells. (E) Airyscan fluorescence imaging of PER2 and BMAL1 in live fixed PER2-EGFP and CRY1-mScarlet-I double KI cells.

Figure 6.. Airyscan imaging analyses of double KI cells showing that most PER2 bodies do not associate with BMAL1/CRY1 and the DNA residence time of BMAL1.

(A) 3D Airyscan imaging result of PER2-EGFP and BMAL1-mScarlet-I double KI cells at hr 10 after synchronization. Representative converted 2D nuclear images were shown. The 3D Airyscan imaging results for other time points are presented in the supplemental data files. Scale bar: 5 μm. (B) 3D imaging of PER2-EGFP and CRY1-mScarlet-I double KI cells at hr 10 after synchronization. Representative converted 2D nuclear images were shown. Scale bar: 5 μm. (C) Numbers of BMAL1 bodies and CRY1 bodies at different time points after dexamethasone synchronization in the double KI cells. Data are means ± SEM (n=8 cells). (D) Numbers of PER2 and BMAL/CRY1 foci that colocalize as a function of time after dexamethasone synchronization. Threshold for colocalization was set < 100 nm. Data are presented as means ± SE (n=8 cells). (E) Percentages of BMAL1 and PER2 foci that colocalize at different time points after synchronization in the PER2 and BMAL1 double KI cells. Data are presented as means ± SEM (n=8 cells). (F) Percentage of CRY1 and PER2 foci that colocalize at different time points after synchronization in the PER2 and CRY1 double KI cells. Data are presented as means ± SEM (n=8 cells). (G) Time series of live cell Airyscan imaging of BMAL1 foci in the BMAL1-mScarlet-I knock-in cells. A selected nuclear region was shown. The arrow indicate an immobile BMAL1 foci. Frame rate: 1.02 s/frame. Scale bar: 5 μm. (H) Summary of duration of immobility of identified immobile BMAL1 foci in the BMAL1-mScarlet-I KI cells. Data are presented as median (n=10 cells).

Figure 7.. PER, BMAL1 and CRY1 microbodies and their interactions revealed by STED microscopy.

(A) STED microscopy images of PER2 and BMAL1 puncta in the PER2-EGFP KI cells at 10 h after synchronization. PER2 and BMAL1 were stained by anti-GFP (Mouse) and anti-BMAL1 antibody (Rabbit), respectively. Scale bar: 5 μm. B: STED microscopy images of PER2 and CRY1 foci in the PER2-EGFP KI cells at 10 h after synchronization. PER2 and CRY1 were stained by anti-GFP and anti-CRY1 antibody, respectively. Scale bar: 5 μm. (C) Percentages of PER2 and BMAL1 bodies that colocalize at 10 h after synchronization. Distance threshold was set < 100 nm. Data are presented as mean ± SEM (n=8 cells). (D) Percentages of PER2 and CRY1 bodies that colocalize at 10 h after synchronization. Distance threshold was set < 100nm. Data are presented as mean ± SEM (n=3 cells). (E) Comparison of PER body sizes estimated by conventional anti-GFP antibody and anti-GFP nanobody. Data are presented as mean ± SEM (n=5 cells). (F) Immunodepletion assay of CRY1 or BMAL1 in extracts prepared from the indicated double KI U2OS cells 10 h after synchronization. After immunodepletion, western blot analyses were performed using the indicated antibodies. (G) Immunodepletion assays of mPER2 in extracts prepared from mice liver tissue (ZT18). (H) Representative immunofluorescence imaging results of mPER2 and mCRY1 in the nucleus using mice SCN tissues harvested at the indicated ZT time points using a rabbit polyclonal PER2 antibody and a mouse CRY1 antibody, respectively. Cells from the core region of SCN were selected here. Scale bar: 5 μm. (I) Quantification of the nuclear PER2/CRY1 foci numbers in the SCN cells from tissues harvested at the indicated time points in ZT or CT.