High Fat Diet and In Utero Exposure to Maternal Obesity Disrupts Circadian Rhythm and Leads to Metabolic Programming of Liver in Rat Offspring

Abstract

The risk of obesity in adulthood is subject to programming beginning at conception. In animal models, exposure to maternal obesity and high fat diets influences the risk of obesity in the offspring. Among other long-term changes, offspring from obese rats develop hyperinsulinemia, hepatic steatosis, and lipogenic gene expression in the liver at weaning. However, the precise underlying mechanisms leading to metabolic dysregulation in the offspring remains unclear. Using a rat model of overfeeding-induced obesity, we previously demonstrated that exposure to maternal obesity from pre-conception to birth, is sufficient to program increased obesity risk in the offspring. Offspring of obese rat dams gain greater body weight and fat mass when fed high fat diet (HFD) as compared to lean dam. Since, disruptions of diurnal circadian rhythm are known to detrimentally impact metabolically active tissues such as liver, we examined the hypothesis that maternal obesity leads to perturbations of core clock components and thus energy metabolism in offspring liver. Offspring from lean and obese dams were examined at post-natal day 35, following a short (2 wk) HFD challenge. Hepatic mRNA expression of circadian (CLOCK, BMAL1, REV-ERBα, CRY, PER) and metabolic (PPARα, SIRT1) genes were strongly suppressed in offspring exposed to both maternal obesity and HFD. Using a mathematical model, we identified two distinct biological mechanisms that modulate PPARα mRNA expression: i) decreased mRNA synthesis rates; and ii) increased non-specific mRNA degradation rate. Moreover, our findings demonstrate that changes in PPARα transcription were associated with epigenomic alterations in H3K4me3 and H3K27me3 histone marks near the PPARα transcription start site. Our findings indicated that offspring from obese rat dams have detrimental alternations to circadian machinery that may contribute to impaired liver metabolism in response to HFD, specifically via reduced PPARα expression prior to obesity development.

About the Authors:

Sarah J. Borengasser

Affiliations Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America, Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Ping Kang

Affiliation: Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America

Jennifer Faske

Affiliation: Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America

Horacio Gomez-Acevedo

Affiliations Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America, Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Michael L. Blackburn

Affiliation: Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America

Thomas M. Badger

Affiliations Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America, Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Kartik Shankar

* E-mail: [email protected]

Affiliations Arkansas Children's Nutrition Center, Little Rock, Arkansas, United States of America, Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, Arkansas, United States of America

Introduction

The epidemic rise in obesity incidence across the age spectrum continues to be a daunting public health challenge. Most concerning is the prevalence of obesity in infants and children as obesity in childhood tracks strongly into adulthood [1]. It is estimated that infants in the United States (birth to 2 y) are twice as likely to be greater than the 95th percentile for weight-to-length measurements, indicating an increasing number of infants weigh more than is considered healthy [2]. Maternal obesity and excessive gestational weight gain have been identified as factors that contribute to increased obesity in the offspring. This is particularly significant, as over 60% of all pregnancies in the United States are in women who are either overweight or obese at conception [3]. Several studies have shown strong associations between maternal OB and obesity in the offspring, both in childhood [4]–[7] and in adulthood [8]–[10]. A recent study of 37,000 individuals showed greater risk of cardiovascular disease and pre-mature death in those born to OB women [11]. Likewise, the diminished risk of obesity in children born to OB women who lost weight prior to pregnancy, also strongly suggests a role for developmental programming [12], [13]. In animal models, exposure to maternal obesity unambiguously influences the risk of obesity in the offspring [14]–[18]. Collectively, the evidence suggests that susceptibility to obesity may begin prior to birth. Consequently, understanding how the intrauterine environment contributes to offspring obesity development is a pertinent health concern.

The Barker hypothesis offers a framework linking the in utero environment to long-term health outcomes in offspring [19]. While originally described in low-birth weight infants, it is now clear that maternal caloric excess and high fat diets also engender similar long-term programming of offspring obesity risk [20]. Both epidemiological [12], [13], [21] and animal [22]–[25] studies demonstrate intrauterine exposure to maternal obesity is associated with increased adiposity and weight gain. We, as well as others, have reported metabolic programming due to maternal nutrition which includes insulin resistance [7], [22], cardiovascular disease [26], [27], nonalcoholic fatty liver disease (NAFLD) [28], [29], decreased energy expenditure [30], and impaired ability to utilize fatty acids for energy in offspring [30]. Notably, our model limits the exposure of maternal obesity exclusively to the intrauterine environment through cross-fostering of offspring to lean surrogate dams at post-natal day (PND) 1 [22]. Our studies have also shown that impaired energy metabolism in obese dam offspring occurs prior to alterations in body weight and adiposity [28], [30]. Although many such advances have furthered our understanding of fetal programming, much of the molecular mechanisms leading to metabolic dysregulation remain to be elucidated.

A close relationship exists between metabolism and diurnal circadian rhythm, in which disruption promotes obesity and metabolic disease [31]–[35]. Perturbations in circadian rhythms are known to impact energetics, particularly in tissues regulating metabolism (considered peripheral clocks) such as liver, skeletal muscle, brown adipose, and white adipose [36]–[38]. Factors such as feeding, fasting, or type of diet can modulate rhythmicity and metabolism. More specifically, recent studies suggest that circadian rhythms are affected by gestational exposure to dietary manipulations such as protein restriction [39] or high fat consumption [40]. We have previously demonstrated that exposure to maternal obesity disrupted targets downstream of peroxisome proliferator-activated receptor (PPAR)α and AMPK which are key regulators responsible for orchestrating fatty acid oxidation prior to obesity development [28]. Since, PPARα expression oscillates in a circadian fashion throughout the day, and also acts as a circadian regulator [41]–[43], it remains unclear whether maternal obesity impairs PPARα signaling via disruption of circadian rhythm.

In the present study, we examined the hypothesis that intrauterine exposure to maternal obesity disrupts offspring circadian rhythm and liver metabolism following a short term high fat diet (HFD) challenge after weaning. First, we investigated mRNA expression of core clock components which included circadian locomotor output cycles kaput (CLOCK), brain and muscle ARNTL-like protein-1 (BMAL1), REV-ERBα, Cryptochromes (CRYs), and Periods (PERs). Second, we applied a mathematical model to identify distinct mechanisms that contribute to changes in PPARα mRNA expression due to maternal obesity and HFD. Lastly, we investigated histone modifications in the PPARα promoter to elucidate epigenetic mechanisms regulating PPARα expression. Our results demonstrate that offspring from obese rat dams have detrimental alterations to circadian machinery that contribute to impaired liver metabolism in response to high fat feeding.

Materials and Methods

Animals and chemicals

