CDKN2A/p16INK4a suppresses hepatic fatty acid oxidation through the AMPK a 2-SIRT1-PPAR a signaling pathway

Mechanisms Tumorigenesis and Target

In addition to their well-known role in the control of cellular proliferation and cancer, cell cycle regulators are increasingly identified as important metabolic modulators. Several GWAS have identified SNPs near CDKN2A, the locus encoding for p16INK4a (p16), associated with elevated risk for cardiovascular diseases and type-2 diabetes development, two pathologies associated with impaired hepatic lipid metabolism. Although p16 was recently shown to control hepatic glucose homeostasis, it is unknown whether p16 also controls hepatic lipid metabolism. Using a combination of in vivo and in vitro approaches, we found that p16 modulates fasting-induced hepatic fatty acid oxidation (FAO) and lipid droplet accumulation. In primary hepatocytes, p16-deficiency was associated with elevated expression of genes involved in fatty acid catabolism. These transcriptional changes led to increased FAO and were associated with enhanced activation of PPARa through a mechanism requiring the catalytic AMPKa2 subunit and SIRT1, two known activators of PPARa. By contrast, p16 overexpression was associated with triglyceride accumulation and increased lipid droplet numbers in vitro, and decreased ketogenesis and hepatic mitochondrial activity in vivo. Finally, gene expression analysis of liver samples from obese patients revealed a negative correlation between CDKN2A expression and PPARA and its target genes. Our findings demonstrate that p16 represses hepatic lipid catabolism during fasting and may thus participate in the preservation of metabolic flexibility.
Cell cycle regulators have been extensively studied in the context of proliferation, cancer development, and aging (1). Progression through the cell cycle requires specific metabolic programs for synthesis of cellular building blocks or ATP production (2). Interestingly, recent work has shown that several cell cycle regulators also modulate metabolism in nonproliferative cells, suggesting new physiological functions of this large family of proteins (3).
P16INK4a (p16) is a cyclin-dependent kinase inhibitor that blocks activation of E2F transcription factors via inhibition of CDK4/6 (4). Interestingly, several GWAS have identified single nucleotide polymorphism near CDKN2A, the locus encoding for p16, as associated with elevated risk for cardiovascular disease and type-2 diabetes (T2D) development (5). In line, we recently reported that p16-deficient mice display elevated hepatic gluconeogenesis during fasting due to activation of a cascade involving CDK4, PKA, CREB, and PGC1a in hepatocytes, suggesting an important role for p16 in metabolic control (6). Interestingly, other cell cycle regulators in the CDK4/Cyclin D-E2F1 pathway, downstream of p16, have also been implicated in the control of hepatic lipid metabolism (7,8). In addition, impaired hepatic lipid metabolism (e.g. as seen during aging or exposure to high fat diet feeding) is associated with increased hepatic expression of p16 (9,10). However, whether p16 directly regulates hepatic lipid homeostasis remains unknown.
The AMP-activated protein kinase (AMPK) is an important regulator of hepatic lipid homeostasis. During prolonged fasting, AMPK senses cellular energetic deficit and activates fatty acid oxidation (FAO) to reestablish normal energy balance (11). AMPK is a heterotrimeric complex composed of one catalytic subunit a and two regulatory subunits b and g, each of which have several different isoforms. There are two catalytic subunit isoforms, AMPKa1 (PRKAA1) and AMPKa2 (PRKAA2), and their phosphorylation at Thr-172 is critical for AMPK activation. The tissue expression of these two isoforms is different, however, whether there are specific roles for each isoform remains unresolved. Interestingly, CDK4 phosphorylates AMPKa2, and not AMPKa1, at several sites (other than Thr-172) thereby suppressing AMPK activity and FAO (12). Moreover, other cell cycle regulators such as Cyclin D1, a CDK4 partner (13), inhibit the activity of PPARa, a master regulator of lipid metabolism gene expression and downstream effector of AMPK activation. Overall, these studies highlight the close interplay between cell cycle and energy balance regulators.
In this study, we assessed the effects of hepatic p16 expression on fasting lipid metabolism, and found that modulation of p16 expression regulates hepatic FAO, mitochondrial function, and ketogenesis. p16-deficiency leads to activation of a cascade involving AMPKa2, SIRT1, and PPARa, which drives enhanced expression of lipid catabolism genes. Interestingly, we found these effects of p16-deficiency to be independent of changes in CDK4 activity. Moreover, adenovirus-mediated overexpression of p16 led to accumulation of LD in vitro, and decreased hepatic mitochondrial activity and ketogenesis in vivo. Our findings highlight a new role for p16 in the hepatic response to fasting and uncover a novel mechanism by which p16 may contribute to the development of metabolic diseases via modulation of hepatic mitochondrial function.

