# MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in Klf10 deficient mice

Pierre Leclère, Déborah Rousseau, Stéphanie Patouraux, Sophie Guérin, Stéphanie Bonnafous, Aline Gréchez-Cassiau, Anthony Ruberto, Carmelo Luci, Malayannan Subramaniam, Albert Tran, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Pierre Leclère, Déborah Rousseau, Stéphanie Patouraux, Sophie Guérin, Stéphanie Bonnafous, et al. MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in Klf10 deficient mice. Scientific Reports, 2020, 10 (1), pp.12139. 10.1038/s41598-020-69085-w. hal-04496680

HAL Id: hal-04496680

https://hal.science/hal-04496680

Submitted on 27 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in Klf10 deficient mice

Pierre Leclère, Déborah Rousseau, Stéphanie Patouraux, Sophie Guérin, Stéphanie Bonnafous, Aline Gréchez-Cassiau, Anthony Ruberto, Carmelo Luci, Malayannan Subramaniam, Albert Tran, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Pierre Leclère, Déborah Rousseau, Stéphanie Patouraux, Sophie Guérin, Stéphanie Bonnafous, et al. MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in Klf10 deficient mice. Scientific Reports, 2020, 10 (1), pp.12139. 10.1038/s41598-020-69085-w. hal-04496680

HAL Id: hal-04496680

https://hal.science/hal-04496680v1

Submitted on 27 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

OPEN

# MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in *Klf10* deficient mice

Pierre S. Leclère<fnref n="1" />,², Déborah Rousseau², Stéphanie Patouraux³, Sophie Guérin<fnref n="1" />,
Stéphanie Bonnafous³, Aline Gréchez-Cassiau<fnref n="1" />, Anthony A. Ruberto<fnref n="1" />, Carmelo Luci²,
Malayannan Subramaniam⁴, Albert Tran³, Franck Delaunay<fnref n="1" />, Philippe Gual²,⁵✉ &
Michèle Teboul<fnref n="1" />,⁵✉

A large number of hepatic functions are regulated by the circadian clock and recent evidence suggests that clock disruption could be a risk factor for liver complications. The circadian transcription factor Krüppel like factor 10 (KLF10) has been involved in liver metabolism as well as cellular inflammatory and death pathways. Here, we show that hepatic steatosis and inflammation display diurnal rhythmicity in mice developing steatohepatitis upon feeding with a methionine and choline deficient diet (MCDD). Core clock gene mRNA oscillations remained mostly unaffected but rhythmic *Klf10* expression was abolished in this model. We further show that *Klf10* deficient mice display enhanced liver injury and fibrosis priming upon MCDD challenge. Silencing *Klf10* also sensitized primary hepatocytes to apoptosis along with increased caspase 3 activation in response to TNFα. This data suggests that MCDD induced steatohepatitis barely affects the core clock mechanism but leads to a reprogramming of circadian gene expression in the liver in analogy to what is observed in other experimental disease paradigms. We further identify KLF10 as a component of this transcriptional reprogramming and a novel hepato-protective factor.

Non-alcoholic fatty liver diseases (NAFLDs) are the most frequent chronic liver pathologies worldwide with a global prevalence ranging from 22 to 28%<fnref n="1" />. NAFLDs cover a spectrum of hepatic abnormalities ranging from reversible steatosis to non-alcoholic steatohepatitis (NASH) and steatofibrosis which may lead to cirrhosis and ultimately hepatocellular carcinoma². NASH is characterized by (i) fatty hepatocytes, (ii) inflammatory foci, (iii) ballooned hepatocytes and (iv) increased hepatocytes apoptosis and necrosis. Liver injury associated with NASH leads to elevated systemic levels of alanine transaminase (ALT) and aspartate transaminase (AST) activity³. Sustained activation of cellular stress in lipid-overloaded hepatocytes and increased activation of death receptor signaling including the Tumor Necrosis Factor alpha (TNFα) pathway, have been associated with hepatocellular injury during NASH⁴. Additional modes of cell death including necroptosis and pyroptosis could also participate to liver injury during NASH⁵,⁶. Hepatocyte injury derived signals such as Danger Associated Molecular Pattern Signals (DAMPs) as well as pro-inflammatory cytokines and chemokines are also strong drivers of NASH progression⁷. Although an increasing number of clinical trials are now in progress, approved pharmacological therapy for burned-out NASH is still missing.

The mammalian circadian timing system (CTS) coordinates most of physiology and behavior with the external light/dark cycle. This system is organized hierarchically with a central pacemaker in the suprachiasmatic

<footnote n="1">¹Université Côte D'Azur, CNRS, INSERM, iBV, Nice, France. ²Université Côte D'Azur, INSERM, C3M, Nice, France. ³Université Côte D'Azur, CHU, INSERM, C3M, Nice, France. ⁴Department of Biochemistry and Molecular Biology, Mayo Clinic, Rochester, MN, USA. ⁵These authors jointly supervised: Philippe Gual and Michèle Teboul. ✉email: Philippe.gual@inserm.fr; Michele.TEBOUL@univ-cotedazur.fr</footnote>
<figure><img src="image_3.png" /></figure>

OPEN

# MCD diet-induced steatohepatitis generates a diurnal rhythm of associated biomarkers and worsens liver injury in *Klf10* deficient mice

Pierre S. Leclère<fnref n="1" />,², Déborah Rousseau², Stéphanie Patouraux³, Sophie Guérin<fnref n="1" />,
Stéphanie Bonnafous³, Aline Gréchez-Cassiau<fnref n="1" />, Anthony A. Ruberto<fnref n="1" />, Carmelo Luci²,
Malayannan Subramaniam⁴, Albert Tran³, Franck Delaunay<fnref n="1" />, Philippe Gual²,⁵✉ &
Michèle Teboul<fnref n="1" />,⁵✉

A large number of hepatic functions are regulated by the circadian clock and recent evidence suggests that clock disruption could be a risk factor for liver complications. The circadian transcription factor Krüppel like factor 10 (KLF10) has been involved in liver metabolism as well as cellular inflammatory and death pathways. Here, we show that hepatic steatosis and inflammation display diurnal rhythmicity in mice developing steatohepatitis upon feeding with a methionine and choline deficient diet (MCDD). Core clock gene mRNA oscillations remained mostly unaffected but rhythmic *Klf10* expression was abolished in this model. We further show that *Klf10* deficient mice display enhanced liver injury and fibrosis priming upon MCDD challenge. Silencing *Klf10* also sensitized primary hepatocytes to apoptosis along with increased caspase 3 activation in response to TNFα. This data suggests that MCDD induced steatohepatitis barely affects the core clock mechanism but leads to a reprogramming of circadian gene expression in the liver in analogy to what is observed in other experimental disease paradigms. We further identify KLF10 as a component of this transcriptional reprogramming and a novel hepato-protective factor.

Non-alcoholic fatty liver diseases (NAFLDs) are the most frequent chronic liver pathologies worldwide with a global prevalence ranging from 22 to 28%<fnref n="1" />. NAFLDs cover a spectrum of hepatic abnormalities ranging from reversible steatosis to non-alcoholic steatohepatitis (NASH) and steatofibrosis which may lead to cirrhosis and ultimately hepatocellular carcinoma². NASH is characterized by (i) fatty hepatocytes, (ii) inflammatory foci, (iii) ballooned hepatocytes and (iv) increased hepatocytes apoptosis and necrosis. Liver injury associated with NASH leads to elevated systemic levels of alanine transaminase (ALT) and aspartate transaminase (AST) activity³. Sustained activation of cellular stress in lipid-overloaded hepatocytes and increased activation of death receptor signaling including the Tumor Necrosis Factor alpha (TNFα) pathway, have been associated with hepatocellular injury during NASH⁴. Additional modes of cell death including necroptosis and pyroptosis could also participate to liver injury during NASH⁵,⁶. Hepatocyte injury derived signals such as Danger Associated Molecular Pattern Signals (DAMPs) as well as pro-inflammatory cytokines and chemokines are also strong drivers of NASH progression⁷. Although an increasing number of clinical trials are now in progress, approved pharmacological therapy for burned-out NASH is still missing.

The mammalian circadian timing system (CTS) coordinates most of physiology and behavior with the external light/dark cycle. This system is organized hierarchically with a central pacemaker in the suprachiasmatic

<footnote n="1">¹Université Côte D'Azur, CNRS, INSERM, iBV, Nice, France. ²Université Côte D'Azur, INSERM, C3M, Nice, France. ³Université Côte D'Azur, CHU, INSERM, C3M, Nice, France. ⁴Department of Biochemistry and Molecular Biology, Mayo Clinic, Rochester, MN, USA. ⁵These authors jointly supervised: Philippe Gual and Michèle Teboul. ✉email: Philippe.gual@inserm.fr; Michele.TEBOUL@univ-cotedazur.fr</footnote>
nuclei (SCN) of the hypothalamus that receives light inputs from the retina and in turn entrains cellular clocks present in virtually all peripheral tissues and extra-hypothalamic brain regions via internal synchronizers⁸. Cellular circadian clocks control a vast array of biological processes in a tissue specific manner and the liver is the organ showing the most extensive circadian regulation⁹. Many features of modern society life style including sleep deprivation, shiftwork, light at night, overfeeding, and night eating contribute to the misalignment of internal body time with the environmental cycle and this may in turn compromise health¹⁰. Consistently, epidemiological studies increasingly support the notion that a poorly synchronized or disrupted CTS is an important risk factor for the development of many chronic metabolic diseases including dyslipidemia, overweight and insulin resistance, and cancer¹¹⁻¹³. Reciprocally, circadian physiology can be altered by diseases¹⁴⁻¹⁶. This is exemplified by the reprogramming of the liver and serum metabolome as well as hepatic transcriptome occurring in obese mice¹⁷⁻²⁰. The interaction of the circadian clock with several mechanisms involved in the onset of hepatic-steatosis and its progression to NASH has been well documented in several reviews²¹⁻²³. Global and tissue specific knock-out or mutation of clock genes (*Rev-erbα/β⁻*, *Cry1/2⁻*, *Bmal1⁻/−* and *ClockΔ19*, *Bmal1Δhbp*) promote hepatic steatosis, and alter blood glucose and lipid homeostasis in lean and high fat diet (HFD) fed obese mice²⁴⁻²⁸. Accordingly, sustained chrono-disruption in mice via experimental chronic jet lag leads to spontaneous development of NASH and hepatocellular carcinoma²³,²⁹. Several transcription factors belonging to the *Krüppel Like Factor* (KLF) family have recently emerged as important circadian regulators involved in liver functions and diseases³⁰⁻³². While some of these KLF family members including KLF4 and KLF6 are potential players in NAFLD development³³⁻³⁵, the significance of KLF10 in liver physiopathology remains unknown although it has been involved in several mechanisms which are dysregulated in NAFLD³⁶⁻³⁸. We previously reported that hepatic *Klf10* is a circadian output of the clock as a result of a direct regulation by the BMAL1/CLOCK heterodimer. Transcriptome profiling of *Klf10⁻/−* mouse liver further identified a differentially expressed gene set that was significantly enriched for lipid and carbohydrate metabolism processes. In addition, KLF10 was shown to regulate hepatic gluconeogenesis through the direct transcriptional repression of *Pepck*³⁰. It has also been recently reported that *Klf10* hepatic expression is strongly increased in fatty liver of obese mice³⁹. KLF10 has also been shown to be involved in inflammation and cell death processes³⁷,⁴⁰. All these features suggest that KLF10 may have a role in the development of steatohepatitis that remains to be explored. In the current study, we use mice fed a methionine and choline deficient diet (MCDD), which is a rodent model of oxidative-stress mediated steatohepatitis. On the one hand, a lack of choline in this diet, hampers the export of triglycerides (TG) via a very low-density lipoprotein (VLDL) packaging from hepatocytes, resulting in hepatic steatosis. On the other hand, the essential amino acid methionine is required for the synthesis of S-adenosylmethionine (SAM) and glutathione, which are both antioxidants. Although this model was widely used and despite the increasing numbers of reports linking circadian disorders and NAFLD, no study investigated the influence of the MCDD challenge on the circadian rhythmicity. In the present study, we examine the steatohepatitis features around the 24-h cycle and the role of the circadian transcription factor KLF10 in the progression of the disease, in mice fed a MCDD.

## Results

**Diurnal rhythmicity of liver steatosis, inflammation and fibrosis priming in MCDD fed mice.** Several experimental paradigms using nutritional challenges in mice have shown a reprogramming of the circadian transcriptome and physiology in many organs suggesting that these pathogenic conditions have the potential to unveil disease related rhythmic patterns¹⁵,¹⁷⁻²⁰. To evaluate whether steatohepatitis features, including hepatic steatosis, inflammation and liver injury, display a circadian rhythmicity, WT mice were fed with a control diet (CD) or a methionine and choline deficient diet (MCDD) for 4 weeks. Blood and liver were collected every 6 h over 24 h at ZT3, 9, 15 and 21, ZT0 and ZT12 being the time when the light is switched on and off, respectively (Fig. 1A). Hepatic steatosis evaluated by the percentage of lipid droplets displayed a robust diurnal rhythmicity with a calculated acrophase during the rest phase at ZT8 (Fig. 1B, Supplementary Table S1) upon MCDD challenge. Consistent with this observation, hepatic expression of *Fsp27* (also known as *Cidec*) and *Pnpla2* which encode for two proteins involved in lipid droplet fusion and unpacking respectively, also showed a rhythmic profile with a peak of expression during the rest phase (Fig. 1C, Supplementary Table S1). The hepatic triglyceride (TG) content displayed rhythmicity in CD fed mice as expected from previous work⁴¹. Expectedly, the TG content was increased upon MCDD, but was no longer rhythmic (Fig. 1B, Supplementary Table S1). These observations suggest that although the elevated content of TG remains constant during the 24 h cycle, the dynamics of CIDE-regulated lipid droplet formation may be a rhythmic process. Hepatic inflammation evaluated by the number of inflammatory foci also exhibited a robust diurnal rhythmicity during steatohepatitis (Fig. 1D, Supplementary Table S1). Interestingly, expression of the inflammatory markers *Tnfα* and *Ccl2*, which are not rhythmic in the CD group, displayed a de novo rhythmicity upon MCDD-induced steatohepatitis (Fig. 1E, Supplementary Table S1). As expected, upon MCDD, both AST and ALT activities were higher in the MCDD fed mice than in the CD fed mice (Fig. 1F, Supplementary Table S1). We observed that both ALT and AST activity varied depending on the time of the day, ALT activity being at ZT3 than at ZT9, and AST activity displaying a 24 h oscillation. (Fig. 1F, Supplementary Table S1). In addition, the fibrosis priming markers *Tgfβ1*, *Col1α1* and *Timp1* were not rhythmic upon CD and, as expected, were upregulated upon MCDD. While *Tgfβ1* did not display circadian rhythmicity, *Col1α1* and *Timp1* both gained rhythmicity upon MCDD (Fig. 1G, Supplementary Table S1). Collectively these results demonstrate that multiple markers that have an arrhythmic expression in chow diet display an up-regulation associated with a time-of-day dependent variation in the MCDD induced NASH mouse model.

