Prenatal metal exposures and childhood gut microbial signatures are associated with depression score in late childhood

BACKGROUND

Childhood depression is a major public health issue worldwide. Previous studies have linked both prenatal metal exposures and the gut microbiome to depression in children. However, few, if any, have studied their interacting effect in specific subgroups of children.

OBJECTIVES

Using an interpretable machine-learning method, this study investigates whether children with specific combinations of prenatal metals and childhood microbial signatures (cliques or groups of metals and microbes) were more likely to have higher depression scores at 9-11 years of age.

METHODS

We leveraged data from a well-characterized pediatric longitudinal birth cohort in Mexico City and its microbiome substudy (n = 112). Eleven metal exposures were measured in maternal whole blood samples in the second and third trimesters of pregnancy. The gut microbial abundances were measured at 9-11-year-olds using shotgun metagenomic sequencing. Depression symptoms were assessed using the Child Depression Index (CDI) t-scores at 9-11 years of age. We used Microbial and Chemical Exposure Analysis (MiCxA), which combines interpretable machine-learning into a regression framework to identify and estimate joint associations of metal-microbial cliques in specific subgroups. Analyses were adjusted for relevant covariates.

RESULTS

We identified a subgroup of children (11.6 % of the sample) characterized by a four-component metal-microbial clique that had a significantly high depression score (15.4 % higher than the rest) in late childhood. This metal-microbial clique consisted of high Zinc in the second trimester, low Cobalt in the third trimester, a high abundance of Bacteroides fragilis, a high abundance of Faecalibacterium prausnitzii. All combinations of cliques (two-, three-, and four-components) were significantly associated with increased log-transformed t-scored CDI (β = 0.14, 95%CI = [0.05,0.23], P < 0.01 for the four-component clique).

SIGNIFICANCE

This study offers a new approach to chemical-microbial analysis and a novel demonstration that children with specific gut microbiome cliques and metal exposures during pregnancy may have a higher likelihood of elevated depression scores.

Prenatal Lead Exposure Is Associated with Reduced Abundance of Beneficial Gut Microbial Cliques in Late Childhood: An Investigation Using Microbial Co-Occurrence Analysis (MiCA)

Many analytical methods used in gut microbiome research focus on either single bacterial taxa or the whole microbiome, ignoring multibacteria relationships (microbial cliques). We present a novel analytical approach to identify microbial cliques within the gut microbiome of children at 9-11 years associated with prenatal lead (Pb) exposure. Data came from a subset of participants (n = 123) in the Programming Research in Obesity, Growth, Environment and Social Stressors cohort. Pb concentrations were measured in maternal whole blood from the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the gut microbiome. Using a novel analytical approach, Microbial Co-occurrence Analysis (MiCA), we paired a machine learning algorithm with randomization-based inference to first identify microbial cliques that were predictive of prenatal Pb exposure and then estimate the association between prenatal Pb exposure and microbial clique abundance. With second-trimester Pb exposure, we identified a two-taxa microbial clique that included Bifidobacterium adolescentis and Ruminococcus callidus and a three-taxa clique that also included Prevotella clara. Increasing second-trimester Pb exposure was associated with significantly increased odds of having the two-taxa microbial clique below the median relative abundance (odds ratio (OR) = 1.03, 95% confidence interval (CI) [1.01-1.05]). Using a novel combination of machine learning and causal inference, MiCA identified a significant association between second-trimester Pb exposure and the reduced abundance of a probiotic microbial clique within the gut microbiome in late childhood.

Cross-sectional associations between phthalates, phenols, and parabens with metabolic syndrome risk during early-to-mid adolescence among a cohort of Mexican youth

BACKGROUND

Epidemiological studies on children and adults have linked toxicants from plastics and personal care products to metabolic disruption. Yet, the impact of endocrine-disrupting chemicals (EDCs) on adolescent metabolic syndrome (MetS) risk during early and mid-adolescence is unclear.

METHODS

To examine the links between exposure to EDCs and MetS risk and its components, cross-sectional data from 344 Mexican youth in early-to-mid adolescence (10-17 years) were analyzed. Urinary biomarker concentrations of phthalates, phenol, and paraben analytes were measured from a single spot urine sample collected in 2015; study personnel obtained anthropometric and metabolic measures. We examined associations between summary phthalates and metabolites, phenol, and paraben analytes with MetS risk z-scores using linear regression, adjusted for specific gravity, sex, age, pubertal status, smoking, alcohol intake, physical activity level, and screen time. As a secondary aim, mediation analysis was conducted to evaluate the role of hormones in the association between summary phthalates with lipids and MetS risk z-scores.