p16-deficiency increases fatty acid catabolic gene expression in hepatocytes
To determine whether p16-deficiency affects hepatic lipid metabolism, we performed transcriptomic profiling of primary hepatocytes isolated from p16 1/1 and p16 2/2 mice cultured under conditions mimicking fasting (as described under "Materials and methods"). There were 3289 differentially expressed genes between p16 1/1 and p16 2/2 cells (false discovery rate , 0.2). Gene Ontology terms enrichment analysis revealed that several metabolic pathways are significantly affected by p16 deficiency, including oxidation reduction, a-amino acid metabolism, and xenobiotic metabolism (Fig. 1A), suggesting that p16 regulates mitochondrial and peroxisomal functions. As expected, principle component analysis of the 3289 differentially expressed genes identified a clear separation based on p16 genotype (Fig. 1B). Interestingly, the loadings for several genes involved in FAO and ketogenesis (e.g. Cpt1a, Ehhadh, Hmgcs2, Bdh1) were highly correlated with PC1 (Fig. 1B), indicating that their expression pattern contributes strongly to the difference between the p16 1/1 and p16 2/2 transcriptomes. We then specifically investigated metabolic pathways associated with fatty acid (FA) mobilization, storage and utilization (Fig. 1C). Interestingly, p16-deficiency was associated with marked changes in these pathways, including activation of both mitochondrial and peroxisomal catabolic pathways (e.g. Ehhadh, Cpt1a, Hmgcs2, and Bdh1, Fig. S1A).
Because primary hepatocytes from p16-deficient mice could be affected by an altered developmental program due to p16deficiency, we next tested whether acute modulation of p16 expression was sufficient to induce similar changes in gene expression in vitro. In murine AML12 cells (Fig. 1D) and immortalized human hepatocytes (IHH, Fig. 1E), siRNA-mediated reduction of p16 led to a strong increase in several genes identified in the transcriptomic analysis, including HMGCS2, CPT1A, and EHHADH. Furthermore, adenovirus-mediated p16-overexpression in HepG2 cells resulted in reduced expression of these same genes (Fig. 1F). Taken together, these results indicate that p16-deficiency in hepatocytes leads to the induction of a transcriptional program indicative of increased fatty acid utilization.
p16 controls mitochondrial fatty acid oxidation and lipid droplet formation in hepatocytes We next sought to investigate whether activation of lipid metabolism gene expression as a result of p16-deficiency is reflected by corresponding functional changes in hepatocyte lipid metabolism and mitochondrial function. Measurement of b-oxidation with 14 C-labeled oleic acid in p16 1/1 and p16 2/2 primary hepatocytes revealed that formation of 14 CO 2 (complete oxidation, Fig. 2A) was increased in p16 2/2 compared with p16 1/1 cells, and associated with a tendency toward higher 14 C-acid soluble metabolites, mostly ketone bodies in hepatocytes (incomplete oxidation, Fig. 2B). To test whether hepatic mitochondrial activity and ketone body production are also modulated by p16 in vivo, we first assessed fasting blood b-hydroxybutyrate (bOHB) levels, which account for ;75% of circulating ketone bodies in mice (14). Unexpectedly, bOHB concentration did not differ between p16 1/1 and p16 2/2 mice after overnight fasting (Fig. 2C). However, acute injection of overnight fasted mice with octanoate, a medium chain fatty acid shown to specifically activate hepatic ketogenesis (15), resulted in a statistically significant augmentation in blood bOHB levels in p16 2/2 mice within 60 min. Moreover, the area under the curve of bOHB was also higher in p16 2/2 mice compared with p16 1/1 littermate controls (Fig. 2D). Finally, we tested whether liver-specific p16 overexpression by adenovirus is sufficient to decrease hepatic mitochondrial function in C57BL6/J mice. Hepatic p16 overexpression (Fig. S2A) resulted in a significant decrease in fasting blood bOHB levels ( Fig. 2E) and hepatic mitochondrial respiratory control ratio (RCR) measured in the presence of palmitoylcarnitine (Fig. 2F) when compared with GFP-overexpressing controls. Importantly, these changes were independent of differences in body weight (Fig. S2B) and plasma free fatty acids (Fig. S2C), the main substrate for ketone body synthesis.
To test which mitochondrial processes were most affected by p16 deficiency, oxygen consumption measurements using the Seahorse Extracellular Flux Analyzer were performed. Interestingly, p16 siRNA-transfected IHH cells displayed increased maximal respiration capacity compared with control siRNAtransfected cells ( Fig. 2G and Fig. S3A). Finally, we assessed the effect of p16-knockdown on intracellular triglyceride and lipid droplet (LD) accumulation by BODIPY staining. P16-knockdown resulted in a significant decrease in LD number per cell (Fig. 2, H and I), which was associated with decreased expression of perilipin 2 (PLIN2, Fig. 2J), an LD-associated protein whose expression correlates with intracellular lipid storage. In line with the decreased mitochondrial activity observed in vivo (Fig. 2, E and F), p16 overexpression in HepG2 cells resulted in a strong increase (.20-fold) in LD per cell compared with control infected cells (Fig. 2, K and L), and associated with increased PLIN2 expression (Fig. 2M). These results indicate that p16 overexpression leads to accumulation of intracellular LD accompanied by reduced mitochondrial FAO.

