Lipocalin-2 (Lcn2) expression is mediated by maternal nutrition during the development of the fetal liver

Abstract

The mechanisms by which maternal protein deficiency programs insulin action in the offspring are poorly understood. The interpretation of transcriptomics is complicated by homeostatic adaptations, for example, changes in amino acid metabolism, which are potentially unrelated to the programming mechanism. The fatty acid composition of the maternal diet modulates the programming of insulin action, offering a possible strategy to circumvent these complications. Fetal livers harvested on d21 of gestation from pregnant rats fed high-protein (18 % w/w) and low-protein (9 % w/w) diets prepared with either corn or soya oil were screened with rat genome microarrays. Although a low-protein maternal diet altered the abundance of more than one hundred mRNAs in the fetal liver, only 40 were changed by the fatty acid composition of the diet (P < 0.05). One of these mRNAs was identified as lipocalin-2 (Lcn2). This pattern of differential expression was confirmed by qRT-PCR. The expression of Lcn2 was decreased by low-protein diets when the diet contained soya oil, whereas the effect of protein was much smaller in the group fed diets prepared with corn oil. The decrease in Lcn2 expression produced by soya oil persisted into adult life. Levels of the Lcn2 protein were closely correlated to the mRNA abundance. The results suggest a possible involvement of Lcn2 in the programming of hepatic function.

Introduction

Despite a body of epidemiological evidence showing that a restricted maternal diet during gestation enhances the risks of metabolic disease in the offspring, molecular markers associated with the underlying mechanisms remain elusive (Gluckman et al. 2009; Tarry-Adkins and Ozanne 2011). The pregnant rat fed a low-protein diet is one of several well-established animal models in which poor fetal growth leads to adverse effects on blood pressure, glucose metabolism, and insulin action in adult life (Langley-Evans et al. 1996; Fernandez-Twinn et al. 2005; Rees et al. 2006; Warner and Ozanne 2010). Restricting the protein content of the maternal diet causes a wide range of physiological changes in both dam and fetus; however, many of these are homeostatic adaptations, with no long-lasting effects on the offspring, for example maintaining the amino acid supply in response to reduced protein intake. Nutrigenomic techniques (proteomics and transcriptomics) cannot differentiate between adaptive homeostatic responses and underlying programming mechanisms, limiting the value of these otherwise powerful tools. Fortunately, the nutritional basis of the low-protein model is more complicated. The final phenotype of the offspring is dependent not only on the protein content but also on other macronutrients in the maternal diet (Langley-Evans 2000). The fatty acid composition of the diet is a particularly important factor; diets prepared with soya oil have been shown to enhance the programming effect of protein deficiency compared to those prepared with corn oil. For example, the amount of insulin released in response to a glucose challenge is reduced in the offspring of dams fed protein-restricted diets prepared with soya oil (low in n-3 fatty acids), whereas when the maternal diet is prepared with corn oil (high in n-3 fatty acids), there is a minimal effect of protein deficiency, i.e., corn oil-based diets suppress the programming of postnatal glucose metabolism (Maloney et al. 2007). This interaction suggests that key pathways and molecules associated with the programming of insulin action should only be sensitive to the protein content of the maternal diet in the fetuses of dams fed the diets prepared with soya oil, whereas homeostatic changes associated with protein and amino acid metabolism should be present in both soya and corn oil groups.

The fetal liver is a key target for developmental programming with a central role in the metabolism of carbohydrates and fats (Desai et al. 1997; Maloney et al. 2003), and as with the glucose tolerance test, the programming of insulin action in the liver also depends on the fatty acid composition of the maternal diet. The levels of mRNAs coding for acyl CoA carboxylase-1 (Acc-1, a regulator of fatty acid synthesis) and liver-type carnitine acyl transferase-1 (L-CPT-1, a regulator of fatty acid beta-oxidation and transport) are higher following insulin stimulation in the offspring of dams fed the low-protein diet compared to those fed the high-protein diet prepared with soya oil; however, when the same diet is prepared with corn oil, the effect of protein restriction is much smaller (Maloney et al. 2007). To identify possible marker genes for developmental events in the liver, we have investigated differentially expressed mRNAs using microarray technology (transcriptomics). We reasoned that by comparing the patterns of gene expression in the fetal liver of rats fed diets with two levels of protein (18 vs. 9 %) prepared with corn (C) or soya oil (S) (C vs. S) in a 2 × 2 factorial design, we would be able to eliminate those mRNAs that responded to protein or oil alone. The remaining mRNAs that were either elevated or suppressed in the livers of the fetuses of dams fed low-protein diets prepared with soya oil compared to corn oil are likely to be of interest as both markers and possible components of the programming mechanism.