RESULTS

The mean (SD) age was 13.2 (1.9) years, and 50.9% were female. Sex-stratified analyses revealed associations between summary phthalates and lipids ratio z-scores, including Σ DEHP [β = 0.21 (95% CI: 0.04, 0.37; p < 0.01)], phthalates from plastic sources (Σ Plastic) [β = 0.22 (95% CI: 0.05, 0.39; p < 0.01)], anti-androgenic phthalates (Σ AA) [β = 0.22 (95% CI: 0.05, 0.39; p < 0.01)], and individual phthalate metabolites (MEHHP, MEOHP, and MECPP) among males. Among females, BPA [β = 0.24 (95% CI: 0.03, 0.44; p < 0.05)] was positively associated with lipids ratio z-score and one phenol (2,5 DCP) [β = 0.09 (95% CI: 0.01, 0.18); p < 0.05)] was associated with increased waist circumference z-score. Results showed no evidence of mediation by hormone concentrations in the association between summary phthalates with lipids ratio or MetS risk z-scores.

CONCLUSION

Higher EDC exposure was positively associated with serum lipids during adolescence, particularly among males.

Prenatal lead exposure is negatively associated with the gut microbiome in childhood

Background

Metal exposures are associated with gut microbiome (GM) composition and function, and exposures early in development may be particularly important. Considering the role of the GM in association with many adverse health outcomes, understanding the relationship between prenatal metal exposures and the GM is critically important. However, there is sparse knowledge of the association between prenatal metal exposure and GM later in childhood.

Objectives

This analysis aims to identify associations between prenatal lead (Pb) exposure and GM composition and function in children 9-11 years old.

Methods

Data come from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) cohort based in Mexico City, Mexico. Prenatal metal concentrations were measured in maternal whole blood drawn during the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the GM. This analysis uses multiple statistical modeling approaches, including linear regression, permutational analysis of variance, weighted quantile sum regression (WQS), and individual taxa regressions, to estimate the association between maternal blood Pb during pregnancy and multiple aspects of the child GM at 9-11 years old, adjusting for relevant confounders.

Results

Of the 123 child participants in this pilot data analysis, 74 were male and 49 were female. Mean prenatal maternal blood Pb was 33.6 (SE = 2.1) ug/L and 34.9 (SE = 2.1) ug/L at second and third trimesters, respectively. Analysis suggests a consistent negative relationship between prenatal maternal blood Pb and the GM at age 9-11, including measures of alpha and beta diversity, microbiome mixture analysis, and individual taxa. The WQS analysis showed a negative association between prenatal Pb exposure and the gut microbiome, for both second and third trimester exposures (2Tβ = -0.17, 95%CI = [-0.46,0.11]; 3Tβ = -0.17, 95%CI = [-0.44,0.10]). Ruminococcus gnavus, Bifidobacterium longum, Alistipes indistinctus, Bacteroides caccae, and Bifidobacterium bifidum all had weights above the importance threshold from 80% or more of the WQS repeated holdouts in association with both second and third trimester Pb exposure.

Discussion

Pilot data analysis suggests a negative association between prenatal Pb exposure and the gut microbiome later in childhood; however, additional investigation is needed.

Prenatal dietary patterns in relation to adolescent offspring adiposity and adipokines in a Mexico City cohort

Maternal diet during pregnancy has been associated with obesity among offspring. The extent to which trimester-specific dietary patterns are associated with markers of adiposity during adolescence remains unclear. We examined associations between prenatal diet patterns with adolescent offspring measures of adiposity and adipokines in 384 mother-adolescent dyads from the Mexico City ELEMENT cohort. Trimester-specific diet patterns were derived from principal component analysis of food frequency questionnaire data. Adolescent anthropometry and serum leptin and adiponectin were measured at 10-17 years. Three maternal diet patterns were identified: Prudent Diet (PD), high in fish and vegetables, the High Meat and Fat Diet (HMFD), high in pork and processed meats, and the Transitioning Mexican Diet (TMD), high in corn tortillas and sugar-sweetened beverages. Multiple linear regression was used to estimate sex-stratified associations among quartiles of diet patterns with adiposity and adipokines, adjusting for maternal marital status, education, and parity. First trimester TMD was associated with greater anthropometric measures and higher leptin in females, while third trimester HMFD was associated higher body fat percentage, triceps thickness, waist circumference, and leptin, but lower adiponectin among males. Contrary to expectation, there were positive associations between the trimester 1 PD pattern and anthropometric measurements in females, and for trimester 2 HMFD and TMD patterns with adipokines among males. Findings suggest maternal diet patterns may influence offspring adiposity markers during adolescence in a sex-specific manner.

