A Western Dietary Pattern during Pregnancy is Associated with Neurodevelopmental Disorders in Childhood and Adolescence
=======================================================================================================================

* David Horner
* Jens Richardt M. Jepsen
* Bo Chawes
* Kristina Aagaard
* Julie B. Rosenberg
* Parisa Mohammadzadeh
* Astrid Sevelsted
* Nilo Følsgaard
* Rebecca Vinding
* Birgitte Fagerlund
* Christos Pantelis
* Niels Bilenberg
* Casper-Emil T. Pedersen
* Anders Eliasen
* Yulu Chen
* Nicole Prince
* Su H. Chu
* Rachel S. Kelly
* Jessica Lasky-Su
* Thorhallur I. Halldorsson
* Marin Strøm
* Katrine Strandberg-Larsen
* Sjurdur F. Olsen
* Birte Y. Glenthøj
* Klaus Bønnelykke
* Bjørn H. Ebdrup
* Jakob Stokholm
* Morten Arendt Rasmussen

## ABSTRACT

Despite the high prevalence of neurodevelopmental disorders, there is a notable gap in clinical studies exploring the impact of maternal diet during pregnancy on child neurodevelopment. This observational clinical study examined the association between pregnancy dietary patterns and neurodevelopmental disorders, as well as their symptoms, in a prospective cohort of 10-year-old children (n=508). Data-driven dietary patterns were derived from self-reported food frequency questionnaires. A Western dietary pattern in pregnancy (per SD change) was significantly associated with attention-deficit / hyperactivity disorder (ADHD) (OR 1.66 [1.21 - 2.27], p=0.002) and autism diagnosis (OR 2.22 [1.33 - 3.74], p=0.002) and associated symptoms (p<0.001). Findings for ADHD were validated in three large (n=59725, n=656, n=348), independent mother-child cohorts. Objective blood metabolome modelling at 24 weeks gestation identified 15 causally mediating metabolites which significantly improved ADHD prediction in external validation. Temporal analyses across five blood metabolome timepoints in two independent mother-child cohorts revealed that the association of Western dietary pattern metabolite scores with neurodevelopmental outcomes was consistently significant in early to mid-pregnancy, independent of later child timepoints. These findings underscore the importance of early intervention and provide robust evidence for targeted prenatal dietary interventions to prevent neurodevelopmental disorders in children.

## INTRODUCTION

Neurodevelopmental disorders, particularly attention-deficit/hyperactivity disorder (ADHD) and autism, are prevalent and an increasing public health concern (1). In 2020, Danish registers reported a cumulative incidence of 5.9% in males and 3% in females for ADHD, and 4.3% in males and 1.8% in females for autism during childhood and adolescence (2). Large meta-analyses support these findings, with global estimates of 5.9% for ADHD (3) and up to 1.1% for autism (4). In addition, neurodevelopmental disorder symptoms manifest as continuous traits within the paediatric general population (5), and can be reliably captured via validated questionnaires (6,7). Studies have linked several prenatal exposures to neurodevelopmental disorders in children, including maternal obesity, metabolic disturbances, stress, and individual nutritional elements (8). Numerous animal studies have shown that high-fat diets can alter brain morphology (9) and behaviour in offspring in ways resembling neurodevelopmental disorders (10,11), however these findings may not extrapolate well to humans. Dietary constituents may contribute directly to the aforementioned aetiologies, providing both the energy substrates and building blocks for the foetal brain (12). Given the potential public health impact, it is vital to robustly investigate the link between maternal diet during pregnancy and childhood neurodevelopmental disorders, particularly those involving Western diets characterised by high consumption of processed meats, sugars, refined grains, and low intakes of fruits and vegetables. These diets, which are prevalent in modern societies, starkly contrast with historical human diets and may influence developmental outcomes (13). Existing literature lacks objective measurements of dietary patterns in pregnancy, and adjustment for confounding dietary patterns in childhood, when assessing this relationship.

The placenta transfers nutrients from the mother to the developing foetus during pregnancy (14). These nutrients, which include the essential n-3 long chain polyunsaturated fatty acids (n3-LCPUFA) and micronutrients such as iron, choline, iodine, zinc, and vitamins B, D, and E, are necessary for foetal brain development and function (12,15), and are obtained as part of a healthy dietary pattern. While the precise mechanisms underlying the influence of maternal diet on offspring neurodevelopment remain to be fully elucidated, various factors may be implicated in developmental processes. Among these factors are the potential impacts of dietary compounds (16) and lipid profiles (17). For example, certain dietary-derived metabolites may interact with developmental pathways in the foetus, potentially affecting neurodevelopmental outcomes (18). Additionally, higher maternal intakes of n3-LCPUFAs, such as docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA), which are associated with better overall diet quality (19), are reported to be neuroprotective (20).

Genetics contribute significantly to neurodevelopmental disorders, with heritability estimates as high as 80% (21). The increasing prevalence of these disorders (22) highlight the need for a better understanding of how environmental factors interact with genetic risk (23). In this context, pregnancy dietary influences may be moderated by the child’s underlying genetic risk for neurodevelopmental disorders, which can be succinctly captured in a polygenic risk score (PRS) (24). While twin studies consider gene-environment interactions when estimating heritability, clinical evidence validating this potential interaction is notably lacking (25).

In this study, we hypothesise that a Western dietary pattern in pregnancy is associated with adverse neurodevelopmental outcomes. We sought to bridge the research gap by leveraging the thorough neurodevelopmental clinical examinations within the Copenhagen Prospective Studies on Asthma in Childhood 2010 (COPSAC2010) mother-child cohort. We employ external validation approaches in three independent mother-child cohorts. Drawing from twelve metabolome datasets from three mother-child cohorts, we juxtapose maternal and child blood profiles to shed light on the temporality of dietary associations and their potential impact on neurodevelopment and delve into which metabolites may be driving dietary associations on neurodevelopment.

## RESULTS

### Cohort Characteristics

During the 10-year clinical visit, a total of 593 children (84.7%) participating in the COpenhagen Prospective Study on Neuro-PSYCHiatric Development study (COPSYCH) underwent clinical examination for neurodevelopmental and other psychiatric disorders. Additionally, 11 children who did not complete the clinical examination, had information regarding their neurodevelopmental symptom loads. There were no significant differences in baseline sociodemographic characteristics between participants and non-participants at the 10-year visit (*Table S1*), nor between male and female children (*Table S2*).

View this table:
[Supplementary Table S1.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T3)

Supplementary Table S1. 
Cohort characteristics stratified by participation in the COPSAC2010 10 year clinical visit. Low, medium, and high income are defined as <50,000, 50,000 - 110,000, and >110,000 euro respectively. Low, medium and high level educational attainment are defined as “primary, secondary, or college graduate”, “tradesman or bachelor’s degree” and “Masters degree”, respectively.

View this table:
[Supplementary Table S2.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T4)

Supplementary Table S2. 
Cohort characteristics stratified by child sex in COPSYCH 10 year clinical visit. Low, medium, and high income are defined as *<50,000, 50,000 - 110,000, and >110,000 euro,* respectively. Low, medium and high level educational attainment are defined as “*primary, secondary, or college graduate*”, *“tradesman or bachelor’s degree*” and “*Masters degree*”, respectively.

We included 508 children with both pregnancy dietary and clinical data in our analysis, 77 (15.2%, 71% male) had any neurodevelopmental disorder. Specifically, 55 (10.8%, 76% male) were diagnosed with ADHD (25 ADHD predominantly inattentive presentation and 30 ADHD combined presentation) and 13 (2.6%, 69% male) with autism.

### Identification of Dietary Patterns

We used principal component analysis (PCA) on 95 nutrient constituents (*Table S3)* from pregnancy food frequency questionnaires (FFQ) assessed at 24 weeks gestation to identify maternal dietary patterns in the COPSAC2010 cohort (*Figure S1*). Principal component 1 (PC1), which explained 44.3% of the variance, had a positive association across all food groups and represents a "Varied dietary pattern”. PC2, which explained 10.7% of the variance, had positive associations with intakes of animal fats, refined grains and high-energy drinks, and negative associations with intakes of fruit, fish, and vegetables, representing a "Western dietary pattern" (*Figure 1*). Regarding macronutrient intake, PC2 (Western dietary pattern) predominantly reflected a higher intake of fats (*Figure S2A*), specifically saturated fatty acids (*Figure S2B)*. Using the maternal PC model, we predicted a child’s Western dietary pattern at 10-years, allowing for a consistent comparison of dietary habits (r=0.22). A Western dietary pattern during pregnancy was negatively associated with social circumstances and positively associated with maternal pre-pregnancy body mass index (BMI), smoking during pregnancy, antibiotic use during pregnancy and a Western dietary pattern in children at 10-years of age (*Table 1*).

![Supplementary Figure S1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F5.medium.gif)

[Supplementary Figure S1.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F5)

Supplementary Figure S1. 
Biplot of the first two principal components (PCs) from the maternal food frequency questionnaire-derived nutrient constituents. PC2 (Western dietary pattern) is significantly associated with any neurodevelopmental disorder (OR 1.53 [1.17 - 2.00], p=0.002), ADHD (OR 1.66 [1.21 - 2.27], p = 0.002), and autism diagnosis (OR 2.22 [1.33 - 3.74], p = 0.002). PC1 is not significantly associated with any neurodevelopmental disorders (p > 0.288). Further details of the loadings for each of the 95 nutrient constituents can be viewed in supplementary table 2.

![Supplementary Figure S2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F6.medium.gif)

[Supplementary Figure S2.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F6)

Supplementary Figure S2. 
Biplot of the first two principal components from the maternal FFQ-derived nutrient constituents **A)** Nutrient constituents are categorised into fatty acids, amino acids, sugars, minerals and vitamins. Fatty acids are a key determinant of PC2 (Western dietary pattern) **B)** Stratified further by fatty acid type, saturated fatty acids are most associated with PC2 (Western dietary pattern).

View this table:
[Supplementary Table S3.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T5)

Supplementary Table S3. 
Descriptive summary statistics for nutrients derived from maternal food frequency questionnaires in pregnancy and principal component loadings for principal component 1 and 2 in the COPSAC2010 cohort.

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F1)

Figure 1. 
Figure showing the pregnancy Varied and Western dietary patterns derived from nutrient constituents (PC1 and PC2) and their associations with food groups. Multivariable analysis shows the Western dietary pattern is highly associated with any neurodevelopmental disorder diagnosis (OR 1.53, [1.17 - 2] (p = 0.002), ADHD diagnosis 1.66 [1.21 - 2.27] (p = 0.002), and autism diagnosis 2.22 [1.33 - 3.74] (p = 0.002), but the Varied dietary pattern is not significantly associated with neurodevelopmental disorders (p>0.288). Yellow represents positive associations with PC1, while red and green denote positive and negative associations with PC2, respectively.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T1)

Table 1. 
Model covariates stratified by low, middle and high (tertiles) for the pregnancy Western dietary pattern, in children who were clinically assessed for neurodevelopmental disorders.

**Footnote:**

† PC2: Principal Component 2, BMI: Body Mass Index, RCT: Randomised Controlled Trial, FFQ: Food Frequency Questionnaire, PRS: Polygenic Risk Score.

### A Western Dietary Pattern during Pregnancy is Associated with Neurodevelopmental Disorders and Symptom Loads in Children

A Western dietary pattern during pregnancy (PC2, per SD change) was significantly associated with any neurodevelopmental disorder (OR 1.53 [1.17-2.00], p=0.002), ADHD (OR 1.66 [1.21-2.27], p=0.002), and autism (OR 2.22 [1.33-3.74], p=0.002) (*Table 2*) in multivariate modelling adjusted for pre-pregnancy maternal BMI, social circumstances, child sex, birth weight, gestational age, pregnancy smoking/antibiotic use, pre-eclampsia and a child Western dietary pattern. It was also significantly associated with symptom loads for ADHD (1.73 [0.98-2.49], p<0.001) and autism (3.21 [1.69-4.74], p<0.001), with consistent associations observed in an earlier assessment of ADHD symptom load at 8 years (1.90 [0.42-4.55], p<0.001). Moreover, ADHD symptoms at 6 years, assessed using the hyperactivity/inattention scale from the Strengths and Difficulties Questionnaire (SDQ), were also significant (0.50 [0.29-0.70], p<0.001). To illustrate the associations between the Western dietary pattern during pregnancy, neurodevelopmental outcomes, and model covariates, we employed Gaussian graphical models (*Figure 2*).

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F2)

