Circadian and Feeding Rhythms Orchestrate the Diurnal Liver ...

31 downloads 0 Views 6MB Size Report
Martin, Laetitia Da Silva, Ivan Montoliu, Sebastiano Collino, Francois-Pierre. Martin, Joanna Ratajczak, Carles Cantó, Martin Kussmann, Felix Naef, and Frédéric.
Cell Reports, Volume 20

Supplemental Information

Circadian and Feeding Rhythms Orchestrate the Diurnal Liver Acetylome Daniel Mauvoisin, Florian Atger, Loïc Dayon, Antonio Núñez Galindo, Jingkui Wang, Eva Martin, Laetitia Da Silva, Ivan Montoliu, Sebastiano Collino, Francois-Pierre Martin, Joanna Ratajczak, Carles Cantó, Martin Kussmann, Felix Naef, and Frédéric Gachon

Supplemental Information Supplemental Figures Figure S1. Experimental design and data quality check, related to Figure 1

A

B 1062

C

236

297

R=0.79

R=0.81

R=0.83

R=0.78

R=0.70

R=0.78

R=0.73

R=0.62

D 0.2

Nucleus Nucleus & Cytoplasm Mitochondria Nucleus & Cytoplasm

TE

NE

Non Nuclear localized

1/12

A. Workflow of the SILAC-based MS analysis of acetylated proteins from total mouse liver extracts (TE) and purified nuclei (NE). The lower graph represents the number of acetylation sites for each time point in TE (black bars, left Y axis) and NE (red bars, right Y axis). B. Number of acetylated sites quantified in TE (red circle, 1298 acetylated sites identified) and NE (green circle, 533 acetylated sites identified). C. Pearson correlation analysis of biological replicates at the eight time points. All values are log2 ratios to the common reference samples. The biological replicates are well correlated (76% average Pearson correlation) and showed similar spread, indicating that the quality of the protein purifications was fairly homogenous. D. Fractions of raw intensity signals quantified with SILAC-MS for nuclear and non-nuclear located proteins in TE (left graph) and NE (right graph).

2/12

Figure S2. Overlap of rhythmic acetylation sites and expression of Acetyl-CoA synthesizing enzymes, related to Figure 2

A

B ZT 0

6

12

18

6

12

18

98

n=34

24

3

HMGCS2 K83; HMGCS2 K310; STIP1 K325 CPS1 K55; CPS1 K527; CPS1 K915; SOD2 K122; SOD2 K68

5 14

C ZT 0

n=17

0

3

20

2

15

1

-1 10

0

3

6

9

-1

12 15 18 21 24

ZT

TE

TE

B M A LK O

A C SS 3

A cs s3

L

K

O

C R Y

TE

TE A M B s2 A

A

cs

C

LY

B

M

A cs s2

A

LK

O

C R Y

TE

TE W T

C R Y

A C LY

A

TE

TE

B M A LK O

A C SS 3

L

K

O

C R Y

TE

TE M B s2 A

A

cs

C

LY

B

M

A cs s2

A

LK

O

W T

C R Y

C R Y

TE

TE

KO/WT log2

L

L

L

A C SS 2

LY

C

A C SS 3 L

A C LY

A cs s3

ACSS3

ACLY

40

TE

TE C R Y

A C SS 3

B M A LK O

TE O K L

A cs s3

A M B s2 cs

O LK

C R Y

TE A cs s2

50 0

A

10

TE

20

A

-3

ACSS2

30

M

0

ACLY

B

-2

-2

0 40

TE

12 15 18 21 24

-1

LY

9

0

C

6

Bmal1 KO

1

W T

3

Cry1/2 DKO

2 -2

A

0

5 50

3 -1

C R Y

-1

ACSS2 p=0.06 Phase=1.32 ACLY p=0.005 Phase=23.7

0

A C LY

SILAC (L/H) ratio log2

0 -2

G

H

L

E

10 L

ZT

ACLY ACSS2 ACSS3

A C SS 3

12 15 18 21 24

L

9

A C SS 3

6

0 15

A C SS 2

3

5 20

A C SS 2

0

-1

10

L

ACLY p=0.0009 Phase=4.9

1

L

0 -2

15

A

-1

-2 2

0

LY

12 15 18 21 24

C

9

-1 3

20

A

6

5

LY

3

0

C

0

Average Absolute MS intensities (relative to ACSS2)

