Non-Alcoholic Fatty Pancreas Disease Pathogenesis: A Role for Developmental Programming and Altered Circadian Rhythms

Abstract

Objectives

Emerging evidence suggests that maternal obesity (MO) predisposes offspring to obesity and the recently described non-alcoholic fatty pancreas disease (NAFPD) but involved mechanisms remain unclear. Using a pathophysiologically relevant murine model, we here investigated a role for the biological clock - molecular core circadian genes (CCG) in the generation of NAFPD.

Design

Female C57BL6 mice were fed an obesogenic diet (OD) or standard chow (SC) for 6 weeks, prior to pregnancy and throughout gestation and lactation: resulting offspring were subsequently weaned onto either OD (Ob_Ob and Con_Ob) or standard chow (Ob_Con and Con_Con) for 6 months. Biochemical, pro-inflammatory and pro-fibrogenic markers associated with NAFPD were then evaluated and CCG mRNA expression in the pancreas determined.

Results

Offspring of obese dams weaned on to OD (Ob_Ob) had significantly increased (p≤0.05): bodyweight, pancreatic triglycerides, macrovesicular pancreatic fatty-infiltration, and pancreatic mRNA expression of TNF-α, IL-6, α-SMA, TGF-β and increased collagen compared to offspring of control dams weaned on to control chow (Con_Con). Analyses of CCG expression demonstrated a phase shift in CLOCK (−4.818, p<0.01), REV-ERB-α (−1.4,p<0.05) and Per2 (3.27,p<0.05) in association with decreased amplitude in BMAL-1 (−0.914,p<0.05) and PER2 (1.18,p<0.005) in Ob_Ob compared to Con_Con. 2-way ANOVA revealed significant interaction between MO and post-weaning OD in expression of CLOCK (p<0.005), PER1 (p<0.005) and PER2 (p<0.05) whilst MO alone influenced the observed rhythmic variance in expression of all 5 measured CCG.

Conclusions

Fetal and neonatal exposure to a maternal obesogenic environment interacts with a post-natal hyper-calorific environment to induce offspring NAFPD through mechanisms involving perturbations in CCG expression.

About the Authors:

Rebeca Carter

Affiliation: Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom

Angelina Mouralidarane

Affiliation: Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom

Junpei Soeda

Affiliation: Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom

Shuvra Ray

Affiliations Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom, Department of Gastroenterology, Guy's and St Thomas' NHS Foundation Trust, London, United Kingdom

Joaquim Pombo

Affiliation: Division of Women's Health, King's College London and King's Health Partners, London, United Kingdom

Ruma Saraswati

Affiliations Division of Women's Health, King's College London and King's Health Partners, London, United Kingdom, Department of Pathology, University College London, London, United Kingdom

Marco Novelli

Affiliation: Department of Pathology, University College London, London, United Kingdom

Giuseppe Fusai

Affiliation: Hepatobiliary and Liver Transplant Unit, Royal Free Hospital, University College London, London, United Kingdom

Francesca Rappa

Affiliations Department of Experimental Biomedicine and Clinical Neurosciences, Section of Human Anatomy, University of Palermo, Palermo, Italy, Euro-Mediterranean Institute of Science and Technology (IEMEST), Palermo, Italy

Chiara Saracino

Affiliation: Department of Medical Sciences, Gastroenterology Unit, IRCCS “Casa Sollievo della Sofferenza” Hospital, San Giovanni Rotondo, Italy

Valerio Pazienza

Affiliation: Department of Medical Sciences, Gastroenterology Unit, IRCCS “Casa Sollievo della Sofferenza” Hospital, San Giovanni Rotondo, Italy

Lucilla Poston

Affiliation: Division of Women's Health, King's College London and King's Health Partners, London, United Kingdom

Paul D. Taylor

Affiliation: Division of Women's Health, King's College London and King's Health Partners, London, United Kingdom

Manlio Vinciguerra

Affiliations Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom, Euro-Mediterranean Institute of Science and Technology (IEMEST), Palermo, Italy, Department of Medical Sciences, Gastroenterology Unit, IRCCS “Casa Sollievo della Sofferenza” Hospital, San Giovanni Rotondo, Italy

Jude A. Oben

* E-mail: [email protected]

Affiliations Institute for Liver and Digestive Health, Royal Free Hospital, University College London, London, United Kingdom, Department of Gastroenterology, Guy's and St Thomas' NHS Foundation Trust, London, United Kingdom

Introduction

Non-alcoholic fatty pancreas disease (NAFPD) is a recently described disease entity associated with an obese and/or dysmetabolic phenotype [1], [2]. NAFPD describes a phenotype ranging from deposition of fat in the pancreas to pancreatic inflammation, and resultant fibrosis. This pancreatic phenotype is similar to that of obesity-induced liver disease, non-alcoholic fatty liver disease (NAFLD), which describes a spectrum from hepatic steatosis through steatohepatitis to cirrhosis, and possible hepatocellular carcinoma [2], [3]. Given that the liver and pancreas have similar embryological origins, it is also plausible, as suggested for NAFLD, that obesity may lead to pancreatic cancer through pancreatic steatosis [4], [5]. This hypothesis is corroborated by studies implicating NAFPD as a risk factor in pancreatic adenocarcinoma [6]–[8].

The prevalence of maternal obesity is increasing worldwide in parallel with adult obesity rates [9], and observational studies suggest an association between maternal obesity, and risk of childhood obesity [10], [11], for which rates are similarly rising with 30% of US and UK children, aged 2 to 15, now classed as overweight or obese [12], [13]. We have previously shown, in a rodent model, that diet induced maternal obesity can play a causative role in the development of NAFPD and that this is exacerbated if the offspring themselves are reared on the same obesogenic diet [1].

