Animal and plant protein intake during infancy and childhood DNA methylation: a meta-analysis in the NutriPROGRAM consortium

ABSTRACT Background Higher early-life animal protein intake is associated with a higher childhood obesity risk compared to plant protein intake. Differential DNA methylation may represent an underlying mechanism. Methods We analysed associations of infant animal and plant protein intakes with DNA methylation in early (2−6 years, N = 579) and late (7̄−12 years, N = 604) childhood in two studies. Study-specific robust linear regression models adjusted for relevant confounders were run, and then meta-analysed using a fixed-effects model. We also performed sex-stratified meta-analyses. Follow-up analyses included pathway analysis and eQTM look-up. Results Infant animal protein intake was not associated with DNA methylation in early childhood, but was associated with late-childhood DNA methylation at cg21300373 (P = 4.27 × 10¯8, MARCHF1) and cg10633363 (P = 1.09 × 10¯7, HOXB9) after FDR correction. Infant plant protein intake was associated with early-childhood DNA methylation at cg25973293 (P = 2.26 × 10−7, C1orf159) and cg15407373 (P = 2.13 × 10−7, MBP) after FDR correction. There was no overlap between the findings from the animal and plant protein analyses. We did not find enriched functional pathways at either time point using CpGs associated with animal and plant protein. These CpGs were not previously associated with childhood gene expression. Sex-stratified meta-analyses showed sex-specific DNA methylation associations for both animal and plant protein intake. Conclusion Infant animal protein intake was associated with DNA methylation at two CpGs in late childhood. Infant plant protein intake was associated with DNA methylation in early childhood at two CpGs. A potential mediating role of DNA methylation at these CpGs between infant protein intake and health outcomes requires further investigation.


Introduction
Dietary composition plays a crucial role in development during the early years of childhood.Protein intake is of particular importance as a source of essential amino acids [1].Recent observational studies and randomized controlled trials confirmed that high protein intake in early life increases early weight gain and the risk of later overweight and obesity, while growth can be restricted at lower intakes [2,3].Although the exact pathophysiological mechanisms are unclear, some studies have shown stronger associations of animal protein as compared to plant protein, especially during the first year of life, with higher body mass index (BMI) in later childhood [4,5] and especially higher body fat mass [4][5][6][7].The most common hypothesis explaining these associations between higher animal protein intake and obesity is the effects of amino acids, especially branchedchained amino acids (BCAA) in addition to other amino acids like arginine [8], assumed to stimulate insulin and insulin-like growth factor 1 (IGF-I) secretion as well as affect preadipocyte metabolism leading to overweight in children [9].
Differential DNA methylation might be a mechanism underlying these associations.Associations of early-life dietary intake with DNA methylation have been reported in several studies [10][11][12].For example, breastfeeding is associated with early-life growth, with DNA methylation being a potential mediating factor [13][14][15].Other studies found associations between the quality and quantity of dietary fat and fatty acids intake and DNA methylation in childhood [16][17][18].
Quality and quantity of infant protein intake have been previously investigated in association with DNA methylation in animal models [19][20][21], but never in humans.One study found that maternal protein restriction led to widespread differential methylation and gene expression in newborn rats [22].Few recent studies found associations between plasma protein metabolites, as an intermediate phenotype of protein intake, and DNA methylation in human adults [23,24].
Better knowledge on possible epigenetic pathways, looking at specific protein sources, might help further understand the underlying mechanisms of early-life dietary programming of later health.Therefore, we examined associations of dietary intake of animal protein in infancy with DNA methylation in early and late childhood.We additionally examined associations of dietary intake of plant protein with DNA methylation at the same time points.

