There is consistent evidence that n-6 and n-3 long-chain (LC) PUFA are important to ensure optimal growth and development, particularly during fetal and early postnatal life( Reference Escolano-Margarit, Campoy and Ramirez-Tortosa 1 , Reference Harslof, Larsen and Ritz 2 ).
Desaturation/elongation from essentials linoleic (LA, 18 : 2n-6) and α-linolenic (ALA, C18 : 3n-3) fatty acids (FA) facilitate the synthesis of important LC-PUFA, including arachidonic acid (AA, 20 : 4n-6), EPA (20 : 5n-3), docosapentaenoic acid (DPA, 22 : 5n-3), and DHA (22 : 6n-3)( Reference Koletzko, Lattka and Zeilinger 3 – Reference Park, Kothapalli and Reardon 5 ). However, the enzymatic machinery involved in this process, including Δ5 and Δ6-desaturases (encoded by fatty acid desaturase (FADS)1 and FADS2 genes, respectively) and elongases (encoded by elongation of very long chain fatty acids (ELOVL) gene), participates in both pathways, triggering a competition between n-3 and n-6 series( Reference de la Garza, Montes and Chisaguano Tonato 6 ). Consequently, inter-individual differences in capacity to form LC-PUFA may depend on the genotype of FADS1 and FADS2 ( Reference Yeates, Love and Engstrom 7 – Reference Lattka, Koletzko and Zeilinger 9 ), which build a gene cluster on chromosome 11 together with a third desaturase gene, FADS3 ( Reference Lattka, Illig and Koletzko 10 , Reference Reardon, Hsieh and Park 11 ).
To date, studies have focused on the effects of FA supplementation during pregnancy on body composition, visual and cognitive development( Reference Campoy, Escolano-Margarit and Anjos 12 ), as well as effects on immune outcomes and cardiovascular function in early infancy( Reference Drewery, Gaitan and Spedale 13 , Reference Rodriguez-Santana, Ochoa and Lara-Villoslada 14 ). Moreover, other authors suggest that fish-oil (FO) or n-3 FA supplementation during pregnancy show a potential role on essential FA status in infants( Reference Rodriguez-Santana, Ochoa and Lara-Villoslada 14 , Reference Glaser, Heinrich and Koletzko 15 ). However, the long-term impact of maternal PUFA status on functionality of the metabolic pathways involved in PUFA desaturation/elongation in children still remains unknown. Interestingly, several studies have demonstrated the association between variants in the human genes FADS1, FADS2 and FADS3 and blood levels of PUFA in European or Caucasian populations( Reference Harslof, Larsen and Ritz 2 , Reference Koletzko, Lattka and Zeilinger 3 , Reference Glaser, Heinrich and Koletzko 15 – Reference Meldrum, Li and Zhang 18 ). It seems that genetic variation in the FADS gene cluster could modify desaturase function supporting that minor alleles increase the concentrations of LA and ALA, and decrease the concentrations of AA, EPA, DPA and DHA in serum, plasma and erythrocytes. Some authors support that 28·5 % of the variability in both PUFA and LC-PUFA levels in human tissues is determined by FADS gene variants( Reference Schaeffer, Gohlke and Muller 19 ). However, diets and genetic population background have also a potential role on associations of PUFA levels with different FADS SNP( Reference Gonzalez-Casanova, Rzehak and Stein 20 , Reference Tanjung, Rzehak and Mansyur 21 ). Glaser et al. ( Reference Glaser, Heinrich and Koletzko 15 ) found that FADS1 and FADS2 gene variants modulate tracking of serum PUFA levels in glycerophospholipids (GPL), which are the main lipid constituents of cell membranes with key functions in cell structure, intracellular signalling, anchoring to membranes and stabilisation of protein structure.
With these considerations in mind, we analysed the cheek cell GPL-FA composition in children from the Nutraceuticals for a Healthier Life (NUHEAL) study( Reference Campoy, Escolano-Margarit and Ramos 22 , Reference Krauss-Etschmann, Shadid and Campoy 23 ) at school age with the following aims:
(1) To analyse if maternal supplementation with DHA and/or 5-methyltetrahydrofolate (5-MTHF) during pregnancy exerted long-term effects on the LC-PUFA concentrations in their offspring cheek cell GPL at 8, 9 and 9·5 years.
(2) To evaluate whether FADS1, FADS2 and FADS3 genotype polymorphisms may modify PUFA concentrations in schoolchildren’s cheek cells.
(3) To explore the possible interaction between the described intervention in the mothers during pregnancy and the offspring FADS polymorphisms on children’s LC-PUFA concentrations in cheek cell GPL.
Methods
Participants
This study is a part of the European randomised multicentre trial NUHEAL, registered at www.clinicaltrials.gov, reference no. NCT01180933. A total of 270 women took part in the study until giving birth; the detailed study design, subject recruitment and population characteristics have been described elsewhere( Reference Campoy, Escolano-Margarit and Ramos 22 , Reference Krauss-Etschmann, Shadid and Campoy 23 ). Briefly, NUHEAL is a multicentre randomised study in healthy pregnant women from Munich, Pécs and Granada, designed to assess the effects of receiving daily supplements of FO (500 mg of DHA+150 mg EPA), 5-MTHF (400 µg), placebo or both (FO+5-MTHF), from gestation week 20 until delivery, on pregnancy outcomes and on offspring long-term growth and development. Study participants were approached again and asked to participate in a new assessment of their children; 147 complied with the request and cheek cell samples were obtained in all of them at 8, 9 and 9·5 years of age. In all, forty-two (28·6 %) were revisited at the Ludwig-Maximilians-University of Munich, Germany, and 105 (71·4 %) at the University of Granada, Spain. This follow-up examination in the children was performed under the NUTRIMENTHE EU Project framework. Dropout rates were similar between intervention groups, mainly for relocation (n 3), loss of contact (n 75) and unwillingness to continue (n 38). In all, four children were born prematurely before the 35th week of pregnancy, one child was born with a congenital left-side anophthalmus, one child developed craniosynostosis, and another was reported to have left-side deafness. There were no differences in maternal mean dietary intake of DHA between intervention groups at the 20th and 30th weeks of pregnancy. Certain baseline characteristics, including maternal parity, BMI, smoking habits, obstetric risk factors, complications at parturition or perinatal morbidity in any intervention group, did not differ between the studied and not-tested populations. Children whose fathers had a high educational level in the placebo and FO+5-MTHF groups showed high drop-out rates (77·9 % of children with fathers with general qualification for university entrance in the placebo group and 72·8 % in the FO+5-MTHF group withdrew at 8 year follow-up, whereas only 41 % of the fathers of followed children in the placebo group and 32·6 % in the FO+5-MTHF group had this qualification (P=0·004)).
The study protocols were approved by the Biomedical Ethical Committees from both participating centres. Written informed consent was obtained from all participants.
Cheek cell sampling procedure
The sampling procedure consisted of mouth cleaning with 200 ml tap water three times before sampling, followed by rubbing of each inner cheek side with a Bio Brush (reference 73506013; Labolan). Subsequently, the mouth was rinsed with 10 ml pure water, and the rinsing solution containing cells washed off the inner cheeks was collected in a tube. The brush used for scraping was inserted in this tube for a few seconds and, after removing the brush, this tube was closed and centrifuged at 1400 g (2710 rpm) for 10 min at 4°C. The supernatant was discarded. The cell pellet was stored at −80°C until further analysis( Reference Klingler, Demmelmair and Koletzko 24 ).
Fatty acid measurement from cheek cell glycerophospholipids
Cheek cell GPL-FA were analysed as described by Klingler et al. ( Reference Klingler, Demmelmair and Koletzko 24 ) and results are presented as weight percentages of all quantified FA. The sample preparation procedure selectively transferred GPL bound FA into their methyl esters, which were quantified by GC. Quality control samples indicated on analytical precision depended on the percentage contribution of the corresponding FA to total GPL-FA and ranged from 2·29 % (LA) to 12·53 % (ALA) (CV). After quality control and plausibility checking, cheek cell GPL-FA percentages were compiled in a MS Excel Table.
FA used for the analyses as control were myristic acid (C14 : 0), palmitic acid (C16 : 0), margaric acid (C17 : 0), stearic acid (C18 : 0), arachidic acid (C20 : 0), behenic acid (C22 : 0), lignoceric acid (C24 : 0), palmitoleic acid (C16 : 1n-7), vaccenic acid (C18 : 1n-7), oleic acid (C18 : 1n-9), eicosaenoic acid (C20 : 1n-9), nervonic acid (C24 : 1n-9), LA (C18 : 2n-6), γ-linolenic acid (18 : 3n-6), dihomo-γ-linolenic acid (C20 : 3n-6), AA (C20 : 4n-6), DPA (C22 : 5n-6), ALA (C18 : 3n-3), EPA (C20 : 5n-3), DPA (C22 : 5n-3) and DHA (C22 : 6n-3).
SNP selection and genotyping
A total of seventeen SNP from the FADS1 (rs174548, rs174556 and rs174561), FADS2 (rs174570, rs174574, rs174575, rs174576, rs174578, rs174579, rs174602, rs498793, rs968567, rs2727271 and rs3834458) and FADS3 (rs174448, rs174449 and rs174455) gene cluster from the children’s umbilical cord were genotyped. Of these, eleven SNP were selected according to its localisation in FADS gene cluster on chromosome 11, while six additional SNP (rs174546, rs174561, rs3834458, rs174548, rs174574 and rs174578) were included due to its relation with FA levels previously reported( Reference Koletzko, Lattka and Zeilinger 3 , Reference Lattka, Koletzko and Zeilinger 9 ). The genotyping efficiency was 10·7-fold, although the two further SNP rs174561 and rs3834458 could not be included in the HapMap database. Genotyping of SNP was carried out with the iPLEX method (Sequenom) by means of matrix-assisted laser desorption ionisation–time of the flight MS method (Mass Array; Sequenom), according to the manufacturer’s instructions. Standard genotyping quality control included 10 % duplicate and negative samples. The genotyping discordance rate was below 0·3 %.
Statistical analyses
The power calculation was calculated by standard approach, setting α value as 0·05 and β value as 0·2. Different studied groups (FO, 5-MTHF, FO+5-MTHF and placebo) were compared with each other for essential FA, including LA, AA, ALA, EPA, n-5 DPA and DHA. A statistical power of 80 % for selected FA levels in children’s cheek cell GPL at school-age was obtained. Sample size required for our genetic study was also calculated using R package ‘powerEQTL version 0.1.3’, which considers different ranges of sample sizes, α value of 0·05 and seventeen selected SNP. A statistical power of 90 % was obtained for analysed SNP, except for rs174570 and rs2727271 whose was of 80 % (see online Supplementary Methods for more details).
Histograms of FA levels were studied and assessed for normality using Shapiro–Wilk tests with the help of SPSS statistical program, version 20.0 (SPSS Inc). Some LC-PUFA were extremely right skewed and thus were transformed using the natural logarithm (ln) and square root. If no normal distribution was achieved, non-parametric tests were used. Descriptive results were expressed as means and standard deviations. Categorical outcomes were assessed with χ 2 tests and continuous variables with Kruskal–Wallis/one-way ANOVA tests to explore the potential association between FADS polymorphisms and type of prenatal supplementation.
Hardy–Weinberg (H–W) equilibrium for the genotypes was calculated using R statistical software (3.2.2 version, ‘genetics’ package). Deviations from H–W equilibrium were studied using Fisher’s exact test. Further, to test the linkage disequilibrium, Lewontin’s D′ and pairwise squared correlations r 2 were performed.
After selecting the most representative SNP of our sample, the data were transformed into longitudinal form to perform a linear mixed model with the aim to study the relationships between FADS and FA concentrations in cheek cells. The outcomes were the FA concentrations, the fixed part was each one of the selected genotypes (homozygous major allele carriers were coded as 0, heterozygous as 1, and homozygous minor allele as 2), and the random effect was the time variable (8, 9 and 9·5 years).
Furthermore, several linear mixed models have been built to assess each FA considering as co-founders the eight selected SNP, group of maternal prenatal supplementation, maternal age at delivery, mother’s BMI at weeks 20 and 30 of pregnancy, mother smoking at weeks 20 and 30 of pregnancy, placental weight and sex. To select those co-founder variables, separate analyses were first performed for each of the eight candidate SNP. Then, linear mixed models were adjusted for all other analysed SNP and those covariates which changed between study groups or by SNP. Those co-founders which showed a P<0·05 in the preliminary analyses were subsequently included in our final adjusted linear mixed models.
Results
The baseline characteristics of the study participants are shown in online Supplementary Table S1. There were no significant differences between study groups with respect to sex of the children, maternal age, BMI, smoking at weeks 20 and 30 of pregnancy and placenta weight between the four groups of prenatal supplementation.
Prenatal supplementation and children fatty acids in cheek cell glycerophospholipids
Depending on mother’s supplementation during pregnancy, percentages of those GPL-FA in cheek cells measured in the NUHEAL children resulted statistically significant were presented in Table 1. Our analysis showed a statistical difference in the GLP AA concentrations in cheek cells from 8 to 9·5 year (P=0·002). As we can see in Fig. 1, levels of GPL AA significantly decrease from 8 to 9·5 year in those children whose mothers were supplemented with FO or FO+5-MTHF. In order to understand our results, we have also explored AA:DHA, AA:EPA and AA:[DHA+EPA] ratios, but not significant differences were observed. Moreover, there was also a decrease in AA:DHA ratio from 8 to 9·5 years which did not resulted significant (P=0·985).
Evolution of arachidonic acid (AA) concentrations in children’s cheek cell glycerophospholipids from 8 to 9·5 years of age, depending on their mothers’ supplementation during pregnancy. Significant decrease in AA concentrations were observed in those children born to mothers who received fish oil (FO) with or without 5-methyltetrahydrofolate (5-MTHF). ***P<0·001, ****P<0·0001 after Bonferroni post hoc correction.
, FO group;
, 5-MTHF group;
, FO + 5-MTHF group;
, placebo group.
Effects of prenatal supplementation on behenic, arachidonic, n-6 docosapentaenoic and oleic acid levels in children’s cheek cell glycerophospholipids at 8, 9 and 9·5 years old† (Mean values and standard deviations)
n, number of cases; P1, level of significance from one-way ANOVA analysis for type of prenatal supplementation in each time point; P2, level of significance from one-way ANOVA analysis between 8, 9 and 9·5 years old; 5-MTHF, 5-methyltetrahydrofolate; FO, fish oil.
* P-values <0·05.
† Percentage of total fatty acids from cheek cell glycerophospholipids.
At 8 years old, children born to mothers who received the FO+5-MTHF supplementation during pregnancy had higher concentrations of behenic (P=0·037) and n-6 DPA acids (P=0·007) in cheek cell GPL than those whose mothers were supplemented with FO, 5-MTHF or Placebo. At 9 years old, oleic acid concentrations in GPL measured in cheek cells were higher in those children born to mothers who received FO during pregnancy respect to those born to mothers from 5-MTHF or FO+5-MTHF groups (P=0·035). Concentrations of other analysed FA, including LA, ALA, EPA, n-3 DPA and DHA, did not show significant differences between the four study groups of children at any individual time point (online Supplementary Table S2).
Effects of fatty acid desaturase polymorphisms on fatty acid composition in cheek cells
Genotype and allele frequencies of seventeen polymorphisms in FADS1, FADS2 and FADS3 genes from NUHEAL children are presented in Table 2. The analysed genotypes were in H–W equilibrium, except for rs174574, rs174576 and rs174578 which were excluded from the following analyses.
Characteristics of the seventeen analysed variants in fatty acid desaturase (FADS) 1/2/3 genes of children in the studied population (Numbers and percentages)
0, homozygous major; 1, heterozygous; 2, homozygous minor; P, level of significance according to Fisher’s exact test.
* P-values <0·05.
Based on the linkage disequilibrium data from Lewontin’s D′ and pairwise-squared correlations r 2 (Fig. 2 and online Supplementary Table S3) eight SNP out of the initial seventeen SNP (rs174556, rs174570, rs174575, rs174579, rs174602, rs498793, rs2727271 and rs174448) were selected for further analyses. For these selected SNP, Lewontin’s D′ ranged between 0·95 and 1·0 and the pairwise-squared correlations r 2 ranged between 0·85 and 1·0.
(a) Pairwise linkage disequilibrium (LD) measured in D′ for the complete seventeen SNP obtained from Nutraceuticals for a Healthier Life (NUHEAL) children; (b) pairwise LD measured in r 2 for the complete seventeen SNP obtained from NUHEAL children.
We analysed the effects of FADS gene cluster polymorphisms on GPL-FA concentrations in children. Table 3 represents the regression analysis between the main seven GPL-FA in children cheek cells and the eight FADS gene cluster polymorphisms selected. The analysed FADS1 and FADS2 polymorphisms showed different associations with LA, ALA, AA, EPA, n-3 DPA and DHA concentrations. The majority of the associations were established between FADS1 and FADS2 SNP with LA and ALA. FADS3 rs174448 did not have any effects on the analysed FA concentrations.
Associations between PUFA concentrations in children’s cheek cell glycerophosphopholipids and selected SNP of fatty acid desaturase (FADS) gene cluster polymorphism (β-Coefficients and 95 % confidence intervals)
Int, intercept (constant from regression model) (mean fatty acid concentration in homozygous carriers of the major allele); β, regression coefficient (linear regression of the fatty acid (outcome) on the respective single SNP); P, level of significance obtained from linear regression.
* P-values <0·05.
The frequencies of major and minor alleles according to genotypes are presented in online Supplementary Table S4. Among all groups of prenatal supplementation, the major allele carriers present highest frequencies for FADS2 and FADS3. No statistical differences were found in the frequencies of major or minor homozygous and heterozygous for FADS1, FADS2 or FADS3 between children depending on their mothers groups of supplementation.
Next, we analysed the effects of SNP on FA concentration depending on the maternal supplementation groups (Table 4).
Relationship between fatty acid glycerophospholipid concentrations in children’s cheek cells and fatty acid desaturase (FADS) gene cluster polymorphisms by type of prenatal supplementation† (Number of samples, mean values and standard deviations)
5-MTHF, 5-methyltetrahydrofolate; FO, fish oil; 0, homozygous major; 1, heterozygous; 2, homozygous minor; P1, level of significance for each prenatal supplementation group considering homozygosis for each selected SNP; P2, level of significance between prenatal supplementation groups; LA, linoleic acid; AA, arachidonic acid; DPA, docosapentaenoic acid; ALA, α-linolenic acid.
* P-value <0·05.
† Random effect=time points (8, 9 and 9·5 years of age).
As shown in Table 4, LA children whose mothers were supplemented with FO, and carrying the minor homozygous of both FADS1 rs174556 and FADS2 rs174575 showed high LA concentrations respect to heterozygosis or major homozygous genotypes. High LA levels were also found in those children born to mothers who received placebo and carrying the heterozygous genotype for FADS1 rs174556, FADS2 rs174575 and rs174579. AA concentrations were increased in children with FADS2 rs174570 and rs174579 major homozygous genotypes in case of maternal FO or placebo supplementation, respectively. No other effects of the different polymorphisms or supplementation type were observed.
Prenatal 5-MTHF or FO+5-MTHF supplementation and specific FADS polymorphisms significantly affected n-6 DPA levels. In fact, those children whose mothers received prenatal 5-MTHF supplementation, and carrying the major homozygous genotypes for FADS3 rs174448, showed high n-6 DPA levels. Regarding prenatal FO+5-MTHF supplementation, we observed higher n-6 DPA levels in children with minor homozygous of FADS2 rs174575 and rs174602. Interestingly, higher n-6 DPA concentrations were found in those children born to mothers who received placebo and carrying the major homozygotic alleles for FADS2 rs174602 and rs2727271 polymorphisms.
Those children, whose mothers received FO supplementation during pregnancy, and having the minor homozygotic alleles for FADS1 rs174556 and FADS2 rs174575 had a higher level of ALA in cheek cell GPL, respect to those carrying the other polymorphisms. Furthermore, children carrying the heterozygous genotypes for rs174579, whose mothers were supplemented with 5-MTHF, showed an increase of ALA concentrations in cheek cell GPL. Higher ALA levels were also found in those children heterozygous for FADS2 rs174575 and rs174579 variants born to mothers who received placebo during pregnancy.
Those children heterozygous for the FADS2 rs174570 polymorphism, and whose mothers received placebo during pregnancy, had a lower concentration of EPA in cheek cell GPL, compared with those children carrying major homozygous; prenatal FO+5-MTHF supplementation was found to be associated with high levels of EPA in those heterozygous children carrying FADS2 rs2727271 or minor homozygosis in FADS2 rs174602 polymorphisms. Children born to mothers who were supplemented with FO had higher DPA concentrations in cheek cell GPL when they carried the major homozygous rs174579 and heterozygous rs2727271. Moreover, prenatal supplementation with FO+5-MTHF determined higher concentrations of DPA when children were heterozygous for FADS2 rs498793. On the other hand, high DHA levels are caused by major homozygous FADS2 rs174602 and heterozygous FADS1 rs174556 and FADS2 rs174575 in those children who mothers received prenatal FO supplementation. No definitive pattern was found in both DPA and DHA concentrations in those children born to mothers who received placebo; in this case, the presence of the different FADS polymorphisms had no effect on aforementioned FA levels.
Finally, we analysed the effects of SNP and different confounders on the FA concentrations in GPL from cheek cells using a lineal mixed model (Table 5). LA concentrations were influenced by rs174556, rs174602, rs498793 and rs174448 genotypes, child’s sex, maternal age and BMI 20 kg/m2 (P<0·0001). ALA concentrations in cheek cell GPL were associated with rs174570, rs174448 and sex. AA was influenced by all the SNP studied, except rs174602, rs498793 and rs174448. Moreover, we did not observe any association between n-6 DPA levels and selected confounders. Regarding n-3 series LC-PUFA, EPA shows a liaison with rs2727271, meanwhile n-3 DPA shows no correlation with any SNP, but with sex and maternal age. Different SNP influenced the DHA concentrations, but the model showed no statistical significance.
Associations between PUFA in glycerophospholipids from children’s cheek cell samples and selected SNP of fatty acid desaturase gene cluster polymorphism adjusted by different co-founders (β-Coefficients and 95 % confidence intervals)
CONS, constant value; MA, maternal age; BMI20, BMI at 20 weeks of pregnancy; BMI30, BMI at 30 weeks of pregnancy; S20, smoking at 20 weeks of pregnancy; S30, smoking at 30 weeks of pregnancy; PWD, placental weight at delivery; GS, group of prenatal supplementation; P, level of significance of each co-founder variable; P-value, unadjusted level of significance considering all analysed co-founders; P-value adjusted, level of significance considering selected co-founders; LA, linoleic acid; AA, arachidonic acid; DPA, docosapentaenoic acid; ALA, α-linoleic acid.
* P-value <0·05.
Discussion
This is the first study suggesting the presence of long-term effects of maternal prenatal supplementation with DHA and/or 5-MTHF on PUFA status and metabolism in the offspring at school age; furthermore, our data also seem to indicate an interaction between maternal supplementation during pregnancy and their offspring polymorphisms in the FADS cluster, determining children’s PUFA concentrations in cheek cell GPL at school age.
Previous results from NUHEAL study demonstrate clear associations between maternal FA concentrations in plasma phospholipids during pregnancy and those levels found in their offspring (umbilical cord)( Reference Escolano-Margarit, Campoy and Ramirez-Tortosa 1 ), as well as changes depending on maternal supplementation during pregnancy. However, the present analysis suggests long-term effects of maternal supplementation with FO or FO+5-MTHF during pregnancy on the FA status in cheek cell GPL in their offspring up to 9·5 years. In fact, supplementation with FO+5-MTHF during pregnancy determined an increase of children’s n-6 DPA concentrations in cheek cell GPL at 8 years old, suggesting an additive effect of both supplements promoting a major activity in the n-6 PUFA series. It is also important to notice that children’s AA and DHA concentrations measured in cheek cell GPL from 8 to 9·5 years are within normal ranges, but a significant decrease of AA concentrations was observed in those children whose mothers were supplemented with FO or FO+5-MTHF. Furthermore, there were no significant differences for AA:DHA, AA:EPA and AA:[DHA+EPA] ratios. Thus, our results suggest that maternal FO prenatal supplementation could determine a beneficial programming effect on the metabolic pathways during fetal life, which would be maintained up to late childhood; specifically, this supplement would program the metabolic n-6 FA pathways with a significant reduction of the synthesis of AA without affecting n-6 concentrations or aforementioned ratios. This apparent beneficial effect could be related to a lower risk of adverse health consequences in the offspring in later life. As our study is not able to demonstrate if this effect is linked to an excess of maternal DHA and EPA intake during pregnancy, our results should be interpreted with caution due to the sample size and the multiple confounding factors along 9·5 years of life, but still are of interest for future studies. In this regard, a recent experimental study in newborn and recently weaned offspring mice supported that changes in specific hepatic microRNA (miRNA) (miR-122 and miR-370, respectively) are independently determined by an excess of maternal lipids during both gestational and lactation periods, leading in turn to impaired lipid metabolism in the offspring( Reference de Paula Simino, de Fante and Figueiredo 25 ). Finally, we suggest that this decline may also be associated with slightly modifications in dietary intake of AA, LA or n-3 LC-PUFA, which have a role in increasing or decreasing the AA concentrations. In fact, we observed an increase of n-3 LC-PUFA intake from 8 to 9 years in both German and Spanish children (data not shown).
Prenatal supplementation with FO alone also determined higher oleic acid concentrations in GPL in the offspring at 9·5 years old, respect to those children born to mothers who received 5-MTHF. Again the interaction between FO and 5-MTHF seems to have long-term effects on the biochemical pathways involved in FA and folate metabolisms.
Our study also tried to evaluate whether FADS1, FADS2 and FADS3 genotype polymorphisms may modify PUFA concentrations in schoolchildren’s cheek cells. To date, findings support that 28·5 % of the variability observed in both PUFA and LC-PUFA levels in human tissues are associated with FADS gene variants( Reference Escolano-Margarit, Campoy and Ramirez-Tortosa 1 , Reference Schaeffer, Gohlke and Muller 19 ); therefore, it is clear that both diet and genetic variation play an important role in determining blood and tissue levels of LA and ALA, as well as their biologically active LC-PUFA derivatives. In fact, this is the first study suggesting the associations of genetic variants in the FADS1, FADS2 and, to a lesser extent, FADS3 gene cluster and children’s cheek cell GPL-FA. Interestingly, our study results seem to indicate a long-term interaction between maternal supplementation during pregnancy and FADS genotypes that extends up to 9·5 years of age; the effect of such potential interaction may determine n-3 and n-6 LC-PUFA concentrations in childhood. Maternal supplementation with FO alone or plus 5-MTHF was associated with both n-3 and n-6 LC-PUFA status in their offspring cheek cell GPL; however, supplementation only with FO or 5-MTHF, upon a possible interaction with FADS polymorphisms mainly causes an increase in the concentration of ALA. This pattern is also seen in infants born to mothers supplemented with placebo, except for n-3 DPA and DHA levels. In the Avon Longitudinal Study of Parents and Children (ALSPAC), child’s FADS genotypes influenced amounts of cord plasma LC-PUFA and FA ratios; furthermore, the authors suggested a specific effect of the minor allele on FA metabolism( Reference Lattka, Koletzko and Zeilinger 9 ).
The present data suggest that maternal FO supplementation in association with minor homozygosis for FADS1 rs174556 and FADS2 rs174575 determine a decrease of FADS activity and subsequent increase of n-3 and n-6 FA precursors. However, other FADS2 polymorphisms seem to be associated with significant decrease of AA, n-3 DPA and DHA in cheek cell GPL. Interestingly, 5-MTHF supplementation during pregnancy did not show any effect on desaturation and elongation of both n-3 and n-6 FA in the offspring, independently of the FADS gene polymorphisms; only, major homozygosis of FADS3 rs174448 was associated with higher n-6 DPA and heterozygosis of FADS2 rs174579 with higher ALA concentrations. To date, it is reported that folic acid improves the re-methylation of homocysteine, which leads to the formation of methionine. Some authors( Reference Sugiyama, Kumazawa and Zhou 26 ) have suggested that methionine stimulates phosphatidylethanolamine methylation, thereby altering the ratio of phosphatidylcholine to phosphatidylethanolamine in liver microsomes; as consequence, Δ5 and Δ6-acyl-coenzyme A desaturases increase their enzymatic activity, which in turn results in higher LC-PUFA concentrations. Our results show that the supplement combining FO+5-MTHF during pregnancy was not associated with the same effect than FO alone on cheek cell GPL-FA in the offspring up to 9·5 years; on the contrary, changes in allele frequencies in FADS2 in the FO+5-MTHF group seem to be associated with an increase in n-6 DPA, n-3 DPA and EPA. With these results we are just speculating that there seems to be an interaction between early life supplementation and offspring FADS polymorphisms determining FA status at childhood. Consequently, future studies are needed to clarify this interaction, evaluating also the potential effects of ELOVL gene polymorphisms, together with FADS gene cluster, to disentangle all actors involved in the potential early programming of LC-PUFA metabolic pathways.
It is important to highlight that FA status is largely based on analysis of FA levels in plasma, erythrocytes and adipose tissue. These approaches, however, are invasive and not well accepted in population studies, particularly in those including infants and young children. In the present study, we measured FA in cheek cells as reference of FA status. Despite its limitations in sample quality and quantity, buccal mucosal cheek cell phospholipids should be proposed as a non-invasive biological marker for LC-PUFA status in children( Reference Klingler, Demmelmair and Koletzko 24 ). In fact, several studies have previously reported that cheek cell FA composition not only reflects dietary FA intake, but also both erythrocyte and plasma levels of essential PUFA( Reference Hoffman, Birch and Birch 27 – Reference Kirby, Woodward and Jackson 30 ). In this regard, we found that daily intake of AA, EPA and folate increase the levels of AA, EPA and n-5 DPA in GPL of children’s cheek cells, while DHA intake decreases the concentration of n-5 DPA (data not shown). Furthermore, several authors had already demonstrated the association between variants in the FADS1, FADS2, FADS3 genes and blood levels of PUFA( Reference Harslof, Larsen and Ritz 2 , Reference Koletzko, Lattka and Zeilinger 3 , Reference Glaser, Heinrich and Koletzko 15 – Reference Meldrum, Li and Zhang 18 ). On one side, PUFA concentrations were strongly associated with higher or lower prevalence of several diseases, including the metabolic syndrome( Reference Guo, Li and Shi 31 ), obesity( Reference Simopoulos 32 , Reference Vaittinen, Mannisto and Kakela 33 ), CVD( Reference Endo and Arita 34 , Reference Takahashi, Ando and Shimada 35 ), atopic sensitisations( Reference Miles and Calder 36 , Reference Rucci, den Dekker and de Jongste 37 ) or rheumatoid arthritis( Reference Brouwers, von and Toes 38 – Reference Zhou, Chen and Hu 40 ). On the other side, FADS polymorphisms have not only shown potential effects on neuropsychological and developmental outcomes in infants( Reference Lauritzen, Sorensen and Harslof 17 , Reference Andersen, Harslof and Schnurr 41 ), but also on FA related disorders, including attention-deficit/hyperactivity disorder, allergic diseases and the metabolic syndrome( Reference Meldrum, Li and Zhang 18 ). With these considerations in mind, we suggest that the profile of FA in cheek cell GPL would be an easy method in children to study whether interaction between FA status and the FADS polymorphisms may be related to the development or evolution of some diseases.
Besides, LC-PUFA status can be considered as a biomarker of the effect of supplementation during pregnancy, which proves again that there is a long-term effect of early nutrition in the offspring( Reference Escolano-Margarit, Campoy and Ramirez-Tortosa 1 ). We demonstrated that prenatal supplementation with FO and 5-MTHF has a long-term effect on cognitive (attention) development in the offspring at 8 years( Reference Catena, Muñoz-Machicao and Torres-Espinola 42 ), but we should seriously consider that may be a potential impact on neurodevelopment depending on children’s own lipid metabolism. Moreover, other authors noticed that specific PUFA may have also general health benefits( Reference Lattka, Klopp and Demmelmair 16 , Reference Chisaguano, Montes and Perez-Berezo 43 ), but may depend on the individual FADS genotypes. So, it will be very interesting to perform further studies to evaluate this hypothesis in different fields of children’s health.
Finally, we believe it is important to mention that the results of our current study have an added value, because to our knowledge, no study has made such a long-term follow-up after a prenatal supplementation. However, in the current study there are limitations that must be acknowledged. First, higher number of cases per group would have been better, but the sample size achieved resulted high enough to detect relevant FA concentrations differences between groups. Furthermore, we did not see a positive or negative trend in our data. Although sample size for genetic study was about 80–90 %, we are all aware that larger sample would have led us to find more compelling results. Having in mind these considerations, we think that the current study is still interesting as a pilot to opening new perspectives focused on the role of both genotype and early nutrition on the long-term development of FA metabolism.
Secondly, the effects of FO and 5-MTHF supplementation during pregnancy on the obtained results at school age should be interpreted with caution, due to the short window of data collection. Current international recommendations for pregnant women indicate that they must achieve at least 200–300 mg of DHA/d( Reference Koletzko, Cetin and Brenna 44 , Reference Koletzko, Lien and Agostoni 45 ). Considering that NUHEAL study included two groups of pregnant women who received daily 500 mg DHA+150 mg EPA (FO) from week 20 until delivery, and the suggested joint action between FO and 5-MTHF supplements, the potential excess of LC-PUFA intake in some women would have long-term unknown consequences. Taken into account these results, we propose that those women who follow a diet rich in both DHA and folic acid should receive an individualised supplement based on their genetic characteristics and real needs, avoiding the excess of these nutrients. In this regard, it is possible that both DHA and folic acid show specific critical windows during pregnancy which do not necessarily coincide. Thus, future studies are needed to further investigate the potential benefits of both DHA and folic acid supplementation at different times of gestation, avoiding overdoses and interaction between both nutrients to favouring fetal development without potential long-term adverse effects.
On the other side, children were programmed to be examined twice within NUTRIMENTHE framework, when they were aged 7·5 and 9 years; due to the long-time required to fulfil the protocol agreed, each assessment was divided into two visits resulting in four visits at 7·5–8 and 9–9·5 years old. So, finally we took cheek cell samples four times during 2 years, coinciding with the four visits performed; these repeated samples’ collection improved the accuracy and reliability of this new methodology. However, cheek cell samples were obtained at 7·5 years with a different type of brush, and therefore were discarded.
In summary, herewith we report that FO or FO+5-MTHF supplementation during pregnancy seems to have long-term effects on the GPL-FA profile in the cheek cells of the offspring. Most of the FADS polymorphisms studied are associated with PUFA measured in the NUHEAL children’s cheek cells. Together, data obtained from this study suggest that long-term FA status during childhood is not only related to early nutrition programming of FA metabolic pathways, but also FADS genetic variations in children. Further research is guaranteed to demonstrate, in a more plausible evolutionary design, the long-term effects of early life nutritional interventions depending on individual FADS genetic variations.
Acknowledgements
The authors are grateful to Pilar Brandi and Gudrun Haile for their collaboration in contacting the families, administering the questionnaires and collecting the samples. The authors thank all people involved in the huge number of tasks of this study. Moreover, the authors thank the participants for their generous collaboration for such a long time and the heads of the institutions (German, Hungarian and Spanish) to ideate this wonderful project, NUHEAL.
This work was supported by the Commission of the European Community’s Seventh Framework Program (FP7/2008–2013), grant agreement no. 212652 (NUTRIMENTHE Project), Marie Curie post-doctoral fellowship (FP7, no. 329812, NutriOmics); within the Sixtth Framework Program, contract no. 007036 (EARNEST Project); and supported in part by the Commission of the European Community within the Fifth Framework Program, contract no. QLK1-CT-1999-00888 (NUHEAL EU project). This publication is the work of the authors and does not necessarily reflect the views of the Commission of the European Community. The work of B. K. is supported by the European Research Council Advanced Grant ERC-2012-AdG – no. 322605 META-GROWTH. Funders had no role in the design, analysis or writing of this article.
All authors helped in the interpretation of results and contributed to manuscript preparation. C. C. coordinated and designed the study, was the coordinator of the study at University of Granada and wrote the manuscript; C. M.-Z. and H. D. performed the cheek cell analysis; C. M.-Z. and J. A. G.-S. wrote the manuscript. E. R. and P. R. have been responsible for the FADS polymorphisms analysis. C. C., C. M.-Z., S. A. and H. A. worked on the data analysis and interpretation; B. K. was the coordinator of the NUHEAL study in Germany and supervised the manuscript.
The results reported in this article are likely to be included in the Doctoral Thesis of H.A. at the University of Granada, PhD Program in Epidemiology and Public Health.
The authors declare that there are no conflicts of interest.
Supplementary materials
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114518002623
Open access
