Introduction

The process of the human brain developing is a protracted process that starts in the third week of gestation with the differentiation of the neural progenitor cells and lasts at least until late adolescence, and possibly throughout life [1]. Early gene expression changes and environmental factors play a crucial role in the normal development of the brain [2]. In particular, placenta is posited as a critical determinant of immediate and long-term neurodevelopmental outcomes in children [3]. The placenta–brain axis refers to this close relationship between the placenta and the brain [4], and, as recently reviewed by Kratimenos et al. [5], compromised placental function can predispose an individual to psychiatric disorders. During pregnancy, placenta is the organ that transfers oxygen and nutrients from the mother’s blood to the fetal brain. Throughout the whole pregnancy, placenta plays a crucial role as it protects the fetus from environmental harm and supports the development of its brain. Additionally, the placenta aids in immune system control and functions as a neuroendocrine organ, generating hormones that act as growth factors and neuropeptides that operate as neurotransmitters [6].

Genomic and epigenomic changes in the placenta have been linked to neurodevelopmental abnormalities and typical neurobehavioral development in epidemiological and animal studies [7]. In particular, the epigenome can be considered a proxy of placental function during pregnancy and a plausible mediator of the association between intrauterine environmental exposures and genetics, and neurodevelopment and brain health trajectories [8]. DNA methylation (DNAm) on cytosine nucleotides followed by guanine (CpG) across the genome is the most studied epigenetic mechanism, having an important role in regulating gene expression [9]. For example, some authors have suggested that differences in placental methylation patterns are associated with differences in brain development, potentially through mechanisms such as the reprogramming of gene expression involved in the hypothalamic-pituitary-adrenal axis [10].

Previous research has revealed that a child’s cognitive capacity is not just influenced by genetics and that there are a number of other factors associated, such as socioeconomic status [11], the mother’s gestational age [12], the mother’s intelligence [13], family context [14], and, most importantly for this study, prenatal exposure to detrimental environmental factors such as cigarette smoke, organochlorine pollutants, maternal malnutrition, obesity, or stress, among others [15,16,17,18,19]. Thus, it is thought that placental function, and in particular placenta DNAm, could be a plausible mechanism contributing to this prenatal environmental influence on cognitive functions later in life. In fact, neurodevelopmental disorders or traits such as autism spectrum disorder or neurobehavior have been investigated in relation to placenta DNAm [19,20,21].

However, to the best of our knowledge, the association between placenta DNAm and child’s cognitive functions has not been addressed before. Notably, a recent epigenome-wide association study (EWAS) in older-aged adults (N = 6809) reported associations of blood DNAm with global cognitive function at one intergenic CpG site on chromosome 12 and with phonemic verbal fluency at one CpG site on chromosome 10 in the INPP5A gene [22]. Moreover, using cord blood data, the Pregnancy and Childhood Epigenetics (PACE) consortium has recently conducted an epigenome-wide association study for cognitive functioning in childhood (N = 3300) [23]. Overall, they did not find robust evidence that cord blood DNAm at the single CpGs investigated could be associated with cognitive functions, either overall, verbal, or perceptive performance. On the other hand, with a much smaller sample size (N = 112) but different neurodevelopmental outcome, Abrishamcar et al. [24] found that methylation at several CpGs in cord blood was significantly mediating the effect of prenatal tobacco exposure and prenatal alcohol exposure on cognition and motor skills at 6 months of age in a south-african cohort. In the same cohort, but another study, the authors found that methylation at 3 CpGs in cord blood was associated with severe neurodevelopmental delay at 2 years of age [25]. In the context of adverse birth outcomes, Camerota et al. [26] assessed neonatal DNAm from buccal swabs, such as in children born very preterm, and found that methylation at 309 CpGs mediated the association between cumulative prenatal risk and child cognitive ability at 3 years of age. The same authors also explored attention problems in another study and found that DNAm at 33 CpG sites was associated with this outcome at 2 years of age [27]. In line with this, other studies have focused more on behavioural traits instead of cognition. Notable contributions, again within the PACE consortium, include work by Neumann et al. [28], who reported that DNAm at 9 CpGs in cord blood were associated with later ADHD symptoms, and by Rijlaarsdam et al. [29] who found that DNAm at 1 CpG and 1 differentially methylated region (DMR) in child blood was associated with general psychopathology factor at school-age.

