Theranostics 2022; 12(4):1589-1606. doi:10.7150/thno.69054 This issue
1. Institute of molecular rhythm and metabolism, Guangzhou University of Chinese Medicine, Guangzhou 510405, China.
2. College of Pharmacy, Jinan University, Guangzhou 510632, China.
3. The First Affiliated Hospital of Jinan University, Guangzhou 510632, China.
# These authors contribute equally to this work.
While growing evidence suggests that circadian clock and obesity are intertwined, the underlying mechanism is poorly understood. Here, we investigate how circadian clock is linked to obesity.
Methods: Metabolomics profiling of WAT (white adipose tissue) samples was performed to identify the metabolites altered in obese model. mRNA levels were analyzed by qPCR assays. Proteins were detected by immunoblotting, immunofluorescence and ELISA. ChIP and luciferase reporter assays were used to investigate epigenetic and transcriptional regulation.
Results: Obesity causes perturbance of circadian clock in WAT in mice and humans, particularly, BMAL1 is markedly reduced. Metabolomic analysis reveals reduced glutamine and methionine in obese WAT. Glutamine metabolism contributes to production of acetyl-CoA, whereas methionine metabolism generates S-adenosyl methionine (SAM). Acetyl-CoA and SAM are the substrates for histone acetylation and methylation, respectively. Reduced glutamine and methionine in obese WAT are associated with decreased H3K27ac and H3K4me3 at Bmal1 promoter. Consistently, glutamine or methionine administration in vitro and in vivo increases H3K27ac or H3K4me3, promoting Bmal1 transcription and expression. A screen of transport and metabolic genes identifies downregulation of the uptake transporter SLC1A5 as a cause of reduced glutamine or methionine in obese WAT. Moreover, we observe impaired expression of PPAR-γ in obese WAT. PPAR-γ trans-activates Slc1a5 via direct binding to a response element in promoter.
Conclusion: Impaired PPAR-γ in obesity provokes downregulation of SLC1A5 and reductions in adipocyte uptake of glutamine and methionine (two epigenetic modulators), leading to disruption of Bmal1. Therefore, PPAR-γ integrates obesity and adipocyte clock, promoting a vicious cycle between circadian disruption and obesity development.
Keywords: PPAR-γ, obesity, circadian clock, BMAL1, SLC1A5
Obesity (defined as abnormal fat accumulation in the body) has become a major public health problem, affecting a growing number of countries. Currently, more than 650 million adults (about 13% of the population) in the world are obese, which has nearly tripled since 1975 . Of note, 33% of Americans are obese, and the prevalence may hit 50% by 2030 . Obesity is closely associated with metabolic complications (e.g., insulin resistance, glucose intolerance, hyperlipidemia, hypertension and fatty livers) and is a risk factor for several deadly diseases, including heart disease, stroke, diabetes and cancers [1,3]. Contributing factors, both genetic and environmental, have been proposed for obesity development, and a combination of excess calorie intake and sedentary lifestyle is regarded as the primary cause . However, calorie restriction coupled with exercise has been largely ineffective in maintaining long-term weight loss, and the obesity pandemic continues unabated . Therefore, there is a desperate need for novel insights into the mechanisms controlling obesity and energy homeostasis.
Circadian clocks orchestrate internal daily (~24 h) rhythms of molecular and cellular processes, leading to circadian variations in physiology and behaviors (e.g., body temperature, heartbeats, feeding behavior and sleep) . In mammals, the clock located in the suprachiasmatic nucleus (SCN) of the brain acts as the central pacemaker coordinating the clocks present in other regions of the brain and in peripheral tissues . These clocks are cell-autonomous oscillators built on several interconnected feedback loops . In the core loop, the transcription factors brain and muscle ARNT-like protein-1 (BMAL1) and circadian locomotor output cycles kaput (CLOCK) associate to form a heterodimer that activates the transcription of clock-controlled genes (CCGs), including periods (PERs) and cryptochromes (CRYs) . As PER and CRY proteins accumulate, they translocate into the nucleus and inhibit the transcriptional activity of BMAL1/CLOCK complex, thereby down-regulating the expression of CCGs in a negative feedback manner . This feedback mechanism generates ~24 h oscillations in circadian gene expression . Extra feedback loops involve the key players REV-REBs, RAR-related orphan receptors (ROR)s, and D-box acting proteins [e.g., albumin D-site-binding protein (DBP) and E4 promoter-binding protein 4 (E4BP4)] that serve to maintain robustness of the clock system by regulating BMAL1 and PER2 expression . In addition, circadian oscillations can be modified by post-transcriptional mechanisms such as phosphorylation of clock proteins by casein kinases .
Disruptions of circadian rhythms, by either genetic or environmental factors, enhance the risk of obesity and metabolic syndrome [12,13-14]. For instance, Clock and Bmal1 mutant mice are obese prone and develop metabolic syndrome [12,15,16]. Social jetlag and shift work (measures of circadian desynchrony in humans) are associated with metabolic disturbances and obesity [14,17]. The mechanisms for circadian dysfunction-induced obesity are complex and multifaceted, including promotion of insulin resistance, inducement of leptin resistance, altered adipogenesis and attenuated feeding rhythm (or increased food intake in rest period) [18-21]. On the other hand, people who are overweight or obese show altered expression of clock genes [22,23]. Likewise, disruptions of clock genes and circadian behaviors are noted in obese mice [24-26]. For instance, Bmal1, Clock and Per2 are altered in fat and liver tissues of high-fat fed mice . These suggest that obesity and circadian system are tightly intertwined (a crosstalk phenomenon), which perhaps is involved in the pathogenesis of obesity and is a contributing factor to the obesity epidemic. However, the molecular mechanisms underlying the crosstalk are poorly defined. In particular, it remains to be answered how circadian clock is linked to obesity, which forms an essential aspect of the clock-obesity crosstalk.
PPARs (peroxisome proliferator-activated receptors) are a family of ligand-responsive nuclear receptors that contain three members, namely, PPAR-α, PPAR-β (also known as PPAR-δ) and PPAR-γ. PPAR family members can be activated by endogenous ligands such as fatty acids and their metabolites . Upon ligand binding, PPARs form heterodimers with retinoid X receptors and bind to a specific DNA response element (PPRE) to regulate gene transcription and expression . PPAR isoforms show tissue-specific differences in their expression and functions . Notably, PPAR-γ is predominantly present in adipose tissues and regarded as a master regulator of adipogenesis and lipid metabolism . Downregulation of PPAR-γ has been noted in obese mice [27,30]. Mutations in PPAR-γ gene lead to partial lipodystrophy in humans [31,32]. In addition, PPAR-γ has a critical role in regulating insulin sensitivity, inflammation, tumor cell growth, apoptosis and differentiation . These specific and crucial functions render PPAR-γ an attractive drug target for metabolic diseases and cancers [34,35]. In fact, the PPAR-γ agonists rosiglitazone and pioglitazone are clinically used as insulin-sensitizing drugs for management of type 2 diabetes.
Here, we show disrupted expression of clock genes in white adipocyte tissue in obese mice and humans, and altered circadian behaviors (feeding rhythm) in obese mice. Mechanistically, PPAR-γ, downregulated in obesity, provokes changes in a cascade of events in the adipocytes, including decreased expression of the uptake transporter SLC1A5, reductions in glutamine and methionine uptake, lower production of acetyl-CoA and S-adenosylmethionine (SAM), and reduced histone acetylation and methylation, thereby repressing Bmal1 transcription and expression. We thus propose that PPAR-γ integrates obesity and adipocyte clock via epigenetic regulation of Bmal1, promoting a vicious cycle between circadian disruption and obesity development.
Disruption of adipocyte clock in obese mice and humans. (A) Schematic diagram of the experimental design for collecting SCN, liver, and WAT from ob/ob, high-fat diet (HFD) and control mice. (B) mRNA expression of clock genes (Bmal1, Rev-erbα, Rev-erbβ, Per1, Cry1 and Dbp) in WAT derived from ob/ob and control mice at 6 circadian times (CT2, CT6, CT10, CT14, CT18 and CT22). CT, Circadian time. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (C) mRNA expression of clock genes (Bmal1, Rev-erbα, Rev-erbβ, Per1, Cry1 and Dbp) in WAT derived from HFD and control mice at 6 circadian time points. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (D) Protein expression of BMAL1, REV-ERBα and DBP in WAT derived from ob/ob, HFD and control mice at 6 circadian times. The middle lane is the pre-stained protein ladder. (E) Immunostaining of BMAL1 protein (green) in WAT derived from ob/ob and control mice at CT10 and CT22. Scale bar = 50 μm. (F) Quantification of green fluorescence intensity in Fig 1E. The fluorescence intensity relative to control (CT10) is shown and expressed as mean ± SD. (G) Protein expression of BMAL1 in omental WAT derived from obese patients (n = 13) and lean individuals (n = 10). The samples for each group were pooled for Western blotting.
To test the effects of obesity on circadian rhythms, we used two mouse models, namely, ob/ob mice and mice with high-fat diet (HFD, 60% kcal from fat)-induced obesity. The two mouse models of obesity were validated by measuring adipocyte diameter and circulating lipid levels (Figure S1). We examined the expression of clock genes in SCN, WAT (white adipose tissue) and liver from the obese and control mice at six circadian times (CT2, CT6, CT10, CT14, CT18 and CT22) (Figure 1A). Obesity resulted in overall disruption of core clock genes (i.e., Bmal1, Clock, Rev-erbα/Nr1d1, Rev-erbβ/Nr1d2, Per1, Per2, Cry1, Cry2, Dbp and Rorα) in WAT (Figure 1B-C & Figure S2). Notably, most clock genes were markedly downregulated, including Bmal1 and its direct targets Rev-erbα/β, Per1/2, Cry1/2, Dbp and Rorα (Figure 1B-C & Figure S2) [36,37]. This suggests a crucial role of Bmal1 in connecting obesity to adipocyte clock machinery. Consistent with mRNA changes, BMAL1, REV-ERBα and DBP proteins were downregulated in WAT of obese mice and their circadian rhythms were blunted (Figure 1D). Deregulation of circadian WAT BMAL1 was further confirmed by immunofluorescence staining (Figure 1E-F). By contrast, core clock genes such as Bmal1, Clock, Rev-erbα and Dbp were unaffected in SCN, and slightly altered in the liver in obese mice compared to control mice (Figure S3). In addition, downregulation of BMAL1 was observed in the omental WAT of obese patients as compared to lean individuals (Figure 1G). Altogether, these findings indicate that obesity is associated with clock dysfunction in WAT.
Because adipocyte clock has a role in regulating feeding rhythms and energy balance , we tested whether obesity affects circadian rhythms related to energy metabolism. Both ob/ob and HFD-induced obese mice showed a decreased locomotor activity with no significant change in circadian period based on actogram analysis (Figure 2A-B & Figure S4). The decreases in daily activities were attributed to reductions in nighttime activities as the daytime activities were increased (Figure 2C). Thus, circadian rhythm in locomotor activity was attenuated in obese mice with an increased percent of activity in the daytime (Figure 2D). The obese mice also had decreased daily energy expenditure (EE) compared to control mice (Figure 2E), consistent with the changes in locomotor activity. Likewise, the EE rhythm was dampened with an increased percent of EE in the light period (Figure 2F). Consistently, circadian rhythm in the daily food intake was also attenuated in obese mice with an increased percent of food intake in the daytime (Figure 2G). Intriguingly, the obesity-associated circadian phenotypes were very similar to those of mice with Bmal1-specifically deleted in the adipocytes . Altogether, obesity in mice attenuates circadian rhythms in locomotor activity, EE, and feeding behavior, that may be linked to dysfunction of adipocyte clock.
Altered circadian behaviors in obese mice. (A) Representative double-plotted actograms for daily locomotor activities of control (left panel), ob/ob (middle panel) and HFD (right panel) mice. Horizontal black and white bars at the top of each actogram represent (subjective) light and dark phases, respectively. (B) Locomotor activities monitored over 3 weeks (1 week under a 12 h light/dark cycle and 2 weeks under constant darkness) for ob/ob, HFD and control mice. ZT, Zeitgeber time. (C) Locomotor activities monitored over 1 week for ob/ob, HFD and control mice kept under a 12 h light/dark condition. Data are mean ± SD (n = 7). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (D) Percent of total activity for the light period. Data are mean ± SD (n = 7). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (E) Energy expenditure (EE) monitored over 3 days for ob/ob, HFD and control mice kept under a 12 h light/dark condition. Data are mean ± SD (n = 3). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (F) Percent of total EE for the light period. Data are mean ± SD (n = 3). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (G) Percent of total food intake for the light period. Data are mean ± SD (n = 3). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test).
Obesity downregulates glutamine and methionine in mouse WAT. (A) Schematic diagram for metabolomic analysis. (B) Heatmap for oscillating metabolites in WAT of ob/ob and control mice at six circadian time points. Red indicates high expression, and blue indicates low expression of metabolites as shown in the scale bar. (C) Venn diagram showing numbers of oscillating metabolites in WAT of ob/ob and control mice. (D) Radar diagram showing phase distribution of metabolites in WAT of ob/ob and control mice. (E) Pie diagram showing percentages of oscillating metabolites in WAT of ob/ob and control mice. (F) Upset plot showing differential metabolites (representative ones are shown in right panel) occurring at more than two circadian times simultaneously. (G) KEGG analysis of obesity-associated differential metabolites in WAT (Top 5 pathways are shown). (H) Glutamine and methionine levels in WAT of ob/ob and control mice at 6 circadian time points. Data are mean ± SD (n = 5). *p < 0.05 (two-way ANOVA and Bonferroni post hoc test). (I) Glutamine and methionine levels in omental WAT of obese and lean humans. Data are mean ± SD (n = 13 for obese & n=10 for lean). *p < 0.05 (t-test).
Because metabolites (particularly oscillating metabolites such as melatonin, glucose, heme and bilirubin) are potential modulators of circadian clock and obesity is associated with metabolic changes in WAT [38-39,40,41], we investigated a role of metabolites in regulation of adipocyte clock by obesity. To this end, we performed global 24-hr metabolomics profiling of WAT from ob/ob and normal mice (Figure 3A). ob/ob mice showed considerable changes in a large number of metabolites in WAT (Figure 3B & Supplementary excel file). JTK_CYCLE algorithm identified 1793 oscillating metabolites in normal mice, and > 35% of them were no longer rhythmic in obese mice (Figure 3C). Analysis of the oscillating metabolites revealed a significantly altered phase distribution in ob/ob mice (Figure 3D). Differential WAT metabolites (based on metaX software calculation) between ob/ob and normal mice were observed at one or more of six circadian times. A majority of differential metabolites were lipids and organic acids (mainly including amino acids) (Figure 3E). We further examined differential metabolites occurring at more than two circadian times simultaneously, which were considered as the ones altered the most (Figure 3F). These metabolites included glutamine, methionine, glutamate and citrate (Figure 3F). KEGG enrichment analysis of differential metabolites identified several specific pathways deregulated in obese mice, including glycerophospholipid metabolism, glutamine and glutamate metabolism, and aminoacyl-tRNA biosynthesis (Figure 3G). Glycerophospholipid mainly serves as a structural component of biological membranes , and its role in regulating circadian clock may be excluded. Mass spectrometric analyses demonstrated downregulation of glutamine and glutamate (a metabolic product of glutamine) in ob/ob mice and loss of their circadian rhythms, consistent with the KEGG findings (Figure 3H & Supplementary excel file). Furthermore, we examined the amino acids (key players) involved in the pathway of aminoacyl-tRNA biosynthesis, and found that glutamine, glutamate and methionine were downregulated to the largest degree (Figure 3H & Supplementary excel file). Supporting the mouse data, both glutamine and methionine were decreased in omental WAT of obese patients as compared to lean individuals (Figure 3I). Glutamine and methionine have been shown to regulate gene transcription and expression [43-44,45], thus may be potential modulators of circadian clock. Altogether, obesity reduces glutamine and methionine, and disrupts their circadian rhythms in mouse WAT.
Obesity impacts epigenetics of Bmal1 in a glutamine- and methionine-dependent manner. (A) Identification of H3K27ac and H3K4me3 as histone marks of Bmal1 gene in adipocytes according to Cistrome database (accession numbers: GSM2433306 and GSM1650056). All Cistrome data have been carefully curated and processed with a streamlined analysis pipeline and evaluated with comprehensive quality control metrics. (B) Protein expression of H3K27ac and H3K4me3 in WAT of ob/ob (upper panel), HFD (lower panel) and control mice (C) Enrichments of H3K27ac and H3K4me3 proteins at Bmal1 promoter in WAT of ob/ob, HFD and control mice at CT6. Pgc-1α and Wnt-6 were used as positive controls for H3K27ac and H3K4me3, respectively. Data are mean ± SD (n = 3). *p < 0.05 (two-way ANOVA and Bonferroni post hoc test). (D) Levels of acetyl-CoA and S-adenosyl methionine (SAM) in WAT of ob/ob, HFD and control mice at 6 circadian time points. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (E) ChIP assays showing enrichments of H3K27ac and H3K4me3 at Bmal1 promoter in WAT of HFD and control mice administered with glutamine, methionine or vehicle. Data are mean ± SD (n = 3). *p < 0.05 as determined by two-way ANOVA and Bonferroni post hoc test.
Metabolism of glutamine in the citric acid cycle generates citrate and contributes to production of acetyl-CoA, the acetyl carrier for histone acetylation . Methionine is an essential metabolite for biosynthesis of SAM which functions as a substrate for histone methylation . Recent studies have proposed glutamine and methionine as critical regulators of gene transcription via epigenetic modifications . We thus investigated whether deregulation of glutamine and methionine in obese mice impacts epigenetic modifications of clock genes. According to Cistrome database , H3K27ac and H3K4me3 (both are associated with active transcription) were histone marks highly enriched at Bmal1 promoter in adipocytes (Figure 4A). Both ob/ob and HFD-induced obese mice had lower levels of total H3K27ac and H3K4me3 proteins in WAT as compared to control mice (Figure 4B). However, other histone marks such as H3K9ac, H3K9me2, H3K9me3 and H3K27me3 were unaffected (Figure S5). ChIP assays showed significantly reduced enrichments of H3K27ac and H3K4me3 at the Bmal1, Pgc-1α and Wnt-6 promoters in WAT of obese mice (Figure 4C). Furthermore, we observed lower levels of acetyl-CoA and SAM in obese WAT (Figure 4D). This was accompanied by reductions in glutamate and citrate (Supplementary excel file). Moreover, glutamine or methionine administration to obese mice (HFD-induced) partially rescued epigenetic changes in Bmal1, Pgc-1α and Wnt-6 (Figure 4E). Additionally, C646 (a histone acetyltransferase inhibitor ) and MM-102 (an H3K4 histone methyltransferase inhibitor ) dose-dependently attenuated the effects of glutamine and methionine on Bmal1 expression as they reduced the level of H3K27ac or H3K4me3 in 3T3-L1 adipocytes (Figure S6). Taken together, these findings indicate that deregulation of glutamine and methionine in obese mice reduces H3K27ac and H3K4me3 at Bmal1 promoter, thereby potentially downregulating Bmal1 transcription and expression.
We next tested whether glutamine and methionine can regulate expression of Bmal1 and related clock genes. Addition of glutamine or methionine enhanced the mRNA expression of Bmal1 and its target clock gene (i.e., Rev-erbα) in synchronized 3T3-L1 adipocytes, whereas deprivation of glutamine or methionine reduced their expression levels (Figure 5A). Consistently, BMAL1 and REV-ERBα proteins in 3T3-L1 adipocytes were increased by the two amino acids (Figure 5B). Likewise, mRNA and protein levels of both BMAL1 and REV-ERBα were elevated in human visceral adipocytes treated with glutamine or methionine (Figure 5C-D). We performed real-time bioluminescence recordings of a BMAL1-luciferase reporter in U2OS cells treated with glutamine or methionine or vehicle (Figure 5E). Glutamine and methionine increased the BMAL1-luciferase reporter activity, suggestive of their abilities to enhance Bmal1 transcription (Figure 5E). We further examined the in vivo effects of glutamine and methionine on adipocyte clock. We initiated the study with 10-week-old mice that had been maintained on a HFD for 6 weeks (Figure 5F). The mice continued on HFD and were treated with glutamine or methionine in drinking water (Figure 5F). We confirmed that glutamine and methionine treatments did increase the concentrations of glutamine/acetyl-CoA and methionine/SAM in mouse WAT, respectively (Figure 5G-H). Glutamine- and methionine-treated mice showed higher adipocyte expression levels of Bmal1 and Rev-erbα as well as increased H3K27ac and H3K4me3 as compared to vehicle-treated mice (Figure 5I-J). Further, glutamine and methionine significantly decreased the body weight, and reduced the mass of WAT and the size of adipocytes, indicating amelioration of mouse obesity (Figure 5K-M). They also protected mice against developing obesity-associated metabolic complications such as hyperlipidemia and hypercholesterolemia (Figure 5N). To examine whether restoration of clock genes is an effect of glutamine and methionine rather than weight loss, we initiated short-term (two weeks) treatment experiments with chow- and HFD-fed mice. Short-term treatment with methionine or glutamine did not impact animal weight, but enhanced the expression of clock genes in WAT, however, obesity-associated metabolic parameters (such as serum triglycerides and cholesterol) were unaffected (Figure S7). Thus, methionine and glutamine had casual effects on adipocyte clock. As the feeding rhythm and EE were normalized in glutamine- and methionine-treated obese mice (Figure 5O-P), we wondered whether Bmal1 enhancement is linked to obesity protection. To test this, we overexpressed Bmal1 in adipocytes in HFD-fed obese mice using recombinant adeno-associated virus serotype 8 (AAV8) vector encoding Bmal1 driven by an adipocyte-specific aP2 promoter (named “AAV8.aP2.Bmal1”), and examined its effects on obesity and associated comorbidities (Figure 5Q-T). We found that adipocyte-specific overexpression of Bmal1 attenuated obesity and metabolic abnormalities (Figure 5R-T). This is consistent with a prior study by Paschos et al in which Bmal1 malfunction in adipocytes results in obesity . Altogether, these findings suggest that glutamine and methionine can restore the function of adipocyte clock to ameliorate obesity in mice.
Glutamine and methionine enhance BMAL1 expression in WAT and protect mice against developing obesity. (A) Effects of glutamine/methionine addition or deprivation on mRNAs of Bmal1 and Rev-erbα in synchronized 3T3-L1 adipocytes. Data are mean ± SD (n = 3). *p < 0.05 as determined by two-way ANOVA and Bonferroni post hoc test. (B) Protein expression of BMAL1 and REV-ERBα in 3T3-L1 adipocytes treated with glutamine or methionine. (C) mRNA expression of BMAL1 and REV-ERBα in human adipocytes treated with glutamine or methionine. Data are mean ± SD (n = 3). *p < 0.05 (t-test). (D) Protein expression of BMAL1 and REV-ERBα in human adipocytes treated with glutamine or methionine. (E) Bioluminescent recordings of BMAL1-Luc U2OS cells after glutamine or methionine treatment. (F) Experimental protocol for assessment of glutamine and methionine on adipocyte clock and obesity development. (G) Glutamine and acetyl-CoA levels in WAT derived from HFD- and chow-fed mice at CT6. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (H) Methionine and SAM in WAT derived from HFD- and chow-fed mice at CT6. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (I) mRNA expression of Bmal1 and Rev-erbα in WAT derived from HFD- and chow-fed mice at CT6. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (J) Protein expression of H3K27ac and H3K4me3 in WAT derived from HFD- and chow-fed mice at CT6. (K) Body weight of HFD- and chow-fed mice. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (L) WAT weight of HFD- and chow-fed mice at CT6. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (M) H&E staining of WAT derived from HFD- and chow-fed mice. Scale bar = 50 μm. (N) Serum triglycerides and cholesterol in HFD and chow mice at CT6. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (O) EE (left panel) and food intake (right panel) monitored over 3 days for HFD- and chow-fed mice kept under a 12 h light/dark condition. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (P) Food intake (day and night) was monitored over 3 days for HFD- and chow-fed mice kept under a 12 h light/dark condition. (Q) Effects of AAV8.aP2.Bmal1 on Bmal1 mRNA in WAT in obese mice induced by HFD. (R) Effects of AAV8.aP2.Bmal1 on body weight of obese mice induced by HFD. (S) Effects of AAV8.aP2.Bmal1 on serum triglycerides and cholesterol in obese mice induced by HFD. (T) Effects of AAV8.aP2.Bmal1 on food intake monitored over 3 days for obese mice kept under a 12 h light/dark condition. Data are mean ± SD (n = 5). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test).
SLC1A5, a dual glutamine and methionine uptake transporter, is downregulated in mouse WAT. (A) mRNA expression of three dual glutamine and methionine uptake transporters (Slc1a5, Slc38a1 and Slc38a2) in WAT of ob/ob, HFD and control mice at 6 circadian times. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (B) Protein expression of SLC1A5 in WAT of ob/ob (top panel), HFD (bottom panel) and control mice at 6 circadian times. (C and D) mRNA (C) and protein (D) expression of SLC1A5 in the omental fat of obese or lean humans. Data are mean ± SD (n = 13 for obese & n = 10 for lean). *p < 0.05 (t-test). (E) Serum glutamine and methionine levels in ob/ob, HFD and control mice at 6 circadian times. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (F) Uptake of [13C5]-glutamine and [2H3]-methionine in 3T3-L1 adipocytes transfected with siSlc1a5 (for gene knockdown/KD) or control with or without Slc1a5 overexpression (OE). Data are mean ± SD (n = 3). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (G) Slc1a5 and Bmal1 mRNAs in synchronized 3T3-L1 adipocytes transfected with siSlc1a5 or control. *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (H) mRNA expression of methionine-related enzymes (Bhmt, Mtr and Mat2a) in WAT derived from ob/ob, HFD and control mice. Data are mean ± SD (n = 5). (I) mRNA expression of glutamine-related enzymes (Glul, Gls and Gls2) in WAT derived from ob/ob, HFD and control mice. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test.
Amino acid levels in the body are tightly regulated by transport and metabolic processes . SLC1A5, SLC38A1 and SLC38A2 are three critical transporters responsible for cellular uptake of both glutamine and methionine . SLC1A5 expression was considerably downregulated, whereas the expression levels of Slc38a1 and Slc38a2 were unaffected in WAT in obese mice compared to normal mice (Figure 6A-B). Consistently, we observed a lower level of SLC1A5 in omental WAT of obese patients (Figure 6C-D). Downregulation of SLC1A5 may account for reduced glutamine and methionine in WAT of obese mice. Supporting this, circulating glutamine and methionine levels were increased in obese mice (Figure 6E). Furthermore, Slc1a5 knockdown in 3T3-L1 adipocytes resulted in reduced uptake of both [13C5]-glutamine and [2H3]-methionine, however, overexpression of Slc1a5 reversed these changes (Figure 6F). Moreover, knockdown of Slc1a5 led to decreased Bmal1 expression and disrupted rhythm in synchronized 3T3-L1 adipocytes (Figure 6G). We analyzed the WAT expression of glutamine- and methionine-related enzymes in obese and control mice. Obesity did not affect the expression of Bhmt, Mtr and Mat2a (Figure 6H). However, Gls and Gls2 (glutaminases, responsible for metabolism of glutamine to glutamate) were downregulated (Figure 6I). We also observed reduced expression of Glul (glutamine synthetase, catalyzing the formation of glutamine from glutamate) in obese mice (Figure 6I), consistent with a prior study . We argued that reduced glutamine and methionine pools in WAT were mainly attributed to impaired uptake transport mediated by SLC1A5. Contributions of the reductions in the three enzymes GLUL, GLS and GLS2 (catalyzing interconversion between glutamine and glutamate, a futile recycling) to a reduced glutamine pool may be negligible because of concurrent reductions in both glutamine and glutamate in obese WAT (Figure 3H & Supplementary excel file).
PPAR-γ is a critical nuclear receptor in obesity development and downregulated in obese conditions . Obesity-induced downregulation of PPAR-γ is likely attributed to induction of the proinflammatory cytokine TNF-a, a repressor of PPAR-γ . We confirmed that PPAR-γ expression was reduced in WAT of both mice and humans with obesity (Figure 7A-C). Also, PPAR-γ was rhythmically expressed in normal mice, but its rhythm was disrupted in obese mice (Figure 7A-B). We thus wondered whether PPAR-γ plays a role in obesity-repression of SLC1A5. ChIP-sequencing revealed PPAR-γ as a potential transcriptional regulator of SLC1A5 in adipose tissue (Figure 7D). We found knockdown of Ppar-γ decreased, whereas overexpression of Ppar-γ increased, the mRNA level of Slc1a5 in 3T3-L1 adipocytes (Figure 7E-F). Consistently, activation of PPAR-γ by rosiglitazone (a specific agonist) increased the mRNA and protein of SLC1A5 in 3T3-L1 adipocytes (Figure 7G). Rosiglitazone can also induce SLC1A5 expression in human adipocytes (Figure 7H). Furthermore, rosiglitazone induced the luciferase activity driven by a -1.4 kb Slc1a5 promoter in a dose-dependent manner (Figure 7I). These data indicated that PPAR-γ transcriptionally regulated the expression of SLC1A5. Moreover, truncation and mutation analyses identified a PPRE element (i.e., -1248/-1231 bp) in Slc1a5 promoter responsible for PPAR-γ action (Figure 7J). ChIP assays confirmed recruitment of WAT PPAR-γ to the PPRE element of Slc1a5 in normal mice (Figure 7K), however, such recruitment was reduced in obese mice (Figure 7K). In addition, PPAR-γ activation enhanced the uptake of [13C5]-glutamine and [2H3]-methionine into 3T3-L1 adipocytes (Figure 7L), accompanied by enhanced enrichments of both H3K27ac and H3K4me3, and increased Bmal1 expression (Figure S8). However, the effect of PPAR-γ activation on Bmal1 expression was attenuated by Slc1a5 knockdown (Figure S9). Taken together, PPAR-γ is down-regulated in obese WAT and is a direct transcriptional activator of Slc1a5 gene. Thereby, impaired PPAR-γ in WAT of obese mice causes reductions in SLC1A5 expression and in adipocyte uptake of both glutamine and methionine.
We have observed perturbance of adipocyte clock (a peripheral circadian clock) in obese conditions. Importantly, we have revealed PPAR-γ as an integrator of obesity and adipocyte clock. Obesity causes a reduction in PPAR-γ, which provokes downregulation of SLC1A5 (a direct target of PPAR-γ and a dual uptake transporter of glutamine and methionine) and reductions in adipocyte glutamine and methionine (two epigenetic activators of Bmal1), leading to disruption of Bmal1 and other clock genes and thus to impaired adipocyte clock function. The mechanisms for obesity-repression of PPAR-γ have been established in the literature. Obese adipose tissue releases the proinflammatory cytokine TNF-α which inhibits PPAR-γ expression via suppressing C/EBPδ (a transcriptional activator of PPAR-γ gene) and enhancing HDAC3 (a PPAR-γ corepressor) activity in the nucleus [54-57]. Additional mechanism involves free fatty acids elevated in obese adipose tissue [58,59]. Free fatty acids inhibit PPAR-γ expression in a TLR4/ER stress-dependent manner . In a prior study, adipocyte-specific deletion of Bmal1 results in obesity in mice due to an attenuated rhythm in food intake (i.e., increased food intake in daytime) caused by temporal changes in circulating polyunsaturated fatty acids and in expression of neurotransmitters for appetite regulation . Therefore, we propose a vicious cycle between circadian disruption and obesity development (Figure S10), in which obesity can be both a cause and an effect of circadian disruption, and vice versa. This vicious cycle appears to perpetuate and intensify obesity, providing an explanation for the increasing prevalence of obesity worldwide.
Impaired PPAR-γ in mouse WAT reduces SLC1A5 expression via a transcriptional mechanism. (A and B) PPAR-γ mRNA (A) and protein (B) in WAT derived from ob/ob, HFD and control mice at 6 circadian times. Data are mean ± SD (n = 5). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (C) PPAR-γ protein in the omental fat of obese and lean subjects. (D) Identification of PPAR-γ as a potential regulator of Slc1a5 gene in adipocytes according to Cistrome database (accession number: GSM2433448). All Cistrome data have been carefully curated and processed with a streamlined analysis pipeline and evaluated with comprehensive quality control metrics . (E) Effects of Ppar-γ knockdown (after 3T3-L1 cells differentiation) on Slc1a5 mRNA in 3T3-L1 adipocytes. Data are mean ± SD (n = 3). *p < 0.05 (t-test). (F) Effects of Ppar-γ overexpression on Slc1a5 mRNA in 3T3-L1 adipocytes. Data are mean ± SD (n = 3). *p < 0.05 (t-test). (G) SLC1A5 mRNA and protein in 3T3-L1 adipocytes treated with rosiglitazone (Rosi, a PPAR-γ agonist). Data are mean ± SD (n = 3). *p < 0.05 (t-test). (H) SLC1A5 mRNA and protein in human adipocytes treated with rosiglitazone. Data are mean ± SD (n = 3). *p < 0.05 (t-test). (I) Effects of rosiglitazone (1 and 2 μM) on the activities of Slc1a5-Luc reporter (-1400/+100 bp) in 3T3-L1 adipocytes. Data are mean ± SD (n = 6). *p < 0.05 (one-way ANOVA and Bonferroni post hoc test). (J) Effects of rosiglitazone (1 μM) on the activities of various versions of Slc1a5-Luc reporters (-1400/+100 bp, -1000/+100 bp and a PPRE-mutated version) in 3T3-L1 adipocytes. Data are mean ± SD (n = 6). *p < 0.05 (t-test). (K) ChIP assays showing recruitment of PPAR-γ protein to the PPRE element of Slc1a5 promoter in WAT of ob/ob, HFD and control mice at CT6. Data are mean ± SD (n = 3). *p < 0.05 at individual time points as determined by two-way ANOVA and Bonferroni post hoc test. (L) Effects of rosiglitazone on [13C5]-glutamine and [2H3]-methionine uptake into 3T3-L1 adipocytes. Data are mean ± SD (n = 3). *p < 0.05 (t-test).
Clock genes such as Rev-erbα/β, Per1/2, Cry1/2, Rorα and Dbp were also disrupted in obese adipose tissue in addition to Bmal1 (Figure 1 & Figure S2). It is most likely that downregulation of these clock genes is secondary to the change in BMAL1 expression as the transcription of these genes is directly driven by BMAL1 and they are BMAL1 target genes [36,37]. It is thus proposed that BMAL1 is a fundamental factor linking obesity to circadian clock (adipocyte clock). In addition to disrupted clock gene expression, obesity results in disturbance of feeding rhythm in mice (Figures 1 & 2), that phenocopies adipocyte-specific knockout of Bmal1 . Therefore, disturbed feeding rhythm in obesity is mainly attributed to impaired expression of adipocyte Bmal1, which has an effect on hypothalamic feeding centers via regulating circulating and hypothalamic polyunsaturated fatty acids . It is noteworthy that obesity is also associated with decreased locomotor activity in mice (Figure 2). However, SCN clock genes were unaffected in obese mice, consistent with prior studies [27,61,62]. We speculate that reduced locomotion in obesity is regulated by the nervous systems outside the SCN (e.g., the dopamine and endocannabinoid systems) [63-64,65].
Both PPAR-γ mRNA and protein oscillate with circadian time in adipose tissue in normal mice, however, the protein is phase-delayed about 4 hours compared to the mRNA (Figure 7A-B). This is probably due to a time lag in the translation from mRNA to protein, as noted previously for hepatic PPAR-γ . The PPAR-γ protein is in phase with the Slc1a5 mRNA, supporting PPAR-γ as a transcriptional activator of Slc1a5 and as a main driver of Slc1a5 rhythm (Figures 6A & 7B). Although the mechanism for circadian oscillations in PPAR-γ remains elusive, nocturnin (Noc, a circadian deadenylase and a post-transcriptional modifier of mRNAs) may regulate the circadian expression pattern of PPAR-γ as Noc-deficient mice show attenuated diurnal rhythm in Ppar-γ mRNA . We analyzed the expression of Noc in obese adipose tissue, and its level was unaltered (Figure S11). Therefore, the blunting of PPAR-γ rhythm in obese adipose tissue is most likely attributed to circadian inhibiting effects of TNF-α, whose level varies with the times of the day [56,68].
We have shown a reduced glutamine level in obese adipose tissue consistent with a recent report by Petrus et al . Reduction in adipocyte glutamine is mainly attributed to impaired expression of SLC1A5, which mediates cellular uptake of glutamine. By contrast, Petrus et al ascribed the reduction of glutamine to the lowered level of GLUL (glutamine synthetase), an enzyme generates glutamine by using glutamate as a substrate. However, glutamine and glutamate are interchangeable, and GLS (glutaminase, also reduced in obese adipose tissue) catalyzes the formation of glutamate from glutamine . If glutamine is decreased due to reduced catabolism of glutamate, one may expect to see an elevation in adipocyte glutamate. This contradicts with the observation that glutamate is considerably reduced in obese adipose tissue (Supplementary excel file). Therefore, we argue that impaired uptake of glutamine from the blood circulation to adipocytes accounts for the reductions in both glutamine and glutamate. Interestingly, glutamine treatment can ameliorate obesity and related comorbidities in mice (Figure 5), consistent with previous observations that glutamine supplementation reduces obesity with an improvement in systemic insulin action [70,71]. Likewise, methionine protects mice against developing obesity and related comorbidities. The main mechanism of their actions refers to the restoration of adipocyte clock which regulates energy balance (feeding rhythm) and obesity development (Figure 5). Since clock genes such as Rev-erbα act as inflammatory repressors , there is a possibility that clock restoration reduces adipocyte inflammation and contributes to the protecting effects of glutamine and methionine on obesity (a condition of chronic low-grade inflammation).
We focused on analysis of changes in clock genes in obese models because we intended to clarify the mechanisms by which obesity regulates circadian clock. We have revealed that obesity causes considerable disruption of clock genes in WAT. It is noteworthy that affected genes should not be limited to clock genes considering that H3K27ac and H3K4me3 are histone marks broadly distributed at active promoters in the genome. In fact, we found that non-clock genes such as Pgc-1α and Wnt-6 (with known H3K27ac or H3K4me3 distribution) are also downregulated in WAT of HFD-fed mice, accompanied by reduced enrichment of H3K27ac or H3K4me3 at the gene promoters (Figures S12, 4C). We have shown that H3K27ac and H3K4me3 are altered in response to the reductions of glutamine and methionine in obese models, while other histone marks potentially affecting the transcription of clock genes are unaffected (Figure S5). The specific effect of methionine change on H3K4me3 is consistent with prior studies in which methionine metabolism and the SAM level affect the expression of H3K4me3, while having no influences on other histone methylation marks such as H3K4me1, H3K4me2, H3K9me2, H3K9me3, H3K27me3, H3K36me3 and H3K79me3 [73-74,75,76]. However, the mechanism by which the methionine effect is specific to H3K4me3 remains unknown. Likewise, why the glutamine effect is specific to H3K27ac is unaddressed.
HFD is known to affect circadian gene expression in the tissues, particularly the liver [77-78,79]. To examine the impact of HFD-induced obesity (an obese condition) on circadian clock, we attempted to minimize the effects of HFD by feeding obese mice (that had been maintained on HFD for 10 weeks) and controls a normal chow diet for one week prior to gene expression analysis. Therefore, the obese condition rather than HFD is the direct cause of disturbances in adipocyte clock and circadian behaviors (Figures 1 & 2). Supporting this, mice fed an HFD for 3 days (in the absence of obesity) showed unaffected expression of clock genes in adipose tissue (Figure S13). In addition, ob/ob mice maintained on a normal chow diet developed obesity and showed the circadian phenotypes similar to those of HFD-induced obese mice (Figures 1 & 2), emphasizing an association of adipocyte cock with obese conditions but not with diets.
In conclusion, obesity disrupts adipocyte clock in mice. Mechanistically, obesity causes a reduction in PPAR-γ, which provokes downregulation of SLC1A5 (a direct target of PPAR-γ and a dual uptake transporter of glutamine and methionine) and reductions in adipocyte glutamine and methionine (two epigenetic activators of Bmal1), contributing to disruption of Bmal1 and other clock genes and thus to impaired adipocyte clock. Therefore, PPAR-γ integrates obesity and adipocyte clock, facilitating a vicious cycle between circadian disruption and obesity development.
Anti-BMAL1 (ab3350), anti-H3K27ac (ab4729), anti-H3K4me3 (ab8580), anti-H3K9me3 (ab8898) and anti-GAPDH (ab8245) antibodies were purchased from Abcam (Cambridge, UK). Anti-REV-ERBα (14506-1-AP), anti-DBP (12662-1-AP), anti-SLC1A5 (20350-1-AP), and anti-Histone H3 (17168-1-AP) antibodies were purchased from Proteintech (Wuhan, China). Anti-PPAR-γ (2443), anti-H3K9ac (9649), anti-H3K9me2 (4658), anti-H3K27me3 (9733) and anti-rabbit IgG (2729) were purchased from Cell Signaling Technology (Danvers, MA). Glutamine, methionine, C646 and MM-102 were purchased from Tsbiochem (Shanghai, China). [13C5]-glutamine and [2H3]-methionine were purchased from Zzbio (Shanghai, China). Dexamethasone was purchased from Sigma-Aldrich (St. Louis, MO). RNAiso Plus reagent, PrimeScript RT Master Mix and ChamQ Universal SYBR qPCR Master Mix were purchased from Vazyme (Nanjing, China). Beetle luciferin was purchased from Promega (Madison, WI). siRNAs for Slc1a5 and Ppar-γ, and negative control were purchase from IGE Biotechnology (Guangzhou, China). pcDNA 3.1, pcDNA 3.1-Ppar-γ, pcDNA 3.1-Slc1a5, pRL-TK, pGL4.0 and Slc1a5 luciferase reporters (-1400/+100 bp, -1000/+100 bp and a PPRE-mutated version) were obtained from Kaile Technologies (Guangzhou, China).
Male wild-type C57BL/6 mice were obtained from HFK Biotechnology (Beijing, China). Male leptin-deficient (ob/ob) mice (on a C57BL/6 background) were purchased from GemPharmatech (Nanjing, China). Mice were housed in a specific pathogen-free facility under a 12 h light/12 h dark cycle with free access to water and food. Diets were normal chow (D12450J, Research Diets, New Brunswick, NJ) and HFD (60% calories from fat, D12492, Research Diets, New Brunswick, NJ). The procedures for animal experiments were approved by the Institutional Animal Care and Use Committee of Guangzhou University of Chinese Medicine.
Eight-week-old ob/ob mice (maintained on a chow diet) were used as one of two obese models. The ob/ob and control wild-type mice were kept under constant darkness one day prior to sample collection. Mice (n = 5 per time point) were sacrificed at 6 circadian times (CT2, CT6, CT10, CT14, CT18 and CT22), and the SCN, serum, liver and perigonadal WAT samples were collected for gene, protein and/or metabolite analyses. For the second model, 4-week-old wild-type mice were fed HFD for 10 weeks to induce obesity, and switched to chow diet feeding for another 1 week. Control group of mice were kept on the chow diet. The two groups of mice were kept under constant darkness one day prior to sample collection. Samples were collected at a 4-h interval around the clock as described above for ob/ob mice. To evaluate the effect of adipocyte Bmal1 on obesity, wild-type mice were fed HFD for 10 weeks. After 8 weeks, mice were treated with AAV8.aP2.Bmal1 (2.5*1011 virus particles per mouse, Viraltherapy Technologies, Wuhan, China) or control vector via intraperitoneal injection at CT6. 4 weeks later, mice were sacrificed, and serum and WAT were collected for biochemical and expression analyses.
Omental fat was obtained from control subjects with a BMI of 18.5-23.9 kg/m2 and from obese subjects with a BMI of 32-60 kg/m2 through abdominal surgery (information provided in Table S1). Individuals with a history of malignancy, chronic inflammatory diseases such as ulcerative colitis and Crohn's disease, active infectious diseases, and drug or alcohol abuse were excluded. Fresh tissue samples were snap-frozen in liquid nitrogen and stored at -80℃ until RNA, protein and metabolite analyses. The individual samples for each group were pooled for Western blotting as the sample amount was limited. Research protocol was approved by the Medical Ethical Committee of the First Affiliated Hospital of Jinan University. An informed consent was obtained from every participant prior to the study.
Tissues were formalin-fixed and embedded in paraffin. 5 μm paraffin-embedded sections were stained with haematoxylin & eosin (H&E). Images were captured using a AXIO Imager M1 microscope (Carl Zeiss, Oberkochen, Germany). Adipocyte size was quantified using the ImageJ software based on six sections per mouse and five mice per group.
Triglycerides and cholesterol in serum samples were measured using the biochemical kits (Jiancheng Bioengineering Institute, Nanjing, China). SAM and acetyl-CoA were measured by using ELISA kits (Meimian Biotechnology, Jiangsu, China).
RNA extraction, cDNA synthesis, and qPCR reactions were performed as described previously . Data were normalized to the housekeeping gene (mouse Gapdh or human GAPDH). Primers are listed in Table S2.
Protein samples were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (10% or 12% gels) and transferred to a polyvinylidene difluoride membrane. The membrane was sequentially incubated with primary antibodies and a horseradish peroxidase-conjugated secondary antibody. Protein bands were visualized by using enhanced chemiluminescence and Omega Lum G imaging system (Aplegen, Pleasanton, CA). Protein levels were quantified with Fluorchem 5500 software (Fisher Scientific, Fair Lawn, NJ) and normalized to GAPDH or Histone H3 (Figure S14).
Adipose tissues were fixed in 4% paraformaldehyde, and then transferred to a grade series of sucrose solution (10%, 20% and 30%). The sections (6 μm thickness) were blocked with 5% BSA and 0.5% Triton X-100 in phosphate-buffered saline (PBS), and then incubated with anti-BMAL1 antibody. After washing with PBS, sections were incubated with Alexa Fluor 488-conjugated anti-rabbit secondary antibody. Sections were washed, mounted and imaged using a laser scanning microscope (Carl Zeiss, Oberkochen, Germany).
Mice were individually housed in running wheel cages (Lafayette Instrument, Lafayette, IN) that were placed in light-tight cabinets under a 12 h light/dark cycle. After acclimation to the system, mice were subjected to 3 weeks of continuous recording (1 week under a 12 h light/dark cycle and 2 weeks under constant darkness). Locomotor activity was analyzed using the ClockLab software (Actimetrics, Wilmette, IL). The data were separated, pooled and averaged on the basis of light/dark cycle.
Mice were individually housed in Promethion Metabolic Cages (Sable Systems, Las Vegas, NV) for one week before data collection. Indirect calorimetric and energy balance parameters including VO2 (oxygen consumption, in l/min), VCO2 (carbon dioxide expiration, in l/min) and food intake were assessed for four days. Energy expenditure (EE, in kJ/min) was calculated by using the Weir's equation (16.3 ×VO2 + 4.6 ×VCO2). Values of EE were normalized to the body weight raised to the power 0.75.
Metabolites in obese and control WAT were extracted with 50% methanol. After centrifugation at 4000 g for 20 min, the supernatant was subjected to analysis by a UHPLC‐QTOF/MS system consisting of an ExionLC UHPLC and a TripleTOF 5600+ mass spectrometer (AB SCIEX, Warrington, UK). Chromatographic separation was performed with an ACQUITY UPLC T3 column (100*2.1 mm, 1.8 μm, Waters, Milford, MA). The mobile phases consisted of 0.1% formic acid (phase A) and acetonitrile in 0.1% formic acid (phase B). The gradient elution program was 5% phase B for 0-0.5 min, 5-100% phase B for 0.5-7 min, 100% phase B for 7-8 min, 100-5% phase B for 8-8.1 min, and 5% phase B for 8.1-10 min. The flow rate was 0.4 ml/min. The capillary voltage was set at 5 kV in positive mode and -4.5 kV in negative mode. LC-MS raw data file were converted into a common format (.mzdata), followed by data processing (including peak picking, peak grouping and retention time correction) with the XCMS software. The HMDB and KEGG databases were used to annotate the metabolites. Metabolite quantitation and differential metabolite identification were performed using metaX software. Differential metabolites between obese and normal WAT were defined as such when Benjamini-Hochberg-adjusted p-value was < 0.05. KEGG pathways fulfilling the criterion of a hypergeometric p < 0.05 were defined as significantly enriched in differential metabolites.
ChIP assays were performed using a SimpleChIP plus Enzymatic Chromatin IP kit (Cell Signaling Technology, Beverly, MA) as described in our previous publication . Briefly, WAT was fixed with formaldehyde, followed by the addition of glycine. Sample was lysed and chromatin was sheared (to a size of 150-900 bp) using the micrococcal nuclease. Sheared chromatin was subjected to immunoprecipitation with anti-H3K4me3, anti-H3K27ac, anti-PPAR-γ or normal IgG, followed by incubation with protein G magnetic beads for 2 h. The immune complex was eluted to collect DNA fragments, and DNA fragments were purified using spin columns. The abundance of precipitated DNA fragments was determined by qPCR using the primers for targeted loci (Table S3).
3T3-L1 adipocytes (3T3-L1 cells were differentiated into adipocytes by treatment with 3-isobutyl-1-methylxanthine, dexamethasone, and insulin) were cultured in Dulbecco's modified Eagle's medium (DMEM). For histone acetylation inhibition, cells were pretreated with C646 (final concentration: 2, 5 or 10 μM) or vehicle for 12 h, and treated with 10 mM glutamine for another 24 h. For histone methylation inhibition, cells were pretreated with MM-102 (final concentration: 5, 20 or 50 μM) or vehicle for 12 h, and treated with 30 μM methionine for another 24 h. At the end of the experiments, cells were harvested for RNA and protein extraction.
Cell synchronization (serum shock) was performed as previously described . Briefly, 3T3-L1 adipocytes were cultured in DMEM containing 10% fetal bovine serum (FBS, Gibco, Waltham, MA). On the next day, the culture medium was replaced with a serum-free medium. 12-h later, 50% FBS was added for 2 h. The medium was then changed back to serum-free medium containing 2 mM glutamine, 30 μM methionine, 1 μM rosiglitazone or vehicle. Cells were harvested for RNA extraction and qPCR at 0, 4, 8, 12, 16, 20 and 24 h after synchronization.
U2OS cells stably overexpressed with BMAL1-Luc (a BMAL1 promoter-driven luciferase reporter) were used to assess transcription activity of BMAL1. U2OS cells were seeded into 35 mm dishes and maintained in DMEM containing 10% FBS. On next day, cells were pre-treated with glutamine- or methionine-free DMEM for 24 h, followed by incubation with a recording medium (10 mM HEPES, 3.5 g/l glucose, 0.35 g/l NaHCO3, 1% penicillin-streptomycin, 10% FBS, 0.1 mM luciferin and 100 nM dexamethasone in DMEM) containing 10 mM glutamine or 30 μM methionine. Luminescence data (counts/s) were collected by using a LumiCycle high-throughput luminometer (Actimetrics, Wilmette, IL).
Four-week-old wild-type mice were fed on HFD for 6 weeks. The mice continued on HFD and we initiated feeding of glutamine (1%, in drinking water) or methionine (0.3%, in drinking water) or pure drinking water for 4 weeks. Body weight was monitored weekly. Mice were sacrificed at CT6, and serum and WAT samples were collected for further experiments. We also performed short-term treatment experiments. Four-week-old wild-type mice were fed on HFD or chow diet for 8 weeks. The mice continued on HFD or chow and we initiated feeding of glutamine (1%, in drinking water) or methionine (0.3%, in drinking water) or pure drinking water for 2 weeks. Body weight was monitored and samples were collected as described above.
Luciferase reporter assays were performed as previously described . Briefly, 3T3-L1 adipocytes were co-transfected with Slc1a5 luciferase reporter (-1400/+100, -1000/+100 bp or PPRE-mutated version) and pRL-TK vector using a JetPRIME transfection kit. 24 h later, the culture medium was replaced with fresh medium containing 1 μM rosiglitazone (a PPAR-γ agonist) or vehicle. After another 24 h, cells were lysed, and luciferase activities were determined by using Dual-Luciferase Reporter Assay System (Promega, Madison, WI). Luciferase activity was expressed as relative luciferase unit (RLU).
Cellular uptake assays were performed with [13C5]-glutamine or [2H3]-methionine following a published protocol . Briefly, 3T3-L1 adipocytes were seeded into 6-well plates and maintained in DMEM containing 10% FBS. Cells were treated with rosiglitazone or siSlc1a5 or control for 24 h with or without following treatment with Slc1a5 plasmid for 24 h. Cells were washed with PBS and incubated with 2 ml isotope-containing medium (glutamine-free DMEM supplemented with 60 μM [13C5]-glutamine or methionine-free DMEM supplemented with 200 μM [2H3]-methionine) for 1 h. Cells were solubilized in 80% methanol. After centrifugation at 14000 g for 20 min, supernatant was collected and dried under nitrogen gas. Residues were resuspended in 80% methanol and subjected to UPLC-QTOF/MS analyses.
Glutamine and methionine (including [13C5]-glutamine and [2H3]-methionine) were quantified using a Waters UPLC-QTOF/MS system consisting an ACQUITY UPLC and a Xevo G2 QTOF mass spectrometer. Chromatographic separation was performed using a BEH Amide column (1.7 μm, 2.1 x 100 mm) (Waters, Milford, MA). The mobile phase consisted of 20 mM ammonium acetate in water (phase A) and acetonitrile in 0.1% formic acid (phase B). The flow rate was 0.2 ml/min. The gradient elution program was 5% phase B for 0-1 min, 5-10% phase B for 1-2 min, 10-50% phase B for 2-3.5 min, 50-95% phase B for 3.5-4 min, 95% phase B for 4-4.5 min, and 95-5% phase B for 4.5-5 min. The mass spectrometer was operated at the positive ion mode.
Data are recorded as mean ± standard deviation (SD). Means of two groups were compared using Student's t-test. Analysis of variance (ANOVA, one-way or two-way) was performed to compare means of more than two groups. In most cases, we performed two-way ANOVA to analyze whether the effect of obesity was significant or not. If so, post-hoc Bonferroni test was used for further group comparisons. The level of significance was set at p < 0.05 (*).
AAV8: adeno-associated virus serotype 8; Ac-CoA: acetyl coenzyme A; BHMT: betaine homocysteine methyltransferase; BMAL1: brain and muscle ARNT-like protein-1; CLOCK: circadian locomotor output cycles kaput; ChIP: chromatin immunoprecipitation; CRY: cryptochrome; CT: circadian time; DBP: albumin D-site-binding protein; DMEM: Dulbecco's Modified Eagle's Medium; E4BP4: E4 promoter-binding protein 4; EE: energy expenditure; FBS: fetal bovine serum; Gln: glutamine; GLUL: glutamine synthetase; GLS: glutaminase; H3K4: histone H3 lysine 4; H3K27: histone H3 lysine 27; HFD: high fat diet; Met: methionine; NPAS2: neuronal PAS domain protein 2; MAT2A: methionine adenosyltransferase II, alpha; MTR: 5-methyltetrahydrofolate-homocysteine methyltransferase; PER: period; qPCR: quantitative polymerase chain reaction; PPAR: peroxisome proliferator-activated receptor; PPRE: PPAR response element; ROR: RAR-related orphan receptor; SAM: S-adenosylmethionine; SLC: solute carrier; WAT: white adipose tissue.
Supplementary figures and tables.
This work was supported by the National Natural Science Foundation of China (No. 81722049 & 82003839) and Natural Science Foundation of Guangdong Province (No. 2019A1515110892 & 2021A1515012189).
B.W., and S.W. designed the study; S.W., Y.L., L.G., Z.Y., J.L., S.R., F.L. and J.C. performed experiments; S.W. and Y.L. collected and analyzed data; Z.D. and Z.W. collected clinical samples; B.W. and S.W. wrote the manuscript.
All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
The authors have declared that no competing interest exists.
1. WHO (2018) - Fact sheet - Obesity and overweight. (https://www.who.int/en/news-room/fact-sheets/detail/obesity-and-overweight)
2. Finkelstein EA, Khavjou OA, Thompson H, Trogdon JG, Pan L, Sherry B. et al. Obesity and severe obesity forecasts through 2030. Am J Prev Med. 2012;42(6):563-70
3. GBD 2017 Causes of Death Collaborators. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018 392,1736-88
4. Ng SW, Popkin BM. Time use and physical activity: a shift away from movement across the globe. Obes Rev. 2012;13(8):659-80
5. Sturm R, Hattori A. Morbid obesity rates continue to rise rapidly in the United States. Int J Obes (Lond). 2013;37(6):889-91
6. Reppert SM, Weaver DR. Coordination of circadian timing in mammals. Nature. 2002;418(6901):935-41
7. Mohawk JA, Green CB, Takahashi JS. Central and peripheral circadian clocks in mammals. Annu Rev Neurosci. 2012;35:445-62
8. King DP, Takahashi JS. Molecular genetics of circadian rhythms in mammals. Annu Rev Neurosci. 2000;23:713-42
9. Shearman LP, Sriram S, Weaver DR, Maywood ES, Chaves I, Zheng B. et al. Interacting molecular loops in the mammalian circadian clock. Science. 2000;288(5468):1013-9
10. Zhang R, Lahens NF, Ballance HI, Hughes ME, Hogenesch JB. A circadian gene expression atlas in mammals: implications for biology and medicine. Proc Natl Acad Sci U S A. 2014;111(45):16219-24
11. Eide EJ, Woolf MF, Kang H, Woolf P, Hurst W, Camacho F. et al. Control of mammalian circadian rhythm by CKIepsilon-regulated proteasome-mediated PER2 degradation. Mol Cell Biol. 2005;25(7):2795-807
12. Turek FW, Joshu C, Kohsaka A, Lin E, Ivanova G, McDearmon E. et al. Obesity and metabolic syndrome in circadian Clock mutant mice. Science. 2005;308(5724):1043-5
13. Touitou Y, Reinberg A, Touitou D. Association between light at night, melatonin secretion, sleep deprivation, and the internal clock: Health impacts and mechanisms of circadian disruption. Life Sci. 2017;173:94-106
14. Parsons MJ, Moffitt TE, Gregory AM, Goldman-Mellor S, Nolan PM, Poulton R. et al. Social jetlag, obesity and metabolic disorder: investigation in a cohort study. Int J Obes (Lond). 2015;39(5):842-8
15. Pan X, Jiang XC, Hussain MM. Impaired cholesterol metabolism and enhanced atherosclerosis in clock mutant mice. Circulation. 2013;128(16):1758-69
16. Lamia KA, Storch KF, Weitz CJ. Physiological significance of a peripheral tissue circadian clock. Proc Natl Acad Sci U S A. 2008;105(39):15172-7
17. Karlsson B, Knutsson A, Lindahl B. Is there an association between shift work and having a metabolic syndrome? Results from a population based study of 27,485 people. Occup Environ Med. 2001;58(11):747-52
18. Shi SQ, Ansari TS, McGuinness OP, Wasserman DH, Johnson CH. Circadian disruption leads to insulin resistance and obesity. Curr Biol. 2013;23(5):372-81
19. Kettner NM, Mayo SA, Hua J, Lee C, Moore DD, Fu L. Circadian Dysfunction Induces Leptin Resistance in Mice. Cell Metab. 2015;22(3):448-59
20. Paschos GK, Ibrahim S, Song WL, Kunieda T, Grant G, Reyes TM. et al. Obesity in mice with adipocyte-specific deletion of clock component Arntl. Nat Med. 2012;18(12):1768-77
21. Shimba S, Ishii N, Ohta Y, Ohno T, Watabe Y, Hayashi M. et al. Brain and muscle Arnt-like protein-1 (BMAL1), a component of the molecular clock, regulates adipogenesis. Proc Natl Acad Sci U S A. 2005;102(34):12071-6
22. Tahira K, Ueno T, Fukuda N, Aoyama T, Tsunemi A, Matsumoto S. et al. Obesity alters the expression profile of clock genes in peripheral blood mononuclear cells. Arch Med Sci. 2011;7(6):933-40
23. Pivovarova O, Gögebakan Ö, Sucher S, Groth J, Murahovschi V, Kessler K. et al. Regulation of the clock gene expression in human adipose tissue by weight loss. Int J Obes (Lond). 2016;40(6):899-906
24. Ando H, Yanagihara H, Hayashi Y, Obi Y, Tsuruoka S, Takamura T. et al. Rhythmic messenger ribonucleic acid expression of clock genes and adipocytokines in mouse visceral adipose tissue. Endocrinology. 2005;146(12):5631-6
25. Grosbellet E, Dumont S, Schuster-Klein C, Guardiola-Lemaitre B, Pevet P, Criscuolo F. et al. Circadian phenotyping of obese and diabetic db/db mice. Biochimie. 2016;124:198-206
26. Kaneko K, Yamada T, Tsukita S, Takahashi K, Ishigaki Y, Oka Y. et al. Obesity alters circadian expressions of molecular clock genes in the brainstem. Brain Res. 2009;1263:58-68
27. Kohsaka A, Laposky AD, Ramsey KM, Estrada C, Joshu C, Kobayashi Y. et al. High-fat diet disrupts behavioral and molecular circadian rhythms in mice. Cell Metab. 2007;6(5):414-21
28. Kersten S, Desvergne B, Wahli W. Roles of PPARs in health and disease. Nature. 2000;405(6785):421-4
29. Spiegelman BM. PPAR-gamma: adipogenic regulator and thiazolidinedione receptor. Diabetes. 1998;47(4):507-14
30. Soccio RE, Li Z, Chen ER, Foong YH, Benson KK, Dispirito JR. et al. Targeting PPARγ in the epigenome rescues genetic metabolic defects in mice. J Clin Invest. 2017;127(4):1451-62
31. Agarwal AK, Garg A. A novel heterozygous mutation in peroxisome proliferator-activated receptor-gamma gene in a patient with familial partial lipodystrophy. J Clin Endocrinol Metab. 2002;87(1):408-11
32. Hegele RA, Cao H, Frankowski C, Mathews ST, Leff T. PPARG F388L, a transactivation-deficient mutant, in familial partial lipodystrophy. Diabetes. 2002;51(12):3586-90
33. Laganà AS, Vitale SG, Nigro A, Sofo V, Salmeri FM, Rossetti P. et al. Pleiotropic Actions of Peroxisome Proliferator-Activated Receptors (PPARs) in Dysregulated Metabolic Homeostasis, Inflammation and Cancer: Current Evidence and Future Perspectives. Int J Mol Sci. 2016;17(7):999
34. Semple RK, Chatterjee VK, O'Rahilly S. PPAR gamma and human metabolic disease. J Clin Invest. 2006;116(3):581-9
35. Peters JM, Shah YM, Gonzalez FJ. The role of peroxisome proliferator-activated receptors in carcinogenesis and chemoprevention. Nat Rev Cancer. 2012;12(3):181-95
36. Hood S, Amir S. The aging clock: circadian rhythms and later life. J Clin Invest. 2017;127(2):437-46
37. Greco CM, Sassone-Corsi P. Circadian blueprint of metabolic pathways in the brain. Nat Rev Neurosci. 2019;20(2):71-82
38. Sánchez DI, González-Fernández B, Crespo I, San-Miguel B, Álvarez M, González-Gallego J. et al. Melatonin modulates dysregulated circadian clocks in mice with diethylnitrosamine-induced hepatocellular carcinoma. J Pineal Res. 2018;65(3):e12506
39. Sulli G, Manoogian ENC, Taub PR, Panda S. Training the Circadian Clock, Clocking the Drugs, and Drugging the Clock to Prevent, Manage, and Treat Chronic Diseases. Trends Pharmacol Sci. 2018;39(9):812-27
40. Wang S, Lin Y, Zhou Z, Gao L, Yang Z, Li F. et al. Circadian Clock Gene Bmal1 Regulates Bilirubin Detoxification: A Potential Mechanism of Feedback Control of Hyperbilirubinemia. Theranostics. 2019;9(18):5122-33
41. Dyar KA, Lutter D, Artati A, Ceglia NJ, Liu Y, Armenta D. et al. Atlas of Circadian Metabolism Reveals System-wide Coordination and Communication between Clocks. Cell. 2018;174(6):1571-85.e11
42. Farooqui AA, Horrocks LA, Farooqui T. Glycerophospholipids in brain: their metabolism, incorporation into membranes, functions, and involvement in neurological disorders. Chem Phys Lipids. 2000;106(1):1-29
43. Ishak Gabra MB, Yang Y, Li H, Senapati P, Hanse EA, Lowman XH. et al. Dietary glutamine supplementation suppresses epigenetically-activated oncogenic pathways to inhibit melanoma tumour growth. Nat Commun. 2020;11(1):3326
44. Roy DG, Chen J, Mamane V, Ma EH, Muhire BM, Sheldon RD. et al. Methionine Metabolism Shapes T Helper Cell Responses through Regulation of Epigenetic Reprogramming. Cell Metab. 2020;31(2):250-66.e9
45. Petrus P, Lecoutre S, Dollet L, Wiel C, Sulen A, Gao H. et al. Glutamine Links Obesity to Inflammation in Human White Adipose Tissue. Cell Metab. 2020;31(2):375-90.e11
46. Mei S, Qin Q, Wu Q, Sun H, Zheng R, Zang C. et al. Cistrome Data Browser: a data portal for ChIP-Seq and chromatin accessibility data in human and mouse. Nucleic Acids Res. 2017;45(D1):D658-62
47. Madiraju P, Pande SV, Prentki M, Madiraju SR. Mitochondrial acetylcarnitine provides acetyl groups for nuclear histone acetylation. Epigenetics. 2009;4(6):399-403
48. Lu C, Thompson CB. Metabolic regulation of epigenetics. Cell Metab. 2012;16(1):9-17
49. Zheng R, Wan C, Mei S, Qin Q, Wu Q, Sun H. et al. Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis. Nucleic Acids Res. 2019;47(D1):D729-35
50. Li QL, Wang DY, Ju LG, Yao J, Gao C, Lei PJ. et al. The hyper-activation of transcriptional enhancers in breast cancer. Clin Epigenetics. 2019;11(1):48
51. Karatas H, Townsend EC, Cao F, Chen Y, Bernard D, Liu L. et al. High-affinity, small-molecule peptidomimetic inhibitors of MLL1/WDR5 protein-protein interaction. J Am Chem Soc. 2013;135(2):669-82
52. Jiang Y, Rose AJ, Sijmonsma TP, Bröer A, Pfenninger A, Herzig S. et al. Mice lacking neutral amino acid transporter B(0)AT1 (Slc6a19) have elevated levels of FGF21 and GLP-1 and improved glycaemic control. Mol Metab. 2015;4(5):406-17
53. Wang W, Zou W. Amino Acids and Their Transporters in T Cell Immunity and Cancer Therapy. Mol Cell. 2020;80(3):384-95
54. Motawi TK, Shaker OG, Ismail MF, Sayed NH. Peroxisome Proliferator-Activated Receptor Gamma in Obesity and Colorectal Cancer: the Role of Epigenetics. Sci Rep. 2017;7(1):10714
55. Hotamisligil GS, Shargill NS, Spiegelman BM. Adipose expression of tumor necrosis factor-alpha: direct role in obesity-linked insulin resistance. Science. 1993;259(5091):87-91
56. Kudo M, Sugawara A, Uruno A, Takeuchi K, Ito S. Transcription suppression of peroxisome proliferator-activated receptor gamma2 gene expression by tumor necrosis factor alpha via an inhibition of CCAAT/ enhancer-binding protein delta during the early stage of adipocyte differentiation. Endocrinology. 2004;145(11):4948-56
57. Ye J. Regulation of PPARgamma function by TNF-alpha. Biochem Biophys Res Commun. 2008;374(3):405-8
58. Ma S, Yada K, Lee H, Fukuda Y, Iida A, Suzuki K. Taheebo Polyphenols Attenuate Free Fatty Acid-Induced Inflammation in Murine and Human Macrophage Cell Lines As Inhibitor of Cyclooxygenase-2. Front Nutr. 2017;4:63
59. Boden G. Obesity and free fatty acids. Endocrinol Metab Clin North Am. 2008;37(3):635-46 viii-ix
60. Nguyen MT, Chen A, Lu WJ, Fan W, Li PP, Oh DY. et al. Regulation of chemokine and chemokine receptor expression by PPARγ in adipocytes and macrophages. PLoS One. 2012;7(4):e34976
61. Ando H, Kumazaki M, Motosugi Y, Ushijima K, Maekawa T, Ishikawa E. et al. Impairment of peripheral circadian clocks precedes metabolic abnormalities in ob/ob mice. Endocrinology. 2011;152(4):1347-54
62. Prasai MJ, George JT, Scott EM. Molecular clocks, type 2 diabetes and cardiovascular disease. Diab Vasc Dis Res. 2008;5(2):89-95
63. Bjursell M, Gerdin AK, Lelliott CJ, Egecioglu E, Elmgren A, Törnell J. et al. Acutely reduced locomotor activity is a major contributor to Western diet-induced obesity in mice. Am J Physiol Endocrinol Metab. 2008;294(2):E251-60
64. Masubuchi S, Honma S, Abe H, Ishizaki K, Namihira M, Ikeda M. et al. Clock genes outside the suprachiasmatic nucleus involved in manifestation of locomotor activity rhythm in rats. Eur J Neurosci. 2000;12(12):4206-14
65. Mieda M, Sakurai T. Bmal1 in the nervous system is essential for normal adaptation of circadian locomotor activity and food intake to periodic feeding. J Neurosci. 2011;31(43):15391-6
66. Deng J, Guo L, Wu B. Circadian Regulation of Hepatic Cytochrome P450 2a5 by Peroxisome Proliferator-Activated Receptor γ. Drug Metab Dispos. 2018;46(11):1538-45
67. Green CB, Douris N, Kojima S, Strayer CA, Fogerty J, Lourim D. et al. Loss of Nocturnin, a circadian deadenylase, confers resistance to hepatic steatosis and diet-induced obesity. Proc Natl Acad Sci U S A. 2007;104(23):9888-93
68. Keller M, Mazuch J, Abraham U, Eom GD, Herzog ED, Volk HD. et al. A circadian clock in macrophages controls inflammatory immune responses. Proc Natl Acad Sci U S A. 2009;106(50):21407-12
69. Bernfeld E, Foster DA. Glutamine as an Essential Amino Acid for KRas-Driven Cancer Cells. Trends Endocrinol Metab. 2019;30(6):357-68
70. Laviano A, Molfino A, Lacaria MT, Canelli A, De Leo S, Preziosa I. et al. Glutamine supplementation favors weight loss in nondieting obese female patients. A pilot study. Eur J Clin Nutr. 2014;68(11):1264-6
71. Abboud KY, Reis SK, Martelli ME, Zordão OP, Tannihão F, de Souza AZZ. et al. Oral Glutamine Supplementation Reduces Obesity, Pro-Inflammatory Markers, and Improves Insulin Sensitivity in DIO Wistar Rats and Reduces Waist Circumference in Overweight and Obese Humans. Nutrients. 2019;11(3):536
72. Wang S, Li F, Lin Y, Wu B. Targeting REV-ERBα for therapeutic purposes: promises and challenges. Theranostics. 2020;10(9):4168-82
73. Li S, Swanson SK, Gogol M, Florens L, Washburn MP, Workman JL. et al. Serine and SAM Responsive Complex SESAME Regulates Histone Modification Crosstalk by Sensing Cellular Metabolism. Mol Cell. 2015;60(3):408-21
74. Singhal NK, Li S, Arning E, Alkhayer K, Clements R, Sarcyk Z. et al. Changes in Methionine Metabolism and Histone H3 Trimethylation Are Linked to Mitochondrial Defects in Multiple Sclerosis. J Neurosci. 2015;35(45):15170-86
75. Shiraki N, Shiraki Y, Tsuyama T, Obata F, Miura M, Nagae G. et al. Methionine metabolism regulates maintenance and differentiation of human pluripotent stem cells. Cell Metab. 2014;19(5):780-94
76. Ding W, Higgins DP, Yadav DK, Godbole AA, Pukkila-Worley R, Walker AK. Stress-responsive and metabolic gene regulation are altered in low S-adenosylmethionine. PLoS Genet. 2018;14(11):e1007812
77. Keleher MR, Zaidi R, Hicks L, Shah S, Xing X, Li D. et al. A high-fat diet alters genome-wide DNA methylation and gene expression in SM/J mice. BMC Genomics. 2018;19(1):888
78. Eckel-Mahan KL, Patel VR, de Mateo S, Orozco-Solis R, Ceglia NJ, Sahar S. et al. Reprogramming of the circadian clock by nutritional challenge. Cell. 2013;155(7):1464-78
79. Stokkan KA, Yamazaki S, Tei H, Sakaki Y, Menaker M. Entrainment of the circadian clock in the liver by feeding. Science. 2001;291(5503):490-3
80. Wang S, Lin Y, Yuan X, Li F, Guo L, Wu B. REV-ERBα integrates colon clock with experimental colitis through regulation of NF-κB/NLRP3 axis. Nat Commun. 2018;9(1):4246
81. Aoyagi T, Shimba S, Tezuka M. Characteristics of circadian gene expressions in mice white adipose tissue and 3T3-L1 adipocytes. J Health Sci. 2005;51:21-32
82. Wang Y, Bai C, Ruan Y, Liu M, Chu Q, Qiu L. et al. Coordinative metabolism of glutamine carbon and nitrogen in proliferating cancer cells under hypoxia. Nat Commun. 2019;10(1):201
Corresponding author: Baojian Wu, Ph.D. E-mail: bj.wucom; Pinghua Sun, Ph.D. E-mail: pinghuasunnycom; Zhiyong Dong, Ph.D. E-mail: dongzy2008edu.cn