p16 expression modulates HMGCS2 expression independently of CDK4
Thus far, we have demonstrated that modulation of p16 expression leads to changes in fatty acid metabolism in both nonproliferative primary hepatocytes and proliferating hepatocyte cell lines, indicating that p16's action on this pathway may be independent of its role in cell cycle regulation. The canonical function of p16 is to block the G0 to S phase transition via inhibition of CDK4-mediated phosphorylation of Rb protein (4). Consequently, we analyzed whether the p16-dependent changes in mitochondrial activity were associated with altered cell cycle progression. To do so, we first analyzed changes in CDK4 protein levels and Rb phosphorylation in p16-silenced IHH cells (Fig. 3A). As expected (16), both CDK4 expression and Rb phosphorylation were reduced within 1-4 h of shifting cells to conditions mimicking fasting (low glucose, no insulin, supplemented with forskolin). Interestingly, p16 silencing did p16INK4a controls hepatic fatty acid oxidation Figure 1. p16-deficiency increases fatty acid catabolism gene expression in hepatocytes. A-C, primary hepatocytes were isolated from p16 1/1 and p16 2/2 mice and incubated for 8 h in 1 mM glucose DMEM supplemented with 10 mM FSK. A, Gene Ontology terms enrichment for differentially expressed genes between p16 1/1 and p16 2/2 primary hepatocytes. B, principle component analysis of differentially expressed genes between p16 1/1 and p16 2/2 primary hepatocytes with selected FAO and ketogenesis gene loadings plotted. C, a heatmap of genes in the major fatty acid metabolic pathways in p16 1/1 and p16 2/2 primary hepatocytes. D-F, mRNA relative expression of P16INK4A, HMGCS2, CPT1A, and EHHADH were measured by qPCR in (D) the murine hepatocyte cell line AML12 transfected with siRNA-CDKN2A (siCDKN2A), (E) the human hepatocyte cell line IHH transfected with siRNA-p16 (sip16), and (F) the human hepatocyte cell line HepG2 infected with p16-adenovirus (ad-p16). All values are expressed as mean 6 S.D. *, compared with siCTR or ad-CTR, unpaired t test: ****, p , 0.0001; ***, p , 0.001; **, p , 0.01; and *, p , 0.05. See also Fig. S1. , and B, incomplete fatty acid oxidation ( 14 C-ASM) measured in primary hepatocytes isolated from p16 1/1 and p16 2/2 mice and incubated with 14 C-labeled oleic acid. C, ketone body (b-OHB) production in p16 1/1 (n = 9) and p16 2/2 (n = 9) 12-week-old male mice fasted overnight and injected with a solution of sodium octanoate (250 mM). D, area under the curve (AUC) of sodium octanoate-stimulated b-OHB production. E and F, 12-week-old male mice were injected with adenovirus GFP (ad-GFP) (n = 10) or adenovirus p16 (ad-p16) (n = 11). E, circulating b-OHB was measured following an overnight fasting. F, RCR was measured on liver homogenates from ad-GFP (n = 5) and ad-p16 (n = 5) injected mice. G, OCR in IHH cells transfected with either siCTR or sip16 for 48 h and incubated for 1 h in Seahorse assay media. Oligomycin, ATP synthase inhibitor. p16INK4a controls hepatic fatty acid oxidation not affect total CDK4 nor Rb phosphorylation under these conditions (Fig. 3A). Moreover, p16 silencing in these cells did not modify cell cycle progression as measured by propidium iodide staining (Fig. 3B). Finally, we assessed whether CDK4 activation was also required for the regulation of HMGCS2 gene expression upon p16 deletion. Pharmacological CDK4 inhibition by PD0332991 (Palbociclib, PD033) ( Fig. 3C) or CDK4 silencing by siRNA ( Fig. 3D and Fig. S4A) did not affect the siRNA-p16mediated increase in HMGCS2 expression.
Because p16 overexpression is associated with a senescent phenotype and senescence increases LD accumulation in hepatocytes (9), we assessed several parameters of senescence in ad-p16 overexpressing HepG2 cells (Fig. 2, K-M). A kinetic study of adenovirus-mediated p16 overexpression in HepG2 cells showed no differences in p21 mRNA expression, a marker of senescence, between ad-p16 and ad-Control cells at 24 and 48 h post-infection. However, an increase of p21 transcript levels was observed in ad-p16-infected HepG2 cells after 72 h (Fig. S4, B and C). Interestingly, the modest increase in p21 mRNA observed at the 72-h time point was not associated with modulation of cell cycle progression nor with an increase of SA-b-galactosidase activity, as measured by C12FDG staining, suggesting that p16 overexpression may not fully induce senescence in this time frame (Fig. S4, D and E). Conversely, BODIPY staining showed that ad-p16-infected HepG2 cells present significantly more lipid droplets per cell as early as 24 h post-infection and this is maintained up to 72 h (Fig. S4, F and G). This effect is associated with a gradual and significant increase of PLIN2 mRNA expression at 48 and 72 h and a significant decrease of HMGCS2 mRNA expression at 72 h post-infection (Fig. S4, H and I). In addition, adenovirus-mediated p16 overexpression in the liver was associated with decreased p21 mRNA expression and no changes in hepatic b-galactosidase activity (Fig. S4, J and K). Taken together, these data suggest that p16 overexpression in hepatocytes may impact lipid metabolism prior to a possible induction of cellular senescence. These results indicate that the effects of p16 on hepatocyte mitochondrial activity are unlikely to be mediated by modification of cell cycle progression and occur independently of CDK4 activity.