Because of all the exposed above, the placenta is probably an appropriate organ for studying the early mechanisms and determinants of cognition. Thus, our goal here is to investigate whether differences in placenta DNAm profile are associated with three different cognitive domains (namely verbal score, perceptive performance score, and general cognitive score) in childhood, accounting for relevant covariates. We also assess the effects of maternal cognitive function, gestational age, and pregnancy complications and adverse birth outcomes on this association. To do this, we conducted epigenome-wide association analyses including data from mother-child pairs from 3 subcohorts within the INMA project (INfancia y Medio Ambiente, acronism of Spanish “Chilhood and Environment) and performed a follow-up functional analysis to help the interpretation of the findings.

Materials/subjects and methods

Study participants

The study population belongs to the INMA project (INfancia y Medio Ambiente-Environment and Childhood), which is a network of birth cohorts in Spain created with the aim to study the role of environmental pollutants in air, water and diet during pregnancy and early childhood in relation to child growth and development [30]. In our study, we used data from participants within the subcohorts of Gipuzkoa, Sabadell and Valencia due to data availability. Of the 2506 mother-child pairs that were followed, a sub-sample of 255 mothers-child pairs with complete data on placenta genome-wide methylation and child’s general cognitive functions at 4 years of age were selected for this study. Mother-child pairs of multiple births or mother-child pairs with a non-European ancestry were excluded. For more details, please see Supplementary Methods.

Placental biopsy, DNA extraction and processing

In INMA, 2506 mother-child pairs were followed until birth and a random selection of 489 placentas were collected. Genomic DNA from placenta was isolated using the DNAeasy® Blood and Tissue Kit, (Qiagen, CA, USA) and stored at −20 °C until further processing. DNAm was assessed with the Infinium MethylationEPIC BeadChip from Illumina, following the manufacturer’s protocol in the Erasmus Medical Centre core facility. For more details on the specific protocols followed, please see Supplementary Methods.

Methylation data acquisition, quality control and normalisation

The methylation data was pre-processed using the PACEAnalysis R package (v.0.1.7) (https.www.epicenteredresearch). The pre-processing pipeline consists of sample quality control, probe quality control, normalisation, batch correction, estimation of cell type proportions, and winsorization of outlier values (the extreme values to the 1% percentile–with 0.5% each at upper and lower ends of the distribution- were winsorized, where percentiles were estimated with the empirical beta-distribution). The final dataset consists of 379 samples (before merging with phenotypic data) and 811,990 probes. DNAm values are expressed as beta values, where 0 means un-methylation and 1 complete methylation. Cell type proportions of six populations (trophoblasts, syncytiotrophoblast, nucleated red blood cell, Hofbauer cells, and endothelial cells) were estimated from DNAm using the placenta reference panel from term placentas implemented in the planet R package [31]. Additional details regarding pre-processing can be found in the supplementary Methods.

Child’s general cognitive functions assessment

A standardised version of the McCarthy Scales of Children’s Functions (MSCA) adapted to the Spanish population by TEA Editions S.A. (official editorial company for adapting tests in Spain) was used to evaluate motor and cognitive functions [32]. In this particular study, we used three out of the five sub-area scales from the original MSCA: verbal scale, perceptive performance scale, and general cognitive scale. The three cognitive continuous scores (verbal, perceptive performance, and general cognitive scores) were standardised across the entire cohort (N = 2150) so that mean = 100 and SD = 15, to facilitate comparison among the different subcohorts and to align with conventional cognitive test scoring [30] (Fig. 1; Table 1). Participants were included if they had replied to all the tests. For more details on this assessment please see Supplementary Methods.

Fig. 1
figure 1

Histograms of the cognitive outcome variables verbal score, perceptive performance score, and general cognitive score.

Full size image
Table 1 Descriptives of the study population (N = 255).
Full size table

Statistical analysis

In this study, we conducted EWAS to test the association between the normalised beta-value at each CpG (dependent variable) and the verbal score, perceptive performance score or general cognitive scores in children (independent variable). For this, we conducted multivariable robust linear regression models to account for potential heteroscedasticity and non-normality, using the R package PACEAnalysis (https.www.epicenteredresearch). We have modelled DNAm as the dependent variable to control for extreme values (with high leverage or influential points) that could skew the association [33], particularly in placental tissue where skewed methylation patterns are common [34]. The main models included the following set of covariates selected among a list of potential a priori confounders: subcohort (Gipuzkoa, Valencia or Sabadell); child’s sex (female/male), obtained from clinical records; gestational age, calculated by ultrasounds; maternal age; maternal education (primary or without education /secondary /university); parity (primiparous vs multiparous); child’s age at MSCA assesment; maternal smoking status during pregnancy (non-smoker /non-sustained smoker – quit smoking in early pregnancy, before the end of first trimester-/sustained smoker throughout pregnancy), obtained through questionnaires administered face to face by trained interviewers in first and third trimester of pregnancy; and five cell type proportions (stromal cells, which corresponds to the lowest proportion of cell types, was removed from the models to avoid collinearity issues). Further, the models were adjusted for contamination score which was created by the PACEanalysis R package, as mentioned in the Supplementary Methods [35]. Multiple testing was corrected using the False discovery rate (FDR) through Benjamini–Hochberg (BH) procedure. CpGs with a nominal p < 1 × 10−4 were considered to show suggestive statistically significant association with the cognitive scores and they were analysed in the downstream analysis. The association coefficients were transformed so that they reflect the difference in DNAm level associated with an interquartile range (IQR) increase in cognitive scores.