Methods

All experimental procedures were approved by the ethical review committee of the Rowett Research Institute and conducted in accordance with the UK Animals (Scientific Procedures) Act, 1986. Four groups of twelve female Rowett Hooded-Lister rats were randomly assigned to one of four experimental diets from approximately 8 weeks of age. The composition of the diets, containing 18 % w/w (18) or 9 % w/w casein (9) and 7 % (w/w) of either C or S has been described previously (Maloney et al. 2007). After feeding the diets for approximately 2 weeks, the animals (body weight 200–250 g) were mated with normal males. The animals continued to be fed the experimental diets for the duration of gestation. On gestation day 21, dams were anesthetized, killed by cardiac puncture, and rapidly dissected. After the fetuses were killed, tissues were rapidly removed, frozen in liquid nitrogen, and stored at −80 °C until required.

Microarray analysis of gene expression

Six dams were chosen from each treatment group, excluding those with large (>12) or small (<8) litters. Female fetuses were used as these were shown to be more affected in adult life (Maloney et al. 2007). The fetal liver from one female fetus was chosen randomly from each litter, and a portion was used for the preparation of total RNA with Trizol reagent (Sigma Poole Dorset UK) followed by DNase treatment as described previously (Maloney et al. 2005). All mRNA samples used had an RNA integrity number >8 (Agilent Bioanalyser), and sexing was verified by PCR. Samples of RNA were hybridized to Affymetrix Rat Genome 230 2.0 microarrays (Affymetrix, Santa Clara CA. USA, provided through the NuGO array pipeline) by ServiceXS (Leiden NL) using their standard protocol.

qRT-PCR

Samples of 200 ng total RNA were reverse transcribed using the TaqMan Reverse Transcription Reagents Kit (Applied Biosystems, Warrington, Cheshire, UK) primed with random hexamers. Primer and probe sets for lipocalin-2 (Lcn2) were supplied by Applied Biosystems (Cat No: Rn00590612_ml) or as described previously (McNeil et al. 2008). The levels of cDNA were measured relative to the 18S ribosomal RNA, and relative target quantity was calculated from the standard curve.

Western blotting

Samples of fetal liver were homogenized in 10 volumes of buffer (150 mM NaCl, 1 % NP-40, 0.5 % sodium deoxycholate, 0.1 % SDS, 10.8 mM sodium phosphate pH 7.4, 0.1 mM PMSF, 5 mg/ml aprotinin). After centrifugation at 10,000×g for 15 min at 4 °C, 100 μg of protein from the supernatant was separated on 15 % SDS-PAGE. The gels were blotted to Immobilon-P membranes (Millipore, Watford, Herts., UK) Antibodies were obtained from R&D Systems (Abingdon, Oxfordshire, UK), and conditions for Western blotting were as recommended by the manufacturer. The primary antibody (lipocalin-2/NGAL, Cat No: AF3508) was reconstituted as per the manufacture’s instructions and used at 1:200 dilution (1 μg/ml). The secondary antibody (Anti-goat IgG-HRP, Cat No: HAF019) was used at 1:1,000 as per the manufacturer’s instructions. The blots were visualized by ECL (Immobilon Western Chemiluminescent HRP substrate Cat No WBKLSO500, Millipore, Watford UK) and imaged with a Fuji Imager LAS3000. Blots were subsequently stripped and reprobed with anti-beta-actin (A5441, Sigma Poole Dorset UK, used at a dilution of 1:10,000 and anti-mouse IgG-HRP, Sigma Cat No A9044, used at a dilution of 1:10,000).

Statistical analysis