p16 protein expression is modulated by nutritional changes in IHH cells
We then wondered if suppression of p16 expression was a physiological mechanism occurring during fasting. p16 mRNA expression was assessed in the liver of mice fasted overnight and refed for 6 or 24 h (Fig. 4A) or in mouse primary hepatocytes (Fig. 4B) and IHH cells (Fig. 4C) incubating in media mimicking fed or fasted conditions. None of these nutritional changes altered the expression of p16 at the mRNA level. Interestingly, p16 protein levels were increased in IHH cells in the refeeding experiment (Fig. 4, D and E), suggesting that p16 is regulated post-translationally. A recent report demonstrated p16INK4a controls hepatic fatty acid oxidation that p16 is degraded by autophagy (17). Accordingly, rapamycin treatment for 8 h in IHH cells led to a strong decrease of p16 protein (Fig. 4F), suggesting that autophagy-mediated p16 degradation likely links nutritional status to p16 activity in hepatocytes.
p16 silencing activates the AMPKa2-SIRT1 pathway To further dissect the mechanism of gene regulation by p16, we focused on HMGCS2 expression. We first assessed whether p16 may inhibit AMPK, an important regulator of hepatic FAO and mitochondrial function, to drive changes in gene expression. Interestingly, p16-silencing in IHH cells resulted in increased AMPK phosphorylation at Thr-172 (Fig. 5A), a marker of AMPK activation. Moreover, siRNA knockdown of the AMPKa2 catalytic subunit (siPRKAA2) (Fig. 5, B and C, and Fig. S5A) abrogated the p16-dependent increase in HMGCS2 expression, whereas knockdown of the AMPKa1 (siPRKAA1) catalytic subunit had no effect (Fig. S5, B-D). Interestingly, we found that knockdown of AMPKa1 increased HMGCS2 expression and decreased PLIN2 expression in IHH cells (Fig. S5H), whereas knockdown of the AMPKa2 increased PLIN2 expression (Fig. S5I), mimicking the effects observed after knockdown or overexpression of p16, respectively. p16 silencing in IHH cells also increased mRNA and protein expression/nuclear localization of AMPKa2, but not AMPKa1 (Fig. 5, B-D, and Fig.  S5, C-G), suggesting that p16 silencing may also increase AMPK activity specifically through activation of the AMPKa2 catalytic subunit. In line, reduced total AMPK and AMPK phosphorylation at Thr-172 were observed in livers of over-night fasted mice overexpressing p16 (Fig. 5E). Finally, p16knockdown in IHH cells was associated with increased mitochondrial ROS, an activating signal for AMPK (18), but not total intracellular ROS (Fig. 5, F and G).
Because AMPK activation increases the NAD 1 /NADH ratio leading to SIRT1 activation (19), the role of SIRT1 in the p16mediated regulation of HMGCS2 gene expression was also assessed. Treatment of IHH cells with EX527 ( Fig. 5H and Fig.  S5J), a pharmacological SIRT1 inhibitor, or siRNA-mediated silencing of SIRT1 (Fig. 5, I and J, Fig. S5K) abrogated the increase of HMGCS2 gene expression upon p16-silencing. Altogether, these data suggest that p16-silencing increases HMGCS2 expression and mitochondrial activity via an AMPKa2 and SIRT1-dependent signaling pathway.

