Next Article in Journal
NMR Fragment-Based Screening against Tandem RNA Recognition Motifs of TDP-43
Next Article in Special Issue
Gene–Environment Interactions on Body Fat Distribution
Previous Article in Journal
Adiponectin and Cardiovascular Risk. From Pathophysiology to Clinic: Focus on Children and Adolescents
Previous Article in Special Issue
Reduced Diet-induced Thermogenesis in Apolipoprotein A-IV Deficient Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maternal and Post-Weaning High-Fat Diets Produce Distinct DNA Methylation Patterns in Hepatic Metabolic Pathways within Specific Genomic Contexts

1
Division of Nutritional Sciences, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
2
Department of Human Genetics, University of California at Los Angeles, CA 90095, USA
3
Department of Food Science and Human Nutrition, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
4
Department of Food Science and Human Nutrition, Division of Nutritional Sciences, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
5
Department of Food Science and Human Nutrition, Division of Nutritional Sciences, Illinois Informatics Institute, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(13), 3229; https://doi.org/10.3390/ijms20133229
Submission received: 28 May 2019 / Revised: 24 June 2019 / Accepted: 26 June 2019 / Published: 30 June 2019
(This article belongs to the Special Issue Nutrition Genomics)

Abstract

:
Calorie-dense high-fat diets (HF) are associated with detrimental health outcomes, including obesity, cardiovascular disease, and diabetes. Both pre- and post-natal HF diets have been hypothesized to negatively impact long-term metabolic health via epigenetic mechanisms. To understand how the timing of HF diet intake impacts DNA methylation and metabolism, male Sprague–Dawley rats were exposed to either maternal HF (MHF) or post-weaning HF diet (PHF). At post-natal week 12, PHF rats had similar body weights but greater hepatic lipid accumulation compared to the MHF rats. Genome-wide DNA methylation was evaluated, and analysis revealed 1744 differentially methylation regions (DMRs) between the groups with the majority of the DMR located outside of gene-coding regions. Within differentially methylated genes (DMGs), intragenic DNA methylation closer to the transcription start site was associated with lower gene expression, whereas DNA methylation further downstream was positively correlated with gene expression. The insulin and phosphatidylinositol (PI) signaling pathways were enriched with 25 DMRs that were associated with 20 DMGs, including PI3 kinase (Pi3k), pyruvate kinase (Pklr), and phosphodiesterase 3 (Pde3). Together, these results suggest that the timing of HF diet intake determines DNA methylation and gene expression patterns in hepatic metabolic pathways that target specific genomic contexts.

1. Introduction

High-fat (HF) diet intake during both the pre-natal and post-natal periods has adverse effects on metabolic health, including higher body weight, reduced insulin sensitivity, systemic inflammation, and dysregulation of hypothalamic satiety cues [1,2]. Furthermore, the metabolic consequences of HF diet consumption have been in part attributed to diet-induced epigenetic modifications. For instance, post-natal HF diet exposure in rats and mice has been associated with differential methylation of the hunger/satiety genes melanocortin 4 receptor (Mc4r), proopiomelanocortin (Pomc), and leptin (Lep) as well as the glycotic enzymes glucokinase (Gck) and pyruvate kinase (Pklr) [3,4,5]. Perinatal HF diet also produces methylation differences in the energy homeostasis gene, peroxisome proliferator-activated receptor α (Ppara) [6], the inflammatory genes, toll-like receptors 1 and 2 (Tlr-1 and Tlr-2) [7], and the hepatic cell cycle inhibitor, cyclin-dependent kinase inhibitor (Cdkn1a) [8].
The perinatal and post-natal periods studied in the aforementioned publications are characterized by separate epigenetic events. De novo DNA methylation is established in the perinatal period while environmental factors during the post-natal period can interfere with the maintenance of such methylation patterns [9,10,11]. Given these differences between perinatal and post-natal epigenetic processes, it is important to understand how each exposure window impacts metabolism and gives rise to a unique methylome. As opposed to using a candidate-gene approach to examine the effects of HF diet on DNA methylation, the present study instead used more comprehensive high-throughput sequencing methodologies. Genome-wide DNA methylation can be measured through using one of the methods that differentiates the methylated cytidines from unmethylated ones: bisulfite conversion, affinity capture, or endonuclease digestion. While bisulfite sequencing provides superior coverage and resolution, it is both labor-intensive and comparatively more expensive [12,13]. Alternatively, Methylated DNA Immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq) is an affinity enrichment method that uses 5mC-specific antibodies to directly quantify methylated DNA sequences. MeDIP-seq offers a more cost-effective option, and integration of Methylation-sensitive Restriction Enzyme digestion followed by sequencing (MRE-seq) compensates for the limitations of MeDIP-seq by having single-base resolution and extensive coverage at regions with low CpG densities. MRE-seq is an endonuclease-based technique in which restriction enzymes identify unmethylated regions [14]. Alone, MRE-seq is confined by the sequence-specificity of existing restriction enzymes, but integrating MRE-seq with MeDIP-seq can improve coverage and resolution of genome-wide methylation and proves to be a viable way of assessing methylation in large samples [15].
The liver changes substantially throughout the early stages of life. Not only does the organ undergo structural transformation, but the functional demands before and after birth are also very different. This is reflected in the altered metabolic profiles including insulin output and hepatic enzyme kinetics in the fetus versus adult [16,17]. Therefore, the present study was designed to combine MRE-seq and MeDIP-seq to survey the effects of HF diet exposure on genome-wide methylation to further elucidate the underlying mechanisms and pathways related to hepatic metabolic adaptation. Male rats were exposed to a HF diet during either the perinatal period (maternal HF, MHF) or at post-weaning (post-weaning HF, PHF). To understand how timing of HF diet exposure impacts the hepatic epigenome, we first compared DNA methylation patterns directly between MHF and PHF rats. Using the identified differentially methylated regions (DMRs), we then compared DNA methylation and gene expression between MHF, PHF, and lifelong control-fed animals (CON) to assess whether any HF diet intake affects DNA methylation. Finally, we further examined the metabolic pathways that were potentially differentially influenced by the HF diet at perinatal and post-natal time points.

2. Results

2.1. Compared to Maternal HF Diet Exposure, Post-Weaning HF Does Not Alter Body Weight but Induces Greater Hepatic Lipid Accumulation

To assess the potential differential effects of HF exposure at different developmental stages, male rats were exposed to HF diet either during gestation and lactation (maternal high fat: MHF) or after weaning (post-weaning high fat: PHF) (Figure 1a). Additionally, a control group (CON) was fed an AIN93G diet during both the perinatal and post-weaning periods. Food intake in grams per week was not significantly different between the three groups, except on week 6 and 7, when CON animals consumed more than PHF rats (Figure 1c). However, given the higher caloric density of the HF diet, the PHF group consistently consumed more calories during the post-weaning period than both CON and MHF (Figure 1d). PHF rats were only significantly heavier than CON, but not MHF, on post-natal weeks 5, 7, 10, 11 and 12 (Figure 1b). Food efficiency (weight gain/kcal intake) was lower in PHF compared to both CON and MHF (Figure S1). ORO staining revealed that lipid accumulation was unchanged in MHF compared to CON, but was significantly higher in PHF (Figure 1e,f). To account for differences in body weight and lipid accumulation, we measured expression of de novo lipogenesis and triglyceride synthesis genes and found decreased levels in PHF and increased levels in MHF (Figure S2).

2.2. Timing HF Diet Exposure Determines Distinct Methylation Patterns

We first compared DNA methylation patterns directly between MHF and PHF to identify regions that differ according to timing of HF diet exposure. Genome-wide DNA methylation was measured in liver using MeDIP-seq and MRE-seq. MeDIP-seq produces extensive CpG coverage while MRE-seq offers high resolution such that combining the two methods has been proven to improve detection quality [18]. Sequencing produced over 49 million filtered MeDIP reads and over 31 million filtered MRE reads for each animal. The R-package MethylMnM was used to integrate reads from both methods to provide a comprehensive analysis of differences in methylation profile between groups. We uncovered 1744 DMRs (FDR q-value < 0.05), of which 990 (57%) were more methylated in the PHF group and 754 (43%) were more methylated in the MHF group (Figure 2a,b). These DMRs were distributed fairly evenly across 21 chromosomes (Data not shown).
Methylation results were validated in multiple animals from each of the groups using both a probe-based PCR method and MSP. First, a DMR located within a gene of interest, pyruvate kinase (Pklr), was verified using a fluorescent probe-based PCR technique (Figure 2, Table 1). Using the fluorescent probe, we confirmed that the MeDIP- and MRE-identified DMR was highly methylated in PHF (74% methylated) when compared to MHF (43% methylated, Figure 2c,d). We used the same method to examine a genomic location upstream of the Il6 gene that was not differentially methylated according to MeDIP and MRE analysis. Probe-based analysis of DNA methylation confirmed this result, as the methylation in the MHF (48% methylated) and PHF (47% methylated) groups were not significantly different (Figure 2). Additionally, we used MSP to confirm that a region within the Acacb gene was differentially methylated between PHF and MHF rats (Table 1). PHF animals had higher methylation than MHF animals (70% and 25%, respectively; Figure 2). MSP also confirmed that a region in the Gpam gene was not differentially methylated, as predicted by MeDIP-seq and MRE-seq (PHF: 57%, MHF: 53%). Given the strict criteria set by the MethylMnM package and the results from the subsequent validation both in this study and in previous work [19], we assume that any MeDIP- and MRE-identified DMRs represent true differential methylation between groups.
Next, we examined DMR location relative to CpG islands. Despite the relative enrichment of DMRs in island regions, most of the differential methylation was observed in the 2000 bp shore sequences flanking the islands. In fact, more than twice as many DMRs were identified in shore regions (Figure 3a, n = 484, 27.8%) than in CpG islands (n = 227, 13.0%). This is consistent with previous reports showing that most treatment-induced changes of DNA methylation does not occur within the conserved CpG islands, but rather is localized primarily to nearby regions with slightly lower CpG densities, including the shores [20,21].
DMRs were then aligned with gene positions. Of the 1744 identified DMRs, 661 (38%) were associated with protein-coding genes (Figure 3b). Among the 661 gene-associated DMRs, only 43 were positioned within the 1500 bp upstream promoter region and 61 were in the 1500 bp downstream regions. The rest of the 557 gene-associated DMRs are located within intragenic regions. Furthermore, 37%, 40%, and 43% were more methylated in the MHF group than the PHF group, and 63%, 60%, and 57% were more methylated in the PHF group than the MHF group within the promoter, gene body, and downstream regions, respectively.
A majority of the DMRs identified (n = 1,083, 62%) were not associated with a protein-coding gene. These regions were further examined for non-coding RNA sequences using the Genome Browser rat dataset, and only one DMR was located in an RNase P. While it is possible that these non-gene-associated DMRs are positioned within regions of “junk DNA”, it may be the case that these regions serve regulatory functions for unknown gene targets.
DMGs were annotated with Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Significant clustering was observed based on functional similarity (Figure S3). When all DMRs were considered, 8 annotation clusters were identified (Enrichment Threshold > 2.0, p < 0.01). Representative annotation terms with the smallest p-value from each cluster included ion binding (GO:0043167, Cluster 1), cell morphogenesis (GO:0000902, Cluster 2), cell fraction (GO:0000267, Cluster 3), cell morphogenesis involved in differentiation (GO:0000904, Cluster 4), stereocilium (GO:0032420, Cluster 5), regulation of synaptic transmission (GO:0050804, Cluster 6), triglyceride biosynthetic process (GO:0019432, Cluster 7), and cAMP catabolic process (GO:0006198, Cluster 8).
Investigation of KEGG pathways demonstrated that seven pathways were DMR-enriched (Table 1, Fisher’s exact p-value < 0.05). Using Benjamini–Hochberg correction, only the phosphatidylinositol (PI) signaling system (rno04070) was significantly enriched (corrected p-value = 0.006). Other pathways of interest included insulin signaling pathway (rno04910), axon guidance (rno04360), pathways in cancer (rno05200), purine metabolism (rno00230), glycerolipid metabolism (rno00561), and Fc gamma R-mediated phagocytosis (rno04666). There were three pathways that were more methylated in the MHF group: phosphatidylinositol signaling system (rno04070), apoptosis (rno04210), and Fc gamma R-mediated phagocytosis (rno04666). On the other hand, the insulin signaling pathway (rno04910) was the only pathway that is more methylated in the PHF group.