Participants
Data from two studies were used in this metaanalysis; the CHOP trial with participants from five European countries, and the Generation R Study with participants from the Netherlands.The CHOP study is a multi-centre double-blind randomized clinical trial with 1,678 children enrolled in Germany, Belgium, Italy, Poland and Spain.The Generation R Study is a prospective populationbased cohort in Rotterdam, the Netherlands.In total, 9,778 mothers were enrolled in the Generation R Study.For CHOP, all children who attended the follow-up visits at ages of 5.5 and 11 years and who agreed on blood withdrawal and DNA methylation measurement were included.In the Generation R Study, a subgroup of children with European ethnic background was selected for DNA methylation measurement at ages of 6 years and 10 years, based on completeness of data.For CHOP, only children who were not breastfed during dietary data collection were included to dietary diary collection.Generation R included only 26 children who were breastfed.In both studies, energy from infant feeding was included in energy intake calculation.Study design and characteristics of both studies have been described previously in detail [25,26].From both studies, we included participants with information about infant animal protein intake and DNA methylation in early (2-6 years) and/or late childhood (7-12 years).Participants with energy intake from animal and plant protein outside ±5 standard deviations (SDs) from the study mean were excluded from the analysis (n = 3 and n = 0 in early childhood for CHOP and Generation R, respectively, n = 1 and n = 0 in late childhood for CHOP and Generation R, respectively).N = 183 and N = 227 were included in both time points for CHOP and Generation R, respectively.Neither of the studies had siblings in the studied sample.Participants with missing data on any of the covariates were excluded from the analysis (complete case analysis).All children in both studies had European ethnicity.Informed consent was obtained for all participants, and both studies received approval from the local ethics committees.

Animal and plant protein intake
Animal protein intake was defined as the average daily protein intake from animal food sources, including animal dairy, meat, eggs, and fish combined.Plant protein intake was defined as the average daily protein intake from non-animal sources.Both intakes were measured at approximately 12 and 14 months of age in CHOP and Generation R, respectively.Because we were interested in relative protein intake and to account for confounding by energy intake, we expressed animal protein and plant protein intake as percentage of total energy intake (E%) and additionally included energy intake in our models.In CHOP, dietary intakes were recorded using prospective 3-day dietary diary protocol following standardized operating procedures [27].Nutrient information was primarily based on the German national food composition database (BLS 3.01).Food items and recipes not identified in the database were added by CHOP dietitians at each study centre according to information from the manufacturers, other databases or ingredients.Food records with energy or macronutrient intakes > 3 SDs of the mean by month and country and those noted by the data entering dietitian to be incomplete or inaccurate or with reported concurrent illness were excluded [28].In the Generation R Study, dietary intake was evaluated using a 211item semi-quantitative Food Frequency Questionnaire (FFQ) covering the previous month.The Dutch Food Composition Table 2006 and standard Dutch portion sizes were used to convert food frequencies into energy and macronutrient intakes [29].

DNA methylation data acquisition, quality control, and normalization
DNA was isolated from peripheral whole blood samples in both studies, using the salting out method.For early childhood, the average age at blood draw was 5.5 and 6 years for CHOP and Generation R, respectively, while it was 11 and 10 years for CHOP and Generation R, respectively, for late childhood.For CHOP and Generation R, 800 and 500 ng of DNA, respectively, per sample underwent bisulphite conversion.DNA methylation in early and late childhood was measured using the Illumina Infinium® HumanMethylation450 BeadChip assay.In both studies, DNA methylation data normalization and quality control were done following the quantile normalization method (CPACOR) by Lehne et.al [30].DNA methylation was analysed as untransformed DNA methylation beta-values ranging from 0 (completely unmethylated) to 1 (completely methylated) expressing the proportion of cells in which the DNA was methylated at a specific cytosine-phosphate-guanine (CpG) site.Nonautosomal probes were excluded.Control probes were removed as part of study-level QC performed on DNA methylation data.The top and bottom 0.5% methylation beta-values were winsorized.

Covariates
Models were adjusted for total energy intake (from FFQ or dietary protocol), age at blood collection, sex, self-reported highest completed maternal educational level (lower/higher), gestational age at birth, selfreported maternal pre-pregnancy BMI, self-reported sustained maternal smoking during pregnancy into at least the second trimester (yes/no), blood cell proportions (B-cells, CD8+ T-cells, CD4+ T-cells, granulocytes, NK-cells, and monocytes) estimated using the reference-based Houseman method [31] in the minfi package in R [32,33], with the Reinius reference [34] in both studies.To control for technical batch effects, both studies adjusted for sample plate number as a technical covariate.For CHOP, the country of study centre was also included in all models.

