Figures
Abstract
Introduction
Breastfeeding effects on cognition are attributed to long-chain polyunsaturated fatty acids (LC-PUFAs), but controversy persists. Genetic variation in fatty acid desaturase (FADS) and elongase (ELOVL) enzymes has been overlooked when studying the effects of LC-PUFAs supply on cognition. We aimed to: 1) to determine whether maternal genetic variants in the FADS cluster and ELOVL genes contribute to differences in LC-PUFA levels in colostrum; 2) to analyze whether these maternal variants are related to child cognition; and 3) to assess whether children's variants modify breastfeeding effects on cognition.
Methods
Data come from two population-based birth cohorts (n = 400 mother-child pairs from INMA-Sabadell; and n = 340 children from INMA-Menorca). LC-PUFAs were measured in 270 colostrum samples from INMA-Sabadell. Tag SNPs were genotyped both in mothers and children (13 in the FADS cluster, 6 in ELOVL2, and 7 in ELOVL5). Child cognition was assessed at 14 mo and 4 y using the Bayley Scales of Infant Development and the McCarthy Scales of Children's Abilities, respectively.
Results
Children of mothers carrying genetic variants associated with lower FADS1 activity (regulating AA and EPA synthesis), higher FADS2 activity (regulating DHA synthesis), and with higher EPA/AA and DHA/AA ratios in colostrum showed a significant advantage in cognition at 14 mo (3.5 to 5.3 points). Not being breastfed conferred an 8- to 9-point disadvantage in cognition among children GG homozygote for rs174468 (low FADS1 activity) but not among those with the A allele. Moreover, not being breastfed resulted in a disadvantage in cognition (5 to 8 points) among children CC homozygote for rs2397142 (low ELOVL5 activity), but not among those carrying the G allele.
Citation: Morales E, Bustamante M, Gonzalez JR, Guxens M, Torrent M, Mendez M, et al. (2011) Genetic Variants of the FADS Gene Cluster and ELOVL Gene Family, Colostrums LC-PUFA Levels, Breastfeeding, and Child Cognition. PLoS ONE 6(2): e17181. https://doi.org/10.1371/journal.pone.0017181
Editor: Maria Moran, Hospital Universitario 12 de Octubre, Spain
Received: September 29, 2010; Accepted: January 22, 2011; Published: February 23, 2011
Copyright: © 2011 Morales et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was funded by grants from Instituto de Salud Carlos III (Red INMA G03/176, CB06/02/0041), Spanish Ministry of Health (FIS-PI041436, FIS- PI081151, 97/0588, 00/0021-2, PI061756 and PS0901958), Generalitat de Catalunya-CIRIT 1999SGR 00241, Beca de la IV convocatoria de Ayudas a la Investigación en Enfermedades Neurodegenerativas de La Caixa, EC Contract No. QLK4-CT-2000-00263, Fundación Roger Torné, and Fundació La Marató de TV3 (Project No. 090430). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Brain development depends on the combined effects of incompletely understood genetic and environmental factors. The majority of dry weight in an adult brain is composed of lipids, 35% of which are long chain polyunsaturated fatty acids (LC-PUFAs) [1]. During development the brain is more highly enriched than most other tissues in LC-PUFAs that are essential for efficient neurogenesis, myelination, neurite outgrowth, dendritic arborisation and neurotransmission [2]–[7]. During pregnancy, the fetus is supplied with preformed maternal LC-PUFAs by placental transfer. After birth breast milk provides a unique supply of crucial LC-PUFAs including eicosapentaenoic acid (EPA), arachidonic acid (AA) and docosahexanoic acid (DHA), which support the accretion of LC-PUFA in the brain growth spurt that takes place during the third trimester of pregnancy and the first two years of life [8], [9]. Higher concentrations of DHA and AA, rapidly incorporated in the developing brain, are found in infants who are breastfed in comparison with infants fed with unsupplemented formulas [10]. However, it remains controversial whether breastfeeding provides a direct nutritional advantage for child cognition– attributed in part to the high levels and range of LC-PUFAs– or simply reflects more favorable socio-environmental factors– no fully controlled in statistical analysis– among women who breast feed [11]–[14].
LC-PUFA synthesis is controlled by key enzymes (desaturases and elongases) encoded by the FADS gene cluster and the ELOVL gene family (see Supporting Information Figure S1). Fatty acid desaturases are encoded by FADS1 and FADS2 genes, which build a gene cluster together with a third desaturase gene, FADS3 [15]. FADS2 is the rate-limiting step on the metabolic pathway leading to AA and DHA production [16]. Genes involved in the elongation of LC-PUFAs include ELOVL4 (specifically expressed in human retina), and ELOVL2 and ELOVL5 that are expressed in brain [17]. Polymorphisms of the FADS gene cluster and ELOVL2 gene have been reported to be associated with plasma variation in LC-PUFA levels [18]–[22]. In addition, in the absence of dietary supplies, desaturation and elongation processes may play a role in ensuring supply of LC-PUFAs such as DHA and AA [23].
Research on the effects of LC-PUFA supply during pregnancy and lactation taking into account maternal and child inter-individual genetic variation in fatty acid desaturase and elongase activities may be critical to elucidate whether breastfeeding influences cognitive development as a unique source of LC-PUFAs, independent of other risk factors [24]. Unlike breastfeeding practises, genotype is not likely to be influenced by social, behavioural or lifestyle variables. Therefore, using the principles of Mendelian randomization, associations between the functional polymorphism of FADS and ELOVL genes and child cognition are unlikely to be confounded by social or behavioral factors and are more likely to be representative of a causal relationship. The aims of this study were: 1) to determine whether maternal genetic variants in the FADS cluster and ELOVL genes contribute to differences in LC-PUFA levels in first breast milk (colostrum); 2) to analyze whether these maternal variants are related to cognitive outcomes in the child; and 3) to assess whether children's variants in these genes modify breastfeeding effects on cognition. To achieve these objectives we used data from two Spanish population-based birth cohorts included in the INMA (INfancia y Medio Ambiente [Environment and Childhood]) Project.
Methods
Ethics statement
Written informed consent was obtained from all participants and the study was approved by the Clinical Research Ethical Committee of the Municipal Institute of Health Care (CEIC-IMAS), Barcelona.
Study design and participants are summarized in Figure 1. Distribution of child and maternal characteristics of participants and non participants are showed in the Online Repository (see Supporting Information Table S5).
The INMA-Sabadell birth cohort
The INMA-Sabadell cohort was established in Sabadell (Barcelona, Spain) between 2004 and 2006 [25]. A total of 657 women were enrolled at the first trimester of pregnancy (participation rate was 60%, 657 of 1097 eligible women), and 622 (94%) followed until the child's birth. Detailed information about child feeding was completed by interviewer-administered questionnaires with mothers at 6 and 14 months. Exclusively breastfeeding was defined as receiving breast milk only but allowing supplementation of non-milk liquids according to the WHO's definition [26]. Mental development was assessed at age 14 months (range 12–17 months) by 2 specially trained psychologists using the Bayley Scales of Infant Development, first edition [27], [28]. The Bayley Scale for Cognitive Development is composed by two scales, the mental development and psychomotor. For the present study, the main outcome was the score in the mental development scale consisting in 163 items that assess age-appropriate cognitive development in areas such as performance ability, memory, and first verbal learning.
The AMICS INMA-Menorca birth cohort
The Menorca Asthma Multicenter Infant cohort Study (AMICS) is a population-based birth cohort established in the island of Menorca (Spain) in 1997–98. A total of 482 children were recruited at birth from 492 pregnant women residing on the island [29]. Breastfeeding was assessed by interviewer-administered questionnaires to mothers three times at ages 6 months, 14 months, and 2 years[30]. Infants not breastfeed received formula feeding at a time before LC-PUFA supplementation of formula became widely available in Spain. Two certified psychologists performed the neuropsychological testing of the children at age 4 years measuring cognitive functioning with the Spanish version of the McCarthy Scales of Children's Abilities (MCSA), consisting of 5 sub-scales assessing several cognitive domains (verbal, perceptual-performance, quantitative, memory, and motor) and a general cognitive index [31]–[33]. For the present study score in the general cognitive index was used as the main outcome.
Genotyping procedure
All genotyping processes were performed at the Barcelona Node of the Spanish National Genotyping Centre (CEGEN). In the INMA-Sabadell cohort maternal DNA samples were extracted from blood obtained at first trimester of pregnancy. Child DNA samples were extracted from cord blood obtained at birth in the INMA-Sabadell cohort (100%) and from blood (85%) obtained at age 4 in the INMA-Menorca cohort (15 % from saliva). Details of the SNP selection process are provided in the supplementary information material (see Supporting information Appendix S1). Rs numbers, gene location, minor allele frequencies and Hardy-Weinberg Equilibrium values are shown in the Online Repository (see Supporting Information Table S1). Genotyping was performed blind to measures of breastfeeding and LC-PUFA exposure, as well as, cognition scores. Finally, twenty-six SNPs out of 36 initially selected were successful genotyped and included in the analysis (13 in the FADS gene cluster, 6 in ELOVL2 gene, and 7 in ELOVL5 gene). Linkage disequilibrium (LD) parameters (D' and r2) were estimated using the Haploview software (see Supporting Information Table S6, Table S7, Table S8).
Measurement of colostrum LC-PUFA
Colostrum samples were obtained from the first 270 willing participants from the INMA-Sabadell cohort. Samples were collected during the first 48–96 hours post-partum, transported to the laboratory in ice boxes less than 2 hours after collection and stored at −80°C until analyzed. The following fatty acids were identified and quantified according to the method developed by Molto-Puigmarti et al. [34]: linoleic acid (LA), α-linolenic acid (ALA); AA; eicosatrienoic acid (DGLA); EPA; docosapentaenoic acid (DPA) and DHA. The ratio between the concentrations of the product-substrate pairs of the desaturase and elongase reactions have been previously used as strong indicators for the efficiency of the enzymes and as a proxy for enzymatic activity [35]. Thus, we further calculated the following ratios as surrogate indexes of FADS and ELOVL5 enzymatic activities: AA:DGLA for FADS1, DGLA:LA and DHA:DPA for FADS2 as previously described [36], [37], and DPA:EPA for ELOVL5 activity. In addition, ratios EPA/AA and DHA/AA, which reflect the proportions of key n-3 to n-6 fatty acids essential for brain growth and function, were calculated.
Statistical analysis
The outcomes of interest, general cognitive index from McCarthy test and the mental development scale score from the Bayley test, were standardized to a mean of 100 corresponding to the mean of the raw scores, and standard deviation of 15. LC-PUFAs were log-transformed to achieve a normal distribution. Bivariate analyses exploring how breastfeeding behaviors and colustrum LC-PUFA profiles were related to either child cognition or to genetic variants used χ2-tests for categorical data and Krukall-Wallis tests for continuous data. For genetic variables, departure from HWE within mothers and children was evaluated using a χ2-test.
Associations between maternal genetic variants and both child cognition and LC-PUFA levels and ratios in colostrum were assessed using multivariable linear models and evaluated using likelihood ratio tests. The resultant regression coefficients represent the mean difference in cognitive scores or LC-PUFA levels between genotype groups. All variables significantly associated with general cognitive scores in the bivariate models (p<0.2) were included in the multivariate model, and only retained if they had an at least marginally significant association (p<0.1) or modified the coefficient by at least 5% [32], [33], [38]. Final multivariate models were adjusted for sex, child age (days), psychologist, quality of neuropsychological test, maternal education, breastfeeding, and use of gas stove at home. For associations between LC-PUFA and maternal genotype, unadjusted results are reported since incorporating potential confounders in multivariable linear models including maternal age, body mass index, smoking habits, and maternal diet intake of fish, seafood, meat and egg intakes measured by a 101-item validated food frequency questionnaire during pregnancy did not meet criteria for confounding (change-in-estimate ≥5%), and finally were excluded for parsimony. Bonferroni correction was applied to take multiple testing into account and p value thresholds were set at 0.002 (0.05/23).
We tested interactions between child's genetic variants and breastfeeding, as well as, between LC-PUFA levels and n-3/n-6 ratios in colostrum. Linear models were also used to assess whether children's genetic variants modified the effects of breastfeeding behaviors on cognitive test scores, adjusting for other covariables. Similarly, linear models were used to assess whether children's genetic variants modified the effects of colostrum LC-PUFA profiles on cognitive scores. For these analyses, selected single LC-PUFA and n-3/n-6 ratios were dichotomized as close as possible to p = 0.50, and analyses were restricted to breastfed children. The resultant regression coefficient gives the mean difference in cognitive scores between the two groups (e.g. children exposed to low vs. high colostrum fatty acid concentrations or children with contrasting genetic variants), having controlled for confounders. Statistical analysis was performed with R statistical package version 2.9.1 using the SNPassoc package [39].
Results
Associations between maternal FADS gene cluster and ELOVL gene family SNPs and LC-PUFA profiles in colostrum
Associations between maternal genetic variants and LC-PUFA levels in colostrum are showed in the Online Repository (see Supporting Information Table S2). The rs174537-T, rs174570-T, rs2072114-G, rs174602-G, rs526126-G, rs174626-C, rs174464-T, and rs174468-G alleles in the FADS cluster were related to lower levels of AA (Table S2). The rs174602-G, and 174464-T in the FADS cluster were associated with lower levels of DHA, and other SNPs in the cluster showed a trend to be associated with DHA levels in the same direction as the observed for AA (Table S2). Regarding elongases rs953413-A and rs3798719-T in ELOVL2 gene were associated with higher EPA levels. None of the SNPs in the ELOVL5 were associated with AA, EPA or DHA levels at a p value of 0.05, although a trend was observed for rs12207094 and EPA levels. Finally, several SNPs in the FADS cluster (rs174464 and rs174468), ELVOL2 (rs3734397, rs953413 and rs3798719) and ELVOL5 (rs17544159, rs9395855 and rs12207094) genes were associated with EPA/AA ratio and/or DHA/AA ratio (Table S2).
Table 1 shows the associations observed between the FADS cluster and ELOVL gene family SNPs and enzymatic indexes estimated based on LC-PUFA profiles. The rs174537-T, rs968567-A, rs2072114-G, rs526126-G, rs174626-C, rs174627-T, rs174464-T, rs174468-G alleles were associated with lower FADS1 index. Three of the previous SNPs (rs174537, rs968567 and rs174627) were inversely associated with FADS2 index. Moreover, the rs174602-G allele was associated with lower FADS2 index (capacity to synthesize DHA from EPA). In ELOVL5, the rs2397142-C allele was nominally associated with lower capacity to metabolize EPA to DPA.
Associations between maternal polymorphisms in the FADS gene cluster and ELOVL genes and child cognition
In the FADS gene cluster, we found that maternal SNPs rs968567, rs174602, rs174627 and rs174464 were nominally associated with child cognition (Table 2). In general the allele associated with higher child cognition scores were the alleles associated with high FADS2 index (the limiting enzyme of the pathway), which correspond to low FADS1 index (as both enzymatic ratios are estimated using DGLA levels). The most significant association found was for the maternal SNP rs174627 (p<0.001). In particular children of mothers carrying at least one copy of the minor allele (T) for this variant had higher cognition scores compared to children of mothers homozygotes for the major allele (Table 2). None of the SNPs in the ELOVL2 gene showed a significant association with child cognition. Regarding the ELOVL5 gene, we found that children of mothers carrying the rs17544159-C allele (associated with higher EPA/AA ratio in colostrums) or carrying the rs12207094-T allele (associated with both higher EPA/AA and DHA/AA ratios in colostrum) had higher cognition scores related to children of mothers homozygote for the major allele.
Associations between child polymorphisms in the FADS gene cluster and ELOVL genes and child cognition
Regarding child genetic variants, after Bonferroni correction no statistically significant associations were found between child polymorphisms and cognition, and we could no replicate the nominal associations found in Menorca cohort (see Supporting Information Table S3).
Child polymorphisms in the FADS cluster and ELOVL5 gene modify breastfeeding effect on cognition
We found two statistically significant gene-breastfeeding interactions that were replicated in both cohorts. We found that not being breastfed conferred a disadvantage in cognition (9 point in the INMA-Menorca cohort and 8 points in the INMA-Sabadell cohort) among children GG homozygotes for rs174468 (low FADS1 index), but not among those carrying at least one copy of the A allele (high FADS1 index) (Figure 2). In addition, not being breastfed resulted in a disadvantage in cognition (8-point in INMA-Menorca and 5-point in the INMA-Sabadell) among children CC homozygotes for rs2397142 (low ELOVL5 index), but not among those carrying at least one copy of the G allele (high ELOVL5 index). In contrast, breastfed children did not differ in cognition score irrespective of their genetic variants in these polymorphisms.
Bars represent standard error. 1Adjusted for sex, school trimester at testing, maternal social class, maternal education, maternal smoking and alcohol consumption in pregnancy, indoor nitrogen dioxide (NO2) levels. 2Adjusted for sex, child age (in days), psychologist, quality of neuropsychological test, maternal education, and use of gas stove at home.
Children' genotypes were not related to potential confounders including maternal social class, maternal education level, reproductive outcomes, diet intakes, or maternal tobacco and alcohol consumption during pregnancy. In addition, the prevalence of breastfeeding was similar across children with different genotypes (see Supporting Information Table S4).
Child polymorphism in ELOVL5 gene modifies the effect of colostrum LC-PUFA levels on cognition
In the INMA-Sabadell cohort we analyzed whether child genetic variants interacted with LC-PUFA levels in colostrum to influence cognition. Among breastfed children carrying at least one G allele for rs2397142 (high ELOVL5 index), there was no association between cognition and LC-PUFA levels in colostrum. In contrast, among children CC homozygotes for rs2397142 (low ELOVL5 index), exposure to high levels of EPA, and high EPA/AA and DHA/AA ratios were associated with a 6.7-, 8.6- and 10.1-points advantage in scores, respectively, relative to children exposed to low levels (Table 3).
Discussion
In the present study, we analysed polymorphisms in genes encoding the key enzymes involved in LC-PUFA synthesis in mothers and their offspring to disentangle their role in modifying potential nutritional advantages of breastfeeding on child cognition. Maternal genetic variants in the FADS gene cluster and ELOVL5 gene were associated with higher colostrum levels of n-3 LC-PUFA (i.e. EPA and DHA) as well as with higher cognitive scores in their children. In infants, we found that genetic variants in the FADS cluster and ELOVL5 gene modified the effects of breastfeeding on cognition. Children with variants associated with lower synthesis of LC-PUFA had higher scores when breastfed, while those with greater capacity to synthesize these fatty acids had higher scores regardless of breastfeeding practises. We confirmed these gene-breastfeeding interactions in two independent population-based birth cohorts. Moreover, we ruled out the possibility that these findings were due to potential confounders including correlations between genes and exposures, intrauterine growth differences, dietary or sociodemographic factors.
Maternal genetic variation in LC-PUFA metabolizing enzymes, colostrum LC-PUFA profile and child cognition
Although fatty acid desaturase activities are present in the mammary gland [40], few studies have assessed whether maternal polymorphisms in key enzymes involved in LC-PUFA synthesis influence the composition of breast milk [20], [22], and the effects on child phenotypes such as cognition have not been investigated so far. To this purpose we used LC-PUFA measurements from early breast milk (colostrum), since its LC-PUFA composition is related to the LC-PUFA composition of brain tissue in the newborn [41]. Our results on the direction of the association between maternal genetic variants in desaturases and LC-PUFA levels in colostrum are consistent with previous findings showing associations between several SNPs in the FADS cluster and plasma, erythrocyte and breast milk LC-PUFA levels in adults [18]–[20], [22], [40]. In our study, the two maternal FADS cluster SNPs most strongly related to disparities in colostrum LC-PUFA levels were rs174537 previously reported to be associated with LC-PUFA levels in genome-wide association studies [18], and rs968567 recently reported to influence FADS2 transcription in an in vitro system [42]. We also identified other FADS cluster SNPs significantly related to variation in colostrum LC-PUFA levels as well as with FADS enzymatic indexes previously described [20], [37]. In general, minor alleles in the FADS cluster (associated with low FADS1 index and high FADS2 indexes) were related to higher child cognitive scores, although no statistically significant associations were found for all of them. FADS2 is considered the limiting enzyme which strengthens the importance of the present results. Since not all the LC-PUFAs were measured, we could not evaluate the same ratios in the n-3 pathway. However, the fact that both types of LC-PUFA compete for the same enzymes system, suggests that similar results are expected in the n-3 pathway. In addition, although tissue LC-PUFA levels are also determined by dietary intake, associations between the analyzed SNPs and both LC-PUFA colostrum levels and desaturase and elongase enzymatic activity estimates did not change after taking into account information on maternal diet habits including fish, seafood, meat and eggs intakes measured by food frequency questionnaire during pregnancy.
AA and DHA have been proposed to play important roles in synaptic transmission and plasticity during early brain development [8]. Furthermore, the imbalance between n-6 and n-3 fatty acids has been suggested to be related to neuropsychiatric disorders in children [8]. Although, infants are capable of synthesizing DHA and AA accumulation of LC-PUFA fatty acids in utero are derived predominantly through placental transfer [43]. Based on the present results, maternal genetic variants influencing fatty acid desaturation and elongation activities– and thus levels of LC-PUFA in colostrum– appear to have functional importance for infant brain development. In this regard, we found that children of mothers with genetically determined lower levels of AA, higher EPA/AA, and higher capacity to synthesize DHA and therefore higher DHA/AA ratios, showed a significant advantage in cognition at age 14 months.
Child polymorphisms in LC-PUFA metabolizing enzymes, breastfeeding and child cognition
After birth, formula fed infants depend on the utilization of body stores or on the endogenous ability to synthesize LC-PUFAs: the latter may help that may compensate for deficient supplies of LC-PUFAs associated with not being breastfed. We found that genetic variants in the FADS gene cluster and ELOVL5 gene in children modified the effect of breastfeeding on cognition; and results were replicated in both cohorts. Not being breastfed conferred a disadvantage in cognition among children with lower FADS1 index, and similarly, not being breastfed conferred a disadvantage among children with lower ELOVL5 index. These results of the FADS cluster variant are in general accordance with those reported by Caspi et al. [44]. As proposed by these authors the genetic modification of breastfeeding effects on cognition is unlikely to be directly caused by the analyzed SNPs, but rather mediated by more efficient LC-PUFA processing.
Biological effects of LC-PUFAs on brain function are assumed to be mediated by the availability of LC-PUFA with >20 carbon atoms and >3 double bonds (i.e. AA, EPA and DHA) [8]. In mammals, ELOVL5 catalyzes the two carbon elongation process required for the synthesis of DHA from EPA [45]. The modification of breastfeeding effects on child cognition by a SNP in ELOVL5 found in the present study could be explained by higher capacity in infants with this variant to synthesize 20 carbon fatty acids from their precursors provided by breast milk. While previous studies have focused primarily on the FADS gene cluster, interestingly we found that breastfed children who were CC homozygotes (low ELOVL5 index) exposed to higher EPA, EPA/AA and DHA/AA in colostrum had a significant advantage in cognition in comparison to those exposed to lower concentrations. These results suggest that homozygosity for the C allele could be related to less efficient synthesis of DHA from EPA, probably due to decreased transcription or to a less enzyme activity. On the contrary, the presence of at least one G allele could be related to more efficient fatty acid metabolism, which could compensate for deficient supplies of LC-PUFAs derived from not being breastfed or being exposed to breast milk containing lower levels of EPA and DHA. Consistent with this mechanism, experimental studies in baboon neonates have revealed that exposure to different levels of DHA in formulas induced global changes in gene expression in cerebral cortex [46], and only a single elongation enzyme was differentially expressed: the analogue to human ELOVL5.
Implications for genetic neuroscience and newborn nutrition
Identifying genes for variation in the range of normal cognitive functioning could help in understanding biological mechanisms and pathways of milder but more-prevalent forms of impaired neurofunctioning, which are often associated with autism, schizophrenia, and attention deficit hyperactivity disorder (ADHD). Complex traits, such as cognition, are the result of combined effects of genetic and environmental factors[47]; however, gene-environment interactions are repeatedly overlooked in genome-wide scans which could result in unexplained heritability [48]. In addition, the results indicate that breastfeeding offers a real nutritional advantage for child cognition, and that this advantage is related to LC-PUFA profiles. The ratios between n-3 and n-6 fatty acids, particularly EPA/AA and DHA/AA, have been emphasized in infant nutrition [8]. Significant improvement in ADHD symptoms have been observed with a significant increase in EPA and DHA supplies, as well as a drop in the AA/EPA ratio [49]. Also a significant association between a SNP in the FADS2 gene and ADHD has been reported (carriers of the genetic variant associated with ADHD have greatly elevated AA contents in their serum phosplipids) [50]. These earlier studies and findings from the present study suggest not only a need to increase n-3 (EPA and DHA) in the diet, including in infant formulae, but also perhaps to decrease the n-6 content (i.e. AA). Furthermore, failure to consider genetic heterogeneity in LC-PUFA metabolism in participants may dilute supplementation effects and as a result inconsistent conclusions are derived from randomized controlled trials.
Strengths and limitations
Strengths of this study include the use of data from two independent population-based birth cohorts. We also had both maternal and child genetic information. Cognitive assessment was measured with validated instruments. The observed associations are unlikely to be due to population stratification since the discovery and replication samples were restricted to Caucasian individuals of European origin recruited in Catalonia and in Menorca. Furthermore, the frequency distribution of minor and major alleles was similar in both cohorts and to that listed in the HapMap CEU sample. In comparison with the study by Caspi et al., we assessed the association of both maternal and child genetic variation in fatty acid desaturases and elongases on cognition. Moreover we assessed the genetic variability in the whole FADS cluster gene and also in elongase genes involved in the LC-PUFA synthesis pathway. In addition, the use of biomarkers such as colostrum levels of LC-PUFA strengthens the findings.
Among the study limitations we acknowledge the small sample size that has probably limited power to detect associations with small effect sizes and associations with rare variants. Since we measured cognitive function at two different ages, phenotypic heterogeneity between cohorts may have limited our ability to detect some direct associations. However, our main aim was to seek for gene-environment interactions, which probably have a more relevant role in complex traits as cognition. In fact we were able to detect and replicate two gene-breastfeeding interactions with a plausible biological justification. Additionally, a proportion of mothers and children in both cohorts were not included in the present study because of non-participation in the neurodevelopment exam or failure in the genotyping. However, the analysis comparing participants and non-participants did not show any difference in general cognition or sociodemographic characteristics (see Supporting Information Table S5). We could no collect maternal DNA in Menorca cohort, thus associations between maternal genetic variants and child cognition could not be replicated and future studies are warranted.
Conclusions
Results of the present study showed that LC-PUFA supplies during pregnancy and lactation, genetically determined by maternal desaturase and elongase activities, appear to have functional importance to the infant brain development. In addition, breastfeeding effects on cognition are also modified by child genetic variants in desaturase and elongase enzymes involved in the control of LC-PUFA pathways. Future research assessing the effects of breastfeeding and LC-PUFA diet supplementation on brain development should take into account genetic heterogeneity in key enzymes of LC-PUFA synthesis.
Supporting Information
Figure S1.
The n-6 and n-3 fatty acid metabolism pathways.
https://doi.org/10.1371/journal.pone.0017181.s002
(DOC)
Table S1.
Characteristics of the SNPs analyzed in the FADS gene cluster, ELOVL2 and ELOVL5 genes. Abbreviations: SNP: single-nucleotide polymorphism; HWE: Hardy-Weinberg equilibrium; MAF: Minor allele frequency; ME: Mendelian errors.
https://doi.org/10.1371/journal.pone.0017181.s003
(DOC)
Table S2.
Associations* between maternal genetic variants in the FADS cluster, ELOLV2 and ELOVL5 genes and levels of LC-PUFA in colostrum. *p-value of the association using additive genetic models. Change in concentrations in the additive genetic models assuming a trend per copy of the minor allele. MAF: Minor allele frequency; AA: Arachidonic acid; EPA: Eicosapentanoic acid; DHA: Docosahexaenoic acid.
https://doi.org/10.1371/journal.pone.0017181.s004
(DOC)
Table S3.
Associations* between child polymorphisms in the FADS cluster, ELOVL2 and ELOVL5 genes and child cognition by cohort. *p values in the additive genetic models assuming a trend per copy of the minor allele. MAF: Minor allele frequency.
https://doi.org/10.1371/journal.pone.0017181.s005
(DOC)
Table S4.
Association of child genetic polymorphisms in the FADS cluster (rs174468) and ELOVL5 (rs2397142) with general cognitive function score and selected covariates by cohort. ¶ Unless otherwise specified, p value derived from chi-2 test; §p value derived from Kruskall-Wallis test.
https://doi.org/10.1371/journal.pone.0017181.s006
(DOC)
Table S5.
Comparison of the distribution of child and maternal characteristics between children included and not included in the study by cohort. ¶ Unless otherwise specified, p value derived from chi-2 test; §p value derived from Kruskall-Wallis test.
https://doi.org/10.1371/journal.pone.0017181.s007
(DOC)
Table S6.
Pair-wise LD parameters (D'r2) in the FADS gene cluster in the INMA Menorca cohort (only SNPs included in the final analysis).
https://doi.org/10.1371/journal.pone.0017181.s008
(DOC)
Table S7.
Pair-wise LD parameters (D'/r2) in the ELOVL2 gene in the INMA Menorca cohort (only SNPs included in the final analysis).
https://doi.org/10.1371/journal.pone.0017181.s009
(DOC)
Table S8.
Pair-wise LD parameters (D'/r2) in the ELVOL5 gene in the INMA Menorca cohort (only SNPs included in the final analysis).
https://doi.org/10.1371/journal.pone.0017181.s010
(DOC)
Acknowledgments
The authors are grateful to Silvia Fochs, Anna Sànchez, Maribel López, Nuria Pey, Muriel Ferrer, and Cristina Capó for their assistance in contacting the families and administering the questionnaires.
The authors would particularly like to thank all the participants for their generous collaboration. A full roster of the INMA Project Investigators can be found at http://www.proyectoinma.org/presentacion-inma/listado-investigadores/en_listado-investigadores.html.
Author Contributions
Conceived and designed the experiments: EM MB JRG MT XE JS. Performed the experiments: MB XE. Analyzed the data: EM JRG. Contributed reagents/materials/analysis tools: MB CMP CLS XE. Wrote the manuscript: EM MB. Critical revision of the manuscript for important intellectual content: MG MT RGE JJ JF MM MV CMP JS.
Referentes
- 1. Haag M (2003) Essential fatty acids and the brain. Can J Psychiatry 48: 195–203.
- 2. Uauy R, Hoffman DR, Peirano P, Birch DG, Birch EE (2001) Essential fatty acids in visual and brain development. Lipids 36: 885–895.
- 3. Innis SM (2000) The role of dietary n-6 and n-3 fatty acids in the developing brain. Dev Neurosci 22: 474–480.
- 4. He C, Qu X, Cui L, Wang J, Kang JX (2009) Improved spatial learning performance of fat-1 mice is associated with enhanced neurogenesis and neuritogenesis by docosahexaenoic acid. Proc Natl Acad Sci U S A 106: 11370–11375.
- 5. Robson LG, Dyall S, Sidloff D, Michael-Titus AT (2010) Omega-3 polyunsaturated fatty acids increase the neurite outgrowth of rat sensory neurones throughout development and in aged animals. Neurobiol Aging 31: 678–687.
- 6. Salvati S, Natali F, Attorri L, Di Benedetto R, Leonardi F, et al. (2008) Eicosapentaenoic acid stimulates the expression of myelin proteins in rat brain. J Neurosci Res 86: 776–784.
- 7. Marza E, Lesa GM (2006) Polyunsaturated fatty acids and neurotransmission in Caenorhabditis elegans. Biochem Soc Trans 34: 77–80.
- 8. Schuchardt JP, Huss M, Stauss-Grabo M, Hahn A (2010) Significance of long-chain polyunsaturated fatty acids (PUFAs) for the development and behaviour of children. Eur J Pediatr 169: 149–164.
- 9. Sauerwald TU, Demmelmair H, Koletzko B (2001) Polyunsaturated fatty acid supply with human milk. Lipids 36: 991–996.
- 10. Makrides M, Neumann MA, Byard RW, Simmer K, Gibson RA (1994) Fatty acid composition of brain, retina, and erythrocytes in breast- and formula-fed infants. Am J Clin Nutr 60: 189–194.
- 11. Anderson JW, Johnstone BM, Remley DT (1999) Breast-feeding and cognitive development: a meta-analysis. Am J Clin Nutr 70: 525–535.
- 12. Mortensen EL, Michaelsen KF, Sanders SA, Reinisch JM (2002) The association between duration of breastfeeding and adult intelligence. JAMA 287: 2365–2371.
- 13. Der G, Batty GD, Deary IJ (2006) Effect of breast feeding on intelligence in children: prospective study, sibling pairs analysis, and meta-analysis. BMJ 333: 945.
- 14.
Simmer K, Patole SK, Rao SC (2008) Longchain polyunsaturated fatty acid supplementation in infants born at term. Cochrane Database Syst Rev. CD000376 p.
- 15. Marquardt A, Stöhr H, White K, Weber BH (2000) cDNA cloning, genomic structure, and chromosomal localization of three members of the human fatty acid desaturase family. Genomics 66: 175–183.
- 16. Sprecher H (2000) Metabolism of highly unsaturated n-3 and n-6 fatty acids. Biochimica Et Biophysica Acta-Molecular and Cell Biology of Lipids 1486: 219–231.
- 17. Jakobsson A, Westerberg R, Jacobsson A (2006) Fatty acid elongases in mammals: their regulation and roles in metabolism. Prog Lipid Res 45: 237–249.
- 18. Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, et al. (2009) Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet 5: e1000338.
- 19. Schaeffer L, Gohlke H, Müller M, Heid IM, Palmer LJ, et al. (2006) Common genetic variants of the FADS1 FADS2 gene cluster and their reconstructed haplotypes are associated with the fatty acid composition in phospholipids. Hum Mol Genet 15: 1745–1756.
- 20. Xie L, Innis SM (2008) Genetic variants of the FADS1 FADS2 gene cluster are associated with altered (n-6) and (n-3) essential fatty acids in plasma and erythrocyte phospholipids in women during pregnancy and in breast milk during lactation. J Nutr 138: 2222–2228.
- 21. Martinelli N, Girelli D, Malerba G, Guarini P, Illig T, et al. (2008) FADS genotypes and desaturase activity estimated by the ratio of arachidonic acid to linoleic acid are associated with inflammation and coronary artery disease. Am J Clin Nutr 88: 941–949.
- 22. Moltó-Puigmartí C, Plat J, Mensink RP, Müller A, Jansen E, et al. (2010) FADS1 FADS2 gene variants modify the association between fish intake and the docosahexaenoic acid proportions in human milk. Am J Clin Nutr 91: 1368–1376.
- 23. Heird WC, Lapillonne A (2005) The role of essential fatty acids in development. Annu Rev Nutr 25: 549–571.
- 24. Smith GD, Ebrahim S (2004) Mendelian randomization: prospects, potentials, and limitations. Int J Epidemiol 33: 30–42.
- 25. Ribas-Fitó N, Ramón R, Ballester F, Grimalt J, Marco A, et al. (2006) Child health and the environment: the INMA Spanish Study. Paediatr Perinat Epidemiol 20: 403–410.
- 26.
World Health Organization (2008) Indicators for assessing infant and young child feeding practices: Part 1 Definitions. Geneva, Switzerland.
- 27.
Bayley N (1977) Escalas Bayley de Desarrollo Infantil. Madrid: TEA Ediciones.
- 28. Vrijheid M, Martinez D, Forns J, Guxens M, Julvez J, et al. (2010) Prenatal exposure to cell phone use and neurodevelopment at 14 months. Epidemiology 21: 259–262.
- 29. Polk S, Sunyer J, Muñoz-Ortiz L, Barnes M, Torrent M, et al. (2004) A prospective study of Fel d1 and Der p1 exposure in infancy and childhood wheezing. Am J Respir Crit Care Med 170: 273–278.
- 30. Ribas-Fitó N, Júlvez J, Torrent M, Grimalt JO, Sunyer J (2007) Beneficial effects of breastfeeding on cognition regardless of DDT concentrations at birth. Am J Epidemiol 166: 1198–1202.
- 31.
McCarthy D (1972) Manual for the McCarthy Scales of Children's Abilities. 1996). New York: Psychological Corporation. (Spanish adaptation: Madrid: TEA Ediciones;.
- 32. Julvez J, Ribas-Fitó N, Torrent M, Forns M, Garcia-Esteban R, et al. (2007) Maternal smoking habits and cognitive development of children at age 4 years in a population-based birth cohort. Int J Epidemiol 36: 825–832.
- 33. Morales E, Julvez J, Torrent M, de Cid R, Guxens M, et al. (2009) Association of early-life exposure to household gas appliances and indoor nitrogen dioxide with cognition and attention behavior in preschoolers. Am J Epidemiol 169: 1327–1336.
- 34. Moltó-Puigmartí C, Castellote AI, López-Sabater MC (2007) Conjugated linoleic acid determination in human milk by fast-gas chromatography. Anal Chim Acta 602: 122–130.
- 35. Gieger C, Geistlinger L, Altmaier E, Hrabé de Angelis M, Kronenberg F, et al. (2008) Genetics meets metabolomics: a genome-wide association study of metabolite profiles in human serum. PLoS Genet 4: e1000282.
- 36. Warensjö E, Rosell M, Hellenius ML, Vessby B, De Faire U, et al. (2009) Associations between estimated fatty acid desaturase activities in serum lipids and adipose tissue in humans: links to obesity and insulin resistance. Lipids Health Dis 8:
- 37. Bokor S, Dumont J, Spinneker A, Gonzalez-Gross M, Nova E, et al. (2010) Single nucleotide polymorphisms in the FADS gene cluster are associated with delta-5 and delta-6 desaturase activities estimated by serum fatty acid ratios. J Lipid Res 51: 2325–2333.
- 38. Ribas-Fitó N, Cardo E, Sala M, Eulàlia de Muga M, Mazón C, et al. (2003) Breastfeeding, exposure to organochlorine compounds, and neurodevelopment in infants. Pediatrics 111: e580–585.
- 39. González JR, Armengol L, Solé X, Guinó E, Mercader JM, et al. (2007) SNPassoc: an R package to perform whole genome association studies. Bioinformatics 23: 654–655.
- 40. Rodriguez-Cruz M, Tovar AR, Palacios-González B, Del Prado M, Torres N (2006) Synthesis of long-chain polyunsaturated fatty acids in lactating mammary gland: role of Delta5 and Delta6 desaturases, SREBP-1, PPARalpha, and PGC-1. J Lipid Res 47: 553–560.
- 41. Xiang M, Alfven G, Blennow M, Trygg M, Zetterstrom R (2000) Longchain polyunsaturated fatty acids in human milk and brain growth during early infancy. Acta Paediatr 89: 142–147.
- 42. Lattka E, Eggers S, Moeller G, Heim K, Weber M, et al. (2010) A common FADS2 promoter polymorphism increases promoter activity and facilitates binding of transcription factor ELK1. J Lipid Res 51: 182–191.
- 43. Koletzko B, Lien E, Agostoni C, Böhles H, Campoy C, et al. (2008) The roles of long-chain polyunsaturated fatty acids in pregnancy, lactation and infancy: review of current knowledge and consensus recommendations. J Perinat Med 36: 5–14.
- 44. Caspi A, Williams B, Kim-Cohen J, Craig IW, Milne BJ, et al. (2007) Moderation of breastfeeding effects on the IQ by genetic variation in fatty acid metabolism. Proc Natl Acad Sci U S A 104: 18860–18865.
- 45. Leonard AE, Kelder B, Bobik EG, Chuang LT, Lewis CJ, et al. (2002) Identification and expression of mammalian long-chain PUFA elongation enzymes. Lipids 37: 733–740.
- 46. Kothapalli KS, Anthony JC, Pan BS, Hsieh AT, Nathanielsz PW, et al. (2007) Differential cerebral cortex transcriptomes of baboon neonates consuming moderate and high docosahexaenoic acid formulas. PLoS One 2: e370.
- 47. Posthuma D, Luciano M, Geus EJ, Wright MJ, Slagboom PE, et al. (2005) A genomewide scan for intelligence identifies quantitative trait loci on 2q and 6p. Am J Hum Genet 77: 318–326.
- 48. Thomas D (2010) Gene-environment-wide association studies: emerging approaches. Nat Rev Genet 11: 259–272.
- 49. Sorgi PJ, Hallowell EM, Hutchins HL, Sears B (2007) Effects of an open-label pilot study with high-dose EPA/DHA concentrates on plasma phospholipids and behavior in children with attention deficit hyperactivity disorder. Nutr J 6: 16.
- 50. Brookes KJ, Chen W, Xu X, Taylor E, Asherson P (2006) Association of fatty acid desaturase genes with attention-deficit/hyperactivity disorder. Biol Psychiatry 60: 1053–1061.