Transcriptional Profiling of Rats Subjected to Gestational Undernourishment: Implications for the Developmental Variations in Metabolic Traits

A link has been established between prenatal nutrition and the development of metabolic and cardiovascular diseases later in life, a process referred to as developmental programming. It has been suggested that the trajectory of development is shifted by alterations in the maternal nutritional state leading to changes in developmental plasticity, in part underpinned by epigenetic changes in gene regulation. However, to date, only candidate gene approaches have been used to assess expression and molecular changes in the offspring of maternally undernourished animals. Furthermore, most work has focused on animals at an age where the programmed phenotype is already manifest and little is known about changes in gene expression in the offspring prior to development of obesity and related metabolic disorders. Gene expression profiles of liver, retroperitoneal white adipose fat, and biceps femoris skeletal muscle tissue from young adult male rats (55 days old) in which nutritional status had been manipulated in utero by maternal undernutrition (UN) were compared to the profiles of offspring of ad libitum fed mothers serving as the control group (AD) (8 offspring/group). The expression profiles were determined using the Illumina RatRef-12 BeadChip. No significant changes in expression were identified for skeletal muscle or white adipose tissue. However, studies of liver tissue showed 249 differentially expressed genes (143 up regulated, 106 down regulated). Although the animals at day 55 have yet to develop obesity they already show biochemical abnormalities and by day 110 express a phenotype characterized by increased adiposity and altered insulin sensitivity. An analysis of pathways affected suggests that intrauterine programming of UN animals to favor fat as an energy source results in mitochondrial dysfunction which initially affects the postnatal hepatic function and subsequently, via the resultant metabolic changes in other organs leads to the evolution of a phenotype similar to that of the metabolic syndrome.


Introduction
Developmental plasticity is the process by which environmental factors acting during critical windows in early life affect the developmental trajectory and thus the mature phenotype. Such changes are often irreversible and are generally considered to be underpinned at least in part by epigenetic processes [1,2]. While adaptive in origin developmental plasticity can have maladaptive consequences [3] and thus in humans can affect health [4].
The developing fetus is completely dependent on its mother and so it follows that maternal nutrition during pregnancy, in particular, has a strong influence on the intrauterine environment. The DOHaD (Developmental Origins of Health and Disease) hypothesis proposes that the metabolic syndrome (i.e. coronary heart disease, type 2 diabetes, stroke, hypertension) originates through developmental plasticity in response to undernutrition during fetal life and infancy [1,5,6]. It has been hypothesized that these diseases of later life arise from a mismatch [7] between the actual postnatal environment and that predicted by the fetus during the phase of developmental plasticity [8,9]. A variety of maternal nutritional manipulations in the rat have induced offspring who develop obesity, insulin resistance, and endothelial dysfunction in later life including global maternal undernutrition as used in this study [10] and maternal low protein by isocaloric diets [11].
The mechanistic basis of such changes is only now being explored. It has been proposed that epigenetic modifications of the fetal genome establish new patterns of gene expression that reset the metabolic state that is maintained into adulthood [12]. Effects on candidate genes using such rat models have focused on genes known to be significant in the metabolic syndrome including peroxisomal proliferator-activated receptor (PPARa) and the glucocorticoid receptor (GR) [11,13,14]. These observations suggest that there are methylation differences in the promoters of these genes and changes in mRNA levels between the offspring of normally and undernourished mothers. A broader picture of the cellular functions and genetic pathways that may be implicated can be obtained using a non-biased, global microarray approach to study gene expression in key tissues between offspring of control and undernourished mothers.
In the present study, an established rat model of balanced maternal undernutrition [10] has been used for investigation of gene expression differences in target tissues (liver, retroperitoneal white adipose fat, and biceps femoris skeletal muscle) between offspring of control and undernourished mothers. We studied young adult rats at day 55 of age when they have not yet developed the abnormal phenotype in order to identify any gene expression changes that may play a role in predisposing these animals to the development of the metabolic syndrome. The findings show that no significant gene expression changes are detected between the two groups in skeletal muscle and white adipose tissue. However, the liver reveals 249 significant gene expression differences (p,0.05) showing significant down regulation of expression of several key genes affecting the entry of glucose into the cell and its metabolism via the glycolytic and tricarboxylic acid pathways. This would result in a substantial impairment in their ability to utilize carbohydrate on restoration of a normal diet. In addition, there was a marked increase in expression of genes associated with the intracellular trafficking of fatty acids, their peroxisomal degradation and their entry into the mitochondrion to undergo b-oxidation. This might suggest that following intrauterine undernutrition these animals have lost the flexibility to switch between carbohydrate and fat as energy sources. The analysis also reveals a substantial impairment of the electron transport chain as evidenced by a down regulation in expression of 3 components of complex I and the ATP synthase that generates ATP from the proton gradient. There is also a deficit in expression of a key enzyme in the biosynthesis of the NAD proton carrier.
We suggest that maternal undernutrition leads to offspring that favor fat as an energy source thus resulting in mitochondrial dysfunction that initially affects postnatal hepatic function and subsequently, other organs.