Statistical analysis
Robust linear regression models were run in both studies to analyse the associations of both infant animal and plant protein intake and DNA methylation at both time points following a pre-specified analysis plan and R code.For animal protein intake, each study ran eight robust linear regression models (four models at each of the two time points) using the rlm() function of the 'MASS' package in R [33,35], with model 4 (fully adjusted model) being the main model.For covariance matrix estimation, 'Heteroscedasticity-Consistent Covariance Matrix Estimation' method was used.For plant protein intake, each study ran the fully adjusted model at each time point.To account for confounding by energy intake, the multivariate nutrient density model was used, expressing animal protein in E% and additionally including energy intake as covariate: (1) DNAm ~ Infant animal/plant protein intake + age at blood collection + sex + technical covariates (2) Model 1 plus adjustment for total energy intake (3) Model 2 plus adjustment for maternal educational level + gestational age at birth + pre-pregnancy BMI + maternal smoking during pregnancy (4) Model 3 plus adjustment for cell counts To identify sex-specific DNA methylation differences, we also performed sex-stratified metaanalyses for both animal and plant protein intake at both time points using the fully adjusted model 4. Before meta-analysis, quality control was done on the results of both studies using the QCEWAS R package [36].Results from both studies were centrally meta-analysed using METAL [37] by MeS, using fixed effects inverse varianceweighted meta-analysis in METAL.The metaanalysis was then independently repeated by JFF, to exclude the possibility of human error.Removing cross-reactive probes [38,39] (N = 43,635 in early childhood and N = 43,702 in late childhood) and probes available only in one study (N = 13,428 in early childhood and N = 14,799 in late childhood) resulted in 415,274 and 414,307 probes included in the meta-analyses for early and late childhood, respectively.We flagged CpGs potentially influenced by a single nucleotide polymorphism (polymorphic sites) [38,39] and CpGs listed as methylation quantitative trait loci (mQTLs) [40].Hartigans' dip test for unimodality was used to check if methylation level at significant CpGs was influenced by a nearby single nucleotide polymorphism.We also flagged CpGs with high between-study heterogeneity, defined as an I 2 >50%.We decided a priori on a false discovery rate (FDR)-adjusted p-value significance level of < 0.05 [41].All statistical analyses were performed using R statistical analysis software [33], except for the meta-analysis in METAL.

Additional analyses
Following the meta-analyses, multiple analyses were performed on the top CpGs (FDR P < 0.05) from the fully adjusted non sex-stratified metaanalyses at both time points, to examine potential functional consequences.First, the EWAS catalog [42] was used to identify previously published associations with other traits for the top CpGs.Second, we checked if the top CpGs were previously reported to be associated with gene expression in childhood blood, measured as expression quantitative trait methylation (eQTM).For this, we used the catalogue of 13.6 million blood autosomal cis-eQTMs in children published by the Human Early Life Exposome (HELIX) project, after cell type adjustment [43].Third, we ran functional enrichment analyses using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) in the MissMethyl R package [44] on all CpGs with P < 1×10¯4 in the fully adjusted models at both time points for both animal and plant protein intakes.

Participants
In total, 579 and 604 children were included in the animal protein meta-analyses and 577 and 604 in the plant protein meta-analysis for early and late childhood, respectively.Study-specific descriptives are presented in Tables 1 and 2. Infant animal and plant protein intakes were normally distributed in both studies.Animal protein had means (SD) of 10.6% (3.2) and 10.7% (3.5) in CHOP for the analyses in early and late childhood DNA methylation, respectively, and 8.1% (2.4) and 8.1% (2.4) in Generation R for analyses in early and late childhood DNA methylation, respectively.Plant protein intake had means (SD) of 4.04 (1.03) and 3.93 (1.06) in CHOP for the analyses in early and late childhood DNA methylation, respectively, and 4.67 (1.31) and 4.76 (1.26) in Generation R for analyses in early and late childhood DNA methylation, respectively (Supplementary Figures S1A-1D and 2A-2D).