Female Sprague-Dawley rats (150–175 g) were obtained from Charles River Laboratories (Wilmington, MA). Animals were housed in an AAALAC-approved animal facility in a temperature and light controlled room (12 h light-12 h dark cycle). All experimental protocols were approved by the Institutional Animal Care and Use Committee at the University of Arkansas for Medical Sciences (Protocol # 2971). Unless specified, all chemicals were obtained from Sigma-Aldrich Chemical Co. (St. Louis, MO).

Experimental protocol

Virgin female Sprague-Dawley rats were intragastrically cannulated to receive total enteral nutrition (TEN) at age 8 wk and allowed to recover for 10 d as previously described [22], [44]–[47]. Rats were fed liquid TEN diets at either 155 kcal/kg3/4·d (referred to as lean dams) or at 220 kcal/kg3/4·d (40% excess calories, referred to as obese dams). We have previously reported body weights and body compositions of lean and obese dams [22]. TEN diets met National Research Council (NRC) nutrient recommendations as used previously by our group [22], [44], [46]–[50] and consisted of 20% protein (casein), 75% carbohydrate (dextrose and maltodextrin), and 5% fat (corn oil) . Infusion of diets was carried out for 3 wk allowing for precise control of both diet composition and caloric intake in a low-stress manner. Animals had ad libitum access to drinking water and body weights were measured three times per week. Following 3 wk of overfeeding to induce obesity in the 220 kcal/kg3/4·d group, lean and obese rats (N = 15/group) were allowed to mate for 1 wk. Each female rat was housed with one lean breeder male and allowed ad libitum access to AIN-93G diet during this period. After mating all female rats (lean and obese) received diets at 220 kcal/kg3/4·d (NRC recommended caloric intake for pregnancy in rats). All rats were allowed to give birth naturally. Numbers and sex of pups, birth weight, and crown-to-rump and anogenital distance were measured for each pup on PND1 as previously described [22], [28]. On PND1, four male and four female pups from each litter were cross-fostered to lean dams that had been previously time-impregnated to give birth on the same day as the obese dams receiving infusion diets. Cross-fostered dams were not cannulated and had ad libitum access to AIN-93G pelleted diets throughout lactation. Using this experimental paradigm, we ensured that offspring's exposure to any effects of maternal obesity was limited exclusively to the intrauterine environment [22]. Female offspring of lean and obese dams were used for separate experiments, and only data from male offspring are reported here. At PND21 male offspring were weaned onto either an AIN-93G (17% kcals fat) or high fat diet (45% kcals fat) for 2 wk. Offspring were weighed and sedated with carbon dioxide and euthanized via exsanguination in the fed condition at PND35 every 4 h over a 24 h period (N = 4/group for all time points except, N = 5/group for Lean-Con 2AM and Obese-HFD 6AM and N = 3/group for Lean-Con 10PM and Obese-Con 10PM, See Figure 1 for detailed offspring grouping). Each offspring group represents 2–4 distinct biological litters. At sacrifice, liver was weighed and immediately frozen in liquid nitrogen and stored at −70°C for later analyses. Serum was obtained by centrifugation of blood samples and stored at −20°C. Animal characteristics and serum parameters are shown in Table S2 in File S1 and Figure S1 in File S1.

[Figure omitted. See PDF.]

Figure 1. Schematic of experimental design, interventions and outcomes.

https://doi.org/10.1371/journal.pone.0084209.g001

Real-time RT-PCR

Total RNA was isolated from liver of offspring at PND35 (N = 3–5 rats per group per time point) using RNeasy mini columns (QIAGEN, Valencia, CA) including on-column DNase digestion. One microgram of total RNA was reverse transcribed using iScript cDNA synthesis kit (BioRad, Hercules, CA). Real-time PCR analysis was performed as described previously using an ABI Prism 7500 Fast instrument (Carlsbad, CA) [50], [51]. Area under the curve (AUC) was calculated using the trapezoidal rule. Gene specific primers were designed using Primer Express Software (Table S1 in File S1). Relative amounts of mRNA were quantified using a standard curve and normalized to the expression of SRP14 mRNA.

Mathematical Modeling

The mammalian circadian cycle has been computationally modeled by Leloup and Golbeter under normal physiological conditions [52], [53]. Here, we extended their model to analyze the effect of maternal obesity and HFD on PPARα mRNA expression. The relative mRNA expression level of PPARα was analyzed with a differential equation that has the form:where X(t) represented the relative gene expression level of PPARα. A more detailed description of the model is described in Table S3 in File S1. Parameter values were obtained by Monte Carlo simulations (n = 10,000) followed by a conjugate-gradient method implemented in Mathematica 8 (Wolfram Research Inc, Champaign, IL) as shown in Table S4 in File S1. To assess the reliability and robustness of the model, we performed a sensitivity analysis using Partial Rank Correlation Coefficient (PRCC) method at two switching points of the dark/light cycle that are critical for the circadian cycle, namely 6AM and 6PM with sample size set to N = 1000 [54].

Chromatin Immunoprecipitation (ChIP)

Histone modifications on the PPARα promoter were assessed using ChIP. Samples were processed according to a previously established protocol [55] and using a ChIP-IT enzymatic kit (Active Motif, Carlsbad, CA) with minor modifications as described previously [56]. Briefly, pools of liver samples (each pool representing 4–5 separate animals) at 6AM and 6PM, from offspring of lean or obese dams fed either control or HFD, were used for analyses. Samples were minced and fixed in a 1% formaldehyde solution and dounce homogenized. Nuclear isolation was performed and chromatin was sheared using a Covaris S220 focused-ultrasonicator (Woburn, MA). MagnaChIP protein A/G beads (Millipore, Billerica, MA) and chromatin were pre-blocked with 0.5% BSA and used to pre-clear the chromatin. Immunoprecipitation was performed using 2.5 µg of ChIP grade antibodies for H3K4me3 (Active Motif, Carlsbad, CA), H3K27me3 (Millipore, Billerica, MA), or matched nonspecific IgG (Millipore, Billerica, MA). Target binding regions on PPARα promoter were amplified via qPCR using an Eco Real-Time PCR System (Illumina, San Diego, CA), ±500 base pairs from the transcription start site (TSS) (See Table S1 in File S1 for sequences).

Statistical Analysis

Data are expressed as means ± SEM, significance was set at p<0.05. Differences in mRNA expression between offspring of lean and obese dams at each time point at PND35 were determined using two-tailed Student's t-test. Differences between offspring of lean and obese dams at PND35 fed AIN-93G or HFD for all other measured parameters were analyzed using two-way analysis of variance (ANOVA). Significant interactions identified by two-way ANOVA were followed by a one-way ANOVA and all pair-wise comparisons by Student-Newman-Keuls. Statistical analyses were performed using SigmaPlot 12.5 software (Systat Software Inc., San Jose, CA).

Results

Maternal obesity and high fat diet disrupt hepatic mRNA amplitude of clock machinery in offspring

We examined mRNA expression of the core clock components (CLOCK, BMAL1, REV-ERBα, Per1, Per2, Per3, Cry1, and Cry2) to evaluate if maternal obesity and/or exposure to HFD altered circadian rhythm over a 24 h period (Figure 2). The clock machinery forms a transcriptional-translational loop (TTL) where BMAL1 dimerizes with CLOCK forming a heterodimer that binds to E-box promoter elements to positively regulate CLOCK and clock controlled genes [32]. As expected, mRNA expression of CLOCK and BMAL1 resembled each other in their cyclic oscillations in which both peaked at 10AM and steadily declined until 6PM, when the lowest levels of expression occurred as depicted in Figure 2A and 2B. Maternal obesity did not alter rhythmicity or levels of mRNA expression of either CLOCK or BMAL1 in control-diet fed offspring of lean and obese dams. However, in HFD-fed offspring AUC for CLOCK and BMAL1 mRNA were lowest in offspring from obese dams indicated by a significant interaction of maternal obesity and HFD consumption (p<0.001) (Table 1).

[Figure omitted. See PDF.]

Figure 2. Hepatic mRNA expression of core circadian genes.

mRNA expression was assessed in lean and obese dam offspring (A–H) fed either control or high fat diets (HFD) from weaning through PND35 (N = 3–5 animals per group, Lean-Con 2AM and Obese-HFD 6AM (N = 5 per group), Lean-Con 10AM and Obese-Con 10PM (N = 3 per group), and all the remaining groups (N = 4)). Gene expression was assessed via real-time RT-PCR. All genes were normalized to SRP14 and expressed relative to Lean-Con at 6AM. Data are presented as mean ± SEM Statistical differences were determined using a Student's t test. * denotes significance, P<0.05.

https://doi.org/10.1371/journal.pone.0084209.g002

[Figure omitted. See PDF.]

Table 1. 24(AUC) (N = 3–5/group).

https://doi.org/10.1371/journal.pone.0084209.t001

Similar to CLOCK/BMAL1, Cry and Per also form a heterodimer, which acts as a negative regulator by impairing the action of CLOCK/BMAL1 to complete the TTL. As anticipated, Cry2 and Per 1, 2, and 3 all exhibited mRNA expression patterns that were antiphasic to that of CLOCK/BMAL1 as these targets peaked at 6PM and reached their lowest levels at 10AM (Figure 2D–G). Per2 mRNA expression of control-fed lean and obese dam offspring displayed slightly shifted expression, peaking occurred at 10PM as opposed to 6PM (Figure 2F). Similar to CLOCK/BMAL1, maternal obesity by itself did not alter mRNA expression of Cry (1, 2) and Per (1,2,3), but only in combination with HFD did maternal obesity show an effect of further reducing expression levels of Cry and Per (Table 1). Cry1 did not follow an anti-phasic pattern of mRNA expression to CLOCK/BMAL1, instead Cry1 was highest at 6AM and 2AM and was lowest at 2PM as shown in Figure 2C.

Two nuclear receptors are central to the regulation of the above mentioned clock components via control of BMAL1 transcription. Retinoic acid receptor-related orphan receptor (ROR)α is thought to activate BMAL1 transcription, while REV-ERBα suppresses it. In the present study, we found that RORα mRNA expression did not differ between groups (data not shown). However, REV-ERBα mRNA expression was robustly induced at 2 PM in lean-CON (∼114-fold) and obese-CON (∼95-fold) relative to 6AM (Figure 2H). However, there effects were severely blunted in high fat-fed offspring via a decrease of ∼81-fold in lean dam offspring and ∼75-fold in obese dam offspring in mRNA expression amplitude as compared to control-fed counterparts (Figure 2H). This dramatic HFD-induced reduction of REV-ERBα mRNA expression did not coincide with substantial changes in BMAL1 transcription in HFD-fed offspring suggesting BMAL1 transcription may be regulated via other mechanisms.

Disruption of circadian machinery alters PPARα in obese dam offspring

PPARα oscillates in a circadian manner regulating multiple pathways relating to hepatic lipid metabolism. Similar to the gene expression patterns of the clock components, PPARα did not differ between lean and obese dam offspring fed control diet (Figure 3A). However, AUC values for PPARα suggested that offspring of obese dams have an exacerbated response when fed HFD by displaying the lowest AUC values among groups, similar to the mRNA expression of core clock machinery (Table 1). These findings were consistent with our previous reports showing decreased activation of PPARα targets in offspring of obese dams [28].

[Figure omitted. See PDF.]

Figure 3. Circadian expression of PPAR-α mRNA.

Hepatic PPARα mRNA was assessed in lean and obese dam offspring fed A) control or B) HFD at PND35 (N = 3–5 animals per group, Lean-Con 2AM and Obese-HFD 6AM (N = 5 per group), Lean-Con 10AM and Obese-Con 10PM (N = 3 per group), and all the remaining groups (N = 4)). Statistical differences were determined using a Student's t test. * denotes significance, P<0.05. Mathematical model derived from relative PPARα mRNA expression of offspring of lean and obese dams fed C) control diet or D) HFD at PND35. Model fitting incorporated eight parameters representing distinct biological mechanisms. E) Optimal values for the parameter νp which represents rates of PPARα mRNA synthesis and F) δp which represents rates of PPARα mRNA degradation were determined by Monte Carlo simulations (n = 10,000) followed by a conjugate-gradient method implemented in Mathematica.

