Breast-feeding is associated with a moderate, but consistent, reduction in the risk of childhood obesity and metabolic disease in later life( Reference Dewey 1 – Reference Arenz, Ruckerl and Koletzko 5 ). Specific aspects of breast-feeding such as feeding behaviour and its unique nutritional composition may contribute to this long-term protective effect. Human milk (HM) contains several biologically active compounds such as oligosaccharides, hormones and growth factors that are not present in infant milk formula (IMF), supporting optimal growth and development( Reference Burrin, Stoll and Moore 6 – Reference Carlson 9 ). Moreover, HM lipids have a distinct physical structure composed of large lipid droplets (average mode diameter of 4 µm) consisting of a TAG core enveloped by a tri-layer of membranes, the native milk fat globule membrane (MFGM)( Reference Gallier, Gragson and Cabral 10 – Reference Rueda 12 ). MFGM components are secreted by the epithelial cells of the mammary gland and consist of phospholipids and sphingolipids, membrane-specific proteins, enzymes, cholesterol, glycoproteins and vitamins. In contrast, the diameter of lipid droplets in standard IMF is about 0·3–1·0 μm and is not surrounded by a membrane, due to the industrial processing procedures to achieve stable and reproducible end products. The IMF fat globules mainly consist of TAG derived from plant oils and have milk proteins adhering to the surface, with or without small amounts phospholipids (PL)( Reference Michalski 13 – Reference Fave, Coste and Armand 16 ).
Several studies have indicated that dietary exposure to a MFGM lipid-containing product may have beneficial effects. Clinical intervention studies with MFGM products in adults have shown long-term anti-carcinogenic effects, hypocholesterolaemic effects, suppression of multiple sclerosis and improved learning and cognitive performance after consuming the MFGM extract (reviewed by Spitsberg( Reference Spitsberg 17 )). Dietary exposure to MFGM lipid product in early life may have beneficial effects on growth and (cognitive) development in rodents and infants( Reference Vickers, Guan and Gustavsson 18 , Reference Timby, Domellöf and Hernell 19 ). Timby et al.( Reference Timby, Domellöf and Hernell 19 , Reference Timby, Lönnerdal and Hernell 20 ) showed that infants fed a low-energy, low-protein, MFGM-enriched formula showed cognitive development and serum cholesterol concentrations more similar to those of breast-fed infants at 12 months of age as compared with a standard IMF. The composition of MFGM has been shown to affect pancreatic gastric lipase activity and lipase interaction, cholesterol availability, lipid absorption and absorption of bile salts( Reference Michalski 13 , Reference Michalski, Briard and Michel 14 , Reference Timby, Domellöf and Hernell 19 ). However, no growth or associated metabolic effects of MFGM have so far been reported in clinical studies.
We hypothesised that in particular the physical structure of dietary lipids in early life contributes to a healthy development of body composition at longer term via effects on growth and metabolic development during infancy. It has been clearly shown that dietary lipid structure affects lipolysis and the metabolic fate of fatty acids (FA) due to distinct differences in lipid digestion and absorption kinetics( Reference Armand, Hamosh and Mehta 21 , Reference Armand, Pasquier and Andre 22 ). Any difference in postprandial lipid handling can impact lipid availability for the development of metabolic organs. This in turn could programme metabolic homoeostasis, energy balance and metabolic response with potential impact on later life health.
We have previously shown that exposing young mice to a concept diet (Nuturis®; Nutricia Research) containing PL derived from the MFGM fraction, coated around large lipid droplets, prevented fat accumulation and improved metabolic profile in adulthood( Reference Oosting, Kegler and Wopereis 23 , Reference Oosting, van Vlies and Kegler 24 ). These data indicated that exposure to an altered dietary lipid structure in early life may be a key determinant of later life metabolic health. Thus far, however, it is not known how specific aspects of the Nuturis® concept, that is, lipid droplet size and/or coating, may contribute to these later life beneficial effects. The objective of the present study was to determine whether large droplet size and/or coating are responsible for the previously found long-term protective effects against adult adiposity. To this end, we used a validated mouse model for nutritional programming( Reference Oosting, Kegler and Boehm 25 ) in which young pre-weaning mice were exposed to intervention diets comprising different features of the Nuturis® concept. After the early diet intervention, all mice were switched to a Western-style diet (WSD) and body composition development was monitored.
Methods
Animals and study design
All experimental procedures were approved by the Animal Experimental Committee (DEC Consult) and complied with the principles of good laboratory animal care. The animals were housed at facilities of the Wageningen University and Research Centre in a 12 h light–12 h dark cycle (light on 06.00 hours=Zeitgeber time 0 h). Room temperature and humidity were controlled (21 (sem 2)°C and 50 (sem 5) %, respectively). During the entire protocol, food and water were available ad libitum. Between postnatal (PN) day 42 and 98, food intake was measured per cage and body weight (BW) was determined per litter before weaning and individually after weaning twice a week. Male and female C57BL/6J mice were time mated. Dams were assigned to the American Institute of Nutrition-93 (AIN93G) diet during pregnancy and lactation. After birth, litters were culled to four male and two female pups and randomly cross-randomised. Litters were randomly allocated to one of the experimental diets (see below) from PN16 until PN42. At weaning on PN21, male pups were housed in pairs and continued on their respective diets until PN42. From PN42 until dissection at PN98, animals were exposed to a WSD (20 %, w/w fat) (Fig. 1). To control for the effects of the WSD challenge in adolescence and adulthood, a reference (REF) group was included, which was fed the control IMF diet from PN15 to PN42 and the AIN93M diet from PN42 to PN98.
Experimental diets
All experimental diets were semi-synthetic (Research Diet Services; Table 1) and had macronutrient and micronutrient compositions according to the AIN formulation of AIN93G-purified diets for laboratory rodents. The experimental diets contained 28·3 %, w/w IMF powder complemented with protein and carbohydrates to meet rodent nutrient requirements. The fat fraction was completely derived from the IMF. The detailed composition and processing procedure of the diets have been described previously( Reference Oosting, van Vlies and Kegler 24 ). The diets were provided as a dough in order to preserve the lipid structure of the products.
S, small droplets; L, large droplets; Scoating, small droplets coated by milk fat globule membrane (MFGM); Lcoating, large droplets coated by MFGM.
* Peak diameter of particle size distribution.
† Sum volume of particles of 2–12 mm (% total volume of all particles).
The experimental diets differed according to lipid droplet size (small (S) v. large (L)) and presence or absence of MFGM as lipid droplet coating: S (control IMF diet), L, Scoating and Lcoating (Nuturis® IMF diet). The MFGM originated from butter serum (SM 2; Corman Food Industry). In the production of anhydrous milk fat from cream, the aqueous streams are a by-product, rich in milk PL, present as MFGM fragments and referred to as butter serum( Reference Vanderghem, Bodson and Danthine 26 ).
The oil blend of IMF was a mixture of oils such as rapeseed, sunflower, fish, coconut and palm oil. After IMF production, protein, lactose, lipid, PL concentrations and droplet size were measured (Table 1). The protein content was analysed according to the Dumas principle; the used N-factor was 6·25 (ISO14891 IDF 185, 2002). The lactose content was determined by an enzymatic method based on Nederlands Normalisatie-Instituut (NEN) 6853, 1989. Lipid content was determined according to the Van Gulik method (ISO11870, 2000). PL content was analysed with P-NMR by Spectral Services. The particle size distribution was determined by laser light scattering. The measurements were performed at 20°C, using a refractive index for fat of 1·46 and 0·01 absorption. The refractive index for the dispersant (water) of 1·33 was used (Mastersizer 2000, Hydro 2000G; Malvern Instruments Limited).
The WSD consisted of 20 %, w/w fat (3 %, w/w soya oil; 17 %, w/w lard and 0·1 %, w/w cholesterol).
Body composition
Body composition was repeatedly measured between PN42 and 98 by dual-energy X-ray absorptiometry scan under general anaesthesia (isoflurane–N2O–O2) using a PIXImus imager (GE Lunar). Relative fat mass (FM) was calculated as the percentage FM (%FM) of total BW. Between PN42 and PN98, BW, lean body mass (LBM), FM and %FM gain (δBW, δLBM, δFM and δ%FM, respectively) were calculated by subtracting PN42 values from PN98.
Blood sampling and dissection
On PN98, mice were anaesthetised (isoflurane–N2O–O2) after fasting overnight and killed by bleeding (eye extraction) and cervical dislocation. Blood samples were collected in K3EDTA-coated 1-ml micro tubes. Plasma was obtained by centrifugation at 1350 g for 12 min at 4°C and stored at −80°C. The liver, pancreas, musculus tibialis and brain as well as epididymal (EPI), retroperitoneal (RP), perirenal (peri) and inguinal (ING) white adipose tissue (WAT) depots were dissected, weighed, snap-frozen in liquid N2 and stored at −80°C.
Plasma analyses
Plasma glucose concentration was measured colorimetrically (glucose oxidase-phenol and 4 aminophenazone (GOD-PAP) method; Roche Diagnostics) and analysed using a micro plate imaging system (Bio-Rad Laboratories, Inc.). Insulin concentration was measured using ELISA (DRG) according to the manufacturer’s protocol. Homoeostasis model assessment of insulin resistance as calculated from fasting plasma glucose and insulin concentrations (glucose (mmol/l)×insulin (pmol/l)/22·5) as an indirect measure of insulin sensitivity. Total cholesterol, HDL-cholesterol, LDL-cholesterol, VLDL-cholesterol and TAG concentrations were determined colorimetrically after enzymatic conversion using a Roche Hitachi 717 analyzer (Reinier de Graaf Laboratory).
Gene expression in retroperitoneal white adipose tissue
Total RNA was isolated from the RP WAT using Trizol–chloroform (Invitrogen) and purified using a RNeasy mini kit (Qiagen Benelux B.V.). Contaminating genomic DNA was removed with the RNase-free DNase set (Qiagen Benelux B.V.). Quality and quantity of RNA were determined by NanoDrop 2000 (Thermo Scientific) and Bioanalyzer (Agilent). Complementary DNA (cDNA) was synthesised using the iScript cDNA synthesis kit (Bio-Rad) according to the manufacturer’s instructions; 9·4 ng of RP fat depot was used as input for each quantitative PCR (qPCR) reaction. 5× Hot FIREPol Evagreen qPCR mix Plus (Bio-Connect) was used according to the manufacturer’s instructions and qPCR was accomplished with a 7900HT Fast Real-Time PCR System (Applied Biosystems). Leptin (marker adiposity), mesoderm-specific transcript/paternal expressed gene 1 (mest/peg1; marker adipocyte size) and preadipocyte factor 1 (Pref1; marker preadipocytes) mRNA expressions were measured. Normalisation of qPCR data was achieved using qBase PLUS (Biogazelle). qBase is a free programme for the analysis of qPCR data, and for this experiment several housekeeping genes (calnexin, ribosomal protein gene L19, ribosomal protein S29) were used to normalise the qPCR data. For detailed description, see the published article of Hellemans et al.( Reference Hellemans, Mortier and De Paepe 27 ). The primer sequences are depicted in Table 2.
Mest/Peg1, mesoderm-specific transcript/paternal expressed gene 1; Pref1, preadipocyte factor 1; Canx, calnexin; Rpl19, ribosomal protein gene L19; Rps29, ribosomal protein S29.
Statistical analyses
Statistical analyses were performed using SPSS 19.0 (SPSS Benelux). Effects of diets on body composition, organ weights, plasma and gene expression parameters were analysed using mixed effects regression models, including droplet size and coating and size by coating interaction as fixed effects. In all regression models, the correlation among animals within the same cage was accounted for by their sharing a common random effect. Pairwise comparisons were adjusted for multiple comparisons, using the least square difference (LSD) approach. Data are shown as mean values with their standard errors. Differences were considered significant when P<0·05 and a trend was defined at P<0·1.
Results
Programming of adult food intake and body composition by dietary lipid structure
Average energy intake from PN42 to PN98 was comparable between groups. To control for the effects of the WSD challenge, a REF group was included. At day 98, BW, LBM, FM and %FM of the REF group were lower than all other groups fed WSD (Fig. 2). However, FM of REF mice was significantly lower (approximately 64 %) compared with S, L and Scoating groups fed the WSD, but not significantly different from FM found in the Lcoating group that showed comparable FM gain with the REF group.
At day 98, mice that received the concept diet containing MFGM coated around large lipid droplets (Lcoating, Nuturis®) showed significantly reduced FM compared with a diet containing small lipid droplets with comparable MFGM surface coating (P=0·05; Fig. 2; analysed using LSD) and a trend when compared with a diet containing large lipid droplets only (FM gain P=0·072 and %FM P=0·059, respectively).
BW was not significantly different between the four groups, although trends supported the above-summarised differences: BW (F=2161, P=0·157), LBM (F=1·074, P=0·312), FM (F=3·620, P=0·072) and %FM (F=3·646, P=0·071) at PN98 (analysed by a mixed effect regression model).
Programming effects of dietary lipid structure in early-life diet on adult organ weights
The weight of the EPI fat depot was lower in mice fed Lcoating in early life compared with small lipid droplets with or without MFGM surface coating (S, P=0·045 and Scoating, P=0·05, respectively; Table 3) and a trend was observed when compared with large lipid droplets only (L; P=0·071). The ING fat depot weight was lower in Lcoating mice (Nuturis®) compared with Scoating (P=0·018; analysed using LSD) and a trend was observed when compared with large lipid droplets only (L; P=0·066). The weight of the peri fat depot was decreased in Lcoating mice compared with L (P=0·049) and a trend was detected when compared with Scoating (P=0·053). The weight of the liver was higher in Scoating and L compared with S (P=0·047 and P=0·021, respectively) and a trend was observed when compared with Lcoating (P=0·073).
S, small droplets; L, large droplets; Scoating, small droplets coated by milk fat globule membrane (MFGM); Lcoating, large droplets coated by MFGM.
a,b Mean values with unlike superscript letters were significantly different (P<0·05).
Effects of dietary lipid structure in early-life diet on adult plasma parameters
At day 98, plasma parameters were only moderately affected by the WSD itself, as indicated by the results of the REF group: only glucose concentrations were 44 % lower in REF mice fed a standard rodent diet compared with the groups fed WSD (P<0·05, data not shown). The multivariate analysis results are depicted in the supplements (online Supplementary Table S1; analysed using LSD).
Programming of markers for adiposity, adipocyte size and cell differentiation
At day 98, leptin expression was significantly lower in the Lcoating group compared with the S control group (P=0·023; Fig. 3(a); analysed using LSD). Besides, mest/peg1 expression in Lcoating (Nuturis®) was decreased compared with L following the WSD challenge (P=0·043; Fig. 3(b)), but only a trend was seen for lower mest/peg1 expression in Lcoating compared with Scoating and S (P=0·053 and P=0·051, respectively).
Discussion
The present study aimed to determine whether providing a diet in early life containing large lipid droplets, MFGM as lipid droplet coating or their combination, was responsible for the long-term protective effect against adiposity in a mouse model. Our results clearly indicate that the combination of a large droplet size and MFGM present as a lipid droplet coating was essential for prevention of excessive body FM accumulation in adult life. Neither an increase in droplet size alone nor MGFM coating of small lipid droplets exerted any beneficial programming of adult body composition. These data suggest that the complex structure of large lipid droplet size consisting of a TAG core surrounded by a trimembrane layer as present in HM may significantly contribute to its reported protective effect against obesity in later life. These outcomes support a role of lipid matrix beyond the FA composition in dietary lipid quality in early life linked to later life health.
The mechanism(s) underlying the reduction in fat accumulation by early feeding with the Nuturis® diet is not completely clear. As food intake between groups was comparable, differences in food intake cannot explain the observed adult phenotype. The structural aspects of the Nuturis® concept, which bring its nutrient matrix closer to the HM structure, may affect digestion and absorption, and thereby modulate the postprandial response and FA metabolism( Reference Fave, Coste and Armand 16 ). We hypothesise that differences in postprandial lipid handling may programme metabolic homoeostasis and energy balance. This could impact the development and function of metabolic organs including intestine, pancreas, liver, brain and WAT, and thereby influence later life metabolic health. Previous studies have indeed shown that the lipid structure impacts gastric emptying, lipid absorption and digestion kinetics as well as the postprandial handling of lipids( Reference Armand, Hamosh and Mehta 21 , Reference Armand, Pasquier and Andre 22 ). For example, in adult rats, HM-derived lipid droplets coated by MFGM alters the TAG appearance compared with small droplets with proteins at the interface( Reference Michalski, Soares and Lopez 28 ). It clearly shows the importance of lipid droplet size as well as the nature of the surface coating, both impacting physical digestion and absorption rate. In general, lipids will either be targeted for β-oxidation or stored in WAT depending on these differences in kinetics and metabolic fate. Clearly, in obese people, lipid handling has been favoured towards lipid storage instead rather than towards β-oxidation( Reference Bourlieu and Michalski 29 ). Absorption and digestion kinetics are different between adults and infants related to the digestive immaturity of the infant( Reference Abrahamse, Minekus and van Aken 30 , Reference Michalski, Genot and Gayet 31 ). The latter provides the possibility for lipids to be a major source of energy as well as education for the developing metabolism, affecting growth and metabolic development. We speculate that in particular the combination of lipid droplet size and its surface composition of PL originating from MFGM fraction as present in Nuturis® during infancy may impact metabolic development and programme later life (metabolic) health. As a consequence, the adult metabolic homoeostasis may be more flexible to handle various lifestyle challenges, for example, a diet higher in (SFA) fats. Fat oxidation could be one of the targets of programming, but others are also possible, including glucose homoeostasis, fat absorption kinetics and postprandial physiology.
Future research is required to assess whether differential digestion and absorption kinetics induced by lipid structure of Nuturis® (Lcoating) indeed underlie the long-term effects found in the present study. Recent imaging experiments comparing Nuturis® lipid droplets indeed confirm the notion that the lipid structure is closer to that observed in HM (Gallier et al., unpublished results).
One of the possible mechanisms that could contribute to the observed protective effect of Nuturis® on body fat accumulation is reduced lipid storage capacity and a decline in lipogenesis in WAT. Previously, we found a reduced expression of PPAR γ (related to lipogenesis) in mice fed Nuturis® compared with mice fed control diet (S) before the WSD challenge( Reference Oosting, Kegler and Wopereis 23 ). Indeed, WAT development during the PN period may determine the size and function of WAT in later life( Reference Budge, Sebert and Sharkey 32 ). The total adipocyte number increases in early life and a constant fat cell number is only established by early adulthood( Reference Knittle, Timmers and Ginsberg-Fellner 33 , Reference Spalding, Arner and Westermark 34 ). FM changes in adulthood are mostly driven by changes in adipocyte cell volume and not cell number( Reference Spalding, Arner and Westermark 34 ). Mest/peg1, a marker for adipocyte size, is highly expressed in WAT and inhibits the Wnt signalling pathway, which stimulates the differentiation of stem cells in preadipocytes( Reference Spalding, Arner and Westermark 34 – Reference Laudes 38 ). We found that early-life feeding with Nuturis® could lower the expression of mest/peg1 in RP WAT, although not significant compared with our control diet containing small lipid droplets with or without MFGM coating. This observation indeed supports the reduced FM seen in Nuturis®-fed mice. FM also strongly corresponded with leptin expression levels in RP WAT, particularly in Nuturis®-fed mice (Lcoating; R 0·732, P=0·000). Plasma leptin concentrations are positively correlated to FM in both humans and rodents and implicated in energy homoeostasis( Reference Maffei, Halaas and Ravussin 39 – Reference Guo, Halo and Leibel 41 ). Detailed analysis of adipose tissue to gain further insights into adipose tissue function is one of the approaches needed to understand the underlying mechanisms for the observed protective effect of Nuturis® feeding in early life. The expression of genes implicated in lipogenesis and lipolysis cannot be easily differentiated from the present study design, because the observed effect is a combination between programming and a WSD effect. Rather, a study with a specific experimental design is required to obtain mechanistic information on lipogenesis and lipolysis. Studies to obtain more detailed information on lipogenesis and lipid storage in WAT as well as in vivo energy expenditure are needed to determine whether these possibilities are indeed the underlying mechanisms.
MFGM supplementation to IMF has been proven to be safe and well tolerated( Reference Billeaud, Puccio and Saliba 42 ). Recent studies by Timby et al.( Reference Timby, Hernell and Vaarala 43 ) suggest beneficial effects of adding MFGM lipid products to IMF in infants on specific immune and cognitive outcomes. Cognitive development and serum cholesterol concentrations in MFGM-enriched formula-fed infants were closer to breast-fed infants at 12 months of age( Reference Timby, Domellöf and Hernell 19 , Reference Timby, Lönnerdal and Hernell 44 ). Although it was hypothesised that these effects could have sustained effects on adult cholesterol metabolism, as well as on the risk for CVD( Reference Timby, Lönnerdal and Hernell 20 ). No effects were reported on growth or long-term health outcomes. In a recently published editorial, Greer & Kleinman( Reference Greer and Kleinman 45 ) raised the question whether adding MFGM combined with lower protein and energy contents to IMF could bring IMF closer to HM composition and health benefits attributed to breast-feeding. Our results clearly show that adding MFGM as a coating to current IMF containing small lipid droplets, only the combination of a large lipid droplet size and the MFGM coating showed beneficial effects on later life FM gain. Notably, in the studies by Timby et al. ( Reference Timby, Domellöf and Hernell 19 , Reference Timby, Lönnerdal and Hernell 20 , Reference Timby, Hernell and Vaarala 43 ) MFGM was added to an IMF that was lower in both energy as well as protein compared with the IMF used to construct the rodent diets. Therefore, the observed effects found in Timby et al. could also, at least in part, be explained by the combination of these nutrient composition and energy differences instead of the MFGM supplementation alone. Clinical studies testing the effect of Nuturis® on growth and body composition development are currently on-going (clinical trial registry NCT01609634 and the Dutch Trial Register NTR3683).
Although we found reduced adiposity in adult life after feeding a diet with an altered lipid structure in early life, we did not observe significant improvements in metabolic parameters as far as tested. This could be explained, at least in part, by the fact that the WSD challenge was rather mild as it did not generate major disturbances in plasma lipid or glucose homoeostasis compared with the unchallenged REF group, nor resulted in hyperphagia or overt obesity. The choice of this relatively mild obesogenic diet was based on the assumption that this would mimick a modern Western lifestyle more closely. Indeed, the WSD challenge in our model did not induce hyperphagia as is commonly seen with more extreme WSD diets in the literature( Reference Wei, Sun and Chen 46 , Reference Duval, Thissen and Keshtkar 47 ). To assess whether the Nuturis® diet may protect against the development of metabolic risk factors, it would be interesting to extend the duration of the mild WSD challenge and/or use a more severe dietary challenge in future studies.
In conclusion, a diet containing the combination of large lipid droplets and a coating with MFGM mimicking the lipid droplet architecture as present in mammalian milk prevents body fat accumulation when challenged with a moderate WSD during adolescence and adulthood. These observations confirm the notion that the complex lipid structure beyond the composition of lipids may be an important factor in the early diet that determines growth and development of body composition and impacts long-term (metabolic) health. Future experimental and clinical studies are required to confirm the protective effect of this concept in humans and to gain more detailed insights about the underlying mechanisms.
Acknowledgements
This study was not supported by external funding. The study was funded by Danone Nutricia Research. A. B., A. O., E. E., D. K., A. K., L. S. and E. M. v. d. B. are employed by Danone Nutricia Research and contributed to the study design, conduct of the study, analysis of the samples and data, interpretation of the findings and preparation of the manuscript. H. J. V. is a consultant for Danone research outside the submitted work, for which his institution (Department of Paediatric Gastroenterology and Hepatology, University Medical Centre Groningen, The Netherlands) is compensated financially.
The authors’ contributions are as follows: A. O. and E. M. v. d. B. designed the research; L. S., E. E. and D. K. conducted the research; A. B., A. O., L. S. and E. E. analysed the data, A. B. wrote the paper, A. O., E. M. v. d. B. and H. J. V. reviewed the manuscript. All the authors read and approved the final version of the manuscript.
Supplementary Material
For supplementary material/s referred to in this article, please visit http://dx.doi.org/doi:10.1017/S0007114516001082