The concept of developmental programming suggests that the early environment from conception through to the early post-natal period can alter gene expression through epigenetic processes in the developing offspring, resulting in a permanent alteration in offspring physiology [14]–[16]. Nutrition is considered a ‘major intrauterine environmental factor that alters expression of the fetal genome’ [17], [18], but mechanistic pathways are unclear.

Circadian clocks are molecular oscillators, which drive daily rhythms of physiology and behaviour [19]. The molecular machinery encoding the biological ‘clock’ involves a transcriptional/translational negative feedback loop. The heterodimer CLOCK (circadian locomotor output kaput cycles) and BMAL1 (brain and muscle anrt-like 1) heterodimer complex is controlled, through a negative feedback loop involving Period and Cryptochrome genes [19]. Regulatory accessory pathways include REV-ERB-α, which modulates the circadian clock through BMAL1 expression [19], [20].

Homeostasis is achieved through interactions between the pace-setting hypothalamic suprachiasmatic nucleus (SCN), also known as the master clock, and the peripheral clocks, located in all body cells. The master clock acts as a pacesetter for all peripheral clocks and is predominantly entrained by light [21]. Peripheral clocks however, have been shown to be capable of acting autonomously, and are sensitive to changes in nutritional status [22]. An association between disruption of CCG and metabolism has been shown in several models [19], [21] and, specifically, is implicated in the regulation of metabolic processes executed by the pancreas, including islet cell growth and development [23].

Mice with CLOCK gene mutations are reported to be hyperphagic and obese [24], and a high fat diet has been reported to alter the expression and rhythmicity of CCG in rodents [25]. Recently, isolated reports have suggested that genes contributing to circadian ‘clocks’ may be vulnerable to modulation by the nutritional environment in early life [26], [27]. We therefore hypothesised that maternal over-nutrition may programme offspring dysmetabolism via an altered expression of CCG, resulting in a permanent disruption of circadian rhythms in the pancreas during development.

Specifically, we interrogated the potential for mechanistic involvement of CCG in the pathogenesis of NAFPD arising from an interaction between maternal obesity and a post-weaning nutritional status.

Materials and Methods

All studies were approved by Local University College London Ethics Committee, and conducted under UK Home Office, Animal in Science Regulation Unit (Scientific Procedures) Act 1986 guidelines. Female C57BL/6J mice (n = 60), proven breeders (one previous litter) and approximately 100-day-old (Charles River Laboratories, UK) were maintained under controlled conditions (22°C, 12-hr light/dark cycle) and fed either a standard chow diet (RM1: 15% crude protein, 50% polysaccharides, 7% simple sugars, 5.3% fibre, 3% lipids, 0.2% methionine, 3.3 kcal/g) ad libitum. They were then randomly allocated to either a control standard chow (RM1, Special Dietary Services, UK) or a semi-synthetic energy-rich and highly palatable obesogenic diet (10% simple sugars, 20% animal lard 28% polysaccharide, 23% protein (w/w), Special Dietary Services, UK, energy 4.5 kcal/g). The pelleted obesogenic diet was supplemented with ad libitum access to sweetened condensed milk (approximately 55% simple sugar, 8% fat, 8% protein, w/w, Nestle, SZ) with added micronutrient mineral mix (AIN93G, Special Dietary Services, UK). Combined intake calculated from measured daily intake of pellets and milk (approximately 16% fat, 33% simple sugars, 15% protein, energy 4.0 kcal/g) [28]. After 6–8 weeks females on the obesogenic diet, achieved a 30% increase in body weight and were then mated with C57BL/6J males (day 0 pregnancy signified by the appearance of a copulation plug). Dams were maintained on the obesogenic or control diet throughout gestation and suckling. Following spontaneous delivery, dams and their litters were left undisturbed for 48 hours when litters were standardized to 6 pups with an equal number of males and females (3♂, 3♀), wherever possible, to standardise milk supply. At 3 weeks, offspring were weaned onto either the control (RM1, standard chow, SC, n = 4–5/group) or obesogenic diet (OD) (n = 4–5/group) until 6 months of age. Offspring had ad libitum access to food and water and were maintained in a 12-hour light/dark cycle in a thermostatically controlled environment (22°C). Each of the four groups included pups randomly selected from litters born to different dams.

Offspring Husbandry

The allocation of maternal and offspring post-weaning diet provided 4 groups:

1. maternal diet of SC followed by a post-weaning diet of SC (Con_Con, n = 4–5).

2. maternal diet of SC followed by a post-weaning OD (Con_Ob, n = 4–5).

3. an obesogenic maternal diet followed by a post-weaning diet of SC (Ob_Con, n = 4–5).

4. an obesogenic maternal diet followed by a post-weaning obesogenic diet (Ob_Ob, n = 4–5).

Glucose tolerance test

An oral glucose tolerance test (OGTT) was performed at 6 months, as previously described, with some modification [29]. Briefly mice were fasted for 5 h and D-glucose at 1.5 g/kg was orally administered. OGTT test was performed using different sets of mice, kept under the same experimental conditions as in the main experiments, to avoid the possible influence of fasting on circadian analyses [29].

Tissue collection

At 6 months of age, offspring being maintained on a 12 h light/12 h dark cycle with food and water available ad libitum, were sacrificed at 4 hourly intervals over a 24-hour period, to allow analyses of CCG expression over a 24-hour period (Zeitgeber Time (ZT) 0, 4, 8, 12, 16, and 20, where ZT0 = light on and ZT12 = light off). Following sacrifice, blood samples were taken and harvested organs appropriately stored until analysed. For sample collection in the dark period, mice were transferred to a lit room and terminated within a few minutes of transfer.