Meta-analysis
Results of the epigenome-wide association study (EWAS) meta-analysis for the fully adjusted models (model 4) at both time points, for CpGs with P < 1×10¯5 are presented in Tables 3 and 4 for animal protein and in Tables 5 and 6 for plant protein.Sex-stratified EWAS results are presented in Supplementary tables S1-S8.Full EWAS meta-analysis results are available via: https://doi.org/10.5281/zenodo.8375454.
Q-Q, Manhattan and Volcano plots for the metaanalysis results of the fully adjusted model at both time points are presented in Figure 1a-f for animal protein and Figure 2a-f for plant protein.Cohort-specific Q-Q plots and lambdas are shown in Supplementary Figures S3A-3D for animal protein and S4A-4D for plant protein.Q-Q plots for sex-stratified meta-analysis results are shown in Supplementary figures S5A-5D for animal protein and Supplementary figures S6A-6D for plant protein.Infant animal protein intake was not associated with DNA methylation at any CpG site in early childhood, after FDR correction.Infant animal protein intake was statistically significantly associated with DNA methylation in late childhood at two CpG sites: cg21300373 (0.16% increase in DNA methylation per each 1 E% increase in infant animal protein intake, SE = 0.0003, P = 4.27×10 −8 , P FDR = 0.018, I 2 = 0) and cg10633363 (0.05% decrease in DNA methylation per each 1 E% increase in infant animal protein intake, SE = 0.0001, P = 1.09×10 −7 , P FDR = 0.022, I 2 = 51) after FDR correction.Forest plots for both CpGs associated with animal protein intake are presented in Supplementary figures S7 and 8 and Supplementary figures S9 and 10 for plant protein.Even though cg10633363 had an I 2 of 50.6, the forest plots did not show clear heterogeneity between studies.Both CpGs were flagged as potentially associated with a SNP (mQTLs).Although the density plots showed a slight indication of multimodality (Supplementary figures S11 and 12), Hartigans' dip test had a P > 0.05 for both CpGs in both studies, indicating no significant deviation from unimodality.
Plant protein intake was not associated with DNA methylation at any CpG site in late childhood, after FDR correction.Neither of the two CpGs associated with animal protein intake were associated with plant protein intake at P < 1×10¯5, or vice versa.
Sex-stratified meta-analyses showed sex-specific DNA methylation differences in both early and late childhood with infant animal protein intake.Infant animal protein intake was associated with early childhood DNA methylation at 16 CpG sites in boys and at 9 CpG sites in girls, after FDR correction.Infant animal protein intake was associated with late childhood DNA methylation at 11 CpG sites in boys and at 1 CpG site in girls, after FDR correction.Cg21300373, one of the CpGs significantly associated with animal protein in late childhood in both sexes combined, was only significant in boys (Effect = 0.002, SE = 0.0004, P = 1.24×10 −6 , P FDR = 0.046, I 2 = 3.9) and not in girls (Effect = 0.001, SE = 0.0004, P = 7.6×10 −4 , P FDR = 0.34, I 2 = 0), after FDR correction.Cg10633363, the other CpG significantly associated with animal protein in late childhood in both sexes combined, was not significant neither in boys (Effect = 0.001, SE = 0.0003, P = 2.41×10 −4 , P FDR = 0.34, I 2 = 0) nor in girls (Effect = 0.001, SE = 0.0003, P = 1.17×10 −5 , P FDR = 0.14, I 2 = 57) (Supplementary tables S1-S4).
For infant plant protein intake, sex-stratified meta-analysis showed associations with DNA methylation in early childhood at 1 CpG site for boys and at 3 CpG sites for girls, respectively.Associations of infant plant protein intake with DNA methylation at late childhood were found at 9 CpG sites for boys and at 1 CpG site for girls, respectively.None of the findings in the sexstratified plant protein meta-analyses was found in the other sex or in any of the other sex-stratified models (Supplementary tables S5-S8).

Functional analysis
Cg21300373 is mapped to MARCHF1 and cg10633363 is mapped to HOXB9, both CpGs are located in the transcription start site.We used the CpGs associated with animal protein for functional analyses.First, a look-up in the EWAS catalog showed that DNA methylation at cg21300373 was previously associated with pancreatic ductal adenocarcinoma and clear cell renal carcinoma [45,46].Cg10633363 was not reported in the EWAS catalog.Second, cg21300373 and cg10633363 were not associated with gene expression in the catalogue of blood autosomal cis-eQTMs in children [47].Third, functional enrichment analysis was done using the CpGs with p-values <1×10¯4.For animal protein models, 120 CpGs for early childhood and 98 CpGs for late childhood models were used as input.For plant protein, 88 CpGs for early childhood and 89 CpGs for late childhood were used as input.There were no enriched GO or KEGG pathways (FDR <0.05) at either animal or plant protein intake at either time point.Cg25973293 associated with infant plant protein intake at early childhood is mapped to C1orf159.Cg25973293 was not reported to be associated with other phenotypes in the EWAS catalog.