Animals
This study exploited a well defined model of developmental programming via maternal undernutrition [10,15]. Briefly, virgin Wistar rats (age, 10065 days) were time mated using a rat estrous cycle monitor to assess the stage of estrous of the animals before introducing the male. After confirmation of mating, rats were housed individually in standard rat cages with free access to water. All rats were kept in the same room with constant temperature maintained at 25uC and a 12-h light/12-h darkness cycle. Animals were assigned to one of two nutritional groups: a) undernutrition (30% of ad libitum) of a standard chow diet throughout gestation (UN group), b) standard chow diet ad libitum throughout gestation (AD group). Food intake and maternal weights were recorded daily until the end of pregnancy. After birth, pups were weighed and litter size was adjusted to 8 pups per litter to assure adequate and standardized nutrition until weaning. Pups from undernourished mothers were cross-fostered onto dams that had received AD feeding throughout pregnancy. At weaning (day 22), animals were weight matched within maternal dietary group to 2 animals per cage. Eight animals per group were culled at day 55 following an overnight fast to represent an age known to precede the development of increased adiposity and altered insulin sensitivity in UN offspring. Tissues (liver, retroperitoneal white adipose fat, and biceps femoris skeletal muscle) were immediately dissected and snap-frozen in liquid nitrogen for molecular analysis. A subcohort of animals (n = 8 per group) were maintained until postnatal day 110 to represent an age where phenotypic changes in growth and metabolism between AD and UN offspring have previously been reported. [10,16] The Animal Ethics Committee of The University of Auckland approved all animal manipulations.

Phenotypic Measurements
Serial dual-energy x-ray absorptiometry (DEXA) analysis was performed on the animals at day 55 and day 110 to establish body composition and assessment of bone parameters. In addition to DEXA, body adiposity was assessed via standard techniques (fat depot dissection and weighing) following post-mortem. Body weight and food and water intake were monitored. Plasma analysis incorporated IGF-1, insulin, C-peptide, and lipid profiles. Statistical analyses were performed with the Knowledge Discovery Environment (InforSense Ltd., London, UK) statistical package. Differences between groups were determined by a T-test for difference in variance with an FDR correction, and data are shown as an adjusted p-value6standard error.

Microarray Hybridization and Analysis
Global gene expression analysis was performed using the Illumina Bead Array Platform. Briefly, total RNA was isolated from cells with TRIzol reagent (Invitrogen, Paisley, Scotland, UK) followed by column purification (RNeasy, Qiagen, UK). RNA quality was assessed with the Lab-on-a-chip system (Agilent, Palo Alto, CA). The RNA samples were amplified following the Illumina TotalPrep RNA amplification protocol. The fluorescent labeling and hybridization of RNA samples were also performed according to custom protocols defined in the manufacturer's instructions. The samples from liver, skeletal muscle, and white adipose fat tissue for male offspring (n = 8) were hybridized to the RatRef-12 BeadChip available from Illumina (sampling 21,290 genes from the rat transcriptome). The hybridization signals were captured and quantified using a Bead Array scanner and the associated Bead Studio software. Bead Studio includes a spectrum of internal controls for determining data quality. These check for consistency in signal intensity, expected background and noise levels, mismatches versus perfect matches, and signal in housekeeping genes (Figure S1A-S1H). Data was exported and normalization was performed in R (Bioconductor) using the Limma [17] and Lumi [18] Bioconductor packages for microarray analysis. Data was filtered for those showing expression on at least one array and a variance-stabilizing transformation [19] was done to take advantage of the multiple technical replicates (on average 20-30 per array) for each gene available on each Illumina chip. Quantile normalization was chosen to remove technical variability between arrays. Density plots of intensity and box plots of amplitude were used to visualize the raw data compared to the filtered and normalized data ( Figure S2A-S2D). Limma was used for pairwise comparisons between AD and UN using a linear model to compute p-values, which were adjusted using the Benjamini and Hochberg multiple testing correction. The derivation of the p-value reflects the degree of variance between biological replicates and thus is a measure of the confidence in assigning significance of small fold changes in gene expression. Genes with a p-value,0.05 were considered to be significantly differentially expressed. Gene ontology analysis was done using DAVID: Database for Analysis, Visualization, Integrated Discovery [20] and Ingenuity Pathway Analysis software (Ingenuity, Stanford USA).

