EMBO/SourceData huggingface dataset:
f2(a) Inhibition of ROCK1 activity ablates the interaction between Beclin1 and ROCK1. HeLa cells were untreated (H2O) or treated with 10 μM ROCK1 inhibitor, Y27632, for 8 h and then cultured in control media (+) or HBSS (−) for 4 h. Cells were collected, lysed and immunoprecipitated with Beclin1 (left panel) or ROCK1 (right panel) antibody. Left panel: immune complexes were divided into two (ROCK1-MOPS buffer and Beclin1-MES buffer) and probed for ROCK1 and Beclin1. Relative values of ROCK1 versus Beclin1 are shown.(b) An increase in ROCK1-mediated phosphorylation of MYPT1 upon starvation. HeLa and EJ cells were starved in HBSS for indicated time points and endogenous ROCK1 was immunoprecipitated. An in vitro kinase assay was performed using MYPT1 as substrate and immune-complexed ROCK1. Proteins were resolved and visualized with Coomassie staining (top), western blot analysis (middle) and autoradiography (bottom). Inputs were examined for ROCK1 and β-actin.(c) HeLa cells incubated in HBSS were collected and endogenous phosphorylation of MYPT1 was determined by western blotting. The same blot was then washed in TBST for 30 min, blocked in 5% milk and was analysed for total MYPT1 as a control. Extracts from this experiment were rerun on 7.5% gel and analysed for ROCK1 and β-actin.(d) ROCK activity is increased upon starvation. Cell lines (HeLa and EJ) grown in control (+Glu) or starved with HBSS for 1 and 3 h were lysed, extracts normalized and ROCK kinase activity was measured by performing enzyme-linked immunosorbent assay. Graphs represent ROCK activity in starved cells plotted as percent of ROCK1 activity in control (glucose-rich) cells. Data were obtained from triplicates performed at the same time (mean±s.d.). P-value was calculated using Student's t-test.(e) ROCK1 activity is not dependent on Beclin1. siControl (siCont.) and siBeclin1 knockdown cells were incubated in starvation media (HBSS) for indicated time. Resulting cell extracts were used to confirm Beclin1 knockdown and ROCK1 expression (left panel). Cell extracts from the same experiment were used for enzyme-linked immunosorbent assay measuring ROCK activity (right panel). Graph represents mean±s.d. of duplicate samples read at the same time (Student's t-test; NS, not significant).
f3(a) The effect of ROCK1 depletion on autophagy marker, LC3II. MEF ROCK1flox/flox cells were generated and treated with Con or Ad-Cre for 72 h. Cells were then incubated in nutrient-rich (+Glucose) or HBSS media for 2 h. Cell extracts were resolved by SDS-PAGE and ROCK1 KO was confirmed by western blotting (left panel). The same blot was used for LC3 and β-actin expression. Cells from the same experiment were fixed with chilled acetone and immunofluroscence was performed against endogenous LC3 (red). Blue represents nuclear DAPI staining. Scale bar, 20 μm.(b) EJ cells (left) with a stable knockdown using shCont. or shROCK1#1 were incubated in control (+glucose) or nutrient-free (HBSS) medium (−glucose) for 24 h and lysed. Lysates were subjected to western blotting with antibodies indicated. HeLa cells (right) with a transient knockdown using siCont. and siROCK1#2 for 30 h were further incubated in control or HBSS media for 8 h. Whole cell extracts were prepared, resolved by SDS-PAGE and analysed for ROCK1, LC3 and β-actin. Relative ratios of LC3II/LC3I are shown.(c) HeLa cells were treated with H2O or 10 μM ROCK1 inhibitor, Y27632, for 8 h. Cells were then cultured in high glucose (0 h) or starved for in HBSS nutrient-free media for 6 h. One hour before collection, cells were incubated in dimethylsulphoxide (control) or 0.1 μM bafilomycin-A1. Whole cell lysates were resolved by SDS-PAGE and were analysed by western blotting with indicated antibodies.(d) HeLa cells incubated with H2O or 10 μM Y27632 for 8 h were washed with PBS and cultured in control (+Glu) or HBSS media (−Glu) for an additional 8 h. Cells were fixed in acetone and immunofluorescence was performed against endogenous LC3 antibody (green). Scale bar, 20 μm. Graph represents % of LC3 punctae/cell (mean±s.d.), n=25 cells; (Student's t-test) **P=0.00041.(e) Transmission electron micrographs of EJ cells with stable shCont. or shROCK1#1 knockdown, cultured in control or HBSS medium for 16 h. Arrow bars indicate autophagic vacuoles. Scale bar, 500 nm.
f4(a) Inhibition of ROCK1 kinase activity prevents metabolic stress-mediated phosphorylation of Beclin1. HeLa cells left untreated or treated with a ROCK1 inhibitor, Y27632 (10 μM), were starved of glucose (4 h) as indicated, and endogenous Beclin1 was immunoprecipitated and blotted against phospho-Ser/Thr antibody. Whole cell lysates were run as inputs for Beclin1 and β-actin.(b) Phosphorylation of Beclin1 by ROCK1. Recombinant ROCK1 was used for in vitro kinase assay in the presence or absence of Y27632, using recombinant His-Beclin1 (left) or Flag-Beclin1 immunopurified from transfected HeLa cells (right). Proteins were resolved by SDS-PAGE; phosphorylated protein was visualized with autoradiography, and Beclin1 by Coomassie staining.(c) Identification of ROCK1-mediated phosphorylation domains in Beclin1. The indicated Flag-Beclin1 fragments were purified from HeLa cells and used as substrates for in vitro ROCK1 kinase assay. 32P-autoradiogram (centre) analysed phosphorylation and western blotting (WB; right) determined protein levels. Schematic representation of Beclin1 domain structure and deletion constructs are shown (left).(d) Identification of the phosphorylation site on Beclin1. Representation of point mutations in the different domains of Beclin1 are shown (left). In vitro kinase assay using Flag-Beclin1 WT and mutants (T38A and T119A), immunoprecipitated from transfected HeLa cells, as substrate and recombinant ROCK1 was performed. Phosphorylation was detected by 32P-autoradiogram, and Flag-Beclin1 levels were examined by WB.(e) Flag-Beclin1 Wt or T119A transfected HeLa cells were incubated in glucose-rich or nutrient-free (HBSS) media, in the presence or absence of ROCK1 inhibitor Y27632. Total cell extracts were used for IP using Flag agarose. Eluted protein was analysed by WB against phospho-T119 (Beclin) antibody and Beclin1. Input for Beclin1 was run on 7.5% gel and immunoblotted.
f5(a) Inhibiting ROCK1 activity increases the association between Beclin1 and Bcl-2. HeLa cells, untreated or treated with 10 μM Y27632 (top panel), were cultured in control or glucose-free (4 h) medium, lysed and cell extracts prepared. Endogenous Bcl-2 was immunoprecipitated with cross-linked Bcl-2 agarose and resulting eluted immune complexes were resolved by SDS-PAGE and blotted against indicated antibodies. EJ shCont. and shROCK1#1 cells (bottom panel) were grown in control or starvation media for 4 h and Bcl-2 IP was performed as above. Inputs were rerun to confirm ROCK1 knockdown.(b) HeLa cells were transfected with full-length constructs of Flag-Beclin WT, Flag-Beclin T119A or Flag-Beclin T119E. Flag-Beclin WT cells were then treated with Y27632 for 8 h. Post transfection and treatment with Y27632, cells were starved in HBSS medium for 4 h, after which cells were collected and whole cell lysates prepared. IP with Flag-agarose was performed for 2 h and proteins bound to beads were eluted and used for western blotting.(c) Flag-Beclin1 Wt or T119A transfected HeLa cells were incubated in glucose-rich or nutrient-free (HBSS) media, in the presence or absence of ROCK1 inhibitor Y27632. Total cell extracts were used for IP using Flag agarose. Eluted protein was analysed by western blotting against Bcl-2 and Beclin1. Inputs and IP were run on separate gels and same exposure is shown.(d) HeLa cells with a transient knockdown using siCont. and siBeclin1-3′-untranslated region for 30 h were further left untreated or treated with 10 μM Y27632. Cells were then incubated in control or HBSS media for 8 h. Whole cell extracts were prepared, resolved by SDS-PAGE and were analysed for Beclin1, LC3 and β-actin.(e) Endogenous LC3 staining in WT- or T119A-Beclin1-transfected cells. HeLa cells transfected with Flag-Beclin WT (±Y27632) and T119A were cultured in HBSS media for 6 h, fixed with cold acetone and endogenous LC3 immunofluorescence was performed. Representative images are shown; scale bar, 20 μm. Arrows indicate punctate LC3 staining. Graph represents mean±s.d., n=25, LC3 puncta/cell.
f6(b) PCR assay using tail genomic DNA from three different genotypes, separated on 0.7% agarose gel (top). Western blot analysis of ROCK1 expression in the heart tissue from three different genotypes using anti-ROCK1 antibody (bottom). β-Actin was used for protein loading control in each lane.(c) Mice were fed or starved for 16 h, and heart samples were prepared in lysis buffer. Resulting proteins were resolved by SDS-PAGE and blotted against ROCK1, LC3, autophagy marker p62 and β-actin.(d) ROCK1-deficient hearts or control hearts from animals fed or starved for 16 h were subjected to electron microscopy to identify autophagosome structures. Arrows indicate autophagic vacuoles. Scale bar, 500 nm.(e) Hearts from fed and starved ROCK1 WT and KO mice were sectioned, fixed in acetone and immunofluorescence assay was performed against endogenous LC3. Data are representative of two independent experiments from each genotype. Scale bar, 50 μm.
Figure 1(E) Confirmation of thalidomide dependency on six hit proteins using an in vitro binding assay. Interaction between bls-CRBN and FLAG-GST-protein in the presence of DMSO or 50 µM thalidomide was detected using the AlphaScreen technology.(F) In vitro binding assay for thalidomide, pomalidomide and lenalidomide. Interaction between bls-CRBN and FLAG-GST-PLZF in the presence of DMSO (3.125, 6.25, 12.5, 25, 50, 100 or 200 µM) thalidomide, pomalidomide or lenalidomide was analysed using the AlphaScreen technology. Data information: All relative AlphaScreen (AS) signals are expressed as relative luminescent signal with luminescent signal of DMSO as one. Error bars denote ± standard deviation (n=3).
Figure 2(E) In vitro binding and ubiquitination assay of AGIA-PLZF. Empty vector, AGIA-PLZF or FLAG-CRBN expressing HEK293T cells were lysed and the lysates were mixed. The first immunoprecipitation with anti-AGIA or anti-FLAG antibodies was performed in the presence of DMSO or 200 µM lenalidomide. The purified AGIA-PLZF or CRL4FLAG-CRBN complex, including AGIA-PLZF and FLAG-CRBN, was incubated with recombinant E1, E2 and HA-ubiquitin in the presence of DMSO or 200 µM lenalidomide. The second immunoprecipitation was performed using anti-AGIA antibody. Ubiquitination of PLZF was analysed by immunoblot.
Figure 3(B) In vitro binding assay using truncated PLZF. Thalidomide-dependent interaction between bls-CRBN and FLAG-GST-PLZF-full length (FL) or truncated FLAG-GST-PLZF was analysed in the presence of DMSO or 50 µM thalidomide using the AlphaScreen technology. Data information: All relative AlphaScreen (AS) signals are expressed as a relative luminescent signal with a luminescent signal of DMSO as one. Error bars denote ± standard deviation (n = 3). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, ****P < 0.0001
Figure 4(D) Immunoblot analysis of endogenous SALL4 or PLZF protein levels in HuH7 cells treated with DMSO, thalidomide or 5-hydroxythalidomide for 24 h. (E) Immunoblot analysis of endogenous PLZF or IKZF1 protein levels in THP-1 cells treated with DMSO, thalidomide or 5-hydroxythalidomide for 24 h. Data information: Data in bar graphs in (D) and (E) were calculated as relative PLZF band intensity with PLZF band intensity of DMSO as 100. Error bars donated ± standard deviation (n = 3) and P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, *P < 0.05, **P < 0.01 and ***P < 0.001).
Figure 5(A, B) Immunoblot analysis of AGIA-MmPLZF (A) or -MmSALL4 (B) in FLAG-MmCRBN-WT, -MmCRBN-I391V or -MmCRBN-V380E/I391V expressing CRBN-/- HEK293T cells treated with DMSO or thalidomide (Tha) for 16 h. Data information: Data in bar graphs were calculated as relative Mm) PLZF or SALL4 band intensity with Mm) PLZF or SALL4, respectively. The band intensity of DMSO was 100. Error bars denote ± standard deviation (n = 3). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001).(C, D) Immunoblot analysis of AGIA-GgPLZF (C) or -GgSALL4 (D) in FLAG-GgCRBN-WT or -GgCRBN-I390V expressing CRBN-/- HEK293T cells treated with DMSO or thalidomide (Tha) for 16 h. (E, F) Immunoblot analysis of AGIA-Mm PLZF (E) or -MmSALL4 (F) in FLAG-MmCRBN-WT expressing CRBN-/- HEK293T cells treated with indicated concentration of DMSO, thalidomide (Tha) or 5-hydroxythalidomide (5-HT) for 16 h. (G, H) Immunoblot analysis of AGIA-GgPLZF (G) or -GgSALL4 (H) in FLAG-GgCRBN-WT expressing CRBN-/- HEK293T cells treated with indicated concentration of DMSO, thalidomide (Tha) or 5-hydroxythalidomide (5-HT) for 16 h. Data information: Data in bar graphs were calculated as relative (Gg or Mm) PLZF or SALL4 band intensity with (Gg or Mm) PLZF or SALL4, respectively. The band intensity of DMSO was 100. Error bars denote ± standard deviation (n = 3). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001).
Figure 6(E) Skeletal patterning of forelimb and hindlimb in E6 chicken embryos infected RCAN virus packaging shControl (n = 10) or shPLZF (n = 9) were analysed by Victoria blue staining. h; humerus, r; radius, u; ulna, fe; femur, fi; fibula, t; tibia. Data in bar graphs were calculated as relative bone length with bone length of shControl as 100. Error bars denote ± standard deviation (forelimb and hindlimb shControl (n = 3), forelimb shPLZF (n = 5) or hindlimb shPLZF (n = 6)). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (**P < 0.01 and ****P < 0.0001).
Figure 7(F) Photographs show E4 chicken embryos treated with HBC (control, n = 10, C1, C3, C5 and C6), 1 µg/µl (3.9 mM) thalidomide (n = 11, T1, T2, T6 and T9) or 1 µg/µl (3.6 mM) 5-hydroxythalidomide (5-HT) (n = 10, H1, H2, H8 and H9) corresponding to immunoblot analysis in Fig 7G. Red arrows show treated regions. Data information Data in bar graphs in (F) were calculated as relative limb bud length with left limb bud length as one. Error bars denote ± standard deviation (F), n = 4). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, **P < 0.01, and ***P < 0.001).(G) Immunoblot analysis of PLZF or SALL4 from tissue of E4 chicken forelimb bud. Endogenous PLZF or SALL4 protein expression was detected by immunoblot using chicken embryos treated with HBC (n =10), 1 µg/µl (3.9 mM) thalidomide (n =11) and 1 µg/µl (3.6 mM) 5-hydroxythalidomide (5-HT) (n =10). Data information: Data in bar graphs in (G) were calculated as relative GgPLZF or GgSALL4 band intensity with GgPLZF or GgSALL4, respectively. The band intensity of HBC was 100. Error bars denote ± standard deviation ((E), n = 3 or (G), n = 4). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (NS = Not Significant, *P < 0.05, and **P < 0.01).
Figure 8E3 chicken forelimb bud electroporated with empty pCAGGS or pCAGGS-GgPLZF was treated with HBC or 1 µg/µl (3.9 mM) thalidomide, and the chicken forelimb bud phenotypes were observed in bright field (BF) or fluorescence of enhanced GFP (EGFP). Fgf10/Fgf8 expression was analysed by section by in situ hybridisation. (Q) Quantitative analysis of phenotypes of chicken limb bud Data in bar graphs in (Q) were calculated as relative limb bud length with limb bud length of control embryo as 100. Error bars denote ± standard deviation (n = 3). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (**P < 0.01).(R) Quantitative analysis of Fgf10/Fgf8 expression Data in bar graphs in (R) were calculated as relative Fgf10/Fgf8 mRNA expressions with Fgf10/Fgf8 mRNA expressions of control embryo as 100. Error bars denote ± standard deviation (n = 3). P-values were calculated by one-way ANOVA with Tukey's post-hoc test (*P < 0.05 and **P < 0.01).
Figure 2(a, b). Decreased Drp1S616 phosphorylation in PINK1-null primary neurons. Lysates of DIV6 neurons generated from PINK1KO and PINK1WT control mice were immune-detected with antibodies against either phospho (Ser616)-Drp1 (pS616) or Drp1. Tuj1 was detected as a loading control (a). Quantitation of pS616/Drp1 for each experiment is shown (b). Student's test. **p<0.01. Data was presented as mean ± SEM of three independent experiments.(c, d). TEM analysis of mitochondria in mouse substantial nigra. Sections of substantial nigra from 18-month-old PINK1KO and PINK1WT control mice were analyzed under TEM (c, left panels, scale bar=2 µm). Amplified images of mitochondria from dot line frames gated region are shown (c, right panels, scale bar=1 µm). Asterisks indicate abnormal elongated mitochondria. A quantitative analysis of mitochondria with perimeter >2µm is shown (d). Over 300 identifiable mitochondria (cristae and/or double membrane) per mice were randomly and blindly selected for perimeter measurement. Student's test. *p<0.05. Data was presented as mean ± SEM of three independent experiments.(e, f). Overexpression of Drp1WT, but not phosphorylation mutant Drp1S616A, suppresses mitochondrial fusion in dendrites of PINK1KO primary neurons. (e). DIV4 neurons generated from PINK1 WT control or PINK1KO mice were co-transfected mito-GFP with either pcDNA3.1 (3.1) or Drp1WT or Drp1S616A. Neurons were immuno-detected with an antibody against myc-tag (Myc) (red, to detect exogenous Drp1). Mito-GFP: green, to detect mitochondria. Representative images were shown, scale bar=25 μm. Magnified dendritic shafts gated by a dot line frame are shown (Zoom, right panels, scale bar=5 μm).(f). Quantification of Mito-index in dendrites of neurons from experiments of (e). 10 neurons were randomly selected and analyzed in each experiment. One-way ANOVA followed with Tukey's test. *p<0.05, ***p<0.001, ns: no significance. Data was presented as mean ± SEM of three independent experiments.
Figure 4(a-e). ATG5 is dispensable for PINK1-mediated mitochondrial fission. MEF cells derived from ATG5-null (ATG5KO) and their wildtype control mice (ATG5WT) were co-transfected Δ110-PINK1-GFP-FKBP /FRB-MTS with either PINK1WT or a kinase-dead mutant PINK1D384N. Cells were induced with 250 nM rapalog (Rapalog) for 2 hours to activate PINK1 kinase, followed by immunodetecting mitochondria (TOM20, red) and Δ110-PINK1-GFP-FKBP /FRB-MTS (PINK1-FKBP, green). Cells treated with solvent (DMSO) were used as a treatment control. Representative images are shown, scale bar=25 μm (a). Higher magnification images are also included (panels beneath the large cell images, scale bar=10 μm). Mitochondrial morphology in different transfections was quantified (b-e). Student's test. *p<0.05, ***p<0.001, ns: no significance. Data was presented as mean ± SEM of three independent experiments. For each condition, >100 cells were analyzed.(f-i). Drp1 is dispensable for recruitment of parkin to mitochondria. Drp1WT and Drp1KO HeLa cells were transfected with GFP-parkin. Cells were treated with 20 μM CCCP for 2 hours (f) followed by immunodetecting mitochondria (TOM20, red) and parkin (GFP-parkin, green). Nuclei were labeled with DAPI (blue). Cells treated with solvent (DMSO) were used as a treatment control. Representative images were shown, scale bar=25 μm. Cell with mitochondrial parkin (%) in different experimental group was quantified (g, i). One-way ANOVA followed with Tukey's test. ***p<0.001, ns: no significance. Data was presented as mean ± SEM of three independent experiments. For each condition, >100 cells were analyzed.(f-i). Drp1 is dispensable for recruitment of parkin to mitochondria. Drp1WT and Drp1KO HeLa cells were transfected with GFP-parkin. Cells were treated with 150 μM actinonin for 6 hours (h), followed by immunodetecting mitochondria (TOM20, red) and parkin (GFP-parkin, green). Nuclei were labeled with DAPI (blue). Cells treated with solvent (DMSO) were used as a treatment control. Representative images were shown, scale bar=25 μm. Cell with mitochondrial parkin (%) in different experimental group was quantified (g, i). One-way ANOVA followed with Tukey's test. ***p<0.001, ns: no significance. Data was presented as mean ± SEM of three independent experiments. For each condition, >100 cells were analyzed.
Figure 6(a, b). Drp1 rescues PINK1KO induced collapsed thorax phenotype in dATG7-deficient Drosophila. Representative SEM images of thorax from PINK1WT, (upper panels) and PINK1KO (lower panels) flies combined with dATP7 deficiency (dATG7KO) or their controls (dATG7WT), followed by expressing either mhc-gal4 (Mhc/+) alone, or Mhc-gal4 driven expression of human Drp1 (Mhc>hDrp1WT), mhc-gal4 driven expression of human Drp1S616A (Mhc>hDrp1S616A) or mhc-gal4 driven expression of human Drp1S616D (Mhc>hDrp1S616D) are shown (a). Quantitation of crushed thorax in flies with indicated genotypes is presented (b). > 100 flies for each experiment are analyzed. Scale bar=200µm. One-way ANOVA followed with Tukey's test. ***p<0.001. Data was presented as mean ± SEM of three independent experiments. (c, d) Drp1 rescues PINK1KO induced mitochondrial defects in dATG7KO Drosophila. Representative TEM images of indirect flight muscle from PINK1WT (upper panels) and PINK1KO (lower panels) flies combined with dATG7KO or their controls (dATG7WT), followed by expressing either mhc-gal4 (Mhc/+) alone, or Mhc-gal4 driven expression of human Drp1(Mhc>hDrp1WT), mhc-gal4 driven expression of human Drp1S616A(Mhc>hDrp1S616A) or mhc-gal4 driven expression of human Drp1S616D(Mhc>hDrp1S616D) are shown (c). Scale bar=1 μm. Normal and abnormal mitochondrial number per mm2 from each experimental group were and quantified. >3 flies/group and 3-6 pictures of different microscopic fields from each fly were analyzed per repeat. Bars represent the mean ± SEM of three independent experiments (d).(e). Drp1 rescues PINK1KO induced ATP reduction in Drosophila. ATP contents of thorax muscle tissues from the indicated genotypes were measured and normalized against the protein levels One-way ANOVA followed with Tukey's test. ***p<0.001, ns: no significance. Data was presented as mean ± SEM of three independent experiments.
Figure 7(g, h). Reduction of Drp1S616 phosphorylation in dermal fibroblasts derived from PD patients with PINK1 mutations. Lysates of cells derived from three unrelated normal control individuals (Control, C) and four familial PD patients with PINK1 mutations (PINK1mutant) were immunoblotted with antibodies recognizing either phosphor-Drp1S616 (pS616) or Drp1 protein. β-actin was detected as a loading control (g). A quantitative analysis of phospho-Drp1S616 (pS616/Drp1) is shown (h). Student's test. **p<0.01, ns: no significance. Data was presented as mean ± SEM of three independent experiments.Drp1S616 phosphorylation in dermal fibroblasts derived from sporadic PD patients. Lysates of cells derived from five unrelated normal control individuals (Control, C) and seven unrelated sporadic PD patients (Sporadic PD, SP) were immunoblotted with antibodies recognizing either phospho-Drp1S616 (pS616) or Drp1 protein. β-actin was detected as a loading control (i).Drp1S616 phosphorylation in dermal fibroblasts derived from sporadic PD patients. A quantitative analysis of phospho-Drp1S616 (pS616/Drp1) in sporadic cases is shown (j). Student's test. *p<0.05. Data was presented as mean ± SEM of three independent experiments.
Figure 2B. Volcano plot showing the magnitude of differential gene up-regulation, as measured by the respective log 2-fold change in expression level, at day (positive) or night (negative) versus statistical significance, as quantified by a negative binomial exact test (red dots, dashed line p<0.05 adjusted for multiple testing using Benjamin-Hochberg). Known circadian-regulated genes are labelled.
Figure 3B. Violin plots depicting the distribution of the quantified average cellular expression of Prok2, ProkR2, Vip, Vipr2, Avp, Avpr1a, Grp and Grpr in each sub-cluster in CT7.5 (left) and CT15.5 (right) datasets.