SILAC (L/H) ratio log2

-2 0

F

A

ACSS3 p=0.23 Phase=7.8 ACSS2 p=0.009 Phase=6.1

Average Absolute MS intensities (relative to ACSS2)

D

ACSS2

30

-3 -1

0

3

6

9

12 15 18 21 24

0 40

30

L A C SS 2

LY C

A

10 50

L

20

-2 0

3/12

A. Venn diagram showing the overlap of rhythmic acetylation sites between the TE dataset (106 rhythmic acetylation sites), the NE dataset (27 rhythmic acetylation sites), and the dataset of (Masri et al., 2013) (19 rhythmic acetylation sites). B-C. Heat maps showing 12 hours rhythmic acetylation sites normalized by their corresponding total protein amount in (B) TE (n=34) and (C) NE (n=17) under light-dark and night-restricted feeding conditions. Data were standardized by rows and gray blocks indicate missing protein data. The polar plots on the right of each heatmap display peak phase distribution of the 12 hours rhythmic acetylation sites in each extract. Colors indicate acetylation sites with a corresponding total protein having a defined mitochondrial (n=17 sites; red) or cytoplasmic (n=20 sites; blue) localization in TE, or acetylation sites from NE (n=17; green). D-E. Temporal protein expression of ACLY (black), ACSS2 (blue), and ACSS3 (orange) in TE (D) and NE (E) (ACSS3 not detected in NE). The values represent the mean ± SEM from two independent biological samples. F-G. Raw intensity signals quantified by SILAC-MS for ACLY (black), ACSS2 (blue), and ACSS3 (orange) in TE (F) and NE (G) (ACSS3 not detected in NE). Data are expressed relative to ACSS2 average signal, error bars represent min to max. H. Differential protein expression level of ACLY, ACSS2, and ACSS3 in TE in Cry1/2 DKO mice (red bars) and Bmal1 KO mice (green bars) compared to respective WT animals. Data show average level ± SEM in four samples collected around the clock. Data of panels D, F, H, and E, G are from (Mauvoisin et al., 2014) and (Wang et al., 2017), respectively.

4/12

Figure S3. Temporal nuclear expression of SIRT1 and SIRT7, related to Figure 3

B

24 hour-Rhythms 12 hour-Rhythms N=106 N=34

SILAC (L/H) ratio log2

A

0.2 0.0 -0.2 -0.4 -0.6

0.2 0.0 -0.2 -0.4 -0.6 -0.8 -1.0

0

3

6

9

12 15 18 21 24

ZT 0.4

SIRT2 p=0.003 0

3

6

9

Phase=0.75 12 15 18 21 24

ZT

SILAC (L/H) ratio log2

SILAC (L/H) ratio log2

D 0.4

SIRT1

-0.8 -1.0

C

RPS19BP1 p=0.0003 Phase=23.3

0.2 0.0 -0.2 -0.4 -0.6 -0.8 -1.0

SIRT7 p=0.0001 Phase=4.22 0

3

6

9

12 15 18 21 24

ZT

A. Proportion of SIRT3 targets (orange) versus non-SIRT3 targets (blue) in rhythmic acetylation sites with 24 hourrhythms and 12 hour-rhythms. B. Temporal nuclear expression of SIRT1 (blue) and its coactivator RPS19BP1 (black). C. Temporal nuclear expression of SIRT2. D. Temporal nuclear expression of SIRT7. For B, C, and D, the values represent the mean ± SEM from two independent biological samples. Data are extracted from (Wang et al., 2017).

5/12

Figure

S4.

Additional

metabolic

pathways

affected

by

rhythmic

acetylation,

related to Figure 5

23 0 1 22 2 21 3 20 19 18 17 16 15 9 14 13 10 12 11

ZT

4 5 6 7 8

°* ATP5B °* ATP5C1 * NDUFA2 NDUFA4 NDUFA5

Oxida/ve Phosphoryla/on

ATP5A1

Qu a Ide n3fie n3 d fie d

Rhythmic

Ac-sites

°* NDUFA10 NDUFAB1 ATP5J °* ATP5H °* NDUFB10 NDUFS1 CYB5;CYB5A UQCRC1 ATP5O * ATP5L °* NDUFV1 NDUFV2 UQCRB °* UQCRC2 MT-ATP8 UQCRCQ

ATP5D ATP5F1

COX6B1