The data generated by the microarrays were analyzed using the MadMax suite of programs (https://madmax.bioinformatics.nl). All of the data met the appropriate quality control criteria and were normalized using the GCRMA (slow) program. Relative changes in mRNA levels were determined using linear models for microarray data, (LIMMA) module of GenePattern (Reich et al. 2006). The corrected gene expression data were analyzed by two-way ANOVA using Genstat (12th/13th edition, VSN International Ltd, Hemel Hempstead, UK). qRT-PCR and Western blot data were analyzed by ANOVA or REML as appropriate.

Results

Maternal live weight gain was reduced in animals fed the protein-restricted diet (P < 0.05), and the effect was greater in corn compared to soya oil diets (interaction protein.oil P < 0.05) as shown in Table 1. The protein and oil in the diet did not affect litter size. Fetal growth was reduced by protein restriction (P < 0.05), but it was unaffected by the oil type. Protein-restricted diets reduced the weight of the fetal liver by approximately 8 %, making the reduction in fetal growth comparable to the previous studies (Rees et al. 1999).

Table 1 Maternal and fetal characteristics at day 21 of gestation

Following normalization, a total of 11,426 mRNAs were found to be expressed above the background level in the fetal liver. In the first instance, we compared the patterns of mRNA abundance in the 18S and 9S groups, as this diet composition produced the greatest difference in hepatic gene expression in adult life (Maloney et al. 2007) and found a total of 110 mRNAs that were different (P < 0.025). These candidate mRNAs were then analyzed by two-way ANOVA using the full dataset of soya and corn oil groups, testing for the effects of protein (Pr), oil (O), and their interactions (Pr·O). This analysis identified 40 mRNAs, the abundance of 6 mRNAs altered by both protein and oil (Pr and O, P < 0.05), and 34 where there was evidence for an interaction between the protein and oil content of the maternal diet (Pr·O, P < 0.05). The remaining mRNAs were either not different or only influenced by the protein content (n = 9 Pr, P < 0.05; O, P > 0.05) or oil in the maternal diet (n = 3 Pr, P > 0.05; O, P < 0.05). The ratio of abundance of the mRNAs that differed in the 18S:9S groups ranged from 4.2:1 to 0.7:1 with 25 of these mRNAs downregulated by the 9 % protein diet. Of these mRNAs, only 12 were increased by more than 1.25-fold or reduced by 0.75-fold by protein restriction as shown in Supplementary Table 1. For comparison, a similar analysis starting with a comparison of the 18C and 9C groups is shown in Supplementary Table 2. This analysis identified 17 mRNAs increased by more than 1.25-fold or reduced by 0.75-fold by protein restriction.

The mRNA identified as Lcn2 showed the largest change (4.2-fold) in response to protein restriction in the fetuses of the dams fed the diet containing soya oil, with little change in the fetuses of dams fed the corn oil diet (Pr, P = 0.045; O, P = 0.005; Pr·O = n.s.), a pattern of expression that closely matched the postnatal phenotype. Since the Lcn2 mRNA was expressed at a high level, it was chosen for further study. Its expression was assessed by qRT-PCR using the mRNA samples used in the initial analysis. Expression was significantly different (P < 0.05) and showed a correlation between the abundance as determined by qRT-PCR and microarray analysis (linear regression P < 0.05, R ² > 0.25). As an additional confirmation, the expression of Lcn2 was also examined in mRNA extracted from the liver of sibling fetuses and the combined data are shown in Fig. 1a. The Lcn2 mRNA abundance is lower in the group fed 9 % protein (P = 0.007) and also lower in the group where the maternal diet contains soya oil (P = 0.008). The effects were additive with no significant interaction between them (P = 0.27). The abundance of the Lcn2 mRNA in the livers of the male fetuses (Fig. 1b) showed a similar numerical pattern with respect to the maternal diet but was more variable than in the females, with the result that there were no significant differences. Lcn2 expression was also analyzed in samples of liver from 6-week-old female offspring from a previous experiment (Fig. 1c). The expression relative to the 18S RNA was substantially lower than at the fetal stages, but the mRNA for Lcn2 was still detectable. Expression in the offspring of dams fed the maternal diet containing soya oil was less than in the offspring of dams fed the corn oil diet (P = 0.035); however, there was no residual effect of protein (P = 0.16).

Fig. 1

Expression of lipocalin-2 (Lcn2) determined by qRT-PCR. Expression is relative to the 18S ribosomal RNA in arbitrary units. Data are mean ± SEM. Columns sharing a superscript letter are not significantly different (P > 0.05). a Expression in fetal liver of female fetuses on d21 of gestation 18S n = 8, 9S n = 10, 18C n = 8, 9C n = 7. Data were analyzed by REML Protein P = 0.007 Oil P = 0.008 Protein·Oil P = 0.273. b Expression in fetal liver of male fetuses on d21 of gestation (n = 6 for all groups). c Expression in liver of female offspring at 6 weeks of age (n = 6 for all groups). Data were analyzed by two-way ANOVA Protein P = 0.167, Oil P = 0.035, Protein·Oil P = 0.94

To confirm that the mRNA differentially expressed was full-length Lcn2, samples of cDNA were amplified using primers covering exons 1–5 of the full-length Lcn2 mRNA. This produced a product of the anticipated size and sequence (not shown). Western blotting with an antibody against Lcn2 detected a single band of 25 kDa (Fig. 2). Lcn2 expression is a part of the innate response to infection (Zhang et al. 2008), and a band of similar size was present in a sample of liver from an adult rat (Sal) infected with Salmonella enterica, which was included as a positive control. Analysis of the band density in fetal liver samples by two-way ANOVA showed expression was reduced in the soya oil group (P < 0.05) and tended (P = 0.069) to be reduced in the group fed the 9 % protein diets. There was a linear correlation (P < 0.001, R ² = 0.32) between the density of the spots on the Western blot and gene expression estimated by qRT-PCR. Immunocytochemistry of sections of fetal liver showed staining in discrete cells; however, this was not associated with specific anatomical features (data not shown).

Fig. 2

Lipocalin-2 protein levels in protein extracted from the liver of female fetuses at d21 of gestation. a The pixel density of replicate 100 μg samples ± SEM (n = 6). Data were analyzed by one-way ANOVA P = 0.017, and columns sharing a superscript letter are not significantly different (P > 0.05). A typical blot is shown in the lower panel Sal = 25 μg protein extracted from the liver of an adult rat infected with Salmonella enterica (Robertson et al. 2003). b Correlation between Lcn2 protein expression (pixel density on Western blots) and mRNA (relative expression by RT-PCR) of a subset of samples, n = 6 per treatment. Line fitted by linear regression (P < 0.001, R ² = 0.32)

Discussion

Numerous studies have made use of transcriptomics to evaluate complex phenotypes programed by nutritional insults and in many cases have identified genes associated with cell proliferation (Mortensen et al. 2010; Swali et al. 2011; Oster et al. 2012; Vijaya et al. 2013) and signaling pathways such as IGF-1 and progesterone (McNeil et al. 2009), a pattern consistent with the known function of these pathways in response to a restricted protein intake. However, it is clear that the programing by protein restriction is not simply a product of limited growth, which in the case of the fetal liver is equally affected by protein-restricted diets prepared with either corn or soya oils. Despite a restriction in hepatic growth, the adult offspring of dams fed a corn oil-based diet exhibit very little programming of insulin action, differences are only apparent in the offspring of dams fed diets prepared with soya oil (Maloney et al. 2007). This study has employed this observation to differentiate between these non-specific growth effects and identify potential gene markers associated with the programming mechanism(s). Maternal protein deficiency alters the abundance of slightly less than one hundred and fifty of the more than eleven thousand mRNAs present in the developing fetal liver irrespective of the maternal diet containing corn or soya oils. Interestingly, these mRNAs show markedly different patterns with very few common hits in soya versus corn. Unfortunately, the small number of targets limited the power of gene ontology analysis.

The mRNA with an abundance that most closely matched the postnatal phenotype was identified as Lcn2, also known as neutrophil gelatinase-associated lipocalin, siderocalin, and 24p3. The expression of the mRNA correlated closely with the abundance of the corresponding protein and was reduced when the dam was fed a low-protein diet, but crucially decreased further when the diet contained n-3 fatty acids from soya oil as opposed to n-6 from corn oil. This pattern of a substantial decrease in the expression of Lcn2 in 9S compared to 18S, and a much smaller change when comparing 9C to 18C, closely resembles the pattern of programming of both glucose tolerance and hepatic gene expression in the offspring at 6 months of age (Maloney et al. 2007). Furthermore, the change in hepatic Lcn2 expression associated with exposure to the soya oil diet in utero persists into adult life. This close correlation strongly suggests that Lcn2 may be part of the mechanism associated with the programming of hepatic glucose metabolism and gene expression.

Lipocalin-2 is a 24 kDa secreted glycoprotein expressed in a variety of tissues. It has been shown to have a role in the sequestration of iron and is an important acute phase protein (Zhang et al. 2008). Lipocalin proteins bind a range of hydrophobic molecules including fatty acids, retinoids, and steroids (Akerstrom et al. 2000; Logdberg and Wester 2000; Wang 2012; Chakraborty et al. 2012). Although the endogenous ligand of Lcn2 has yet to be identified, there is evidence that it may be involved in the metabolism of fatty acids (Law et al. 2010) and may play a role in remodeling the plasma membrane by redistributing polyunsaturated fatty acids among different types of phospholipids (Yang et al. 2012). Fatty acid analysis (McNeil and Rees unpublished) showed that n-3 fatty acids were lower (with DHA reduced by approximately 65 %) in animals fed the corn oil compared to the soya oil diets, and that there was also a corresponding small but significant (P < 0.05) increase in arachidonic acid (AA). Protein deficiency had no major effect on the abundance of DHA or AA. These differences in the ratios of n-3 and n-6 polyunsaturated fatty acids could be an important driver of changes in Lcn2. The promoter region of the Lcn2 gene contains multiple transcription factor binding sites, which include glucocorticoid (Garay-Rojas et al. 1996) and estrogen response elements (Seth et al. 2002), suggesting that endocrine factors may also play a role in regulating Lcn2 expression. A number of studies have shown that exposure to excessive concentrations of glucocorticoids during fetal development programs hyperglycemia in the offspring (Drake et al. 2007); however, there was no correlation between corticosterone concentrations in fetal plasma and Lcn2 expression (McNeil and Rees unpublished data). Gender is also an important factor as the maternal diet programs hepatic lipid metabolism in female offspring, whereas in males, parameters like blood pressure are more sensitive. As Lcn2 has been shown to modify estradiol production and estrogen receptor expression in metabolic tissues (Guo et al. 2012), gender-specific effects of a persistent change in its expression may become apparent later in life.

The Lcn2 mRNA is present in the liver, kidney, and spleen during organogenesis in the mouse (Garay-Rojas et al. 1996); however, animals lacking the Lcn2 gene develop normally with no apparent defects (Berger et al. 2006). In rodents, elevated serum Lcn2 concentrations are a characteristic of several models of obesity (Yan et al. 2007), and mice with a targeted deletion of the Lcn2 gene have disturbances in glucose and lipid metabolism when they are fed high-fat diets (Guo et al. 2010). In vitro, Lcn2 regulates apoptosis, growth, and development in a variety of cell types (Devireddy et al. 2005), suggesting that it has a secondary role in development. Such small non-critical changes in the structure and organization of organs and tissues are ideally suited to the regulation of metabolic plasticity. In this study, we have concentrated on the latter stages of hepatic development in the rat as this has been shown to be a critical period for the programming of insulin action in the low-protein model (Langley-Evans et al. 1996; Martin-Gronert and Ozanne 2010). Just before birth, there is a rapid increase in the weight of the fetal liver, which appears to be very sensitive to the maternal diet (Rees et al. 1999; Gruppuso et al. 2005). This suggests a possible involvement of Lcn2 in the changes in hepatic function, which take place around the time of birth and are associated with the acquisition of a more mature hepatocyte function. Lcn2 has been shown to be a regulator of apoptosis in the hematopoietic compartment (Liu et al. 2011b), and the surge in Lcn2 expression may be associated with the transfer of hematopoietic cells from the liver to the bone marrow late in gestation.

In adult humans, elevated plasma concentrations of the homologous protein LCN2 have been reported in obese and insulin-resistant individuals, although these findings have not been replicated in all populations (Catalan et al. 2009; Wang et al. 2007; Wallenius et al. 2011; Liu et al. 2011a). It is interesting that LCN2 concentrations were increased in the plasma of pregnant women, with a high body mass index and furthermore were correlated with markers of insulin resistance (Cesur et al. 2012). Serum LCN2 was also higher in women with gestational diabetes and preeclampsia (D’Anna et al. 2008, 2009). These findings suggest that LCN2 may also be of value as a possible marker for fetal programming in humans.