Several sensitivity analyses were run consisting of: (1) additionally adjusting for maternal cognitive function score (WAIS-III Similarity subtest conducted on the participant mothers when the child was 4 years old); (2) unadjusting for gestational age; and (3) excluding infants that were born preterm (gestational age < 37 weeks), small for gestational age (neonate’s birth weight below the percentile 10 taking into account the reference table of the nearest week of gestation [36]), and/or mother-child pairs that experienced pregnancy complications (preeclampsia, gestational diabetes). A summary of the models and the corresponding sample sizes is available in Supplementary Table 1. Pearson’s correlation and GGally R package were used to compare the resulting estimates between models. The difference in the effect size among main models and alternative models was calculated as (effect size main model − effect size alternative model)/effect size main model × 100. Venn diagrams to visualise the overlap between CpGs identified in the different models were conducted using the webtool https://bioinformatics.psb.ugent.be/webtools/Venn/. We used the Quality Control module from the EASIER R package to perform the quality control of EWAS results, annotate CpGs, and estimate FDRs to account for multiple testing [37].

Differentially methylated regions

DMRs were identified using the dmrff package in R [38]. This method first identifies candidate DMRs by screening the meta-level EWAS results for genomic regions each covered by a sequence of CpG sites with EWAS effects in the same direction, EWAS p < 0.05, and <500 bp gaps between consecutive CpG sites. Then, summary statistics are calculated for each candidate DMR by meta-analysing the EWAS summary statistics of the CpG sites in the region (regions were considered DMRs if containing more than 1 CpG site). Multiple testing was corrected using the FDR.

Biological interpretation