**Steatohepatitis induces discrete alterations of the liver molecular clock and abolishes *Klf10* oscillation.** To investigate whether the observed de novo rhythmicity of steatohepatitis features could be
nuclei (SCN) of the hypothalamus that receives light inputs from the retina and in turn entrains cellular clocks present in virtually all peripheral tissues and extra-hypothalamic brain regions via internal synchronizers⁸. Cellular circadian clocks control a vast array of biological processes in a tissue specific manner and the liver is the organ showing the most extensive circadian regulation⁹. Many features of modern society life style including sleep deprivation, shiftwork, light at night, overfeeding, and night eating contribute to the misalignment of internal body time with the environmental cycle and this may in turn compromise health¹⁰. Consistently, epidemiological studies increasingly support the notion that a poorly synchronized or disrupted CTS is an important risk factor for the development of many chronic metabolic diseases including dyslipidemia, overweight and insulin resistance, and cancer¹¹⁻¹³. Reciprocally, circadian physiology can be altered by diseases¹⁴⁻¹⁶. This is exemplified by the reprogramming of the liver and serum metabolome as well as hepatic transcriptome occurring in obese mice¹⁷⁻²⁰. The interaction of the circadian clock with several mechanisms involved in the onset of hepatic-steatosis and its progression to NASH has been well documented in several reviews²¹⁻²³. Global and tissue specific knock-out or mutation of clock genes (*Rev-erbα/β*⁻, *Cry1/2*⁻, *Bmal1*⁻/*and ClockΔ19*, *Bmal1*¹*ahcp*) promote hepatic steatosis, and alter blood glucose and lipid homeostasis in lean and high fat diet (HFD) fed obese mice²⁴⁻²⁸. Accordingly, sustained chrono-disruption in mice via experimental chronic jet lag leads to spontaneous development of NASH and hepatocellular carcinoma²³,²⁹. Several transcription factors belonging to the *Krüppel Like Factor* (KLF) family have recently emerged as important circadian regulators involved in liver functions and diseases³⁰⁻³². While some of these KLF family members including KLF4 and KLF6 are potential players in NAFLD development³³⁻³⁵, the significance of KLF10 in liver physiopathology remains unknown although it has been involved in several mechanisms which are dysregulated in NAFLD³⁶⁻³⁸. We previously reported that hepatic *Klf10* is a circadian output of the clock as a result of a direct regulation by the BMAL1/CLOCK heterodimer. Transcriptome profiling of *Klf10*⁻/⁻ mouse liver further identified a differentially expressed gene set that was significantly enriched for lipid and carbohydrate metabolism processes. In addition, KLF10 was shown to regulate hepatic gluconeogenesis through the direct transcriptional repression of *Pepck*³⁰. It has also been recently reported that *Klf10* hepatic expression is strongly increased in fatty liver of obese mice³⁹. KLF10 has also been shown to be involved in inflammation and cell death processes³⁷,⁴⁰. All these features suggest that KLF10 may have a role in the development of steatohepatitis that remains to be explored. In the current study, we use mice fed a methionine and choline deficient diet (MCDD), which is a rodent model of oxidative-stress mediated steatohepatitis. On the one hand, a lack of choline in this diet, hampers the export of triglycerides (TG) via a very low-density lipoprotein (VLDL) packaging from hepatocytes, resulting in hepatic steatosis. On the other hand, the essential amino acid methionine is required for the synthesis of S-adenosylmethionine (SAM) and glutathione, which are both antioxidants. Although this model was widely used and despite the increasing numbers of reports linking circadian disorders and NAFLD, no study investigated the influence of the MCDD challenge on the circadian rhythmicity. In the present study, we examine the steatohepatitis features around the 24-h cycle and the role of the circadian transcription factor KLF10 in the progression of the disease, in mice fed a MCDD.

## Results

**Diurnal rhythmicity of liver steatosis, inflammation and fibrosis priming in MCDD fed mice.** Several experimental paradigms using nutritional challenges in mice have shown a reprogramming of the circadian transcriptome and physiology in many organs suggesting that these pathogenic conditions have the potential to unveil disease related rhythmic patterns¹⁵,¹⁷⁻²⁰. To evaluate whether steatohepatitis features, including hepatic steatosis, inflammation and liver injury, display a circadian rhythmicity, WT mice were fed with a control diet (CD) or a methionine and choline deficient diet (MCDD) for 4 weeks. Blood and liver were collected every 6 h over 24 h at ZT3, 9, 15 and 21, ZT0 and ZT12 being the time when the light is switched on and off, respectively (Fig. 1A). Hepatic steatosis evaluated by the percentage of lipid droplets displayed a robust diurnal rhythmicity with a calculated acrophase during the rest phase at ZT8 (Fig. 1B, Supplementary Table S1) upon MCDD challenge. Consistent with this observation, hepatic expression of *Fsp27* (also known as *Cidec*) and *Pnpla2* which encode for two proteins involved in lipid droplet fusion and unpacking respectively, also showed a rhythmic profile with a peak of expression during the rest phase (Fig. 1C, Supplementary Table S1). The hepatic triglyceride (TG) content displayed rhythmicity in CD fed mice as expected from previous work⁴¹. Expectedly, the TG content was increased upon MCDD, but was no longer rhythmic (Fig. 1B, Supplementary Table S1). These observations suggest that although the elevated content of TG remains constant during the 24 h cycle, the dynamics of CIDE-regulated lipid droplet formation may be a rhythmic process. Hepatic inflammation evaluated by the number of inflammatory foci also exhibited a robust diurnal rhythmicity during steatohepatitis (Fig. 1D, Supplementary Table S1). Interestingly, expression of the inflammatory markers *Tnfα* and *Ccl2*, which are not rhythmic in the CD group, displayed a de novo rhythmicity upon MCDD-induced steatohepatitis (Fig. 1E, Supplementary Table S1). As expected, upon MCDD, both AST and ALT activities were higher in the MCDD fed mice than in the CD fed mice (Fig. 1F, Supplementary Table S1). We observed that both ALT and AST activity varied depending on the time of the day, ALT activity being at ZT3 than at ZT9, and AST activity displaying a 24 h oscillation. (Fig. 1F, Supplementary Table S1). In addition, the fibrosis priming markers *Tgfβ1*, *Col1α1* and *Timp1* were not rhythmic upon CD and, as expected, were upregulated upon MCDD. While *Tgfβ1* did not display circadian rhythmicity, *Col1α1* and *Timp1* both gained rhythmicity upon MCDD (Fig. 1G, Supplementary Table S1). Collectively these results demonstrate that multiple markers that have an arrhythmic expression in chow diet display an up-regulation associated with a time-of-day dependent variation in the MCDD induced NASH mouse model.

**Steatohepatitis induces discrete alterations of the liver molecular clock and abolishes *Klf10* oscillation.** To investigate whether the observed de novo rhythmicity of steatohepatitis features could be
<figure><img src="image_4.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>F</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>G</figcaption></figure>

<figure><img src="image_9.png" /></figure>

**Figure 1.** Hepatic steatosis and inflammation display circadian rhythmicity during steatohepatitis. (A) WT male mice were fed on a control diet (CD) (n=3-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and the blood and liver were sampled around the clock every 6 h. (B) Quantification of hepatic steatosis from H&E stained liver sections (%) and total liver triglyceride contents (mg/g of tissue). (C) Hepatic gene expression of *Pnpla2* and *Fsp27*. (D) quantification of inflammatory foci from the H&E staining of liver tissue section samples. (E) Gene expression of *Tnfa* and *Ccl2*. (F) Serum ALT and AST activity (IU/L). (G) Gene expression of *Tgfβ*, *Col1α1* and *Timp1*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the liver complications and related gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor) (Supplementary Table S1). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle.
*, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.

linked to a modification of the hepatic circadian clock gene network, we profiled clock and clock-controlled gene expression in the same set of mice. We observed robust oscillations in both normal and pathogenic livers for all
<figure><img src="image_4.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>E</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>F</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>G</figcaption></figure>

<figure><img src="image_10.png" /></figure>

**Figure 1.** Hepatic steatosis and inflammation display circadian rhythmicity during steatohepatitis. (A) WT male mice were fed on a control diet (CD) (n=3-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and the blood and liver were sampled around the clock every 6 h. (B) Quantification of hepatic steatosis from H&E stained liver sections (%) and total liver triglyceride contents (mg/g of tissue). (C) Hepatic gene expression of *Pnpla2* and *Fsp27*. (D) quantification of inflammatory foci from the H&E staining of liver tissue section samples. (E) Gene expression of *Tnfa* and *Ccl2*. (F) Serum ALT and AST activity (IU/L). (G) Gene expression of *Tgfβ*, *Col1α1* and *Timp1*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the liver complications and related gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor) (Supplementary Table S1). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle.
*, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.

linked to a modification of the hepatic circadian clock gene network, we profiled clock and clock-controlled gene expression in the same set of mice. We observed robust oscillations in both normal and pathogenic livers for all
analyzed core clock genes including *Bmal1*, *Rev-erbα* and β, *Per2*, *Cry1*, *Rorc* (Fig. 2A, Supplementary Table S2). However, we noticed that the phase of *Bmal1* and *Cry1* was advanced by approximately 3 h and that amplitude of *Rorc* was augmented by twofolds (Fig. 2A, Supplementary Table S2). In contrast the profile of the clock-controlled transcription factor *Dbp* was more severely altered as its amplitude was dampened by twofolds in MCDD fed mice (Fig. 2A). We also found that *Klf10* expression became arrhythmic upon this challenge (Fig. 2B, Supplementary Table S2). To address whether the observed changes in clock and clock-controlled gene expression are restricted to the liver, we also determined the expression of the same set of genes in the kidney as it has been described that steatohepatitis and development of chronic kidney disease are linked<sup>42</sup> (Fig. 3A). Interestingly, clock genes were also rhythmic in both groups of mice with discrete modifications of their phase and amplitude (Fig. 3B, Supplementary Table S3). However, among all the evaluated genes, *Klf10* was also the only gene becoming arrhythmic upon MCDD in the kidney (Fig. 3C, Supplementary Table S3). We next determined whether this loss of rhythmicity was also observed for other *Klf5* known to be linked to liver physiology and disease. As shown in Fig. 2B and Supplementary Table S2, hepatic expression of *Klf4*, previously described in the M2 polarization of Kupffer cells during NASH<sup>34</sup>, was increased upon MCDD. Similarly, *Klf6*, which was associated with NAFLD and fibrosis<sup>33,43</sup>, was also up regulated and displayed a rhythmic pattern. In contrast, expression of the *Klf10* paralog *Klf11*, which has been previously reported to be important in hepatic lipid metabolism<sup>35</sup>, was not affected with the development of steatohepatitis in liver. Collectively, these gene expression profiles suggest that the MCDD challenge has no dramatic effects on the core circadian clock dynamics. Conversely the clock output *Klf10* appears to be depending on MCDD induced pathways that overcome its physiological circadian regulation normally occurring in healthy organs.

**Lack of KLF10 does not impact on the development of MCDD induced steatosis and inflammation.** Since the hepatic circadian expression of *Klf10* was strongly impacted upon the MCDD challenge, we then determined its significance during steatohepatitis. Wild type and *Klf10*-/- mice were fed with MCDD for 4 weeks and livers were sampled at two different ZTs during the rest phase as the steatosis and inflammatory foci are more important during this time window of the 24-h cycle (ZT3 and ZT9). Following MCDD challenge, the number of lipid droplets and inflammatory foci were similar in both genotypes at both ZTs (Fig. 4A,B, Supplementary Fig. S1a, b). In accordance with the quantification of the inflammatory foci, the recruitment of granulocytes (Ly6C+/Ly6G+) and inflammatory monocytes (Ly6C<sup>high</sup>) into the liver were also similar in both genotypes at ZT3 (Fig. 4C). In addition, the hepatic expression of inflammatory markers *Tnfα* and *Ccl2*, was also similar in both genotypes at both ZTs (Fig. 4D and Supplementary Fig. S1c). This data indicates that at this stage of the MCDD induced steatohepatitis, KLF10 does not appear to play a key role in the development of hepatic steatosis and inflammation.

**Lack of KLF10 aggravates MCDD induced liver injury and TNFα induced hepatocyte death.** As liver injury is a hallmark of steatohepatitis, we compared serum ALT activity and liver TUNEL assay in WT and *Klf10*-/- MCDD fed mice. ALT activity measurements were done in samples taken during the ZT21-ZT3 time window, as this corresponds to the highest activity (Fig. 1F). *Klf10*-/- mice displayed significantly higher circulating ALT activity as compared to their WT controls suggesting that more hepatocyte death occurred in *Klf10*-/- MCDD fed mice (Fig. 5A). In addition, liver sections of *Klf10*-/- mice tended to contain more TUNEL positive nuclei than WT mice (Supplementary Fig. S2a) (evaluated at ZT3). Taken together, these results suggest that more hepatocyte death occurred in *Klf10*-/- MCDD fed mice. To determine whether KLF10 could be a protective factor against hepatocyte death, we challenged *Klf10* silenced vs control mouse primary hepatocytes with TNFα and actinomycin D (ActD). TNFα is a well-known factor mediating hepatocyte death in steatohepatitis and is strongly upregulated in MCDD mice (Figs. 1E, 4D and Supplementary Fig. S1c). Unchallenged *Klf10*-silenced primary hepatocytes did not display any change regarding the expression of inflammatory markers (*Tnfα*, *Il6* and *Ikb*) (Supplementary Fig. S2b). In response to TNFα/ActD, *Klf10*-silenced primary hepatocytes displayed reduced cell viability (Fig. 5B) associated with a higher percentage of dead cells and early apoptotic cells. (Fig. 5C). These cells also exhibited higher levels of activated (cleaved) caspase 3 (p17) in response to TNFα/ActD compared to control primary hepatocytes (Fig. 5D). These results suggest that hepatocytes lacking *Klf10* are more prone to apoptosis mediated by TNFα, which could support the aggravated liver injury observed in *Klf10*-/- mice upon MCDD challenge.

**Lack of KLF10 enhanced the priming of hepatic fibrosis induced by 4 weeks MCDD challenge.** Since KLF10 deficiency was associated with increased liver injury, we then evaluated the hepatic fibrosis. Both genotypes developed no or few hepatic fibrosis after only 4 weeks MCDD feeding as assessed by Sirius red-stained liver sections (Fig. 6A,B). The absence of collagen deposition in WT mice was expected since we have previously reported that this requires a long lasting MCDD-challenge<sup>44</sup>. Interestingly, *Klf10*-/- mice displayed significantly higher hepatic expression of fibrosis markers including *Tgfβ*, *Timp1* and *Col1α1* (Fig. 6C) compared to control animals, even before liver collagen accumulation. This suggests that these mice could be more prone to develop severe liver complications in response to long-lasting MCDD challenge by enhancing liver injury and fibrogenesis.

Discussion

Recent studies have shown that nutritional challenges including high fat and ketogenic diets can reprogram extensively and in a tissue-specific manner the mouse circadian physiology<sup>15,17-19</sup>. This reprogramming is also observed in case of diseases such as rheumatoid arthritis<sup>16</sup>. In the present work we used the MCDD experimental paradigm to study the link between circadian timing and the development of steatohepatitis in the mouse. We
analyzed core clock genes including *Bmal1*, *Rev-erbα* and β, *Per2*, *Cry1*, *Rorc* (Fig. 2A, Supplementary Table S2). However, we noticed that the phase of *Bmal1* and *Cry1* was advanced by approximately 3 h and that amplitude of *Rorc* was augmented by twofolds (Fig. 2A, Supplementary Table S2). In contrast the profile of the clock-controlled transcription factor *Dbp* was more severely altered as its amplitude was dampened by twofolds in MCDD fed mice (Fig. 2A). We also found that *Klf10* expression became arrhythmic upon this challenge (Fig. 2B, Supplementary Table S2). To address whether the observed changes in clock and clock-controlled gene expression are restricted to the liver, we also determined the expression of the same set of genes in the kidney as it has been described that steatohepatitis and development of chronic kidney disease are linked<sup>42</sup> (Fig. 3A). Interestingly, clock genes were also rhythmic in both groups of mice with discrete modifications of their phase and amplitude (Fig. 3B, Supplementary Table S3). However, among all the evaluated genes, *Klf10* was also the only gene becoming arrhythmic upon MCDD in the kidney (Fig. 3C, Supplementary Table S3). We next determined whether this loss of rhythmicity was also observed for other *Klf5* known to be linked to liver physiology and disease. As shown in Fig. 2B and Supplementary Table S2, hepatic expression of *Klf4*, previously described in the M2 polarization of Kupffer cells during NASH<sup>34</sup>, was increased upon MCDD. Similarly, *Klf6*, which was associated with NAFLD and fibrosis<sup>33,43</sup>, was also up regulated and displayed a rhythmic pattern. In contrast, expression of the *Klf10* paralog *Klf11*, which has been previously reported to be important in hepatic lipid metabolism<sup>35</sup>, was not affected with the development of steatohepatitis in liver. Collectively, these gene expression profiles suggest that the MCDD challenge has no dramatic effects on the core circadian clock dynamics. Conversely the clock output *Klf10* appears to be depending on MCDD induced pathways that overcome its physiological circadian regulation normally occurring in healthy organs.

**Lack of KLF10 does not impact on the development of MCDD induced steatosis and inflammation.** Since the hepatic circadian expression of *Klf10* was strongly impacted upon the MCDD challenge, we then determined its significance during steatohepatitis. Wild type and *Klf10*-/- mice were fed with MCDD for 4 weeks and livers were sampled at two different ZTs during the rest phase as the steatosis and inflammatory foci are more important during this time window of the 24-h cycle (ZT3 and ZT9). Following MCDD challenge, the number of lipid droplets and inflammatory foci were similar in both genotypes at both ZTs (Fig. 4A,B, Supplementary Fig. S1a, b). In accordance with the quantification of the inflammatory foci, the recruitment of granulocytes (Ly6C+/Ly6G+) and inflammatory monocytes (Ly6C<sup>high</sup>) into the liver were also similar in both genotypes at ZT3 (Fig. 4C). In addition, the hepatic expression of inflammatory markers *Tnfα* and *Ccl2*, was also similar in both genotypes at both ZTs (Fig. 4D and Supplementary Fig. S1c). This data indicates that at this stage of the MCDD induced steatohepatitis, KLF10 does not appear to play a key role in the development of hepatic steatosis and inflammation.

**Lack of KLF10 aggravates MCDD induced liver injury and TNFα induced hepatocyte death.** As liver injury is a hallmark of steatohepatitis, we compared serum ALT activity and liver TUNEL assay in WT and *Klf10*-/- MCDD fed mice. ALT activity measurements were done in samples taken during the ZT21-ZT3 time window, as this corresponds to the highest activity (Fig. 1F). *Klf10*-/- mice displayed significantly higher circulating ALT activity as compared to their WT controls suggesting that more hepatocyte death occurred in *Klf10*-/- MCDD fed mice (Fig. 5A). In addition, liver sections of *Klf10*-/- mice tended to contain more TUNEL positive nuclei than WT mice (Supplementary Fig. S2a) (evaluated at ZT3). Taken together, these results suggest that more hepatocyte death occurred in *Klf10*-/- MCDD fed mice. To determine whether KLF10 could be a protective factor against hepatocyte death, we challenged *Klf10* silenced vs control mouse primary hepatocytes with TNFα and actinomycin D (ActD). TNFα is a well-known factor mediating hepatocyte death in steatohepatitis and is strongly upregulated in MCDD mice (Figs. 1E, 4D and Supplementary Fig. S1c). Unchallenged *Klf10*-silenced primary hepatocytes did not display any change regarding the expression of inflammatory markers (*Tnfα*, *Il6* and *Ikb*) (Supplementary Fig. S2b). In response to TNFα/ActD, *Klf10*-silenced primary hepatocytes displayed reduced cell viability (Fig. 5B) associated with a higher percentage of dead cells and early apoptotic cells. (Fig. 5C). These cells also exhibited higher levels of activated (cleaved) caspase 3 (p17) in response to TNFα/ActD compared to control primary hepatocytes (Fig. 5D). These results suggest that hepatocytes lacking *Klf10* are more prone to apoptosis mediated by TNFα, which could support the aggravated liver injury observed in *Klf10*-/- mice upon MCDD challenge.

**Lack of KLF10 enhanced the priming of hepatic fibrosis induced by 4 weeks MCDD challenge.** Since KLF10 deficiency was associated with increased liver injury, we then evaluated the hepatic fibrosis. Both genotypes developed no or few hepatic fibrosis after only 4 weeks MCDD feeding as assessed by Sirius red-stained liver sections (Fig. 6A,B). The absence of collagen deposition in WT mice was expected since we have previously reported that this requires a long lasting MCDD-challenge<sup>44</sup>. Interestingly, *Klf10*-/- mice displayed significantly higher hepatic expression of fibrosis markers including *Tgfβ*, *Timp1* and *Col1α1* (Fig. 6C) compared to control animals, even before liver collagen accumulation. This suggests that these mice could be more prone to develop severe liver complications in response to long-lasting MCDD challenge by enhancing liver injury and fibrogenesis.

Discussion

Recent studies have shown that nutritional challenges including high fat and ketogenic diets can reprogram extensively and in a tissue-specific manner the mouse circadian physiology<sup>15,17-19</sup>. This reprogramming is also observed in case of diseases such as rheumatoid arthritis<sup>16</sup>. In the present work we used the MCDD experimental paradigm to study the link between circadian timing and the development of steatohepatitis in the mouse. We
<figure><img src="image_10.png" /><figcaption>**Figure 2.** Steatohepatitis in mice leads to *Klf10* loss of rhythmicity in the liver. WT male mice were fed on a control diet (CD) (n=3-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and the blood and liver were sampled around the clock every 6 h. (A) Gene expression of *Bmal1*, *Rev-erba*, *Per2*, *Cry1*, *Rev-erb-β*, *Rorc* and *Dbp* in mice fed a CD or a MCDD. (B) Gene expression of *Klf10*, *Klf11*, *Klf4* and *Klf6*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the hepatic gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor) (Supplementary Table S2). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle. *, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.</figcaption></figure>

found that the hepatic molecular clock is mostly resilient to this metabolic challenge and, that markers of steatosis and lobular inflammation, as well as serum transaminases activity displayed a diurnal rhythmicity. The finding that multiple challenges can consistently unveil pathology-related diurnal oscillations has several important basic mechanistic and translational implications.
<figure><img src="image_11.png" /><figcaption>**Figure 2.** Steatohepatitis in mice leads to *Klf10* loss of rhythmicity in the liver. WT male mice were fed on a control diet (CD) (n=3-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and the blood and liver were sampled around the clock every 6 h. (A) Gene expression of *Bmal1*, *Rev-erba*, *Per2*, *Cry1*, *Rev-erb-β*, *Rorc* and *Dbp* in mice fed a CD or a MCDD. (B) Gene expression of *Klf10*, *Klf11*, *Klf4* and *Klf6*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the hepatic gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor) (Supplementary Table S2). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle. *, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.</figcaption></figure>

found that the hepatic molecular clock is mostly resilient to this metabolic challenge and, that markers of steatosis and lobular inflammation, as well as serum transaminases activity displayed a diurnal rhythmicity. The finding that multiple challenges can consistently unveil pathology-related diurnal oscillations has several important basic mechanistic and translational implications.
<figure><img src="image_11.png" /><figcaption>**Figure 3.** Steatohepatitis in mice leads to *Klf10* loss of rhythmicity in the kidney. (A) WT male mice were fed on a control diet (CD) (n=2-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and kidneys were sampled around the clock every 6 h. (B) Gene expression of *Bmal1*, *Rev-erba*, *Per2*, *Cry1*, *Rev-erb-β*, *Rorc* and *Dbp*. (C) Gene expression of *Klf10* and *Klf11*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the renal gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor). (Supplementary Table S3). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle. *, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 3.** Steatohepatitis in mice leads to *Klf10* loss of rhythmicity in the kidney. (A) WT male mice were fed on a control diet (CD) (n=2-4 mice per ZT) or a methionine and choline deficient diet (MCDD) (n=5-6 mice per ZT) for 4 weeks and kidneys were sampled around the clock every 6 h. (B) Gene expression of *Bmal1*, *Rev-erba*, *Per2*, *Cry1*, *Rev-erb-β*, *Rorc* and *Dbp*. (C) Gene expression of *Klf10* and *Klf11*. All data are expressed as mean ± SEM. Gene expressions are normalized to *B2m* and expressed relative to the CD ZT3 level. Rhythmicity of the renal gene expressions was evaluated by nonlinear regression cosine fitting analysis (cosinor). (Supplementary Table S3). The rhythmicity of each parameter is indicated on their respective graph with the corresponding color. ZT3 values were double-plotted to complete the 24 h cycle. *, rhythmic (p < 0.05), nsr, non-significantly rhythmic, nd, not detected.</figcaption></figure>
A

MCDD

<figure><img src="image_12.png" /></figure>

B

Hepatic
steatosis

<figure><img src="image_13.png" /></figure>

C

Granulocytes

<figure><img src="image_14.png" /></figure>

D

Monocytes

<figure><img src="image_15.png" /></figure>

Tnfa

<figure><img src="image_16.png" /></figure>

Ccl2

**Figure 4.** *Klf10* deficiency does not impact the development of hepatic steatosis and inflammation during steatohepatitis in mice. WT and *Klf10*-/- male mice were fed on a methionine and choline deficient diet (MCDD) for 4 weeks and sacrificed at ZT 3. (A) Representative images, showing the presence of lipid droplets and inflammatory foci (black arrows), from H&E stained liver sections of WT and *Klf10*-/- upon MCDD. (B) Quantification of hepatic steatosis and inflammatory foci from H&E stained liver sections (n = 8–11 mice/group). (C) Number of hepatic granulocytes and monocytes assessed by flow cytometry on hepatic non parenchymal cells stained for CD45, Ly6C and Ly6G. Data represents the number of cells/g of tissue (n = 5–6 mice/group). (D) Hepatic expression of *Tnfa* and *Ccl2* (n = 8–12 mice/group). All data are expressed as mean ± SEM. The mRNA levels are normalized to *B2m* and expressed relative to the CD ZT3 level (n = 12). Statistically significance was analyzed using the Mann–Whitney test vs control group WT CD, \* p < 0.05, \*\* p < 0.0001.

The emergence of robust circadian rhythms in a pathologic setting may seem paradoxical at first glance, but this phenomenon is likely to reflect the fact that the CLOCK:BMAL1 heterodimer acts primarily as a pioneer factor regulating rhythmically chromatin accessibility to other transcription factors rather than being a necessary and sufficient driver of circadian transcription⁴⁵,⁴⁶. The extensive de novo rhythmicity induced by the high fat diet regimen was for example explained by the circadian nuclear accumulation of the lipogenic transcription factor PPARγ¹⁸,²⁰. It is therefore conceivable that transcriptional responses to acute or chronic metabolic challenges involve genes which are normally repressed or not significantly expressed but which contain CLOCK:BMAL1 binding sites in their enhancers so that their pathogenic activation, by other transcription factors such as PPARγ, CREBH, NF-κB, gets temporally constrained by the circadian clock¹⁸,⁴⁷,⁴⁸. In line with this model, the *Fsp27/Cidec* mRNA coding for a protein promoting lipid droplet formation that serves as a marker of liver steatosis displayed a high amplitude rhythm with a phase preceding the peak of droplets in MCDD fed mice, and this gene was previously reported to be bound by the CLOCK:BMAL1 heterodimer despite not being expressed in healthy liver¹⁶. Interestingly, the α and β isoforms of mouse *Fsp27* are regulated by the peroxisome proliferator-activated receptor gamma (PPARγ) and Cyclic-AMP-responsive-element-binding protein H (CREBH) respectively⁴⁹,⁵⁰, which both become rhythmic. *Fsp27* could therefore become rhythmic in the in the MCDD model through a mechanism similar to that described by Murakami et al.²⁰. An additional and not exclusive mechanism may also be attributed to a putative disruption of homeostatic microbiota rhythmicity by the MCDD model since this challenge was reported to alter the gut microbiota⁵¹. Importantly, disruption of homeostatic microbiome rhythmicity was reported to lead to genome-wide de novo oscillations in both intestine and liver¹⁹,²⁰,⁵²,⁵³. Unexpectedly, the diurnal pattern of MCDD induced steatosis was not paralleled by a similar rhythm of hepatic TG. In high fat
<figure><img src="image_13.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Inflammatory foci</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

**Figure 4.** *Klf10* deficiency does not impact the development of hepatic steatosis and inflammation during steatohepatitis in mice. WT and *Klf10*⁻/⁻ male mice were fed on a methionine and choline deficient diet (MCDD) for 4 weeks and sacrificed at ZT 3. (A) Representative images, showing the presence of lipid droplets and inflammatory foci (black arrows), from H&E stained liver sections of WT and *Klf10*⁻/⁻ upon MCDD. (B) Quantification of hepatic steatosis and inflammatory foci from H&E stained liver sections (n = 8–11 mice/group). (C) Number of hepatic granulocytes and monocytes assessed by flow cytometry on hepatic non parenchymal cells stained for CD45, Ly6C and Ly6G. Data represents the number of cells/g of tissue (n = 5–6 mice/group). (D) Hepatic expression of *Tnfa* and *Ccl2* (n = 8–12 mice/group). All data are expressed as mean ± SEM. The mRNA levels are normalized to *B2m* and expressed relative to the CD ZT3 level (n = 12). Statistically significance was analyzed using the Mann–Whitney test vs control group WT CD, \* p < 0.05, \*\* p < 0.0001.

The emergence of robust circadian rhythms in a pathologic setting may seem paradoxical at first glance, but this phenomenon is likely to reflect the fact that the CLOCK:BMAL1 heterodimer acts primarily as a pioneer factor regulating rhythmically chromatin accessibility to other transcription factors rather than being a necessary and sufficient driver of circadian transcription⁴⁵,⁴⁶. The extensive de novo rhythmicity induced by the high fat diet regimen was for example explained by the circadian nuclear accumulation of the lipogenic transcription factor PPARγ¹⁸,²⁰. It is therefore conceivable that transcriptional responses to acute or chronic metabolic challenges involve genes which are normally repressed or not significantly expressed but which contain CLOCK:BMAL1 binding sites in their enhancers so that their pathogenic activation, by other transcription factors such as PPARγ, CREBH, NF-κB, gets temporally constrained by the circadian clock¹⁸,⁴⁷,⁴⁸. In line with this model, the *Fsp27/Cidec* mRNA coding for a protein promoting lipid droplet formation that serves as a marker of liver steatosis displayed a high amplitude rhythm with a phase preceding the peak of droplets in MCDD fed mice, and this gene was previously reported to be bound by the CLOCK:BMAL1 heterodimer despite not being expressed in healthy liver¹⁶. Interestingly, the α and β isoforms of mouse *Fsp27* are regulated by the peroxisome proliferator-activated receptor gamma (PPARγ) and Cyclic-AMP-responsive-element-binding protein H (CREBH) respectively⁴⁹,⁵⁰, which both become rhythmic. *Fsp27* could therefore become rhythmic in the in the MCDD model through a mechanism similar to that described by Murakami et al.²⁰. An additional and not exclusive mechanism may also be attributed to a putative disruption of homeostatic microbiota rhythmicity by the MCDD model since this challenge was reported to alter the gut microbiota⁵¹. Importantly, disruption of homeostatic microbiome rhythmicity was reported to lead to genome-wide de novo oscillations in both intestine and liver¹⁹,²⁰,⁵²,⁵³. Unexpectedly, the diurnal pattern of MCDD induced steatosis was not paralleled by a similar rhythm of hepatic TG. In high fat
A

ALT

<figure><img src="image_17.png" /></figure>

B

Klf10

<figure><img src="image_18.png" /></figure>

C

<figure><img src="image_19.png" /></figure>

D

<figure><img src="image_20.png" /></figure>

**Figure 5.** *Klf10* deficiency aggravates liver injury during steatohepatitis in mice and hepatocyte death in vitro. (A) Serum ALT activity of WT and *Klf10*/- male mice were fed a methionine and choline deficient diet (MCDD) for 4 weeks and sacrificed at ZT 3 and 21 (n=8-20 mice/group). (B-D) Mouse primary hepatocytes were transfected with control (siCtrl) or *Klf10* (si*Klf10*) siRNA and were stimulated with TNFα (20 ng/mL) and actinomycin D (0.1 µg/mL) for 16 h (B-C) or 12 h (D). (B) Hepatocyte *Klf10* gene expression (n=6). Gene expression is normalized to *B2m* and expressed as relative expression of the untreated siCtrl condition, 72 h after transfection. Hepatocytes cell viability assessed by MTT assay (n=6). (C) Hepatocyte cell death and apoptosis assessed by flow cytometry following 7aaD and annexinV staining (n=5). (D) Western blot of cleaved CASPASE 3 fragment (17 kDa) (c-CASP3) and HSP90 protein (left) and densitometry quantification of cleaved CASPASE 3 was normalized to HSP90 in four independent experiments (right). All the data are expressed as mean ± SEM. Statistical significance was tested using the Mann-Whitney test; \* and $p < 0.05$, \*\* and $$p < 0.01$, \*\*\* and $$$p < 0.001$, vs control group (A)/cells (B, C, D) (stars) or the indicated groups (A)/cells (B, C, D) (dollar symbols).

diet challenged mice, metabolomics revealed that hepatic glycerol levels are elevated and not rhythmic¹⁷. While the rhythm of intermediated products of fatty acid esterification and TG hydrolysis (DAG and MAG) have not been well established in fatty liver, our results could also suggest that lipid droplets dynamics rather than TG
A

<figure><img src="image_19.png" /></figure>

ALT

B

<figure><img src="image_20.png" /></figure>

Klf10

C

<figure><img src="image_21.png" /></figure>

D

<figure><img src="image_22.png" /><figcaption>**Figure 5.** *Klf10* deficiency aggravates liver injury during steatohepatitis in mice and hepatocyte death in vitro. (A) Serum ALT activity of WT and *Klf10*/- male mice were fed a methionine and choline deficient diet (MCDD) for 4 weeks and sacrificed at ZT 3 and 21 (n=8-20 mice/group). (B-D) Mouse primary hepatocytes were transfected with control (siCtrl) or *Klf10* (si*Klf10*) siRNA and were stimulated with TNFα (20 ng/mL) and actinomycin D (0.1 µg/mL) for 16 h (B-C) or 12 h (D). (B) Hepatocyte *Klf10* gene expression (n=6). Gene expression is normalized to *B2m* and expressed as relative expression of the untreated siCtrl condition, 72 h after transfection. Hepatocytes cell viability assessed by MTT assay (n=6). (C) Hepatocyte cell death and apoptosis assessed by flow cytometry following 7aaD and annexinV staining (n=5). (D) Western blot of cleaved CASPASE 3 fragment (17 kDa) (c-CASP3) and HSP90 protein (left) and densitometry quantification of cleaved CASPASE 3 was normalized to HSP90 in four independent experiments (right). All the data are expressed as mean ± SEM. Statistical significance was tested using the Mann-Whitney test; \* and $p < 0.05$, \*\* and $$p < 0.01$, \*\*\* and $$$p < 0.001$, vs control group (A)/cells (B, C, D) (stars) or the indicated groups (A)/cells (B, C, D) (dollar symbols).</figcaption></figure>

diet challenged mice, metabolomics revealed that hepatic glycerol levels are elevated and not rhythmic¹⁷. While the rhythm of intermediated products of fatty acid esterification and TG hydrolysis (DAG and MAG) have not been well established in fatty liver, our results could also suggest that lipid droplets dynamics rather than TG
<figure><img src="image_21.png" /><figcaption>**Figure 6.** *Klf10* deficient mice display enhanced fibrosis priming upon MCDD challenge. (A) Representative images of liver sections from WT and *Klf10*-/- animals fed a CD or MCDD for 4 weeks (ZT3). (B) Quantification of hepatic fibrosis using the clinical fibrosis score (Ishak) (n=7-12 mice/group). (C) Hepatic gene expression of *Tgfβ1*, *Col1α1* and *Timp1* in WT and *Klf10*-/- mice sacrificed at ZT3 (n=8-12 mice/group). All the data are expressed as mean ± SEM. The mRNA levels are normalized to *B2m* and expressed as relative to the CD level ZT3 (n=12). Statistical significance was tested using the Mann-Whiney test vs the control group (stars) or the indicated groups (dollar symbols), **** p < 0.0001; ** p < 0.01; * p < 0.05.</figcaption></figure>

metabolism could be rhythmic in the MCDD hepatic steatosis model. The mechanism underlying such rhythmicity could operate at the level of droplet biogenesis and/or degradation, or alternatively be the result of interaction with clock-controlled cellular processes³⁴ such as mitochondrial metabolism⁵⁵, autophagy⁵⁶ or the unfolded protein response⁵⁷,⁵⁸. One important role of lipid droplets besides the storage of TGs, steryl esters and retinyl esters, is to safely sequester otherwise toxic lipids such as overabundant non esterified fatty acids (NEFA). This protective function against lipotoxicity is probably the reason for the abundant accumulation of lipid droplets in many disease states such as NAFLD⁵⁹. This implies that the intracellular toxicity of harmful lipids such as NEFAs is likely to be rhythmic.

Immune cells trafficking displays diurnal variations under basal conditions in multiple organs, including the liver<sup>60,61</sup>. Our data provides the evidence that there is a clear diurnal variation of inflammatory foci during the rest phase in the liver of mice fed a MCDD. Interestingly, this rhythmicity is close to what was observed in the aortic plaque during atherosclerosis<sup>62</sup>. This rhythm could be the result of a general inflammatory circadian homing code. Hepatic *Tnfa* and *Ccl2* gene expression gained rhythmicity with a peak of expression during the day (ZT2 and 4, respectively). TNFα and CCL2 play a crucial role in hepatic inflammation and leukocyte recruitment into the liver during steatohepatitis, respectively<sup>63</sup>. These observations are in line with previous studies in mouse models of inflammation mediated by *Plasmodium* infection and diet-mediated atherosclerosis<sup>62,64</sup>. In the latter, CCL2 was involved in the circadian recruitment of granulocytes and monocytes at the site of inflammation in a clock dependent fashion<sup>62</sup>.

In humans, NASH occurs in association with insulin resistance and obesity. In contrast, MCDD-fed mice lose weight, are hypoglycemic and do not have systemic insulin resistance, however they develop hepatic insulin resistance<sup>65</sup>. Our data obtained with the MCDD paradigm should stimulate further studies with models
<figure><img src="image_23.png" /><figcaption>**Figure 6.** *Klf10* deficient mice display enhanced fibrosis priming upon MCDD challenge. (A) Representative images of liver sections from WT and *Klf10*-/- animals fed a CD or MCDD for 4 weeks (ZT3). (B) Quantification of hepatic fibrosis using the clinical fibrosis score (Ishak) (n=7-12 mice/group). (C) Hepatic gene expression of *Tgfβ1*, *Col1α1* and *Timp1* in WT and *Klf10*-/- mice sacrificed at ZT3 (n=8-12 mice/group). All the data are expressed as mean ± SEM. The mRNA levels are normalized to *B2m* and expressed as relative to the CD level ZT3 (n=12). Statistical significance was tested using the Mann-Whiney test vs the control group (stars) or the indicated groups (dollar symbols), **** p < 0.0001; ** p < 0.01; *p* < 0.05.</figcaption></figure>

metabolism could be rhythmic in the MCDD hepatic steatosis model. The mechanism underlying such rhythmicity could operate at the level of droplet biogenesis and/or degradation, or alternatively be the result of interaction with clock-controlled cellular processes³⁴ such as mitochondrial metabolism⁵⁵, autophagy⁵⁶ or the unfolded protein response⁵⁷,⁵⁸. One important role of lipid droplets besides the storage of TGs, steryl esters and retinyl esters, is to safely sequester otherwise toxic lipids such as overabundant non esterified fatty acids (NEFA). This protective function against lipotoxicity is probably the reason for the abundant accumulation of lipid droplets in many disease states such as NAFLD⁵⁹. This implies that the intracellular toxicity of harmful lipids such as NEFAs is likely to be rhythmic.

Immune cells trafficking displays diurnal variations under basal conditions in multiple organs, including the liver<sup>60,61</sup>. Our data provides the evidence that there is a clear diurnal variation of inflammatory foci during the rest phase in the liver of mice fed a MCDD. Interestingly, this rhythmicity is close to what was observed in the aortic plaque during atherosclerosis<sup>62</sup>. This rhythm could be the result of a general inflammatory circadian homing code. Hepatic *Tnfa* and *Ccl2* gene expression gained rhythmicity with a peak of expression during the day (ZT2 and 4, respectively). TNFα and CCL2 play a crucial role in hepatic inflammation and leukocyte recruitment into the liver during steatohepatitis, respectively<sup>63</sup>. These observations are in line with previous studies in mouse models of inflammation mediated by *Plasmodium* infection and diet-mediated atherosclerosis<sup>62,64</sup>. In the latter, CCL2 was involved in the circadian recruitment of granulocytes and monocytes at the site of inflammation in a clock dependent fashion<sup>62</sup>.

In humans, NASH occurs in association with insulin resistance and obesity. In contrast, MCDD-fed mice lose weight, are hypoglycemic and do not have systemic insulin resistance, however they develop hepatic insulin resistance<sup>65</sup>. Our data obtained with the MCDD paradigm should stimulate further studies with models
recapitulating metabolic disorders such as obesity and insulin resistance to investigate if the de novo oscillations we observe with the MCDD challenge are also features of NASH.

We previously showed that *Klf10* is a circadian clock-controlled gene that regulates metabolic genes in non-pathogenic liver³⁰. Interestingly, this transcription factor loses rhythmicity with steatohepatitis. This could indicate that upon MCDD challenge, *Klf10* could be regulated by additional mechanisms independent of the CLOCK:BMAL1 heterodimer as previously described³⁰,⁶⁶. In addition, our data demonstrates that loss of KLF10 is associated with higher liver injury despite hepatic steatosis and inflammation similar to that seen in WT mice upon MCDD challenge. In accordance with the greater liver injury observed in *Klf10*⁻/⁻ mice, we show that these mutant mice display enhanced hepatic fibrosis priming, even before collagen accumulation. This indicates that KLF10 is required to attenuate hepatocyte death and subsequently disease progression. We provide evidence that *Klf10* protects mouse primary hepatocytes against TNFα induced apoptosis. This is in contrast to many previous studies showing that KLF10 acts as a tumor suppressor through TGF signaling by playing an important role in induction of apoptosis³⁶,⁶⁷. However, all these studies that supported the pro-apoptotic role of KLF10 were conducted in cancer cells while we used primary cells. So far the role of KLF10 in humans has been only investigated in pancreatic, breast and prostate cancer cells where it has been described to inhibit proliferation and to induce apoptosis thus having a potential tumor suppressor function⁴⁰. The role of KLF10 in human liver remains unknown and our work should stimulate studies to investigate its role in the development of NAFLDs in obese patients.

Despite a growing body of evidence supporting the relevance of circadian timing for the prevention, diagnosis and treatment of diseases, this remains mostly overlooked in the current general medical, occupational and clinical practice. Our data adds a new piece of evidence to suggest that taking into account circadian timing may be critical for both the study and the management of NAFLD. Most preclinical studies using NAFLD models are performed without any circadian timing consideration and our results suggest that this could lead to important biases in the interpretation of the data. Along the same line, as the diagnosis of NALFD patients is highly constrained by logistical and organizational considerations, its accuracy may be questioned in the light of the disease related circadian variations seen in mice if this would also hold true for humans. This remains to be evaluated, yet a rhythmicity of ALT activity has already been reported in cirrhotic patients⁶⁸. Finally, although no approved treatment is yet available for NAFLD, future trials designed to tackle this current limitation may also be improved by using a timed delivery of the assessed drugs. This is the basis for chronotherapy which already has proven to be effective for improving the tolerability and/or efficacy of treatments targeting cancer, asthma, infectious agents or inflammation⁶⁹,⁷⁰. This innovative therapeutics approach would be particularly relevant to hepatopathy as the liver is probably the most rhythmic organ.

## Materials and methods

**Animals and study design.** Animal experiment procedures were carried out in accordance with the CNRS and INSERM institutional guidelines. The local ethical committee (Comité Institutionnel d’Éthique Pour l’Animal de Laboratoire CIEPAL-AZUR, national agreement No. 28 PEA No. NCE 20-355) specifically approved this study. All the animals were housed under a 12 h light/12 h dark cycle in a temperature and humidity controlled facility. WT and *Klf10*⁻/⁻ mice in the C57BL/6j genetic background were generated by crossing of heterozygous animals. Genotyping of the mice was performed by polymerase chain reaction (PCR) of genomic DNA. For experiments, 3–6 month-old WT (Janvier-Labs), and WT / *Klf10*⁻/⁻³⁰ male mice were fed ad libitum either a methionine- and choline-deficient diet (MCDD) or control diet (CD) for 4 weeks [MCDD # E15653-947; Ctrl D # E15654-047 (SSNIFF, Soest, Germany)]. We used male mice because it has been shown that male mice are more susceptible to the progression of NAFLD in response to MCDD than female mice ⁷¹. Blood was sampled via retro-orbital puncture and mice were sacrificed by cervical dislocation. For the circadian time series, blood sampling and animal sacrifice were performed every 6 h over 24 h at Zeitgeber (ZT) 3, ZT9, ZT15 and ZT21, ZT0 being the time where light is switched on (7.00 am). Following sacrifice, liver and kidney were sampled as described in the following section. For experiments comparing the WT and *Klf10*⁻/⁻ genotypes, blood and liver were sampled at ZT3, ZT9 and ZT21.

**Tissue sampling and preparation.** Upon sampling, kidneys and pieces (50 to 100 mg) of the liver were immediately frozen in liquid nitrogen and stored at – 80 °C until analysis. Another piece of the liver (200 to 250 mg) was fixed in buffered formalin, paraffin-embedded, sectioned, and stained with Hematoxylin–Eosin. A third piece (400 to 500 mg) was used for immune cells extraction and flow cytometry analysis.

**Triglyceride extraction and quantification.** Lipids were extracted from 50 to 100 mg of frozen liver tissue using the methanol:chloroform phase extraction method. Livers were first sonicated and dissolved in methanol (Sigma-Aldrich, St Quentin Fallavier, FR). Then 2 vol of CHCl₃ were added and lysates were incubated at 4 °C overnight. Samples were rinsed with 0.05% CaCl₂ and the organic phase was collected and evaporated using speed-vacuum centrifuge. Dry pellets were resuspended in a PBS-BSA solution (5%). TG quantification was then performed using the Triglycerides FS kit (# CQN: KS, Dia Sys) according to the manufacturer recommendations.

**Real-time quantitative PCR analysis.** Real-time quantitative PCR analysis have been performed as previously described⁴⁴,⁷² and described in the Supplementary Methods. Briefly, total liver RNA was extracted from 20 to 25 mg of tissue, using the RNeasy Mini Kit (74104, Qiagen, Hilden, Germany) and treated with Turbo DNA-free DNase (AM 1907, Thermo Fisher Scientific Inc., Waltham, MA, USA) according to the manufacturer's protocol. The quantity and quality of the RNA samples were determined using the Agilent 2100 Bioanalyzer with
recapitulating metabolic disorders such as obesity and insulin resistance to investigate if the de novo oscillations we observe with the MCDD challenge are also features of NASH.

We previously showed that *Klf10* is a circadian clock-controlled gene that regulates metabolic genes in non-pathogenic liver³⁰. Interestingly, this transcription factor loses rhythmicity with steatohepatitis. This could indicate that upon MCDD challenge, *Klf10* could be regulated by additional mechanisms independent of the CLOCK:BMAL1 heterodimer as previously described³⁰,⁶⁶. In addition, our data demonstrates that loss of KLF10 is associated with higher liver injury despite hepatic steatosis and inflammation similar to that seen in WT mice upon MCDD challenge. In accordance with the greater liver injury observed in *Klf10*⁻/⁻ mice, we show that these mutant mice display enhanced hepatic fibrosis priming, even before collagen accumulation. This indicates that KLF10 is required to attenuate hepatocyte death and subsequently disease progression. We provide evidence that *Klf10* protects mouse primary hepatocytes against TNFα induced apoptosis. This is in contrast to many previous studies showing that KLF10 acts as a tumor suppressor through TGF signaling by playing an important role in induction of apoptosis³⁶,⁶⁷. However, all these studies that supported the pro-apoptotic role of KLF10 were conducted in cancer cells while we used primary cells. So far the role of KLF10 in humans has been only investigated in pancreatic, breast and prostate cancer cells where it has been described to inhibit proliferation and to induce apoptosis thus having a potential tumor suppressor function⁴⁰. The role of KLF10 in human liver remains unknown and our work should stimulate studies to investigate its role in the development of NAFLDs in obese patients.

Despite a growing body of evidence supporting the relevance of circadian timing for the prevention, diagnosis and treatment of diseases, this remains mostly overlooked in the current general medical, occupational and clinical practice. Our data adds a new piece of evidence to suggest that taking into account circadian timing may be critical for both the study and the management of NAFLD. Most preclinical studies using NAFLD models are performed without any circadian timing consideration and our results suggest that this could lead to important biases in the interpretation of the data. Along the same line, as the diagnosis of NALFD patients is highly constrained by logistical and organizational considerations, its accuracy may be questioned in the light of the disease related circadian variations seen in mice if this would also hold true for humans. This remains to be evaluated, yet a rhythmicity of ALT activity has already been reported in cirrhotic patients⁶⁸. Finally, although no approved treatment is yet available for NAFLD, future trials designed to tackle this current limitation may also be improved by using a timed delivery of the assessed drugs. This is the basis for chronotherapy which already has proven to be effective for improving the tolerability and/or efficacy of treatments targeting cancer, asthma, infectious agents or inflammation⁶⁹,⁷⁰. This innovative therapeutics approach would be particularly relevant to hepatopathy as the liver is probably the most rhythmic organ.

## Materials and methods

**Animals and study design.** Animal experiment procedures were carried out in accordance with the CNRS and INSERM institutional guidelines. The local ethical committee (Comité Institutionnel d’Éthique Pour l’Animal de Laboratoire CIEPAL-AZUR, national agreement No. 28 PEA No. NCE 20-355) specifically approved this study. All the animals were housed under a 12 h light/12 h dark cycle in a temperature and humidity controlled facility. WT and *Klf10*⁻/⁻ mice in the C57BL/6j genetic background were generated by crossing of heterozygous animals. Genotyping of the mice was performed by polymerase chain reaction (PCR) of genomic DNA. For experiments, 3–6 month-old WT (Janvier-Labs), and WT / *Klf10*⁻/⁻³⁰ male mice were fed ad libitum either a methionine- and choline-deficient diet (MCDD) or control diet (CD) for 4 weeks [MCDD # E15653-947; Ctrl D # E15654-047 (SSNIFF, Soest, Germany)]. We used male mice because it has been shown that male mice are more susceptible to the progression of NAFLD in response to MCDD than female mice ⁷¹. Blood was sampled via retro-orbital puncture and mice were sacrificed by cervical dislocation. For the circadian time series, blood sampling and animal sacrifice were performed every 6 h over 24 h at Zeitgeber (ZT) 3, ZT9, ZT15 and ZT21, ZT0 being the time where light is switched on (7.00 am). Following sacrifice, liver and kidney were sampled as described in the following section. For experiments comparing the WT and *Klf10*⁻/⁻ genotypes, blood and liver were sampled at ZT3, ZT9 and ZT21.

**Tissue sampling and preparation.** Upon sampling, kidneys and pieces (50 to 100 mg) of the liver were immediately frozen in liquid nitrogen and stored at – 80 °C until analysis. Another piece of the liver (200 to 250 mg) was fixed in buffered formalin, paraffin-embedded, sectioned, and stained with Hematoxylin–Eosin. A third piece (400 to 500 mg) was used for immune cells extraction and flow cytometry analysis.

**Triglyceride extraction and quantification.** Lipids were extracted from 50 to 100 mg of frozen liver tissue using the methanol:chloroform phase extraction method. Livers were first sonicated and dissolved in methanol (Sigma-Aldrich, St Quentin Fallavier, FR). Then 2 vol of CHCl₃ were added and lysates were incubated at 4 °C overnight. Samples were rinsed with 0.05% CaCl₂ and the organic phase was collected and evaporated using speed-vacuum centrifuge. Dry pellets were resuspended in a PBS-BSA solution (5%). TG quantification was then performed using the Triglycerides FS kit (# CQN: KS, Dia Sys) according to the manufacturer recommendations.

**Real-time quantitative PCR analysis.** Real-time quantitative PCR analysis have been performed as previously described⁴⁴,⁷² and described in the Supplementary Methods. Briefly, total liver RNA was extracted from 20 to 25 mg of tissue, using the RNeasy Mini Kit (74104, Qiagen, Hilden, Germany) and treated with Turbo DNA-free DNase (AM 1907, Thermo Fisher Scientific Inc., Waltham, MA, USA) according to the manufacturer's protocol. The quantity and quality of the RNA samples were determined using the Agilent 2100 Bioanalyzer with
RNA 6000 Nano Kit (#5067-1511, Agilent Technologies, Santa Clara, CA, USA). Total RNA (1 µg) was reverse transcribed with the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA). Real-time quantitative PCR was performed for each sample using the StepOne Plus Real-Time PCR System (Thermo Fisher Scientific Inc., Waltham, MA, USA). Gene expression was normalized to the *B2m* gene and calculated based on the comparative cycle threshold CT method (2⁻ΔΔCT). Clock genes were determined using the same methods with sybergreen detection as previously described ⁷². All primers and assays used are described in Supplementary Methods.

**Serum transaminases activity measurement.** After coagulation, serum was prepared by a centrifugation 10,000 for 5 min. Determination of serum transaminases (AST/ALT) was performed using in vitro test with pyridoxal phosphate activation on Roche/Hitachi cobas c systems (ASTPM, ALTPM, cobas, Meylan, France). Roche/Hitachi cobas c systems automatically calculated the concentration of each sample.

**Immune cells staining and flow cytometry analysis.** A fraction of the mouse liver (400 to 500 mg) was cut into very small pieces before incubation in RPMI 1640 medium containing 0.01% DNase I (Roche Diagnostics) and 0.05% collagenase type IV (Sigma-Aldrich, St Quentin Fallavier, FR) for 1 h at 37 °C under agitation. Then the pieces were crushed on a 100 µm cell strainer and washed with RPMI 1640 medium. The cellular suspension was centrifuged at 630 g for 5 min at room temperature (RT). The cell pellet was resuspended and centrifuged on a density cushion of Percoll 40% and 70% at 780 g for 20 min at RT. The NPC fraction located between 40 and 70% Percoll was collected, centrifuged at 630 g for 5 min at 4 °C and resuspended in PBS containing 3% FBS, 5 mM EDTA to perform flow cytometric analysis. The cells were surface-stained with antibodies against mouse CD45 (clone 30-F11), Ly6G (clone 1A8), Ly6C (clone AL-21) purchased from BD Biosciences (Le pont de Claix, FR) and eBioscience (Paris, France). Cells were acquired using a fluorescence activated cell sorting (FACS) Canto II flow cytometer, and the data were analyzed using FlowJo (BD Biosciences).

**Cellular models and treatments.** Primary hepatocytes were isolated as previously described ⁷³ and described in Supplementary Methods. Klf10 siRNAs (#11320001 MSS238499, Life Technologies, St Aubin, France) or control siRNA (#12935400, Life Technology, St Aubin, FR) were transfected into 200,000 (MTT) or 400,000 (western blot and flow cytometry) viable primary hepatocytes using lipofectamine RNAiMAX (#13778075, Life Technologies) according to the manufacturer's recommendations. After 48-h transfection, cells were treated with TNFα (#315-01A Preprotech, Neuilly-sur-Seine, France) (20 ng/ml) and actinomycin D (#A9415 Sigma-Aldrich, St Quentin Fallavier, FR) (0.1 µg/ml) in a William medium containing 0.5% BSA for 12 h or 16 h as indicated.

**MTT and cell death assays.** Cell survival was assessed using the MTT assay using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (#M2003-1G Sigma-Aldrich, St Quentin Fallavier, FR) according to the manufacturer's recommendations. Cell death was analyzed using a double fluorescent staining annexinV-PE and 7-AAD according to the manufacturer's instructions (AnnexinV-PE apoptosis detection kit I, BD Biosciences, Le pont-de-claix, France). Single cell fluorescence was analyzed by flow cytometry (BD Canto II) and FlowJo software (BD Biosciences).

**Immunoblotting.** Cells were solubilized in lysis buffer (20 mM Tris, pH 7.4, 150 mM NaCl, 10 mM EDTA, 150 mM NaF, 2 mM sodium orthovanadate, 10 mM pyrophosphate, protease inhibitor cocktail, and 1% Triton X-100) for 45 min at 4 °C. Lysates were cleared (14,000 rpm, 15 min). Proteins were quantified (BCA Protein assay kit, 23225 Thermo Fisher Scientific Inc, Waltham, MA, USA.), separated by SDS-PAGE and immunoblotted as previously described ⁷⁴. The proteins were probed with anti-caspase 3 (#9662, Cell Signaling) and anti-HSP90 (#4877, Cell Signaling) antibodies at 1 µg/mL as indicated. Detection was performed using electrochemio-luminescent method (Western Lightning Plus NEL105001EA) exposure to Amersham Hyperfilm ECL (#28906837 GE life sciences, Pittsburgh, PA). Densitometry analysis was performed using Fiji software. Densitometry of the detected c-casp3 fragment was normalized by the densitometry of the detected HSP90 band of the corresponding sample.

**TUNEL assay.** Liver paraffin embedded sections were deparaffinized, rehydrated and stained for 3' hydroxyl termini of DNA double strand breaks using the Apop Tag plus Peroxidase in situ Apoptosis kit (S7101, Merck, Meyzileu, FR), and nuclei were counter-colored using hematoxylin. All the steps were performed according to the manufacturer's instructions. The whole sections were then imaged using the Vectra Polaris digital slide scanner (CLS143455, Akoya Biosciences Inc, Marlborough, MA, USA). Image analysis were performed via IHC-based cell segmentation and TUNEL positive nuclei quantification using the HALO software (Incidia labs, Albuquerque, NM, USA). Percentage of positive nuclei was then normalized by the total analyzed surface.

**Statistical analysis.** Circadian rhythmicity of time series was assessed using the non-linear regression analysis based cosinor method implemented in R as described previously ³⁰,⁷². Statistical significance was determined using a boot-strap analysis to compare confidence intervals between oscillation parameters (mean level, amplitude and phase). Significance threshold was set at p < 0.05. For Single time point data, statistical significances between two groups were determined using the non-parametric Mann-Whitney test. Differences were considered significant for p < 0.05.
RNA 6000 Nano Kit (#5067-1511, Agilent Technologies, Santa Clara, CA, USA). Total RNA (1 µg) was reverse transcribed with the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA). Real-time quantitative PCR was performed for each sample using the StepOne Plus Real-Time PCR System (Thermo Fisher Scientific Inc., Waltham, MA, USA). Gene expression was normalized to the B2m gene and calculated based on the comparative cycle threshold CT method (2⁻ΔΔCT). Clock genes were determined using the same methods with sybergreen detection as previously described ⁷². All primers and assays used are described in Supplementary Methods.

**Serum transaminases activity measurement.** After coagulation, serum was prepared by a centrifugation 10,000 for 5 min. Determination of serum transaminases (AST/ALT) was performed using in vitro test with pyridoxal phosphate activation on Roche/Hitachi cobas c systems (ASTPM, ALTPM, cobas, Meylan, France). Roche/Hitachi cobas c systems automatically calculated the concentration of each sample.

**Immune cells staining and flow cytometry analysis.** A fraction of the mouse liver (400 to 500 mg) was cut into very small pieces before incubation in RPMI 1640 medium containing 0.01% DNase I (Roche Diagnostics) and 0.05% collagenase type IV (Sigma-Aldrich, St Quentin Fallavier, FR) for 1 h at 37 °C under agitation. Then the pieces were crushed on a 100 µm cell strainer and washed with RPMI 1640 medium. The cellular suspension was centrifuged at 630 g for 5 min at room temperature (RT). The cell pellet was resuspended and centrifuged on a density cushion of Percoll 40% and 70% at 780 g for 20 min at RT. The NPC fraction located between 40 and 70% Percoll was collected, centrifuged at 630 g for 5 min at 4 °C and resuspended in PBS containing 3% FBS, 5 mM EDTA to perform flow cytometric analysis. The cells were surface-stained with antibodies against mouse CD45 (clone 30-F11), Ly6G (clone 1A8), Ly6C (clone AL-21) purchased from BD Biosciences (Le pont de Claix, FR) and eBioscience (Paris, France). Cells were acquired using a fluorescence activated cell sorting (FACS) Canto II flow cytometer, and the data were analyzed using FlowJo (BD Biosciences).

**Cellular models and treatments.** Primary hepatocytes were isolated as previously described ⁷³ and described in Supplementary Methods. Klf10 siRNAs (#11320001 MSS238499, Life Technologies, St Aubin, France) or control siRNA (#12935400, Life Technology, St Aubin, FR) were transfected into 200,000 (MTT) or 400,000 (western blot and flow cytometry) viable primary hepatocytes using lipofectamine RNAiMAX (#13778075, Life Technologies) according to the manufacturer's recommendations. After 48-h transfection, cells were treated with TNFα (#315-01A Preprotech, Neuilly-sur-Seine, France) (20 ng/ml) and actinomycin D (#A9415 Sigma-Aldrich, St Quentin Fallavier, FR) (0.1 µg/ml) in a William medium containing 0.5% BSA for 12 h or 16 h as indicated.

**MTT and cell death assays.** Cell survival was assessed using the MTT assay using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (#M2003-1G Sigma-Aldrich, St Quentin Fallavier, FR) according to the manufacturer's recommendations. Cell death was analyzed using a double fluorescent staining annexinV-PE and 7-AAD according to the manufacturer's instructions (AnnexinV-PE apoptosis detection kit I, BD Biosciences, Le pont-de-claix, France). Single cell fluorescence was analyzed by flow cytometry (BD Canto II) and FlowJo software (BD Biosciences).

**Immunoblotting.** Cells were solubilized in lysis buffer (20 mM Tris, pH 7.4, 150 mM NaCl, 10 mM EDTA, 150 mM NaF, 2 mM sodium orthovanadate, 10 mM pyrophosphate, protease inhibitor cocktail, and 1% Triton X-100) for 45 min at 4 °C. Lysates were cleared (14,000 rpm, 15 min). Proteins were quantified (BCA Protein assay kit, 23225 Thermo Fisher Scientific Inc, Waltham, MA, USA.), separated by SDS-PAGE and immunoblotted as previously described ⁷⁴. The proteins were probed with anti-caspase 3 (#9662, Cell Signaling) and anti-HSP90 (#4877, Cell Signaling) antibodies at 1 µg/mL as indicated. Detection was performed using electrochemio-luminescent method (Western Lightning Plus NEL105001EA) exposure to Amersham Hyperfilm ECL (#28906837 GE life sciences, Pittsburgh, PA). Densitometry analysis was performed using Fiji software. Densitometry of the detected c-casp3 fragment was normalized by the densitometry of the detected HSP90 band of the corresponding sample.

**TUNEL assay.** Liver paraffin embedded sections were deparaffinized, rehydrated and stained for 3' hydroxyl termini of DNA double strand breaks using the Apop Tag plus Peroxidase in situ Apoptosis kit (S7101, Merck, Meyzileu, FR), and nuclei were counter-colored using hematoxylin. All the steps were performed according to the manufacturer's instructions. The whole sections were then imaged using the Vectra Polaris digital slide scanner (CLS143455, Akoya Biosciences Inc, Marlborough, MA, USA). Image analysis were performed via IHC-based cell segmentation and TUNEL positive nuclei quantification using the HALO software (Incidia labs, Albuquerque, NM, USA). Percentage of positive nuclei was then normalized by the total analyzed surface.

**Statistical analysis.** Circadian rhythmicity of time series was assessed using the non-linear regression analysis based cosinor method implemented in R as described previously ³⁰,⁷². Statistical significance was determined using a boot-strap analysis to compare confidence intervals between oscillation parameters (mean level, amplitude and phase). Significance threshold was set at p < 0.05. For Single time point data, statistical significances between two groups were determined using the non-parametric Mann-Whitney test. Differences were considered significant for p < 0.05.
Received: 8 January 2020; Accepted: 3 July 2020

Published online: 22 July 2020

References

1. Younossi, Z. M. et al. Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. *Hepatology* **64**, 73–84. https://doi.org/10.1002/hep.28431 (2016).

2. Anstee, Q. M. et al. From NASH to HCC: current concepts and future challenges. *Nat Rev Gastroenterol Hepatol* **16**, 411–428. https://doi.org/10.1038/s41575-019-0145-7 (2019).

3. Kew, M. C. Serum aminotransferase concentration as evidence of hepatocellular damage. *Lancet* **355**, 591–592. https://doi.org/10.1016/S0140-6736(99)00219-6 (2000).

4. Schwabe, R. F. & Luedde, T. Apoptosis and necroptosis in the liver: A matter of life and death. *Nat. Rev. Gastroenterol. Hepatol.* **15**, 738–752. https://doi.org/10.1038/s41575-018-0065-y (2018).

5. Gautheron, J. et al. A positive feedback loop between RIP3 and JNK controls non-alcoholic steatohepatitis. *EMBO Mol. Med.* **6**, 1062–1074. https://doi.org/10.15252/emmm.201403856 (2014).

6. Xu, B. et al. Gasdermin D plays a key role as a pyroptosis executor of non-alcoholic steatohepatitis in humans and mice. *J. Hepatol.* **68**, 773–782. https://doi.org/10.1016/j.jhep.2017.11.040 (2018).

7. Schuster, S. et al. Triggering and resolution of inflammation in NASH. *Nat. Rev. Gastroenterol. Hepatol.* **15**, 349–364. https://doi.org/10.1038/s41575-018-0009-6 (2018).

8. Panda, S. Circadian physiology of metabolism. *Science* **354**, 1008–1015. https://doi.org/10.1126/science.aah4967 (2016).

9. Zhang, R. et al. A circadian gene expression atlas in mammals: Implications for biology and medicine. *Proc. Natl. Acad. Sci. U S A* **111**, 16219–16224. https://doi.org/10.1073/pnas.1408886111 (2014).

10. Roenneberg, T. & Merrow, M. The Circadian clock and human health. *Curr. Biol.* **26**, R432-443. https://doi.org/10.1016/j.cub.2016.04.011 (2016).

11. Schernhammer, E. S. et al. Night-shift work and risk of colorectal cancer in the nurses' health study. *J. Natl. Cancer Inst.* **95**, 825–828. https://doi.org/10.1093/jnci/95.11.825 (2003).

12. Tucker, P. et al. Shiftwork and metabolic dysfunction. *Chronobiol. Int.* **29**, 549–555. https://doi.org/10.3109/07420528.2012.675259 (2012).

13. Pan, A. et al. Rotating night shift work and risk of type 2 diabetes: Two prospective cohort studies in women. *PLoS Med.* **8**, e1001141. https://doi.org/10.1371/journal.pmed.1001141 (2011).

14. Masri, S. et al. Lung Adenocarcinoma distally rewires hepatic circadian homeostasis. *Cell* **165**, 896–909. https://doi.org/10.1016/j.cell.2016.04.039 (2016).

15. 15Tognini, P. et al. Distinct Circadian signatures in liver and gut clocks revealed by ketogenic diet. *Cell Metab* **26**, 523–538 e525, https://doi.org/10.1016/j.cmet.2017.08.015 (2017).

16. Poolman, T. M. et al. Rheumatoid arthritis reprograms circadian output pathways. *Arthritis Res. Ther.* **21**, 47. https://doi.org/10.1186/s13075-019-1825-y (2019).

17. 17Dyar, K. A. et al. Atlas of Circadian metabolism reveals system-wide coordination and communication between clocks. *Cell* **174**, 1571–1585 e1511, https://doi.org/10.1016/j.cell.2018.08.042 (2018).

18. Eckel-mahan, K. L. et al. Reprogramming of the Circadian clock by nutritional challenge. *Cell* **155**, 1464–1478. https://doi.org/10.1016/j.cell.2013.11.034 (2013).

19. Leone, V. et al. Effects of diurnal variation of gut microbes and high-fat feeding on host circadian clock function and metabolism. *Cell Host Microbe* **17**, 681–689. https://doi.org/10.1016/j.chom.2015.03.006 (2015).

20. Murakami, M. et al. Gut microbiota directs PPAR γ-driven reprogramming of the liver circadian clock by nutritional challenge. *EMBO Rep.* **17**, 1292–1303. https://doi.org/10.15252/embr.201642463 (2016).

21. Gnocchi, D. et al. Circadian rhythms: A possible new player in non-alcoholic fatty liver disease pathophysiology. *J. Mol. Med. (Berl)* **97**, 741–759. https://doi.org/10.1007/s00109-019-01780-2 (2019).

22. Mazzoccoli, G., De Cosmo, S. & Mazza, T. The Biological Clock: A Pivotal Hub in Non-alcoholic Fatty Liver Disease Pathogenesis. *Front Physiol* **9**, 193. https://doi.org/10.3389/fphys.2018.00193 (2018).

23. Mukherji, A. et al. The circadian clock and liver function in health and disease. *J. Hepatol.* **71**, 200–211. https://doi.org/10.1016/j.jhep.2019.03.020 (2019).

24. Bugge, A. et al. Rev-erbalpha and Rev-erbbeta coordinately protect the circadian clock and normal metabolic function. *Genes Dev.* **26**, 657–667. https://doi.org/10.1101/gad.186858.112 (2012).

25. Delezie, J. et al. The nuclear receptor REV-ERBalpha is required for the daily balance of carbohydrate and lipid metabolism. *FASEB J.* **26**, 3321–3335. https://doi.org/10.1096/fj.12-208751 (2012).

26. Barclay, J. L. et al. High-fat diet-induced hyperinsulinemia and tissue-specific insulin resistance in Cry-deficient mice. *Am. J. Physiol. Endocrinol. Metab.* **304**, E1053-1063. https://doi.org/10.1152/ajpendo.00512.2012 (2013).

27. Rudic, R. D. et al. BMAL1 and CLOCK, two essential components of the circadian clock, are involved in glucose homeostasis. *PLoS Biol.* **2**, e377. https://doi.org/10.1371/journal.pbio.0020377 (2004).

28. Turek, F. W. et al. Obesity and metabolic syndrome in circadian clock mutant mice. *Science* **308**, 1043–1045. https://doi.org/10.1126/science.1108750 (2005).

29. Kettner, N. M. et al. Circadian homeostasis of liver metabolism suppresses hepatocarcinogenesis. *Cancer Cell* **30**, 909–924. https://doi.org/10.1016/j.ccell.2016.10.007 (2016).

30. Guillaumond, F. et al. Krüppel-like factor KLF10 is a link between the circadian clock and metabolism in liver. *Mol. Cell Biol.* **30**, 3059–3070. https://doi.org/10.1128/MCB.01141-09 (2010).

31. Jeyaraj, D. et al. Klf15 orchestrates circadian nitrogen homeostasis. *Cell Metab.* **15**, 311–323. https://doi.org/10.1016/j.cmet.2012.01.020 (2012).

32. Han, S. et al. Circadian control of bile acid synthesis by a KLF15-Fgf15 axis. *Nat. Commun.* **6**, 7231. https://doi.org/10.1038/ncomm s8231 (2015).

33. Bechmann, L. P. et al. Post-transcriptional activation of PPAR alpha by KLF6 in hepatic steatosis. *J. Hepatol.* **58**, 1000–1006. https://doi.org/10.1016/j.jhep.2013.01.020 (2013).

34. Han, Y. H. et al. RORalpha induces KLF4-mediated M2 polarization in the liver macrophages that protect against nonalcoholic steatohepatitis. *Cell Rep.* **20**, 124–135. https://doi.org/10.1016/j.celrep.2017.06.017 (2017).

35. Zhang, H. et al. Mouse KLF11 regulates hepatic lipid metabolism. *J. Hepatol.* **58**, 763–770. https://doi.org/10.1016/j.jhep.2012.11.024 (2013).

36. Ribeiro, A. et al. The transforming growth factor beta(1)-inducible transcription factor TIEG1, mediates apoptosis through oxidative stress. *Hepatology* **30**, 1490–1497. https://doi.org/10.1002/hep.510300620 (1999).

37. Venuprasad, K. et al. The E3 ubiquitin ligase Itch regulates expression of transcription factor Foxp3 and airway inflammation by enhancing the function of transcription factor TIEG1. *Nat. Immunol.* **9**, 245–253. https://doi.org/10.1038/ni1564 (2008).
Received: 8 January 2020; Accepted: 3 July 2020

Published online: 22 July 2020

References

1. Younossi, Z. M. et al. Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. *Hepatology* **64**, 73–84. https://doi.org/10.1002/hep.28431 (2016).

2. Anstee, Q. M. et al. From NASH to HCC: current concepts and future challenges. *Nat Rev Gastroenterol Hepatol* **16**, 411–428. https://doi.org/10.1038/s41575-019-0145-7 (2019).

3. Kew, M. C. Serum aminotransferase concentration as evidence of hepatocellular damage. *Lancet* **355**, 591–592. https://doi.org/10.1016/S0140-6736(99)00219-6 (2000).

4. Schwabe, R. F. & Luedde, T. Apoptosis and necroptosis in the liver: A matter of life and death. *Nat. Rev. Gastroenterol. Hepatol.* **15**, 738–752. https://doi.org/10.1038/s41575-018-0065-y (2018).

5. Gautheron, J. et al. A positive feedback loop between RIP3 and JNK controls non-alcoholic steatohepatitis. *EMBO Mol. Med.* **6**, 1062–1074. https://doi.org/10.15252/emmm.201403856 (2014).

6. Xu, B. et al. Gasdermin D plays a key role as a pyroptosis executor of non-alcoholic steatohepatitis in humans and mice. *J. Hepatol.* **68**, 773–782. https://doi.org/10.1016/j.jhep.2017.11.040 (2018).

7. Schuster, S. et al. Triggering and resolution of inflammation in NASH. *Nat. Rev. Gastroenterol. Hepatol.* **15**, 349–364. https://doi.org/10.1038/s41575-018-0009-6 (2018).

8. Panda, S. Circadian physiology of metabolism. *Science* **354**, 1008–1015. https://doi.org/10.1126/science.aah4967 (2016).

9. Zhang, R. et al. A circadian gene expression atlas in mammals: Implications for biology and medicine. *Proc. Natl. Acad. Sci. U S A* **111**, 16219–16224. https://doi.org/10.1073/pnas.1408886111 (2014).

10. Roenneberg, T. & Merrow, M. The Circadian clock and human health. *Curr. Biol.* **26**, R432-443. https://doi.org/10.1016/j.cub.2016.04.011 (2016).

11. Schernhammer, E. S. et al. Night-shift work and risk of colorectal cancer in the nurses' health study. *J. Natl. Cancer Inst.* **95**, 825–828. https://doi.org/10.1093/jnci/95.11.825 (2003).

12. Tucker, P. et al. Shiftwork and metabolic dysfunction. *Chronobiol. Int.* **29**, 549–555. https://doi.org/10.3109/07420528.2012.675259 (2012).

13. Pan, A. et al. Rotating night shift work and risk of type 2 diabetes: Two prospective cohort studies in women. *PLoS Med.* **8**, e1001141. https://doi.org/10.1371/journal.pmed.1001141 (2011).

14. Masri, S. et al. Lung Adenocarcinoma distally rewires hepatic circadian homeostasis. *Cell* **165**, 896–909. https://doi.org/10.1016/j.cell.2016.04.039 (2016).

15. 15Tognini, P. et al. Distinct Circadian signatures in liver and gut clocks revealed by ketogenic diet. *Cell Metab* **26**, 523–538 e525, https://doi.org/10.1016/j.cmet.2017.08.015 (2017).

16. Poolman, T. M. et al. Rheumatoid arthritis reprograms circadian output pathways. *Arthritis Res. Ther.* **21**, 47. https://doi.org/10.1186/s13075-019-1825-y (2019).

17. 17Dyar, K. A. et al. Atlas of Circadian metabolism reveals system-wide coordination and communication between clocks. *Cell* **174**, 1571–1585 e1511, https://doi.org/10.1016/j.cell.2018.08.042 (2018).

18. Eckel-mahan, K. L. et al. Reprogramming of the Circadian clock by nutritional challenge. *Cell* **155**, 1464–1478. https://doi.org/10.1016/j.cell.2013.11.034 (2013).

19. Leone, V. et al. Effects of diurnal variation of gut microbes and high-fat feeding on host circadian clock function and metabolism. *Cell Host Microbe* **17**, 681–689. https://doi.org/10.1016/j.chom.2015.03.006 (2015).

20. Murakami, M. et al. Gut microbiota directs PPAR γ-driven reprogramming of the liver circadian clock by nutritional challenge. *EMBO Rep.* **17**, 1292–1303. https://doi.org/10.15252/embr.201642463 (2016).

21. Gnocchi, D. et al. Circadian rhythms: A possible new player in non-alcoholic fatty liver disease pathophysiology. *J. Mol. Med. (Berl)* **97**, 741–759. https://doi.org/10.1007/s00109-019-01780-2 (2019).

22. Mazzoccoli, G., De Cosmo, S. & Mazza, T. The Biological Clock: A Pivotal Hub in Non-alcoholic Fatty Liver Disease Pathogenesis. *Front Physiol* **9**, 193. https://doi.org/10.3389/fphys.2018.00193 (2018).

23. Mukherji, A. et al. The circadian clock and liver function in health and disease. *J. Hepatol.* **71**, 200–211. https://doi.org/10.1016/j.jhep.2019.03.020 (2019).

24. Bugge, A. et al. Rev-erbalpha and Rev-erbbeta coordinately protect the circadian clock and normal metabolic function. *Genes Dev.* **26**, 657–667. https://doi.org/10.1101/gad.186858.112 (2012).

25. Delezie, J. et al. The nuclear receptor REV-ERBalpha is required for the daily balance of carbohydrate and lipid metabolism. *FASEB J.* **26**, 3321–3335. https://doi.org/10.1096/fj.12-208751 (2012).

26. Barclay, J. L. et al. High-fat diet-induced hyperinsulinemia and tissue-specific insulin resistance in Cry-deficient mice. *Am. J. Physiol. Endocrinol. Metab.* **304**, E1053-1063. https://doi.org/10.1152/ajpendo.00512.2012 (2013).

27. Rudic, R. D. et al. BMAL1 and CLOCK, two essential components of the circadian clock, are involved in glucose homeostasis. *PLoS Biol.* **2**, e377. https://doi.org/10.1371/journal.pbio.0020377 (2004).

28. Turek, F. W. et al. Obesity and metabolic syndrome in circadian clock mutant mice. *Science* **308**, 1043–1045. https://doi.org/10.1126/science.1108750 (2005).

29. Kettner, N. M. et al. Circadian homeostasis of liver metabolism suppresses hepatocarcinogenesis. *Cancer Cell* **30**, 909–924. https://doi.org/10.1016/j.ccell.2016.10.007 (2016).

30. Guillaumond, F. et al. Krüppel-like factor KLF10 is a link between the circadian clock and metabolism in liver. *Mol. Cell Biol.* **30**, 3059–3070. https://doi.org/10.1128/MCB.01141-09 (2010).

31. Jeyaraj, D. et al. Klf15 orchestrates circadian nitrogen homeostasis. *Cell Metab.* **15**, 311–323. https://doi.org/10.1016/j.cmet.2012.01.020 (2012).

32. Han, S. et al. Circadian control of bile acid synthesis by a KLF15-Fgf15 axis. *Nat. Commun.* **6**, 7231. https://doi.org/10.1038/ncomm s8231 (2015).

33. Bechmann, L. P. et al. Post-transcriptional activation of PPAR alpha by KLF6 in hepatic steatosis. *J. Hepatol.* **58**, 1000–1006. https://doi.org/10.1016/j.jhep.2013.01.020 (2013).

34. Han, Y. H. et al. RORalpha induces KLF4-mediated M2 polarization in the liver macrophages that protect against nonalcoholic steatohepatitis. *Cell Rep.* **20**, 124–135. https://doi.org/10.1016/j.celrep.2017.06.017 (2017).

35. Zhang, H. et al. Mouse KLF11 regulates hepatic lipid metabolism. *J. Hepatol.* **58**, 763–770. https://doi.org/10.1016/j.jhep.2012.11.024 (2013).

36. Ribeiro, A. et al. The transforming growth factor beta(1)-inducible transcription factor TIEG1, mediates apoptosis through oxidative stress. *Hepatology* **30**, 1490–1497. https://doi.org/10.1002/hep.510300620 (1999).

37. Venuprasad, K. et al. The E3 ubiquitin ligase Itch regulates expression of transcription factor Foxp3 and airway inflammation by enhancing the function of transcription factor TIEG1. *Nat. Immunol.* **9**, 245–253. https://doi.org/10.1038/ni1564 (2008).
38. Papadakis, K. A. et al. Kruppel-like factor KLF10 deficiency predisposes to colitis through colonic macrophage dysregulation. *Am. J. Physiol. Gastrointest. Liver Physiol.* **309**, G900-909. https://doi.org/10.1152/ajpgi.00309.2015 (2015).

39. Kim, J. K. et al. Progression of diet induced nonalcoholic steatohepatitis is accompanied by increased expression of Kruppel-like-factor 10 in mice. *J. Transl. Med.* **12**, 186. https://doi.org/10.1186/1479-5876-12-186 (2014).

40. Subramaniam, M. et al. Functional role of KLF10 in multiple disease processes. *BioFactors* **36**, 8–18. https://doi.org/10.1002/biof.67 (2010).

41. Aviram, R. et al. Lipidomics analyses reveal temporal and spatial lipid organization and uncover daily oscillations in intracellular organelles. *Mol. Cell* **62**, 636–648. https://doi.org/10.1016/j.molcel.2016.04.002 (2016).

42. Targher, G. & Byrne, C. D. Non-alcoholic fatty liver disease: An emerging driving force in chronic kidney disease. *Nat. Rev. Nephrol.* **13**, 297–310. https://doi.org/10.1038/nrneph.2017.16 (2017).

43. 43Miele, L. et al. The Kruppel-like factor 6 genotype is associated with fibrosis in nonalcoholic fatty liver disease. *Gastroenterology* **135**, 282–291 e281, https://doi.org/10.1053/j.gastro.2008.04.004 (2008).

44. Patouraux, S. et al. CD44 is a key player in non-alcoholic steatohepatitis. *J Hepatol* **67**, 328–338. https://doi.org/10.1016/j.jhep.2017.03.003 (2017).

45. Trott, A. J. & Menet, J. S. Regulation of circadian clock transcriptional output by CLOCK:BMAL1. *PLoS Genet* **14**, e1007156. https://doi.org/10.1371/journal.pgen.1007156 (2018).

46. Menet, J. S., Pescatore, S. & Rosbash, M. CLOCK:BMAL1 is a pioneer-like transcription factor. *Genes Dev.* **28**, 8–13. https://doi.org/10.1101/gad.228536.113 (2014).

47. Zheng, Z. et al. CREBH couples circadian clock with hepatic lipid metabolism. *Diabetes* **65**, 3369–3383. https://doi.org/10.2337/db16-0298 (2016).

48. Spengler, M. L. et al. Core circadian protein CLOCK is a positive regulator of NF-kappaB-mediated transcription. *Proc Natl Acad Sci U S A* **109**, E2457-2465. https://doi.org/10.1073/pnas.1206274109 (2012).

49. Xu, X. et al. Transcriptional activation of Fsp27 by the liver-enriched transcription factor CREBH promotes lipid droplet growth and hepatic steatosis. *Hepatology* **61**, 857–869. https://doi.org/10.1002/hep.27371 (2015).

50. Matsusue, K. et al. Hepatic steatosis in leptin-deficient mice is promoted by the PPARgamma target gene Fsp27. *Cell Metab.* **7**, 302–311. https://doi.org/10.1016/j.cmet.2008.03.003 (2008).

51. Ye, J. Z. et al. Dynamic alterations in the gut microbiota and metabolome during the development of methionine-choline-deficient diet-induced nonalcoholic steatohepatitis. *World J. Gastroenterol.* **24**, 2468–2481. https://doi.org/10.3748/wj.v24.i23.2468 (2018).

52. Thaiss, C. A., Nobs, S. P. & Elinav, E. NFIL-trating the host circadian rhythm-microbes fine-tune the epithelial clock. *Cell Metab.* **26**, 699–700. https://doi.org/10.1016/j.cmet.2017.10.014 (2017).

53. 53Weger, B. D. et al. The mouse microbiome is required for sex-specific diurnal rhythms of gene expression and metabolism. *Cell Metab* **29**, 362–382 e368, https://doi.org/10.1016/j.cmet.2018.09.023 (2019).

54. Chaix, A., Zarrinpar, A. & Panda, S. The circadian coordination of cell biology. *J Cell Biol* **215**, 15–25. https://doi.org/10.1083/jcb.201603076 (2016).

55. Jacobi, D. et al. Hepatic Bmal1 regulates rhythmic mitochondrial dynamics and promotes metabolic fitness. *Cell Metab.* **22**, 709–720. https://doi.org/10.1016/j.cmet.2015.08.006 (2015).

56. Ma, D., Panda, S. & Lin, J. D. Temporal orchestration of circadian autophagy rhythm by C/EBPbeta. *EMBO J* **30**, 4642–4651. https://doi.org/10.1038/embj.2011.322 (2011).

57. Bu, Y. et al. A PERK-miR-211 axis suppresses circadian regulators and protein synthesis to promote cancer cell survival. *Nat Cell Biol* **20**, 104–115. https://doi.org/10.1038/s41556-017-0006-y (2018).

58. 58Zhu, B. et al. A cell-autonomous mammalian 12 hr clock coordinates metabolic and stress rhythms. *Cell Metab* **25**, 1305–1319 e1309, https://doi.org/10.1016/j.cmet.2017.05.004 (2017).

59. Gluchowski, N. L. et al. Lipid droplets and liver disease: From basic biology to clinical implications. *Nat Rev Gastroenterol Hepatol* **14**, 343–355. https://doi.org/10.1038/nrgastro.2017.32 (2017).

60. Godinho-Silva, C. et al. Light-entrained and brain-tuned circadian circuits regulate ILC3s and gut homeostasis. *Nature* **574**, 254–258. https://doi.org/10.1038/s41586-019-1579-3 (2019).

61. 61He, W. et al. Circadian expression of migratory factors establishes lineage-specific signatures that guide the homing of leukocyte subsets to tissues. *Immunity* **49**, 1175-1190 e1177, https://doi.org/10.1016/j.immuni.2018.10.007 (2018).

62. 62Winter, C. et al. Chrono-pharmacological targeting of the CCL2-CCR2 axis ameliorates atherosclerosis. *Cell Metab* **28**, 175–182 e175, https://doi.org/10.1016/j.cmet.2018.05.002 (2018).

63. Baeck, C. et al. Pharmacological inhibition of the chemokine CCL2 (MCP-1) diminishes liver macrophage infiltration and steatohepatitis in chronic hepatic injury. *Gut* **61**, 416–426. https://doi.org/10.1136/gutjnl-2011-300304 (2012).

64. 64Hirako, I. C. et al. Daily Rhythms of TNFalpha expression and food intake regulate synchrony of plasmodium stages with the host circadian cycle. *Cell Host Microbe* **23**, 796–808 e796, https://doi.org/10.1016/j.chom.2018.04.016 (2018).

65. Leclercq, I. A. et al. Intrahepatic insulin resistance in a murine model of steatohepatitis: effect of PPARgamma agonist pioglitazone. *Lab Invest* **87**, 56–65. https://doi.org/10.1038/labinvest.3700489 (2007).

66. Iizuka, K., Takeda, J. & Horikawa, Y. Kruppel-like factor-10 is directly regulated by carbohydrate response element-binding protein in rat primary hepatocytes. *Biochem Biophys Res Commun* **412**, 638–643. https://doi.org/10.1016/j.bbrc.2011.08.016 (2011).

67. Jiang, L. et al. Transactivation of the TIEG1 confers growth inhibition of transforming growth factor-beta-susceptible hepatocellular carcinoma cells. *World J Gastroenterol* **18**, 2035–2042. https://doi.org/10.3748/wj.v18.i17.2035 (2012).

68. Cordoba, J. et al. Diurnal variation of serum alanine transaminase activity in chronic liver disease. *Hepatology* **28**, 1724–1725. https://doi.org/10.1002/hep.510280640 (1998).

69. Ruben, M. D. et al. Dosing time matters. *Science* **365**, 547–549. https://doi.org/10.1126/science.aax7621 (2019).

70. Cederroth, C. R. et al. Medicine in the fourth dimension. *Cell Metab.* **30**, 238–250. https://doi.org/10.1016/j.cmet.2019.06.019 (2019).

71. Lee, Y. H. et al. Sex-specific metabolic interactions between liver and adipose tissue in MCD diet-induced non-alcoholic fatty liver disease. *Oncotarget* **7**, 46959–46971. https://doi.org/10.18632/oncotarget.10506 (2016).

72. Feillet, C. et al. Sexual dimorphism in circadian physiology is altered in LXRalpha deficient mice. *PLoS ONE* **11**, e0150665. https://doi.org/10.1371/journal.pone.0150665 (2016).

73. Patouraux, S. et al. Osteopontin deficiency aggravates hepatic injury induced by ischemia-reperfusion in mice. *Cell Death Dis* **5**, e1208. https://doi.org/10.1038/cddis.2014.174 (2014).

74. Sans, A. et al. The differential expression of cide family members is associated with Nafld progression from steatosis to steatohepatitis. *Sci. Rep.* **9**, 7501. https://doi.org/10.1038/s41598-019-43928-7 (2019).

## Acknowledgements

We thank Angelo Bernabeo and the technical staff from the iBV animal facility, Samah Rekima from the iBV histology core facility, Frederic Larbret from the C3M flow cytometry core facility and Marie Irondelle from the C3M imaging core facility for their skillful assistance and support. This work was supported by grants from INSERM (France), charities (Association Française pour l'Etude du Foie (AFEF) to PG). This work was also
38. Papadakis, K. A. et al. Kruppel-like factor KLF10 deficiency predisposes to colitis through colonic macrophage dysregulation. *Am. J. Physiol. Gastrointest. Liver Physiol.* **309**, G900-909. https://doi.org/10.1152/ajpgi.00309.2015 (2015).

39. Kim, J. K. et al. Progression of diet induced nonalcoholic steatohepatitis is accompanied by increased expression of Kruppel-like-factor 10 in mice. *J. Transl. Med.* **12**, 186. https://doi.org/10.1186/1479-5876-12-186 (2014).

40. Subramaniam, M. et al. Functional role of KLF10 in multiple disease processes. *BioFactors* **36**, 8–18. https://doi.org/10.1002/biof.67 (2010).

41. Aviram, R. et al. Lipidomics analyses reveal temporal and spatial lipid organization and uncover daily oscillations in intracellular organelles. *Mol. Cell* **62**, 636–648. https://doi.org/10.1016/j.molcel.2016.04.002 (2016).

42. Targher, G. & Byrne, C. D. Non-alcoholic fatty liver disease: An emerging driving force in chronic kidney disease. *Nat. Rev. Nephrol.* **13**, 297–310. https://doi.org/10.1038/nrneph.2017.16 (2017).

43. 43Miele, L. et al. The Kruppel-like factor 6 genotype is associated with fibrosis in nonalcoholic fatty liver disease. *Gastroenterology* **135**, 282–291 e281, https://doi.org/10.1053/j.gastro.2008.04.004 (2008).

44. Patouraux, S. et al. CD44 is a key player in non-alcoholic steatohepatitis. *J Hepatol* **67**, 328–338. https://doi.org/10.1016/j.jhep.2017.03.003 (2017).

45. Trott, A. J. & Menet, J. S. Regulation of circadian clock transcriptional output by CLOCK:BMAL1. *PLoS Genet* **14**, e1007156. https://doi.org/10.1371/journal.pgen.1007156 (2018).

46. Menet, J. S., Pescatore, S. & Rosbash, M. CLOCK:BMAL1 is a pioneer-like transcription factor. *Genes Dev.* **28**, 8–13. https://doi.org/10.1101/gad.228536.113 (2014).

47. Zheng, Z. et al. CREBH couples circadian clock with hepatic lipid metabolism. *Diabetes* **65**, 3369–3383. https://doi.org/10.2337/db16-0298 (2016).

48. Spengler, M. L. et al. Core circadian protein CLOCK is a positive regulator of NF-kappaB-mediated transcription. *Proc Natl Acad Sci U S A* **109**, E2457-2465. https://doi.org/10.1073/pnas.1206274109 (2012).

49. Xu, X. et al. Transcriptional activation of Fsp27 by the liver-enriched transcription factor CREBH promotes lipid droplet growth and hepatic steatosis. *Hepatology* **61**, 857–869. https://doi.org/10.1002/hep.27371 (2015).

50. Matsusue, K. et al. Hepatic steatosis in leptin-deficient mice is promoted by the PPARgamma target gene Fsp27. *Cell Metab.* **7**, 302–311. https://doi.org/10.1016/j.cmet.2008.03.003 (2008).

51. Ye, J. Z. et al. Dynamic alterations in the gut microbiota and metabolome during the development of methionine-choline-deficient diet-induced nonalcoholic steatohepatitis. *World J. Gastroenterol.* **24**, 2468–2481. https://doi.org/10.3748/wj.v24.i23.2468 (2018).

52. Thaiss, C. A., Nobs, S. P. & Elinav, E. NFIL-trating the host circadian rhythm-microbes fine-tune the epithelial clock. *Cell Metab.* **26**, 699–700. https://doi.org/10.1016/j.cmet.2017.10.014 (2017).

53. 53Weger, B. D. et al. The mouse microbiome is required for sex-specific diurnal rhythms of gene expression and metabolism. *Cell Metab* **29**, 362–382 e368, https://doi.org/10.1016/j.cmet.2018.09.023 (2019).

54. Chaix, A., Zarrinpar, A. & Panda, S. The circadian coordination of cell biology. *J Cell Biol* **215**, 15–25. https://doi.org/10.1083/jcb.201603076 (2016).

55. Jacobi, D. et al. Hepatic Bmal1 regulates rhythmic mitochondrial dynamics and promotes metabolic fitness. *Cell Metab.* **22**, 709–720. https://doi.org/10.1016/j.cmet.2015.08.006 (2015).

56. Ma, D., Panda, S. & Lin, J. D. Temporal orchestration of circadian autophagy rhythm by C/EBPbeta. *EMBO J* **30**, 4642–4651. https://doi.org/10.1038/embj.2011.322 (2011).

57. Bu, Y. et al. A PERK-miR-211 axis suppresses circadian regulators and protein synthesis to promote cancer cell survival. *Nat Cell Biol* **20**, 104–115. https://doi.org/10.1038/s41556-017-0006-y (2018).

58. 58Zhu, B. et al. A cell-autonomous mammalian 12 hr clock coordinates metabolic and stress rhythms. *Cell Metab* **25**, 1305–1319 e1309, https://doi.org/10.1016/j.cmet.2017.05.004 (2017).

59. Gluchowski, N. L. et al. Lipid droplets and liver disease: From basic biology to clinical implications. *Nat Rev Gastroenterol Hepatol* **14**, 343–355. https://doi.org/10.1038/nrgastro.2017.32 (2017).

60. Godinho-Silva, C. et al. Light-entrained and brain-tuned circadian circuits regulate ILC3s and gut homeostasis. *Nature* **574**, 254–258. https://doi.org/10.1038/s41586-019-1579-3 (2019).

61. 61He, W. et al. Circadian expression of migratory factors establishes lineage-specific signatures that guide the homing of leukocyte subsets to tissues. *Immunity* **49**, 1175–1190 e1177, https://doi.org/10.1016/j.immuni.2018.10.007 (2018).

62. 62Winter, C. et al. Chrono-pharmacological targeting of the CCL2-CCR2 axis ameliorates atherosclerosis. *Cell Metab* **28**, 175–182 e175, https://doi.org/10.1016/j.cmet.2018.05.002 (2018).

63. Baeck, C. et al. Pharmacological inhibition of the chemokine CCL2 (MCP-1) diminishes liver macrophage infiltration and steatohepatitis in chronic hepatic injury. *Gut* **61**, 416–426. https://doi.org/10.1136/gutjnl-2011-300304 (2012).

64. 64Hirako, I. C. et al. Daily Rhythms of TNFalpha expression and food intake regulate synchrony of plasmodium stages with the host circadian cycle. *Cell Host Microbe* **23**, 796–808 e796, https://doi.org/10.1016/j.chom.2018.04.016 (2018).

65. Leclercq, I. A. et al. Intrahepatic insulin resistance in a murine model of steatohepatitis: effect of PPARgamma agonist pioglitazone. *Lab Invest* **87**, 56–65. https://doi.org/10.1038/labinvest.3700489 (2007).

66. Iizuka, K., Takeda, J. & Horikawa, Y. Kruppel-like factor-10 is directly regulated by carbohydrate response element-binding protein in rat primary hepatocytes. *Biochem Biophys Res Commun* **412**, 638–643. https://doi.org/10.1016/j.bbrc.2011.08.016 (2011).

67. Jiang, L. et al. Transactivation of the TIEG1 confers growth inhibition of transforming growth factor-beta-susceptible hepatocellular carcinoma cells. *World J Gastroenterol* **18**, 2035–2042. https://doi.org/10.3748/wj.v18.i17.2035 (2012).

68. Cordoba, J. et al. Diurnal variation of serum alanine transaminase activity in chronic liver disease. *Hepatology* **28**, 1724–1725. https://doi.org/10.1002/hep.510280640 (1998).

69. Ruben, M. D. et al. Dosing time matters. *Science* **365**, 547–549. https://doi.org/10.1126/science.aax7621 (2019).

70. Cederroth, C. R. et al. Medicine in the fourth dimension. *Cell Metab.* **30**, 238–250. https://doi.org/10.1016/j.cmet.2019.06.019 (2019).

71. Lee, Y. H. et al. Sex-specific metabolic interactions between liver and adipose tissue in MCD diet-induced non-alcoholic fatty liver disease. *Oncotarget* **7**, 46959–46971. https://doi.org/10.18632/oncotarget.10506 (2016).

72. Feillet, C. et al. Sexual dimorphism in circadian physiology is altered in LXRalpha deficient mice. *PLoS ONE* **11**, e0150665. https://doi.org/10.1371/journal.pone.0150665 (2016).

73. Patouraux, S. et al. Osteopontin deficiency aggravates hepatic injury induced by ischemia-reperfusion in mice. *Cell Death Dis* **5**, e1208. https://doi.org/10.1038/cddis.2014.174 (2014).

74. Sans, A. et al. The differential expression of cide family members is associated with Nafld progression from steatosis to steatohepatitis. *Sci. Rep.* **9**, 7501. https://doi.org/10.1038/s41598-019-43928-7 (2019).

## Acknowledgements

We thank Angelo Bernabeo and the technical staff from the iBV animal facility, Samah Rekima from the iBV histology core facility, Frederic Larbret from the C3M flow cytometry core facility and Marie Irondelle from the C3M imaging core facility for their skillful assistance and support. This work was supported by grants from INSERM (France), charities (Association Française pour l'Etude du Foie (AFEF) to PG). This work was also
funded by the French Government (National Research Agency, ANR): #ANR-18-CE14-0019-02 and through the "Investments for the Future" LABEX SIGNALIFE (#ANR-11-LABX-0028-01) and the UCAJEDI Investments in the Future project (#ANR-15-IDEX-01). This study was supported in part by research funding from the ANR HyClock and iCycle projects (ANR-14-CE09-0011 and ANR-16-CE0016) the Canceropôle PACA, Institut National du Cancer and Région Sud.

**Author contributions**

P.S.L., F.D., P.G. and M.T. designed the research, analyzed and interpreted the data and wrote the paper. P.S.L., D.R., S.P., S.G., S.B., A.G.C., A.A.R., and C.L. performed the experiments. M.S. generated the mutant mice. A.T. contributed to histology analysis. All the other authors edited and approved the final submitted draft.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-69085-w.

**Correspondence** and requests for materials should be addressed to P.G. or M.T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_22.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020
funded by the French Government (National Research Agency, ANR): #ANR-18-CE14-0019-02 and through the "Investments for the Future" LABEX SIGNALIFE (#ANR-11-LABX-0028-01) and the UCAJEDI Investments in the Future project (#ANR-15-IDEX-01). This study was supported in part by research funding from the ANR HyClock and iCycle projects (ANR-14-CE09-0011 and ANR-16-CE0016) the Canceropôle PACA, Institut National du Cancer and Région Sud.

**Author contributions**

P.S.L., F.D., P.G. and M.T. designed the research, analyzed and interpreted the data and wrote the paper. P.S.L., D.R., S.P., S.G., S.B., A.G.C., A.A.R., and C.L. performed the experiments. M.S. generated the mutant mice. A.T. contributed to histology analysis. All the other authors edited and approved the final submitted draft.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-69085-w.

**Correspondence** and requests for materials should be addressed to P.G. or M.T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_24.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020