Figure 2. 
Graphical models illustrate associations (partial correlations p<0.05) between the Western dietary pattern during pregnancy, neurodevelopmental outcomes A: ADHD diagnosis, B: ADHD Symptom Load (ADHD-RS), C: Autism diagnosis, D: Autism Symptom Load (SRS-2), and model covariates. Models include child ADHD PRS (A&B) and autism PRS (C&D) as additional covariates.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T2)

Table 2. 
Regression analysis of the Western dietary pattern during pregnancy and neurodevelopmental outcomes. Odds ratios and estimates are interpreted as the association of 1 standard deviation of the pregnancy Western dietary pattern. **Model A**: Univariate analysis. **Model B**: Multivariable analysis (multivariable model adjusted for pre-pregnancy maternal BMI, social circumstances (household income, maternal education level, and maternal age), child sex, birth weight, gestational age, pregnancy smoking/antibiotic use, pre-eclampsia and a child Western dietary pattern at 10-years of age).

Our findings remained significant after further adjustment for maternal genetic risk for ADHD, total maternal energy intake, genomic principal components, intake of prenatal supplements (*Table S4*), and two prenatal nutrition supplementation randomised controlled trials (n3-LCPUFA and high-dose D vitamin) (*Table S5*). Associations of the pregnancy Western dietary pattern were significant for both male (1.42 [1.02-1.99], p=0.039) and female children (1.86 [1.15-3.07], p=0.012) for any neurodevelopmental disorder diagnosis (*Table S6*). Neurodevelopmental outcomes showed no significant differences related to child sex (p-interaction>0.12). Associations were comparable for both ADHD-predominantly inattentive presentation (OR 1.71 [1.13-2.56], p=0.009) and ADHD-combined presentation (OR 1.51 [1.03-2.20], p=0.033). In a subanalysis in children without a neurodevelopmental disorder diagnosis (n=428), the pregnancy Western dietary pattern association with ADHD symptom load remained significant (0.63 [0.01-1.25], p=0.048) (*Table S7*).

View this table:
[Supplementary Table S4.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T6)

Supplementary Table S4. 
Daily estimated intake of prenatal supplements assessed at 24 weeks gestation, with 1 month recall, including vitamins, minerals, and fatty acids, with mean (SD) values and number of participants taking each supplement.

View this table:
[Supplementary Table S5.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T7)

Supplementary Table S5. 
Subanalysis with extra adjustments for the pregnancy Western dietary patterns effect on neurodevelopmental outcomes (multivariable models). **Model 1:** Adjusted for prenatal RCT n-3 LCPUFA. **Model 2:** Adjusted for prenatal RCT high-dose Vitamin D. **Model 3:** Adjusted for prenatal supplements of vitamins, minerals and fatty acids (P1-PC6). **Model 4:** Adjusted for food frequency derived estimated energy intake. **Model 5:** Adjusted for maternal ADHD genetic polygenic risk score **Model 6:** Adjusted for Genomic Principal Components.

View this table:
[Supplementary Table S6.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T8)

Supplementary Table S6. 
Sex stratification and statistical interaction effect of the pregnancy Western dietary pattern effect on neurodevelopmental outcomes (multivariable models).

View this table:
[Supplementary Table S7.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T9)

Supplementary Table S7. 
The pregnancy Western dietary patterns effect on ADHD diagnosis subtypes, and subscores of ADHD and autism symptoms loads (multivariable models). **Column 1**: The entire COPSAC2010 **Column 2:** Subanalysis with any neurodevelopmental disorder diagnosis excluded.

### Validating our findings for ADHD in the Danish National Birth Cohort

We sought replication of our findings in the Danish National Birth Cohort (DNBC), which had an identical FFQ assessment undertaken at a comparable gestational age, neurodevelopmental diagnoses derived from national registries, and ADHD symptoms from the hyperactivity/inattention scale of the SDQ questionnaire. In our replication analysis including 59,725 mother-child pairs, we observed prevalences of 6.8% for any neurodevelopmental disorder, 4.6% for ADHD diagnosis, and 2.7% for autism diagnosis. We found that the association of the Western dietary pattern in pregnancy, trained in the COPSAC2010 cohort, was significantly associated with ADHD diagnosis (OR 1.07 [1.02-1.11], p=0.002), but not for any neurodevelopmental disorder diagnosis (OR 1.03 [1.00-1.07], p=0.087) or autism diagnosis (OR 1.04 [0.98-1.09], p=0.351), in multivariable analysis (*Table S8*). Reinforcing our method, we observed PC2 trained within the DNBC cohort was also significantly associated with ADHD diagnosis (OR 1.05 [1.01-1.10], p=0.018), with a strong correlation between PC2 trained in COPSAC2010 and the DNBC (r=0.91). Moreover, ADHD symptoms at 7 years, assessed in a subset of 37,608 children, were also significantly associated with the Western dietary pattern in pregnancy, trained in COPSAC2010 (0.11 [0.09-0.13] p<0.001), and within the DNBC (0.11 [0.09-0.13] p<0.001), in multivariable analysis.

View this table:
[Supplementary Table S8.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T10)

Supplementary Table S8. 
External validation of the association between the pregnancy Western dietary pattern (PC2), trained within COPSAC2010, and within the DNBC and registry derived neurodevelopmental diagnosis in the DNBC cohort.

### A Blood Metabolomics-Derived Pregnancy Western Dietary Pattern Score Is Associated with Neurodevelopmental Disorders in COPSAC2010

We used the blood metabolome from week 24 gestation to objectively measure the pregnancy Western dietary pattern and validate our self-reported FFQ findings. 43.0% of 760 metabolites were significantly associated with the pregnancy Western dietary pattern, with 34.5% surviving false discovery rate (q<0.05). In blood metabolome modelling of the pregnancy Western dietary pattern, 43 metabolites emerged as biomarkers (RMSECV=0.87, R2CV=0.24). Loadings for these metabolites can be seen in *Figure 3* and *Table S9*.

View this table:
[Supplementary Table S9.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T11)

Supplementary Table S9. 
The 43 metabolites surviving regularisation in the pregnancy Western dietary patterns sparse partial least squares model in COPSAC2010. Metabolites name, super pathway, sub pathways,loadings for component 1 and whether or not the metabolite significantly mediated the association between the Western dietary pattern and any Neurodevelopmental Disorder Diagnosis.

![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F3)

Figure 3. 
The 43 metabolites, selected by the sparse partial least squares model, represent those associated with the Western dietary pattern during pregnancy at 24 weeks gestation in COPSAC2010. Positive metabolite scores indicate a positive association, while negative scores suggest an inverse relationship with this dietary pattern. A systematic backward elimination pinpointed 15 metabolites as mediators between the Western dietary pattern and any neurodevelopmental disorders. Notably, dietary-derived compounds, like ergothioneine, suggest potential protective roles, while certain lipid-associated metabolites hint at possible detrimental impacts on neurodevelopment.

A Western dietary pattern-metabolite score (WDP-MS) was derived from this model and was associated (per SD change) with any neurodevelopmental disorder diagnosis (OR 1.42 [1.07-1.90], p=0.017), and ADHD diagnosis (OR 1.39 [1.00-1.95], p<0.05), but not significantly with autism diagnosis (OR 1.09 [0.59-2.00], p=0.781).

To discern whether any of the 43 metabolites acted as biomarkers, or genuine mediators in the association between the Western dietary pattern and neurodevelopmental disorders, we deployed a backward elimination strategy. Our multivariable mediation analysis, found that 15 metabolites (*Figure 3, Table S9*) mediated the association between the Western dietary pattern and any neurodevelopmental disorder diagnosis. The Western dietary pattern had a total effect of 5.5% on any neurodevelopmental disorder diagnosis, of which 80.3% (p=0.003) was mediated through the selected 15 metabolites. It became evident that a predominant portion of the metabolites acting as mediators were derived from dietary components and lipids. Prominently, plant-based metabolites such as ergothioneine, carotene diol, and tartrate indicate potential protective effects. Conversely, many lipid-associated metabolites, such as medium chain fatty acids like caprate and caprylate, predominantly exhibited positive loadings, underscoring their possible detrimental impact on neurodevelopment.

We explored potential time-dependent associations by predicting an additional four WDP-MS using comparable blood metabolome datasets in mothers (one week postpartum) and children (6 months, 18 months, and 6 years). In a sensitivity analysis, which compared the five WDP-MS without including other covariates, significant associations were observed solely with the 24 week gestational WDP-MS for any neurodevelopmental disorder diagnosis (OR 1.50 [1.07-2.11], p=0.020), ADHD diagnosis (OR 1.58 [1.08-2.33], p=0.020) and directionally for autism diagnosis (OR 1.74 [0.84-3.57], p=0.131). This was the case despite the strong correlations noted between maternal dietary metabolite scores (r=0.68, *Table S10*). Findings were consistent for ADHD (1.43 [0.41-2.45], p=0.006) and autism symptom loads (2.86 [0.84-4.89], p=0.006). We substantiated these temporal associations by employing linear mixed models. Through this approach, we discerned that the initial magnitude of the Western dietary pattern during pregnancy, reflected by the model’s intercept, likewise served as a predictor for any neurodevelopmental disorder (OR 1.56 [1.24-1.97], p<0.001), ADHD diagnosis (OR 1.63 [1.25-2.12], p<0.001) but not autism diagnosis (OR 1.34 [0.813-2.18], p=0.242). Additionally, there were likewise significant associations with ADHD (1.69 [0.99-2.40], p<0.001) and autism symptom loads (3.65 [2.24-5.06], p<0.001).

View this table:
[Supplementary Table S10.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T12)

Supplementary Table S10. 
Correlation matrices of the respective maternal and child Western dietary pattern Metabolite Scores in the COPSAC2010 and VDAART cohorts.

### Western Dietary Pattern-Metabolite Score Predicts ADHD Diagnosis in the VDAART Mother-Child Cohort

As our blood metabolome modelling of the Western dietary pattern during pregnancy was associated with ADHD diagnosis and symptom load in COPSAC2010, we sought replication of our findings in the Vitamin D Antenatal Asthma Reduction Trial cohort (VDAART), a large US-based mother-child cohort. We included participants with corresponding pregnancy blood metabolome and ADHD outcome data (*n=656*). Among the children in the VDAART cohort, caregivers reported a diagnosis for ADHD or attention deficit disorder for 18 children (2.7%) by age 6 years and 57 children (8.7%) by age 8 years. Significant cohort characteristic differences were observed between the COPSAC2010 and VDAART cohorts (*Table S11*), encompassing age, income, maternal education, maternal age, gestational age, maternal smoking, race, maternal pre-pregnancy BMI, birthweight, and caesarean section delivery (p<0.001). Supplementary Figures 3 shows the metabolomic variability in the COPSAC2010 and VDAART, illustrating closely matched variances across pregnancy samples.

View this table:
[Supplementary Table S11.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T13)

Supplementary Table S11. 
Cohort characteristics stratified by cohorts. P-values summarise cohort characteristic differences between COPSAC2010/VDAART and COPSAC2010/COPSAC2000 respectively. Low, medium, and high income are defined as *<50,000, 50,000 - 110,000, and >110,000 Euro* in COPSAC2010/COPSAC2000, and *<30,000, 30,000 - 99,999, and >$100,000 USD* in VDAART, respectively. Low, medium and high level educational attainment are defined as “*primary, secondary, or college graduate*”, *“tradesman or bachelor’s degree*” and “*Masters degree*” in COPSAC2010/COPSAC2000 and “*Did not graduate from high school/Graduated from high school*”, “*Technical school/Junior college/some college*” and “*Graduate school/College graduate*” in VDAART.

View this table:
[Supplementary Table S12.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T14)

Supplementary Table S12. 
Maternal pre-pregnancy BMI and child sex associations on neurodevelopmental outcomes, with the Western dietary pattern included in multivariable models, shows the independent associations of maternal pre-pregnancy BMI and child sex on neurodevelopmental diagnosis and symptom loads. Furthermore, models show the associations with child ADHD and autism polygenic risk scores (PRS), likewise also adjusted for the Western dietary pattern, on neurodevelopmental diagnosis and symptom loads.