Quantitative Real-Time PCR
The comparative expression levels in AD and UN RNA of seven genes identified as differentially expressed from our microarray data were verified by qRT-PCR, together with the expression levels of five genes not present on the array. Cyclophilin was used for the generation of a standard curve and normalization of the RNA concentration. Each qRT-PCR was performed in triplicate. Real-time RT-PCR was performed in 96-well white plates (Abgene) using the Verso SYBR Green 2-Step qRT-PCR Fluorescein Kit (Thermo Scientific, UK) according to the manufacturer's protocols; the resulting fluorescence was quantified using an iCycler system (Bio-Rad). The threshold cycle value (C T ) was obtained for each well as the cycle number at which the measured fluorescence crossed the arbitrary threshold of 150 units (all values are in log phase). The average C T was calculated for each gene in each sample. Data were normalized to cyclophilin, with DC T calculated as follows: DC T = C T (test)2C T (cyclophilin). DDC T values were then calculated as the change in DC T for the UN sample relative to the DC T value for the AD sample using the equation: DDC T = 2 CT(test) 22 CT(cyclophilin) . This data is plotted in Figure 1 and primer sequences are shown in Table S1.

Phenotypic Assessment
For each of the 8 animals in each of the AD and UN groups, a series of phenotypic measurements relevant to the metabolic syndrome was recorded at day 55 (Table 1), the age at which gene expression analysis was performed, and in a subcohort of animals at day 110 ( Table 2). From these measurements, it is clear that at day 55 there are very few differences in key indicators (insulin resistance, glucose intolerance, elevated blood lipids and increased adiposity) as might be expected at this early stage in the progression of the metabolic syndrome. However, by day 110 there are clear indications of abnormalities in insulin sensitivity and markers of adiposity including total fat mass and plasma leptin concentrations.

Expression Analysis
Expression profiling was done for liver, white adipose fat, and skeletal muscle tissue from young adult (55-day-old) male rats born to either AD or UN dams. In each group there were eight biological replicates, all of which were raised postnatally on an ad libitum diet. Expression profiles were determined using the Illumina RatRef-12 BeadChip that permits interrogation of 21,290 rat genes. Genes with a p-value,0.05 were considered to show significant differential expression between the two groups. Using these criteria, no significant gene expression changes between the two groups were identified for skeletal muscle or white adipose tissue in the day 55 animals. However, the comparison for liver revealed a list of 249 differentially expressed genes. Table S2 summarizes the expression values for all genes represented on the array and Table S3 the list of 249 differentially expressed genes with associated p-values and fold changes between the two groups.  Quantitative RT-PCR (qRT-PCR) Analysis of Selected Genes. In order to confirm the expression data derived from array analysis, seven genes ( , and RBP1 (XM_230637)) found to be differentially expressed were selected for qRT-PCR on the pooled RNA samples for the biological replicates in each treatment group. Cyclophilin was used as an invariant control for normalization. Figure 1A summarizes the data and illustrates that the differences in DDCt values between the AD and UN RNA populations for each gene are consistent with the expression changes determined by array analysis. These values were then converted to fold change to directly compare the qRT-PCR values to the expression values from the microarray ( Table 3). Previous analysis has provided evidence that the expression of the PPARa, GR and AOX (acyl-CoA oxidase) genes is altered in the liver of offspring from normally nourished mothers compared to mothers who had either a low protein diet or balanced undernutrition during pregnancy [11,21]. It has also been suggested that PEPCK (phosphoenolpyruvate carboxykinase) and 11b-HSD2 (11b-hydroxysteroid dehydrogenase type 2) may also be genes influenced by diet [22,23]. As these genes were not represented by any of the oligonucleotides on the RatRef-12 array, they were assessed by qRT-PCR in the AD and UN RNA populations using cyclophilin as the normalization control. The findings are shown in figure 1B for the 8 biological replicates in each treatment group. For PPARa and PEPCK, there is no evidence for any differences in expression levels between the two treatment groups. However, GR, AOX and 11b-HSD2 show between 1.5 to 2-fold increases in expression in the UN treatment group.

