# Maternal high-protein and low-protein diets perturb hypothalamus and liver transcriptome and metabolic homeostasis in adult mouse offspring

Lisa J. Martin, Qingying Meng, Montgomery Blencowe, Sandrine Lagarrigue, Sheila Xiao, Calvin Pan, Julian Wier, William C. Temple, Sherin U. Devaskar, Aldons J. Lusis, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Lisa J. Martin, Qingying Meng, Montgomery Blencowe, Sandrine Lagarrigue, Sheila Xiao, et al., Maternal high-protein and low-protein diets perturb hypothalamus and liver transcriptome and metabolic homeostasis in adult mouse offspring. *Frontiers in Genetics*, 2018, 9, 10.3389/fgene.2018.00642. hal-01977616

HAL Id: hal-01977616

https://hal.science/hal-01977616

Submitted on 10 Jan 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for
updates

# Maternal High-Protein and Low-Protein Diets Perturb Hypothalamus and Liver Transcriptome and Metabolic Homeostasis in Adult Mouse Offspring

Lisa J. Martin<fnref n="1" />, Qingying Meng², Montgomery Blencowe², Sandrine Lagarrigue³, Sheila Xiao<fnref n="1" />, Calvin Pan<fnref n="1" />, Julian Wier², William C. Temple<fnref n="1" />, Sherin U. Devaskar⁴, Aldons J. Lusis<fnref n="1" />,<fnref n="5" />,⁶* and Xia Yang²*

OPEN ACCESS

**Edited by:**
Katia Aquilano,
Università degli Studi di Roma Tor
Vergata, Italy

**Reviewed by:**
Simone Carotti,
Università Campus Bio-Medico, Italy
Wimal Pathmasiri,
University of North Carolina at
Chapel Hill, United States

*Correspondence:
Aldons J. Lusis
jlusis@mednet.ucla.edu
Xia Yang
xyang123@ucla.edu

**Specialty section:**
This article was submitted to
Nutrigenomics,
a section of the journal
Frontiers in Genetics

**Received:** 09 October 2018
**Accepted:** 27 November 2018
**Published:** 11 December 2018

**Citation:**
Martin LJ, Meng Q, Blencowe M, Lagarrigue S, Xiao S, Pan C, Wier J, Temple WC, Devaskar SU, Lusis AJ and Yang X (2018) Maternal High-Protein and Low-Protein Diets Perturb Hypothalamus and Liver Transcriptome and Metabolic Homeostasis in Adult Mouse Offspring. Front. Genet. 9:642. doi: 10.3389/fgene.2018.00642

Early life nutritional imbalances are risk factors for metabolic dysfunctions in adulthood, but the long term effects of perinatal exposure to high versus low protein diets are not completely understood. We exposed C57BL/6J offspring to a high protein/low carbohydrate (HP/LC) or low protein/high carbohydrate (LP/HC) diet during gestation and lactation, and measured metabolic phenotypes between birth and 10 months of age in male offspring. Perinatal HP/LC and LP/HC exposures resulted in a decreased ability to clear glucose in the offspring, with reduced baseline insulin and glucose concentrations in the LP/HC group and a reduced insulin response post-glucose challenge in the HP/LC group. The LP/HC diet group also showed reduced birth and weanling weights, whereas the HP/LC offspring displayed increased weanling weight with increased adiposity beyond 5 months of age. Gene expression profiling of hypothalamus and liver revealed alterations in diverse molecular pathways by both diets. Specifically, hypothalamic transcriptome and pathway analyses demonstrated perturbations of MAPK and hedgehog signaling, processes associated with neural restructuring and transmission, and phosphate metabolism by perinatal protein imbalances. Liver transcriptomics revealed changes in purine and phosphate metabolism, hedgehog signaling, and circadian rhythm pathways. Our results indicate maternal protein imbalances perturbing molecular pathways in central and peripheral metabolic tissues, thereby predisposing the male offspring to metabolic dysfunctions.

**Keywords:** maternal diet, high protein, low protein, gene expression, metabolic dysfunction, glucose intolerance, hypothalamus, liver

<footnote n="1">¹ Department of Medicine, Division of Cardiology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, United States, ² Department of Integrative Biology and Physiology, University of California, Los Angeles, Los Angeles, CA, United States, ³ INRA, Agrocampus Ouest, UMR1348 PEGASE, Rennes, France, ⁴ Department of Pediatrics, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, United States,</footnote>

<footnote n="5">⁵ Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, United States, ⁶ Department of Microbiology, Immunology and Molecular Genetics, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA, United States</footnote>
# INTRODUCTION

Metabolic syndrome and type 2 diabetes mellitus (T2DM) have reached epidemic proportions worldwide (Saklayen, 2018). A major risk factor for the development of T2DM is the quantity and quality of intra-uterine nutrient exposure (Gluckman et al., 2008). Studies have shown that the maternal diet can modify the development of regulatory systems in utero and postnatally. Individual studies of high and low protein maternal diets have shown their influence on body weight, metabolic phenotypes and programming of food intake in offspring (Oken and Gillman, 2003; Frias and Grove, 2012). Our previous mouse study on intrauterine growth restriction (IUGR) induced by a low protein maternal diet demonstrated sex-specific phenotypic differences in offspring at 9 months of age (Bhasin et al., 2009). While the male offspring displayed glucose intolerance and increased adiposity, the female offspring were glucose tolerant, without catch-up growth or adiposity. However, the mechanisms underlying these metabolic changes are yet to be elucidated. In addition, previous studies have not systematically compared the effect of low protein/high carbohydrate (LP/HC) versus high protein/low carbohydrate (HP/LC) maternal diets on modulating tissue-specific gene expression profiles in key metabolic tissues and the associated pathological metabolic variations in offspring. As gene expression differences involved in metabolically relevant pathways in the hypothalamus (the control center of energy balance and metabolism) and liver (critical for lipid and glucose homeostasis) may underlie the differing metabolic phenotypes in adulthood, we aim to investigate how early-life nutritional imbalances contribute to metabolic syndrome outcomes in adulthood.