View this table:
[Supplementary Table S13.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/T15)

Supplementary Table S13. 
Analysis of missing metabolites and their association with neurodevelopmental outcomes. Missingness rates of all metabolites and the selected 43 metabolites in the pregnancy Western dietary pattern metabolite score, alongside their odds ratios and confidence intervals for neurodevelopmental disorders, symptom loads and the Western dietary pattern.

In the VDAART cohort, which includes blood metabolome profiling at two pregnancy timepoints during early (10-18 week) and late gestation (32-38 weeks), we used COPSAC2010 cohort-trained models to predict two WDP-MS to assess their association with ADHD diagnosis. The metabolite scores at both time points were positively associated with intakes of deep-fried foods, processed meats and margarine, and negatively associated with intakes of various vegetables, whole grain foods and fruit in the VDAART cohort, as assessed by independent food frequency questionnaires (*Figure S4*). There was a strong correlation between the VDAART FFQ’s principal component 1 and the predicted WDP-MS during early (r=0.48, p=1.4×10⁻⁴⁶) and late (r=0.46, p=2.3×10⁻⁴²) pregnancy, underscoring the WDP-MS validity in predicting dietary patterns across diverse populations.

![Supplementary Figure S3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F7.medium.gif)

[Supplementary Figure S3.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F7)

Supplementary Figure S3. 
Comparison of maternal blood metabolomes at three different pregnancy timepoints from two mother-child cohorts. Top panels show a Principal Component Analysis (PCA) scoreplot on all metabolites, as well as the selected metabolite scores for the COPSAC2010 vs VDAART 10-18 weeks (left) and COPSAC2010 vs VDAART 32-38 weeks (right). Bottom panels show the relative variation per metabolite computed by the ratio of sums of squares (SSQtime / SSQresidual) from a one-way anova model with Time/Cohort at predictor and Comparison of per metabolite standard deviation within cohort relative to 24 week gestation pregnancy time point from COPSAC2010 for the COPSAC2010 vs VDAART 10-18 weeks (left) and COPSAC2010 vs VDAART 32-38 weeks (right).

![Supplementary Figure S4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F8.medium.gif)

[Supplementary Figure S4.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F8)

Supplementary Figure S4. 
Associations between dietary intake from food frequency questionnaires (FFQs) during pregnancy and Western Dietary Pattern Metabolite Scores (WDP-MS) at two distinct timepoints (10-18 weeks and 32-38 weeks) in the VDAART cohort. The depicted associations, based on the VDAART cohort, were assessed using COPSAC2010 cohort-trained models that shared overlapping metabolites.

Despite significant differences in cohort characteristics, our blood metabolomics analysis successfully replicated the univariate association between the WDP-MS and ADHD diagnosis at age 6 in both early (OR 2.64 [1.56-4.64], p<0.001) and late pregnancy (OR 1.85 [1.12-3.16], p=0.020), in the VDAART cohort. After adjusting for similar covariates to those used in the COPSAC2010 cohort, a significant association remained between the WDP-MS and ADHD diagnosis in early pregnancy (OR 2.33 [1.26-4.42], p=0.008) but not in late pregnancy (OR 1.67 [0.90-3.20], p=0.113). For further robustness of our findings, we assessed children ever having had a ADHD diagnosis up to age 8. The WDP-MS successfully univariately replicated, albeit with attenuated estimates in early (OR 1.71 [1.28-2.32], p<0.001) and late pregnancy (OR 1.37 [1.04-1.83], p=0.030). Associations remained significant after multivariable adjustment in early pregnancy (OR 1.45 [1.02-2.06], p=0.039) but not late pregnancy (OR 1.11 [0.79-1.58], p=0.544).

In sensitivity analysis, we independently employed the backward elimination strategy, previously outlined for the COPSAC2010 cohort, on the metabolites overlapping with the VDAART 10-18 weeks and 32-38 weeks time points. In doing so, we identified the mediating metabolites in the association between the Western dietary pattern in COPSAC2010 to ADHD diagnosis (*Figure S5*). Leveraging COPSAC2010’s mediating metabolites, we performed a sensitivity analysis to contrast the performance of a WDP-MS using only these potentially causally mediating metabolites against the WDP-MS with all metabolites, for predicting ADHD in the VDAART cohort. At the 10-18 week timepoint, in a multivariable model, the subset of mediating metabolites improved ADHD diagnosis prediction at 6 years (OR 2.85 [1.60-5.20], p<0.001), significantly moreso (nested model chi-squared test p=0.012) than our model with all 18 metabolites. However, at the 32-38 week timepoint the subset of mediating metabolites (OR 1.81 [0.95-3.58], p=0.078) did not perform significantly better than the WDP-MS with all 19 metabolites (nested model chi-squared test p=0.423). To replicate the potential time-dependent association of the Western dietary pattern, we predicted three additional WDP-MS for children in VDAART (1 year, 3 years, and 6 years). When incorporating the three child WDP-MS in a model with the two pregnancy WDP-MS timepoints, we observed significant associations only with the early pregnancy WDP-MS for ADHD diagnosis at age 6 (OR 1.97 [1.01-3.80], p = 0.044) and ADHD diagnosis up to age 8 (OR 1.65 [1.11-2.45], p=0.013).

![Supplementary Figure S5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F9.medium.gif)

[Supplementary Figure S5.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F9)

Supplementary Figure S5. 
Metabolites associated with a Western dietary pattern at VDAART gestational timepoints (10-18 weeks and 32-38 weeks), based on models from the COPSAC2010 cohort. The selected metabolite scores are depicted by bars, distinguished by their metabolic pathway. Striped bars indicate metabolites that mediate the association with ADHD Diagnosis in COPSAC2010, while solid bars signify non-mediating metabolites. The directionality of the bars represents the positive or negative metabolite score. Metabolite scores from VDAART are colour-coded for each timepoint: 10-18 weeks (green) and 32-38 weeks (purple).

### Western Dietary Pattern-Metabolite Score Predicts ADHD Symptoms in COPSAC2000

We sought further replication in the COPSAC2000 mother-child cohort, where we included 328 participants with both neonatal dried blood spots (DBS) metabolic profiling and ADHD symptom load from the Adult ADHD Self-Report Scale (ASRS) at age 18. We developed a Western dietary pattern dry blood spot (WDP-DBS) score using overlapping DBS metabolites from the COPSAC2010 cohort to predict the association with ADHD symptom load. Supplementary Figure 6 illustrates the metabolomic variability between the COPSAC2010 and COPSAC2000 newborn samples, demonstrating comparable variance across cohorts.

To substantiate the WDP-DBS score’s relevance, we first confirmed its significance within the COPSAC2010 cohort, showing a strong association with ADHD symptom load, in multivariable modelling (1.19 [0.41-1.98], p=0.003), before applying it to COPSAC2000 data. In COPSAC2000, we successfully replicated this association between the WDP-DBS and ADHD symptom load (0.44 [0.01-0.86], p=0.045), but this finding didn’t survive covariate adjustment (0.34 [-0.08-0.77], p=0.114). In sensitivity analysis, we enhanced cross-correlations between COPSAC2010 and COPSAC2000 cohorts from >0.4 to >0.6, and restricted the analysis to 626 metabolites with similar distributions between cohorts (Kolmogorov-Smirnov p ≥ 0.05). This approach led to more pronounced univariate (0.49 [0.06-0.91], p=0.027) and multivariable associations (0.41 [-0.02-0.83], p=0.062).

### Moderating Effects of Genetics, Maternal BMI, and Child Sex on Diet-Neurodevelopment Associations

In a subanalysis, we investigated if established risk factors (child genetic risk (ADHD and autism PRS), maternal pre-pregnancy BMI, and child sex) moderate the association of the Western dietary pattern during pregnancy with neurodevelopmental outcomes in COPSAC2010. Supplementary Table 12 provides more details on the independent associations of these factors on neurodevelopmental outcomes.

ADHD (11.8%) and autism (2.7%) diagnosis, had limited variability in our dataset, prompting us to assess the association of the Western dietary pattern using strata of median split maternal pre-pregnancy BMI (above and below 23.7) and PRS scores (*Figure 4A/B*). In multivariable analyses, children born to mothers with a higher BMI and with a greater genetic predisposition for the respective disorder showed the strongest associations with the Western dietary pattern. For ADHD diagnosis, this group had an odds ratio of 2.18 [1.30-3.74], p=0.004, and for autism diagnosis, it was 4.59 [2.33-9.78], p<10⁻4. When further stratified by child sex, findings continued to show consistent significance across both sexes (*Figure S7A/B*). Findings were again consistent when considering the objectively measured WDP-MS (*Figure S8A/B*).

![Supplementary Figure S6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F10.medium.gif)

[Supplementary Figure S6.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F10)

Supplementary Figure S6. 
Comparative analysis of metabolomic data across newborn dry blood spots from the COPSAC2010 and COPSAC2000 mother-child cohorts. Principal Component Analysis (PCA) plots illustrate the distribution and separation of metabolomic data between cohorts. Panels A - C represent overlapping metabolomes used in our analysis (n=951, n=774, n=626), showcasing subsets of metabolites selected based on varying degrees of correlation and distributional similarity between cohorts. A: displays the PCA for metabolites with correlations greater than 0.4 (n=951), B: Correlations greater than 0.6 (n=774), and C: Metabolites with a correlation greater than 0.6 passing the two-sample Kolmogorov-Smirnov test with a threshold of >=0.05 (n=626).

![Supplementary Figure S7.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F11.medium.gif)

[Supplementary Figure S7.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F11)

Supplementary Figure S7. 
Modulation of the Western Dietary Pattern, Stratified by Child Sex, on Neurodevelopmental Outcomes. Odds ratio estimates of ADHD (**A**) and autism diagnoses (B) based on interactions of the Western dietary pattern, maternal pre-pregnancy BMI (split at median value 23.7), and child’s polygenic risk score (PRS) (median split) for ADHD and autism, stratified by child sex. Estimates for ADHD and autism symptom loads based on the Western dietary pattern, considering tertiles of maternal pre-pregnancy BMI (<22.2, 22.2-25.4, >25.4) and child’s PRS for ADHD and autism, stratified by child sex. Stars represent significance levels: * indicates p < 0.05, ** indicates p < 0.01, and \***| indicates p < 0.001, with "NS" denoting non-significant results (p ≥ 0.05). Further details, including the individual associations of these modulating factors, can be found in Table S11.

![Supplementary Figure S8.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F12.medium.gif)

[Supplementary Figure S8.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F12)

Supplementary Figure S8. 
Modulation of the Western Dietary Pattern Metabolite Score on Neurodevelopmental Outcomes in COPSAC2010 Cohort. Odds ratio estimates for ADHD (A) and autism diagnoses (B) based on interactions of the Western dietary pattern metabolite score, maternal pre-pregnancy BMI (split at median value 23.7), and child’s polygenic risk score (PRS) (median split) for ADHD and autism. Estimates for ADHD (C) and autism (D) symptom loads are based on the Western dietary pattern metabolite score, considering tertiles of maternal pre-pregnancy BMI and child’s PRS for ADHD and autism. Stars represent significance levels: * indicates p < 0.05, ** indicates p < 0.01, and \***| indicates p < 0.001, with "NS" denoting non-significant results (p ≥ 0.05). Further details, including the individual associations of these modulating factors, can be found in Table S11.

![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F4)

Figure 4. 
Modulation of the Western Dietary Pattern on Neurodevelopmental Outcomes in COPSAC2010 Cohort. Odds ratio estimates for ADHD (A) and autism diagnoses (B) based on interactions of the Western dietary pattern, maternal pre-pregnancy BMI (split at median value 23.7), and child’s polygenic risk score (PRS) (median split) for ADHD and autism. Estimates for ADHD (C) and autism (D) symptom loads based on the Western dietary pattern, considering tertiles of maternal pre-pregnancy BMI and child’s PRS for ADHD and autism. Odds ratios and estimates are in relation to a change of 1 SD of the Western dietary pattern. Stars represent significance levels: * indicates p < 0.05, ** indicates p < 0.01, and \***| indicates p < 0.001, with "NS" denoting non-significant results (p ≥ 0.05). Further details, including the individual associations of these modulating factors, can be found in Table S11.