Prenatal Pb exposure is associated with reduced abundance of beneficial gut microbial cliques in late childhood: an investigation using Microbial Co-occurrence Analysis (MiCA)

Background

Many analytical methods used in gut microbiome research focus on either single bacterial taxa or the whole microbiome, ignoring multi-bacteria relationships (microbial cliques). We present a novel analytical approach to identify multiple bacterial taxa within the gut microbiome of children at 9-11 years associated with prenatal Pb exposure.

Methods

Data came from a subset of participants (n=123) in the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) cohort. Pb concentrations were measured in maternal whole blood from the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the gut microbiome. Using a novel analytical approach, Microbial Co-occurrence Analysis (MiCA), we paired a machine-learning algorithm with randomization-based inference to first identify microbial cliques that were predictive of prenatal Pb exposure and then estimate the association between prenatal Pb exposure and microbial clique abundance.

Results

With second-trimester Pb exposure, we identified a 2-taxa microbial clique that included Bifidobacterium adolescentis and Ruminococcus callidus, and a 3-taxa clique that added Prevotella clara. Increasing second-trimester Pb exposure was associated with significantly increased odds of having the 2-taxa microbial clique below the 50th percentile relative abundance (OR=1.03,95%CI[1.01-1.05]). In an analysis of Pb concentration at or above vs. below the United States and Mexico guidelines for child Pb exposure, odds of the 2-taxa clique in low abundance were 3.36(95%CI[1.32-8.51]) and 6.11(95%CI[1.87-19.93]), respectively. Trends were similar with the 3-taxa clique but not statistically significant.

Discussion

Using a novel combination of machine-learning and causal-inference, MiCA identified a significant association between second-trimester Pb exposure and reduced abundance of a probiotic microbial clique within the gut microbiome in late childhood. Pb exposure levels at the guidelines for child Pb poisoning in the United States, and Mexico are not sufficient to protect against the potential loss of probiotic benefits.

Prenatal Lead Exposure is Negatively Associated with the Gut Microbiome in Childhood

Background

Metal exposures are associated with gut microbiome (GM) composition and function, and exposures early in development may be particularly important. Considering the role of the GM in association with many adverse health outcomes, understanding the relationship between prenatal metal exposures and the GM is critically important. However, there is sparse knowledge of the association between prenatal metal exposure and GM later in childhood.

Objectives

This analysis aims to identify associations between prenatal lead (Pb) exposure and GM composition and function in children 9-11 years old.

Methods

Data come from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) cohort based in Mexico City, Mexico. Prenatal metal concentrations were measured in maternal whole blood drawn during the second and third trimesters of pregnancy. Stool samples collected at 9-11 years old underwent metagenomic sequencing to assess the GM. This analysis uses multiple statistical modeling approaches, including linear regression, permutational analysis of variance, weighted quantile sum regression (WQS), and individual taxa regressions, to estimate the association between maternal blood Pb during pregnancy and multiple aspects of the child GM at 9-11 years old, adjusting for relevant confounders.

Results

Of the 123 child participants in this pilot data analysis, 74 were male and 49 were female. Mean prenatal maternal blood Pb was 33.6(SE=2.1) ug/L and 34.9(SE=2.1) ug/L at second and third trimesters, respectively. Analysis suggests a consistent negative relationship between prenatal maternal blood Pb and the GM at age 9-11, including measures of alpha and beta diversity, microbiome mixture analysis, and individual taxa. The WQS analysis showed a negative association between prenatal Pb exposure and the gut microbiome, for both second and third trimester exposures (2Tβ=-0.17,95%CI=[-0.46,0.11]; 3Tβ=-0.17,95%CI=[-0.44,0.10]). Ruminococcus gnavus, Bifidobacterium longum, Alistipes indistinctus, Bacteroides caccae, and Bifidobacterium bifidum all had weights above the importance threshold from 80% or more of the WQS repeated holdouts in association with both second and third trimester Pb exposure.