Fibrogenic markers and circadian genes

Pancreatic tissue mRNA was assayed by quantitative 2-step PCR for expression of pro-inflammatory markers (Interleukin-6 (IL-6) and tumour necrosis factor-α (TNFα), pro-fibrogenic markers, (collagen type 1-α2; α-smooth muscle actin, αSMA and transforming growth factor–β, TGF-β1) and CCG (Clock, BMAL-1, Period 1, Period 2 and REVERB- α) using pups from all time points. Primer sequences were as shown in Table 1. GAPDH was used as reference gene. The reference were determined by the software program geNORM (Primerdesign, UK).

[Figure omitted. See PDF.]

Table 1. Primer sequences were obtained for the relevant genes of interest from GenBank (except for the primers for the circadian study arm).

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

Western Blotting

Pancreata from the four experimental mice groups were lysed and processed for immunoblotting analysis with specific CLOCK and BMAL1 antibodies as previously described [30], [31]. We pooled together (n = 4) and analyzed equal amounts of proteins lysates of pancreata for each group/time point, as previously reported [32].

Tissue Triglyceride Content

An adaptation of the Folch Method and triglyceride assay reagents (Roche Diagnostics) [33] was used to determine murine whole pancreas tissue triglyceride content.

Histology

Following sacrifice, sections of pancreas were fixed in formalin with the spleen attached to aid orientation in histological analysis. Sections were embedded in paraffin and stained with both H&E and Masson's Trichrome to determine adipocyte infiltration and pericellular fibrosis respectively. Fat infiltration and extent of fibrosis was graded as previously described [34].

Statistical Analysis

Data are expressed as mean ± standard error of the mean (SEM). Means of each group were compared using both one-way and two-way ANOVA, as indicated. Statistical significance was assumed as p<0.05. Cosinor analysis was employed, in addition to ANOVA, to determine rhythmicity of circadian gene expression within a 24-hour period. Cosinor analysis evaluates the ‘mesor’ (circadian rhythm adjusted mean, based on the parameters of a cosine function), timing of the oscillatory crest and amplitude, p<0.05 regarded as significant. Two-way ANOVA was also used to determine the effect of gender and nutritional group on the observed variance between offspring and the relative influence of maternal verses postnatal diet on offspring phenotype. Since no main effect of gender was observed for the reported biomarkers (Table 2) male and female data were combined in the cosinor analysis. The statistical unit ‘n’ used throughout the analysis is number of dams, not number of pups. Statistical analysis was conducted using GraphPad Prism 5 and Stata version 11.2 (StataCorp, College Station, Texas).

[Figure omitted. See PDF.]

Table 2. Two-Way ANOVA results when comparing the influence of gender and group on the overall variance seen between offspring.

https://doi.org/10.1371/journal.pone.0089505.t002

Results

Dysmetabolic Phenotype and NAFPD

Body Weights.

As described in detail above, the allocation of maternal and offspring post-weaning diets provided 4 groups: Con_Con; Con_Ob; Ob_Con; Ob_Ob. Offspring exposed to a post-weaning OD (Con_Ob and Ob_Ob) were significantly heavier than control (Con_Con; p<0.001; p<0.001). Two-way ANOVA further revealed 83% of the variance seen between the groups to be attributable to the post-weaning diet (p<0.0001), with no apparent main effect of MO on body weight, and no interaction between MO and post-weaning diet attributed to variance (Figure 1A).

[Figure omitted. See PDF.]

Figure 1. Anthropometric Data: Con_Con, Con_Ob, Ob_Con, Ob_Ob.

A = Offspring Body Weights; B = Offspring Pancreas Weights; C = Pancreatic Triglyceride Concentrations; D = Macrovesicular Fat Infiltration; E = Oral Glucose Tolerance Test. * = p<.05; ** = p<.001; *** = p<.0001 (n = 4–5/group). ANOVA with Tukey post hoc test.

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

Pancreas Weights.

At 6 months, only offspring exposed to both interventions (Ob_Ob) were found to have significantly heavier pancreas weights compared to control (Con_Con, p<0.01). Moreover, these offspring also displayed heavier pancreata compared to offspring exposed to MO alone (Ob_Con, p<0.01). Two-way ANOVA attributed the post-weaning diet to 64.18% of the observed variance between the groups (p<0.001), with no statistical main effect of MO, nor an interaction between the variables as responsible for the variance (Figure 1B).

Triglyceride Content.

As for pancreatic triglyceride content, this was significantly greater in Ob_Ob offspring than both Ob_Con offspring (p<0.001) and Con_Con offspring (p<0.001). There was no significant difference between offspring exposed to a post-weaning OD alone (Con_Ob) and control. Two-way ANOVA attributed 54.38% of the observed variance to the post-weaning diet (p<0.01) but no statistical power was attributed to the effect of MO (Figure 1C).

Pancreatic Macrovesicular Fat Infiltration was significantly greater in Ob_Ob compared to all other groups: Ob_Con (p<0.05), Con_Ob (p<0.05) and Con_Con (p<0.05) (Figure 1D, n = 4–5). Two-way ANOVA attributed 23% (p<0.05) of observed variance to a post-weaning diet, with no apparent main effect of MO. However, 24.61% (p<0.05) of observed variance was attributed to an interaction between MO and post-weaning OD on two-way ANOVA. (Figure 1D).

Oral Glucose Tolerance Test.

The influence of maternal obesity on offspring glucose tolerance was evident following an oral glucose tolerance test (OGTT), with both Ob_Con offspring (p<0.05) and Ob_Ob offspring (p<0.0001) displaying a significant increase in AUC compared with Con_Con offspring (Figure 1E and Figure S1). Moreover, the OB_Ob AUC was significantly greater than the response of Con_Ob (p<0.05, Figure 1E), indicating the influence of maternal obesity on offspring glucose metabolism. An influence of the postnatal diet was also apparent with Con_Ob demonstrating a significant increase in AUC compared with Con_Con offspring (p<0.05).

Histology.

Histological analysis (H&E staining) of pancreata revealed a marked increase in macrovesicular adipocyte infiltration in offspring exposed to both obesogenic interventions (i.e. Ob_Ob) (Figure 2A, Ob_Ob compared to the other groups), confirming the biochemical results above (Figure 1). Upon scoring for intra-lobular fat Ob_Ob offspring displayed higher values (0.73±0.02), compared to Con_Con (0.02±0.003), Con_Ob (0.2±0.01), or Ob_Con (0.3±0.02) [33]. Masson's Trichrome staining for pericellular fibrosis highlighted areas of occasional increased fibrosis in Ob_Ob offspring (Figure 2B, Ob_Ob compared to the other panels) but on scoring fibrosis grade was not significantly different between the groups.

[Figure omitted. See PDF.]

Figure 2. Pancreatic Histology: Con_Con = Top Left; Ob_Ob = Top Right; Ob_Con = Bottom Left; Ob_Ob = Bottom Right.

A = Haematoxylin and Eosin Stain for Pancreas Histology (×40); B = Masson's Trichrome Stain for Pancreas Histology (×40) (n = 4–5/group).

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

Pancreatic Inflammation and fibrosis

We then sought to analyze the gene expression levels of the major inflammatory cytokines (IL-6, TNF-α) and the fibrotic indicators (α-SMA, collagen, TGF-β) in the pancreata of the 4 mice groups, to test the potential effects of MO and post-weaning diets on these pathogenic processes. All the marker tested in PCR showed same trend across the gropus in all time points, therefore, ZT 8 was chosen for further analysis in all groups.

IL-6.

IL-6 mRNA expression in Ob_Ob offspring was higher when compared to all other groups: versus Ob_Con group (p<0.05), versus Con_Ob group (p<0.05), and versus Con_Con (p<0.001). There was no significant difference between offspring exposed to a post-weaning OD alone (Con_Ob) or MO alone (Ob_Con) compared to controls (Con_Con), indicating the need for the combined presence of both interventions to induce upregulation of IL-6 expression. Two-Way ANOVA further corroborated these results, revealing MO and the post-weaning OD to be attributable to 43.71% (p<0.001) and 34.53% (p<0.01) of the overall variance in IL-6 expression respectively, but also, that there was a significant interaction between these variables, contributing 10.18% (p<0.05) of the overall variance (Figure 3A).

[Figure omitted. See PDF.]

Figure 3. mRNA Expression of markers of Pancreatic Injury: Con_Con, Con_Ob, Ob_Con, Ob_Ob.

A = IL-6; B = TNF-α; C = α-SMA; D = Collagen; E = TGF- β. * = p<.05; ** = p<.001; *** = p<.0001 (n = 4–5/group). ANOVA with Tukey post hoc test.

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

TNF-α.

The influence of the interventions on IL-6 was mirrored in TNF-α, another marker of pancreatic inflammation, with Ob_Ob offspring demonstrating a significantly increased expression of this cytokine compared to the other offspring groups: versus Ob_Con (p<0.001), versus Con_Ob (p<0.0001) and versus Con_Con (p<0.0001). In addition, MO was found to have an independent effect with Ob_Con offspring exhibiting significant upregulation of TNF-α compared with control (p<0.05). Two-way ANOVA attributed 61.22% and 18.35% of the overall variance to MO (p<0.0001) and the post-weaning diet to (p<0.001), respectively. An overall significant interaction between the variables attributed 4.63% (p<0.05) of the total variance. As with IL-6, there was no significant upregulation of TNF-α by a post-weaning OD in isolation, compared with control (Figure 3B).

Markers of Pancreatic Fibrosis

α-SMA.

The influence of MO on markers of offspring pancreatic fibrotic injury was evident in α-SMA expression, with Ob_Ob displaying upregulation compared to the other offspring groups: versus Ob_Con (p<0.001); versus Con_Ob (p<0.001) and versus Con_Con (p<0.0001). Two-way ANOVA further supported these results, attributing 63.77% (p<0.0001) and 13.08% (p<0.01) of the observed variance to MO and post-weaning OD respectively, and 5.53% (p<0.05) of the variance to an overall interaction between the variables. No significant difference was noted between offspring exposed to a post-weaning OD alone (Con_Ob) and control (Figure 3C).

Collagen.

Similarly, Ob_Ob collagen mRNA expression was significantly greater than the other 3 offspring groups: versus Ob_Con (p<0.001), versus Con_Ob (p<0.0001), and versus Con_Con (p<0.0001). Two-way ANOVA attributed 61.26% (p<0.0001) and 28.39% (p<0.001) of observed variance to MO and post-weaning OD respectively and 7.14% (p<0.05) of the variance to an overall interaction between the variables. Again, no significant difference was observed between Con_Oband Con_Con (Figure 3D).

TGF-β.

TGF-β mRNA expression was significantly greater in Ob_Ob compared with Ob_Con (p<0.05) and Con_Con (p<0.001) indicating an additive effect of the postweaning obesogenic diet. No significant difference was noted between Ob_Ob and Con_Ob, nor was a significant difference noted between Con_Ob and control. Two-way ANOVA revealed MO contributing 64% (p = 0.0009) of the observed variance, indicating MO as having the main effect on the observed variance in results (Figure 3E).

Core Circadian Genes

Results are divided into ANOVA and cosinor analysis of CCG of offspring over 6 separate time points (Figure 4):

[Figure omitted. See PDF.]

Figure 4. mRNA Expression of Circadian Genes: Con_Con, Con_Ob, Ob_Con, Ob_Ob.

A = Clock; B = Per-2; C = BMAL-1; D = Per-1; E = REV-ERB-α. %%% = p<.0001; ∞∞∞ = p<.0001; * = p<.05; ** = p<.001; *** = p<.0001; ## = p<.001; § = p<.05. (% = Con_Con versus Con_Ob; ∞ = Con_Con versus Ob_Con; * = Con_Con versus Ob_Ob; # = Con_Ob versus Ob_Ob; § = Ob_Con versus Ob_Ob) (n = 4–5/group). ANOVA with Tukey post hoc test.

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

ANOVA Analysis of Circadian Gene mRNA Expression

Both CLOCK and Per2 gene expression displayed differences as a result of MO and a post-weaning OD. Analysis of CLOCK gene expression revealed a significant difference between Con_Con expression and the 3 other groups at ZT4: vs Ob_Ob (p<0.001), vs Ob_Con (p<0.0001) and vs Con_Ob (p<0.05). Two-way ANOVA revealed that MO accounted for 64% of the total variance seen and an interaction between the variables as accountable for 23% (p<0.05) (Figure 4A). Per-2 gene expression revealed a different pattern, with a decreased expression in Ob_Ob compared with the three other groups at ZT8: vs Ob_Con (p<0.05), vs Con_Ob (p<0.001) and vs Con_Con (p<0.05). Two-way ANOVA further revealed an overall interaction as being attributable to 23% of total variance (p<0.05), with MO attributing a further 39% (p<0.05), similarly indicating that MO had the main effect on the observed variance (Figure 4B).

At ZT16, a perturbation was also observed for BMAL-1 as noted for Per-2 expression at ZT8 (Ob_Ob vs all other groups, p<0.05). Two-way ANOVA revealed that MO and an overall interaction between MO and a post-weaning OD accounted for 18%, and 54% variance respectively, p<0.05. Additionally, the post-weaning diet contributed a further 44% (p<0.05) of the overall variance (Figure 4C).

Further interactions between the sequential circadian genes were observed with the same patterns of gene expression noted for CLOCK at ZT4 and Per-1 at ZT16. As observed with CLOCK at ZT4, Per-1 gene expression was different between Con_Con and Con_Ob (p<0.0001), Ob_Con (p<0.0001) and Ob_Ob (p<0.0001). Two-way ANOVA revealed an overall interaction between MO and OD as contributing 10% of total variance (p<0.05), MO alone contributing 49% (p<0.0001) and the post-weaning OD 24% (p<0.001), (Figure 4D). Two-way ANOVA revealed that MO attributed to 31% (p<0.05) of the variance in gene expression.

Protein levels of core clock genes CLOCK and BMAL1

To determine if variations observed at the mRNA level were mirrored at the protein level, we now analyzed protein expression levels of the core components of the circadian clock machinery, CLOCK and BMAL1, in the pancreata of the 4 experimental groups at ZT 0 h, 8 h and 16 h by immunoblotting (Figure S2). Strikingly higher protein levels of CLOCK and BMAL1 were observed at ZT16, with a peak of BMAL1 expression present in the Ob_Ob offspring, while an increased expression of CLOCK was detected in the Con_Ob group. Highest levels of CLOCK and BMAL1 proteins at ZT16 are consistent with observed mRNA expression (Figure 4A and C).

Cosinor Analysis of mRNA Expression

Clock.

Cosinor analysis revealed a phase shift of −4.818005 Hours (p<0.01) when comparing Ob_Ob with Ob_Con, suggesting a main effect of the postweaning diet. No other significant differences were seen between the groups in relation to CLOCK gene expression (Figure 5A; Table 3).

[Figure omitted. See PDF.]

Figure 5. Cosinor Analysis of Circadian Genes: Con_Con, Con_Ob, Ob_Con, Ob_Ob.

A = Clock; B = Bmal-1; C = Per-1; D = Per-2; E = REV-ERB-α. (% = Con_Con versus Con_Ob; * = Con_Con versus Ob_Ob; § = Ob_Con versus Ob_Ob) (n = 4–5/group).

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

[Figure omitted. See PDF.]

Table 3. Results of Cosinor Analysis of Circadian Gene mRNA Expression.

https://doi.org/10.1371/journal.pone.0089505.t003

BMAL-1.

Cosinor analysis revealed significantly reduced amplitude, with respect to Con_Con, in Ob_Ob offspring (−0.914, p<0.05) and Con_Ob offspring (−0.972, p<0.05), suggesting main effect of postweaning diet. Furthermore, maxima were calculated to be significantly reduced in Ob_Ob offspring when compared with Con_Con (−1.86 p<0.05). No other significant differences were seen between the groups in relation to BMAL-1 gene expression (Figure 5B; Table 3).

Period 1.

Cosinor analysis yielded no statistically significant difference between the groups in Per 1 (Figure 5C; Table 3).

Period 2.

There was a phase shift of 3.27 Hours (p<0.05) between Ob_Ob and Ob_Con offspring suggesting main effect of MO. Amplitude was also calculated as increased in Ob_Ob (1.03, p<0.05) and Ob_Con (1.18, p<0.005), when compared with Con_Con suggesting evidence of both MO and an interaction between the variables as having a significant effect. Furthermore, maxima were increased in Ob_Ob compared with Con_Con (1.68, p<0.05), and minima decreased Ob_Con compared with Con_Con (−1.11, p<0.05). No other significant differences were seen between the groups in relation to Per-2 gene expression (Figure 5D; Table 3).

REVERB-α.

There was a phase shift of −1.4 Hours (p<0.05) when comparing offspring exposed to a post-weaning OD (Con_Ob, Ob_Ob) with COn_Con offspring. No other significant differences were seen between the groups in relation to REVERB-α gene expression (Figure 5E; Table 3).

Discussion

NAFPD is an emerging clinical entity associated with dysmetabolism and characterized by pancreatic fat deposition, inflammation, fibrosis and pancreatitis. This pancreatic phenotype is strikingly similar to that of obesity-induced NAFLD, which describes a spectrum characterized by hepatic steatosis, steatohepatitis and cirrhosis. The similarity in phenotypes could be due to the common embryonic origins of the liver and pancreas. We have previously shown that offspring exposure to MO throughout pregnancy and lactation induces a significant increase in markers indicative of a NAFPD phenotype, when compared to offspring exposed to a normal intrauterine and perinatal environment [1]. Here, we corroborate these results using a more pathophysiologically relevant model not involving cross-fostering which has since been shown to influence metabolic phenotype [35]. In addition to the previous findings of increased pancreatic triglycerides, collagen and TGF-β expression [1], our results show that MO in conjunction with a post-weaning OD (Ob_Ob) significantly impacts upon pancreas weight, pancreatic triglycerides concentration and macrovesicular fat concentration.

Moreover, offspring glucose tolerance was impaired not only in Ob_Ob but also in offspring exposed to MO in isolation (Ob_Con), accentuating an independent effect of MO. Thus, MO appears to have a deleterious effect on offspring phenotype, particularly when coupled with a post-weaning OD, which in isolation only significantly influenced offspring bodyweight. Since the offspring of obese mothers, weaned onto an OD (Ob_Ob) demonstrated the highest pancreatic triglyceride (TG) concentration, with no change in TGs in Con_Ob, this was most likely due to an interaction between the two interventions (maternal and post weaning OD diet) rather than the post-weaning diet alone (Figure 1C). The TG profile was mirrored in the infiltration of macrovesicular fat.

However, as Ob_Ob fat infiltration was higher than Con_Ob, this exposed an independent influence of MO in these offspring (Figure 1D). Thus, MO appears to have a priming or ‘first hit’ effect, which only become evident when offspring are subsequently challenged by a hypercalorific diet, exacerbating disease phenotype as a result. OGTT observations indicated impaired glucose tolerance in offspring of obese dams, which builds on a recent report [36], demonstrating that diet-induced obesity in female dams, whose offspring were weaned onto standard chow, led to insulin resistance at 8–11 months. Here we show this occurs at 6 months, and that it is exacerbated by a post-weaning obesogenic diet (Figure 1E).

Moreover, we found that control offspring exposed to a post-weaning OD showed less evidence of inflammation than offspring of the obese dams, on the same diet adding to evidence for a priming effect of MO. The higher TNF-α expression in the MO offspring, strengthens a role for MO in the programming of offspring disease. This is not the first study to demonstrate a role for MO in activation of inflammatory cytokines in offspring [1], [37], but is the first to demonstrate it in pancreatic tissue, in the absence of cross-fostering. Similar patterns were observed in markers of pancreatic fibrosis, further supporting the priming effect of MO. To summarize, these phenotypic results consistently highlight the importance of the maternal environment on offspring phenotype. This is emphasised by the two-way ANOVA analysis, which consistently implicate MO in the variance amongst the markers of pancreatic injury.

Support for the concept of programming of clock genes by nutrition is provided by a recent report in which adult offspring of protein-restricted rat dams demonstrated permanently altered expression in a functional network of hypothalamic nuclear receptors and co-regulators of the circadian clock involved in lipid metabolism [38]. Moreover, rats exposed to chronic phase shifts of the photoperiod during gestation showed altered energy balance with lifelong consequences for the metabolic homeostasis in the offspring which develop increased obesity and abnormal glucose metabolism [39].

Our results further report a significant disruption in CCG, which could have a mechanistic role for CCG in developmentally programmed NAFPD. Moreover, analysis of CCG expression at different stages of daylight revealed corresponding genes being affected in the same pattern but at opposite times of day (CLOCK and Per-2, BMAL-1 and Per-1; Figure 4A–D). We were able to confirm high protein expression levels of core clock component CLOCK and BMAL1 at ZT16, consistent with their mRNA levels (Figure S2). In addition, two-way ANOVA demonstrated a consistent role both for MO, and an interaction between MO and a post-weaning OD, in the observed variance in CCG. Cosinor analysis build on these observations, describing an array of perturbations in the normal diurnal rhythms, with differences predominantly found between Ob_Ob offspring and Con_Con offspring.

Overall, this evidence suggests that MO, in conjunction with a post-weaning OD, results in a significant disruption of the biological diurnal rhythms. Moreover, as no significant difference is seen between Con_Ob and Con_Con in isolation, in any of the CCG, further credence is given to role of MO in the perturbation of diurnal rhythms. As disruption of CCG has been implicated in altered regulation of metabolic pancreatic processes [40] and in hyperphagia and obesity [24], and the circadian rhythms are programmable [38], it is plausible that expression of CCG could be irreversibly altered following exposure to MO in utero. Thus, perturbation of CCG may be involved mechanistically in the priming effects of MO through epigenetic mechanisms, when the offspring is challenged post weaning with an OD.

In conclusion, our results, in a novel and pathophysiologically relevant NAFPD model, consistently show that offspring exposed to both MO and a post-weaning OD, have a more pronounced dysmetabolic and NAFPD phenotype than control since offspring exposed to a post-weaning OD in isolation were rarely significantly different to control and suggests, that MO has a ‘priming’ effect on offspring phenotype which in the presence of a post-weaning OD results in an exacerbated dysmetabolic offspring phenotype. These findings support the proposal that MO, through developmental programming, may be exacerbating the obesity epidemic and its associated metabolic disorders such as NAFPD [41]. Mechanistically, since CCG are entrainable by nutritional stimuli [19], [21], [41], [42], and in light of the significant differences between Con_Con and Ob_Con/Ob_Ob on cosinor and ANOVA analyses, we propose that the intra-uterine period may through perturbation of CCG be involved in programmed NAFPD.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Oral glucose tolerance tests at 6 months. Blood glucose concentrations, after oral glucose administration (1.5 g/kg), are represented. Ob_Ob showed significantly higher glucose levels at 15 minutes than Con_Con and Con_Ob (p<0.001 and p<0.05, respectively). **p<0.001 vs. Con_Ob, # p<0.05 vs. Con_Con.

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

(TIF)

Figure S2.

Upper panel: pancreatic tissues from 4 different animals per group/time points were pooled together as previously described [32], lysed and equal amounts of proteins (45 µg) were loaded on a 10% polyacrylamide gel, separated by electrophoresis and immunoblotted with specific CLOCK and BMAL1 primary antibodies. β-actin expression served as loading control. Lower panel: densitometric quantification of BMAL1 and CLOCK proteins normalized to β-actin expression.

https://doi.org/10.1371/journal.pone.0089505.s002

(TIF)

Acknowledgments

We gratefully acknowledge Paul Seed and David Sugden for their help with the statistical analysis.

Author Contributions

Conceived and designed the experiments: RC AM JS SR LP PT JO. Performed the experiments: RC AM JS SR FR CS VP. Analyzed the data: RC AM JS JP RS MV MN GF FR CS VP. Contributed reagents/materials/analysis tools: JP RS MV MN GF. Wrote the paper: RC AM SR MV JO.

Citation: Carter R, Mouralidarane A, Soeda J, Ray S, Pombo J, Saraswati R, et al. (2014) Non-Alcoholic Fatty Pancreas Disease Pathogenesis: A Role for Developmental Programming and Altered Circadian Rhythms. PLoS ONE 9(3): e89505. https://doi.org/10.1371/journal.pone.0089505

References

1. Oben JA, Patel T, Mouralidarane A, Samuelsson AM, Matthews P, et al. (2010) Maternal obesity programmes offspring development of non-alcoholic fatty pancreas disease. Biochem Biophys Res Commun 394: 24–28.

2. Fraulob JC, Ogg-Diamantino R, Fernandes-Santos C, Aguila MB, Mandarim-de-Lacerda CA (2010) A Mouse Model of Metabolic Syndrome: Insulin Resistance, Fatty Liver and Non-Alcoholic Fatty Pancreas Disease (NAFPD) in C57BL/6 Mice Fed a High Fat Diet. J Clin Biochem Nutr 46: 212–223.

3. Smits MM, van Geenen EJ (2011) The clinical significance of pancreatic steatosis. Nat Rev Gastroenterol Hepatol 8: 169–177.

4. Pitt HA (2007) Hepato-pancreato-biliary fat: the good, the bad and the ugly. HPB (Oxford) 9: 92–97.

5. Patel AV, Rodriguez C, Bernstein L, Chao A, Thun MJ, et al. (2005) Obesity, recreational physical activity, and risk of pancreatic cancer in a large U.S. Cohort. Cancer Epidemiol Biomarkers Prev 14: 459–466.

6. Larsson SC, Orsini N, Wolk A (2007) Body mass index and pancreatic cancer risk: A meta-analysis of prospective studies. Int J Cancer 120: 1993–1998.

7. Bracci PM (2012) Obesity and pancreatic cancer: overview of epidemiologic evidence and biologic mechanisms. Mol Carcinog 51: 53–63.

8. Renehan AG, Tyson M, Egger M, Heller RF, Zwahlen M (2008) Body-mass index and incidence of cancer: a systematic review and meta-analysis of prospective observational studies. Lancet 371: 569–578.

9. Organization WH (2012) Obesity and overweight. WHO. Available: http://www.who.int/mediacentre/factsheets/fs311/en/index.html. Accessed 2014 Feb 28.

10. Nelson SM, Matthews P, Poston L (2010) Maternal metabolism and obesity: modifiable determinants of pregnancy outcome. Hum Reprod Update 16: 255–275.

11. Norman JE, Reynolds RM (2011) The consequences of obesity and excess weight gain in pregnancy. Proc Nutr Soc 70: 450–456.

12. Hurt RT, Frazier TH, McClave SA, Kaplan LM (2011) Obesity epidemic: overview, pathophysiology, and the intensive care unit conundrum. JPEN J Parenter Enteral Nutr 35: 4S–13S.

13. 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.

14. Hinde K, Capitanio JP (2010) Lactational programming? Mother's milk energy predicts infant behavior and temperament in rhesus macaques (Macaca mulatta). Am J Primatol 72: 522–529.

15. Alfaradhi MZ, Ozanne SE (2011) Developmental programming in response to maternal overnutrition. Front Genet 2: 27.

16. Peaston AE, Whitelaw E (2006) Epigenetics and phenotypic variation in mammals. Mamm Genome 17: 365–374.

17. Wu G, Bazer FW, Cudd TA, Meininger CJ, Spencer TE (2004) Maternal nutrition and fetal development. J Nutr 134: 2169–2172.

18. Gluckman PD, Hanson MA, Cooper C, Thornburg KL (2008) Effect of in utero and early-life conditions on adult health and disease. N Engl J Med 359: 61–73.

19. Mazzoccoli G, Pazienza V, Vinciguerra M (2012) Clock genes and clock-controlled genes in the regulation of metabolic rhythms. Chronobiol Int 29: 227–251.

20. Cho H, Zhao X, Hatori M, Yu RT, Barish GD, et al. (2012) Regulation of circadian behaviour and metabolism by REV-ERB-alpha and REV-ERB-beta. Nature 485: 123–127.

21. Tevy MF, Giebultowicz J, Pincus Z, Mazzoccoli G, Vinciguerra M (2013) Aging signaling pathways and circadian clock-dependent metabolic derangements. Trends Endocrinol Metab 24: 229–237.

22. Schibler U, Ripperger J, Brown SA (2003) Peripheral circadian oscillators in mammals: time and food. J Biol Rhythms 18: 250–260.

23. Lamia KA, Evans RM (2010) Metabolism: Tick, tock, a beta-cell clock. Nature 466: 571–572.

24. Turek FW, Joshu C, Kohsaka A, Lin E, Ivanova G, et al. (2005) Obesity and metabolic syndrome in circadian Clock mutant mice. Science 308: 1043–1045.

25. Kohsaka A, Laposky AD, Ramsey KM, Estrada C, Joshu C, et al. (2007) High-fat diet disrupts behavioral and molecular circadian rhythms in mice. Cell Metab 6: 414–421.

26. Orozco-Solis R, Matos RJ, Lopes de Souza S, Grit I, Kaeffer B, et al. (2011) Perinatal nutrient restriction induces long-lasting alterations in the circadian expression pattern of genes regulating food intake and energy metabolism. Int J Obes (Lond) 35: 990–1000.

27. 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.

28. 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.

29. Deblon N, Bourgoin L, Veyrat-Durebex C, Peyrou M, Vinciguerra M, et al. (2012) Chronic mTOR inhibition by rapamycin induces muscle insulin resistance despite weight loss in rats. Br J Pharmacol 165: 2325–2340.

30. Benegiamo G, Mazzoccoli G, Cappello F, Rappa F, Scibetta N, et al. (2013) Mutual Antagonism between Circadian Protein Period 2 and Hepatitis C Virus Replication in Hepatocytes. PLoS One 8: e60527.

31. Clement S, Peyrou M, Sanchez-Pareja A, Bourgoin L, Ramadori P, et al. (2011) Down-regulation of phosphatase and tensin homolog by hepatitis C virus core 3a in hepatocytes triggers the formation of large lipid droplets. Hepatology 54: 38–49.

32. Field MD, Maywood ES, O'Brien JA, Weaver DR, Reppert SM, et al. (2000) Analysis of clock proteins in mouse SCN demonstrates phylogenetic divergence of the circadian clockwork and resetting mechanisms. Neuron 25 (2) 437–47.

33. Veyrat-Durebex C, Montet X, Vinciguerra M, Gjinovci A, Meda P, et al. (2009) The Lou/C rat: a model of spontaneous food restriction associated with improved insulin sensitivity and decreased lipid storage in adipose tissue. Am J Physiol Endocrinol Metab 296: E1120–1132.

34. Mathur A, Marine M, Lu D, Swartz-Basile DA, Saxena R, et al. (2007) Nonalcoholic fatty pancreas disease. HPB (Oxford) 9: 312–318.

35. Matthews PA, Samuelsson AM, Seed P, Pombo J, Oben JA, et al. (2011) Fostering in mice induces cardiovascular and metabolic dysfunction in adulthood. J Physiol 589: 3969–3981.

36. Nivoit P, Morens C, Van Assche FA, Jansen E, Poston L, et al. (2009) Established diet-induced obesity in female rats leads to offspring hyperphagia, adiposity and insulin resistance. Diabetologia 52: 1133–1142.

37. Grabiec K, Milewska M, Grzelkowska-Kowalczyk K (2012) [Maternal obesity and the development of skeletal muscle in offspring–fetal origin of metabolic disorders]. Postepy Hig Med Dosw (Online) 66: 1–10.

38. Orozco-Solis R, Matos RJ, Guzman-Quevedo O, Lopes de Souza S, Bihouee A, et al. (2010) Nutritional programming in the rat is linked to long-lasting changes in nutrient sensing and energy homeostasis in the hypothalamus. PLoS One 5: e13537.

39. Varcoe TJ, Boden MJ, Voultsios A, Salkeld MD, Rattanatray L, et al. (2013) Characterisation of the maternal response to chronic phase shifts during gestation in the rat: implications for fetal metabolic programming. PLoS One 8: e53800.

40. Kohsaka A, Bass J (2007) A sense of time: how molecular clocks organize metabolism. Trends Endocrinol Metab 18: 4–11.

41. 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.

42. Vinciguerra M, Tevy MF, Mazzoccoli G (2013) A ticking clock links metabolic pathways and organ systems function in health and disease. Clin Exp Med Apr 4. [Epub ahead of print].