To assess the association of the Western dietary pattern on ADHD and autism symptom loads, we employed a tertile split. Concordantly, we found that both ADHD and autism consistently showed that the association of the Western dietary pattern was most pronounced in children born to mothers in the highest pre-pregnancy BMI group (>25.4) and in children with the highest genetic risk for the respective disorders (*Figure 4C/D*). For ADHD symptom load, this group had an estimate of 5.07 [3.12-7.03], p<10⁻⁶, and for autism symptom load it was 13.5 [9.51-17.5], p<10⁻¹. When further stratified by child sex, this pronounced association persisted significantly only for male children (ADHD: 8.24 [5.19-11.3], p<10⁻⁶, autism: 24.6 [19.0-30.2], p<10⁻⁴ (*Figure S7C/D*). Findings were again consistent when considering the objectively measured WDP-MS (*Figure S8A/B*). A further graphic is provided (*Figure S9*) for ADHD (*A-C*) and autism symptom loads (D-F) showing model estimates when considering maternal pre-pregnancy BMI across the categories <25, 25-30, and >30.

![Supplementary Figure S9.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2024/06/25/2024.03.07.24303907/F13.medium.gif)

[Supplementary Figure S9.](http://medrxiv.org/content/early/2024/06/25/2024.03.07.24303907/F13)

Supplementary Figure S9. 
The Western dietary patterns metabolite scores association with ADHD and autism symptom loads, considering standard clinical classifications of maternal pre-pregnancy BMI. Panels A and D represent ADHD and autism symptom loads respectively, stratified by low/high genetic risk (median cut). Panels B and E further contextualise these associations based on maternal pre-pregnancy BMI categories. Panels C and F again contextualise these associations by stratifying by child sex (male sex shown).

## DISCUSSION

In our large prospective general population COPSAC2010 mother-child cohort, we observed strong associations between a Western dietary pattern during pregnancy with ADHD and autism diagnoses, as well as related symptom loads for these disorders in 10-year-old children. Our findings for ADHD were corroborated in three additional cohorts, which together encompass more than 60,000 mother-child pairs. From the 43 metabolites associated with a Western dietary pattern, a subset of 15 mediating metabolites notably strengthened the association with ADHD diagnosis in external validation. In the COPSAC2010 and VDAART cohorts, we consistently found dietary metabolite score associations in early-mid pregnancy had the strongest associations with neurodevelopmental outcomes, suggesting that this may be a particularly sensitive period of neurodevelopment to dietary influences. Finally, we observed that the association of the Western dietary pattern during pregnancy with ADHD and autism diagnoses, as well as symptom loads, was markedly stronger in children with higher genetic predisposition and maternal pre-pregnancy BMI, particularly in male children.

Expanding on prior research, our study, reinforced by clinically assessed outcomes, presents compelling evidence indicating a potential association between maternal dietary habits during pregnancy and child neurodevelopment. The persistence of the Western dietary pattern’s association on ADHD symptoms, even in children without a neurodevelopmental disorder diagnosis, suggests that the pregnancy Western dietary pattern may also affect the subclinical population. A Norwegian study of 77,768 mother-child pairs found that higher maternal adherence to dietary guidelines was inversely correlated with ADHD diagnosis and symptom severity, while ultra-processed food was positively associated with ADHD symptom severity (26). This study revealed no link between the diet of 3-year-olds and ADHD outcomes, aligning with our results. Conversely, the Raine Study associated a Western dietary pattern in 1,799 adolescents with ADHD but did not explore maternal dietary influences (27). In a cross-cohort study of 800 pregnant women in the US, weak associations were found between established reference dietary patterns and autism-related traits and diagnosis (28). Our study adds to this existing literature by identifying potential temporal associations of a Western dietary pattern on child neurodevelopmental outcomes, using both maternal and child self-reported dietary intakes, and the objective blood metabolome (29).

Our findings suggest a potential interplay between the prenatal environment and child genetics that may influence child neurodevelopment. Males with high genetic risk born from mothers with high pre-pregnancy BMI had the strongest associations with the pregnancy Western dietary pattern. Twin studies support the hypothesis that such interactions play a significant role in the development of these disorders (30,31). Whilst heritability for ADHD and autism is high, with genetic influences accounting for approximately 70-80% of the phenotypic variance (32), current estimates for heritability encompass gene-environment interaction effects (33). Our findings highlight the potential significant role of the environment in the aetiology of ADHD and autism. It is worth noting that much of the heritability for ADHD and autism is explained by SNPs in regulatory regions rather than coding regions (34), providing further evidence for the role of the environment in these disorders. This may aid researchers in understanding the underlying mechanisms contributing to the development of these disorders.

The main strength of our study is the prospectively assessed COPSAC2010 cohort followed from pregnancy through childhood, with a deep psychiatric and dimensional psychopathological characterization provided by the COPSYCH project. We use pregnancy and early-life samples taken before onset of disease and symptoms, which allows for the identification of potential risk factors and the directionality of associations. The breadth of data from the COPSAC2010 cohort enabled the inclusion of relevant covariates, and the depth of phenotyping allowed for genetic characterisation, making assessments of gene-environment interactions possible. By utilising a validated semiquantitative food frequency questionnaire and identifying dietary patterns via data-driven approaches, our findings more accurately reflect the dietary patterns that are existing in the population (35). We validate this methodology by identifying a highly correlated dietary pattern within the DNBC. Although our analysis lacked information on morning sickness and hyperemesis gravidarum, our FFQ was registered at 24 weeks gestation with 1 month recall, outside the period when these symptoms are commonly experienced. Furthermore, universal accessibility of a uniform prenatal care in Denmark enhances the reliability of our findings. The pregnancy Western dietary pattern was derived from estimates of nutrient constituents, this approach is not novel (36), and is supported by a recent meta-analysis across 123 studies (37), which reported the reliably of FFQs to capture nutrient constituents with correlation coefficients consistently exceeding 0.5 for most nutrients. This nutrient-centric methodology may offer accurate representation of dietary patterns across diverse populations, since nutrients align more closely with biological and physiological processes than food group categories.

Acknowledging the limitations of FFQs, including recall bias and quantification challenges, our study incorporates blood metabolomics to provide objective biomarkers for dietary intake. This approach, while not without its own limitations, such as the transient nature of dietary-derived biomarkers, has validated the Western dietary pattern both internally, through a strong metabolomic imprint, and externally, in two independent FFQs in the VDAART cohort. Many studies have likewise found strong dietary signatures in the metabolome, suggesting this may be a promising avenue for personalised nutrition and targeted therapeutic strategies (29).

Our findings for ADHD were strengthened by external validation in three independent mother-child cohorts. This validation enhances the reliability and generalisability of our results, despite the presence of notable differences in cohort characteristics. The validation for ADHD diagnosis in the large-scale DNBC cohort, encompassing 59,725 mother-child pairs, underscores the significant relevance and utility of our findings in the general population. Moreover, validation from the US-based VDAART cohort demonstrates the broad applicability of our findings across diverse ethnic and racial groups, dietary habits, and socio-economic spectra. In VDAART, our external validation was significantly strengthened by focusing on a subset of metabolites identified in causal mediation analysis, implying these metabolites may be crucial intermediaries in the associations observed. Additionally, the VDAART findings enhance our understanding from COPSAC2010, with stronger associations of the Western dietary pattern for ADHD in early pregnancy (10-18 weeks), which may reflect a more sensitive period of neurodevelopment. Given this context, it may not be surprising to see weaker associations from neonatal DBS in COPSAC2000, potentially reflecting dietary intakes in late pregnancy.

We cannot rule out potential confounding by maternal phenotypic factors that relate to child neurodevelopmental disorders, however supplemental analysis adjusting for maternal genetic risk did not change our main findings. Given the observational nature of our data, we cannot elucidate the critical window of these dietary associations. Moreover, it is possible that there are discrepancies in the metabolomics platforms concerning metabolite detection, annotations, and cross-platform compatibility, which may introduce bias. Nevertheless, the same metabolomics platform was employed for all maternal and child timepoints in the COPSAC2010 and VDAART cohorts (HD4, Metabolon, Inc.), as well as for the DBS metabolomics profiles in both the COPSAC2010 and COPSAC2000 cohorts (38). A further limitation is that metabolite scores for children, derived from modelling on pregnant mothers, may not reflect dietary patterns in childhood. However, correlations between pregnancy and child metabolite scores were consistently ≥ 0.33 in both COPSAC2010 and VDAART. Whilst dietary intakes remain stable across pregnancy trimesters (39), our ability to distinguish potential acute effects of a pregnancy Western dietary pattern or longer-term effects is limited by our reliance on a single FFQ assessment.

Moving forward, future research should focus on establishing causality to further elucidate the relationship between Western dietary patterns in pregnancy and neurodevelopmental outcomes. While animal studies have provided valuable insights (40), further clinical investigations are warranted. Our findings underscore the potential for targeted prevention strategies, and conducting randomised controlled trials would be instrumental in establishing causality. Additionally, genetic studies could strengthen causal inference in our observational data (41).

In conclusion, a Western dietary pattern in pregnancy was strongly associated with ADHD and autism, and with symptom loads for these disorders. Our findings suggest that early to mid-pregnancy may be a particularly sensitive window during which dietary influences may impact child neurodevelopment. Highlighting specific subgroup effects, our research points towards the potential for tailored dietary guidelines during pregnancy, especially for at-risk populations. Ultimately, our research underscores the importance of developing targeted dietary interventions for pregnant women to potentially mitigate the risk of neurodevelopmental disorders in children.

## Data Availability

Participant-level personally identifiable data are protected under the Danish Data Protection Act and European Regulation 2016/679 of the European Parliament and of the Council (GDPR) that prohibit distribution even in pseudo-anonymized form. However, participant-level data can be made available under a data transfer agreement as part of a collaboration effort.

[https://copsac.com/](https://copsac.com/) 

## Authors Contributions

DH has written the first draft of the manuscript. All co-authors have provided important intellectual input and contributed considerably to the analyses and interpretation of the data. All authors guarantee that the accuracy and integrity of any part of the work have been appropriately investigated and resolved and all have approved the final version of the manuscript. The corresponding author had full access to the data, with the exception of the DNBC cohort, where analysis was run by TH. The corresponding author had final responsibility for the decision to submit for publication.

## Competing Interests Statement

BE is part of the Advisory Board of Eli Lilly Denmark A/S, Janssen-Cilag, Lundbeck Pharma A/S, and Takeda Pharmaceutical Company Ltd; and has received lecture fees from Bristol-Myers Squibb, Boehringer Ingelheim, Otsuka Pharma Scandinavia AB, Eli Lilly Company, and Lundbeck Pharma A/S. BYG has been the leader of a Lundbeck Foundation Centre of Excellence for Clinical Intervention and Neuropsychiatric Schizophrenia Research (CINS) (January 2009 – December 2021), which was partially financed by an independent grant from the Lundbeck Foundation based on international review and partially financed by the Mental Health Services in the Capital Region of Denmark, the University of Copenhagen, and other foundations. All grants are the property of the Mental Health Services in the Capital Region of Denmark and administered by them. She has no other conflicts to disclose. JL-S is a scientific advisor for Precion Inc and a consultant to Tru Diagnostic, Inc.

Remaining authors declare no potential, perceived, or real conflict of interest regarding the content of this manuscript. The funding agencies did not have any role in design and conduct of the study; collection, management, and interpretation of the data; or preparation, review, or approval of the manuscript. No pharmaceutical company was involved in the study.

## METHODS

### Study Design

The primary analysis of our research was conducted within the Copenhagen Prospective Studies on Asthma in Childhood (COPSAC2010) mother-child cohort, supplemented by validation analyses in three independent mother-child cohorts: the Danish National Birth Cohort (DNBC), the Vitamin D Antenatal Asthma Reduction Trial (VDAART), and COPSAC2000. The COPSAC2010 cohort includes 700 mother-child pairs with extensive phenotyping from 14 clinical visits and exposure assessments since birth, including infants born within a gestational range of 30 to 42 weeks, thus covering both preterm and full-term births (42). Outcome data for this study was derived from the COpenhagen Prospective Study on Neuro-PSYCHiatric Development (COPSYCH) study (43)nested within the COPSAC2010 mother-child cohort, at age 10 years.

In validation analysis, we utilised the DNBC, a large mother-child cohort that recruited over 100,000 pregnant mothers from early gestation, approximately 35% of pregnancies during the recruitment period, to validate our findings of the food frequency questionnaire derived dietary pattern (44,45). Furthermore, we used the large independent US-based VDAART mother-child cohort ([NCT00920621](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT00920621&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom)) (*n=656*) (46), to validate our blood metabolome modelling (47), and further validated our findings in the COPSAC2000 cohort (n=348), using metabolic profiles from neonatal dried blood spots (DBS) (48).

### The COPSYCH 10-year visit

The COPSYCH study comprised a 2-day clinical visit that focused on neurodevelopment, reflected by neurocognition and psychopathology (43). Clinical examinations were carried out between January 2019 and December 2021.

We used both the International Classification of Disorders 10th Revision (ICD-10) of Mental and Behavioural Disorders: Clinical descriptions and diagnostic guidelines and Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria for research diagnoses based on clinical information from the 10-year visit. ADHD ICD-10 diagnostic codes assigned at the COPSYCH visit included: DF90.0 (*Disturbance of activity and attention*), DF90.8 (*Other hyperkinetic disorders*) and DF98.8C (*Attention deficit disorder without hyperactivity*). For secondary analysis, we further specified *ADHD combined presentation* (DSM-5) (fulfilling ICD-10 DF90.0 or DF90.8) and *ADHD predominantly inattentive presentation* (DSM-5) (fulfilling ICD-10 DF98.8*C*). Autism ICD-10 diagnostic codes assigned were: DF84.0 (*Childhood autism*), DF84.5 (*Asperger’s syndrome*) and DF84.8 (*Other pervasive developmental disorders*). For our analysis, we also included additional neurodevelopmental disorder codes under the umbrella of ’any neurodevelopmental disorder’, including DF95.1 (*chronic motor or vocal tic disorder*), DF95.2 (*Tourette’s syndrome*), DF88 (*Other disorders of psychological development*), and DF89 (*Unspecified disorder of psychological development*), in addition to ADHD and autism as previously described.

The ADHD Rating Scale (ADHD-RS) questionnaire was completed by parents to assess the severity of inattentive and impulsive-hyperactive symptoms (6), and the Social Responsiveness Scale-2 (SRS-2) was used to assess the severity of autism traits in the whole cohort (7). Symptom loads are derived from the total scores of the ADHD-RS (Q1-Q18) and SRS-2 questionnaires. Subscales for ADHD-RS and SRS-2 are utilised in subanalysis. In addition, at 6 years, we evaluated hyperactivity/inattention problems using the subscale from the Strengths and Difficulties Questionnaire (SDQ), which scores range from 0 to 10 (49,50).

The Kiddie Schedule for Affective Disorders and Schizophrenia for School-Age Children Present and Lifetime Version (K-SADS-PL) (25) was used to assess current and lifetime psychopathology through a general screening interview and supplementary interview for relevant disorders (51). All K-SADS-PL interviews were conducted by a trained clinician. If certain thresholds are reached, a supplementary interview is conducted to assign relevant diagnoses based on all available information. Consensus diagnoses were made between senior researcher and psychologist with specialty in child and adolescent psychiatry (JRMJ) and at least two examiners according to both the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (52) and the International Classification of Disorders 10th Revision (ICD-10) of Mental and Behavioural Disorders: Clinical descriptions and diagnostic guidelines (53). Diagnostic classification was further validated by a professor in child and adolescent psychiatry (NB). An estimation of interrater reliability on the symptom level was performed based on video recordings of 10 participants with JRJM as the gold standard. The overall agreement, including both currently present and not present symptoms, was 99.48% (95%CI 99.25-99.66). The agreement specifically on present symptoms was 88.48% (95%CI 82.60-92.92).

### Food Frequency Questionnaires and Pregnancy Dietary Patterns

At 24 weeks of gestation, upon their recruitment to the COPSAC2010 cohort, mothers were asked to complete a validated semi-quantitative FFQ with a one-month recall period, encompassing 360 items detailing their dietary intake over the previous month. Mothers in COPSAC2010 and DNBC were subsequently excluded from the analysis for reporting unrealistic energy intakes, either below 4,200 kJ/day or exceeding 16,700 kJ/day. Utilised in both the DNBC validation cohort and the COPSAC2010 cohort, this FFQ has been validated in both a group of younger women (54) and for pregnant women (55,56). The timing of FFQ completion at week 25 in the DNBC cohort, closely aligned with the COPSAC2010 cohort’s at 24 weeks, and thus facilitates the examination of dietary pattern associations on neurodevelopmental outcomes (45). Mothers reported their intake of foods with natural units like apples, while for items without natural units, such as lasagne, they indicated frequency of consumption. Nutrient intakes were estimated using standard portion sizes and recipes, with frequency categories converted to daily intake. Intakes of various food items were calculated based on standard portion size assumptions, measured in grams per day. Nutrient constituents were then determined using the National Food Institute’s Food Composition Databank, version 7 (57). Nutrient constituents were selected for analysis due to its association with physiological and metabolic functions, which provides a perspective on biological mechanisms applicable across various populations. During the COPSYCH clinical assessment at the age of 10, we collected data on the dietary habits of children. The FFQ, encompassing 145 distinct food items, aimed to capture the consistent dietary habits of children. It has been previously validated for use among adolescents aged 12-15 in the DNBC study (58). Intakes of food groups and nutrient components were measured using the same method as the pregnancy questionnaire mentioned above.

In COPSAC2010, principal component analysis (PCA) was used to identify maternal dietary patterns based on the calculated estimates of energy, macronutrient, and micronutrient intake in 613 mothers. In the DNBC cohort, dietary pattern validation was achieved using the principal component model originally trained in COPSAC2010 data, ensuring the identification of directly comparable dietary patterns across the two cohorts. To validate our methodology, we also assessed for associations of a comparable dietary pattern trained within the DNBC. In COPSAC2010, the derived dietary patterns were associated with estimates for food groups for interpretability. The sum totals of nutrient constituents were excluded to avoid redundancy in the PCA.

In the VDAART cohort, mothers were asked to complete a semi-quantitative food frequency questionnaire (FFQ) during the clinical visit corresponding to when the blood metabolome samples were taken: 10-18 weeks and 32-38 weeks gestation (59). The FFQ captured the average weekly consumption of each food item(s) over the past year using the following ordinal categories:"Less than once per week" as 0, "Once per week" as 1, "2-4 times per week" as 2, "Nearly daily or daily" as 3, and "Twice or more per day" as 4.

### DNBC, VDAART and COPSAC2000 Neurodevelopmental Outcomes

In the DNBC cohort, neurodevelopmental outcomes were derived from National registers. Specifically, outcome data was extracted from the Danish National Patient Register (60), the Danish Psychiatric Central Research Register (61), the Danish National Prescription Registry (62). Children were defined as having an ADHD diagnosis if they had received an ADHD diagnosis or redeemed a prescription for ADHD medication after reaching 3 years of age. The diagnoses were obtained from the Danish National Patient Register and Danish Psychiatric Central Research Register. These were recorded according to the International Classification of Diseases Tenth revision (ICD-10) and we defined a case by using the diagnostic codes from F90.x. Prescriptions of ADHD medication were identified from the Danish National Prescription Registry using the Anatomical Therapeutic Chemical Classification codes for methylphenidate (N06BA04), modafinil (N06BA07), and atomoxetine (N06BA09). Similarly, autism diagnosis was defined as children aged 2 years or older with a diagnosis code of ICD-10: F84.0, F84.1, F84.5, F84.8, or F84.9, derived from the Danish National Patient Register, or the Danish Psychiatric Central Research Register. Any neurodevelopmental diagnosis disorder was defined as having an ADHD, or autism diagnosis, as defined above, or any of the following ICD-10 codes: F80-83 (other developmental disorder), F84 (Pervasive developmental disorders), (F88-89 (other/unspecified disorders of psychological development), F95 (tic disorders), F98·8 (ADHD – inattentive type) and F20-F29 (schizophrenia spectrum disorder). In addition, at 7 years, we evaluated hyperactivity/inattention problems using the subscale from the SDQ, which scores range from 0 to 10 (49,50).

Parents of VDAART children at the six-, seven- and eight-year follow-up visits were asked if their child “*had ever received a physician diagnosis of ADD or ADHD*” or “*had ever received a physician diagnosis of autism*”. Our main analysis was on children at the age of 6. Further for an analysis of robustness we assessed associations on whether children had ever had a diagnosis of ADD/ADHD or autism up to the age of eight.

In COPSAC2000, participants were asked to complete the Adult ADHD Self-Report Scale v1.1 (ASRS) at the 18-year clinical cohort visit (63). The ASRS has been validated and demonstrated good reliability and validity in screening for ADHD in adults. In our analysis, participants received a score of 1 for any question that deemed a symptom of ADHD (such that total scores could range from 0 to 18) (64). As a subanalysis we stratified inattention-based questions (Q1-Q4, Q7-Q11) and hyperactivity-impulsivity questions (Q5-Q6, Q12-Q18).

### Blood Metabolome

Untargeted plasma metabolomics data was collected from mothers at mid-pregnancy (24 weeks gestation) and one week postpartum, and from children at ages 6 months, 18 months, and 6 years in the COPSAC2010 cohort. A blood sample was collected using an EDTA tube during the visit to the research clinic site, then centrifuged for 10 minutes at 4000 rpm to extract the plasma. The supernatant was secured and stored at −80 °C for future analysis. The untargeted metabolomic analysis of the plasma samples for both COPSAC2010 and the VDAART cohorts were carried out by Metabolon, Inc. (NC, USA).

For sample preparation, we employed the MicroLab STAR® system from Hamilton Company, which is an automated process. Each sample was enhanced with recovery standards to ensure quality control (QC) before metabolites were extracted using methanol. This extraction process involved vigorous shaking for 2 minutes using a Glen Mills GenoGrinder 2000, followed by centrifugation, to precipitate proteins. The extract was then divided into four aliquots for further analysis on four different LC-MS/MS platforms. The aliquots were dried using a TurboVap® (Zymark) to evaporate the organic solvent, then stored under nitrogen overnight prior to LC-MS/MS preparation.

The analysis through LC-MS/MS was conducted using an ACQUITY Ultra-Performance Liquid Chromatography (UPLC) system by Waters, Milford, USA, coupled with a Q Exactive™ Hybrid Quadrupole-Orbitrap™ mass spectrometer, which includes a heated electrospray ionisation (HESI-II) source, provided by ThermoFisher Scientific, Waltham, Massachusetts, USA. We prepared the sample extracts in solvent mixes specifically chosen for each of the four LC methods applied: two methods were reverse phase UPLC-ESI(+) MS/MS for analysing both hydrophilic and hydrophobic molecules; another was a reverse phase UPLC-(−) MS/MS; and the last was HILIC/UPLC-(−) MS/MS. The mass spectrometry analysis alternated between full scan MS and data-dependent MSn scans with dynamic exclusion, covering a scan range from 70 to 1000 m/z for both ion modes.

In terms of data collection and quality control, the raw data was subjected to extraction, peak identification and followed by QC procedures. The semi-quantification of samples was based on the area-under-the-curve method. Further details are discussed in our previously published work (65).

Data preprocessing was done by excluding metabolites with >33% missingness. Remaining missing data was imputed using random forest imputation (missForest package, v1.5) (66) and the metabolome data was log-transformed, centred and scaled, prior to analysis. A pattern of missingness identified in the metabolite subset (0.11 [0.02-0.20], p=0.021) suggested under detection in the Western dietary pattern, with no corresponding bias for neurodevelopmental outcomes or across the full metabolome (*Table S13*), reinforcing the validity of our imputation approach. The metabolome data for mothers at mid-pregnancy (24 weeks gestation) had a total of 760 annotated metabolites for analysis. Of these, 744 metabolites overlapped between the 24-week gestation period and one week postpartum, 516 metabolites overlapped with 6-month data, 707 metabolites overlapped with 18-month data and 540 metabolites overlapped with 6-year data.

We used two VDAART blood metabolome pregnancy time points (10-18 weeks and 32-38 weeks), measured on the same platform as the COPSAC2010 samples, to replicate our findings for ADHD. The 24-week pregnancy time point in COPSAC2010 was compared with those for early pregnancy (10-18 weeks gestation) and late pregnancy (32-38 weeks gestation) in VDAART, identifying 640 and 689 overlapping metabolites, respectively. We further predicted metabolite scores for the VDAART children, with 523 overlapping metabolites at 1 year, 530 at 3 years, and 677 overlapping metabolites at 6 years. The metabolite scores were predicted using the model trained on COPSAC2010 data; applied to log-transformed, centred and scaled VDAART data. The metabolite score predictions were subsequently scaled, and thus results should be interpreted as per SD change in the given population.

### Neonatal Dried Blood Spot samples

DBS samples from COPSAC2010 (*n=677*) and COPSAC2000 (*n=387*) were respectively collected at age 2-3 days and 1-12 days after birth, and stored at -20°C at the Danish National Biobank until analysis. Metabolic profiles of the DBS samples were acquired by liquid chromatography mass spectrometry (LC-MS). Data preprocessing was executed using both XCMS and MZmine, with the quality assessed based on the distribution of pooled samples (38). Metabolite extraction from DBS samples, a 3.2mm punch in diameter, was conducted using a Microlab STAR automated liquid handler from Hamilton Bonaduz AG, Bonaduz, Switzerland. These samples were extracted onto 96-well plates in batches with 100 μL of 80% methanol, specifically designed for analytical processes. Subsequently, 75 μL of the supernatants were shifted to fresh plates, air-dried under a nitrogen atmosphere, and reconstituted in 75 μL of 2.5% methanol. Before injection, 65 μL was transferred onto the final plates. For the analysis, samples from both COPSAC2000 and COPSAC2010 were divided, with the former distributed across six batches and the latter across ten. Each batch or plate comprised eight water blanks, an internal standard mix, four external controls, three paper blanks, four pooled samples, two diluted pools, and 74 cohort samples. All extractions utilised LC-MS grade solvents provided by Thermo Fisher Scientific, Waltham, MA, USA. An addition of 24 isotope-labelled internal standards, encompassing amino acids and acylcarnitines, was made to the extraction solvents. These standards, consistent in concentration across samples, were employed to ensure quality and evaluate the efficacy of batch correction and normalisation that followed. To further gauge the variability in sample preparation, acquisition, and preprocessing, quality control pooled samples (PSs) were periodically assessed throughout the analysis run. Metabolic profiles were generated using a high-resolution Thermo Scientific Q-Exactive Orbitrap mass spectrometer, linked to a Dionex Ultimate 3000 UPLC. This complex procedure involved a series of specific injection, gradient operation, and parameters for mass spectrometric analysis. Preprocessing of the acquired data involved the use of software like XCMS, MZmine, and MSconvert, with specific protocols established for peak detection, alignment, deconvolution, and filtering. Ensuring data quality was paramount, and this was maintained through the use of internal standards and a series of checks, including principal component analysis (PCA) to confirm the clustering of pooled samples. For a more detailed insight into the equipment settings, data preprocessing steps, and quality controls, refer to the specific sections on mass spectrometry settings, LC-MS/MS data preprocessing, and quality control in the supplementary information. Further in-depth description of sample preparation, LC-MS metabolic profiling, and data pre-processing for both cohorts can be found in previous work (38).There were 2313 features detected in COPSAC2010 and 2363 features in COPSAC2000. As the majority of features were unnamed, we merged both datasets using an inexact matching criterion (mass ±0.01, retention time ± 0.1) using the fuzzyjoin R package (v0.1.6) (67) and subsequently found 1253 overlapping metabolites between cohorts. For further robustness, we applied a further filter, only retaining features that had a cross-correlation between cohorts of >0.4, ending with 951 overlapping features between cohorts.

### Information on Covariates

In our multivariable analysis, we included the following covariates for COPSAC2010: pre-pregnancy maternal BMI, social circumstances (the first principal component of household income, maternal education level, and maternal age at birth), child sex, birth weight, gestational age, smoking during pregnancy, antibiotic use during pregnancy, pre-eclampsia and a child Western dietary pattern, assessed by a FFQ from children at 10 years old. Missing covariate data was imputed, based on the available covariate information available using the imputePCA function with one component from the missMDA (v1.18) R package.

Additionally, we conducted sensitivity analysis in COPSAC2010 to explore the influence of the following variables: the results of two randomised controlled trials conducted within the COPSAC2010 cohort (n3-LCPUFA supplementation and high-dose vitamin D during pregnancy) (42), total energy intake derived from the maternal FFQ, maternal ADHD genetic risk, intake of prenatal supplements, and genomic PCs. Specifically, adjustment of intake of prenatal supplements was conducted by adjusting for the first seven principal components (explaining 78.6% of the total variance) of the estimated intakes of 31 vitamins, minerals and fatty acids, recorded at the time of the 24 week FFQ, and whether prenatal supplements were taken (yes (94.1%) or no (5.9%)). Missing data were imputed using the same approach described above from the missMDA R package.

Covariates used in multivariable analysis for the DNBC included pre-pregnancy maternal BMI, parental social group (high, medium, skilled, student, unskilled, unemployed), maternal age, child sex, birth weight, gestational age and smoking during pregnancy. Missing covariate data were imputed based on existing covariates using multiple imputation as implemented by the PROC MI procedure in SAS (version 9.1, SAS Institute, Cary, NC).

Covariates used in multivariable analysis for the VDAART included pre-pregnancy maternal BMI, social circumstances (the first principal component of household income, maternal education level, and maternal age at birth), child sex, birth weight, gestational age, smoking during pregnancy, pre-eclampsia, and race. Missing covariate data was imputed, based on the available covariate information using the imputePCA function with one component from the missMDA (v1.18) R package.

Covariates used in multivariable analysis for COPSAC2000 included social circumstances (the first principal component of household income, maternal education level, and maternal age at birth), child sex, birth weight, antibiotics during pregnancy and smoking during pregnancy.

The child sex variable in COPSAC2010, the DNBC, COPSAC2000 and VDAART was interpreted as sex assigned at birth.

### Polygenic Risk Scores

Maternal and child genotypes in the COPSAC2010 cohort were analysed using the Illumina Infinium HumanOmniExpressExome BeadChip. To calculate PRS for ADHD and autism, we used data from genome-wide association meta-analyses of diagnosed ADHD cases (*n=20,183*) and controls (*n=35,191*) (68) and autism cases (*n=18,381*) and controls (*n=27,969)* (69).

The software package Polygenic risk score continuous shrinkage (CS) was used for PRS construction through regularising SNP effects using a shrinkage prior (70). PRS-CS uses a linkage disequilibrium (LD) reference panel constructed via European ancestry samples from the 1000 genomes project. We used the automatic estimation of the polygenicity parameter (phi). Subsequent to SNP effect adjustment with PRS-CS, we used the software PLINK2 (71) to aggregate SNP effects into PRSs for each individual. The scores were scaled to mean 0 and SD 1 for each phenotype.

### Statistical Analysis

Logistic regression models were used to determine the associations of maternal dietary patterns on neurodevelopmental disorders, and linear regression models were used to assess the associations of maternal dietary patterns on parentally reported symptom loads for ADHD and autism. All multivariable models were adjusted for the previously mentioned covariates. Dietary principal components and metabolite scores were scaled to enhanced interpretation (Odds ratios and estimates are interpreted as per SD change). We used linear interaction models to examine the modulating effects of genetics, child sex and maternal BMI on neurodevelopmental outcomes. All interaction models were multivariable (see ’Information on covariates’), additionally including genetic risk (see ’Polygenic Risk Scores’). We excluded five individuals from the analysis. Specifically, we removed the second twin from twin pairs since the exposure of interest was not independent. However, non-twin sibling pairs were retained in the analysis, as the exposure and outcome measures were independently assessed for each sibling.

We used Gaussian graphical network models via the framework described by Williams and Rast (72). Gaussian graphical models show the non-zero relationships (95% Cl) between the dietary pattern, neurodevelopmental outcomes and model covariates, controlling for the linear effects of all covariates expressed as partial correlations. We integrated gestational age and birthweight to refine our Gaussian graphical network models, as these variables are highly collinear (73). From these models, significant associations in the resulting precision matrices were visualised as interconnected network diagrams utilising the ggraph (v2.1.0), igraph (v1.2.11), and qgraph (v1.9.2) libraries in R.

We used the caret package (v6.0.90) to establish Western dietary pattern metabolite scores (WDP-MS), via sparse partial least squares regression on the COPSAC2010 metabolomics datasets, with the pregnancy Western dietary pattern as the response. Individual models were created for predicting WDP-MS at other timepoints, using the subset of overlapping metabolites (See ’Blood Metabolomè). To improve interpretability and reduce the risk of overfitting, we used single-component models, applying cross-validated predictions (repeated cv, number of segments=5, repeats=10). After reviewing models that varied in sparsity (incremented by 0.1 from 0 to 1), we opted for the model with the lowest root-mean-standard error for cross-validation (RMSECV) as our selection criteria. For external validation between COPSAC2010 and VDAART, we adjusted our model parameters to balance against overfitting whilst ensuring optimal performance, leading to the selection of models that were within +0.01 RMSECV of the best-performing model. We similarly trained a sPLS model in the COPSAC2010 cohort on the Western dietary pattern from DBS metabolic profiles that overlapped with the COPSAC2000 cohort.

To elucidate the role of metabolites as potential biomarkers or mediators in the association between a Western dietary pattern and neurodevelopmental disorders, we employed a systematic backward elimination strategy. This approach was designed to iteratively exclude metabolites, aiming to identify those with the most profound mediating influence on the outcome. We employed 10,000 iterations for each step of this elimination. We used two linear models: one linking a composite metabolite score to the Western dietary pattern and covariates, and another assessing its mediating role between dietary patterns and neurodevelopmental outcome. This causal mediation analysis was conducted with the mediation package in R (v4.5.0).

When comparing the associations of the mother and child WDP-MS in logistic regression, we utilised the imputePCA function from the missMDA package (v1.18) (74) to impute missing data with one component (12.5% in COPSAC2010). This was done using the available metabolite scores as a reference. Linear mixed models with variable slopes and intercepts were applied to the WDP-MS data, integrating gestational age alongside chronological age. The slopes and intercepts were treated as random effects, allowing for modelling of the longitudinal association of the identified maternal dietary pattern on neurodevelopmental outcomes. Accordingly, the random effects were extracted for inference, the "lme4" R-package (1.1.28) was used for modelling (75), both slopes and intercepts were scaled (estimates are interpreted as per SD change).

A significance level of 0.05 was used in all analyses and false-discovery rate control (FDR) applied where relevant (<0.05). All data analyses were performed with the statistical software R version 4.1.1. Other R packages utilised in this analysis include tidyverse (v1.3.1), dplyr (v1.0.10), broom (v0.7.12), lubridate (v1.8.0), ggpubr (v0.4.0) and tableone (v0.13.0).

## Data Availability

Participant-level personally identifiable data are protected under the Danish Data Protection Act and European Regulation 2016/679 of the European Parliament and of the Council (GDPR) that prohibit distribution even in pseudo-anonymized form. However, participant-level data can be made available under a data transfer agreement as part of a collaboration effort.

## Code Availability

The custom code employed in this research is freely accessible to the public for transparency and reproducibility purposes.

## Acknowledgements

We express our deepest gratitude to the children and families of the COPSAC2010, COPSAC2000 and VDAART cohorts for all their support and commitment. We acknowledge and appreciate the unique efforts of the COPSAC research team. We acknowledge all funding received by COPSAC, listed on [www.copsac.com](http://www.copsac.com). The Lundbeck Foundation (Grant no R16-A1694 and R269-2017-5); The Ministry of Health (Grant no 903516); Danish Council for Strategic Research (Grant no 0603-00280B) and The Capital Region Research Foundation have provided core support to the COPSAC research centre. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 946228) (BC). MAR is funded by the Novo Nordisk Foundation (Grant no NNF21OC0068517). CP was supported by an Australian National Health and Medical Research Council (NHMRC) L3 Investigator Grant (1196508) and by a grant from the Lundbeck Foundation (ID: R246-2016-3237). JL-S (R01HL123915, R01HL155742, and R01HL141826) and SHC (K01HL153941) are funded through the National Institute of Heart, Lung, and Blood Institute.

The Danish National Birth Cohort was established with a significant grant from the Danish National Research Foundation. Additional support was obtained from the Danish Regional Committees, the Pharmacy Foundation, the Egmont Foundation, the March of Dimes Birth Defects Foundation, the Health Foundation and other minor grants. The DNBC Biobank has been supported by the Novo Nordisk Foundation and the Lundbeck Foundation. Follow-up of mothers and children have been supported by the Danish Medical Research Council (SSVF 0646, 271-08-0839/06-066023, O602-01042B, 0602-02738B), the Lundbeck Foundation (195/04, R100-A9193), The Innovation Fund Denmark 0603-00294B (09–067124), the Nordea Foundation (02-2013-2014), Aarhus Ideas (AU R9-A959-13-S804), University of Copenhagen Strategic Grant (IFSV 2012), and the Danish Council for Independent Research (DFF – 4183-00594 and DFF - 4183-00152).

## Footnotes

*   We have added a new validation mother-child cohort of 59725 (Danish National Birth Cohort) to validate our association between an FFQ-derived pregnancy Western dietary pattern and ADHD. Furthermore we have changed the term "Unhealthy" dietary pattern to "Western" dietary pattern throughout the manuscript, and title.

*   Received March 7, 2024.
*   Revision received June 25, 2024.
*   Accepted June 25, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References (Main Text)

1.  1.Bougeard C, Picarel-Blanchot F, Schmid R, Campbell R, Buitelaar J. Prevalence of Autism Spectrum Disorder and Co-morbidities in Children and Adolescents: A Systematic Literature Review. Front Psychiatry. 2021 Oct 27;12:744709.
    
    

2.  2.Dalsgaard S, Thorsteinsson E, Trabjerg BB, Schullehner J, Plana-Ripoll O, Brikell I, et al. Incidence Rates and Cumulative Incidences of the Full Spectrum of Diagnosed Mental Disorders in Childhood and Adolescence. JAMA Psychiatry. 2020 Feb 1;77(2):155–64.
    
    

3.  3.Willcutt EG. The prevalence of DSM-IV attention-deficit/hyperactivity disorder: a meta-analytic review. Neurotherapeutics. 2012 Jul;9(3):490–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s13311-012-0135-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22976615&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308826800002&link_type=ISI) 

4.  4.Talantseva OI, Romanova RS, Shurdova EM, Dolgorukova TA, Sologub PS, Titova OS, et al. The global prevalence of autism spectrum disorder: A three-level meta-analysis. Front Psychiatry. 2023 Feb 9;14:1071181.
    
    

5.  5.Larsson H, Anckarsater H, Råstam M, Chang Z, Lichtenstein P. Childhood attention-deficit hyperactivity disorder as an extreme of a continuous trait: a quantitative genetic study of 8,500 twin pairs. J Child Psychol Psychiatry. 2012 Jan;53(1):73–80.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-7610.2011.02467.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21923806&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298003300010&link_type=ISI) 

6.  6.Makransky G, Bilenberg N. Psychometric properties of the parent and teacher ADHD Rating Scale (ADHD-RS): measurement invariance across gender, age, and informant. Assessment. 2014 Dec;21(6):694–705.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1073191114535242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24852496&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

7.  7.Constantino JN, Davis SA, Todd RD, Schindler MK, Gross MM, Brophy SL, et al. Validation of a brief quantitative measure of autistic traits: comparison of the social responsiveness scale with the autism diagnostic interview-revised. J Autism Dev Disord. 2003 Aug;33(4):427–33.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1025014929212&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12959421&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184587800007&link_type=ISI) 

8.  8.Cortés-Albornoz MC, García-Guáqueta DP, Velez-van-Meerbeke A, Talero-Gutiérrez C. Maternal Nutrition and Neurodevelopment: A Scoping Review. Nutrients [Internet]. 2021 Oct 8;13(10). Available from: doi:10.3390/nu13103530
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu13103530&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34684531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

9.  9.Fernandes DJ, Spring S, Roy AR, Qiu LR, Yee Y, Nieman BJ, et al. Exposure to maternal high-fat diet induces extensive changes in the brain of adult offspring. Transl Psychiatry. 2021 Mar 2;11(1):149.
    
    

10. 10.Giriko CÁ, Andreoli CA, Mennitti LV, Hosoume LF, Souto TDS, Silva AV da, et al. Delayed physical and neurobehavioral development and increased aggressive and depression-like behaviors in the rat offspring of dams fed a high-fat diet. Int J Dev Neurosci. 2013 Dec;31(8):731–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijdevneu.2013.09.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24071008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

11. 11.Mitchell AJ, Khambadkone SG, Dunn G, Bagley J, Tamashiro KLK, Fair D, et al. Maternal Western-style diet reduces social engagement and increases idiosyncratic behavior in Japanese macaque offspring. Brain Behav Immun [Internet]. 2022 Jul 6; Available from: doi:10.1016/j.bbi.2022.07.004
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2022.07.004&link_type=DOI) 

12. 12.Prado E, Dewey K. Nutrition and Brain Development in Early Life [Internet]. Prenatal and Childhood Nutrition. 2015. p. 79–126. Available from: doi:10.1201/b18040-6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1201/b18040-6&link_type=DOI) 