Discussion

Pilot data analysis suggests a negative association between prenatal Pb exposure and the gut microbiome later in childhood; however, additional investigation is needed.

Associations between exposure to phthalates, phenols, and parabens with objective and subjective measures of sleep health among Mexican women in midlife: a cross-sectional and retrospective analysis

Endocrine-disrupting chemicals (EDCs) may impact sleep during the menopausal transition by altering sex hormones. However, these studies are scarce among Latin American women. This investigation utilized cross-sectional and retrospective data from midlife women enrolled in the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) study to examine associations between exposure to EDCs (phthalates, phenols, and parabens) and sleep health measures. For cross-sectional analyses, single spot urine samples were collected between 2017-2019 from a pilot sample of women (N = 91) of midlife age to estimate the urinary concentration of individual phthalates, phenols, and parabens and to calculate the summary concentration of phthalate mixtures. Seven-day nightly sleep duration, midpoint, and fragmentation were obtained from wrist-actigraphy devices and estimated from the actigraphy data using a pruned dynamic programming algorithm. Self-reported poor sleep quality was assessed by one item from the Pittsburgh Sleep Quality Index (PSQI). We examined associations between urinary summary phthalate mixtures, phthalate metabolites, phenol, and paraben analytes with each sleep measure using linear or logistic (to compute odds of poor sleep quality only) regression models adjusted for specific gravity, age, and socioeconomic status. We ran similar regression models for retrospective analyses (N = 74), except that urine exposure biomarker data were collected in 2008 when women were 24-50 years old. At the 2017-2019 midlife visit, 38% reported poor sleep quality. Cross-sectionally, EDCs were associated with longer sleep duration, earlier sleep timing, and more fragmented sleep. For example, every 1-unit IQR increase in the phenol triclosan was associated with a 26.3 min per night (95% CI: 10.5, 42.2; P < 0.05) longer sleep duration and marginally associated with 0.2 decimal hours (95% CI: -0.4, 0.0; P < 0.10) earlier sleep midpoint; while every 1-unit IQR increase in the phthalate metabolite MEHP was associated with 1.1% higher sleep fragmentation (95% CI: 0.1, 2.1; P < 0.05). Retrospective study results generally mirrored cross-sectional results such that EDCs were linked to longer sleep duration, earlier sleep timing, and more fragmented sleep. EDCs were not significantly associated with odds of self-reported poor sleep quality. Results from cross-sectional and retrospective analyses revealed that higher exposure to EDCs was predictive of longer sleep duration, earlier sleep timing, and more fragmented sleep among midlife women.

DNA Methylation Is a Potential Biomarker for Cardiometabolic Health in Mexican Children and Adolescents

DNA methylation (DNAm) is a plausible mechanism underlying cardiometabolic abnormalities, but evidence is limited among youth. This analysis included 410 offspring of the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) birth cohort followed up to two time points in late childhood/adolescence. At Time 1, DNAm was quantified in blood leukocytes at long interspersed nuclear elements (LINE-1), H19, and 11β-hydroxysteroid dehydrogenase type 2 (11β-HSD-2), and at Time 2 in peroxisome proliferator-activated receptor alpha (PPAR-α). At each time point, cardiometabolic risk factors were assessed including lipid profiles, glucose, blood pressure, and anthropometry. Linear mixed effects models were used for LINE-1, H19, and 11β-HSD-2 to account for the repeated-measure outcomes. Linear regression models were conducted for the cross-sectional association between PPAR-α with the outcomes. DNAm at LINE-1 was associated with log glucose at site 1 [β = -0.029, p = 0.0006] and with log high-density lipoprotein cholesterol at site 3 [β = 0.063, p = 0.0072]. 11β-HSD-2 DNAm at site 4 was associated with log glucose (β = -0.018, p = 0.0018). DNAm at LINE-1 and 11β-HSD-2 was associated with few cardiometabolic risk factors among youth in a locus-specific manner. These findings underscore the potential for epigenetic biomarkers to increase our understanding of cardiometabolic risk earlier in life.

Duration, timing, and consistency of sleep in relation to inflammatory cytokines in Mexican adolescents