p16 regulates fatty acid catabolism via activation of PPARa in hepatocytes
As PPARa is a master regulator of fasting-induced ketogenesis and a key effector of AMPK activation (20), we next assessed the effects of p16-silencing on PPARa activation in vitro. PPARa activation using the specific agonist GW7647 increased HMGCS2 mRNA expression in primary murine hepatocytes, AML12, IHH, and HepG2 cells and p16-silencing potentiated this induction (Fig. 6, A-C, Fig. S6, A-D) in all models. Conversely, PPARa knockdown partially blocked the induction of HMGCS2 mRNA upon p16-silencing in AML12 and IHH cells (Fig. 6, D and E, Fig. S6, E-H). To assess whether p16 may regulate FAO and mitochondrial function in human liver, expression of CDKN2A, PPARA, and HMGCS2 was measured in liver

Discussion
In the present study, we describe a novel function of the cell cycle inhibitor p16 as a regulator of hepatic fatty acid metabolism. p16 deletion in both murine and human hepatocytes enhanced fatty acid oxidation and ketogenesis by a mechanism involving activation of AMPKa2, SIRT1, and consequently, PPARa. Moreover, our data indicate that regulation of these target genes by p16 is independent of its action on CDK4 and independent of alterations in cell cycle progression.
Cell cycle-independent regulation of metabolism is increasingly becoming an established function of many cell cycle regulators. The CDK4-E2F1 pathway was shown to increase hepatic lipogenesis and decrease FAO in nonproliferating cells under feeding conditions (7,8), thus contributing to steatosis and NAFLD development. In contrast, Cyclin D1, a CDK4 partner, was demonstrated to inhibit hepatic lipogenesis through repression of ChREBP, HNF4 (22), and PPARg (23). More recently, Cyclin D1 was found to inhibit PPARa transcriptional activity in a CDK4-independent manner (13), highlighting the complexity of metabolic regulation by cell cycle regulators. Our data indicate that p16 inhibits fasting-induced FAO through PPARa modulation independently of CDK4 inhibition, suggesting that p16 and CDK4 may exert opposing actions on FA metabolism, perhaps depending on the nutritional context. Further studies are necessary to better understand the interplay of p16 and CDK4 on hepatic metabolic functions.
Our results show that p16 acts via AMPK and SIRT1 to modulate PPARa transcriptional activity in hepatocytes. AMPK and SIRT1 are important nutrient sensors regulating FAO during fasting. SIRT1 deacetylates PGC1a, thereby increasing its association with PPARa and subsequent transcriptional activity (24). Previous studies suggest that the AMPKa2 catalytic isoform controls the balance between lipogenesis and FAO in the liver (25). More recently, AMPKa2 was also shown to be required for efficient FAO in mouse embryonic fibroblasts (12). Several studies have shown that the AMPKa2 subunit could be nuclear localized (26,27), and previous work from our group has shown that AMPK is recruited to chromatin and participates in transcriptional regulation with glucocorticoid receptor and PPARa (28). In the present study, we show that p16silencing is associated with increased AMPK Thr-172 phosphorylation and a specific increase in the mRNA and protein expression of the AMPKa2 (but not a1). Furthermore, only knockdown of AMPKa2 was sufficient to reverse the effects of p16-silencing on HMGCS2 gene expression. In line, adenovirusmediated p16 overexpression in the liver in vivo decreases FAO and AMPK expression and phosphorylation. Together, these results suggest that p16 could specifically activate the AMPKa2 isoform to modulate PPARa transcriptional activity.
The mechanism by which AMPKa2 activity is regulated in hepatocytes is not fully understood. AMPK is an energetic sensor that is also highly sensitive to the cellular redox state. Indeed, AMPK can be activated by ROS-induced ATP depletion (18) and various mechanisms of direct ROS interaction (29,30). Increased mtROS production, as observed upon p16knockdown, could thus explain increased AMPK activity (Fig.  5G). However, further studies are required to better characterize the signaling cascade leading to a specific activation of AMPKa2 in p16-deficient hepatocytes.
We observed that the p16 protein level is increased in IHH cells during refeeding conditions, without modulation of p16 mRNA expression by fasting. In line, the p16 protein was recently demonstrated to be degraded by autophagy (17), a process induced by fasting and inhibited by refeeding (31,32). Moreover, partial deletion of ATG5 in mice, resulting in  Table S1. inhibition of autophagy, increased p16 protein levels in the liver (33). We propose that p16 protein degradation during fasting, likely via an autophagy-dependent mechanism, could participate in the physiological activation of FAO by PPARa in hepatocytes.
In line with our previous study showing that p16-deficiency enhances gluconeogenesis through the CDK4-PKA-CREB pathway (6), the increased PPARa activity upon p16-deficiency could also contribute to this increase via regulation of the expression of several genes involved in the conversion of glycerol to glucose (34) or modulation of substrate utilization. Indeed, both pyruvate (via expression of Pdk4) (35) and glycerol (36) metabolism are sensitive to PPARa modulation.
Finally, we observed that whereas overexpression of p16 in HepG2 cells increases lipid accumulation, p16 impacts neither the cell cycle nor senescence in hepatocytes in the conditions studied. In line, adenovirus-mediated overexpression of p16 in the mice livers led to a decrease of circulating ketone bodies and mitochondrial activity in the presence of palmitoylcarnitine without induction of senescence. Aging and senescence are associated with decreased fatty acid oxidation resulting in increased hepatic fat accumulation (9). Although p16 expression is a widely used marker of senescence, recent evidence suggests that a senescence-associated increase in liver p16 expression may not be related to its expression in hepatocytes (37,38). It is also possible that increased p16 may impact these metabolic pathways prior to the onset of outright senescence in a chronic setting.
Overall, our results provide further evidence that p16 plays an important role in the regulation of hepatic glucose and lipid metabolism and again support the genetic link between p16/ CDKN2A and metabolic disorders such as T2D and cardiovascular disease. Lowering p16 expression may help to restore hepatic FAO and lipid homeostasis by activating the AMPKa2-SIRT1-PPARa signaling pathway.