13. 13.Cordain L, Eaton SB, Sebastian A, Mann N, Lindeberg S, Watkins BA, et al. Origins and evolution of the Western diet: health implications for the 21st century. Am J Clin Nutr. 2005 Feb;81(2):341–54.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI4MS8yLzM0MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA2LzI1LzIwMjQuMDMuMDcuMjQzMDM5MDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Devarshi PP, Grant RW, Ikonte CJ, Hazels Mitmesser S. Maternal Omega-3 Nutrition, Placental Transfer and Fetal Brain Development in Gestational Diabetes and Preeclampsia. Nutrients [Internet]. 2019 May 18;11(5). Available from: doi:10.3390/nu11051107
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu11051107&link_type=DOI) 

15. 15.Taylor MR, Chuang C, Carrasco KD, Nagatomo S, Rucklidge JJ. Dietary and Micronutrient Treatments for Children with Neurodevelopment Disorders [Internet]. Vol. 5, Current Developmental Disorders Reports. 2018. p. 243–52. Available from: doi:10.1007/s40474-018-0150-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40474-018-0150-5&link_type=DOI) 

16. 16.Halliwell B, Cheah IK, Tang RMY. Ergothioneine - a diet-derived antioxidant with therapeutic potential. FEBS Lett. 2018 Oct;592(20):3357–66.
    
    