OBJECTIVE

To evaluate whether sleep duration, timing, and variability were associated with inflammatory cytokines in a cohort of Mexico City adolescents.

METHODS

The analytic sample comprised >500 adolescents who were part of an ongoing longitudinal study in Mexico City. At two time points during mid-to-late puberty (average age 14, n = 391) and late-to-post puberty (average age 16, n = 345), adolescents completed a follow-up visit that included 7-day wrist actigraphy and clinical assessment of plasma inflammatory cytokines (high-sensitivity C-reactive protein, Interleukin 1β, Interleukin 6, and Tumor Necrosis Factor ɑ). Sleep characteristics included weekday and weekend sleep duration and midpoint (median of bed and wake time), as well as sleep variability (SD of sleep duration across 7 days) and social jetlag (midpoint difference from weekdays to weekends). At each time point, multivariable linear regression models were run with log inflammatory levels as the outcome and categories of sleep characteristics as predictors, while adjusting for potential confounders (specific to each model). Analyses were run unstratified and sex-stratified.

RESULTS

In the mid-to-late pubertal visit, weekday sleep duration was inversely associated with natural log hs-CRP after adjustment (Q4 vs Q1: β = -0.41, 95% Confidence Interval (CI) -0.81 to -0.01) and later sleep midpoint was positively associated with log hs-CRP (Q4 vs Q1: β = 0.55, 95% CI 0.13 to 0.97). Sleep duration variability was associated with higher IL-1β among boys, while in girls social jetlag was associated with higher IL-1β and weekend sleep duration was inversely associated with IL-6. At the late-to-post pubertal visit, there were few associations except for a positive association between weekday sleep duration and hs-CRP among boys (β = 0.60, 95% CI 0.04 to 1.16) and a non-linear positive association between social jetlag and hs-CRP among girls (β = 0.80, 95% CI 0.22 to 1.37 comparing 2 to 3 h of social jetlag vs <1 h).

CONCLUSION

Later timing, shorter duration, and inconsistency of sleep were related to higher levels of inflammatory biomarkers, but associations were more evident at the mid-to-late pubertal visit than the late-to-post pubertal visit.

Biomarkers of maternal lead exposure during pregnancy using micro-spatial child deciduous dentine measurements

BACKGROUND

Lead is a toxic chemical of public health concern, however limited biomarkers are able to reconstruct prior lead exposures in early-life when biospecimens are not collected and stored. Although child tooth dentine measurements accurately assess past child perinatal lead exposure, it has not been established if they reflect maternal exposure in pregnancy.

AIM

To assess the prenatal relationship between child tooth dentine and maternal blood lead measurements and to estimate maternal lead exposure during the 2nd and 3rd trimesters of pregnancy from weekly child dentine profiles.

METHODS

We measured early-life lead exposure in child tooth dentine and maternal blood from 419 child-mother dyads enrolled in the Programming Research in Obesity, Growth, Environment and Social Stress (PROGRESS) cohort. We employed the Super-Learner algorithm to determine the relationship of dentine lead data with maternal blood lead concentrations and to predict maternal lead from child dentine lead data in blinded analyses. We validated and quantified the bias of our results internally.

RESULTS

Mothers had moderate blood lead levels (trimesters: 2nd = 29.45 ug/L, 3rd = 31.78 ug/L). Trimester-averaged and weekly child dentine lead measurements were highly correlated with maternal blood levels in the corresponding trimesters. The predicted trimester-specific maternal lead levels were significantly correlated with actual measured blood values (trimesters: 2nd = 0.83; 3rd = 0.88). Biomarkers of maternal lead exposure discriminated women highly exposed to lead (>mean) with 85 % and 96 % specificity in the 2nd and 3rd trimesters, respectively, with 80 % sensitivity.

DISCUSSION

Weekly child dentine lead levels can serve as biomarkers of past child and maternal lead exposures during pregnancy.

Sleep Difficulties among Mexican Adolescents: Subjective and Objective Assessments of Sleep

OBJECTIVE/BACKGROUND

Self-reported sleep difficulties, such as insomnia symptoms, have been reported among adolescents. Yet, studies of their prevalence and correlates are scarce among Latin Americans. This study sought (1) to describe associations between sociodemographic and lifestyle factors with self-reported sleep difficulties and (2) to examine associations between self-reported sleep difficulties and actigraphy-based sleep.