2.3. Post-Weaning HF Increases Gene Expression While Maternal HF Decreases Gene Expression in the PI and Insulin Signaling Pathways

From the pathway analysis, we chose to further investigate the PI signaling and insulin signaling pathways. Not only are these two pathways interconnected, but they also serve metabolic roles that are known to be perturbed by high calorie diets. Furthermore, the PI and insulin signaling pathways were statistically the most DMR-enriched according to the DAVID analysis (p = 0.006 and p = 0.4, respectively, Table 2). Additionally, we compared methylation and gene expression results to the CON group in order to assess not only differences in timing of the HF exposure, but also to examine whether any HF diet intake might disrupt metabolic pathways compared to control diet (Figure 4).
We combined and condensed the pathways by focusing on 27 genes, 17 of which were differentially methylated and 10 of which were not differentially methylated but served crucial metabolic functions (Figure 4). We first examined the methylation in each of the genes relative to the CON group (Figure 5a,b). Generally, methylation of Dgkg, Ip3k, Pik3c2b, Pde3, and Lar was lower in PHF than CON while methylation of Pklr was higher in PHF than CON. Additionally, methylation was lower in Inpp5, Pik3c2b, Cbl, Lar, and Pklr in MHF than CON while methylation of Ip3k and Pde3 was higher in MHF than in CON. We then performed quantitative PCR for all 27 genes and found 10 that were differentially expressed (p < 0.05). Among the differentially methylated genes, diacylglycerol kinase gamma (Dgkg), inositol-triphosphate 3-kinase (Ip3k), inositol polyphosphate-5-phosphatase (Inpp5), PI3 kinase (Pi3kc2b), phosphodiesterase 3 (Pde3), Cbl proto-oncogene (Cbl), and leukocyte antigen-related protein tyrosine phosphatase (Lar), were up-regulated in PHF rats, while pyruvate kinase (Pklr) was up-regulated in the MHF group (Figure 5c). Of the 10 genes that were not differentially methylated, phosphorylase kinase (Phk) and phosphatidylinositol synthase (PI synthase) were more highly expressed in the PHF group. From this differential gene expression, we were able to visually represent changes within glucose and lipid homeostatic processes between perinatal and post-weaning HF diet (Figure 4).

2.4. Genomic Context of DNA Methylation Is Indicative of Gene Expression

Among the 17 DMGs of interest, we examined the distance of the gene-associated DMR to the nearest CpG island as well as their intragenic locations (Figure 6). Of the 8 differentially expressed DMGs, only one contained a DMR that was located within a CpG island (Pklr, Figure 6a). Interestingly, the remaining 7 DMRs fell outside the 2000 bp shore region (Figure 6a). We also found that except for one, all DMRs were intragenic, with a majority located in introns (Figure 6b).
Lastly, we investigated the intragenic locations of these differentially methylated and differentially expressed regions. Based on qPCR results, a gene expression t-score was calculated for every gene in every animal. Values that were four standard deviations away from the mean were chosen. Relative intragenic position of a given DMR was calculated by dividing the distance of the DMR from the transcription start site (TSS) by the total length of the gene (Figure 7 x-axis). The kernel density estimation was plotted to visualize the distribution of intragenic DMR positions. Interestingly, none of the differentially expressed DMR was in the promoter regions (data not shown). Genes highly expressed in the MHF group tended to be less methylated, and this low methylation was found almost exclusively in the middle of the intragenic region (blue curves, Figure 7a). Furthermore, higher gene expression in the PHF group was correlated with lower methylation near the TSS and higher methylation in downstream regions closer to the transcription end site (TES, yellow vs red curves, Figure 7b). However, we highlight the idea that gene body methylation may play an activating or repressive role in transcriptional regulation [22].

3. Discussion