17. 17.Prince N, Chu SH, Chen Y, Mendez KM, Hanson E, Green-Snyder L, et al. Phenotypically driven subgroups of ASD display distinct metabolomic profiles. Brain Behav Immun. 2023 Jul;111:21–9.
    
    

18. 18. Poornimai Abirami GP, Radhakrishnan RK, Johnson E, Roshan SA, Yesudhas A, Parveen S, et al. The Regulation of Reactive Neuroblastosis, Neuroplasticity, and Nutraceuticals for Effective Management of Autism Spectrum Disorder. Adv Neurobiol. 2020;24:207–22.
    
    

19. 19.Naveed S, Lakka T, Haapala EA. An Overview on the Associations between Health Behaviors and Brain Health in Children and Adolescents with Special Reference to Diet Quality [Internet]. Vol. 17, International Journal of Environmental Research and Public Health. 2020. p. 953. Available from: doi:10.3390/ijerph17030953
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijerph17030953&link_type=DOI) 

20. 20.Basak S, Duttaroy AK. Maternal PUFAs, Placental Epigenetics, and Their Relevance to Fetal Growth and Brain Development. Reprod Sci [Internet]. 2022 Jun 8; Available from: doi:10.1007/s43032-022-00989-w
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s43032-022-00989-w&link_type=DOI) 

21. 21.Pettersson E, Lichtenstein P, Larsson H, Song J, Attention Deficit/Hyperactivity Disorder Working Group of the iPSYCH-Broad-PGC Consortium, Autism Spectrum Disorder Working Group of the iPSYCH-Broad-PGC Consortium, Bipolar Disorder Working Group of the PGC, Eating Disorder Working Group of the PGC, Major Depressive Disorder Working Group of the PGC, Obsessive Compulsive Disorders and Tourette Syndrome Working Group of the PGC, Schizophrenia CLOZUK, Substance Use Disorder Working Group of the PGC, Agrawal A, et al. Genetic influences on eight psychiatric disorders based on family data of 4 408 646 full and half-siblings, and genetic data of 333 748 cases and controls - CORRIGENDUM. Psychol Med. 2019 Jan;49(2):351.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291718002945&link_type=DOI) 

22. 22.Atladottir HO, Gyllenberg D, Langridge A, Sandin S, Hansen SN, Leonard H, et al. The increasing prevalence of reported diagnoses of childhood psychiatric disorders: a descriptive multinational comparison. Eur Child Adolesc Psychiatry. 2015 Feb;24(2):173–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00787-014-0553-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24796725&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

23. 23.Tran NQV, Miyake K. Neurodevelopmental Disorders and Environmental Toxicants: Epigenetics as an Underlying Mechanism. Int J Genomics Proteomics [Internet]. 2017 [cited 2023 Jan 7];2017. Available from: [https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5439185/](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5439185/)
    
    

24. 24.Keers R, Coleman JRI, Lester KJ, Roberts S, Breen G, Thastum M, et al. A Genome-Wide Test of the Differential Susceptibility Hypothesis Reveals a Genetic Predictor of Differential Response to Psychological Treatments for Child Anxiety Disorders. Psychother Psychosom. 2016 Apr 5;85(3):146–58.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

25. 25.Faraone SV, Banaschewski T, Coghill D, Zheng Y, Biederman J, Bellgrove MA, et al. The World Federation of ADHD International Consensus Statement: 208 Evidence-based conclusions about the disorder. Neurosci Biobehav Rev. 2021 Sep;128:789–818.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neubiorev.2021.01.022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

26. 26.Borge TC, Biele G, Papadopoulou E, Andersen LF, Jacka F, Eggesbø M, et al. The associations between maternal and child diet quality and child ADHD - findings from a large Norwegian pregnancy cohort study. BMC Psychiatry. 2021 Mar 8;21(1):139.
    
    

27. 27.Howard AL, Robinson M, Smith GJ, Ambrosini GL, Piek JP, Oddy WH. ADHD is associated with a “Western” dietary pattern in adolescents. J Atten Disord. 2011 Jul;15(5):403–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1087054710365990&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20631199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

28. 28.Vecchione R, Wang S, Rando J, Chavarro JE, Croen LA, Fallin MD, et al. Maternal Dietary Patterns during Pregnancy and Child Autism-Related Traits: Results from Two US Cohorts. Nutrients [Internet]. 2022 Jun 30;14(13). Available from: doi:10.3390/nu14132729
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu14132729&link_type=DOI) 

29. 29.Noerman S, Landberg R. Blood metabolite profiles linking dietary patterns with health-Toward precision nutrition. J Intern Med. 2023 Apr;293(4):408–32.
    
    

30. 30.Faraone SV, Larsson H. Genetics of attention deficit hyperactivity disorder. Mol Psychiatry. 2019 Apr;24(4):562–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41380-018-0070-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

31. 31.Santos JX, Rasga C, Marques AR, Martiniano H, Asif M, Vilela J, et al. A Role for Gene-Environment Interactions in Autism Spectrum Disorder Is Supported by Variants in Genes Regulating the Effects of Exposure to Xenobiotics. Front Neurosci. 2022 May 19;16:862315.
    
    

32. 32.Ma SL, Chen LH, Lee CC, Lai KYC, Hung SF, Tang CP, et al. Genetic Overlap Between Attention Deficit/Hyperactivity Disorder and Autism Spectrum Disorder in Gene. Front Neurosci. 2021 Apr 27;15:649588.
    
    

33. 33.Bateson P, Gluckman P. Plasticity, Robustness, Development and Evolution. Cambridge University Press; 2011. 166 p.
    
    

34. 34.Gusev A, Lee SH, Trynka G, Finucane H, Vilhjálmsson BJ, Xu H, et al. Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. Am J Hum Genet. 2014 Nov 6;95(5):535–52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2014.10.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25439723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

35. 35.Newby PK, Tucker KL. Empirically derived eating patterns using factor or cluster analysis: a review. Nutr Rev. 2004 May;62(5):177–203.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1301/nr.2004.may.177-203&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15212319&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227434100001&link_type=ISI) 

36. 36.Zhao J, Li Z, Gao Q, Zhao H, Chen S, Huang L, et al. A review of statistical methods for dietary pattern analysis. Nutr J. 2021 Apr 19;20(1):37.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12937-021-00692-7&link_type=DOI) 

37. 37.Cui Q, Xia Y, Wu Q, Chang Q, Niu K, Zhao Y. A meta-analysis of the reproducibility of food frequency questionnaires in nutritional epidemiological studies. Int J Behav Nutr Phys Act. 2021 Jan 11;18(1):12.
    
    

38. 38.Gürdeniz G, Ernst M, Rago D, Kim M, Courraud J, Stokholm J, et al. Neonatal metabolome of caesarean section and risk of childhood asthma. Eur Respir J [Internet]. 2022 Jun;59(6). Available from: doi:10.1183/13993003.02406-2021
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjEyOiI1OS82LzIxMDI0MDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNi8yNS8yMDI0LjAzLjA3LjI0MzAzOTA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

39. 39.Savard C, Lemieux S, Weisnagel S, Fontaine-Bisson B, Gagnon C, Robitaille J, et al. Trimester-Specific Dietary Intakes in a Sample of French-Canadian Pregnant Women in Comparison with National Nutritional Guidelines [Internet]. Vol. 10, Nutrients. 2018. p. 768. Available from: doi:10.3390/nu10060768
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu10060768&link_type=DOI) 

40. 40.Urbonaite G, Knyzeliene A, Bunn FS, Smalskys A, Neniskyte U. The impact of maternal high-fat diet on offspring neurodevelopment [Internet]. Vol. 16, Frontiers in Neuroscience. 2022. Available from: doi:10.3389/fnins.2022.909762
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2022.909762&link_type=DOI) 

41. 41.Pingault JB, O’Reilly PF, Schoeler T, Ploubidis GB, Rijsdijk F, Dudbridge F. Using genetic data to strengthen causal inference in observational research. Nat Rev Genet. 2018 Sep;19(9):566–80.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

## Methods-only References

1.  42.Bisgaard H, Vissing NH, Carson CG, Bischoff AL, Følsgaard NV, Kreiner-Møller E, et al. Deep phenotyping of the unselected COPSAC 2010 birth cohort study [Internet]. Vol. 43, Clinical & Experimental Allergy. 2013. p. 1384–94. Available from: doi:10.1111/cea.12213
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cea.12213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24118234&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

2.  43.Mohammadzadeh P, Rosenberg JB, Vinding R, Møllegaard Jepsen JR, Lindberg U, Følsgaard N, et al. Effects of prenatal nutrient supplementation and early life exposures on neurodevelopment at age 10: a randomised controlled trial - the COPSYCH study protocol. BMJ Open. 2022 Feb 1;12(2):e047706.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiMTIvMi9lMDQ3NzA2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDYvMjUvMjAyNC4wMy4wNy4yNDMwMzkwNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  44.Olsen J, Melbye M, Olsen SF, Sørensen TI, Aaby P, Andersen AM, et al. The Danish National Birth Cohort--its background, structure and aim. Scand J Public Health. 2001 Dec;29(4):300–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/140349401317115268&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11775787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172566800011&link_type=ISI) 

4.  45.Olsen SF, Mikkelsen TB, Knudsen VK, Orozova-Bekkevold I, Halldórsson TI, Strøm M, et al. Data collected on maternal dietary exposures in the Danish National Birth Cohort. Paediatr Perinat Epidemiol. 2007 Jan;21(1):76–86.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-3016.2007.00777.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17239183&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244524300010&link_type=ISI) 

5.  46.Litonjua AA, Lange NE, Carey VJ, Brown S, Laranjo N, Harshfield BJ, et al. The Vitamin D Antenatal Asthma Reduction Trial (VDAART): rationale, design, and methods of a randomized, controlled trial of vitamin D supplementation in pregnancy for the primary prevention of asthma and allergies in children. Contemp Clin Trials. 2014 May;38(1):37–50.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cct.2014.02.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24614387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000337203900005&link_type=ISI) 

6.  47.Kelly RS, Boulin A, Laranjo N, Lee-Sarwar K, Chu SH, Yadama AP, et al. Metabolomics and Communication Skills Development in Children; Evidence from the Ages and Stages Questionnaire. Metabolites [Internet]. 2019 Mar 5;9(3). Available from: doi:10.3390/metabo9030042
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/metabo9030042&link_type=DOI) 

7.  48.Bisgaard H. The Copenhagen Prospective Study on Asthma in Childhood (COPSAC): design, rationale, and baseline data from a longitudinal birth cohort study. Ann Allergy Asthma Immunol. 2004 Oct;93(4):381–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1081-1206(10)61398-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15521375&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224655700013&link_type=ISI) 

8.  49.Obel C, Heiervang E, Rodriguez A, Heyerdahl S, Smedje H, Sourander A, et al. The Strengths and Difficulties Questionnaire in the Nordic countries. Eur Child Adolesc Psychiatry. 2004;13 Suppl 2:II32–9.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15243784&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

9.  50.Sveen TH, Berg-Nielsen TS, Lydersen S, Wichstrøm L. Detecting psychiatric disorders in preschoolers: screening with the strengths and difficulties questionnaire. J Am Acad Child Adolesc Psychiatry. 2013 Jul;52(7):728–36.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaac.2013.04.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23800486&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

10. 51.Kaufman J, Birmaher B, Brent D, Rao U, Flynn C, Moreci P, et al. Schedule for Affective Disorders and Schizophrenia for School-Age Children-Present and Lifetime Version (K-SADS-PL): initial reliability and validity data. J Am Acad Child Adolesc Psychiatry. 1997 Jul;36(7):980–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00004583-199707000-00021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9204677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XF81700021&link_type=ISI) 

11. 52.American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders: DSM-5-TR. 2022. 1120 p.
    
    

12. 53.World Health Organization. The International Statistical Classification of Diseases and Health Related Problems ICD-10: Tenth Revision. Volume 2: Instruction Manual. World Health Organization; 2004. 198 p.
    
    

13. 54.Friis S, Kruger Kjaer S, Stripp C, Overvad K. Reproducibility and relative validity of a self-administered semiquantitative food frequency questionnaire applied to younger women. J Clin Epidemiol. 1997 Mar;50(3):303–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0895-4356(96)00379-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9120530&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WQ59200012&link_type=ISI) 

14. 55.Mikkelsen TB, Olsen SF, Rasmussen SE, Osler M. Relative validity of fruit and vegetable intake estimated by the food frequency questionnaire used in the Danish National Birth Cohort. Scand J Public Health. 2007;35(2):172–9.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17454921&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

15. 56.Mikkelsen TB, Osler M, Olsen SF. Validity of protein, retinol, folic acid and n-3 fatty acid intakes estimated from the food-frequency questionnaire used in the Danish National Birth Cohort. Public Health Nutr. 2006 Sep;9(6):771–8.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16925883&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239972500015&link_type=ISI) 

16. 57.Paperpile [Internet]. [cited 2023 Oct 26]. Danish Food Composition Data, Version 7. Available from: [https://paperpile.com/app/p/49f46071-353d-0183-9823-b8b1c9c5898d](https://paperpile.com/app/p/49f46071-353d-0183-9823-b8b1c9c5898d)
    
    

17. 58.Bjerregaard AA, Halldorsson TI, Kampmann FB, Olsen SF, Tetens I. Relative validity of a web-based food frequency questionnaire for Danish adolescents. Nutr J. 2018 Jan 12;17(1):9.
    
    

18. 59.Kachroo P, Kelly RS, Mirzakhani H, Lee-Sarwar K, Chawes BL, Blighe K, et al. Fish oil supplementation during pregnancy is protective against asthma/wheeze in offspring. J Allergy Clin Immunol Pract. 2020 Jan;8(1):388–91.e2.
    
    

19. 60.Schmidt M, Schmidt SAJ, Sandegaard JL, Ehrenstein V, Pedersen L, Sørensen HT. The Danish National Patient Registry: a review of content, data quality, and research potential. Clin Epidemiol. 2015 Nov 17;7:449–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2147/CLEP.S91125&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26604824&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

20. 61.Mors O, Perto GP, Mortensen PB. The Danish Psychiatric Central Research Register. Scand J Public Health. 2011 Jul;39(7 Suppl):54–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1403494810395825&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21775352&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292984700013&link_type=ISI) 

21. 62.Pottegård A, Schmidt SAJ, Wallach-Kildemoes H, Sørensen HT, Hallas J, Schmidt M. Data Resource Profile: The Danish National Prescription Registry. Int J Epidemiol. 2017 Jun 1;46(3):798–798f.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27789670&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

22. 63.Kessler RC, Adler L, Ames M, Demler O, Faraone S, Hiripi E, et al. The World Health Organization Adult ADHD Self-Report Scale (ASRS): a short screening scale for use in the general population. Psychol Med. 2005 Feb;35(2):245–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291704002892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15841682&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231499500009&link_type=ISI) 

23. 64.Adler LA, Faraone SV, Sarocco P, Atkins N, Khachatryan A. Establishing US norms for the Adult ADHD Self-Report Scale (ASRS-v1.1) and characterising symptom burden among adults with self-reported ADHD. Int J Clin Pract. 2019 Jan;73(1):e13260.
    
    

24. 65.Rago D, Rasmussen MA, Lee-Sarwar KA, Weiss ST, Lasky-Su J, Stokholm J, et al. Fish-oil supplementation in pregnancy, child metabolomics and asthma risk. EBioMedicine. 2019 Aug;46:399–410.
    
    

25. 66.Stekhoven DJ, Bühlmann P. MissForest--non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012 Jan 1;28(1):112–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr597&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22039212&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298383800017&link_type=ISI) 

26. 67.Website [Internet]. Available from: Robinson D Bryan J Elias. 2020. “Fuzzyjoin: Join Tables Together on Inexact Matching. 2020” Version 0.1.6 (May). [https://cran.r-project.org/web/packages/fuzzyjoin/index.html](https://cran.r-project.org/web/packages/fuzzyjoin/index.html).
    
    

27. 68.Demontis D, Walters RK, Martin J, Mattheisen M, Als TD, Agerbo E, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. 2019 Jan;51(1):63–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0269-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30478444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

28. 69.Grove J, Ripke S, Als TD, Mattheisen M, Walters RK, Won H, et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. 2019 Mar;51(3):431–44.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0344-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30804558&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

29. 70.Ge T, Chen CY, Ni Y, Feng YCA, Smoller JW. Polygenic prediction via Bayesian regression and continuous shrinkage priors. Nat Commun. 2019 Apr 16;10(1):1776.
    
    

30. 71.PLINK 2.0 [Internet]. [cited 2023 Mar 29]. Available from: [https://www.cog-genomics.org/plink/2.0/](https://www.cog-genomics.org/plink/2.0/)
    
    

31. 72.Williams DR, Rast P. Back to the basics: Rethinking partial correlation network methodology. Br J Math Stat Psychol. 2020 May;73(2):187–212.
    
    

32. 73.Marsál K, Persson PH, Larsen T, Lilja H, Selbing A, Sultan B. Intrauterine growth curves based on ultrasonically estimated foetal weights. Acta Paediatr. 1996 Jul;85(7):843–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1651-2227.1996.tb14164.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8819552&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UX51600014&link_type=ISI) 

33. 74.Josse J, Husson F. missMDA: A Package for Handling Missing Values in Multivariate Data Analysis [Internet]. Vol. 70, Journal of Statistical Software. 2016. Available from: doi:10.18637/jss.v070.i01
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v070.i01&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom) 

34. 75.Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using**lme4** [Internet]. Vol. 67, Journal of Statistical Software. 2015. Available from: doi:10.18637/jss.v067.i01
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v067.i01&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23757445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F06%2F25%2F2024.03.07.24303907.atom)