PARTICIPANTS

Participants included 477 Mexican adolescents from the ELEMENT cohort.

METHODS

Over 7 days, self-reported sleep measures (hard time falling asleep, overall sleep difficulties, and specific types of sleep difficulties) were obtained from daily sleep diaries. Actigraphy-based sleep measures (duration, i.e. sleep onset to morning wake, midpoint, and fragmentation) were concurrently assessed using a wrist actigraph.

RESULTS

Mean (SD) age was 15.9 (2.2) years, and 53.5% were females. Mean (SD) sleep duration was 8.5 (1.2) h/night. Half reported a hard time falling asleep at least 3 days, and 25% had sleep difficulties at least 3 days over 7 days. The 3 types of sleep difficulties commonly reported among the entire cohort were insomnia/restlessness (29%), environmental (27%), and mental/emotional difficulties (19%). Female sex, smoking behavior, and socioeconomic indicators were among the most consistent factors associated with sleep difficulties. Subjective sleep difficulties were associated with shorter sleep duration (β = -20.8 [-35.3, -6.2] min), while subjective hard time falling asleep was associated with longer sleep duration (β = 11.3 [4.6, 27.2] min).

CONCLUSION

A high proportion of Mexican adolescents in the sample reported sleep difficulties. Findings demonstrate the importance of obtaining subjective and objective sleep measures for a more comprehensive assessment of adolescent sleep.

Diet Quality Scores and Cardiometabolic Risk Factors in Mexican Children and Adolescents: A Longitudinal Analysis

There is limited evidence for the effects of diet on cardiometabolic profiles during the pubertal transition. We collected repeated measures of diet quality and cardiometabolic risk factors among Mexican youth. This analysis included 574 offspring of the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) birth cohort followed up to three time points. Dietary Approaches to Stop Hypertension (DASH), alternate Mediterranean Diet (aMedDiet), and Children’s Dietary Inflammatory Index (C-DIITM) scores were computed from food frequency questionnaires. Higher DASH and aMedDiet scores reflect a higher diet quality, and lower C-DII scores reflect an anti-inflammatory diet. Cardiometabolic risk factors were lipid profile, glucose homeostasis, blood pressure, and waist circumference. Linear mixed models were used between quartiles of each diet score and outcomes. Compared to the first quartile, the fourth DASH quartile was inversely associated with log serum insulin (μIU/mL) [β = −0.19, p = 0.0034] and log-Homeostatic Model Assessment of Insulin Resistance [β = −0.25, p = 0.0008]. Additionally, log serum triglycerides (mg/dL) was linearly associated with aMedDiet score [β = −0.03, p = 0.0022]. Boys in the highest aMedDiet quartile had higher serum high-density lipoprotein cholesterol (mg/dL) [β = 4.13, p = 0.0034] compared to the reference quartile. Higher diet quality was associated with a better cardiometabolic profile among Mexican youth.

Early Gestational Exposure to High-Molecular-Weight Phthalates and Its Association with 48-Month-Old Children's Motor and Cognitive Scores

In utero phthalate exposure has been associated with neurodevelopmental disorders, nevertheless, trimester-specific susceptibility remains understudied. Our aim was to identify susceptible windows to the effects of gestational High-Molecular-Weight Phthalates (HMWP) exposure on 48 months’ neurodevelopment. We measured six HMWP metabolites (MEHP, MEHHP, MEOHP, MECPP, MBzP and MCPP) in urine samples collected during each trimester from women in the Early Life Exposure in Mexico to Environmental Toxicants (ELEMENT) cohort (n = 218). We assessed children’s motor (MS), cognitive (GCI) and memory (MeS) abilities using McCarthy Scales of Children’s Abilities (MSCA). We used linear regression models to examine associations between trimester-specific phthalate metabolites and MSCA scores, adjusted for sex, gestational age, breastfeeding, and maternal IQ. Although phthalate concentrations were similar across trimesters, first and second trimester phthalates were inversely associated with MS and GCI, with first trimester associations with MS being the strongest and statistically significant. Stronger associations were seen with MS and GCI among boys compared to girls, however interaction terms were not statistically significant. Our results suggest that early gestation is a sensitive window of exposure to HMWP for neurodevelopment, particularly in boys. Regulations on phthalate content in food as well as pregnancy consumption guidelines are necessary to protect future generations.