Materials and methods
Genetic terminology: CDKN2A versus p16 The use of CDKN2A in the manuscript refers to the entire CDKN2A locus (P16 and P14ARF in human or P19arf in mouse), when analysis of P16 mRNA expression alone or its modulation was not possible (human liver transcriptomic data or siRNA-mediated knockdown of p16 in the murine AML12 cell line). Conversely, when using the terminology P16, this means that analysis of P16 mRNA expression alone or its modulation was technically possible (siRNA in human cell lines, adenovirus-mediated overexpression of P16 or measure of P16 mRNA expression by qPCR).

Mouse model and diets
p16 2/2 and littermate control (p16 1/1 ) mice on a C57Bl6/J background (provided by P. Krimpenfort and backcrossed more than 8 generations) were housed under SPF conditions (22 6 2°C) in conventional cages with free access to water and food unless indicated otherwise. Experimental procedures were conducted with the approval of the ethics committee for animal experimentation of the Nord Pas-de-Calais region (APAFIS number 12317-2015121612289958 V15).

Ketogenesis assay in vivo
Following an overnight fast, 12-week-old p16 1/1 (n = 9) and p16 2/2 (n = 9) male mice were injected with 250 mM sodium octanoate and b-hydroxybutyrate was measured with a ketometer and b-ketone test strips (Abbott Diabetes Care, Freestyle Optium Neo) before the injection and every 30 min following the injection.

Adenoviral-mediated overexpression of p16
16-Week-old male mice were injected via the tail vein with a solution containing 1 3 10 8 genome copies of ad-GFP (n = 10) or ad-p16 (n = 11). 4 days after adenovirus injections, mice were fasted overnight, plasma b-hydroxybutyrate was measured as described above, and livers were harvested for mitochondrial function analysis.

Primary hepatocyte isolation, culture, and treatments
Primary hepatocytes were isolated from 12-week-old male mice using the two-step collagenase perfusion method, essentially as described previously (6). Following isolation, cells were seeded on collagen-coated plates in Williams E media (Life Technologies) supplemented with 0.1% BSA, 2 mM glutamine, 100 nM bovine insulin, 100 nM dexamethasone, and penicillin and streptomycin. The following morning, cells were washed with PBS and incubated in DMEM, 1 mM glucose, supplemented with 0.1% BSA and antibiotics for 8 h, described as fasting conditions.