Discussion
Animals at postnatal day 55 were chosen for gene expression analysis because they had completed puberty but had not begun to exhibit the features of the metabolic syndrome phenotype as observed later in adulthood and, as shown in the present study, as early as day 110 postnatally [10,24]. Our results presented show that while there are differences in gene expression in hepatic tissue in day 55 males, there are no observed differences in white adipose fat or skeletal muscle. This is an intriguing finding, as both fat and skeletal muscle is adversely affected in the adult. This suggests that there is a cascade of biological processes that amplifies over time to lead to the phenotypic abnormalities. Previously Ozanne et al. [25] has shown, in a different but related experimental model of maternal low protein intake, that the expression of a limited number of genes in skeletal muscle is minimal early in life and accumulates later in life. It would appear that major changes in gene expression in the liver precede those in other metabolically relevant tissues.
The data for all three tissues is derived from 8 biological replicates in each treatment group with technical replication on each array varying from 20-30 fold. This means that for each gene represented on the array, there are between 160 and 240 measurements of expression level. This permits an accurate estimate of the variance that is used to derive the final p-value and provides confidence in the detection of small gene expression changes. Thus the changes found in liver are based on many measurements and are highly significant. This is further underscored by the absence of any significant changes in gene expression between the two groups in skeletal muscle or white adipose fat tissue subjected to the same analysis. The qRT-PCR analysis of candidate genes (PPARa and PEPCK) showed very little difference in expression levels between the two groups in the day 55 male offspring. This contrasts with the 10-fold, 3-fold and 3-fold increase for PPARa, GR and AOX, respectively, in the day 34 offspring born to nutritionally restricted mothers as reported by Lillycrop et al. [11]. Four factors may explain this difference. First, the restricted diet is different (50% protein restriction only) whereas the diet in the study reported here is 30% of normal chow for all nutrients and may impact differently on gene expression. Second, in the Lillycrop et al. study the animals were weaned on day 28 and spent only 6 days in the post-weaning period compared to 33 days in this present study. Third, the composition of the postweaning diet in the present study is grain-based as opposed to a purified diet containing 50% sucrose in the Lillycrop et al. study. Finally, the male and female tissues were pooled in the Lillycrop et al. study and there may be gender differences we have not explored. Similar to previous reports, GR, 11b-HSD2 and AOX show increases of between 1.5 to 2-fold in the offspring of undernourished mothers used in this study [11].  Analysis of the gene categories and associated pathways associated with differentially expressed genes (Table S3) has revealed some relevant changes that potentially impact on the metabolic phenotype. Table 4 summarizes the relevant categories to emerge from analysis of gene functions (some of which emerged from DAVID and Ingenuity analyses) associated with the list of differentially expressed genes. Figure 2A and 2B summarize the interrelationships between pathways involved in glucose, fat, and energy metabolism and highlight differentially regulated genes that influence these pathways.