https://doi.org/10.1371/journal.pone.0084209.g003

Mathematical modeling of PPARα time changes

We utilized mathematical modeling of PPARα mRNA expression to identify mechanistic parameters that drive alterations in PPARα rhythm [52], [53]. The model consisted of a differential equation with eight parameters that could contribute to the expression or regulation of PPARα and are defined in Table S3 in File S1. The fitted model for PPARα relative mRNA expression is shown in Figure 3C and 3D. An optimal value for each parameter was determined by Monte Carlo simulations (n = 10,000) followed by a conjugated gradient method (Table S4 in File S1). Partial Rank Correlation Coefficients (PRCC) was performed to determine the predominant contribution of each parameter in the model.

PRCC revealed three parameters were consistently significant for the PPARα model (see Table S5 in File S1 for the PRCC values), which included vp, Kp, and δp (p<0.001 for all parameters). Figure 3E and 3F show the estimated values for vp and δp where vp represented the maximum rate of mRNA synthesis and δp non-specific rates of mRNA degradation. Decreased νp values, shown in Figure 3E, indicated that the maximum rate of PPARα mRNA synthesis is decreased due to maternal obesity. Figure 3F showed a stepwise increase of δp, which suggests that obese-HFD offspring have the highest levels of non-specific PPARα mRNA degradation. Kp represented the activation constant for enhancement of gene expression by BMAL1 (Table S5 in File S1).

Histone modifications near the PPARα transcription start site

The TTL, as mentioned previously, involves dynamic transcription factor binding to gene promoters of core clock machinery and metabolism-related targets over the course of the day. Assessing chromatin states offers insight into whether gene transcription was activated or repressed. Here, we focused on PPARα transcription at the start of the light (6AM) and dark (6PM) cycles. Accordingly, ChIP was performed on the PPARα promoter to assess whether the effects of maternal obesity and HFD on PPARα mRNA expression were related to epigenetic changes, specifically histone modifications. ChIP was performed for an active gene transcription mark (H3K4me3) and a silencing transcription mark (H3K27me3) in liver at 6AM and 6PM. Quantitative RT-PCR was performed ±500 base pairs from the PPARα TSS to assess the histone modification status of H3K4me3 and H3K27me3 as these two histone marks are bivalently located near the TSS as shown in Figure 4 and Figure 5.

[Figure omitted. See PDF.]

Figure 4. H3K4me3 enrichment on the PPARα promoter.

Chromatin immunoprecipitation for H3K4me3 was carried out and A) upstream (−500 base pairs) and B) downstream (+500 base pairs from TSS) regions of the PPAR-α gene were amplified in pools of liver samples (each pool represents 3–5 separate animals) from offspring of lean or obese dams fed either control or HFD (run in triplicate). Enrichment was determined by real time RT-PCR and normalized to input levels. Data are presented as mean ± SEM. Statistical differences were determined using a two-way ANOVA to examine the effects of maternal obesity and post-weaning HFD. Significant interactions were followed by one way ANOVA and Student-Newman-Keuls post hoc analyses (P<0.05). Bold values represent significant main effects and interactions and values with different letter superscripts are significantly different from each other (P<0.05).

https://doi.org/10.1371/journal.pone.0084209.g004

[Figure omitted. See PDF.]

Figure 5. H3K27me3 enrichment on the PPARα promoter.