Individual probes and probes included in the DMRs showing FDR significance were looked up in the EWAS catalogue [39] and EWAS atlas [40] to examine potential associations with exposures and health outcomes based on existing studies. To further characterise potential genetic influences on these sites, we used the placenta mQTL lists available from Delahaye et al. [41], the fetal brain mQTL list from Hannon et al. [42] (http://epigenetics.essex.ac.uk/mQTL/), the brain mQTL resource from the Brain xQTL Serve (https://mostafavilab.stat.ubc.ca/xqtl/) [43], and the whole blood mQTL list from GoDMC (http://mqtldb.godmc.org.uk/) [44]. To assess whether methylation levels of CpGs were associated with the expression levels of nearby genes, we consulted the placenta Expression Quantitative Trait Methylation (eQTM) lists available from Delahaye et al. [41], and Deyssenroth et al. [45], and the blood eQTM list from the eQTM HELIX catalogue (https://helixomics.isglobal.org/), generated from childhood blood samples.

Then, we did several functional enrichment analyses with the suggestive CpGs (p < 1 × 10−4), using the Functional Enrichment module of the EASIER R package [37]. The University of California-Santa Cruz (UCSC) genome browser was used to further explore the genomic context of the FDR-significant CpGs and DMRs. For more details on the functional enrichment analyses performed, please see Supplementary Methods.

Finally, we also explored whether genomic regions of the identified DMRs overlap with the 242, 226, and 1271 genome-wide significant loci previously reported in GWAS on intelligence [46], cognitive performance, and educational attainment [47] (0.5 Mb window centred to the genomic 250 locus indicated in the original studies).

Results

Study descriptives

Our study includes 255 mother-child pairs from three subcohorts of the INMA project. Among the participant 255 mother-child pairs, 42 were from INMA Valencia, 111 from INMA Sabadell and 102 from INMA Gipuzkoa subcohorts. Table 1 shows the main sociodemographic characteristics, reproductive/birth and cognitive scores of the included mother-child pairs. Around half of the children were females and half were males (49.4% and 50.6%), and all were of European ancestry. For half of the mothers, it was their first pregnancy (56.86%), and the majority of them had a secondary education or superior. Almost 75% of the mothers did not smoke throughout the whole pregnancy and 14.12% smoked only during the first trimester. Mean age of the mother at the delivery date was 31.34 years (SD = 3.88) and the mean gestational age was 39.72 weeks (SD = 1.31). Average cell type proportions were similar to those reported by previous authors [48].

The child’s cognitive scores were assessed through the MSCA test at 4 years of age (mean = 4.68, SD = 0.49). The mean of the main study variables (verbal score, perceptive performance score, and general cognitive score) is close to 100 because they were standardised in the whole cohort population (N = 2150) to facilitate comparison among the different subcohorts [30] (Fig. 1; Table 1). Verbal and perceptive performance score weakly correlated with each other (r = 0.366), whereas general cognitive score strongly correlated with both verbal (r = 0.848) and perceptive performance score (r = 0.769) (Fig. 2). This is not unexpected since general cognitive score is calculated using the cognitive tasks included in both the verbal score and the perceptive performance score, while verbal score and perceptive performance score do not share any MSCA cognitive subtest.

Fig. 2
figure 2

Correlation between the verbal score, perceptive performance score, and general cognitive score.

Full size image

Epigenome-wide association studies

We run three separate EWAS, one for each cognitive domain, adjusting for covariates. We did not observe any relevant genomic inflation in the main three models: lambdas were 0.89, 0.97, and 0.92 for general cognitive score, verbal score, and perceptive performance score models, respectively. QQplots are shown in Supplementary Fig. 1. The analysis assessing the association between placenta DNAm and general cognitive score revealed 4 significant CpGs after correcting for FDR (p < 7.0611 × 10−8): cg00866476, which was negatively associated with the general cognitive score, and cg15480200, cg02986379, and cg14113931, which were positively associated with the general cognitive score (Fig. 3). Cg1548200 is annotated to an open sea region within chromosome 10, cg02986379 to a south shore region annotated to the DAB2 gene in chromosome 5, cg00866476 to a north shore region not annotated to any gene in chromosome 2, and cg14113931 to a north shore region annotated to CEP76 and PSMG2 genes in chromosome 18 (Table 2). When running the analysis separately by subcohort, we observed a consistent effect in all subcohorts (see Supplementary Fig. 2 for forest plots), showing heterogeneity (i2) values ranging from 0 to 0.2 and p > 0.05. None of the probes reached statistical significance in the analysis assessing the other two cognitive domains.

Fig. 3: Manhattan plots of the EWAS of DNAm and MSCA scores.
figure 3

a general cognitive score; b verbal score; and c perceptive performance score. Each dot represents one of the 708,105 CpGs. The y axis represents the negative logarithm of the association p values, and the x axis shows the genomic coordinates of the CpGs. The red abline indicates the FDR significance, while green abline indicates a suggestive cutoff of p < 1 × 10–4.

Full size image
Table 2 Statistically significant associations between DNAm and childhood general cognitive score after FDR correction.
Full size table

When exploring results at a suggestive significance threshold of p < 1 × 10−4, methylation at 103, 172, and 112 CpGs was associated with the general cognitive, the verbal, and the perceptive performance score scores, respectively (Fig. 3). See the annotated results of all probes that reached a significance of p < 1 × 10−4 or each of the main models (Supplementary Table 2). Complete results are available at Zenodo platform (https://zenodo.org/records/13379121).

When comparing all the CpGs that where significant at a suggestive p-value (p < 1 × 10−4) in the three main models, we observed that of the 172 suggestive CpGs for verbal score, 33 were common with the suggestive CpGs for general cognition. Of the 112 suggestive CpGs for perceptive performance score, 22 were common with CpGs for the general cognitive score. None of the suggestive CpGs were shared between models assessing the perceptive performance score and the verbal score. A Venn diagram summarising these findings is shown in Supplementary Fig. 3.

Results of the sensitivity analyses without adjusting for gestational age indicate that gestational age is not significantly affecting our main findings, which is consistent with the fact that gestational age is not associated with any of the 3 cognitive domains in our dataset (p > 0.05) (Lambdas ranging from 0.89 to 0.97) (summary statistics are shown in Supplementary Table 3, correlation between models is shown in Supplementary Fig. 4). Results of the sensitivity analyses with additional adjusting for maternal cognitive function indicate that the main findings are not influenced by the cognitive function of the mother (Lambdas ranging from 0.88 to 0.97) (descriptives are shown in Supplementary Table 4, correlation between models is shown in Supplementary Fig. 4, summary statitics are shown in Supplementary Table 5). Finally, results of the sensitivity analyses excluding pregnancy complications and adverse birth outcomes indicate that the main findings are not driven by the subset of mother-child pairs that experienced pregnancy complications and adverse birth outcomes (Lambdas ranging from 1.0 to 1.1) (descriptives are shown in Supplementary Table 6, correlation between models is shown in Supplementary Fig. 5, summary statitics are shown in Supplementary Table 7). Hence, when looking at the FDR-significant CpGs for the general cognitive score in the main model, we found that direction and magnitude of the effect was rather consistent between main and all the sensitivity models (with % change ranging from 0.1 to 10.18%).

Differentially methylated regions analysis

In the DMRs analysis, two DMRs containing 5 and 2 CpGs were associated with general cognitive score, two DMRs containing 7 and 3 CpGs were associated with perceptive performance score, and two other DMRs containing 8 and 2 CpGs were associated with verbal score (Table 3). One of the DMRs (chr18:71879836–71879850) was associated with both verbal and general cognitive scores. Notably, none of the DMRs included the CpGs that were significantly associated with the general cognitive score after FDR correction. As shown in Table 3 and Supplementary Fig. 6, DMR on chromosome 3 (chr3:168864101–168864417) and DMR on chromosome 6 (chr6:31803111–31803199), both associated with perceptive performance score, were annotated to a CpG island and a south shore within the MDS1 and EVI1 Complex Locus (MECOM) and Small Nucleolar RNA Host Gene 32 (SNHG32/C6orf48) genes, respectively; DMR on chromosome 19 (chr19:9785727–9786077), associated with verbal score, was annotated to a CpG island and a south shore within the Zinc Finger Protein 562 (ZNF562) gene. Moreover, according to ENCODE data on several cell lines, including different blood and embryonic stem cell types, these three DMRs present H3K27Ac histone marks and overlap with DNAse hypersensitive areas, which are usually associated to active regulatory elements. DMRs on chr6 and 19 overlap with active chromatin states (active promoters). In contrast, DMR on chr3 overlaps with polycomb-repressed states. The rest of DMRs were annotated to open sea regions in the genome.

Table 3 Statistically significant associations between placenta differentially methylated regions (DMRs) and childhood cognitive scores after FDR correction.
Full size table

Follow-up analysis

We conducted a set of downstream analysis aimed to facilitate biological interpretation of the main findings. First, we searched the FDR-significant CpGs of the general cognitive score single-probe analysis and CpGs located within the FDR-significant DMRs in the EWAS catalogue and EWAS atlas (see Supplementary Table 8a). According to these databases, methylation at 3 of the 4 CpGs associated with the general cognitive score have been previously associated with gestational age (in cord blood and fetal brain), age (in blood), exposure to environmental factors (smoking, high maternal plasma glucose, in blood), or associated to traits or health conditions (preterm birth, eosinophilia, in blood). Regarding the CpGs included in the DMRs (27 CpGs), several of them have been associated with age in blood (only one with gestational age), and traits or health conditions (rheumatoid arthritis, B acute lymphoblastic leukaemia, type 2 diabetes—in blood-, or renal, colorectal or oral cancer—in tumoral tissues) (Supplementary Table 8b).

Of the top FDR-significant CpGs, cg14113931, corresponds to a brain and blood methylation quantitative trailt locus (mQTL), and cg00866476 to a blood mQTL (Supplementary Table 8a). Notably, none of these 4 top CpGs were identified as FDR-significant in previous studies exploring the association between DNAm in cord blood or whole blood and cognitive functions in childhood or adulthood [22,23,24,25,26]. To further explore these tissue differences, we run the same analyses with cord blood DNAm data in the INMA cohort individuals that have methylation data on both placenta and cord blood (n = 148). They showed very different effect directions/sizes across the two tissues. Hence, methylation level means were also very different across the two tissues for cg15480200 and cg14113931, which were hemymethylated in placenta (Supplementary Table 9).

Of the FDR-significant DMRs CpGs, cg07751698 and cg22034155, associated with perceptive performance, correspond to fetal brain mQTLs (and a blood mQTL in the case of the latter). Four CpGs from the same DMR associated with general cognitive score (cg02712303, cg23190366, cg09433114, cg05720891), correspond to 4 blood mQTLs, one of them also being a blood eQTL regulating Brain-specific angiogenesis inhibitor 1-associated protein 2 –BAIAP2-AS1- gene expression (cg02712303), and another (cg23190366) also being a brain mQTL. Finally, regarding the longest DMR associated with verbal score, 3 out of the 8 included CpGs (cg24874111, cg20631204, and cg14073063) correspond to fetal brain mQTLs (cg24874111 and cg14073063 also being blood eQTLS regulating Zinc Finger Protein 562 –ZNF562- gene expression), while another CpG (cg17704570) corresponds to a blood mQTL (Supplementary Table 8b).

We also explored whether genomic regions of the FDR-significant hits overlapped with the 242, 226, and 1271 genome-wide significant loci previously reported in GWAS of intelligence [46], cognitive performance, and educational attainment [47] (0.5 Mb window centred to the CpG position). We found that the genomic regions of CpGs cg15480200 and cg14113931 include SNPs previously associated with educational attainment (rs12765185, located on genomic coordinate chr10: 134977077, and rs8097125, located on genomic coordinate chr18: 13005473, respectively [47]).

Hence, the genomic regions of the two DMRs associated with perceptive performance score include SNPs previously associated with educational attainment (rs13099165, located on genomic coordinate chr3: 168740425; rs9267658, located on genomic coordinate chr6: 31845985; rs9267677, located on genomic coordinate chr6: 31892641 [47]) and with cognitive performance (rs707938, located on genomic coordinate chr6: 31729359 [47]); and the genomic region of the previously mentioned DMR associated with verbal score includes SNPs previously associated with educational attainment (rs2287838, located on genomic coordinate chr19: 9959014 [47]).

The gene-set enrichment analysis of the suggestive CpGs (p < 1 × 10−4) conducted with Missmethyl did not show any significant results after multiple-testing correction (Supplementary Table 10; Supplementary Table 11). However, within the top biological pathways and processes there were some related with the formation of the nervous system or with the development of the fetus, as for example, the ventral spinal cord formation or the Notch pathway, and others related to the immune system, reparation, or biosynthesis. The enrichment analysis using the Consensus Path tool revealed an enrichment of 19, 9, and 5 gene-sets significant at FDR (p < 0.05) for the general cognitive score, verbal score, and perceptive performance score, respectively (Supplementary Table 12). The most recurrent gene-sets were related to Hedgehog and Notch signalling pathways, both involved in placenta development, fetus formation, and brain growth [49,50,51].

Finally, we checked whether the suggestive CpGs were over-represented in the placenta germline DMRs (gDMRs) (regions of the genome that show differential methylation according to the parental origin of the genome, associated with the expression of imprinted genes, which are expressed in a parent-of-origin-specific manner) [52], or in placenta partially methylated domains (PMDs) (large regions that cover ~37% of the human genome, which are partially methylated and contain placenta-specific repressed genes) [53] (Supplementary Table 13). Notably, there was a highly significant enrichment of suggestive CpGs of the general cognitive score in both gDMRs and PMDs, while there was not an enrichment of suggestive CpGs from the verbal score and perceptive performance score scores in these regions.

Discussion

In this study, we conducted, to the best of our knowledge, the first epigenome-wide-analysis investigating whether placenta DNAm profile is associated with childhood cognitive functions measured in three main domains of the McCarthy Scales, such as verbal score, perceptive performance score, and general cognitive score, at the age of 4 years. Moreover, we adjusted the models for confounders and explored the effect of maternal cognitive function, gestational age, pregnancy complications, and adverse birth outcomes in this association. The main analysis, at a single-probe level, did not identify significant associations between placenta DNAm and perceptive performance score or verbal score. However, we did find some evidence of FDR-corrected significant association between DNAm and the general cognitive score in 4 CpGs. Overall, the magnitude of the beta coefficients of the top CpGs associated with the general cognitive score is rather small (the difference in DNAm level % associated with an IQR increase in cognitive scores ranges from −0.8 to 1), and therefore, its functional implications remain unclear. Nonetheless, there are examples in the literature (for example, in blood samples) where small effects persist and have been replicated across populations and across time, as summarised by Breton and Marsit [54]. These small effect sizes, while challenging to interpret in isolation, may collectively contribute to significant biological and clinical outcomes. Another explanation for these small changes could be that DNAm is altered only in a specific cell type of the placenta, and this effect is diluted when analysing the bulk tissue.

Interestingly, one of the FDR-significant CpGs, cg02986379, is located in the promoter region of the DAB2 gene, which encodes a clathrin and cargo binding endocytic adaptor protein, playing a role in cellular trafficking of a number of transmembrane receptors and signalling proteins, cell differentiation, and proliferation. This gene is highly expressed in the placenta, where it seems to play a role in placenta growth and development [55], and its compromised abundance is associated with fetal growth restriction and low birth weight pathology. In the nervous system, DAB2 seems to have a role in neuronal development and nerve growth factor-mediated neurite outgrowth [56, 57]. Similarly, another of the FDR-significant CpGs, cg14113931, is located in the body and promoter of CEP76 and PSMG2 genes, respectively. It has not been described a particular role for CEP76 gene in the placenta, but this gene encodes an important centrosomal protein that regulates cell cycle and mitotic progression [58, 59]. PSMG2 encodes an adaptor protein involved in chaperone-mediated protein complex assembly, with no specific function in placenta, but associated with cognitive impairment and Alzheimer’s disease [60]. The other 2 top CpGs indentified, cg15480200 and cg00866476, are annotated to genomic open sea regions. Notably, none of these 4 top CpGs were identified as FDR-significant in previous studies exploring the association between DNAm in cord blood or whole blood and cognitive functions in childhood or adulthood [22,23,24,25,26]. This discrepancy in findings could be attributed to differences in tissue type (supported by the very different effect directions/sizes and methylation level means across cord blood and placenta DNAm at these CpGs in the same individuals from the INMA cohort), as well as other factors such as study design, population characteristics, or sample variability. Interestingly, genomic regions of 2 of these 4 CpGs include SNPs previously associated with educational attainment, a trait highly associated with intelligence test scores [61]. According to EWAS atlas and EWAS catalogue, 3 of these 4 top CpGs (cg15480200, cg00866476, and cg14113931) have been associated with gestational age and preterm birth in whole blood or cord blood studies [62,63,64]. Since preterm birth and reduced gestational age are associated with cognitive functions later in life, this may explain an association between DNAm at these particular sites and child’s cognition if mediated by adverse birth outcomes such as preterm birth. However, this is not supported by the high correlation between the main models and models unadjusted by gestational age, or models excluding pregnancy complications and adverse birth outcomes (preterm). It might be that gestational age affects this association at earlier stages than 4 years of age, and further investigation is required.

Similarly, maternal heritability did not seem to influence our results, given that our models adjusted by maternal cognitive function showed a high correlation with the main models. On the other hand, 3 of these 4 CpGs (cg15480200, cg00866476, cg14113931) were associated with blood mQTLs according to GoDMC database [44], and cg14113931 corresponds to a brain mQTL according to Brain xQTL Serve resource [43]. None of the CpGs correspond to previously reported placenta mQTLs. However, mQTL studies in the placenta are scarce and smaller than in the other tissues, so it is unclear whether genetics might be, at least partially, underlying these associations.

When comparing all the suggestive CpGs in the three main models, we observed that a small proportion were common between general cognitive score and verbal score, a slightly smaller proportion were common between cognitive general score and perceptive performance score, and none of the suggestive CpGs were shared between models assessing the perceptive performance score and verbal score scores. This is not unexpected since the general cognitive score is calculated using the cognitive tasks included in both the verbal score and the perceptive performance score, while the verbal score and perceptive performance score do not share any MSCA cognitive subtests. Actually, verbal and perceptive performance scores are weakly correlated with each other, and general cognitive scores are strongly correlated with both verbal and perceptive performance scores. The fact that the general cognitive score includes more tasks than the other two cognitive domains could be underlying the more significant results identified in this model compared to the other 2 for the single-probe analyses.

Functional enrichment of the suggestive CpGs in the three main models revealed gene sets involved in placenta development, fetus formation, and brain growth, such as Hedgehog and Notch pathways [49,50,51]. Epigenetic changes in these pathways might be associated with changes in placental function, fetal development, and specifically brain development, having a potential impact on later-in-life cognitive functions. Finally, the significant enrichment of suggestive CpGs from the general cognitive score model in both PMDs and gDMRs regions, suggests a specific regulatory role of these CpGs on placenta genes that are usually repressed in this tissue [53], and a potential involvement of these CpGs in genomic imprinting [65]. In fact, imprinted genes are strongly over-represented in the placenta, and it has been stated that changes in their expression in this tissue may also have consequences for brain in both the offspring and the mother, possibly through an imbalance in nutrients supply/demand or by alterations in the endocrine signalling from the placenta to the mother [66, 67]. However, the functional implications of this need further exploration, and functional and genomic enrichment analysis have to be interpreted with caution given the relaxed suggestive p-value used in these particular analyses (p < 1 × 10−4).

When exploring the association between child cognitive functions and DMRs in the placenta, we found 2 DMRs negatively associated with general cognitive score, 2 DMRs negatively associated with perceptive performance score, and 2 DMRs negatively associated with verbal score. Interestingly, three of them are annotated to regulatory regions of the following genes: MECOM, C6orf48/SNHG32, and ZNF562. First, MECOM or EVI1 encodes a zinc finger transcription factor, expressed in the placenta among other tissues, known to be involved in hematopoiesis, cell cycle regulation, differentiation, proliferation, and embryogenesis [68, 69]. This gene is also known as a negative regulator of NF-κB-dependent inflammation [70]. Neurodevelopmental disorders such as ASDs or cognitive impairment are associated with inflammation during development [71, 72]. This can be, in part, caused by a deregulation of the maternal-fetal immune environment [73, 74]. Infections during pregnancy are one of the main triggers of maternal–fetal immune dysregulation [73, 74]. In line with this, methylation levels at three CpGs included in this DMR (cg07751698, cg06238409, and cg13147822) were previously associated with the presence of several microorganisms in the placenta [75]. Second, ZNF562 is involved in transcription regulation and, according to the STRING webtool [76], appears to interact with SMARCAD1 and TRIM28, both of which are involved in the regulation of the inflammatory response [77, 78]. Interestingly, two CpGs included in this DMR are reported fetal brain mQTLs and eQTMs in childhood blood that regulate ZNF562 gene expression. However, few data have been published to date regarding this gene, especially in relation to placenta or fetal development. Third, the non-coding small nucleolar RNA C6orf48/SNHG32, of unknown function, is downregulated in cytotrophoblasts from severe preeclampsia cases [79]. Interestingly, a DMR within this gene in the placenta was associated with following a vegetarian diet before pregnancy [80]. Moreover, eight CpGs within this region in the placenta were previously associated with fetal growth restriction in mothers with insufficient gestational weight gain [81]. Overall, the fact that these DMRs overlap with regulatory elements of the associated genes supports the potential functional relevance of these regions, especially regarding the 2 DMRs associated with inflammation-related genes. Hence, similar to what was found in the single-CpG analysis, genomic regions of 3 of these DMRs include SNP’s previously associated with cognitive performance and with educational attainment (which is highly associated with intelligence test scores [61]). Finally, even though the DMR in chromosome 17 associated with general cognitive score is located in an open sea region, it includes 4 CpGs reported as mQTLs in blood, 1 of which is also a reported brain mQTL, and another an eQTM in child blood that regulates BAIAP2-AS1 gene expression, a gene that is altered in the brain of schizophrenia patients [82].

The current findings should be interpreted in the context of several limitations. The reduced sample size warrants future genome-wide studies with larger sample sizes to replicate and extend our findings, and to explore sex and cell type influences on the assessed associations. Since this study is based on a European population sample, this should be replicated in other ancestries and in more diverse settings. Additionally, eventhough we used a genome-wide array to measure DNAm, this only accounts for only 2–3% of the CpG sites in the genome. We could only compare our results with the other two EWAS on cognitive functions performed in cord blood and whole blood (rather than placenta) because of the novelty of the study. Finally, although our findings are based on placenta tissue, we cannot discard the fact that the methylation status of these marks in the placenta could also reflect the same status in other target tissues, such as the fetal brain, which are, similarly, or even more relevant for this association. In line with this, we detected that 6 of the CpGs included in the FDR-significant DMRs correspond to mQTLs of adult or fetal brains [42, 43].

Strengths of this study include investigating the methylome in relation to future cognitive functions in a crucial tissue for embryonic fetal development, and in particular for brain development, such as the placenta; and the use of detailed information on factors that influence neurodevelopment to adjust the models and to perform sensitivity analyses. Moreover, to the best of our knowledge, this is the first time that a study on placenta has attempted to establish a link between DNAm and childhood cognitive function, which may provide new avenues for other researchers to investigate this topic in the future.

In summary, we found some evidence of placenta DNAm being associated with childhood cognitive functions. This association was not affected by maternal cognitive function, gestational age, pregnancy complications, and adverse birth outcomes in our cohort. Top CpGs and DMRs are mapped to genes involved in placenta fetal and brain development, and inflammation. Moreover, some are located close to loci of cognitive performance and educational attainment. These findings suggest that placental DNAm could be a mechanism contributing to the alteration of important pathways in the placenta that have a consequence on the offspring’s brain development and cognitive function. However, the limitations of the study, such as the reduced sample size, call for further research in other cohorts (including other ancestries, social contexts, sex-stratified analyses, etc), and in vivo and in vitro experiments to achieve a more comprehensive understanding of this association and its biological implications.