Glucose Metabolism
In utero the fetus relies mainly on maternal glucose for its energy supply and only in late gestation develops the capacity for glucose storage and mobilization and for gluconeogenesis in preparation for postnatal life. Glucose and glycogen metabolism appear to have developed inappropriately in the livers of UN animals [26,27]. Many of the changes in hepatic gene expression we detected are relatively modest and this may reflect alterations in hepatic development. Burns et al. reported changes in the balance of periportal to periarterial hepatocytes which have different metabolic profiles in the offspring of low protein fed dams [28].
The livers of UN animals show changes in the anomalous expression of RNA for several enzymes involved in glycogenesis, glycogenolysis and glycolysis. There is a 54% increase in the expression of the glycogen synthase gene (GYS2), however, this is counterbalanced by a marked decrease in expression of glucan branching enzyme (GBE1) that is essential for the formation and storage of glycogen. Furthermore, there is a modest decrease in phosphoglucomutase (PGM1) that converts the glucose-1-phosphate resulting from glycogenolysis to glucose-6-phosphate destined for glycolysis. There is also a deficit of triose phosphate isomerase that could result in an inefficient conversion of fructose 1-6 bisphosphate to glyceraldehyde-3 phosphate in the glycolytic process. The reduced expression of triose phosphate isomerase would therefore affect the efficiency of glucose conversion to pyruvate and might indicate an increased accumulation of dihydroxyacetone phosphate. This can readily be reduced to glycerol-3-P that is the backbone for phosphatidates and signaling phospholipids. Overproduction of these could disrupt intracellular function. The efficiency of the citric acid cycle (TCA) would be compromised by the reduction in dehydrogenase E1 (DHTKD1). This enzyme catalyses the conversion of a-ketoglutarate to succinyl CoA-a key step in the TCA. Its subsequent conversion to succinate is the only step in glucose catabolism that produces high-energy phosphate (GTP) directly, without the need for oxidative phosphorylation. The GTP can be used for nucleotidespecific metabolic reactions or it can convert ADP to ATP. Production of oxaloacetate further on in the cycle would also be affected. This would slow the rate of formation of citric acid and the entry into the TCA of acetyl-CoA from oxidative decarboxylation of pyruvate. Supplementation of oxaloacetate from metabolism of branched-chain amino acids would also be reduced. This means that most of the oxaloacetate is required as an intermediate for TCA function and less is available for gluconeogenesis. It is interesting to note that there is a reduction in the transcripts of other enzymes involved in the catabolism of amino acids that supply intermediates to the TCA cycle. In the fetus, this impasse may be circumvented to an extent by the presence of the glyoxalate cycle in which acetate can be converted to either glucose or succinate. However, this pathway disappears after birth and hence the inability for gluconeogenesis would be exacerbated postnatally.

Insulin and IGF-1 Signalling
The insulin and IGF-1 receptors are closely related structurally and share many post-receptor signaling mechanisms. The livers of UN animals show changes in expression of genes whose products are known to attenuate the PI3 kinase signaling pathway that is common to both receptors. Both IGFBP2 and Rbp1 are known to interact with components of the PI3Kinase/ Akt signaling pathway IGFBP2 is likely to act by sequestering IGF-1 and preventing binding to its receptor and hence inhibiting downstream activation of PI3K [29]. Rbp1 has been shown to interfere with the heterodimerization of the p85/p110 subunits of Akt [30]. This appears to be a retinoic acid independent effect of Rbp1. An increase in expression of Tribbles homolog 2 (TRB2) has been shown to interfere with the phosphorylation step by PI3K that activates Akt, the effector kinase of this signaling cascade. Conversely increased expression of the protein tyrosine phosphatase (PTPRB) will increase the rate of dephosphorylation of Akt to its inactive form. The concerted effect of these gene expression changes is likely to have the effect of blunting insulin signaling and increasing insulin insensitivity. It is interesting to note that while there is not a significant increase in the levels of circulating insulin, there is a discernable upward trend in UN animals that could indicate the initial stages in the establishment of insulin resistance. One of the major targets of Akt is the insulin-dependent glucose transporter GLUT4 that on activation migrates from the ER to the plasma membrane and facilitates glucose entry into the cell. An additional perturbation of this process in UN animals is suggested by the increased expression of p45/STRAP that has been shown to interfere with GLUT4 [26,31] trafficking between the ER and the cell membrane. Taken together these findings suggest that the livers of UN animals have altered insulin sensitivity and glucose homeostasis.
A primary function of IGF-1 is to drive postnatal cell proliferation and growth, principally through the MAP kinase signaling pathway. It is only present in low levels in the fetus (where IGF-2 predominates) and increases dramatically in the perinatal period. The liver is the major source of circulating IGF-1. The actions of IGFs are modulated by the presence of a number of binding proteins (IGFBPs), one of which IGFBP2, shows a .2 fold increase in expression in liver of UN animals. IGFBP2 is secreted into the blood and is thought to have a negative effect on IGF-1 activity. In this study, the highly significant reductions in body length and body, spleen, and liver weights lend support to this role. The levels of IGFBP2 have been shown to be high in the early to mid-term fetus and then drop markedly towards term. The persistence of elevated IGFBP2 in UN animals is therefore an abnormality that could have a wide range of effects on postnatal growth and organ development.