Chromatin immunoprecipitation for H3K27me3 was carried out and A) upstream (−500 base pairs) and B) downstream (+500 base pairs from TSS) regions of the PPAR-α gene were amplified in pools of liver samples (each pool represents 3–5 separate animals) from offspring of lean or obese dams fed either control or HFD (run in triplicate). Enrichment was determined by real time RT-PCR and normalized to input levels. Data are presented as mean ± SEM. Statistical differences were determined using a two-way ANOVA to examine the effects of maternal obesity and post-weaning HFD (P<0.05). Significant interactions were followed by one way ANOVA and Student-Newman-Keuls post hoc analyses (P<0.05). Bold values represent significant main effects and interactions and values with different letter superscripts are significantly different from each other (P<0.05).

https://doi.org/10.1371/journal.pone.0084209.g005

H3K4me3 showed a divergent pattern in control-fed offspring of lean and obese dams, displaying higher enrichment occurring at 6AM and lower at 6PM, both upstream (−500 base pairs) and downstream (+500 base pairs) from the PPARα TSS (Figure 4). Lean and obese dam offspring fed HFD lacked differential H3K4me3 enrichment differences between 6AM and 6PM upstream from the TSS (Figure 4A). The downstream region from the TSS was similar to the upstream region as lean and obese dam offspring fed the control diet displayed differential enrichment (Figure 4A and 4B). These findings would suggest that control-fed offspring of lean and obese dams would demonstrate more rhythmic PPARα mRNA expression while HFD-fed groups would lack oscillatory expression. In fact, PPARα mRNA expression was reflective of these histone changes as control-fed offspring of lean and obese dams showed rhythmic expression that was substantially blunted in HFD-fed offspring of lean and obese dams (Figure 3A and 3B).

We also assessed the repressive histone mark, H3K27me3. There was a significant effect due to maternal obesity at 6PM +500 base pairs from the TSS, but this effect was driven by the obese-HFD group as shown in Figure 5B. Promoter chromatin enrichment of H3K27me3 was significantly higher in HFD-fed offspring of lean and obese dams at 6AM (Figure 5A). At 6PM, the obese-HFD group had higher enrichment levels as compared to all other groups (p<0.001) as shown in Figure 5B. In addition, at 6AM, enrichment of H3K27me3 at both ±500 base pairs from the PPARα TSS was increased in lean and obese dam offspring fed HFD (Figure 5A and 5B). This would suggest that obese-HFD should display the lowest PPARα mRNA expression as enrichment of H3K27me3 as elevated at both 6AM and 6PM and at both TSS locations which is concordant with the PPARα mRNA expression shown in Figure 3A and 3B.

Enhancer of zeste homolog (EZH2) is recognized as the primary histone methyltransferase within the polycomb repressive complex (PRC)-2 that mediates H3K27 trimethylation and has also been shown to co-immunoprecipitate with CLOCK and BMAL1 in mouse liver [57]. EZH2 mRNA expression is shown in Figure 6A and 6B. There were no differences in levels of expression between lean-con, obese-con, or lean-HFD. However, there was an increase in EZH2 expression at nearly every time point in obese-HFD offspring. AUC values in Table 1 further support that obese-HFD had the highest EZH2 levels. In addition, EZH2 has been shown to be associated with sirtuin 1 (SIRT1) [58], [59]; SIRT1 deletion leads to increased EZH2 expression levels. Our data coincided with this relationship as SIRT1 levels were reduced (Figure 6C and 6D) and EZH2 levels were increased (Figure 6A and 6B). Notably, EZH2 mRNA expression was also reflective of H3K27me3 at 6PM (Figure 5B) and was consistent with PPARα mRNA expression (Figure 3A and 3B).

[Figure omitted. See PDF.]

Figure 6. Circadian expression of EZH2 and SIRT1 mRNA.

Hepatic mRNA expression of EZH2 in A) Control-fed and B) HFD-fed offspring of lean and obese dams at PND35 (N = 3–5 animals per group, Lean-Con 2AM and Obese-HFD 6AM (N = 5 per group), Lean-Con 10AM and Obese-Con 10PM (N = 3 per group), and all the remaining groups (N = 4)). C) Hepatic mRNA expression of SIRT1 in Control-fed and D) HFD-fed offspring of lean and obese dams at PND35 (N = 3–5 animals per group, Lean-Con 2AM and Obese-HFD 6AM (N = 5 per group), Lean-Con 10AM and Obese-Con 10PM (N = 3 per group), and all the remaining groups (N = 4)). Data are presented as mean ± SEM. Statistical differences were determined using a Student's t test. * denotes significance, P<0.05.

https://doi.org/10.1371/journal.pone.0084209.g006

Discussion

The maintenance of diurnal circadian rhythm is necessary for normal metabolic function. Misalignment of activities of the molecular clock such as oscillatory disruption in the period, phase, or amplitude of core clock components are associated with the development of metabolic diseases, including obesity, diabetes, and cardiovascular dysfunction [37], [60]. Moreover, alterations in diet composition (high fat diets), timing of eating (night eating), or life-style factors (sleep deprivation or night-shift work) are known to alter circadian rhythms and impair metabolism. However, the influence of gestational experiences (such as maternal obesity) on circadian rhythm in the offspring remains poorly understood. Since, maternal diet and obesity during pregnancy have been shown to influence offspring metabolism, appetite, and adiposity; altered circadian rhythms could be causal mediators.

Although circadian disruptions in periodicity or phase were not apparent in the present report, our findings unequivocally demonstrated that oscillatory amplitude of both core clock regulators and circadian metabolic regulators were altered in offspring when challenged with HFD that was exacerbated in offspring of obese dams. In particular, exposure to both maternal obesity and 2 week post-weaning HFD challenge resulted in the most impairment of gene transcription of core clock machinery and metabolic and epigenetic targets. Changes in PPARα mRNA expression were specifically linked to decreased rates of mRNA synthesis and increased rates of degradation as determined by mathematical modeling. Additionally, epigenomic changes were evident via differences in enrichment of H3K4me3 and H3K27me3 histone marks on the PPARα promoter. In agreement with mRNA expression, the coupling of maternal obesity and HFD exposure led to the most dramatic changes in our mathematical model and histone mark enrichment.

In addition to its central role as a pleiotropic regulator of lipid metabolism, PPARα directly interacts with core clock components in a circadian fashion. Dietary challenges such as fasting and HFD lead to increased PPARα as a normal adaptive response aimed at shunting lipids toward oxidation and away from storage. In an ongoing study, where offspring were fasted for 24 h (at PND22), we observed an attenuated induction in PPARα mRNA and nuclear protein expression in obese dam offspring compared to lean dam counterparts (data not shown). These findings clearly indicate that maternal obesity impairs the underlying mechanisms regulating transcriptional induction of PPARα. Moreover, we have previously demonstrated that exposure to maternal obesity leads to reduced mRNA expression of PPARα target genes prior to the development of obesity or adiposity gains [28]. Consequently, offspring of obese dams are unable to adequately mount a response to metabolic demands that require mobilization of lipids (viz. fasting and high fat diets); specifically, via an inability to induce PPARα and its downstream targets. We have previously shown that another critical regulator of hepatic fatty acid oxidation, SIRT3, was reduced at both the transcript and mitochondrial protein levels in offspring of obese dams [30]. Hence, our studies indicate that complementary pathways involved in hepatic lipid metabolism (PPARα and SIRT3) are impaired by maternal obesity. These may, in concert, contribute to increased offspring susceptibility to nonalcoholic fatty liver disease at weaning [28], [30] and development of obesity in later life [22], [61]. Nevertheless, our previous work did not investigate mechanisms related to changes in PPARα transcription. Indeed, there are limited studies regarding control of PPARα transcription and most research has focused on PPARα-mediated transactivation of downstream targets.