Discussion
In this meta-analysis, infant animal protein intake was associated with DNA methylation in late childhood at the age of 7-12 years at two CpG sites; cg21300373 and cg10633363, but not with DNA methylation in early childhood.No associations with specific functional pathways or gene expression were identified.Infant plant protein intake was associated with DNA methylation in early childhood at the age of 2-6 years at two CpG sites; cg25973293 and cg15407373.Cg21300373 was mapped to the transcription start site of MARCHF1 which is a member of the MARCH family of membrane-bound E3 ubiquitin ligases and was linked to glucose-tolerance and lipid storage [48].E3 ubiquitin ligases play a pivotal role in obesity-induced insulin resistance in humans [49].DNA methylation of MARCHF1 gene was found to be associated with adiposity in a previous EWAS in multi-ethnic Asian adults [50].MARCHF1 genetic variants have also been associated with adiposity in adults [51,52].Therefore, this may be a first link between protein intake-related DNA methylation and obesity in children.
Cg10633363 was mapped to the transcription start site of HOXB9.HOXB9 is one of the HOX genes, a group of related genes responsible for mapping body organs along the head-tail axis during embryonic development.HOXB9 has not been linked to body composition.DNA methylation levels at CpG sites in or close to HOXB9 have been found to be associated with epigenetic ageing of liver tissue [53], a process that was previously found to be accelerated by obesity [54].
Cg25973293, associated with plant protein, is mapped to C1orf159 gene in an open sea.C1orf159 is a protein coding gene whose increased expression was found to be an unfavourable prognostic marker in renal [55] and liver [56] cancer.Despite being significant after FDR correction, cg15407373 had heterogeneity I 2 = 95.6 denoting high heterogeneity between results from both analysed studies results.Cg15407373, associated with plant protein, is mapped to MBP gene in an island.MBP gene is a protein coding gene responsible for encoding Myelin Basic Protein, a protein playing a pivotal role in nervous system development [57].
No associations of genetic variants of HOXB9, C1orf159 and MBP genes with metabolic phenotypes have been reported in the GWAS catalog [58].Although all CpGs discovered in animal protein and plant protein EWAS are listed to be associated with a genetic variant, the level of methylation at these sites did not seem to be strongly influenced by a nearby polymorphic site in our populations, as the distribution of the methylation values for both sites in both studies did not differ significantly from unimodality.Using CpGs associated with infant animal and plant protein intakes with P < 1×10¯4 for early and late childhood models, no enriched functional pathways or biological processes were found.Further research is needed to confirm any mediating role for DNA methylation in the associations of animal protein intake and health in later life.The infant plant protein intake EWAS showed no association for the two CpGs discovered in infant animal protein EWAS with DNA methylation at either of the two time points investigated.Sexspecific DNA methylation differences were identified at both time points for both infant animal and plant protein intakes.
Previous studies showed associations between early life protein intake, especially animal protein, and body composition, mainly higher fat mass [59].Although studies in animal models have reported that quantity and quality of protein intake in pups are associated with DNA methylation [21,60], whether this may be a potential underlying mechanism in the association of protein with adiposity has not been investigated in humans [4,5,61].
Perinatal exposures, such as early nutrition, are important for early health programming.It is also becoming clear that nutritional associations with DNA methylation are not limited to perinatal period, but further extend to postnatal life [62].Infant breastfeeding was found to be associated with child's growth and development, with DNA methylation being a potential mediating factor [14].Our study contributes to this evidence, linking infant animal protein intake with DNA methylation.These associations appeared only later in childhood, suggesting an association with extended exposure.Furthermore, different associations were found for infant animal and plant protein intakes, suggesting a role for the difference in metabolism of animal and plant protein.The quality of protein and availability of methyl group donors might play a role in differential methylation of some CpG sites.Animal protein as a nutritional source is rich in methionine [63], a precursor to the universal methyl-donor S-adenosylmethionine, which is assumed to change the DNA methylation pattern, especially at specific loci rather than genome-wide [64].
Sex-specific differences were observed for infant animal and plant protein intakes which might be attributed to disparities in protein metabolism between the sexes [59,65].Further research including more studies at multiple time points might better reveal the association of animal protein intake and DNA methylation throughout the life course and whether associations of animal protein with DNA methylation indeed mediate associations with body composition.
This study has a number of strengths.We defined a strict age window for infant animal and plant protein intakes, because dietary composition in early life can change drastically over time [66].Both cohorts included in this meta-analysis have detailed dietary data, as well as information on relevant confounders.The analyses were run using a pre-specified analysis plan.We included motheroffspring data from multiple European countries.The multivariate nutrient density model was used to account for total energy intake [67].However, this study has a limited sample size, which may have limited our ability to find smaller effect sizes.As macronutrient intake is known to be correlated with total energy intake, analyses may be confounded by total energy intake [67], to overcome this, we adjusted for total energy intake in the main models while using E% from animal protein sources as the exposure.With this adjustment we can interpret the results as the effect of isocaloric replacement of animal protein with any other macronutrient.
Children from both cohorts were of European ethnicity which might limit generalizability of the findings to other ethnic groups, however it reduces the heterogeneity between the two studies included.Both cohorts relied on self-reported dietary assessment questionnaires, which might have introduced measurement error.Despite both cohorts using different dietary assessment methods, forest plots for significant CpGs showed consistent results from both cohorts with slightly different effect sizes, which might be a reason for the marginally high I 2 (50.6%) for cg10633363 [68].
In conclusion, this meta-analysis showed associations of infant animal protein intake with DNA methylation at cg21300373 and cg10633363 in late childhood.No associations were found between infant animal protein intake and DNA methylation in early childhood.Association of infant plant protein intake with DNA methylation was found at cg25973293 and cg15407373 in early childhood.No associations were found between plant protein and DNA methylation in late childhood.Potential mediating epigenetic pathways between infant protein intake and health outcomes require further investigation using larger sample sizes.