Microarray analysis
Mouse primary hepatocytes transcriptomics-Analysis was performed using Mouse Gene 2.0ST arrays (Affymetrix). Array data processing was performed using Bioconductor in the Renvironment (r-project.org). Gene expression changes were calculated after signal normalization using robust multichip averaging in the oligo package (39). Differential gene expression was assessed using the limma package (40) with a threshold of 5% false discovery rate. Gene Ontology (GO) terms enrichment of selected clusters was performed using the clusterProfiler package (41,42).
Human liver transcriptomics-Liver biopsies were obtained from morbidly obese patients during abdominal surgery. Patients underwent a full clinical workup to assess metabolic characteristics and surgical risks. All experimental procedures abide by the Declaration of Helsinki Principles and written informed consent was obtained from each participant as part of the ABOS protocol (ID: NCT01129297). The protocol was also reviewed and approved by the Lille University Hospital ethical committee. RNA quality was assessed using the Agilent Bioanalyzer and only samples with RIN . 5 were included for analysis. Transcriptomic analysis of human liver biopsies was performed using Affymetrix HTA 2.0 arrays (21). Array pre-processing and gene expression analysis were performed as described for p16INK4a controls hepatic fatty acid oxidation murine arrays. Basic clinical characteristics of the patient population are provided in Table S1.
Cell culture, siRNA, and adenovirus experiments Alpha mouse liver 12 (AML12), IHH, and the human hepatoma cell line HepG2 were maintained as described under in details, under a humidified incubator at 37°C under 5% CO 2.
Human hepatoma cell line HepG2-HepG2 cells (catalog number HB-8065; ATCC) were seeded on 24-well-plates and starved during 24 h for cell cycle synchronization. Starvation medium was replaced by medium supplemented with 2% FCS containing p16-adenovirus (kind gift of Dr. A. Mazo, University of Barcelona, Spain) or empty adenovirus with a multiplicity of infection of 60. 24 h after infection, medium was replaced by fresh medium with 10% FCS. 72 h after infection, cells were collected for BODIPY staining or mRNA analysis.

Reverse transcription real time quantitative PCR analysis
Total RNA was isolated from cells using EXTRACT ALL ® (Eurobio). cDNA was generated using the high capacity reverse transcription kit (Applied Biosystems, Life Technologies, USA). Real-time qPCR was performed on a Stratagene Mx3005P system (Agilent Technologies, Santa Clara, CA) using specific primers (see supporting "Methods"). mRNA levels were normalized to cyclophilin A mRNA and fold-induction was calculated using the DDC t method.

Protein extraction and Western blotting analysis
Cells were scraped in protein lysis buffer (50 mM Tris-HCl, pH 8, 137 mM NaCl, 5 mM EDTA, 2 mM EGTA, 1% Triton X-100) on ice, transferred to 1.5-ml Eppendorf tubes and rotated for 20 min at 4°C, followed by sonication for 5 min and centrifugation at 13,000 3 g for 10 min at 4°C. The resulting supernatants were stored in aliquots at 220°C until analysis. Cell lysate protein concentrations were determined using a BCA protein assay kit (Pierce). The protein concentration of samples was equalized and cell lysates were mixed with 43 LDS nonreducing loading buffer (Life Technologies). Samples were heated at 95°C for 5 min before loading and separated on precast 4-12% BisTris polyacrylamide gels (Life Technologies). Proteins were transferred to a nitrocellulose membrane using the iBlot2 Dry Blotting System (Life Technologies). Membranes were then incubated overnight at 4°C with various primary antibodies in blocking buffer containing 5% nonfat milk at the dilution specified by the manufacturers. The following antibodies were used: p16 (BD Biosciences, number 554079), CDK4 (Santa Cruz Biotechnology, number sc-260), p-Rb (Abcam, ab10921), AMPK (Cell Signaling, number 2532), and p-AMPK T172 (Cell Signaling, number 2531).

Cell cycle analysis
Cells were washed with PBS, harvested, and fixed with cold 100% ethanol for 10 min at 4°C. Cells were again washed with PBS, centrifuged, and incubated with a solution of 50 mg/ml of propidium iodide (PI) (Sigma-Aldrich) with 200 mg/ml of RNase A (Qiagen) and NucBlue ® (Invitrogen) for 30 min at 37°C. Cells were then washed, resuspended in PBS, and analyzed by flow cytometry with BD LSR Fortessa X-20 cell analyzer (BD Biosciences). PI fluorescence was analyzed with blue laser 488-and 575-nm emission filters. NucBlue ® fluorescence was analyzed with UV laser 355-and 496-nm emission filters. Acquisition was performed with FACS DIVA software (BD Biosciences). PI fluorescence was analyzed on live singlecell population. The percentage of cells in each cell cycle phase was analyzed by using the FlowJo TM software algorithm (BD Biosciences).

Immunofluorescence analysis
IHH cells seeded on coverslips were washed with cold PBS and fixed in 4% paraformaldehyde for 10 min at 4°C. Cells were permeabilized with 0.2% Triton X-100 (Euromedex) in PBS for 20 min at room temperature and blocked for 30 min using blocking solution (BSA 1%, 0.2% Triton X-100 in PBS). Cells were then incubated overnight at 4°C with AMPKa1 (Abcam, ab32047) and AMPKa2 (Abcam, ab3760) antibodies in blocking solution before being washed and incubated 1 h at room temperature with secondary antibodies coupled to Alexa 594 (Molecular Probes by Life Technologies) and incubated for 10 min with Hoechst (Invitrogen). Cells were mounted in fluorescence mounting medium (GBl Laboratory). Acquisition was performed using an inverted confocal microscope (LSM 710, Zeiss) with a 340 oil-immersion lens (NA 1.3 with an optical resolution of 176 nm). Alexa 594 and NucBlue were imaged p16INK4a controls hepatic fatty acid oxidation using UV, argon 488 nm and 561 nm lasers. Images were processed with Zen software and analyzed with Icy software.