In the current study, we tested the hypothesis that intrauterine HP/LC and LP/HC diets would subsequently affect adult metabolic phenotypes and long-term gene expression levels in key metabolic sites. Using HP/LC or LP/HC dietary interventions during gestation and lactation, we followed male offspring until 10 months of age and assessed body weight, adiposity, insulin production, glucose clearance, and global gene expression changes in the hypothalamus and liver. We discovered that mice expressed differential programming by perinatal exposure to protein imbalances, and these changes persist into adulthood. Importantly, we found that gene expression changes brought about by perinatal nutritional perturbations revealed enrichment of genes involved in diverse pathways including those important for systemic energy homeostasis.

# MATERIALS AND METHODS

## Study Approval

This study was approved by the UCLA Animal Research Committee and was performed in accordance with the National Institutes of Health guidelines for the use of experimental animals. C57BL/6J mice were purchased from the Jackson Laboratory (Bar Harbor, ME, United States).

## Diets

The LP/HC diet (D02041002, Research Diets Inc., New Brunswick, NJ, United States) contains 9% protein, 4.4% fat, and 77% carbohydrates. The HP/LC diet (D02041001, Research Diets Inc., New Brunswick, NJ, United States) has 23% protein, 4.4% fat and 64% carbohydrate by weight. Standard chow diet (TD 7013, Harlan Teklad, Placentia, CA, United States), containing 18% protein, 6% fat, and 45% carbohydrate by weight, was used as the control diet. The diets were matched for total caloric content, making them isocaloric.

## Maternal Diet Exposure Studies

As most first litters did not survive, all the offspring were from second litters from dams 10 to 12 weeks old. C57BL/6 females were mated overnight with males between 10 and 16 weeks of age. Gestational day 0 (GD0) was determined by the detection of a vaginal plug in the morning. On GD8, pregnant females were placed on a LP/HC, HP/LC, or chow diet. Pups were weaned at 28 days of age into cages of four animals per cage, separated by sex and maternal dietary environment. A total of n = 22, n = 27, and n = 20 male offspring from the LP/HC, HP/LC, or chow diet group, respectively, were examined.

## Body Weight and Body Composition

Body weight of individual mice was measured on a scale, starting day 2 (birth weight) and thereafter monthly until sacrifice. Body composition was measured monthly until sacrifice by a rodent Nuclear Magnetic Resonance (NMR) scanner (Bruker Biospin, Billerica, MA, United States) that was standardized to an internal control provided by the manufacturer. Adiposity was determined as adipose tissue mass per unit (g) body weight.

## Glucose Tolerance Tests and Insulin Measurement