Mitochondrial Activity
The array analysis indicates a number of changes in the expression of genes that could affect the efficiency of electron transport and ATP generation in the mitochondrion in liver from UN animals. Transcripts of NADH dehydrogenase subunits (NDUFC2, NDUFB6 and NDUFA8) of complex I of the electron transport chain are all reduced as is quinolate phosphoribosyltransferase (QPRT) required for the synthesis of NAD, the electron carrier. In addition, the ATP synthase (ATP5o) that generates ATP from the electron gradient is reduced by a similar amount. Changes are also observed in peptidylprolyl isomerase D (PPID) that controls the mitochondrial permeability transition pore and TIM14, a component of the mitochondrial protein import motor. Taken together these findings suggest a degree of     Table 4. Cont.
mitochondrial dysfunction in the livers of UN animals that could significantly affect the energy balance of the tissue and may be the precursor of later adiposity. Since the liver acts as the initial nutrient sensor, this defect could affect not only its own metabolism but also that of other organs that respond to its systemic signals.

Fat
In this well-characterized model, it has been shown that undernourishment of dams during pregnancy predisposes offspring to later adiposity and hyperinsulinaemia when fed a normal diet [21] but only at an age later than we have undertaken gene expression analysis. In the present study, even as early as day 110, male offspring of UN dams are showing features of the metabolic syndrome related to increased adiposity and hyperinsulinemia and hyperleptinemia. The absence of any differences in fat expression profiles between control and UN animals at day 55 may suggest that the development of adiposity is secondary to early metabolic perturbations in postnatal life. It is interesting to note that at day 35 UN rats show the onset of a sedentary phenotype (and hence reduced energy expenditure) before the presence of obesity [32]. This may reflect early changes in neurological gene expression modifying behaviour.
The elevated expression of RNA for fatty acid binding protein (FABP7) suggests that in livers from UN animals there is an increase in the intracellular trafficking of fatty acids. This is paralleled by increases in expression of AOX (from qRT-PCR analysis), which catalyses the first step in the peroxisomal degradation of fatty acids, and carnitine acyltransferase (CRAT) that shuttles fatty acids into the mitochondrial matrix. It is interesting to note that the expression of the first enzyme in the biosynthetic pathway of carnitine, trimethyllysine hydroxylase (Tmlhe) is down regulated and this may reflect a homeostatic response to excess fatty acid trafficking. These findings suggest that more energy is derived from b-oxidation of fatty acids in these animals. However, this can only be achieved if the NADH from boxidation regenerates NAD and forms ATP. The observed defects (especially in complex I) indicate a disruption in this process. The provenance of these fatty acids is not clear since there were no significant differences between treatment groups in the circulating levels of either triglycerides or free fatty acids although LDL levels are elevated by maternal undernutrition. Other possible sources of fatty acids could be from increased endogenous synthesis or from turnover of cellular components. There is no evidence from the microarray analysis for an increase in enzymes for the latter; indeed the modest decrease in levels of diacylglycerol acyltransferase (DFAT2) in UN animals would indicate a blunting of triglyceride synthesis in liver.
The increased expression of SREBP cleavage activating protein (SCAP) suggests that there may be higher levels of sterol synthesis taking place in livers of undernourished animals. The formation of a transcriptionally active sterol response element binding protein (SREPB) involves the proteolytic action of SCAP [33]. In support of these observations the levels of circulating LDL cholesterol in the UN group is significantly elevated. Interestingly a 60% increase in 11b-HSD2 transcripts was seen in the livers of UN animals accompanied by a similar increase in the expression levels of GR. 11b-HSD2 appears in late gestation/postnatal tissues in human fetal tissues. [34] Elevated levels of this enzyme in the livers of UN animals would be a mechanism to moderate the effects of high levels of circulating corticosterone that may arise in response to nutritional deprivation/stress.