BODIPY 493/503 staining
Cells seeded on coverslips were fixed in 4% paraformaldehyde for 10 min at 4°C and incubated with BODIPY 493/503 (Molecular Probes) at 200 ng/ml and NucBlue ® for 20 min at room temperature. Cells were mounted in fluorescent mounting medium (GBl Laboratory). Acquisitions were performed using an inverted confocal microscope (LSM 710, Zeiss) with an airyscan detector and 340 oil-immersion lens (NA 1.3 with an optical resolution of 176 nm) and 1.83 digital zoom. BOD-IPY and NucBlue fluorescence were imaged using argon 488 nm and UV lasers, respectively. Images were processed with Zen and analyzed with Icy programs respectively.

Mitochondrial activity measurement
Fatty acid oxidation test-Fatty acid oxidation was measured essentially as described previously (44). Freshly isolated primary hepatocytes were plated in DMEM-high glucose (25 mM) 1 10% FBS for 4-6 h and changed to DMEM-low glucose (5 mM) overnight. 24 h after isolation, cells were incubated with 500 mM [ 14 C]palmitic acid:BSA complex. After 3-4 h, media was transferred to a CO 2 trap and 14 CO 2 was captured in 1 N NaOH in adjacent wells. Captured 14 CO 2 (complete oxidation) and 14 C-acid soluble metabolites (ASMs, incomplete oxidation) were subjected to scintillation counting.
Seahorse analysis-Measurements of oxygen consumption rate (OCR) were performed using the XF24 apparatus (Seahorse Bioscience, North Billerica, MA). Briefly, cells were plated into XF24 (V7) polystyrene cell culture plates (Seahorse Bioscience, North Billerica). IHH cells were seeded at 70,000/ well (XF24 plate). The cells were incubated for 24 h and transfected with siRNA for 48 h. Sensor cartridges were calibrated prior to each assay. IHH cells were incubated for 1 h in complete Seahorse assay medium (10 mM glucose, 2 mM glutamine) supplemented or not with 0.5 mM Etomoxir (Eto) in a 37°C/ non-CO 2 incubator for 60 min prior to the start of an assay. All experiments were performed at 37°C. Each measurement cycle consisted of a mixing time of 3 min, a waiting time of 2 min, and a data acquisition period of 3 min. OCR data points represent to the average rates during the measurement cycles. Oligomycin, FCCP, and Rotenone/Antimycin A were prepared at, respectively, 2 and 1 and 0.5 mM concentrations in the desired assay medium and adjusted to pH 7.4. A volume of 50 to 60 ml of compound was added to each injection port. In a typical experiment, 3 baseline measurements were taken prior to the addition of any compound, and 3 response measurements were taken after the addition of each compound. OCR were normalized to protein concentration (pmol/min/mg of proteins).

ROS measurement
Total ROS and mitochondrial superoxide generation was determined using the CellROX Deep Red (Invitrogen) and MitoSOX red (Invitrogen), respectively. Briefly, IHH cells were seeded on a 12-well-plate and transfected with siRNA-CTR or siRNA-p16. 48 h after transfectection, IHH cells were harvested and resuspended in culture media containing 5 mM CellROX or MitoSOX and incubating for 20 min at 37°C in the dark. Cells were then washed, resuspended in PBS, and analyzed by flow cytometry with BD LSR Fortessa X-20 cell analyzer (BD Biosciences). CellROX fluorescence was analyzed with extinction and emission at 640 and 665 nm. MitoSOX fluorescence was analyzed with extinction and emission at 510 and 580 nm.
HTRF ® P-AMPK and total AMPK assays A piece of liver (;50 mg) was mechanically disrupted with a Polytron in ;500 ml of lysis buffer. Homogenates were centrifuged and supernatant was collected for protein quantification. Samples containing 0.15 or 0.5 mg/ml of proteins were, respectively, added into 96-well HTRF plates for P-AMPK and total AMPK analysis, according to manufacturer's instructions (Cisbio).

Statistics
All data were expressed as mean 6 S.D. or S.E. of measurement, as indicated in the figure legends. Student's unpaired t test was used for assessing statistical differences between two groups, whereas comparison more than two groups was performed using one-way ANOVA. Differences between two groups under two different factors were examined by a twoway ANOVA. In all analyses, a p values of p , 0.05 was considered significant.

Data availability
Raw datasets for mouse primary hepatocytes transcriptomics have been deposited in the GEO database through GEO Series accession number GSE134625.