Glucose tolerance tests were performed as previously described (Ganguly and Devaskar, 2008) in 9-month-old male offspring. After an overnight fast, glucose was measured at 9AM in blood collected from saphenous venous puncture, after which 2 mg/g body weight of glucose was administered intraperitoneally. Blood glucose was measured at half hour intervals for 2 h. The One Touch Ultra (LifeScan) glucometer was used to measure whole blood glucose concentrations (Ganguly and Devaskar, 2008). Insulin was measured with an ELISA kit (ALPCO, Salem, NH #80-INSMSU-E10) using 5 µL of plasma separated from blood collected at baseline and at 30 min post-glucose challenge.

## Phenotypic Data Analysis

All figures demonstrate data displayed as means ± SEM, with the exception of insulin concentrations that were expressed as a ratio of 30 min/baseline values in **Figure 2**, indicating medians. The two-way ANOVA model was employed for growth curves and glucose tolerance test curves, to simultaneously compare the effect of diets (treatment) and time points within each group. GraphPad Prism version 7.0 software was used for graphing and for statistical analysis of phenotypic data.
## Microarray Analysis

Hypothalamus and liver were dissected between 10 AM and 12 PM from male mice at 10 months of age, placed immediately into liquid nitrogen, and stored at -80°C until RNA was extracted with the RNAeasy kit (Qiagen). A total of 40 males from three different *in utero* exposure groups were analyzed for gene expression changes in hypothalamus and liver (15, 15, and 10 for LP/HC, HP/LC, and Control, respectively) using the Affymetrix Mouse 430 2.0 microarray chip. The chip contained 45,000 probes representing 35,000 genes. *P*-values for differential gene expression were determined using the Student's *t*-test for pairwise comparisons, each treatment compared with control, followed by multiple testing correction using the *q*-value approach (Storey, 2003). For pathway analysis, differentially expressed probes with *p* ≤ 0.01 were analyzed for KEGG pathway and Gene Ontology enrichment using DAVID (Huang da et al., 2009).

The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO series accession number GSE120548.

## RESULTS

### Effects of Maternal HP/LC and LP/HC Diets on Birth Weights, Perinatal Weights, and Adult Growth Curves

*In utero* LP/HC exposure resulted in IUGR as reflected by significantly lower birth weights, whereas HP/LC litter birth weights were not different from control mice (Figure 1A). At weaning, however, HP/LC offspring weighed significantly more than controls (Figure 1B), whereas IUGR LP/HC pups still weighed less at weaning than controls (Figure 1B). IUGR LP/HC weighed significantly less than controls throughout the study, but the higher body weights of HP/LC at weaning did not persist over time and after 3 months HP/LC body weights normalized and were not different from controls (Figure 1C).

### Effects of Maternal HP/LC and LP/HC Diets on Adiposity in Adulthood

In order to compare the impact of *in utero*/postnatal exposure on adult fat storage, we measured adiposity using NMR in the offspring at 5 months of age. HP/LC male offspring had significantly increased adiposity (Figure 1D), which supports previous findings (Thone-Reineke et al., 2006), while LP/HC offspring showed no difference from controls in adiposity (Figure 1D), also supporting previous findings (Daenzer et al., 2002).

### Effects of Maternal HP/LC and LP/HC Diets on Glucose Tolerance and Insulin Response

A key metabolic syndrome indicator is glucose intolerance, which is a measure of baseline glucose, baseline insulin, and the ability

to respond to an exogenous glucose load with adequate insulin production. We therefore carried out an Intraperitoneal Glucose Tolerance Test (IPGTT) at 9 months of age.

Early life exposure to LP/HC and HP/LC impaired the ability of offspring to clear glucose in adulthood. Both LP/HC and HP/LC mice had a significantly higher area under the curve (AUC) for blood glucose levels measured at 0.5, 1.0, and 1.5 h post glucose administration (Figure 2A). For LP/HC offspring, both baseline glucose and insulin concentrations were reduced when compared to control (Figures 2B,C), the latter of which likely contributed to the higher AUC for blood glucose (Figure 2D). In contrast, baseline glucose and insulin concentrations of HP/LC offspring were similar to controls; however, at 30 min post glucose injection, the insulin concentrations were significantly lower (Figure 2E), indicating an inability to mount an effective insulin response to the glucose challenge, and resulting in the higher AUC. We also examined the insulin response per individual offspring by using an insulin response ratio (IRR), which was the insulin concentration value at 30 min post glucose challenge divided by the baseline value (Figure 2F). This comparison provided further evidence that the higher AUC in HP/LC male offspring was likely due to an insufficient insulin response, with half of HP/LC animals having an IRR < 1 (lower insulin after the challenge than at baseline) (Figure 2F). Overall, the median IRR was 2.2 for control, 2.3 for LP/HC, and 1.0 for HP/LC (Figure 2F).

## Gene Expression and Pathway Alterations in the Liver and Hypothalamus

The metabolic phenotypic differences in adult offspring were accompanied by significant differences in gene expression in two metabolically relevant tissues, the liver and hypothalamus collected at final dissection conducted at 10 months of age. For the hypothalamus, we identified 1,441 and 843 differentially expressed genes (DEGs) for HP/LC and LP/HC groups, respectively, by Student's *t*-test using a *p* < 0.01 as the cutoff. For liver, 2014 and 926 DEGs were identified for HP/LC and LP/HC groups, respectively. At a false discovery rate (FDR) < 5%, very few DEGs were identified (2, 2, 14, and 2 genes altered in the LP/HC hypothalamus, HP/LC hypothalamus, LP/HC liver, and HP/LC liver, respectively; Supplementary Tables S5-S8).

We queried the DEGs reaching *p* < 0.01 against canonical pathways curated in KEGG and GO, and found significant over-representation of genes in diverse pathways ranging from metabolic, neural structure, protein modification, to signaling pathways such as MAPK (Table 1). Pathway analysis does not rely on the reliability of individual DEGs but the aggregate activities of multiple DEGs in coherent pathways, thus reducing false positive detection. Compared to the LP/HC group, the HP/LC group had more significant pathways that were expressed in both liver and hypothalamus at Bonferroni-corrected *p* < 0.05. Although the majority of the pathways found in the LP/HC group only reached uncorrected *p* < 0.05 and are considered suggestive pathways here, many overlapped with those from the HP/LC group.

In both HP/LC and LP/HC conditions, shared KEGG/GO pathways within the hypothalamus include phosphate
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Birth weights, perinatal weights, and growth rates of pups exposed to different *in utero*/perinatal conditions. (A) Significantly lower birth weights resulted from *in utero*/perinatal exposure to LP/HC maternal diet (*** indicates *p* < 0.0001 by Welch's *t*-test). HP/LC and vehicle were not different from control. Pups were weighed at postnatal day 2. (B) Body weights at weaning (postnatal day 28) for male LP/HC and HP/LC, versus control (* indicates *p* < 0.001 by Welch's *t*-test). (C) Male control, LP/HC and HP/LC, adult growth curves from 1 to 10 months of age. (D) Adiposity at 5 months of age showing increased adiposity in HP/LC males (* indicates *p* < 0.05 by Welch's *t*-test). Sample sizes range from 15 to 27, 15 to 22, and 9 to 20 for HP/LC, LP/HC, and Control groups, respectively.</figcaption></figure>

metabolism, oxidative phosphorylation, neural restructuring, apoptosis, and protein transport. Within the HP/LC condition, KEGG/GO pathways were specifically enriched for synaptic transmission, nucleotide processing, biopolymer modification, MAPK signaling pathway, circadian rhythm, long-term potentiation, and gap junction signaling in the hypothalamus. Within the LP/HC condition, hedgehog signaling was uniquely enriched (**Supplementary Tables S1–S4**). At FDR < 0.05, only *Ebf1* was significant for the HP/LC condition within hypothalamic tissue. *Ebf1* encodes Early B Cell Factor 1, a transcription factor regulating peroxisome proliferator-activated receptor gamma. Importantly, it has been linked to childhood obesity and neuronal differentiation (Moreno and Bronner-Fraser, 2005; Jimenez et al., 2007; Comuzzie et al., 2012).

The liver transcriptome pathway analysis showed differential gene expression for RNA processing, phosphate metabolism, cell cycle, and circadian rhythm in both nutritionally modified states. Purine metabolism, apoptosis, response to DNA damage stimulus and DNA repair within the liver were uniquely enriched within the HP/LC diets. Several pathways enriched among the genes altered by the LP/HC diet within

the liver include biopolymer modification, cellular protein localization/transport, embryonic development, and hedgehog signaling (**Supplementary Tables S1–S4**). A notable proportion of genes within the enriched pathways in liver are related to the biological clock, including *Wee1*, *Per2* and *Slc5a6*, each previously reported to be associated with the development of metabolic syndrome (Englund et al., 2009; Bechtold et al., 2010; He et al., 2016). Both *Wee1* and *Per2* had an FDR of <5%, while *Slc5a6* had an FDR of 0.056, all representing top DEGs in liver. Additional top liver DEGs at FDR < 5% are *Ndrg1* and *Rory* (**Supplementary Tables S5–S8**). *Ndrg1* (*N*-Myc Downstream Regulated 1) has been shown to promote the differentiation of adipocytes from precursor cells through the induction of *Ppary* (peroxisome proliferator-activated receptor) expression (Cai et al., 2017). A single nucleotide polymorphism (SNP) within *Rory* (*RAR*-related orphan receptor gamma) has been implicated in potentially contributing to T2DM susceptibility (Wang et al., 2003). Another study also demonstrated several SNPs within *Rory* to be associated with obesity in cattle (Barendse et al., 2006).

Across the dietary groups and the two tissues, consistent pathways include circadian rhythm, phosphate metabolism, RNA processing, cell cycle/apoptosis, and embryonic development.
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

C

<figure><img src="image_7.png" /></figure>

D

<figure><img src="image_8.png" /></figure>

E

<figure><img src="image_9.png" /></figure>

F

<figure><img src="image_10.png" /></figure>

**FIGURE 2** | Intraperitoneal Glucose Tolerance Test, baseline glucose, baseline insulin, and insulin response at 9 months of age for males LP/HC and HP/LC vs. control. (A) Intraperitoneal Glucose Tolerance Test. (B) Average baseline and 30 min insulin levels Baseline Glucose. (C) Baseline Insulin. (D) Baseline Glucose. (E) Individual IRR. (F) AUC. Median ratio is indicated by the bar in each treatment group. Sample size ranges are 11–12, 12–13, and 11–15 for HP/LC, LP/HC, and Control groups, respectively. (* indicates *p* < 0.05 by Welch's *t*-test).

DISCUSSION

Our study showed that maternal diets with abnormal nutritional compositions such as HP/LC or a LP/HC diet led to developmental differences, which persisted into adulthood affecting glucose metabolism, body weight, adiposity and gene expression levels in the hypothalamus and liver. Male offspring from mothers exposed to LP/HC diet showed a significantly

reduced birth weight, which persisted at weaning (day 28) and up to 10 months of age. Male offspring from the HP/LC group trended toward a higher birth weight, and had a significantly higher body weight at weaning and increased adiposity at 5 months. Additionally, both sets of offspring were found to have an impaired ability to clear glucose in adulthood. The IPGTT conducted in HP/LC and LP/HC litters resulted in a significantly higher AUC for both groups compared with control.
<table><caption>TABLE 1 | KEGG and GO pathway enrichment of differentially expressed genes in the hypothalamus, liver or both under different nutritional exposures in early life.</caption><thead><tr><th rowspan="2">Pathways</th><th colspan="2">Liver</th><th colspan="2">Hypothalamus</th></tr><tr><th>HP/LC</th><th>LP/HC</th><th>HP/LC</th><th>LP/HC</th></tr></thead><tbody><tr><td>Phosphate metabolism</td><td>4.10E-08</td><td>3.40E-04</td><td>1.50E-09*</td><td>5.10E-04</td></tr><tr><td>Biopolymer modification</td><td>-</td><td>5.80E-04</td><td>1.70E-07*</td><td>-</td></tr><tr><td>Cellular localization of protein</td><td>-</td><td>3.00E-04</td><td>-</td><td>7.00E-05</td></tr><tr><td>RNA processing</td><td>3.70E-06</td><td>7.80E-04</td><td>-</td><td>7.10E-05</td></tr><tr><td>Nucleotide processing</td><td>-</td><td>-</td><td>3.80E-05*</td><td>-</td></tr><tr><td>Neural restructuring</td><td>-</td><td>-</td><td>7.90E-05*</td><td>3.00E-03</td></tr><tr><td>Protein transport</td><td>-</td><td>2.00E-04</td><td>8.70E-06*</td><td>5.00E-05</td></tr><tr><td>Cell cycle</td><td>1.00E-05</td><td>2.10E-04</td><td>3.30E-03</td><td>-</td></tr><tr><td>Synaptic transmission</td><td>-</td><td>-</td><td>1.00E-05*</td><td>-</td></tr><tr><td>MAPK signaling pathway</td><td>-</td><td>-</td><td>3.60E-05*</td><td>-</td></tr><tr><td>Gap junction</td><td>-</td><td>-</td><td>7.70E-05*</td><td>-</td></tr><tr><td>Apoptosis</td><td>8.70E-04</td><td>-</td><td>5.00E-04</td><td>1.50E-03</td></tr><tr><td>Purine metabolism</td><td>5.50E-06*</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Hedgehog signaling pathway</td><td>-</td><td>4.50E-03</td><td>-</td><td>5.20E-04</td></tr><tr><td>Colorectal cancer</td><td>1.70E-05*</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Embryonic development</td><td>-</td><td>3.60E-03</td><td>3.90E-03</td><td>1.60E-03</td></tr><tr><td>Oxidative phosphorylation</td><td>-</td><td>-</td><td>3.80E-03</td><td>4.30E-05*</td></tr><tr><td>Circadian rhythm related pathways</td><td>6.30E-05</td><td>2.40E-04</td><td>1.20E-03</td><td>-</td></tr><tr><td>Long-term potentiation</td><td>-</td><td>-</td><td>7.90E-05*</td><td>-</td></tr><tr><td>Response to DNA damage stimulus</td><td>1.80E-06*</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Pathways shown are those that passed Bonferroni-corrected *p* < 0.05 in at least one group or those that are consistent in at least two groups at uncorrected *p* < 0.05. All *p*-values are based on one-sided Fisher's exact test and Bonferroni correction was carried out on KEGG and GO pathways separately. * indicates passing the Bonferroni-corrected *p* < 0.05.

Interestingly, the LP/HC diet presented offspring with reduced baseline insulin and glucose levels, whereas the HP/LC diet had a reduced insulin response 30 min after glucose exposure. This led us to investigate what genes and pathways might underlie these phenotypic differences in the offspring. We found diet- and/or tissue-specific genes and pathways (e.g., gap junction in hypothalamus and purine metabolism in liver for HP/LC, neuronal pathways for both diets in hypothalamus, hedgehog signaling for LP/HC in both tissues), as well as shared pathways across tissues and diets (e.g., cell cycle and circadian rhythm).

There is a well-cataloged history of studies exploring the effects of maternal diet alterations on diseases in adult offspring, specifically within the context of metabolic perturbations (Langley-Evans et al., 1994; Rees et al., 2000; Daenzer et al., 2002; Zambrano et al., 2006). Previous reports have associated both maternal high- and low-protein diets with increased blood pressure, while maternal LP diet results in reduction of body weight in offspring (Jahan-Mihan et al., 2015). Our results indicate that the inadequacy of the low protein maternal diet induces low birth weight and lower body weights throughout life but no difference in adiposity. Additionally, consistent with our previous findings, LP/HC offspring had diminished glucose clearance rates (Bhasin et al., 2009). The HP/LC maternal diet outcomes indicated that too much stimulation of the nutrient-sensing pathways in utero, while not resulting in higher birth weights, did result in higher body weight at weaning, and higher adiposity and glucose intolerance in adulthood. There

are conflicting results regarding the physiological impact of increased maternal protein intake on offspring in previous studies. One study exploring high protein maternal diet on offspring reported no change in body weight, degree of adiposity, or glucose intolerance (Thone-Reineke et al., 2006), while another reported a reduction in body weight with an increase in adiposity (Daenzer et al., 2002). Differences in results may be due to variations in maternal dietary composition, specifically the ratio of protein to complementary nutrients and the total energetic intake; one diet had a protein to carbohydrate ratio of 1.12 while the other was 0.91. Our study utilized a protein-carbohydrate ratio of 0.35, which is significantly lower than the other two studies and may contribute to the variation in results recorded by different groups.

In our study, we essentially found a U-shaped curve in terms of the relationship between maternal protein intake and glucose tolerance in the offspring's adult life. It is predicted that optimal levels for many essential nutrients follow a U-shaped curve: levels that are too low as well as too high result in detrimental health outcomes. Sufficient vitamin A, for instance, is essential in pregnancy to prevent diaphragmatic hernia in the fetus (Clugston et al., 2010), but excess vitamin A results in cleft palate (Ackermans et al., 2011). High protein diets have been advocated as a way to achieve weight loss or prevent weight gain in humans, and our findings suggest that this type of diet, while effective in preventing weight gain, would not be beneficial to the developing offspring. Furthermore, the HP/LC
diet resulted in high baseline glucose levels compared with controls at 10 months of age. Therefore, the current trend of providing high protein to low birth weight infants in Neonatal Intensive Care Units (Hay and Thuren, 2010) to combat postnatal growth restriction, may have far reaching implications to these infants that may prove detrimental toward subsequently developing glucose intolerance, and needs close monitoring upon follow-up.

Our hypothalamic transcriptome and pathway analyses in the offspring indicate alterations in tissue gene programming, likely through the developmental consequences resulting from embryonic nutritional imbalances, that connects maternal diet with long-term physiological outcomes in offspring. Key hypothalamic pathways such as synaptic transmission, neural restructuring and neural development were particularly altered by the HP/LC maternal diet suggesting that these pathways may lead to abnormal neural circuitry in the hypothalamus, which plays a key role in energy balance, ultimately contributing to dysregulation of metabolic homeostasis including increased fat storage and glucose intolerance in adulthood. These findings corroborate a recent study where a high fat maternal diet led to significant alterations in hypothalamic gene expression, resulting in a disruption of regular hypothalamic circuitry through ECM and dopamine availability pathway perturbations (Barrand et al., 2017). While the pathways are not the same between maternal diets, they all reflect changes in neuronal signaling due to maternal dietary alterations. In addition, mitochondrial dysfunction, specifically with regard to oxidative phosphorylation within the hypothalamus has been reported to be associated with neurodegeneration and anorexia in anx/anx mice (Lindfors et al., 2011). This seems to link with our data showing enrichment of oxidative phosphorylation pathways within hypothalamic tissue for both diets and the development of feeding aberrations observed in a previous study (Lindfors et al., 2011). Finally, our finding that *Ebf1* expression is significantly altered by HP/LC maternal diet is suggestive of a metabolic disorder and dysregulation of neural development, as previous studies have linked this gene with childhood obesity and neuronal differentiation (Moreno and Bronner-Fraser, 2005; Comuzzie et al., 2012).

In the liver, circadian rhythm and metabolic pathways such as phosphate metabolism and metabolic regulation were prominent, indicating that these pathways may be major effectors of metabolic syndrome outcomes due to altered maternal diets. Our findings that a number of circadian rhythm related pathways (Table 1), including rhythmic processes and transcriptional regulation via *Per2*, are altered within hepatic tissue under a LP/HC maternal diet is an interesting outcome of the current study and is supported by previous work done in intrauterine growth restricted male rat models (Freije et al., 2015). Additionally, physiological evidence directly supports the role of circadian rhythm dysregulation in increasing the risk for obesity, type 2 diabetes, and other metabolic syndrome traits (Blattler et al., 2012), thus supporting our physiological findings. Notably, around 10% of transcripts oscillate diurnally within the liver (Panda et al., 2002), many of which are involved in energy homeostasis and metabolism, further supporting the

potential pathophysiological repercussions of hepatic circadian dysregulation. Furthermore, the expression of *Slc5a6*, a biotin transporter, is dependent on cell autonomous circadian clock genes *Clock* and *Bmal1* (He et al., 2016). Notably, biotinylation is an essential posttranslational modification of a number of mammalian carboxylases, including acetyl-CoA carboxylase, which has been linked to non-alcoholic fatty liver disease (Tessari et al., 2009).

We also report the alteration of MAPK (for HP/LC) and hedgehog (for LP/HC) signaling pathways, both of which have been previously linked with the facilitation of insulin signaling. Disruption of insulin production or transduction of its signal is intimately associated with generation of metabolic syndrome (Guo, 2013). More specifically, our data suggests that the HP/LC maternal diet disrupts the MAPK activity within the hypothalamus. MAPK facilitates hypothalamic insulin signaling involved in determining the organism's energy and nutrient homeostasis (Mayer and Belsham, 2009; Schultze et al., 2012). The hedgehog-signaling pathway was predominantly affected by the LP/HC condition for both the liver and hypothalamus. The activation of this pathway is highly correlated with the severity of liver damage in patients with non-alcoholic fatty liver disease and is closely linked with hepatic repair mechanisms. Additionally, hedgehog perturbation is involved in both congenital and adult onset metabolic syndrome (King et al., 2008; Choi et al., 2011; Guy et al., 2012).

We acknowledge that this is a discovery study and the findings provide the framework for future avenues of investigation. In particular, pathohistological and structural analysis of both the hypothalamus and liver tissues are necessary to show association between the highlighted pathways/genes and tissue pathology. Additionally, protein level assays and perturbation experiments are necessary to validate the functional and physiological outcomes of the potential causal genes highlighted by our current study. Moreover, it would be interesting to explore further the classification of adiposity and the relevant abundance of white versus brown adipose tissue. This would help characterize more precisely the effects of nutritional perturbation on adipogenesis. A potential limitation of this study is a lack of further validation of our Affymetrix data via RT-PCR, although the validity of microarrays has been highly supported in the literature (Morey et al., 2006; Wang et al., 2006). Additionally, our study mainly focused on identifying the pathways perturbed by diets instead of individual DEGs. The reliability for the significant pathways is high because the random chance to have multiple significant DEGs from the same pathway is extremely low.

Overall, the prevalence of metabolic disorders such as T2DM and metabolic syndrome in modern societies is growing and the etiology is clearly complex, involving diet, physical activity, and possibly endocrine disrupting chemicals in the environment (De Coster and van Larebeke, 2012). Here, we demonstrate that either a deficiency or an excess of protein in early life can result in adverse effects on metabolic outcomes in adulthood, especially with respect to glucose intolerance and T2DM-like phenotypes. Importantly, these early life perturbations altered gene pathways in adulthood, which may underlie the observed metabolic dysregulation, although the cause and effect relationship remains
to be further explored. Our unique study design directly
comparing the effects between low and high protein maternal
intake on both metabolic phenotypes of the offspring and multi-
tissue transcriptomes offers a comprehensive understanding of
the role of maternal protein imbalance in metabolic health. Based
on our findings, a balanced approach to maternal and postnatal
nutrition toward achieving optimal outcomes is necessary for a
normal physiological phenotype in the offspring.

AUTHOR CONTRIBUTIONS

LJM designed and conducted the experiments, researched the
data, analyzed all the data, and drafted the manuscript. QM, SL,
and CP analyzed the gene expression data. MB analyzed and
summarized all the data, and wrote the manuscript. SX and WT
conducted the experiments. JW researched and analyzed the data.
SD developed experimental design and assisted in phenotypic

experiments and interpretation. AJL developed the experimental
design. XY supervised the data analysis and manuscript writing.
All authors reviewed and edited the manuscript.

FUNDING

AJL was supported by the National Institutes of Health (NIH) grants HL28481 and HL30568. LJM was supported by NIH HD07228. SD was supported by NIH HD-81206 and HD-41230. XY was supported by NIH DK104363.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: https://www.frontiersin.org/articles/10.3389/fgene.
2018.00642/full#supplementary-material

REFERENCES

Ackermans, M. M., Zhou, H., Carels, C. E., Wagener, F. A., and Von den Hoff, J. W. (2011). Vitamin A and clefting: putative biological mechanisms. *Nutr. Rev.* 69, 613–624. doi: 10.1111/j.1753-4887.2011.00425.x

Barendse, W., Bunch, R. J., Kijas, J. W., and Thomas, M. B. (2006). The effect of genetic variation of the retinoic acid receptor-related orphan receptor C (RORC) gene on fatness in cattle. *Genetics* **175**, 843–853. doi: 10.1534/genetics.106.064535

Barrand, S., Crowley, T. M., Wood-Bradley, R. J., De Jong, K. A., and Armitage, J. A. (2017). Impact of maternal high fat diet on hypothalamic transcriptome in neonatal Sprague Dawley rats. *PLoS One* 12:e0189492. doi: 10.1371/journal.pone.0189492

Bechtold, D. A., Gibbs, J. E., and Loudon, A. S. (2010). Circadian dysfunction in disease. *Trends Pharmacol. Sci.* **31**, 191–198. doi: 10.1016/j.tips.2010.01.002

Bhasin, K. K., van Nas, A., Martin, L. J., Davis, R. C., Devaskar, S. U., and Lusis, A. J. (2009). Maternal low-protein diet or hypercholesterolemia reduces circulating essential amino acids and leads to intrauterine growth restriction. *Diabetes* **58**, 559–566. doi: 10.2337/db07-1530

Blattler, S. M., Cunningham, J. T., Verdeguer, F., Chim, H., Haas, W., Liu, H., et al. (2012). Yin Yang 1 deficiency in skeletal muscle protects against rapamycin-induced diabetic-like symptoms through activation of insulin/IGF signaling. *Cell Metab.* **15**, 505–517. doi: 10.1016/j.cmet.2012.03.008

Cai, K., El-Merahbi, R., Loeffler, M., Mayer, A. E., and Sumara, G. (2017). Ndrg1 promotes adipocyte differentiation and sustains their function. *Sci. Rep.* **7**:7191. doi: 10.1038/s41598-017-07497-x

Choi, S. S., Omenetti, A., Syn, W. K., and Diehl, A. M. (2011). The role of hedgehog signaling in fibrogenic liver repair. *Int. J. Biochem. Cell Biol.* **43**, 238–244. doi: 10.1016/j.biocel.2010.10.015

Clugston, R. D., Zhang, W., Alvarez, S., de Lera, A. R., and Greer, J. J. (2010). Understanding abnormal retinoid signaling as a causative mechanism in congenital diaphragmatic hernia. *Am. J. Respir. Cell Mol. Biol.* **42**, 276–285. doi: 10.1165/rcmb.2009-0076OC

Comuzzi, A. G., Cole, S. A., Laston, S. L., Voruganti, V. S., Haack, K., Gibbs, R. A., et al. (2012). Novel genetic loci identified for the pathophysiology of childhood obesity in the hispanic population. *PLoS One* 7:e51954. doi: 10.1371/journal.pone.0051954

Daenzer, M., Ortmann, S., Klaus, S., and Metges, C. C. (2002). Prenatal high protein exposure decreases energy expenditure and increases adiposity in young rats. *J. Nutr.* **132**, 142–144. doi: 10.1093/jjn/132.2.142

De Coster, S., and van Larebeke, N. (2012). Endocrine-disrupting chemicals: associated disorders and mechanisms of action. *J. Environ. Public Health* 2012:713696. doi: 10.1155/2012/713696

Englund, A., Kovanen, L., Saarikoski, S. T., Haukka, J., Reunanen, A., Aromaa, A., et al. (2009). NPAS2 and PER2 are linked to risk factors of the metabolic syndrome. *J. Circadian Rhythms* **7**:5. doi: 10.1186/1740-3391-7-5

Freije, W. A., Thamotharan, S., Lee, R., Shin, B. C., and Devaskar, S. U. (2015). The hepatic transcriptome of young suckling and aging intrauterine growth restricted male rats. *J. Cell. Biochem.* **116**, 566–579. doi: 10.1002/jcb.25008

Frias, A. E., and Grove, K. L. (2012). Obesity: a transgenerational problem linked to nutrition during pregnancy. *Semin. Reprod. Med.* **30**, 472–478. doi: 10.1055/s-0032-1328875

Ganguly, A., and Devaskar, S. U. (2008). Glucose transporter isoform-3-null heterozygous mutation causes sexually dimorphic adiposity with insulin resistance. *Am. J. Physiol. Endocrinol. Metab.* **294**, E1144–E1151. doi: 10.1152/ajpendo.90251.2008

Gluckman, P. D., Hanson, M. A., Cooper, C., and Thornburg, K. L. (2008). Effect of in utero and early-life conditions on adult health and disease. *N. Engl. J. Med.* **359**, 61–73. doi: 10.1056/NEJMra0708473

Guo, S. (2013). Insulin signaling, resistance, and the metabolic syndrome: insights from mouse models to disease mechanisms. *J. Endocrinol.* **220**, T1–T23. doi: 10.1530/JOE-13-0327

Guy, C. D., Suzuki, A., Zdanowicz, M., Abdelmalek, M. F., Burchette, J., Unalp, A., et al. (2012). Hedgehog pathway activation parallels histologic severity of injury and fibrosis in human nonalcoholic fatty liver disease. *Hepatology* **55**, 1711–1721. doi: 10.1002/hep.25559

Hay, W. W., and Thureen, P. (2010). Protein for preterm infants: how much is needed? How much is enough? How much is too much? *Pediatr. Neonatol.* **51**, 198–207. doi: 10.1016/S1875-9572(10)60039-3

He, L., Hamm, J. A., Reddy, A., Sams, D., Peliciari-Garcia, R. A., McGinnis, G. R., et al. (2016). Biotinylation: a novel posttranslational modification linking cell autonomous circadian clocks with metabolism. *Am. J. Physiol. Heart Circ. Physiol.* **310**, H1520–H1532. doi: 10.1152/ajpheart.00959.2015

Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nat. Protoc.* **4**, 44–57. doi: 10.1038/nprot.2008.211

Jahan-Mihan, A., Rodriguez, J., Christie, C., Sadeghi, M., and Zerbe, T. (2015). The role of maternal dietary proteins in development of metabolic syndrome in offspring. *Nutrients* **7**, 9185–9217. doi: 10.3390/nu7115460

Jimenez, M. A., Åkerblad, P., Sigvardsson, M., and Rosen, E. D. (2007). Critical role for Ebf1 and Ebf2 in the adipogenic transcriptional cascade. *Mol. Cell. Biol.* **27**, 743–757. doi: 10.1128/MCB.01557-06

King, P. J., Guasti, L., and Laufer, E. (2008). Hedgehog signalling in endocrine development and disease. *J. Endocrinol.* **198**, 439–450. doi: 10.1677/JOE-08-0161
Langley-Evans, S. C., Phillips, G. J., and Jackson, A. A. (1994). In utero exposure to maternal low protein diets induces hypertension in weanling rats, independently of maternal blood pressure changes. *Clin. Nutr. Edinb. Scotl.* 13, 319-324. doi: 10.1016/0261-5614(94)90056-6

Lindfors, C., Nilsson, I. A., Garcia-Roves, P. M., Zuberi, A. R., Karimi, M., Donahue, L. R., et al. (2011). Hypothalamic mitochondrial dysfunction associated with anorexia in the anx/anx mouse. *Proc. Natl. Acad. Sci. U.S.A.* 108, 18108-18113. doi: 10.1073/pnas.1114863108

Mayer, C. M., and Belsham, D. D. (2009). Insulin directly regulates NPY and AgRP gene expression via the MAPK MEK/ERK signal transduction pathway in mHypoE-46 hypothalamic neurons. *Mol. Cell. Endocrinol.* 307, 99-108. doi: 10.1016/j.mce.2009.02.031

Moreno, T. A., and Bronner-Fraser, M. (2005). Noelins modulate the timing of neuronal differentiation during development. *Dev. Biol.* 288, 434-447. doi: 10.1016/j.ydbio.2005.09.050

Morey, J. S., Ryan, J. C., and Dolah, F. M. (2006). Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. *Biol. Proced. Online* 8, 175-193. doi: 10.1251/bpo126

Oken, E., and Gillman, M. W. (2003). Fetal origins of obesity. *Obes. Res.* 11, 496-506. doi: 10.1038/oby.2003.69

Panda, S., Antoch, M. P., Miller, B. H., Su, A. I., Schook, A. B., Straume, M., et al. (2002). Coordinated transcription of key pathways in the mouse by the circadian clock. *Cell* 109, 307-320. doi: 10.1016/S0092-8674(02)00722-5

Rees, W. D., Hay, S. M., Brown, D. S., Antipatis, C., and Palmer, R. M. (2000). Maternal protein deficiency causes hypermethylation of DNA in the livers of rat fetuses. *J. Nutr.* 130, 1821-1826. doi: 10.1093/jn/130.7.1821

Saklayen, M. G. (2018). The global epidemic of the metabolic syndrome. *Curr. Hypertens. Rep.* 20:12. doi: 10.1007/s11906-018-0812-z

Schultze, S. M., Hemmings, B. A., Niessen, M., and Tschopp, O. (2012). PI3K/AKT, MAPK and AMPK signalling: protein kinases in glucose homeostasis. *Expert Rev. Mol. Med.* 14:e1. doi: 10.1017/S1462399411002109 doi: 10.1159/000243170

Storey, J. D. (2003). The positive false discovery rate: a Bayesian interpretation and the q-value. *Ann. Statist.* 31, 2013-2035. doi: 10.1214/aos/1074290335

Tessari, P., Coracina, A., Cosma, A., and Tiengo, A. (2009). Hepatic lipid metabolism and non-alcoholic fatty liver disease. *Nutr. Metab. Cardiovasc. Dis.* 19, 291-302. doi: 10.1016/j.numecd.2008.12.015

Thone-Reineke, C., Kalk, P., Dorn, M., Klaus, S., Simon, K., Pfab, T., et al. (2006). High-protein nutrition during pregnancy and lactation programs blood pressure, food efficiency, and body weight of the offspring in a sex-dependent manner. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* 291, R1025-R1030. doi: 10.1152/ajpregu.00898.2005

Wang, H., Chu, W., Das, S. K., Zheng, Z., Hasstedt, S. J., and Elbein, S. C. (2003). Molecular screening and association studies of retinoid-related orphan receptor γ (RORC): a positional and functional candidate for type 2 diabetes. *Mol. Genet. Metab.* 79, 176-182. doi: 10.1016/S1096-7192(03)00096-9

Wang, Y., Barbacioru, C., Hyland, F., Xiao, W., Hunkapiller, K. L., Blake, J., et al. (2006). Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. *BMC Genomics* 7:59. doi: 10.1186/1471-2164-7-59

Zambrano, E., Bautista, C. J., Deas, M., Martinez-Samayoa, P. M., Gonzalez-Zamorano, M., Ledesma, H., et al. (2006). A low maternal protein diet during pregnancy and lactation has sex- and window of exposure-specific effects on offspring growth and food intake, glucose metabolism and serum leptin in the rat. *J. Physiol.* 571(Pt 1), 221-230. doi: 10.1113/jphysiol.2005.100313

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Martin, Meng, Blencowe, Lagarrigue, Xiao, Pan, Wier, Temple, Devaskar, Lusis and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.