Figure 1 .
Figure 1.Q-Q, Manhattan and volcano plots for associations between infant animal protein intake (E%) and early and late childhood DNA methylation.Significance line in Manhattan and volcano plots is set to represent the Bonferroni-corrected p-value threshold of P < 1.1 × 10¯7.

Figure 2 .
Figure 2. Q-Q, Manhattan and volcano plots for associations between infant plant protein intake (E%) and early and late childhood DNA methylation.Significance line in Manhattan and volcano plots is set to represent the Bonferroni-corrected p-value threshold of P < 1.1 × 10¯7.

Table 1 .
Study-specific descriptives for animal protein analysis.
*In both cohorts, maternal educational level was defined as 'Low' if no, primary or secondary school education obtained, and as 'High' if college or university degree obtained.

Table 2 .
Study-specific descriptives for plant protein analysis.

Table 3 .
Meta-analysis results for associations of infant animal protein intake (E%) and early childhood whole blood DNA methylation levels at P < 1 × 10¯5 . .Effect sizes and SE are presented as a percent change in DNA methylation per 1 E% increase in energy from animal protein intake.
a b Studies are arranged from left to right as follows: CHOP then Generation R. c methylation quantitative trait loci (mQTL) utilized from GoDMC database.EWAS model was adjusted for total energy intake, age at blood collection, sex, maternal educational level, gestational age, pre-pregnancy BMI, maternal smoking during pregnancy, technical covariates, cell counts.

Table 4 .
Meta-analysis results for associations of infant animal protein intake (E%) and late childhood whole blood DNA methylation levels at P < 1 × 10¯5 . .
a Effect sizes and SE are presented as a percent change in DNA methylation per 1 E% increase in energy from animal protein intake.b Studies are arranged from left to right as follows: CHOP then Generation R. c methylation quantitative trait loci (mQTL) utilized from GoDMC database.EWAS model was adjusted for total energy intake, age at blood collection, sex, maternal educational level, gestational age, pre-pregnancy BMI, maternal smoking during pregnancy, technical covariates, cell counts.

Table 5 .
Meta-analysis results for associations of infant plant protein intake (E%) and early childhood whole blood DNA methylation levels at P < 1 × 10¯5 . .
a Effect sizes and SE are presented as a percent change in DNA methylation per 1 E% increase in energy from animal protein intake.b Studies are arranged from left to right as follows: CHOP then Generation R. c methylation quantitative trait loci (mQTL) utilized from GoDMC database.EWAS model was adjusted for total energy intake, age at blood collection, sex, maternal educational level, gestational age, pre-pregnancy BMI, maternal smoking during pregnancy, technical covariates, cell counts.

Table 6 .
Meta-analysis results for associations of infant plant protein intake (E%) and late childhood whole blood DNA methylation levels at P < 1 × 10¯5 . .CpG Effect a SE a P value FDR Direction b I 2 Polymorphic mQTL c Chr Position Effect sizes and SE are presented as a percent change in DNA methylation per 1 E% increase in energy from animal protein intake.Studies are arranged from left to right as follows: CHOP then Generation R. methylation quantitative trait loci (mQTL) utilized from GoDMC database.EWAS model was adjusted for total energy intake, age at blood collection, sex, maternal educational level, gestational age, pre-pregnancy BMI, maternal smoking during pregnancy, technical covariates, cell counts.