The use of a mathematical model offered valuable insight into the magnitude that distinct biological mechanisms played in modifying PPARα mRNA expression. We identified that mRNA synthesis and degradation rates were mechanisms involved in modulating mRNA expression of PPARα (Figure 3E and 3F). Notably, obese-HFD offspring had the highest rates of degradation and the lowest PPARα mRNA expression. These findings suggest an imbalance between transcription and mRNA degradation. There are specific mechanisms associated with the coupling of initiation and decay of mRNA which include mRNA imprinting, alternative transcription start site, alternative splicing, alternative poly-adenylation, and promoter-regulated decay as recently described in a review by Haimovich et al. (2013) [62]. Each mechanism has different factor(s) that regulate mRNA stability and mRNA decay according to location, binding sites, and timing. It is likely that one or more of these specialized mechanisms are influenced by maternal obesity. The mathematical model was limited to the eight parameters that were chosen a priori; it is likely that other mechanisms, not included in the model, may have also contributed to changes in PPARα mRNA expression. Nevertheless, this model was chosen due to its capacity to detect parameters contributing to mRNA expression due to alterations in circadian rhythmicity. A limitation of the study, but not the model, was that we did not experimentally evaluate rates of PPARα mRNA synthesis and decay. While identifying the precise machinery involved with mRNA stability and degradation was outside of the scope of the current study, these certainly warrant further investigation.

The prevailing hypothesis in developmental programming is that gestational events alter epigenetic and epigenomic regulation of key genes [63]–[66]. While the initiating signals whereby maternal obesity influences epigenetic regulation of genes remains unclear, increasing information about epigenetically regulated metabolic targets genes is forthcoming. Embryonic development and differentiation are associated with large changes in the epigenomic landscape [67]. In fact, recent evidence suggests that key physiological patterns such as circadian rhythms are associated with highly coordinated alterations in epigenetic histone marks and transcription factors [60], [68]–[70]. Consequently, via ChIP, we assessed if maternal obesity and HFD led to differences in histone demarcation (H3K4me3 and H3K27me3) on the PPARα promoter that were associated with changes in gene expression. H3K4me3 is known to follow a rhythmic pattern of enrichment at TSSs and a lag time between H3K4me3 and permissive gene transcription has been reported [68] and our findings supported this notion (Figures 4A and 4B). Although our H3K4me3 findings offer valuable insight into the HFD-induced reduction in PPARα mRNA amplitude, they do not fully explain the exacerbated suppression of mRNA expression in the obese-HFD group. It has been reported that H3K4me3 may be more strongly correlated with increasing gene expression in low CpG promoters as compared to high CpG promoters [71]. PPARα is considered a high CpG promoter which may be a possible explanation for the reason H3K4me3 was not as indicative of HFD-induced changes in PPARα mRNA expression. However, the rhythmic pattern of PPARα mRNA expression can, in part, be explained by H3K4me3 enrichment near the TSS.

Polycomb repressive complexes (PRC)1 and PRC2 mediate the developmentally important repressive mark H3K27me3. PRC1 and PRC2 have crucial roles in pattern specification, organ development, and cellular proliferation and differentiation [72]–[74]. The histone methyltransferase EZH2 is part of PRC2 and the sole enzyme that methylates H3K27me3, but it has also been shown to be involved in the maintenance of circadian clock machinery. Knockdown of EZH2 disrupts circadian rhythm and has been linked to changes in H3K27me3 binding to clock gene promoters [57]. More recently, there is evidence that reduced SIRT1 leads to increased stability and expression of EZH2 [59]. Our data are consistent with this model, where in response to HFD, a decline in SIRT1 was associated with increased expression of EZH2 and thus increased H3K27me3.

Here we showed, H3K27me3 enrichment directly corresponded to the highly suppressive HFD effect on PPARα mRNA expression in both lean-HFD and obese-HFD offspring (Figure 5A and 5B). Interestingly, obese-HFD groups showed high levels of enrichment at both 6AM and 6PM as opposed to the lean-HFD rats which only had increased enrichment at 6AM (Figure 5A and 5B) suggesting PPARα mRNA expression may be more suppressed by HFD in obese-dam offspring. Indeed, PPARα mRNA expression was lowest in the obese-HFD group (Figure 3B and Table 1). These findings suggest that maternal obesity and HFD are more strongly related to epigenetic changes associated with gene silencing, rather than the repression of gene activation. These results are also in agreement with our modeling data which indicated increased rates of mRNA degradation. In physiological terms, this relationship makes sense as well, as increased EZH2 is associated with adipogenesis and reduced SIRT1 is linked to decreased energy expenditure. Although we did not measure markers of adipogenesis or energy expenditure in this study, we have recently demonstrated increased adipogenesis in white adipose tissue [61] and decreased energy expenditure in obese dam offspring [30]. Our epigenetic changes are in agreement with other reports that found maternal HFD [40], [75] and maternal protein restrictions [76], [77] led to histone modifications.

There are a few limitations to the current study. We did not monitor feeding behavior or cage activity which are both known to be associated with circadian rhythmicity [78]–[81]. We have previously shown that exposure to maternal obesity does not alter food intake in offspring on either a control or HFD [22]; however, it is possible that the timing of food intake or activity is altered which is associated with circadian disruption. In addition, the present studies focused solely on the effects of maternal obesity or HFD on offspring liver and did not examine effects on the central clock. Each of these certainly warrants further investigation in future studies. It is also worthy to note that recent studies have examined the influence of maternal HFD-induced obesity on gene regulation of circadian rhythms. These studies conducted in non-human primates also found a robust effect of maternal and post-natal HFD on expression of the Clock paralog, Npas2 [40] and SIRT1 expression and acetylation [82]. In contrast to our findings, these effects were mainly driven by maternal HFD and not by maternal obesity per se, as revealed by diet-reversal during gestation. While the precise reasons are unclear, differences in diet composition to produce obesity, degree of maternal adiposity, exposure to maternal diets during lactation, and species differences may certainly contribute to differences in findings to those observed in our model.

Collectively, our results indicated that the combination of maternal obesity and HFD led to the greatest disruption of core clock machinery and reductions of PPARα mRNA expression. Mathematical modeling revealed that exposure to maternal obesity led to decreased mRNA synthesis rates and a stepwise increase in mRNA degradation, with obese-HFD displaying the highest rates of mRNA degradation. Epigenomic changes in H3K4me3 appeared to contribute to the rhythmic expression of PPARα in control-fed groups and blunted rhythmic expression in HFD-fed groups while H3K27me3 appeared to play a greater role in HFD-induced effects in the presence or absence of maternal obesity. In conclusion, exposure to maternal obesity in utero and post-weaning HFD appear to “poise” offspring to be hyper-responsive to HFD promoting development of obesity in later life.

Supporting Information

[Figure omitted. See PDF.]

File S1.

SUPPLEMENTARY MATERIALS AND METHODS. FIGURE S1. Serum Parameters at 6AM. TABLE S1. Primers Sequences for Real-time RT-PCR Analyses. TABLE S2. Animal Characteristics of Offspring of Lean and Obese Dam Offspring at 6AM. TABLE S3. Model Parameters and Definitions. TABLE S4. Estimated PPARα Parameter Values. TABLE S5. Most Highly Correlated Model Parameter PRCC Values.

https://doi.org/10.1371/journal.pone.0084209.s001

(DOCX)

Acknowledgments

We thank Matt Ferguson, Hoy Pittman, Bobby Fay and other members of the ACNC-Animal Research Core Facility for their assistance with animal studies.

Author Contributions

Conceived and designed the experiments: SJB TMB KS. Performed the experiments: SJB PK JF . Analyzed the data: SJB HGA MLB. Wrote the paper: SJB HGA MLB TMB KS.