The present study used deep sequencing technologies to examine DNA methylation profiles in animals exposed to HF diet at different time points. First, to address whether timing of HF diet exposure induced methylation changes, we compared between MHF and PHF rats. Over 1700 DMRs were identified, of which a large portion were in CpG shores and intragenic regions. Gene-associated DMRs also clustered in the insulin and PI signaling pathways. We then compared DNA methylation in the identified regions to lifelong control-fed animals to uncover the impact of HF intake regardless of timing. In these metabolic pathways, we conclude that any HF diet decreases DNA methylation, as both the MHF and PHF groups generally had lower methylation than CON; however, expression of key metabolic regulators tended to be higher in PHF and lower in MHF. Finally, we examined DMR location relative to CpG islands and gene features. We found the greatest number of DMRs located in CpG shores and showed an inverse correlation between gene expression and distance of intragenic DNA methylation from the TSS.
The current study used two diets. The control and HF diets had equal amounts of soybean oil, but the HF diet contained higher amounts of lard. Lard provides more saturated and monounsaturated fat compared to soybean oil. This has been shown to increase adiposity and markers of inflammation [23]. Chronic obesity-induced inflammation also underlies insulin resistance [24]. Interestingly, we found differential methylation in the insulin signaling pathway, so inflammatory factors may be a potential link between HF diet and epigenetic modifications. This can be tested in future experiments using HF diets with different amounts of lard, soybean oil, or fish oil to distinguish the effects of inflammation versus obesity. All other dietary components were kept constant, except for the major carbohydrate source, which was cornstarch in the control diet and sucrose in the HF diet. Sucrose and cornstarch have been shown to differentially impact body weight gain, circulating lipid profile, insulin sensitivity, and hepatic gene expression [25,26,27]. Thus, it is possible that the observed effects could have resulted from dietary carbohydrate differences rather than fat intake, so further investigation should examine individual dietary components to better understand their effects on DNA methylation.
Overall, we observed higher body weight and hepatic lipid accumulation in PHF rats. At post-natal week 12, the PHF group weighed more than CON, but there was no difference in body weight between MHF and PHF. Given the large difference in caloric intake, this finding was surprising. However, while HF diet administration has been shown to induce weight gain during adulthood, high energy intake during the post-weaning period results in only moderate body weight differences [28,29,30,31]. During the post-weaning period, rapid growth occurs, which may mask any additional diet-induced weight gain. We found greater hepatic lipid deposition in the PHF group compared to both CON and MHF. To understand the potential mechanisms mediating these differences, we measured expression of genes that were related to de novo lipogenesis and triglyceride synthesis. We found higher expression in MHF animals and lower expression in PHF rats. The down-regulation of these pathways in PHF animals is not surprising and is likely a direct effect of abundant dietary fat and a reduced need for de novo synthesis. However, the influx of dietary fat likely leads to the increased quantity of hepatic lipid that we observed. On the other hand, MHF rats had higher expression of de novo lipogenesis genes than controls, despite the eating the same post-weaning diet. Previous reports have similarly found increased acetyl-CoA carboxylase (Acc) expression in liver and adipose tissue of maternal HF-fed animals [32,33]. Moreover, disturbances in maternal diet have been shown to drive metabolic outcomes via epigenetic programming [34,35]. In our study, it is unclear whether this was due to maternal programming, considering that there were no methylation differences at these loci at post-natal week 12. Moreover, there were no differences in hepatic lipid accumulation between MHF and CON. Thus, while gene expression might suggest that the MHF group is more efficient at converting excess energy to fat, there are other factors that prevent hepatic accumulation. Future work should quantify lipid export and size of different adipose depots to find whether changes in the liver could affect other organs.
Among 1744 DMRs identified across 21 chromosomes, ~28% of which were located within CpG shores, more than twice the number located in CpG islands (13%). This is in line with previous studies illustrating that DMRs within the same tissue type occur more frequently in shores versus islands [20,21]. We also found many DMRs in intragenic and intergenic regions as opposed to promoter regions, which is expected, as CpG islands tend to be localized in promoters [36]. Although a majority of the DMRs resided in intragenic regions, only one corresponded with an annotated non-coding RNA. We attribute this finding to the fact that information regarding the function of intergenic regions is sparse, especially for the rat. For example, the microRNA database, mirBase, annotates 4694 miRNA sequences in humans and 3232 in mouse, but only 1318 in rat. As non-coding regions become more widely studied and compiled in other model organisms, we suspect that our identified intergenic regions will show more regulatory functions. Over the past decade, ENCODE and similar projects have gradually uncovered the intricacies of the vast genetic regulatory network [37]. Most of these advances have focused on human transcriptional regulators, but as the data becomes publicly available for other model organisms we might be able to assign function to previously unnamed genomic regions.
Functional clustering and pathway analysis revealed DMR enrichment in the insulin and PI signaling pathways. Both are nutrient-sensing pathways important in carbohydrate and lipid metabolism, suggesting that maternal and post-weaning HF dietary exposure may cause distinct metabolic outcomes. Indeed, PHF animals had greater lipid accumulation in the liver than did MHF animals. In the MHF group, we also observed higher methylation and lower gene expression of Pi3k and Pde3, two genes involved in a phosphorylation cascade that inhibits lipolysis. Three separate PI3K genes were more highly methylated in the MHF compared to the PHF group, including PI3K regulatory subunit 1 (Pi3kr1), PI3K regulatory subunit 3 (Pi3kr3), and PI3K catalytic subunit type 2 beta (Pi3kc2b). While Pi3kr1 and Pi3kr3 expression was modestly decreased in the MHF group, only Pi3kc2b expression was significantly reduced. PI3K is a key component of the insulin signaling pathway, but its regulation is also controlled in large part by PI signaling [38]. Previous studies have shown that when members of the PI signaling system are knocked down, PI3K activity and the insulin signaling pathway are dysregulated. For example, Inositol Polyphosphate Phosphatase-Like 1 (Inppl1) knock out mice display insulin hypersensitivity, increased levels of phosphorylated AKT, and protection against diet-induced obesity [39,40]. Moreover, knocking down inositol polyphosphate-4-phosphatase type II (Inpp4b), a component of the PI signaling pathway, results in the deregulation the PI3K-AKT pathway such that the magnitude and duration of insulin-stimulated AKT activation are significantly altered [41]. In our experiment, we observed four PI signaling genes (Ip3k, Inpp5, DG kinase, and PI synthase) that were more highly expressed in the PHF than MHF group. This difference in PI signaling genes may reflect the distinct regulation of PI3K between the two groups. The observed differences in methylation of PI3K as well as the gene expression changes of its regulators in the PI signaling pathway indicate that PI3K’s function as a catalytic protein may be altered in PHF compared to MHF diet exposure.
In the MHF group, we observed higher methylation within the first intron of Pde3, which was associated with the decreased expression of Pde3. Inhibition of Pde3 sufficiently blocks the antilipolytic action of insulin via its role in hydrolyzing cAMP and cGMP [42]. In fact, an adipocyte Pde3 knockout has been shown to result in higher levels of insulin-stimulated lipogenesis [43]. Similarly, knockout mice had altered phosphorylation state of insulin- and cAMP-signaling components as well as higher lipolysis upon catecholamine stimulation. Thus, the decrease in Pde3 expression in the MHF group suggests that insulin may have less of an inhibitory effect on the catabolism of triacylglercerol and contribute to the lower hepatic lipid accumulation in the MHF group.
In addition to identifying genes that have annotated functions in PI and insulin signaling pathways, our results corroborate GWAS analyses that have investigated the association between novel loci and traits such as body–mass index (BMI), waist circumference, triglyceride levels, cholesterol and diabetes [44,45,46]. Specifically, previous studies have found that SNPs within Cdh23 were correlated with BMI [45], SNPs in Tph2 were correlated with obesity [44], SNPs in Prox1 were correlated with glucose metabolism [47], SNPs in Lipc were correlated with cholesterol [48,49,50], and SNPs in Pnpla3 were correlated with hepatic lipid content [51]. We found DMRs associated with each of these genes: Cdh23, Tph2, and Lipc were more methylated in MHF while Prox1 and Pnpla3 were more methylated in PHF. Although reproducibility across GWAS studies is often lacking, we provide evidence that strengthens past findings. Furthermore, we suggest that future research may benefit from combining multiple sequencing modalities to identify robust biomarkers of metabolic health.
Beyond analyzing the role of HF diet in controlling metabolic processes, we also attempted to characterize the differences in epigenetic landscape between MHF and PHF and to examine how DNA methylation is associated with gene expression. We identified that most DMRs were in CpG shores. CpG island methylation has been studied extensively in the context of islands at TSSs where higher island methylation represses gene expression. Reports have shown that intragenic CpG islands are more likely to contain alternative promoters and other regulatory features [14,52]. Differential DNA methylation in CpG shores has also been shown to be associated with gene expression [20], but it is unknown whether CpGs outside CG-dense regions perform similar regulatory functions. In our study, we also found that out of the DMRs associated with insulin and PI signaling, five were located in CpG islands, one was in a CpG shore, and the remaining 11 were more than 2 kb away from the nearest island in the open sea. CpGs in the open sea have been shown to be generally hypomethylated [53,54,55,56]. While in cancer, this hypomethylation has been associated with chromosomal instability, open sea methylation also naturally changes with during development. From birth to age 10, CpG islands tend to gain methylation whereas open sea CpGs tend to lose methylation [53,57,58,59,60,61]. Thus, it may be the case that the identified DMRs are a result of diet-mediated alterations in developmental processes. Future work should focus on defining the role of DNA methylation outside of CpG islands to understand how metabolic disturbances impact such loci.
Through the investigation of the relationship between intragenic DMR location and gene expression, we observed that all differentially expressed DMGs contained DMRs in the gene body. Five out of 8 differentially expressed DMGs contained a DMR exclusively in an intronic region while the remaining 3 genes had a DMR spanning both an intron and an exon. The role of intronic DNA methylation is not well understood. Past investigation has hypothesized that it controls alternative splicing [62,63]. Others have reported regulatory elements within introns that when methylated have either repressive or activating role in gene expression [64,65]. Additionally, reports have speculated that the balance between intronic and exonic DNA methylation to be important in determining nucleosome spacing and subsequent Pol II binding [66,67,68,69]. Thus, it is difficult to discern which mechanism is at play in each of the identified genes.
There is currently no clear consensus as to the function of gene body DNA methylation. It has been proposed to mediate chromatin structure, splicing, and transcriptional kinetics; however, it is clear that gene body methylation does not necessarily act in the same way as traditionally studied promoter methylation [70]. We observed that DMR distance from the TSS was associated with gene expression. Higher expression was associated with higher methylation in more distal regions. These findings are supported by previous reports that gene expression increased with higher gene body methylation, but did not distinguish among different intragenic locations [71,72]. Other studies showed that highly expressed genes were characterized by low methylation in the gene body near the TSS and high intragenic methylation closer to the 3’ TES [22,73]. Similarly, low gene expression is correlated with high methylation in the first exon [74].
HF diet may affect the epigenetic landscape via numerous mechanisms. HF diet can incite inflammation and hormonal changes that have been associated with altered DNA methylation patterns [75]. However, it may be the case that HF diet acts directly on epigenetic modifiers and methylation processes. For instance, short-chain fatty acids such as butyrate and acetate have been shown to inhibit histone deacetylases (HDACs) and alter DNA methylation [76]. Furthermore, acetyl-CoA, the end-product of beta oxidation, is a rate-limiting cofactor in histone acetylase (HAT) activity [77]. Due to the bidirectional recruitment and interaction between histone modifications and DNA methylation, HF diet has the potential to perturb chromatin state and severely impact transcription.
The observed methylation discrepancies that occur during the perinatal versus post-weaning periods may be explained by the drastic differences in methylation events between the two groups. During the perinatal period, an initial wave of demethylation is followed by de novo methylation by methyltransferases DNMT3a and DNMT3b [9]. However, during the post-weaning period, methylation patterns have already been established, and further upkeep relies solely on the maintenance DNA methyltransferase DNMT1 [78]. HF diet may interfere differently with certain types of methyltransferases. Previous reports have shown that altering dietary folate and choline has different effects on the expression of DNMT1 and DNMT3 and that these changes vary with length and timing of exposure [79,80,81,82]. Another study showed that HF diet differentially affected the expression of DNMT1, DNMT3a, and DNMT3b, as well as their binding to the DNA [83]. After eight weeks of post-weaning HF diet intake, DNMT1 binding at the leptin receptor promoter significantly decreased, and DNMT3b binding significantly increased. The study also investigated methyl binding domain protein 2 (MBD2), which may play a role in DNA demethylation. In addition to DNA methyltransferases, it is possible that HF diet differentially affects mediators of active DNA demethylation. In the future, next generation sequencing will allow for a genome-wide view of epigenetic modifier binding in response to dietary challenges.
We provide novel insight into epigenetic programming by HF diet; however, the current study has limitations that should be addressed in future investigation. First, we examine only male offspring in our study. Previously, we have found sex-specific physiological and molecular changes resulting from HF diet intake [84,85,86]. In particular, we have observed large differences in gene expression and DNA methylation in male rats [19,87,88]. In the current study, we build upon prior results in male offspring; however, it is unknown whether these results would be broadly applicable to females. Thus, subsequent studies should examine whether DNA methylation patterns are robust across sexes. Furthermore, we found molecular changes, but we do not provide evidence that these changes induce functional consequences. As this was a genome-wide exploratory study, it was unclear whether certain pathways would be more enriched for changes in DNA methylation. However, now that PI and insulin signaling have been identified, follow-up studies should measure insulin sensitivity and glucose tolerance. Finally, we measured DNA methylation and gene expression at one time point at post-natal week 12, but it is unknown whether DNA methylation is labile between birth, weaning, and further into the post-natal period. It could be the case that younger animals are better equipped than older animals to combat disturbances in the methylome to normalize gene expression. Further work should explore molecular changes at multiple time points to discern whether epigenetic alterations have a similar impact at different life stages.

4. Materials and Methods

4.1. Animals and Diets

Timed-pregnant Sprague Dawley rat dams (Charles River Laboratories) were divided into two groups. One group of dams was fed a standard AIN93G diet (n = 12; Research Diets, Inc.; 16%, 64%, 20% calories from fat, carbohydrate, and protein, respectively) while the other group was fed a HF diet (n = 12; Research Diets, Inc.; 45%, 35%, 20% calories from fat, carbohydrate, and protein, respectively) during gestation and lactation [85]. On post-natal day 21, male pups from control diet-fed dams were weaned onto control diet (n = 10) or HF diet (n = 10). Male pups from HF-fed dams were weaned onto the AIN93G control diet (n = 10). This created three groups of pups: control (CON), post-weaning high fat (PHF), and maternal high fat (MHF). Diet composition is detailed in Table S1.
Rats were individually housed in standard polycarbonate cages in a humidity- and temperature-controlled room on a 12-hour light-dark cycle with ad libitum access to food and drinking water. Body weight and food intake was measured weekly. Each treatment group was kept on their respective post-weaning diet until 12 weeks of age. Before sacrifice, animals were fasted for 12 h and received free access to water. Euthanasia was performed via CO2 followed by decapitation. The median lobe of the liver was immediately frozen in liquid nitrogen and stored at −70 °C. Institutional and governmental regulations regarding the ethical use of animals were followed during the study. The protocol for ethical use of animals for this study was approved by the University of Illinois Institutional Animal Care and Use Committee (IACUC protocol no. 09112).

4.2. Histological Analysis