Ribosomal Proteins and Protein Turnover
An intriguing finding of the analysis is an increase (in the range 7-24%) in transcripts in UN animals of 13 of 80 recognized ribosomal proteins (79 of which are represented on the array by multiple oligonucleotides) that occur in the large and small ribosomal subunits. Eight of the proteins are components of the mitochondrial ribosome. Similar increases in expression are also seen in brix (BXDC1), involved in ribosome biogenesis, NPG1, a nucleolar GTPase required for the maturation and nuclear export of pre-ribosomes, and Cebpz, a gene important in the processing of rDNA involved in the 60 s subunit [35]. These findings suggest a greater degree of metabolic flux in ribosomes from UN animals, however it is difficult to understand the significance of these observations given the fact that on the basis of the current annotation only about 16% of the ribosomal proteins show changes in expression levels. It is possible that the increased turnover of ribosomal proteins is contributing to the higher levels of urea observed in the day 55 UN animals.

Effects on Gene Transcription
It is notable that a group of genes that influence transcription show upregulation to varying extents. Some (such as the Gtf and MED factors) have a general effect on the transcriptional apparatus. The most notable increase is in the Rpb1 gene known to be expressed highly in the liver, particularly in hepatic stellate cells [36,37]. In liver Rpb1 has an important role in the production of retinoid derivatives that activate the RAR and RXR retinoic acid receptors [38]. RXR can heterodimerise with the PPAR proteins to form active nuclear receptors that target genes involved in fat metabolism and energy homeostasis; for example Crat and AOX are shown to be upregulated in the UN offspring. Thus elevated or abnormal Rbp1 expression may have a major role in establishing the metabolic syndrome phenotype.
Two other genes show significant changes in gene expression that may impact hepatic function. The first is Col14a1 (a member of the collagen gene family) that was originally called undulin [39,40]. This gene is expressed in hepatic stellate cells that have a role in the formation of the hepatic extra cellular matrix [41] and shows a 20% increase in gene expression. This may be a marker of some underlying pathology as Col14a1 has been associated with the rearrangement of connective tissue occurring in hepatic fibrosis [39]. It is of note that Rbp1 expression changes have also been associated with hepatic fibrosis [37] and thus the gene may be involved in the regulation of Col14a1. The second gene is Adra1b (adrenergic receptor alpha 1b showing a 79% upregulation) that has been associated with carbohydrate metabolism in the liver [42] and disruption of glucose homeostasis when the gene is inactivated [43,44]. Given that the role of the receptor is to stimulate glycogenolysis and the process of glycogen formation is impaired in UN animals, the significance of increased expression is not clear. This may reflect an irreversible change in response to the intrauterine nutritional state or a postnatal compensatory mechanism to the altered metabolic profile set in utero.

Conclusions
Studies to date have primarily focused on a candidate gene approach in animals where the classical features of the metabolic syndrome in the programmed phenotype are already evident. There is a relative paucity of data on gene expression, either by candidate or an array-based approach, on key tissues related to the metabolic syndrome in animals at an age preceding development of the metabolic phenotype. The gene expression changes in appropriate pathways in the livers of the male offspring of maternally undernourished dams at day 55 would suggest that these animals may be predisposed to a persistent and perturbed ability to coordinate fat and carbohydrate metabolism with a shift to a use of fatty acids as an energy source. We have also shown that these animals proceed to develop a phenotype similar to that of the metabolic syndrome as early as postnatal day 110. Array analysis of tissues at the day 110 time point was beyond the scope of the present trial. However, although the present study cannot directly correlate the observed gene expression changes at day 55 with the phenotype at day 110, it provides clear evidence of disturbed hepatic function in a number of key genes related to lipid oxidation and mitochondrial function at a pre-phenotypic age. There are parallels with the observations of Koves et al. [27] and Sparks et al. [45] where diet induced obesity in adult life induces oxidative stress associated with increased b-oxidation of fat metabolism, impaired switching of carbohydrate substrates, depletion of TCA cycle intermediates, and decreased expression of components of the electron transport chain. They propose that mitochondrial dysfunction may be driving the pre-diabetic insulin resistant state. The observations in our rat model indicate that the liver is irreversibly programmed to respond to a nutritionally restricted environment and that this persists into early adulthood (day 55). At this age, similar changes are not observed in adipose tissue or skeletal muscle, suggesting that the liver can meet the immediate energy requirements of these peripheral tissues and that it manifests metabolic abnormalities in advance of the full metabolic syndrome phenotype. These data have uncovered potential candidate genes and pathways that when perturbed lead to the development of the metabolic syndrome in older animals and therefore provide a focus for more detailed gene-specific studies.