Citation: Borengasser SJ, Kang P, Faske J, Gomez-Acevedo H, Blackburn ML, Badger TM, et al. (2014) High Fat Diet and In Utero Exposure to Maternal Obesity Disrupts Circadian Rhythm and Leads to Metabolic Programming of Liver in Rat Offspring. PLoS ONE 9(1): e84209. https://doi.org/10.1371/journal.pone.0084209

References

1. Deckelbaum RJ, Williams CL (2001) Childhood obesity: the health issue. Obes Res 9 (Suppl 4) 239S–243S 10.1038/oby.2001.125 [doi].

2. Ogden CL, Carroll MD, Kit BK, Flegal KM (2012) Prevalence of obesity and trends in body mass index among US children and adolescents, 1999–2010. JAMA 307: 483–490 jama.2012.40 [pii];10.1001/jama.2012.40 [doi].

3. Flegal KM, Carroll MD, Kit BK, Ogden CL (2012) Prevalence of obesity and trends in the distribution of body mass index among US adults, 1999–2010. JAMA 307: 491–497 jama.2012.39 [pii];10.1001/jama.2012.39 [doi].

4. Whitaker RC (2004) Predicting preschooler obesity at birth: the role of maternal obesity in early pregnancy. Pediatrics 114: e29–e36.

5. Olson CM, Strawderman MS, Dennison BA (2009) Maternal weight gain during pregnancy and child weight at age 3 years. Matern Child Health J 13: 839–846 10.1007/s10995-008-0413-6 [doi].

6. Catalano PM, Farrell K, Thomas A, Huston-Presley L, Mencin P, et al. (2009) Perinatal risk factors for childhood obesity and metabolic dysregulation. Am J Clin Nutr 90: 1303–1313 ajcn.2008.27416 [pii];10.3945/ajcn.2008.27416 [doi].

7. Catalano PM, Presley L, Minium J, Hauguel-de MS (2009) Fetuses of obese mothers develop insulin resistance in utero. Diabetes Care 32: 1076–1080 dc08-2077 [pii];10.2337/dc08-2077 [doi].

8. Reynolds RM, Osmond C, Phillips DI, Godfrey KM (2010) Maternal BMI, parity, and pregnancy weight gain: influences on offspring adiposity in young adulthood. J Clin Endocrinol Metab 95: 5365–5369 jc.2010-0697 [pii];10.1210/jc.2010-0697 [doi].

9. Lawlor DA, Smith GD, O'Callaghan M, Alati R, Mamun AA, et al. (2007) Epidemiologic evidence for the fetal overnutrition hypothesis: findings from the mater-university study of pregnancy and its outcomes. Am J Epidemiol 165: 418–424 kwk030 [pii];10.1093/aje/kwk030 [doi].

10. Jaaskelainen A, Pussinen J, Nuutinen O, Schwab U, Pirkola J, et al. (2011) Intergenerational transmission of overweight among Finnish adolescents and their parents: a 16-year follow-up study. Int J Obes (Lond) 35: 1289–1294 ijo2011150 [pii];10.1038/ijo.2011.150 [doi].

11. Reynolds RM, Allan KM, Raja EA, Bhattacharya S, McNeill G, et al. (2013) Maternal obesity during pregnancy and premature mortality from cardiovascular event in adult offspring: follow-up of 1 323 275 person years. BMJ 347: f4539.

12. Kral JG, Biron S, Simard S, Hould FS, Lebel S, et al. (2006) Large maternal weight loss from obesity surgery prevents transmission of obesity to children who were followed for 2 to 18 years. Pediatrics 118: e1644–e1649 118/6/e1644 [pii];10.1542/peds.2006-1379 [doi].

13. Smith J, Cianflone K, Biron S, Hould FS, Lebel S, et al. (2009) Effects of maternal surgical weight loss in mothers on intergenerational transmission of obesity. J Clin Endocrinol Metab 94: 4275–4283 jc.2009-0709 [pii];10.1210/jc.2009-0709 [doi].

14. Ainge H, Thompson C, Ozanne SE, Rooney KB (2011) A systematic review on animal models of maternal high fat feeding and offspring glycaemic control. Int J Obes (Lond) 35: 325–335 ijo2010149 [pii];10.1038/ijo.2010.149 [doi].

15. Nathanielsz PW, Poston L, Taylor PD (2007) In utero exposure to maternal obesity and diabetes: animal models that identify and characterize implications for future health. Clin Perinatol 34: 515–26, v S0095-5108(07)00073-5 [pii];10.1016/j.clp.2007.09.005 [doi].

16. Heerwagen MJ, Miller MR, Barbour LA, Friedman JE (2010) Maternal obesity and fetal metabolic programming: a fertile epigenetic soil. Am J Physiol Regul Integr Comp Physiol 299: R711–R722 ajpregu.00310.2010 [pii];10.1152/ajpregu.00310.2010 [doi].

17. Alfaradhi MZ, Ozanne SE (2011) Developmental programming in response to maternal overnutrition. Front Genet 2: 27 10.3389/fgene.2011.00027 [doi].

18. Li M, Sloboda DM, Vickers MH (2011) Maternal obesity and developmental programming of metabolic disorders in offspring: evidence from animal models. Exp Diabetes Res 2011: 592408 10.1155/2011/592408 [doi].

19. Barker DJ (2007) The origins of the developmental origins theory. J Intern Med 261: 412–417 JIM1809 [pii];10.1111/j.1365-2796.2007.01809.x [doi].

20. Catalano PM (2003) Obesity and pregnancy–the propagation of a viscous cycle? J Clin Endocrinol Metab 88: 3505–3506.

21. Villamor E, Cnattingius S (2006) Interpregnancy weight change and risk of adverse pregnancy outcomes: a population-based study. Lancet 368: 1164–1170 S0140-6736(06)69473-7 [pii];10.1016/S0140-6736(06)69473-7 [doi].

22. Shankar K, Harrell A, Liu X, Gilchrist JM, Ronis MJ, et al. (2008) Maternal obesity at conception programs obesity in the offspring. Am J Physiol Regul Integr Comp Physiol 294: R528–R538 00316.2007 [pii];10.1152/ajpregu.00316.2007 [doi].

23. Levin BE, Govek E (1998) Gestational obesity accentuates obesity in obesity-prone progeny. Am J Physiol 275: R1374–R1379.

24. Samuelsson AM, Matthews PA, Argenton M, Christie MR, McConnell JM, et al. (2008) Diet-induced obesity in female mice leads to offspring hyperphagia, adiposity, hypertension, and insulin resistance: a novel murine model of developmental programming. Hypertension 51: 383–392 HYPERTENSIONAHA.107.101477 [pii];10.1161/HYPERTENSIONAHA.107.101477 [doi].

25. Bayol SA, Simbi BH, Bertrand JA, Stickland NC (2008) Offspring from mothers fed a ‘junk food’ diet in pregnancy and lactation exhibit exacerbated adiposity that is more pronounced in females. J Physiol 586: 3219–3230 jphysiol.2008.153817 [pii];10.1113/jphysiol.2008.153817 [doi].

26. Barker DJ (1995) Fetal origins of coronary heart disease. BMJ 311: 171–174.

27. Singhal A, Lucas A (2004) Early origins of cardiovascular disease: is there a unifying hypothesis? Lancet 363: 1642–1645 10.1016/S0140-6736(04)16210-7 [doi];S0140-6736(04)16210-7 [pii].