Frozen liver samples were embedded in Tissue-Tek OCT compound (VWR, cat. #25608-930) and sectioned, stained, imaged, and quantified using a previously published protocol [86]. Briefly, OCT-embedded liver tissues were sectioned to 7 µm and stained with hematoxylin and eosin (H&E) and Oil Red O solution (Newcomer Supply, cat. #1277A). The slides were imaged using the Nanozoomer imaging system at the Carl R. Woese Institute for Genomic Biology core facilities at the University of Illinois (Hamamatsu Photonics, Hamamatsu City, Japan). Quantification of lipid accumulation was normalized to total protein as previously described [87].

4.3. Genomic DNA Isolation

Ten mg of liver tissue was ground in liquid nitrogen and genomic DNA was extracted in 600 μL of Extraction Buffer (50 mM Tris, pH 8.0, 1 mM EDTA, pH 8.0, 0.5% SDS, 1 mg/mL Proteinase K) at 55 °C overnight. Lysate were centrifuged and supernatant was collected and mixed with phenol/chloroform. The mixture was transferred to a phase lock gel tube (Fisher Scientific, cat. #NC1092951) and centrifuged at 16,000× g for 5 min for phase separation. The upper phase was transferred and incubated with 1 μL of RNase (Roche, 10 mg/mL) for 1 h at 37 °C. The extraction was then repeated, and the resulting purified DNA was precipitated with 1/10 vol of 3 M Na Acetate (pH 5.2) and 2.5 vol of 100% ethanol. The DNA pellet was washed with 70% ethanol and resuspended in TE. DNA electrophoresis was used to confirm the integrity of the extracted genomic DNA.

4.4. MeDIP-Seq and MRE-seq Sequencing

High-throughput sequencing was conducted using liver genomic DNA from each group by complementary MeDIP and MRE methods to measure methylated and unmethylated DNA, respectively. Animals were chosen through an extensive screening process in which gene expression and histology were measured and the best representatives from each group were used for sequencing. Moreover, additional precautions were taken to avoid false positives by using two complementary detection techniques (MeDIP-seq and MRE-seq), setting stringent statistical criteria (details follow), and validating in multiple animals using a probe-based method (details follow).
MeDIP-seq and MRE-seq were conducted in the laboratory of Ting Wang at the University of Washington (St. Louis, MO, USA) using previously established protocols [89]. Briefly, the MeDIP-seq method requires sonication of the DNA into ~500 bp fragments followed by immunoprecipitation using a 5mC-specific antibody. High-throughput sequencing is then used to quantify enriched methylated DNA sequences. MRE-seq method uses restriction enzymes that are sensitive to CpG methylation to cut and enrich unmethylated DNA fragments, thus providing additional coverage and resolution. For MRE-seq, five restriction enzymes were used which have the following recognition sequences CCGG, CCGC, GCGC, ACGT, and CGCG.

4.5. Probe-Based Analysis for Validation of DNA Methylation

A previously established, probe-based method, Quantitative Analysis of Methylated Alleles (QAMA), was modified and used to validate methylation measured by MeDIP-seq and MRE-seq [90]. The method uses two different fluorescently labeled probes, one for methylated DNA and the other for unmethylated DNA, to detect methylated versus unmethylated region of DNA. Genomic DNA was isolated from 16 individual animals using the DNeasy Tissue Kit (Zymo Research, Irvine, CA, USA) and treated with sodium bisulfite reagent using the EZ Methylation Gold kit (Zymo Research). After bisulfite conversion, the DNA sample was diluted to 10 ng/μL for quantitative PCR analysis. The PCR was performed in a 96-well optical plate with a final reaction volume of 20 μL. Each reaction contained 10 μL of Taqman® Universal PCR Master Mix without AmpErase®, 4 μL of nuclease-free water, 1.2 μL of 5 μM each of the forward and reverse primers, and 0.8 μL of 2.5 μM each of the fluorescently labeled methylated and unmethylated probes. The reaction was as follows: 95 °C for 15 min to activate AmpliTaq Gold DNA polymerase followed by 45 cycles of 95 °C for 15 s and 60 °C for 1 min in StepOnePlus real-time PCR system (Applied Biosystems, Foster City, CA, USA). Primer and probe were designed using Primer Express Software 3.0.1 (Applied Biosystems; Table 2).
Forward and reverse primers were designed to amplify the bisulfite converted DMR sequence. Primers had a melting temperature (Tm) of 58–60 °C and were associated with regions free of CpGs and nucleotide polymorphisms. Two different probes were selected, each with a Tm of 68–70 °C. The methylated probe binds to a sequence with protected (methylated) CpGs while the unmethylated probe binds to a sequence with unprotected (bisulfite converted) CpGs. The methylated probe was labeled with VIC at 5′ and the unmethylated probe was labeled with FAM at 5′, two fluorophores that have absorption/emission spectrums of 538/554 nm and 494/518 nm, respectively. Both probes were labeled with an MGB quencher at the 3′-end. The relative amount of VIC and FAM fluorescence in each reaction was measured and quantities of methylated and unmethylated DNA were analyzed based on VIC and FAM standard curves, respectively. Percentage of methylated DNA was calculated with the following equation:
(quantity of methylated DNA)/(quantity of methylated DNA + quantity of unmethylated DNA)

4.6. Methylation Specific PCR

In addition to the probe-based method discussed above, methylation specific PCR (MSP) was used to validate sequencing results. Protocols from previous experimentation were used for primer design, genomic DNA isolation, bisulfite conversion, and qPCR [19]. The percentage of methylated DNA was calculated as the ratio discussed above. All MSP primer information can be found in Table 3.

4.7. Gene Expression Analysis

Fifty mg of the liver tissue was homogenized and prepared for mRNA isolation using TRIzol® Reagent (Life Technologies, Carlsbad, CA, USA). Direct-zol™ RNA MiniPrep columns (Zymo Research) with in-column DNase I digestion were used to extract total RNA. RNA concentration and purity were measured using Nanodrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA). Reverse transcription was performed using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). The mRNA mixed with reverse transcriptase was incubated in the 2720 Thermal Cycler (Applied Biosystems) at 25 °C for 10 min, 37 °C for 2 h and 85 °C for 5 min. A serially diluted standard curve and all cDNA samples were amplified using Power SYBR® Green Master Mix (Life Technologies) in a StepOnePlus™ Real-Time PCR System (Life Technologies) at 95 °C for 10 min, followed by 35 cycles of 95 °C for 15 s and 60°C for 1 min. A melting curve was included to check the purity of the PCR product. Primers used for quantitative PCR are shown in Table S2. All quantities reported are normalized to β-actin, as it has been shown as a reliable internal control gene and has been used extensively in the past [84,91,92].

4.8. DMR Identification

Given the resolution of MeDIP-seq, we partitioned the genome into bins of 500 bp in size and ran MethylMnM to test for statistically different methylation status within each bin (FDR q-value < 0.05). The software also outputs MeDIP and MRE RPKM for each animal which was used in subsequent data visualization and analysis.

4.9. Genomic Location of DMRs

DMR location within the genomic context was examined to identify DMR location relative to CpG islands and genes. CpG islands satisfied three criteria: (i) sequences were of length 200 bp or longer, (ii) had guanine-cytosine (GC) content of 50% or greater, and (iii) had a ratio of observed to expected CpGs of 0.6 or higher. Shores were defined as the 2000 bp regions flanking either side of an island [20]. Differentially methylated genes (DMGs) were classified as intragenic, promoter, or downstream regions, where promoter and downstream regions were defined as the 1500 bp sequence upstream of the TSS and the 1500 bp region downstream of the transcription end site, respectively [88].

4.10. Pathway Analysis and Functional Clustering