ATP5J2

Ac-level different in Cry 1/2 DKO° and Bmal1 KO*mice Lysine AASS

I

V

III

°*

IV

II

Succinyl-CoA Synthase Succinate SUCLG2 SUCLA2 SUCLG1 dehydrogenase

°*

α-aminoadipate-6-semialdehyde

°*

SDHA

a-ketoglutarate dehydrogenase DLST DLD Succinyl-CoA

°*

Succinate

AASDH

Lysine degrada/on

COX4I1 COX5A

Fumarate

Fumarate hydratase

°*

α-ketoglutarate

FH

Isocitrate dehydrogenase

α-aminoadipate

°* IDH3G

TCA cycle

IDH2 AADAT

Malate dehydrogenase

Malate

D-isocitrate

*

MDH2

Aconitase ACO2

°*

α-ketoadipate BCKDK BCKDHA

Oxaloacetate

Citrate

Pyruvate carboxylase PC

°*

Citrate Synthase Acetyl-CoA CS

°*

Glutaryl-CoA

Pyruvate Pyruvate dehydrogenase complex

GCDH

*

PDHA1 PDHB DLAT

°*

DLD

Crotonyl-CoA ECHS1

β-Hydroxybutyryl-CoA

*

HADH

acetoacetyl-CoA

to ketogenesis

Rhythmic protein acetylation and their regulation by the circadian clock for lysine degradation, oxidative phosphorylation, and TCA cycle pathways. Each dot represents a unique acetylation site within the protein of interest. Grey and black dots represent non-rhythmic identified and quantified acetylation sites, respectively, whereas phase of rhythmic acetylation sites are color-coded. Superscripted red dots and green asterisks indicate protein acetylation level significantly different in Cry1/2 DKO and Bmal1 KO mice, respectively.

6/12

Figure S5. Differential acetylation of the enzymes involved in the different metabolic pathways in Cry1/2 and Bmal1 KO mice, related to Figure 5 and S4

A

Fa3y acid oxida6on 2.0

§ §

1.0

§§ * § *

0.5 0.0

§

§

§

* ** *

§

ACSF2

-2.0

1.0

ECI1

ACADL

ACAD12

EHHADH

HADH

* *

*

§

*

*

*

**

* *

*

*

**

*

*



*

**

*

§

*

*

§ HADHA

§

§ *

0.5

* *

*

*

0.0

-0.5 -1.0

ARG1 *

ASS1

§

*

*

*

Ketogenesis 1.0

*

§

§*

ACAA2

C §

*

-1.0

2.0

ACAD10

*

-0.5

D

*

*

Urea cycle

0.5 § 0.0

*

§

-1.0

B

*

*

§

-0.5 -1.5

KO/WT (log2)

§

KO/WT (log2)

KO/WT (log2)

1.5

OTC

GLUD1

CPS1

ACAT1

GOT2

HMGCS2

Ethanol metabolism

KO/WT (log2)

1.5 1.0

-0.5

*

-1.0

E 1.5

ALDH1A1 §

ADH1

-2.0

KO/WT (log2)

* *

*

§

*

* §

§

*

0.0

-1.5

§

§ § § ALDH1A7

§

ALDH1L1

§

§

§

ALDH2

ALDH4A1

§

ALDH6A1

§

ALDH9A1

Oxida6ve phosphoryla6on §

*

1.0

*

0.5

§

§

§

§§ *

* §

*

§

*

§ §

§

§

*

§

§

0.0

§

-0.5 -1.0

F

*

ATP5A1

ATP5B

ATP5C1

§

* ATP5F1

ATP5H

ATP5J

* ATP5L

ATP5O

UQCRB §

TCA cycle 1.0

KO/WT (log2)

§

§

0.5

0.5

*

*

*

§

§ § *

§

* § *

§ §

*

* * *

* *

*

§

*

§

*

*

* *

* § §

0.0

*

-0.5

§

§

**

*

*

§

§

*

-1.0

G 1.5 KO/WT (log2)

1.0 0.5

DLD

PDHA1

CS

IDH2

SUCLG2

SDHA

FH

MDH2

PC

Lysine degrada6on § § §

*

*

§

0.0 -0.5 -1.0 -1.5

ACO2

AASS

HDAH

*

ZT 0 6 12 18 0 6 12 18 Cry1/2 DKO Bmal1 KO paired t-test KO-WT § p