28. Shankar K, Kang P, Harrell A, Zhong Y, Marecki JC, et al. (2010) Maternal overweight programs insulin and adiponectin signaling in the offspring. Endocrinology 151: 2577–2589 en.2010-0017 [pii];10.1210/en.2010-0017 [doi].

29. McCurdy CE, Bishop JM, Williams SM, Grayson BE, Smith MS, et al. (2009) Maternal high-fat diet triggers lipotoxicity in the fetal livers of nonhuman primates. J Clin Invest 119: 323–335 32661 [pii];10.1172/JCI32661 [doi].

30. Borengasser SJ, Lau F, Kang P, Blackburn ML, Ronis MJ, et al. (2011) Maternal obesity during gestation impairs fatty acid oxidation and mitochondrial SIRT3 expression in rat offspring at weaning. PLoS One 6: e24068 10.1371/journal.pone.0024068 [doi];PONE-D-11-05564 [pii].

31. Shi SQ, Ansari TS, McGuinness OP, Wasserman DH, Johnson CH (2013) Circadian Disruption Leads to Insulin Resistance and Obesity. Curr Biol S0960-9822(13)00083-3 [pii];10.1016/j.cub.2013.01.048 [doi].

32. Albrecht U (2012) Timing to perfection: the biology of central and peripheral circadian clocks. Neuron 74: 246–260 S0896-6273(12)00332-7 [pii];10.1016/j.neuron.2012.04.006 [doi].

33. Asher G, Schibler U (2011) Crosstalk between components of circadian and metabolic cycles in mammals. Cell Metab 13: 125–137 S1550-4131(11)00007-6 [pii];10.1016/j.cmet.2011.01.006 [doi].

34. Delezie J, Challet E (2011) Interactions between metabolism and circadian clocks: reciprocal disturbances. Ann N Y Acad Sci 1243: 30–46 10.1111/j.1749-6632.2011.06246.x [doi].

35. Coomans CP, van den Berg SA, Houben T, van Klinken JB, van den Berg R, et al. (2013) Detrimental effects of constant light exposure and high-fat diet on circadian energy metabolism and insulin sensitivity. FASEB J fj.12-210898 [pii];10.1096/fj.12-210898 [doi].

36. Richards J, Gumz ML (2012) Advances in understanding the peripheral circadian clocks. FASEB J fj.12-203554 [pii];10.1096/fj.12-203554 [doi].

37. Bass J, Takahashi JS (2010) Circadian integration of metabolism and energetics. Science 330: 1349–1354 330/6009/1349 [pii];10.1126/science.1195027 [doi].

38. Yang X, Downes M, Yu RT, Bookout AL, He W, et al. (2006) Nuclear receptor expression links the circadian clock to metabolism. Cell 126: 801–810 S0092-8674(06)00978-0 [pii];10.1016/j.cell.2006.06.050 [doi].

39. Sutton GM, Centanni AV, Butler AA (2010) Protein malnutrition during pregnancy in C57BL/6J mice results in offspring with altered circadian physiology before obesity. Endocrinology 151: 1570–1580 en.2009-1133 [pii];10.1210/en.2009-1133 [doi].

40. Suter M, Bocock P, Showalter L, Hu M, Shope C, et al. (2011) Epigenomics: maternal high-fat diet exposure in utero disrupts peripheral circadian gene expression in nonhuman primates. FASEB J 25: 714–726 fj.10-172080 [pii];10.1096/fj.10-172080 [doi].

41. Oishi K, Shirai H, Ishida N (2008) PPARalpha is involved in photoentrainment of the circadian clock. Neuroreport 19: 487–489 10.1097/WNR.0b013e3282f7968f [doi];00001756-200803050-00020 [pii].

42. Oishi K, Shirai H, Ishida N (2005) CLOCK is involved in the circadian transactivation of peroxisome-proliferator-activated receptor alpha (PPARalpha) in mice. Biochem J 386: 575–581 BJ20041150 [pii];10.1042/BJ20041150 [doi].

43. Canaple L, Rambaud J, Dkhissi-Benyahya O, Rayet B, Tan NS, et al. (2006) Reciprocal regulation of brain and muscle Arnt-like protein 1 and peroxisome proliferator-activated receptor alpha defines a novel positive feedback loop in the rodent liver circadian clock. Mol Endocrinol 20: 1715–1727 me.2006-0052 [pii];10.1210/me.2006-0052 [doi].

44. Baumgardner JN, Shankar K, Hennings L, Badger TM, Ronis MJ (2008) A new model for nonalcoholic steatohepatitis in the rat utilizing total enteral nutrition to overfeed a high-polyunsaturated fat diet. Am J Physiol Gastrointest Liver Physiol 294: G27–G38 00296.2007 [pii];10.1152/ajpgi.00296.2007 [doi].

45. Shankar K, Harrell A, Kang P, Singhal R, Ronis MJ, et al. (2010) Carbohydrate-responsive gene expression in the adipose tissue of rats. Endocrinology 151: 153–164.

46. Badger TM, Crouch J, Irby D, Hakkak R, Shahare M (1993) Episodic excretion of ethanol during chronic intragastric ethanol infusion in the male rat: continuous vs. cyclic ethanol and nutrient infusions. J Pharmacol Exp Ther 264: 938–943.

47. Badger TM, Ronis MJ, Lumpkin CK, Valentine CR, Shahare M, et al. (1993) Effects of chronic ethanol on growth hormone secretion and hepatic cytochrome P450 isozymes of the rat. J Pharmacol Exp Ther 264: 438–447.

48. Korourian S, Hakkak R, Ronis MJ, Shelnutt SR, Waldron J, et al. (1999) Diet and risk of ethanol-induced hepatotoxicity: carbohydrate-fat relationships in rats. Toxicol Sci 47: 110–117.

49. Shankar K, Liu X, Singhal R, Chen JR, Nagarajan S, et al. (2008) Chronic ethanol consumption leads to disruption of vitamin D3 homeostasis associated with induction of renal 1,25 dihydroxyvitamin D3-24-hydroxylase (CYP24A1). Endocrinology 149: 1748–1756 en.2007-0903 [pii];10.1210/en.2007-0903 [doi].

50. Shankar K, Hidestrand M, Liu X, Xiao R, Skinner CM, et al. (2006) Physiologic and genomic analyses of nutrition-ethanol interactions during gestation: Implications for fetal ethanol toxicity. Exp Biol Med (Maywood) 231: 1379–1397 231/8/1379 [pii].

51. Saben J, Zhong Y, Gomez-Acevedo H, Thakali KM, Borengasser SJ, et al. (2013) Early growth response protein-1 mediates lipotoxicity-associated placental inflammation: role in maternal obesity. Am J Physiol Endocrinol Metab 305: E1–14 10.1152/ajpendo.00076.2013 [doi].

52. Leloup JC, Goldbeter A (2004) Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms. J Theor Biol 230: 541–562 10.1016/j.jtbi.2004.04.040 [doi];S0022519304002437 [pii].

53. Leloup JC, Goldbeter A (2003) Toward a detailed computational model for the mammalian circadian clock. Proc Natl Acad Sci U S A 100: 7051–7056 10.1073/pnas.1132112100 [doi];1132112100 [pii].

54. Marino S, Hogue IB, Ray CJ, Kirschner DE (2008) A methodology for performing global uncertainty and sensitivity analysis in systems biology. J Theor Biol 254: 178–196 S0022-5193(08)00189-6 [pii];10.1016/j.jtbi.2008.04.011 [doi].

55. Schmidt D, Wilson MD, Spyrou C, Brown GD, Hadfield J, et al. (2009) ChIP-seq: using high-throughput sequencing to discover protein-DNA interactions. Methods 48: 240–248 S1046-2023(09)00047-4 [pii];10.1016/j.ymeth.2009.03.001 [doi].

56. Singhal R, Shankar K, Badger TM, Ronis MJ (2008) Estrogenic status modulates aryl hydrocarbon receptor–mediated hepatic gene expression and carcinogenicity. Carcinogenesis 29: 227–236 bgm288 [pii];10.1093/carcin/bgm288 [doi].

57. Etchegaray JP, Yang X, DeBruyne JP, Peters AH, Weaver DR, et al. (2006) The polycomb group protein EZH2 is required for mammalian circadian clock function. J Biol Chem 281: 21209–21215 M603722200 [pii];10.1074/jbc.M603722200 [doi].

58. Lu L, Li L, Lu X, Wu XS, Liu DP, et al. (2011) Inhibition of SIRT1 increases EZH2 protein level and enhances the repression of EZH2 on target gene expression. Chin Med Sci J 26: 77–84 9736 [pii].

59. O'Hagan HM, Mohammad HP, Baylin SB (2008) Double strand breaks can initiate gene silencing and SIRT1-dependent onset of DNA methylation in an exogenous promoter CpG island. PLoS Genet 4: e1000155 10.1371/journal.pgen.1000155 [doi].

60. Feng D, Lazar MA (2012) Clocks, metabolism, and the epigenome. Mol Cell 47: 158–167 S1097-2765(12)00551-5 [pii];10.1016/j.molcel.2012.06.026 [doi].

61. Borengasser SJ, Zhong Y, Kang P, Lindsey F, Ronis MJ, et al. (2013) Maternal Obesity Enhances White Adipose Tissue Differentiation and Alters Genome-Scale DNA Methylation in Male Rat Offspring. Endocrinology en.2012-2255 [pii];10.1210/en.2012-2255 [doi].

62. Haimovich G, Choder M, Singer RH, Trcek T (2013) The fate of the messenger is pre-determined: A new model for regulation of gene expression. Biochim Biophys Acta S1874-9399(13)00008-4 [pii];10.1016/j.bbagrm.2013.01.004 [doi].

63. Smallwood SA, Tomizawa S, Krueger F, Ruf N, Carli N, et al. (2011) Dynamic CpG island methylation landscape in oocytes and preimplantation embryos. Nat Genet 43: 811–814 ng.864 [pii];10.1038/ng.864 [doi].

64. Maltepe E, Bakardjiev AI, Fisher SJ (2010) The placenta: transcriptional, epigenetic, and physiological integration during development. J Clin Invest 120: 1016–1025 41211 [pii];10.1172/JCI41211 [doi].

65. Meissner A (2010) Epigenetic modifications in pluripotent and differentiated cells. Nat Biotechnol 28: 1079–1088 nbt.1684 [pii];10.1038/nbt.1684 [doi].

66. Hawkins RD, Hon GC, Lee LK, Ngo Q, Lister R, et al. (2010) Distinct epigenomic landscapes of pluripotent and lineage-committed human cells. Cell Stem Cell 6: 479–491 S1934-5909(10)00147-5 [pii];10.1016/j.stem.2010.03.018 [doi].

67. Zhu J, Adli M, Zou JY, Verstappen G, Coyne M, et al. (2013) Genome-wide chromatin state transitions associated with developmental and environmental cues. Cell 152: 642–654 S0092-8674(12)01555-3 [pii];10.1016/j.cell.2012.12.033 [doi].

68. Koike N, Yoo SH, Huang HC, Kumar V, Lee C, et al. (2012) Transcriptional architecture and chromatin landscape of the core circadian clock in mammals. Science 338: 349–354 science.1226339 [pii];10.1126/science.1226339 [doi].

69. Gallou-Kabani C, Vige A, Junien C (2007) Lifelong circadian and epigenetic drifts in metabolic syndrome. Epigenetics 2: 137–146 4897 [pii].

70. Masri S, Zocchi L, Katada S, Mora E, Sassone-Corsi P (2012) The circadian clock transcriptional complex: metabolic feedback intersects with epigenetic control. Ann N Y Acad Sci 1264: 103–109 10.1111/j.1749-6632.2012.06649.x [doi].

71. Karlic R, Chung HR, Lasserre J, Vlahovicek K, Vingron M (2010) Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A 107: 2926–2931 0909344107 [pii];10.1073/pnas.0909344107 [doi].

72. Margueron R, Reinberg D (2011) The Polycomb complex PRC2 and its mark in life. Nature 469: 343–349 nature09784 [pii];10.1038/nature09784 [doi].

73. Ku M, Koche RP, Rheinbay E, Mendenhall EM, Endoh M, et al. (2008) Genomewide analysis of PRC1 and PRC2 occupancy identifies two classes of bivalent domains. PLoS Genet 4: e1000242 10.1371/journal.pgen.1000242 [doi].

74. O'Meara MM, Simon JA (2012) Inner workings and regulatory inputs that control Polycomb repressive complex 2. Chromosoma 121: 221–234 10.1007/s00412-012-0361-1 [doi].

75. Aagaard-Tillery KM, Grove K, Bishop J, Ke X, Fu Q, et al. (2008) Developmental origins of disease and determinants of chromatin structure: maternal diet modifies the primate fetal epigenome. J Mol Endocrinol 41: 91–102 JME-08-0025 [pii];10.1677/JME-08-0025 [doi].

76. Lillycrop KA, Phillips ES, Torrens C, Hanson MA, Jackson AA, et al. (2008) Feeding pregnant rats a protein-restricted diet persistently alters the methylation of specific cytosines in the hepatic PPAR alpha promoter of the offspring. Br J Nutr 100: 278–282 S0007114507894438 [pii];10.1017/S0007114507894438 [doi].

77. Lillycrop KA, Phillips ES, Jackson AA, Hanson MA, Burdge GC (2005) Dietary protein restriction of pregnant rats induces and folic acid supplementation prevents epigenetic modification of hepatic gene expression in the offspring. J Nutr 135: 1382–1386 135/6/1382 [pii].

78. Stokkan KA, Yamazaki S, Tei H, Sakaki Y, Menaker M (2001) Entrainment of the circadian clock in the liver by feeding. Science 291: 490–493 10.1126/science.291.5503.490 [doi];291/5503/490 [pii].

79. Damiola F, Le MN, Preitner N, Kornmann B, Fleury-Olela F, et al. (2000) Restricted feeding uncouples circadian oscillators in peripheral tissues from the central pacemaker in the suprachiasmatic nucleus. Genes Dev 14: 2950–2961.

80. Schroeder AM, Truong D, Loh DH, Jordan MC, Roos KP, et al. (2012) Voluntary scheduled exercise alters diurnal rhythms of behaviour, physiology and gene expression in wild-type and vasoactive intestinal peptide-deficient mice. J Physiol 590: 6213–6226 jphysiol.2012.233676 [pii];10.1113/jphysiol.2012.233676 [doi].

81. Wolff G, Esser KA (2012) Scheduled exercise phase shifts the circadian clock in skeletal muscle. Med Sci Sports Exerc 44: 1663–1670 10.1249/MSS.0b013e318255cf4c [doi].

82. Suter MA, Chen A, Burdine MS, Choudhury M, Harris RA, et al. (2012) A maternal high-fat diet modulates fetal SIRT1 histone and protein deacetylase activity in nonhuman primates. FASEB J 26: 5106–5114 fj.12-212878 [pii];10.1096/fj.12-212878 [doi].