DMGs were annotated using Gene Ontology (GO) terms (Biological Processes (BP), Cellular Component (CC), and Molecular Function (MF)) as well as KEGG pathway enrichment. Functional annotation of DMGs was then used to compute Pearson correlation coefficients (PCC) for pairs of genes to reflect functional similarity between genes. Analysis was performed with the Guide for Association Index for Networks (GAIN) tool available online through the University of Minnesota (http://franklin-umh.cs.umn.edu/similarity_index/index.php) [93]. Given a bipartite network of independent sets, X (genes) and Y (annotation terms), in which interactions are either present or absent, GAIN calculates PCC by considering the likelihood of observing overlap between two vertices in set X given their vertex degree and the total number of Y vertices in the network. A PCC value of 1 represents a perfect functional similarity between two genes, whereas a PCC of 0 indicates no functional similarity.
Further pathway analysis of DMGs was conducted using DAVID Bioinformatics Resources 6.7 (http://david.abcc.ncifcrf.gov/ [94]). Clusters were formed based on the same GO categories and KEGG pathway enrichments used in the PCC calculations. DAVID analysis generated a hierarchical organization in which related genes were categorized into annotation terms and annotation terms were categorized into annotation clusters. Functional annotation clustering was performed using the following criteria: Similarity Threshold = 0.85 (the minimum similarity value to be considered biological significant), Similarity Term Overlap = 3 (the minimum number of annotation terms overlapped between two genes), Count Threshold = 3 (the minimum number of genes to constitute an annotation term), Minimum Annotation Terms = 3 (the minimum number of annotation terms to constitute an annotation cluster), modified Fisher’s exact p-value < 0.01 (lower value indicates higher gene-enrichment within an annotation term), and Enrichment Score > 2.0 (based on the geometric mean of annotation term p-values within an annotation cluster).

4.11. General Statistical Analysis

Body weight and food intake were analyzed using ANOVA with repeated-measures followed by post-hoc Tukey HSD test. All gene expression and methylation validation qPCR data were analyzed with R Statistical Software version 3.1.2. Methylation validation was analyzed using a two-tailed student’s t-test. Gene expression was analyzed via one-way ANOVA followed by Tukey HSD test.

5. Conclusions

To our knowledge, this is the first study to show widespread epigenetic changes in the insulin and PI signaling pathways in response to HF diet. Moreover, we demonstrated that HF diet exposure during either the perinatal or post-weaning period is crucial in determining hepatic DNA methylation profile. Given our results, we predict that perinatal HF diet alters maternal physiology which changes de novo epigenetic modifications in the fetus that persist into adulthood. Post-weaning HF diet, on the other hand, likely works through a more direct mechanism in which HF diet alters methylation maintenance in the mature liver. These distinct processes ultimately result in unique hepatic methylation patterns in the PI and insulin signaling pathways. We also illustrated that genomic location of DMRs, namely their position relative to CpG islands and TSSs, was indicative of gene expression. Our data suggest that maternal and post-weaning HF exposure differentially affect the epigenome within specific genomic contexts.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/13/3229/s1.

Author Contributions

Conceptualization, L.M. and Y-X.P.; methodology, L.M and H.C.; formal analysis, L.M. and H.W.; investigation, L.M. and P.M.J.; resources, Y-X.P.; data curation, L.M.; writing—original draft preparation, L.M.; writing—review and editing, L.M., H.C.; supervision, Y-X.P.; project administration, H.C.; funding acquisition, Y-X.P.

Funding

This project was supported by UIUC Research Board grant #12192, the USDA Cooperative State Research, Education and Extension Service, Hatch project numbers # ILLU-698-923 and ILLU-698-391, and the Data Purchase Program from University Library.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

HFhigh-fat diet
MeDIP-seqmethylated DNA immunoprecipitation with high-throughput sequencing
MRE-seqmethylation-sensitive restriction enzyme digestion sequencing
DMRdifferentially methylation regions
DMGdifferentially methylated genes
Mc4rmelanocortin 4 receptor
Pparaperoxisome proliferator-activated receptor α
Tlr-1 and Tlr-2toll-like receptors 1 and 2
Cdkn1acyclin-dependent kinase inhibitor
H&Ehematoxylin and eosin
QAMAQuantitative Analysis of Methylated Alleles
MSPmethylation specific PCR
GOGene Ontology
KEGGKyoto Encyclopedia of Genes and Genomes
PCCPearson correlation coefficients

References

  1. Chang, G.Q.; Gaysinskaya, V.; Karatayev, O.; Leibowitz, S.F. Maternal high-fat diet and fetal programming: Increased proliferation of hypothalamic peptide-producing neurons that increase risk for overeating and obesity. J. Neurosci. 2008, 28, 12107–12119. [Google Scholar] [CrossRef] [PubMed]
  2. Srinivasan, M.; Katewa, S.D.; Palaniyappan, A.; Pandya, J.D.; Patel, M.S. Maternal high-fat diet consumption results in fetal malprogramming predisposing to the onset of metabolic syndrome-like phenotype in adulthood. Am. J. Physiol. Endocrinol. Metab. 2006, 291, E792–E799. [Google Scholar] [CrossRef] [PubMed]
  3. Milagro, F.I.; Campion, J.; García-Díaz, D.F.; Goyenechea, E.; Paternain, L.; Martínez, J.A. High fat diet-induced obesity modifies the methylation pattern of leptin promoter in rats. J. Physiol. Biochem. 2009, 65, 1–9. [Google Scholar] [CrossRef] [PubMed]
  4. Widiker, S.; Kärst, S.; Wagener, A.; Brockmann, G.A. High-fat diet leads to a decreased methylation of theMc4r gene in the obese BFMI and the lean B6 mouse lines. J. Appl. Genet. 2010, 51, 193–197. [Google Scholar] [CrossRef] [PubMed]
  5. Jiang, M.; Zhang, Y.; Liu, M.; Lan, M.S.; Fei, J.; Fan, W.; Gao, X.; Lu, D. Hypermethylation of hepatic glucokinase and L-type pyruvate kinase promoters in high-fat diet-induced obese rats. Endocrinology 2011, 152, 1284–1289. [Google Scholar] [CrossRef] [PubMed]
  6. Ge, Z.J.; Luo, S.M.; Lin, F.; Liang, Q.X.; Huang, L.; Wei, Y.C.; Hou, Y.; Han, Z.M.; Schatten, H.; Sun, Q.Y. DNA methylation in oocytes and liver of female mice and their offspring: Effects of high-fat-diet-induced obesity. Environ. Health Perspect. 2014, 122, 159–164. [Google Scholar] [CrossRef] [PubMed]
  7. Ding, Y.; Li, J.; Liu, S.; Zhang, L.; Xiao, H.; Li, J.; Chen, H.; Petersen, R.B.; Huang, K.; Zheng, L. DNA hypomethylation of inflammation-associated genes in adipose tissue of female mice after multigenerational high fat diet feeding. Int. J. Obes. 2014, 38, 198–204. [Google Scholar] [CrossRef]
  8. Dudley, K.J.; Sloboda, D.M.; Connor, K.L.; Beltrand, J.; Vickers, M.H. Offspring of mothers fed a high fat diet display hepatic cell cycle inhibition and associated changes in gene expression and DNA methylation. PLoS ONE 2011, 6, e21662. [Google Scholar] [CrossRef]
  9. Okano, M.; Bell, D.W.; Haber, D.A.; Li, E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell 1999, 99, 247–257. [Google Scholar] [CrossRef]
  10. Feil, R.; Fraga, M.F. Epigenetics and the environment: Emerging patterns and implications. Nat. Rev. Genet. 2011, 13, 97–109. [Google Scholar] [CrossRef]
  11. Feil, R.; Walter, J.; Allen, N.D.; Reik, W. Developmental control of allelic methylation in the imprinted mouse Igf2 and H19 genes. Development 1994, 120, 2933–2943. [Google Scholar] [PubMed]
  12. Lister, R.; Ecker, J.R. Finding the fifth base: Genome-wide sequencing of cytosine methylation. Genome Res. 2009, 19, 959–966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lister, R.; Gregory, B.D.; Ecker, J.R. Next is now: New technologies for sequencing of genomes, transcriptomes, and beyond. Curr. Opin. Plant Biol. 2009, 12, 107–118. [Google Scholar] [CrossRef] [PubMed]
  14. Maunakea, A.K.; Nagarajan, R.P.; Bilenky, M.; Ballinger, T.J.; D’Souza, C.; Fouse, S.D.; Johnson, B.E.; Hong, C.; Nielsen, C.; Zhao, Y.; et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature 2010, 466, 253–257. [Google Scholar] [CrossRef] [PubMed]
  15. Stevens, M.; Cheng, J.B.; Li, D.; Xie, M.; Hong, C.; Maire, C.L.; Ligon, K.L.; Hirst, M.; Marra, M.A.; Costello, J.F.; et al. Estimating absolute methylation levels at single CpG resolution from methylation enrichment and restriction enzyme sequencing methods. Genome Res. 2013, 23, 1541–1553. [Google Scholar] [CrossRef] [PubMed]
  16. Aldoretta, P.W.; Carver, T.D.; Hay, W.W., Jr. Maturation of glucose-stimulated insulin secretion in fetal sheep. Neonatology 1998, 73, 375–386. [Google Scholar] [CrossRef] [PubMed]
  17. Hay, W.W., Jr. Recent observations on the regulation of fetal metabolism by glucose. J. Physiol. 2006, 572, 17–24. [Google Scholar] [CrossRef] [PubMed]
  18. Harris, R.A.; Wang, T.; Coarfa, C.; Nagarajan, R.P.; Hong, C.; Downey, S.L.; Johnson, B.E.; Fouse, S.D.; Delaney, A.; Zhao, Y.; et al. Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications. Nat. Biotechnol. 2010, 28, 1097–1105. [Google Scholar] [CrossRef]
  19. Moody, L.; Pan, Y.-X.; Chen, H. Postnatal diet remodels hepatic DNA methylation in metabolic pathways established by a maternal high-fat diet. Epigenomics 2017, 9, 1387–1402. [Google Scholar] [CrossRef] [PubMed]
  20. Irizarry, R.A.; Wu, H.; Feinberg, A.P. A species-generalized probabilistic model-based definition of CpG islands. Mamm. Genome 2009, 20, 674–680. [Google Scholar] [CrossRef] [Green Version]
  21. Doi, A.; Park, I.-H.; Wen, B.; Murakami, P.; Aryee, M.J.; Irizarry, R.; Herb, B.; Ladd-Acosta, C.; Rho, J.; Loewer, S.; et al. Differential methylation of tissue- and cancer-specific CpG island shores distinguishes human induced pluripotent stem cells, embryonic stem cells and fibroblasts. Nat. Genet. 2009, 41, 1350–1353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ball, M.P.; Li, J.B.; Gao, Y.; Lee, J.-H.; Leproust, E.M.; Park, I.-H.; Xie, B.; Daley, G.Q.; Church, G.M. Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat. Biotechnol. 2009, 27, 361–368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wang, X.; Cheng, M.; Zhao, M.; Ge, A.; Guo, F.; Zhang, M.; Yang, Y.; Liu, L.; Yang, N. Differential effects of high-fat-diet rich in lard oil or soybean oil on osteopontin expression and inflammation of adipose tissue in diet-induced obese rats. Eur. J. Nutr. 2013, 52, 1181–1189. [Google Scholar] [CrossRef] [PubMed]
  24. Xu, H.; Barnes, G.T.; Yang, Q.; Tan, G.; Yang, D.; Chou, C.J.; Sole, J.; Nichols, A.; Ross, J.S.; Tartaglia, L.A.; et al. Chronic inflammation in fat plays a crucial role in the development of obesity-related insulin resistance. J. Clin. Investig. 2003, 112, 1821–1830. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, X.; Ogawa, H.; Kishida, T.; Ebihara, K. The effect of high-amylose cornstarch on lipid metabolism in OVX rats is affected by fructose feeding. J. Nutr. Biochem. 2010, 21, 89–97. [Google Scholar] [CrossRef] [PubMed]
  26. Blé-Castillo, J.L.; Aparicio-Trápala, M.A.; Juárez-Rojop, I.E.; Torres-Lopez, J.E.; Méndez, J.D.; Aguilar-Mariscal, H.; Olvera-Hernández, V.; Palma-Cordova, L.C.; Diaz-Zagoya, J.C. Differential Effects of High-Carbohydrate and High-Fat Diet Composition on Metabolic Control and Insulin Resistance in Normal Rats. Int. J. Environ. Res. Public Health 2012, 9, 1663–1676. [Google Scholar] [CrossRef] [Green Version]
  27. Hallfrisch, J.; Cohen, L.; Reiser, S. Effects of Feeding Rats Sucrose in a High-Fat Diet. J. Nutr. 1981, 111, 531–536. [Google Scholar] [CrossRef]
  28. Chambers, T.J.G.; Morgan, M.D.; Heger, A.H.; Sharpe, R.M.; Drake, A.J.; Heger, A. High-fat diet disrupts metabolism in two generations of rats in a parent-of-origin specific manner. Sci. Rep. 2016, 6, 31857. [Google Scholar] [CrossRef] [Green Version]
  29. Ghibaudi, L.; Cook, J.; Farley, C.; Van Heek, M.; Hwa, J.J. Fat intake affects adiposity, comorbidity factors, and energy metabolism of sprague-dawley rats. Obes. Res. 2002, 10, 956–963. [Google Scholar] [CrossRef]
  30. Loche, E.; Blackmore, H.L.; Carpenter, A.A.; Beeson, J.H.; Pinnock, A.; Ashmore, T.J.; Aiken, C.E.; De Almeida-Faria, J.; Schoonejans, J.M.; Giussani, D.A.; et al. Maternal diet-induced obesity programmes cardiac dysfunction in male mice independently of post-weaning diet. Cardiovasc. Res. 2018, 114, 1372–1384. [Google Scholar] [CrossRef] [Green Version]
  31. Ibáñez, C.A.; Erthal, R.P.; Ogo, F.M.; Peres, M.N.C.; Vieira, H.R.; Conejo, C.; Tófolo, L.P.; Francisco, F.A.; Silveira, S.D.S.; Malta, A.; et al. A High Fat Diet during Adolescence in Male Rats Negatively Programs Reproductive and Metabolic Function Which is Partially Ameliorated by Exercise. Front. Physiol. 2017, 8, 807. [Google Scholar] [CrossRef] [PubMed]
  32. Almeida, M.M.; Dias-Rocha, C.P.; Souza, A.S.; Muros, M.F.; Mendonca, L.S.; Pazos-Moura, C.C.; Trevenzoli, I.H. Perinatal maternal high-fat diet induces early obesity and sex-specific alterations of the endocannabinoid system in white and brown adipose tissue of weanling rat offspring. Br. J. Nutr. 2017, 118, 788–803. [Google Scholar] [CrossRef] [Green Version]
  33. Oliveira, L.S.; Souza, L.L.; Souza, A.F.P.; Cordeiro, A.; Kluck, G.E.G.; Atella, G.C.; Trevenzoli, I.H.; Pazos-Moura, C.C. Perinatal maternal high-fat diet promotes alterations in hepatic lipid metabolism and resistance to the hypolipidemic effect of fish oil in adolescent rat offspring. Mol. Nutr. Food Res. 2016, 60, 2493–2504. [Google Scholar] [CrossRef] [PubMed]
  34. Keleher, M.R.; Zaidi, R.; Shah, S.; Oakley, M.E.; Pavlatos, C.; El Idrissi, S.; Xing, X.; Li, D.; Wang, T.; Cheverud, J.M. Maternal high-fat diet associated with altered gene expression, DNA methylation, and obesity risk in mouse offspring. PLoS ONE 2018, 13, e0192606. [Google Scholar] [CrossRef] [PubMed]
  35. Martínez, D.; Pentinat, T.; Ribó, S.; Daviaud, C.; Bloks, V.W.; Cebria, J.; Villalmanzo, N.; Kalko, S.G.; Ramon-Krauel, M.; Diaz, R.; et al. In utero undernutrition in male mice programs liver lipid metabolism in the second-generation offspring involving altered Lxra DNA methylation. Cell Metab. 2014, 19, 941–951. [Google Scholar] [CrossRef] [PubMed]
  36. Illingworth, R.S.; Bird, A.P. CpG islands-’a rough guide’. FEBS Lett. 2009, 583, 1713–1720. [Google Scholar] [CrossRef] [PubMed]
  37. Consortium, E.P. An integrated encyclopedia of DNA elements in the human genome. Nature 2012, 489, 57–74. [Google Scholar] [CrossRef]
  38. Shepherd, P.R.; Withers, D.J.; Siddle, K. Phosphoinositide 3-kinase: The key switch mechanism in insulin signalling. Biochem. J. 1998, 333, 471–490. [Google Scholar] [CrossRef]
  39. Sleeman, M.W.; Wortley, K.E.; Lai, K.-M.V.; Gowen, L.C.; Kintner, J.; Kline, W.O.; Garcia, K.; Stitt, T.N.; Yancopoulos, G.D.; Wiegand, S.J.; et al. Absence of the lipid phosphatase SHIP2 confers resistance to dietary obesity. Nat. Med. 2005, 11, 199–205. [Google Scholar] [CrossRef]
  40. Clément, S.; Krause, U.; Desmedt, F.; Tanti, J.-F.; Behrends, J.; Pesesse, X.; Sasaki, T.; Penninger, J.; Doherty, M.; Malaisse, W.; et al. The lipid phosphatase SHIP2 controls insulin sensitivity. Nature 2001, 409, 92–97. [Google Scholar] [CrossRef]
  41. Gewinner, C.; Wang, Z.C.; Richardson, A.; Teruya-Feldstein, J.; Etemadmoghadam, D.; Bowtell, D.; Barretina, J.; Lin, W.M.; Rameh, L.; Salmena, L.; et al. Evidence that inositol polyphosphate 4-phosphatase type II is a tumor suppressor that inhibits PI3K signaling. Cancer Cell 2009, 16, 115–125. [Google Scholar] [CrossRef] [PubMed]
  42. Eriksson, H.; Ridderstråle, M.; Degerman, E.; Ekholm, D.; Smith, C.J.; Manganiello, V.C.; Belfrage, P.; Tornqvist, H. Evidence for the key role of the adipocyte cGMP-inhibited cAMP phosphodiesterase in the antilipolytic action of insulin. Evidence for the key role of the adipocyte cGMP-inhibited cAMP phosphodiesterase in the antilipolytic action of insulin. Biochim. Biophys. Acta 1995, 1266, 101–107. [Google Scholar] [CrossRef]
  43. Choi, Y.H.; Park, S.; Hockman, S.; Zmuda-Trzebiatowska, E.; Svennelid, F.; Haluzík, M.; Gavrilova, O.; Ahmad, F.; Pépin, L.; Napolitano, M.; et al. Alterations in regulation of energy homeostasis in cyclic nucleotide phosphodiesterase 3B-null mice. J. Clin. Investig. 2006, 116, 3240–3251. [Google Scholar] [CrossRef] [PubMed]
  44. Cheverud, J.M.; Ehrich, T.H.; Hrbek, T.; Kenney, J.P.; Pletscher, L.S.; Semenkovich, C.F. Quantitative trait loci for obesity- and diabetes-related traits and their dietary responses to high-fat feeding in LGXSM recombinant inbred mouse strains. Diabetes 2004, 53, 3328–3336. [Google Scholar] [CrossRef] [PubMed]
  45. Dobrin, R.; Chudin, E.; Suver, C.; Beaulaurier, J.; Zhang, B.; Castro, V.; Zhu, J.; Wang, S.; Molony, C.; Lum, P.Y.; et al. A survey of the genetics of stomach, liver, and adipose gene expression from a morbidly obese cohort. Genome Res. 2011, 21, 1008–1016. [Google Scholar]
  46. Welter, D.; MacArthur, J.; Morales, J.; Burdett, T.; Hall, P.; Junkins, H.; Klemm, A.; Flicek, P.; Manolio, T.; Hindorff, L.; et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014, 42, D1001–D1006. [Google Scholar] [CrossRef] [PubMed]
  47. Dupuis, J.; Langenberg, C.; Prokopenko, I.; Saxena, R.; Soranzo, N.; Jackson, A.U.; Wheeler, E.; Glazer, N.L.; Bouatia-Naji, N.; Gloyn, A.L.; et al. New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat. Genet. 2010, 42, 105–116. [Google Scholar] [CrossRef]
  48. Aulchenko, Y.S.; Ripatti, S.; Lindqvist, I.; Boomsma, D.; Heid, I.M.; Pramstaller, P.P.; Penninx, B.W.J.H.; Janssens, A.C.J.W.; Wilson, J.F.; Spector, T.; et al. Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat. Genet. 2009, 41, 47–55. [Google Scholar]
  49. Kathiresan, S.; Willer, C.J.; Peloso, G.M.; Demissie, S.; Musunuru, K.; Schadt, E.E.; Kaplan, L.; Bennett, D.; Li, Y.; Tanaka, T.; et al. Common variants at 30 loci contribute to polygenic dyslipidemia. Common variants at 30 loci contribute to polygenic dyslipidemia. Nat. Genet. 2009, 41, 56–65. [Google Scholar] [CrossRef]
  50. Sabatti, C.; Service, S.K.; Hartikainen, A.-L.; Pouta, A.; Ripatti, S.; Brodsky, J.; Jones, C.G.; Zaitlen, N.A.; Varilo, T.; Kaakinen, M.; et al. Genome-wide association analysis of metabolic traits in a birth cohort from a founder population. Nat. Genet. 2009, 41, 35–46. [Google Scholar] [CrossRef]
  51. DiStefano, J.K.; Kingsley, C.; Wood, G.C.; Chu, X.; Argyropoulos, G.; Still, C.D.; Doné, S.C.; Legendre, C.; Tembe, W.; Gerhard, G.S. Genome-wide analysis of hepatic lipid content in extreme obesity. Acta Diabetol. 2015, 52, 373–382. [Google Scholar] [CrossRef] [PubMed]
  52. Medvedeva, Y.; Fridman, M.V.; Oparina, N.J.; Malko, D.B.; Ermakova, E.O.; Kulakovskiy, I.V.; Heinzel, A.; Makeev, V.J. Intergenic, gene terminal, and intragenic CpG islands in the human genome. BMC Genom. 2010, 11, 48. [Google Scholar] [CrossRef] [PubMed]
  53. Pérez, R.F.; Santamarina, P.; Tejedor, J.R.; Urdinguio, R.G.; Álvarez-Pitti, J.; Redon, P.; Fernández, A.F.; Fraga, M.F.; Lurbe, E. Longitudinal genome-wide DNA methylation analysis uncovers persistent early-life DNA methylation changes. J. Transl. Med. 2019, 17, 15. [Google Scholar] [CrossRef] [PubMed]
  54. Rodriguez, J.; Frigola, J.; Vendrell, E.; Risques, R.-A.; Fraga, M.F.; Morales, C.; Moreno, V.; Esteller, M.; Capellà, G.; Ribas, M.; et al. Chromosomal instability correlates with genome-wide DNA demethylation in human primary colorectal cancers. Cancer Res. 2006, 66, 8462–9468. [Google Scholar] [CrossRef] [PubMed]
  55. Sandoval, J.; Heyn, H.A.; Serra-Musach, J.; Bibikova, M.; Moran, S.; Pujana, M.A.; Esteller, M. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics 2011, 6, 692–702. [Google Scholar] [CrossRef] [PubMed]
  56. Timp, W.; Bravo, H.C.; McDonald, O.G.; Goggins, M.; Umbricht, C.; Zeiger, M.; Feinberg, A.P.; Irizarry, R.A. Large hypomethylated blocks as a universal defining epigenetic alteration in human solid tumors. Genome Med. 2014, 6, 61. [Google Scholar] [CrossRef] [PubMed]
  57. Acevedo, N.; Reinius, L.E.; Vitezic, M.; Fortino, V.; Söderhäll, C.; Honkanen, H.; Veijola, R.; Simell, O.; Toppari, J.; Ilonen, J.; et al. Age-associated DNA methylation changes in immune genes, histone modifiers and chromatin remodeling factors within 5 years after birth in human blood leukocytes. Clin. Epigenet. 2015, 7, 34. [Google Scholar] [CrossRef]
  58. Alisch, R.S.; Barwick, B.G.; Chopra, P.; Myrick, L.K.; Satten, G.A.; Conneely, K.N.; Warren, S.T. Age-associated DNA methylation in pediatric populations. Genome Res. 2012, 22, 623–632. [Google Scholar] [CrossRef] [Green Version]
  59. Christensen, B.C.; Houseman, E.A.; Marsit, C.J.; Zheng, S.; Wrensch, M.R.; Wiemels, J.L.; Nelson, H.H.; Karagas, M.R.; Padbury, J.F.; Bueno, R.; et al. Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet. 2009, 5, e1000602. [Google Scholar] [CrossRef]
  60. Lurbe, E.; Urdinguio, R.G.; Torró, M.I.; Bayón, G.F.; Álvarez-Pitti, J.; Fernández, A.F.; Redón, P.; Fraga, M.F. Longitudinal study of DNA methylation during the first 5 years of life. J. Transl. Med. 2016, 14, 160. [Google Scholar]
  61. Day, K.; Waite, L.L.; Thalacker-Mercer, A.; West, A.; Bamman, M.M.; Brooks, J.D.; Myers, R.M.; Absher, D. Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape. Genome Biol. 2013, 14, R102. [Google Scholar] [CrossRef] [PubMed]
  62. Gelfman, S.; Ast, G. When epigenetics meets alternative splicing: The roles of DNA methylation and GC architecture. Epigenomics 2013, 5, 351–353. [Google Scholar] [CrossRef] [PubMed]
  63. Gelfman, S.; Cohen, N.; Yearim, A.; Ast, G. DNA-methylation effect on cotranscriptional splicing is dependent on GC architecture of the exon-intron structure. Genome Res. 2013, 23, 789–799. [Google Scholar] [CrossRef] [PubMed]
  64. Hoivik, E.A.; Bjanesoy, T.E.; Mai, O.; Okamoto, S.; Minokoshi, Y.; Shima, Y.; Morohashi, K.-I.; Boehm, U.; Bakke, M. DNA methylation of intronic enhancers directs tissue-specific expression of steroidogenic factor 1/adrenal 4 binding protein (SF-1/Ad4BP). Endocrinology 2011, 152, 2100–2112. [Google Scholar] [CrossRef] [PubMed]
  65. Unoki, M.; Nakamura, Y. Methylation at CpG islands in intron 1 of EGR2 confers enhancer-like activity. FEBS Lett. 2003, 554, 67–72. [Google Scholar] [CrossRef]
  66. Portella, G.; Battistini, F.; Orozco, M. Understanding the connection between epigenetic DNA methylation and nucleosome positioning from computer simulations. PLoS Comput. Biol. 2013, 9, e1003354. [Google Scholar] [CrossRef] [PubMed]
  67. Liu, Y.; Liang, G.; Kelly, T.K.; Lay, F.D.; Berman, B.P.; Jones, P.A. Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules. Genome Res. 2012, 22, 2497–2506. [Google Scholar] [Green Version]
  68. Chodavarapu, R.K.; Feng, S.; Bernatavichute, Y.V.; Chen, P.-Y.; Stroud, H.; Yu, Y.; Hetzel, J.A.; Kuo, F.; Kim, J.; Cokus, S.J.; et al. Relationship between nucleosome positioning and DNA methylation. Nature 2010, 466, 388–392. [Google Scholar] [CrossRef] [Green Version]
  69. Pennings, S.; Allan, J.; Davey, C.S. DNA methylation, nucleosome formation and positioning. Brief. Funct. Genom. 2005, 3, 351–361. [Google Scholar] [CrossRef]
  70. Jones, P.A. Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 2012, 13, 484–492. [Google Scholar] [CrossRef]
  71. Rauch, T.A.; Pfeifer, G.P. The MIRA method for DNA methylation analysis. Methods Mol. Biol. 2009, 507, 65–75. [Google Scholar] [PubMed]
  72. Yang, Y.; Yang, C.; Zhu, Y.; Chen, H.; Zhao, R.; He, X.; Tao, L.; Wang, P.; Zhou, L.; Zhao, L.; et al. Intragenic and extragenic disruptions of FOXL2 mapped by whole genome low-coverage sequencing in two BPES families with chromosome reciprocal translocation. Genomics 2014, 104, 170–176. [Google Scholar] [CrossRef] [PubMed]
  73. Lister, R.; Pelizzola, M.; Dowen, R.H.; Hawkins, R.D.; Hon, G.; Tonti-Filippini, J.; Nery, J.R.; Lee, L.; Ye, Z.; Ngo, Q.-M.; et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 2009, 462, 315–322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Brenet, F.; Moh, M.; Funk, P.; Feierstein, E.; Viale, A.J.; Socci, N.D.; Scandura, J.M. DNA methylation of the first exon is tightly linked to transcriptional silencing. PLoS ONE 2011, 6, e14524. [Google Scholar] [CrossRef] [PubMed]
  75. Burdge, G.C.; Hanson, M.A.; Slater-Jefferies, J.L.; Lillycrop, K.A. Epigenetic regulation of transcription: A mechanism for inducing variations in phenotype (fetal programming) by differences in nutrition during early life? Br. J. Nutr. 2007, 97, 1036–1046. [Google Scholar] [CrossRef] [PubMed]
  76. Waldecker, M.; Kautenburger, T.; Daumann, H.; Busch, C.; Schrenk, D. Inhibition of histone-deacetylase activity by short-chain fatty acids and some polyphenol metabolites formed in the colon. J. Nutr. Biochem. 2008, 19, 587–593. [Google Scholar] [CrossRef] [PubMed]
  77. Takahashi, H.; McCaffery, J.M.; Irizarry, R.A.; Boeke, J.D. Nucleocytosolic acetyl-coenzyme a synthetase is required for histone acetylation and global transcription. Mol. Cell 2006, 23, 207–217. [Google Scholar] [CrossRef] [PubMed]
  78. Pradhan, S.; Bacolla, A.; Wells, R.D.; Roberts, R. Recombinant human DNA (cytosine-5) methyltransferase. I. Expression, purification, and comparison of de novo and maintenance methylation. J. Biol. Chem. 1999, 274, 33002–33010. [Google Scholar] [CrossRef]
  79. Pogribny, I.P.; Tryndyak, V.P.; Bagnyukova, T.V.; Melnyk, S.; Montgomery, B.; Ross, S.A.; Latendresse, J.R.; Rusyn, I.; Beland, F.A. Hepatic epigenetic phenotype predetermines individual susceptibility to hepatic steatosis in mice fed a lipogenic methyl-deficient diet. J. Hepatol. 2009, 51, 176–186. [Google Scholar] [CrossRef] [Green Version]
  80. Ding, Y.B.; He, J.L.; Liu, X.Q.; Chen, X.M.; Long, C.L.; Wang, Y.X. Expression of DNA methyltransferases in the mouse uterus during early pregnancy and susceptibility to dietary folate deficiency. Reproduction 2012, 144, 91–100. [Google Scholar] [CrossRef]
  81. Ghoshal, K.; Li, X.; Datta, J.; Bai, S.; Pogribny, I.; Pogribny, M.; Huang, Y.; Young, D.; Jacob, S.T. A folate- and methyl-deficient diet alters the expression of DNA methyltransferases and methyl CpG binding proteins involved in epigenetic gene silencing in livers of F344 rats. J. Nutr. 2006, 136, 1522–1527. [Google Scholar] [CrossRef] [PubMed]
  82. Kovacheva, V.P.; Mellott, T.J.; Davison, J.M.; Wagner, N.; Lopez-Coviella, I.; Schnitzler, A.C.; Blusztajn, J.K. Gestational choline deficiency causes global and Igf2 gene DNA hypermethylation by up-regulation of Dnmt1 expression. J. Biol. Chem. 2007, 282, 31777–31788. [Google Scholar] [CrossRef]
  83. Xia, L.; Wang, C.; Lu, Y.; Fan, C.; Ding, X.; Fu, H.; Qi, K. Time-specific changes in DNA methyltransferases associated with the leptin promoter during the development of obesity. Nutr. Hosp. 2014, 30, 1248–1255. [Google Scholar] [PubMed]
  84. Strakovsky, R.S.; Wang, H.; Engeseth, N.J.; Flaws, J.A.; Helferich, W.G.; Pan, Y.-X.; Lezmi, S. Developmental bisphenol A (BPA) exposure leads to sex-specific modification of hepatic gene expression and epigenome at birth that may exacerbate high-fat diet-induced hepatic steatosis. Toxicol. Appl. Pharmacol. 2015, 284, 101–112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Zhang, X.; Zhou, D.; Pan, Y.-X.; Strakovsky, R.S.; Pan, Y. Gestational high fat diet programs hepatic phosphoenolpyruvate carboxykinase gene expression and histone modification in neonatal offspring rats. J. Physiol. 2011, 589, 2707–2717. [Google Scholar]
  86. Zhou, D.; Wang, H.; Cui, H.; Chen, H.; Pan, Y.-X. Early-life exposure to high-fat diet may predispose rats to gender-specific hepatic fat accumulation by programming Pepck expression. J. Nutr. Biochem. 2015, 26, 433–440. [Google Scholar] [CrossRef] [PubMed]
  87. Zhang, X.; Strakovsky, R.; Zhou, D.; Zhang, Y.; Pan, Y.-X. A maternal high-fat diet represses the expression of antioxidant defense genes and induces the cellular senescence pathway in the liver of male offspring rats. J. Nutr. 2011, 141, 1254–1259. [Google Scholar] [CrossRef]
  88. Zhang, Y.; Wang, H.; Zhou, D.; Moody, L.; Lezmi, S.; Chen, H.; Pan, Y.-X. High-fat diet caused widespread epigenomic differences on hepatic methylome in rat. Physiol. Genom. 2015, 47, 514–523. [Google Scholar] [CrossRef]
  89. Li, D.; Zhang, B.; Xing, X.; Wang, T. Combining MeDIP-seq and MRE-seq to investigate genome-wide CpG methylation. Methods 2015, 72, 29–40. [Google Scholar] [CrossRef]
  90. Claussen, U.; Michel, S.; Boehringer, S.; Gillessen-Kaesbach, G.; Buiting, K.; Horsthemke, B.; Nazlican, H.; Zeschnigk, M. Somatic mosaicism in patients with Angelman syndrome and an imprinting defect. Hum. Mol. Genet. 2004, 13, 2547–2555. [Google Scholar]
  91. Suzuki, T.; Higgins, P.J.; Crawford, D.R. Control Selection for RNA Quantitation. Biotechniques 2000, 29, 332–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Wang, H.; Wilson, G.J.; Zhou, D.; Lezmi, S.; Chen, X.; Layman, D.K.; Pan, Y.-X. Induction of autophagy through the activating transcription factor 4 (ATF4)-dependent amino acid response pathway in maternal skeletal muscle may function as the molecular memory in response to gestational protein restriction to alert offspring to maternal nutrition. Br. J. Nutr. 2015, 114, 519–532. [Google Scholar] [PubMed] [Green Version]
  93. Bass, J.I.F.; Diallo, A.; Nelson, J.; Soto, J.M.; Myers, C.L.; Walhout, A.J.M. Using networks to measure similarity between genes: Association index selection. Nat. Methods 2013, 10, 1169–1176. [Google Scholar] [CrossRef] [PubMed]
  94. Dennis, G.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4, R60. [Google Scholar] [CrossRef]
Figure 1. Physiological differences between perinatal and post-weaning HF diet exposure. (a) The experimental design involved pregnant dams that were fed either a HF diet (45% kcal from fat) or a standard AIN93G control diet throughout gestation and lactation. At weaning, male pups from control diet-fed dams were weaned onto either the same control diet or HF diet, and male pups from HF-fed dams were weaned onto the control diet, creating three groups of pups, control (CON; n = 8), post-weaning high fat (PHF; n = 9) and maternal high fat (MHF; n = 7), respectively. (b) Body weight was measured weekly during the post-weaning period. (c) Weekly food intake was monitored during the post-natal period. (d) Caloric intake as kcal/week was calculated based on the grams of food intake. (e) Liver cross-sections were stained with hematoxylin and eosin (H&E) (top) and Oil Red O (ORO; bottom). (f) PHF animals had greater hepatic lipid accumulation than MHF and CON rats. Lipid accumulation was normalized to total protein. Data points represent mean ± standard error of the mean (SEM). * denotes a significant difference (p < 0.05) between PHF and CON. ** denotes a significant difference (p < 0.01) between PHF and CON as well as between PHF and MHF.
Figure 1. Physiological differences between perinatal and post-weaning HF diet exposure. (a) The experimental design involved pregnant dams that were fed either a HF diet (45% kcal from fat) or a standard AIN93G control diet throughout gestation and lactation. At weaning, male pups from control diet-fed dams were weaned onto either the same control diet or HF diet, and male pups from HF-fed dams were weaned onto the control diet, creating three groups of pups, control (CON; n = 8), post-weaning high fat (PHF; n = 9) and maternal high fat (MHF; n = 7), respectively. (b) Body weight was measured weekly during the post-weaning period. (c) Weekly food intake was monitored during the post-natal period. (d) Caloric intake as kcal/week was calculated based on the grams of food intake. (e) Liver cross-sections were stained with hematoxylin and eosin (H&E) (top) and Oil Red O (ORO; bottom). (f) PHF animals had greater hepatic lipid accumulation than MHF and CON rats. Lipid accumulation was normalized to total protein. Data points represent mean ± standard error of the mean (SEM). * denotes a significant difference (p < 0.05) between PHF and CON. ** denotes a significant difference (p < 0.01) between PHF and CON as well as between PHF and MHF.
Ijms 20 03229 g001
Figure 2. Identification of differentially methylated regions in maternal and post-weaning HF diet groups. (a) Heatmap shows the weighted MeDIP and MRE reads in the MHF and PHF conditions (black and white), where darker shades represent more reads. Overall methylation is shown in the far-right column. Blue represents higher methylation in the MHF group compared to the PHF group and red represents higher in the PHF group compared to the MHF group (FDR q-value < 0.05). (b) Pie chart shows overall numbers of DMRs that are more methylated in MHF animals (blue, n = 754) and that are more methylated in PHF animals (red, n = 990). (cf) Probe-based analysis of DNA methylation was used to validate results from MeDIP and MRE-seq. (c) Representative tracks from the Genome Browser showing a DMR within the Pklr gene that was identified by MeDIP-seq and MRE-seq analysis to be more methylated in the PHF group. The region of interest used for validation is outlined in orange. The top tracks show the MRE reads in black. The red track indicates the MeDIP reads for each group. The blue track represents annotated genes and indicates the DMR position relative to exons (blue bar) and introns (blue hashed line). (d) A probe-based analysis of DNA methylation was used to validate the DMR. The PHF group had a higher percentage of methylated DNA (% of amplified DNA) and a lower percentage of unmethylated DNA compared to the MHF group (* p < 0.00001). (e) Representative tracks from the Genome Browser showing a genomic region upstream of the Il6 gene identified by the joint MeDIP-seq and MRE-seq analysis as being not differentially methylated in the computational analysis. (f) Methylation in the Il6 locus was measured using probe-based analysis of DNA methylation. No difference was detected between the MHF and PHF groups in this region (p = 0.2). (g) MSP was used to confirm methylation status of a region in the Acacb gene identified as significantly differentially methylated by MeDIP-seq and MRE-seq (p = 0.007). (h) MSP was used to validate no change in methylation status of a DMR within the Gpam gene that was not identified to be significant by MeDIP-seq and MRE-seq (p = 0.7). All bar graphs show mean ± SEM.
Figure 2. Identification of differentially methylated regions in maternal and post-weaning HF diet groups. (a) Heatmap shows the weighted MeDIP and MRE reads in the MHF and PHF conditions (black and white), where darker shades represent more reads. Overall methylation is shown in the far-right column. Blue represents higher methylation in the MHF group compared to the PHF group and red represents higher in the PHF group compared to the MHF group (FDR q-value < 0.05). (b) Pie chart shows overall numbers of DMRs that are more methylated in MHF animals (blue, n = 754) and that are more methylated in PHF animals (red, n = 990). (cf) Probe-based analysis of DNA methylation was used to validate results from MeDIP and MRE-seq. (c) Representative tracks from the Genome Browser showing a DMR within the Pklr gene that was identified by MeDIP-seq and MRE-seq analysis to be more methylated in the PHF group. The region of interest used for validation is outlined in orange. The top tracks show the MRE reads in black. The red track indicates the MeDIP reads for each group. The blue track represents annotated genes and indicates the DMR position relative to exons (blue bar) and introns (blue hashed line). (d) A probe-based analysis of DNA methylation was used to validate the DMR. The PHF group had a higher percentage of methylated DNA (% of amplified DNA) and a lower percentage of unmethylated DNA compared to the MHF group (* p < 0.00001). (e) Representative tracks from the Genome Browser showing a genomic region upstream of the Il6 gene identified by the joint MeDIP-seq and MRE-seq analysis as being not differentially methylated in the computational analysis. (f) Methylation in the Il6 locus was measured using probe-based analysis of DNA methylation. No difference was detected between the MHF and PHF groups in this region (p = 0.2). (g) MSP was used to confirm methylation status of a region in the Acacb gene identified as significantly differentially methylated by MeDIP-seq and MRE-seq (p = 0.007). (h) MSP was used to validate no change in methylation status of a DMR within the Gpam gene that was not identified to be significant by MeDIP-seq and MRE-seq (p = 0.7). All bar graphs show mean ± SEM.
Ijms 20 03229 g002
Figure 3. DMR distributions in relevance to CpG islands and to gene locations and structures. (a) DMR distribution around CpG islands. CpG shore was defined as the 2000 base pair sequence flanking either side of a CpG island. DMRs in the shores (dark gray, n = 484) outnumbered the DMRs located within CpG islands (black, n = 227). (b) DMR locations relative to gene structures. Promoter and downstream regions were defined as ± 1500 bp upstream of the transcription start site and downstream of the transcription end site, respectively. Top panel depicts the number of DMRs not associated with a gene (blue, n = 1083) and DMRs located within the gene body (red, n = 557), within the promoter (gray, n = 43), and within the downstream region (yellow, n = 61). A majority of DMRs are not associated with protein-coding genes. Bottom panels classify the DMRs within the promoter, gene body, and downstream regions. DMRs that are more methylated in MHF (MHF > PHF) are in darker shades, while those that more methylated in PHF (PHF > MHF) are in lighter shades.
Figure 3. DMR distributions in relevance to CpG islands and to gene locations and structures. (a) DMR distribution around CpG islands. CpG shore was defined as the 2000 base pair sequence flanking either side of a CpG island. DMRs in the shores (dark gray, n = 484) outnumbered the DMRs located within CpG islands (black, n = 227). (b) DMR locations relative to gene structures. Promoter and downstream regions were defined as ± 1500 bp upstream of the transcription start site and downstream of the transcription end site, respectively. Top panel depicts the number of DMRs not associated with a gene (blue, n = 1083) and DMRs located within the gene body (red, n = 557), within the promoter (gray, n = 43), and within the downstream region (yellow, n = 61). A majority of DMRs are not associated with protein-coding genes. Bottom panels classify the DMRs within the promoter, gene body, and downstream regions. DMRs that are more methylated in MHF (MHF > PHF) are in darker shades, while those that more methylated in PHF (PHF > MHF) are in lighter shades.
Ijms 20 03229 g003
Figure 4. Summary of epigenetic and gene expression differences in the insulin and PI signaling pathways. Genes in each pathway are boxed and filled with two colors. The color in the left side of the box indicates methylation levels, while the color in the right side of the box denotes gene expression levels. White color signifies no change, blue color signifies higher methylation or expression in MHF animals, and red signifies higher methylation or expression in PHF animals. Yellow and purple boxes represent metabolic functions. A solid line represents a direct relationship, a dashed line represents a multi-step process, and compounds with no outline are other important signaling molecules.
Figure 4. Summary of epigenetic and gene expression differences in the insulin and PI signaling pathways. Genes in each pathway are boxed and filled with two colors. The color in the left side of the box indicates methylation levels, while the color in the right side of the box denotes gene expression levels. White color signifies no change, blue color signifies higher methylation or expression in MHF animals, and red signifies higher methylation or expression in PHF animals. Yellow and purple boxes represent metabolic functions. A solid line represents a direct relationship, a dashed line represents a multi-step process, and compounds with no outline are other important signaling molecules.
Ijms 20 03229 g004
Figure 5. Differential gene expression in the insulin and PI signaling pathways. (a) Methylated DNA reads were measured using MeDIP-seq. † denotes a ≥15% difference between CON and either PHF or MHF. Data are presented as mean ±SEM across each DMR. (b) Unmethylated DNA was measured using MRE-seq. † denotes a ≥15% difference between CON and either PHF or MHF. Data are presented as mean ±SEM across each DMR. (c) Eight DMGs were significantly differentially expressed between MHF and PHF groups. The bar graph represents the mRNA quantity relative to B Actin. Data are presented as mean ±SEM. * p < 0.05.
Figure 5. Differential gene expression in the insulin and PI signaling pathways. (a) Methylated DNA reads were measured using MeDIP-seq. † denotes a ≥15% difference between CON and either PHF or MHF. Data are presented as mean ±SEM across each DMR. (b) Unmethylated DNA was measured using MRE-seq. † denotes a ≥15% difference between CON and either PHF or MHF. Data are presented as mean ±SEM across each DMR. (c) Eight DMGs were significantly differentially expressed between MHF and PHF groups. The bar graph represents the mRNA quantity relative to B Actin. Data are presented as mean ±SEM. * p < 0.05.
Ijms 20 03229 g005
Figure 6. Genomic context of DMRs. (a) DMGs within the PI and insulin signaling pathways were mapped according to their distance from the nearest CpG island. Gene names in red and blue denote those that were more highly expressed in the PHF and MHF group, respectively. (b) Intragenic DMR position is shown for each DMG in the PI and insulin signaling pathways. Blue circles represent DMRs in which DNA methylation was higher in the MHF group. Red circles represent DMRs in which DNA methylation was higher in the PHF group. Gene names in red and blue denote those that were more highly expressed in the PHF and MHF group, respectively.
Figure 6. Genomic context of DMRs. (a) DMGs within the PI and insulin signaling pathways were mapped according to their distance from the nearest CpG island. Gene names in red and blue denote those that were more highly expressed in the PHF and MHF group, respectively. (b) Intragenic DMR position is shown for each DMG in the PI and insulin signaling pathways. Blue circles represent DMRs in which DNA methylation was higher in the MHF group. Red circles represent DMRs in which DNA methylation was higher in the PHF group. Gene names in red and blue denote those that were more highly expressed in the PHF and MHF group, respectively.
Ijms 20 03229 g006
Figure 7. Gene expression is associated with intragenic DMR position. X-axis indicates intragenic DMR position relative to the gene’s transcription start site (x = 0) and transcription end site (x = 100). Y-axis represents the density of DMRs. (a) the density plot depicts DMGs that were more highly expressed in the MHF group. Higher gene expression in MHF was associated with lower DNA methylation in the middle of the intragenic region. (b) The density plot shows DMGs that were more highly expressed in the PHF group. Higher PHF expression was associated with higher downstream DNA methylation and lower DNA methylation near the TSS.
Figure 7. Gene expression is associated with intragenic DMR position. X-axis indicates intragenic DMR position relative to the gene’s transcription start site (x = 0) and transcription end site (x = 100). Y-axis represents the density of DMRs. (a) the density plot depicts DMGs that were more highly expressed in the MHF group. Higher gene expression in MHF was associated with lower DNA methylation in the middle of the intragenic region. (b) The density plot shows DMGs that were more highly expressed in the PHF group. Higher PHF expression was associated with higher downstream DNA methylation and lower DNA methylation near the TSS.
Ijms 20 03229 g007
Table 1. QAMA Primers and Probes.
Table 1. QAMA Primers and Probes.
GenePositionSequence
Pklr (+ chr2)Forward Primer + 43335′-TGGTGTTATTTAGATGTTGGAGAGTATGA-3′
Reverse Primer + 45575′-AACATAATACAATCAACCCCATCCA-3′
Methylated Probe + 44675′-VIC-AGGTTCGATTAATTCGGGCG-MGB-3′
Unmethylated Probe + 44705′-FAM-TGATTAATTTGGGTGGAGATAA-MGB-3′
IL-6Forward Primer − 2516GTGAGTAAGGGATTTAGTTTGAGTATGGT
Reverse Primer − 2441CTTATTCCTAAATATCTAATACCCTCTTATAACCTC
Methylated Probe − 24805′-VIC-CGTGTGTGAATGTGCGTTA-3′
Unmethylated Probe − 24835′-FAM-TTTGTGTGTGAATGTGTGTT-MGB-3′
Table 2. Pathway Analysis.
Table 2. Pathway Analysis.
PathwayDMGsFold Enrichmentp-Value
Phosphatidylinositol signaling system115.20.00004
Insulin signaling pathway112.80.006
Axon guidance102.60.01
Pathways in cancer181.90.01
Purine metabolism102.10.04
Glycerolipid metabolism53.80.04
Fc gamma R-mediated phagocytosis72.60.05
Maternal versus post-weaning HF diet exposure results in seven differentially methylated pathways. DMGs column refers to the number of differentially methylated genes in the pathway. Fold Enrichment signifies the degree of enrichment within each pathway given the overall list of DMGs. A lower p-value indicates higher gene-enrichment within a pathway.
Table 3. MSP Primers.
Table 3. MSP Primers.
GenePositionSequence
AcacbForward M + 75,0965′-TTGGGTTCGGTTTTTAGTTTCG-3′
Reverse M + 75,2225′-ACGTATATCCCTATAATCCAACTCGC-3′
Forward UM + 75,0955′-TTTGGGTTTGGTTTTTAGTTTTGAA-3′
Reverse UM + 75,2255′-CACATATATCCCTATAATCCAACTCACTCT -3′
GpamForward M + 72895′-AGTCGTAGTGGTCGGGTAATCG-3′
Reverse M + 73575′-CCGCTTATTTTAAACAACATCGAA-3′
Forward UM + 72885′-AAGTTAAGTTGTAGTGGTTGGGTAATTG-3′
Reverse UM + 73605′-CCCACTTATTTTAAACAACATCAAACC-3′

Share and Cite

MDPI and ACS Style

Moody, L.; Wang, H.; Jung, P.M.; Chen, H.; Pan, Y.-X. Maternal and Post-Weaning High-Fat Diets Produce Distinct DNA Methylation Patterns in Hepatic Metabolic Pathways within Specific Genomic Contexts. Int. J. Mol. Sci. 2019, 20, 3229. https://doi.org/10.3390/ijms20133229

AMA Style

Moody L, Wang H, Jung PM, Chen H, Pan Y-X. Maternal and Post-Weaning High-Fat Diets Produce Distinct DNA Methylation Patterns in Hepatic Metabolic Pathways within Specific Genomic Contexts. International Journal of Molecular Sciences. 2019; 20(13):3229. https://doi.org/10.3390/ijms20133229

Chicago/Turabian Style

Moody, Laura, Huan Wang, Paul M. Jung, Hong Chen, and Yuan-Xiang Pan. 2019. "Maternal and Post-Weaning High-Fat Diets Produce Distinct DNA Methylation Patterns in Hepatic Metabolic Pathways within Specific Genomic Contexts" International Journal of Molecular Sciences 20, no. 13: 3229. https://doi.org/10.3390/ijms20133229

APA Style

Moody, L., Wang, H., Jung, P. M., Chen, H., & Pan, Y. -X. (2019). Maternal and Post-Weaning High-Fat Diets Produce Distinct DNA Methylation Patterns in Hepatic Metabolic Pathways within Specific Genomic Contexts. International Journal of Molecular Sciences, 20(13), 3229. https://doi.org/10.3390/ijms20133229

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop