Next Article in Journal
Comprehensive Metabolic Profiling of Euphorbiasteroid in Rats by Integrating UPLC-Q/TOF-MS and NMR as Well as Microbial Biotransformation
Previous Article in Journal
Maternal High-Fat Diet Programs White and Brown Adipose Tissues In Vivo in Mice, with Different Metabolic and Microbiota Patterns in Obesity-Susceptible or Obesity-Resistant Offspring
Previous Article in Special Issue
Urinary Untargeted Metabolic Profile Differentiates Children with Autism from Their Unaffected Siblings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maternal Serum and Placental Metabolomes in Association with Prenatal Phthalate Exposure and Neurodevelopmental Outcomes in the MARBLES Cohort

by
Mariana Parenti
1,
Rebecca J. Schmidt
2,3,
Sally Ozonoff
3,4,
Hyeong-Moo Shin
5,
Daniel J. Tancredi
6,
Paula Krakowiak
2,
Irva Hertz-Picciotto
2,3,
Cheryl K. Walker
3,7 and
Carolyn M. Slupsky
1,8,*
1
Department of Nutrition, University of California, Davis, CA 95616, USA
2
Department of Public Health Sciences, University of California, Davis, CA 95616, USA
3
Medical Investigation of Neurodevelopmental Disorders (MIND) Institute, University of California, Davis, CA 95616, USA
4
Department of Psychiatry and Behavioral Sciences, University of California, Davis, CA 95616, USA
5
Department of Environmental Science, Baylor University, Waco, TX 76798, USA
6
Department of Pediatrics, School of Medicine, University of California, Davis, CA 95616, USA
7
Department of Obstetrics & Gynecology, School of Medicine, University of California, Davis, CA 95616, USA
8
Department of Food Science and Technology, University of California, Davis, CA 95616, USA
*
Author to whom correspondence should be addressed.
Submission received: 12 August 2022 / Revised: 26 August 2022 / Accepted: 27 August 2022 / Published: 2 September 2022
(This article belongs to the Special Issue Metabolomics of Autism Spectrum Disorder)

Abstract

:
Prenatal exposure to phthalates, a family of endocrine-disrupting plasticizers, is associated with disruption of maternal metabolism and impaired neurodevelopment. We investigated associations between prenatal phthalate exposure and alterations of both the maternal third trimester serum metabolome and the placental metabolome at birth, and associations of these with child neurodevelopmental outcomes using data and samples from the Markers of Autism Risk in Babies Learning Early Signs (MARBLES) cohort. The third trimester serum (n = 106) and placental (n = 132) metabolomes were investigated using 1H nuclear magnetic resonance spectroscopy. Children were assessed clinically for autism spectrum disorder (ASD) and cognitive development. Although none of the urinary phthalate metabolite concentrations were associated with maternal serum metabolites after adjustment for covariates, mixture analysis using quantile g-computation revealed alterations in placental metabolites with increasing concentrations of phthalate metabolites that included reduced concentrations of 2-hydoxybutyrate, carnitine, O-acetylcarnitine, glucitol, and N-acetylneuraminate. Child neurodevelopmental outcome was not associated with the third trimester serum metabolome, but it was correlated with the placental metabolome in male children only. Maternal phthalate exposure during pregnancy is associated with differences in the placental metabolome at delivery, and the placental metabolome is associated with neurodevelopmental outcomes in males in a cohort with high familial ASD risk.

1. Introduction

Phthalates are a family of plasticizers that can easily leach into the environment [1]. They are manufactured in products that humans interact with every day, such as vinyl building materials, food packaging, and personal care products, but are also found in house dust and food. Phthalates are endocrine-disrupting compounds that have been linked with differences in hormone levels in humans, including insulin, sex hormones, and thyroid hormones [2,3,4,5,6]. Moreover, phthalate exposure during gestation has been associated with impaired neurodevelopmental outcomes in humans [7]. Phthalate exposure—whether through transdermal absorption, inhalation, or diet—is usually measured based on urinary concentrations of phthalate monoesters and/or secondary oxidation products [8]. The fact that individual phthalate metabolites have been detected in greater than 89% of samples from 11,071 National Health and Nutrition Examination Survey (NHANES) participants suggests that phthalate exposure is ubiquitous in the U.S. [9]. In an analysis of chemical exposure among 268 pregnant women in NHANES, at least four phthalate metabolites were measured in each pregnant woman, with a median of nine phthalate metabolites detected in each subject [10]. More recent studies have detected some phthalate metabolites in 100% of samples from 205 pregnancies in California [11] and in 171 pregnancies in the United States and Puerto Rico [12]. Furthermore, the level of exposure tends to be higher in historically marginalized communities, including non-Hispanic Black and Hispanic women [13,14]. In addition to direct maternal exposure, phthalates and their metabolites have also been measured in amniotic fluid, cord blood, and breastmilk, indicating that exposure begins in utero and continues through early life [15,16,17]. In addition to endocrine disruption, phthalate exposure is known to induce oxidative stress in humans [18,19,20,21]. A delicate balance of oxidation is necessary for normal development [22], where low levels of reactive oxygen species act as signaling molecules. However, upsetting this balance has been associated with neurodevelopmental disorders including autism spectrum disorder (ASD) [23].
Urinary phthalate metabolites measured during gestation, including metabolites of di-2-ethylhexyl phthalate (DEHP), di-n-butyl phthalate (DnBP), and di-n-octyl phthalate (DnOP), have been associated with increased autistic traits [24,25,26]; however, some studies have shown no correlation between phthalate exposure and ASD [27,28]. There could be several reasons for this. For instance, ASD requires clinical training to assess, and many studies do not use gold-standard tools, such as the ASD Diagnostic Observation Schedule. Symptoms of autism are variable and can be difficult to diagnose, particularly in young children [29,30,31,32]. It could also be that phthalate exposure itself does not contribute to ASD and that it may serve as a proxy for other exposures and/or maternal behaviors related to ASD risk [33], such as nutrition [34], and other environmental exposures, such as air pollution, pesticides, and other plastic additives [35]. On the other hand, mixed evidence for the association between prenatal phthalate exposure and ASD may be related to how the exposure is measured. Investigating exposures to phthalate mixtures could be important to understanding the relationship between prenatal phthalate exposure and ASD risk [25]; so too could investigating the timing of the exposure, as sensitive windows of phthalate exposure for ASD risk may exist during gestation and postnatal life [25,26]. It has also been hypothesized that phthalate exposure could impact ASD risk through genetic susceptibility to oxidative stress since phthalates have been shown to elicit oxidative stress [36]. While it is difficult to investigate ASD as a binary variable given its multifactorial nature [30], a metabolomics approach allows us to investigate intermediate phenotypes associated with phthalate exposure that could be relevant to ASD development.
We have previously investigated the association between prenatal phthalate exposure and neurodevelopmental outcomes and determined that individual phthalate metabolites alone were not associated with ASD in this same cohort of children at increased risk for ASD as part of the Markers of Autism Risk in Babies, Learning Early Signs (MARBLES) study [27]. However, this study did not investigate the subtle changes in metabolism associated with increased susceptibility to oxidative stress. Metabolomics, the study of small molecules, provides an unbiased view of metabolic disturbances that represent the downstream effects of an exposure’s influence on gene expression and protein activity. Phthalate exposure has been shown to induce changes in serum metabolism, as well as placental methylation and transcription in humans [37,38,39].
Despite the heterogeneity in ASD presentation, previous research has identified metabolic disturbances in both children with ASD and their mothers during pregnancy [40,41,42,43]. Among children with ASD, evidence of altered mitochondrial dysfunction, oxidative stress, one-carbon metabolism, and amino acid metabolism have been reported [40]. Multiple metabolic subgroups among children with ASD have also been described [41]. Disruptions in metabolites involved in one-carbon metabolism and the transsulfuration pathway reflected in maternal blood have also been observed during gestation [42], and ASD risk has also been associated with disturbances in biosynthetic and metabolic pathways related to glycosphingolipids, phospholipids, bile acids, pyrimidines, N-glycans, and steroid hormones in maternal serum collected in mid-pregnancy [43]. This is the first study to the best of our knowledge to investigate the third trimester maternal serum metabolome using proton nuclear magnetic resonance (1H NMR)-based metabolomics in relation to ASD. Additionally, the cohort includes children who did not develop ASD, but who were otherwise non-typically developing. A recent study of maternal third trimester plasma collected as part of the MARBLES study has identified alterations in the prostaglandin pathway related to non-typical development [44].
In the present analysis, we aimed to investigate the associations between (a) maternal phthalate exposure during pregnancy and maternal serum metabolites in the third trimester and placental metabolites at birth, and (b) the serum and placental metabolomes and neurodevelopmental outcomes using 1H NMR-based metabolomic analyses of samples collected in the MARBLES cohort. Indeed, maternal serum and placenta metabolite concentrations could indicate the availability of substrates necessary for fetal growth and development. For instance, reduced concentrations of amino acids might reflect reduced availability for building proteins necessary for growth and development, for both the fetus and the placenta. Although phthalate exposure is associated with changes in placental growth and development [45,46], this is also the first study to the best of our knowledge to investigate the associations between phthalate exposure and the placental metabolome. Furthermore, with 14 phthalate metabolites from 8 phthalate diesters, we had the opportunity to investigate the effects of exposure to different mixtures of phthalates.

2. Materials and Methods

2.1. Study Population

This study included participants from the MARBLES cohort [47]. In this ongoing prospective observational study, families with one or more previous children with ASD are enrolled prior to or early on during a subsequent pregnancy. Children from these pregnancies are followed and their behavior and cognitive development are assessed clinically. At approximately 3 years of age, the children are classified into three groups of neurodevelopmental outcomes: typically developing, ASD, and non-typically developing. MARBLES enrollment began in 2006, with placenta samples used in this analysis collected as early as 2007. Maternal serum sample collection as part of MARBLES began in 2009. Maternal serum and placenta samples from pregnancies where maternal urinary phthalate metabolites were measured during the second and third trimesters were selected [27]. Of 200 analyzed serum samples from 119 pregnancies (1 twin pregnancy), 108 samples (1 from a twin pregnancy) were collected in the third trimester. Only 22 samples were collected in the first trimester and 69 were collected in the second trimester. Given that serum samples from the second and third trimesters were only available for 59 subjects, we opted to investigate the association between phthalate exposure and the serum metabolome using only third trimester samples. Furthermore, this allowed us to maintain a temporal relationship between earlier phthalate exposure and later serum metabolism since phthalate exposure was measured during the second and/or third trimester. While there were 108 third trimester serum samples (from 106 maternal subjects, including 1 twin pregnancy) from unique pregnancies, we selected the sample from the first available pregnancy in these analyses or the first twin. Of the 144 placenta samples connected to 133 mothers, we selected the placenta sample from the first available pregnancy to ensure independence between samples. One of these samples was excluded because it was connected to data from twins of different sexes and the placenta is fetal tissue. Of these 106 serum samples and 132 placenta samples, 78 mother-child pairs provided both serum and placenta samples, 54 pairs provided only a placenta sample, and 28 provided only a serum sample.

2.2. Child Neurodevelopmental Assessment

At 3 years of age, the children of interest were clinically assessed for ASD using the gold-standard ASD Diagnostic Observation Schedule (ADOS) [48]. Cognitive, language, and motor development were also assessed using the Mullen Scales of Early Learning (MSEL) [49]. Scores on both the ADOS and MSEL at approximately 36 months were used to classify neurodevelopmental outcomes as previously described [50]. Children with ADOS scores equal to or greater than the ASD cutoff were classified as ASD (serum: n = 27; placenta: n = 33). Children without ASD but who had ADOS scores within 3 points of the ASD cutoff or MSEL subdomain scores 1.5 to 2 standard deviations below the mean were classified as non-typically developing (Non-TD; serum: n = 14; placenta: n = 17). The remaining participants who did not meet the criteria for ASD or Non-TD were classified as typically developing (TD; serum: n = 66; placenta: n = 83).

2.3. Urinary Phthalate Metabolite Analysis

Phthalate metabolite concentrations were measured in pooled urine samples collected during the second and/or third trimesters as previously described [27] using online solid phase extraction coupled with high-performance liquid chromatography with isotope dilution-tandem mass spectrometry. The 14 phthalate metabolites measured included: monoethyl phthalate (MEP), mono-isobutyl phthalate (MiBP), monohydroxy-isobutyl phthalate (MHiBP), mono-n-butyl phthalate (MBP), monohydroxy-n-butyl phthalate (MHBP), monobenzyl phthalate (MBzP), mono(2-ethylhexyl) phthalate (MEHP), mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), mono(2-ethyl-5-oxohexyl) phthalate (MEOHP), mono(2-ethyl-5-carboxypentyl) phthalate (MECPP), mono(3-carboxypropyl) phthalate (MCPP), mono-isononyl phthalate (MNP), mono-carboxyisooctyl phthalate (MCOP), and mono-carboxyisononyl phthalate (MCNP). Concentrations were adjusted for specific gravity to correct for dilution and reported as µg/L. Machine-observed concentrations were used for values below the limit of detection. Urinary concentrations of DEHP metabolites were divided by their molar weights and summed to calculate the molar sum of DEHP metabolites as a marker of total DEHP exposure (ΣDEHP = MEHP + MEHHP + MEOHP + MECPP; µmol/L). The distribution of phthalate metabolite concentrations is presented in Table S1.

2.4. Serum Sample Preparation and 1H NMR Metabolomics Analysis

Serum samples collected during the third trimester were used in this study. Whole blood was centrifuged, and serum was stored at −80 °C prior to extraction. Samples were thawed on ice and filtered using Amicon Ultra-0.5 mL 3000 MW centrifugal filters (Millipore, Burlington, MA, USA) which were prepared by washing three times with ultra-pure water to remove glycerol. The resulting filtrate containing the water-soluble (polar) metabolites was collected and its volume was adjusted with ultra-pure water if an insufficient sample was collected. An internal standard (Chenomx, Edmonton, AB, USA) containing 5.0 mM 3-(trimethylsilyl)-1-propanesulfonic acid-d6 (DSS-d6), 0.2% NaN3, and 99.8% D2O was added. Each sample’s pH was adjusted to 6.8 ± 0.1 and 180 µL was loaded into a 3 mm NMR tube (Bruker, Billerica, MA, USA). Samples were stored at 4 °C until spectral acquisition within the same day.
1H-NMR spectra were collected using a Bruker Avance 600 MHz spectrometer (Bruker, Billerica, MA, USA) using the noesypr1d pulse sequence. Spectra were acquired at 25 °C, with water saturation of 2.5 s during the prescan delay, a mixing time of 100 ms, a sweep width of 12 ppm, an acquisition time of 2.5 s, 8 dummy scans, and 32 transients. All spectra were zero-filled to 128 K data points and Fourier transformed with a 0.5 Hz line broadening applied. Spectra were manually phase- and baseline-corrected using Chenomx NMR Processor (v. 8.1, Chenomx, Edmonton, AB, USA), and metabolite concentrations were quantified using Chenomx Profiler (v. 8.1). The quantification relies on the internal standard (DSS-d6) to determine each metabolite’s concentration using a library of compound spectral signatures which allows for the absolute quantification of many compounds within a spectrum.

2.5. Placenta Sample Preparation and 1H NMR Metabolomics Analysis

Placentas were collected at delivery and full-thickness tissue sections were stored at −80 °C. All samples were partially-thawed and a 6 mm biopsy punch was used to aliquot samples for metabolic analyses which were again stored at −80 °C. Placenta samples were then cryoground to a uniform fine powder using mortars, pestles, and liquid nitrogen. The polar metabolites from each sample were extracted in a two-step CHCl3:MeOH:H2O extraction [51]. Approximately 80 mg of each cryoground sample was used and the upper layer containing polar metabolites was collected, measured, and frozen before being dried using a miVac concentrator system (Genevac, Warminster, PA, USA). The dried samples were stored at −80 °C until preparation for 1H NMR spectroscopy. Samples were reconstituted in 10 mM potassium phosphate buffer. The internal standard DSS-d6 was added and the samples were prepared for spectral acquisition and metabolite identification and quantification as described above for serum samples.

2.6. Statistical Analysis

All statistical analyses were conducted using the R statistical language (R Foundation for Statistical Computing, Vienna, Austria) with RStudio. Demographic differences between neurodevelopmental outcome groups were assessed using chi-squared tests for categorical variables, Fisher’s Exact test for categorical variables with empty cells, and Kruskal–Wallis tests for continuous variables. Serum metabolites were corrected for dilution and reported in µmol/L. Placenta metabolites were corrected and are reported in nmol/g tissue. Phthalate, placenta, and serum metabolites were log10-transformed to improve skewed distributions.
We identified and quantified 59 serum metabolites and 62 placental metabolites. However, metabolites that were identified in samples but may have been introduced during sample preparation (serum: ethanol, glycerol, and isopropanol; placenta: ethanol, isopropanol, and methanol) and metabolites that were not quantified in at least 20% of samples (serum: acetaminophen, ethyl-β-D-glucoronide, fructose, inosine, N-phenylacetylglycine, propylene glycol, trimethylamine N-oxide, and valproate; placenta: acetone, ascorbate, propylene glycol, and urea) were excluded from these analyses. Thus, we included 48 serum metabolites and 54 placental metabolites in this analysis, in conjunction with 14 urinary phthalate metabolites previously measured and described elsewhere [27]. Correlations between placental metabolite concentrations and serum metabolite concentrations are presented in Figure S1.
Univariate normality was assessed using the Shapiro–Wilks test and multivariate normality with Mardia’s test (package MNV). We used two approaches in evaluating the effects of phthalate exposure on the serum and placental metabolomes. In a substance-by-substance approach, correlations between urinary phthalate metabolites and either placenta metabolites or serum metabolites were evaluated using Spearman’s rank correlation coefficient and corrected for false discovery rate (FDR) using the Benjamini–Hochberg procedure [52]. Estimates of the total and direct effects of each phthalate metabolite on each placenta and serum metabolite were evaluated independently and adjusted for covariates using multiple linear regression (MLR) and FDR-corrected. In a mixture approach, the total and direct effects of phthalate mixture exposure on each placenta and serum metabolite were evaluated independently and adjusted for covariates using quantile g-computation (implemented using the package qgcomp) and FDR-corrected [53]. Quantile g-computation is useful for investigating the effect of an exposure mixture by estimating the effect of increasing all exposures by one quantile, adjusted for covariates [53]. In this analysis, phthalate exposure was split into quartiles, so each model estimates the change in the log10-transformed metabolite concentration after increasing all phthalate metabolites by 1 quartile (ψ).
The associations between neurodevelopmental outcome and either the serum metabolome or placental metabolome adjusted for covariates were assessed using permutation multivariate analysis of variance (PERMANOVA). The association between neurodevelopmental outcome and each serum or placenta metabolite was assessed using single-response PERMANOVA adjusted for covariates and FDR-corrected [54]. PERMANOVA was implemented using the package vegan [55]. We considered p-values (for multivariate analyses) less than 0.05 and FDR-corrected p-values less than 0.1 to be statistically significant.
Because the associations between birth year and phthalate exposure were closely tied in the placenta, we tested the relationship between birth year and placental metabolites in those subjects with low DEHP exposure by conducting MLR for associations between birth year and log10-transformed metabolite concentration after adjustment for gestational age at delivery, maternal metabolic condition, homeownership status, maternal race/ethnicity, and maternal education in subjects with a ΣDEHP less than the median (n = 66). As normality assumptions were not met, significance was tested by permuting the residuals under the reduced model (containing only the covariates) and regressing them on birth year [56].
Individual directed acyclic graphs (DAGs) were constructed in R (package dagitty) [57] to identify potential confounders. These potential covariates included maternal prenatal vitamin use during the first month of pregnancy, birth year, birth weight, maternal age, maternal educational attainment, maternal race/ethnicity, gestational age at sampling and/or at gestation, method of delivery, maternal metabolic condition, and child’s sex. Information about self-reported maternal prenatal vitamin use during the 6 months prior and throughout pregnancy was collected in the MARBLES cohort. Maternal prenatal vitamin use in the first month of pregnancy was used here because was inversely associated with ASD risk in the MARBLES cohort [58]. The maternal metabolic condition was coded to integrate maternal pre-pregnancy body mass index (BMI), hypertensive disorders (including preeclampsia), and diabetes during pregnancy with five levels: BMI ≤ 25 and no metabolic conditions [reference], 25 < BMI < 30 and no metabolic conditions, BMI ≥ 30 and no other metabolic conditions, any hypertensive disorder (without any diabetes) at any BMI, and diabetes at any BMI [59]. From the DAGs (Figure S2), we identified adjustment sets to reduce bias in the estimates of associations between (a) urinary phthalate metabolites and serum metabolites and (b) serum metabolites and neurodevelopmental outcomes. The minimal adjustment set to reduce bias in the estimates of the total effect of phthalate exposure on each serum metabolite included birth year, maternal education, maternal race/ethnicity, and homeownership status. We also identified the maternal metabolic condition as a potential mediator on the path between prenatal phthalate exposure and serum metabolite concentrations. The minimal adjustment set to reduce bias in the estimates of the direct effect of phthalate exposure on each serum metabolite included birth year, maternal race/ethnicity, homeownership status, and maternal metabolic condition. The minimal adjustment set to reduce bias in the estimates of the total effect of the serum metabolome on neurodevelopmental outcome included birth year, maternal race/ethnicity, homeownership status, and maternal metabolic condition. Additionally, fasted time and gestational age at sample collection were included in all adjusted models.
From the DAG (Figure S3), we identified adjustment sets to reduce bias in the estimates of associations between (a) urinary phthalate metabolites and placental metabolites and (b) placental metabolites and neurodevelopmental outcomes. The minimal adjustment set to reduce bias in the estimates of the total effect of phthalate exposure on each placental metabolite included birth year, maternal education, maternal race/ethnicity, and homeownership status. We also identified maternal metabolic condition, delivery mode, and gestational age at delivery as potential mediators on the path between prenatal phthalate exposure and placental metabolite concentrations. The minimal adjustment set to reduce bias in the estimates of the direct effect of phthalate exposure on each placental metabolite included birth year, maternal education, maternal race/ethnicity, homeownership status, maternal metabolic condition, delivery mode, and gestational age at delivery. The minimal adjustment set to reduce bias in the estimates of the total effect of the placental metabolome on neurodevelopmental outcome included birth year, birth weight, maternal race/ethnicity, homeownership status, maternal metabolic condition, maternal age, prenatal vitamin use in the first month of pregnancy, delivery mode, and gestational age at delivery.

3. Results

Table 1 summarizes the demographic characteristics of the mothers and children studied grouped by the child’s neurodevelopmental outcome. No significant differences were observed in the subpopulation in which the serum metabolome was examined. However, among subjects providing placenta samples, the child’s neurodevelopmental outcome was significantly associated with prenatal vitamin use in the first month of pregnancy and gestational age at delivery (p < 0.05).

3.1. Associations between Phthalate Metabolites and Serum Metabolites

Spearman correlations were used to explore the associations between urinary phthalate metabolites and serum metabolites. After FDR correction, no associations remained significant (Figure 1). MLR was also conducted to estimate the total and direct effects of individual urinary phthalate metabolites on individual serum metabolites after adjustment for covariates. The minimal adjustment set to reduce bias in the estimates of the total effect of phthalate exposure on each serum metabolite included birth year, maternal education, maternal race/ethnicity, and homeownership status. The minimal adjustment set to reduce bias in the estimates of the direct effect of phthalate exposure on each serum metabolite included birth year, maternal race/ethnicity, homeownership status, and maternal metabolic condition. Additionally, fasted time and gestational age at sample collection were included in all adjusted models. We observed no significant total nor direct estimated effects of individual phthalate metabolites on individual serum metabolites following FDR correction (Table S2). Similarly, mixture analysis using quantile g-computation did not reveal significant total or direct estimated effects of phthalates on serum metabolites after FDR correction (Table S3).

3.2. Associations between Serum Metabolites and Neurodevelopmental Outcomes

PERMANOVA was used to explore the association between neurodevelopmental outcome and the serum metabolome, and no association was observed (R2 = 0.0165, p = 0.6129 under 9999 permutations), nor was an association observed after adjustment for covariates, which included birth year, maternal race/ethnicity, homeownership status, maternal metabolic condition, fasted time, and gestational age at sample collection (Table S4, R2 = 0.0144, p = 0.7125 under 9999 permutations). Furthermore, in single-response PERMANOVAs exploring the association between neurodevelopmental outcome and individual serum metabolites, only cystine was elevated in mothers carrying a child later diagnosed with ASD. Post hoc pairwise comparisons revealed that cystine concentrations were elevated in the ASD group compared to the TD group (p = 0.0357). However, this association did not survive FDR correction (Table S5). Stratifying analysis on the child of interest’s sex similarly showed null results across both sexes (Table S4).

3.3. Associations between Phthalate Metabolites and Placenta Metabolites

Spearman correlations were used to explore the associations between urinary phthalate metabolites and placenta metabolites. After FDR correction, only the DEHP metabolites MEHP, MEHHP, MEOHP, and MECPP remained significantly associated with multiple placenta metabolites (Figure 1). They were weakly to moderately positively correlated with cystine, methionine, pyroglutamate, and tryptophan and weakly to moderately negatively correlated with betaine, glutathione, inosine, NAD+, O-acetylcarnitine, O-phosphocholine, and O-phosphoethanolamine. Additionally, MEP was also negatively associated with inosine. We also used MLR to estimate the total and direct effects of individual urinary phthalate metabolites on individual placental metabolites after adjustment for covariates. The minimal adjustment set to reduce bias in the estimates of the total effect of phthalate exposure on each placental metabolite included birth year, maternal education, maternal race/ethnicity, and homeownership status. The minimal adjustment set to reduce bias in the estimates of the direct effect of phthalate exposure on each placental metabolite included birth year, maternal education, maternal race/ethnicity, homeownership status, maternal metabolic condition, delivery mode, and gestational age at delivery. However, we did not observe significant estimates of the total or direct effects of individual phthalate metabolites on individual placental metabolites after FDR correction (Table S6).
We employed mixture analysis using quantile g-computation to investigate the total and direct estimated effects of phthalates on placental metabolites using the same covariates (Table 2). After FDR correction, we observed a significant, negative total estimated effect of phthalates on placental glucitol (also known as sorbitol), carnitine, O-acetylcarnitine, N-acetylneuraminate (also known as sialic acid), and 2-hydroxybutyrate (Table 2). We observed significant, negative direct effect estimates of the phthalate mixture on placental glucitol, carnitine, O-acetylcarnitine, N-acetylneuraminate, and 2-hydroxybutyrate (Table 2).

3.4. Associations between Placenta Metabolites and Birth Year

Birth year and phthalate exposure, especially to DEHP, are inversely related in this cohort. Since adjustment for birth year attenuated the associations between phthalate exposure and placental metabolites, we tested the association between birth year and placental metabolites in those subjects with low DEHP exposure. In subjects with a molar sum of DEHP metabolites less than or equal to the median (n = 66), we conducted MLR to test the association between birth year and log10-transformed metabolite concentration after adjustment for gestational age at delivery, homeownership status, maternal metabolic condition, maternal race/ethnicity, and maternal education (Table 3 and Table S7). After FDR correction, creatine, glutathione, inosine, NAD+, O-phosphocholine, and O-phosphoethanolamine were positively associated with birth year, while arginine, cystine, isoleucine, leucine, lysine, methionine, ornithine, phenylalanine, pyroglutamate, serine, tryptophan, tyrosine, and uracil were significantly negatively associated with the birth year (Table 3, FDR p < 0.1). Additional associations with the birth year for metabolites with FDR p > 0.1 are reported in the Supplementary Materials (Table S7).

3.5. Associations between Placenta Metabolites and Neurodevelopmental Outcomes

PERMANOVA was used to explore the association between neurodevelopmental outcomes and the placental metabolome. In an unadjusted analysis, a small correlation was observed (R2 = 0.0289, p = 0.0303 under 9999 permutations). However, no association was observed after adjustment for covariates, including birth year, birth weight, maternal race/ethnicity, homeownership status, maternal metabolic condition, maternal age, prenatal vitamin use in the first month of pregnancy, delivery mode, and gestational age at delivery (Table 4, R2 = 0.0201, p = 0.1270 under 9999 permutations). When stratifying the analysis by the child’s sex, we observed that the placental metabolome was associated with neurodevelopmental outcome after adjustment for covariates in males (Table 4, R2 = 0.0440, p = 0.0451 under 9999 permutations). Post hoc pairwise PERMANOVAs of the placental metabolomes from male children corrected for multiple testing showed that the ASD and TD groups tended to differ (R2 = 0.0361, p = 0.0705), but that the Non-TD group did not differ significantly from the others (p < 0.1).
We also conducted single-response PERMANOVAs to investigate the association between placenta metabolites and neurodevelopmental outcomes in all children and stratified by sex. Overall, concentrations of niacinamide, taurine, and uridine differed by the neurodevelopmental outcome after adjustment for covariates (Table S8). While the post hoc pairwise testing revealed no difference in niacinamide concentrations between groups, taurine and uridine concentrations were lower in the ASD group than in the TD group (taurine, p = 0.0567; uridine, p = 0.0255). However, these associations did not survive FDR correction (Table S8). In female children, concentrations of choline, glycerol, and taurine were associated with different neurodevelopmental outcomes after adjustment for covariates (Table S8). In the post hoc pairwise comparisons, we found that choline and glycerol concentrations were lower in ASD compared to TD (choline: p = 0.0447; glycerol, p = 0.0492). Taurine concentrations were lower in ASD compared to both the TD (p = 0.0228) and Non-TD groups (p = 0.0419). However, these associations did not survive FDR correction (Table S8). In male children, concentrations of alanine, cystine, ethanolamine, fumarate, hypoxanthine, methionine, NAD+, pyroglutamate, and uridine differed by the neurodevelopmental outcome after adjustment for covariates (Table S8). No pairwise group differences in ethanolamine, hypoxanthine, or pyroglutamate concentrations were observed in the post hoc pairwise comparisons. Alanine, cystine, fumarate, methionine and uridine concentrations were lower in the ASD group than in TD group (alanine, p = 0.0702; cystine, p = 0.0645; fumarate, p = 0.0222; methionine, p = 0.0978; and uridine, p = 0.0096), while NAD+ concentrations were higher in the ASD group than the TD group (p = 0.0045). However, these associations also did not survive FDR correction (Table S8).

4. Discussion

In this study, we investigated the associations between prenatal phthalate exposure and the maternal serum metabolome during the third trimester and the placental metabolome at birth. In the placenta, we observed that phthalates were negatively associated with 2-hydroxybutyrate, carnitine, glucitol, O-acetylcarnitine, and N-acetylneuraminate concentrations. We also observed that neurodevelopmental outcome was associated with the placental metabolome in male children, though associations with individual metabolites did not survive FDR correction. However, we did not observe that the maternal serum metabolome was related to either phthalate exposure or neurodevelopmental outcome.

4.1. Associations with Phthalate Exposure

In the U.S. generally, and this cohort specifically, phthalate exposure patterns have changed over the course of the study (2007–2014) [9,11]. In the MARBLES cohort, exposure to the phthalate diesters butyl benzyl phthalate, dibutyl phthalate, DEHP, and diethyl phthalate have decreased, while exposures to di-isobutyl phthalate, di-isononyl phthalate, and DnOP have increased [11]. Thus, there are differing profiles of phthalate exposure over the course of the study. While many studies report only individual associations or investigate the effects of phthalate exposure substance-by-substance, phthalate exposure generally occurs as a combination of several phthalates together with other environmental toxicants. These profiles are important to consider because chemical mixtures may exert different biological effects compared to each component substance alone [60]. Furthermore, the pattern of urinary phthalate metabolites may point to differences in phthalate metabolism and excretion between individuals. In the case of DEHP, a higher ratio of its metabolically active monoester MEHP to its secondary oxidized metabolites (MEHHP, MEOPH, and MECPP) is associated with greater physiological effects [61]. Thus, analysis of the mixture provides information about both environmental phthalate exposure and individual phthalate metabolism relevant to the physiological effects.
We took two approaches to analyze the relationship between phthalate exposure and the placental metabolome. First, in a substance approach, we did not observe significant effects of individual phthalate metabolites on individual placental metabolites after adjusting for covariates. However, when we applied mixture analysis, we found that phthalates were negatively associated with 2-hydroxybutyrate, carnitine, glucitol, O-acetylcarnitine, and N-acetylneuraminate concentrations after adjustment for covariates.
Carnitine and acylcarnitines are classically involved in lipid transport into the mitochondria and subsequent metabolism [62]. Exposure to DEHP and its metabolites has been associated with elevated lipid concentrations in rat trophoblasts [63] and altered placental lipid metabolism [64]. Thus, lower carnitine and O-acetylcarnitine concentrations with increasing exposure to the phthalate mixture could reflect impaired placental lipid transport and might suggest, in conjunction with evidence of oxidative stress, impaired mitochondrial function. Moreover, carnitine increases cell proliferation and amino acid transporter expression and abundance through insulin-like growth factor I signaling in rat trophoblasts [65]. Carnitine is essential to fetal development [66]. During pregnancy, maternal circulating carnitine levels drop in part due to placental transfer to the fetus via carnitine/organic cation transporter 2 (OCTN2) [67]. The placenta cannot synthesize carnitine, so OCTN2-mediated transport is critical [68]. However, under hypoxic conditions such as preeclampsia, placental OCTN2 abundance and carnitine concentrations decrease [69]. Prenatal phthalate exposure has been associated with an increased risk of pregnancy-induced hypertensive disorders, including preeclampsia [70,71]. The lower placental concentrations of O-acetylcarnine and carnitine with increasing exposure to phthalates could point to altered placental function.
Carnitine can also act as an osmolyte [72]. Glucitol is a polyol and an osmolyte produced by the placenta [73]. In early development, the polyol pathway is an important source of energy in a low-oxygen environment [74]. However, placental glucitol concentrations increase with gestational age [75] after the placenta becomes the site of fetal gas exchange. Similar to other highly metabolic organs, the placenta has relatively high osmolarity [75]. While we did not observe changes in the concentrations of other osmolytes (including taurine and myo-inositol) in response to the phthalate mixture, changes in glucitol concentrations could indicate changes in metabolism.
N-Acetylneuraminate, also known as sialic acid, has a role in neurodevelopment, immune function, and angiogenesis. Sialic acid is an essential nutrient in neurodevelopment [76]. In rats, maternal prenatal supplementation with sialic acid improved offspring cognition and resulted in elevated concentrations of sialic acid in the brain [77]. However, we did not observe that placental sialic acid was associated with ASD or Non-TD risk. In the placenta, sialic acid is essential starting in early pregnancy for protecting the conceptus from detection and destruction by the maternal innate immune response [78]. Decreased levels of placental sialyated oligosaccharides have been reported in both spontaneous and recurrent miscarriages [79] and preeclampsia [80]. Furthermore, sialic acid has been implicated in regulating angiogenesis. The angiogenic vascular endothelial growth factor must be sialyated in order to bind to its receptor [81], which is expressed in the placenta [68]. In the placenta, angiogenesis continues throughout pregnancy and is critical for increasing blood flow (and thus nutrients) to the fetus in the third trimester [82]. Phthalates have been associated with markers of impaired placental angiogenesis in human and animal studies [83]. In our study, sialic acid concentrations decreased as exposure to phthalates increased, suggesting that phthalate exposure may be related to impaired placental function.
2-Hydroxybutyrate is produced through methionine and threonine catabolism, as well as through the production of cysteine via the transsulfuration pathway [84]. Under conditions of oxidative stress, glutathione biosynthesis increases, and cysteine may become limiting. Thus, as a by-product of cysteine production through glutamate-cysteine ligase, 2-hydroxybutyrate is a sign of glutathione synthesis and oxidative stress [85]. DEHP has been shown to reduce glutathione concentrations in mouse kidneys through the depletion of transcription factor nuclear factor erythroid 2-related factor 2 (Nrf2) and glutamate-cysteine ligase abundance [86]. Nrf2 has been identified as a potential target to improve antioxidant defense in exposure to DEHP and other phthalates [87]. In our analysis, placental glutathione concentrations were negatively associated with increased phthalate concentrations prior to FDR correction. Thus, the negative association between the phthalates and 2-hydroxybutyrate may reflect reduced glutathione biosynthesis.
While concentrations of these metabolites were not significantly or strongly associated with birth year, we identified birth year as a potential confounder in our analysis. Our samples were stored at −80 °C, though some samples had been stored for over a decade before extraction. Cold storage, even at −80 °C, is associated with plasma sample degradation over time [88]. While the human EDTA plasma metabolome is relatively stable in samples stored for up to 7 years, amino acids, lipids, and energy metabolites have been shown to be sensitive to additional storage time at −80 °C [88], and thus we adjusted for the birth year in our serum analyses. While similar analyses have not been carried out in placenta samples, we observed similar effects when we reframed the birth year as a marker of storage time. With an earlier birth year, and thus increased storage time, we observed decreasing levels of the antioxidants inosine and glutathione, energy metabolite NAD+, and O-phosphocholine and O-phosphoethanolamine, as well as increasing levels of the glutathione breakdown-product pyroglutamate and cystine, the oxidized form of cysteine. We also observed increasing concentrations of many amino acids and uracil, potentially indicating the breakdown of macromolecules. It is also possible that over the course of the study, sample processing time improved. Previously, time delay in sample processing was shown to alter metabolite concentrations in placenta samples [89,90]. It is possible that these differences in metabolite concentrations are a sign of oxidation or degradation of the samples during processing and/or storage. So, while the birth year is associated with changes in phthalate exposure, it is also associated with changes in the placental metabolome. Thus, especially for metabolites related to oxidative stress, the potential effects of phthalate exposure could not be disentangled from the potential effects of birth year and long-term storage.
In this analysis, we did not observe an effect of phthalate exposure on the maternal third trimester serum metabolome. In our analysis of the serum metabolome, samples were collected throughout the third trimester. Over the course of pregnancy, serum concentrations of some amino acids increase, such as glutamine and methionine, while others decrease, such as leucine and valine [91]. Indeed, in this cohort, we observed a significant effect of gestational age at sample collection on the serum metabolome. Furthermore, participants did not fast before the serum was collected which contributed to variation in the serum metabolome. While we were able to adjust for time since the last meal or snack, we were not able to account for the potential effects of the size or nutrient composition of the most recent intake.

4.2. Associations with Neurodevelopmental Outcome

While phthalate exposures have been changing, ASD incidence has increased, due in part to increased awareness and evolving diagnostic criteria. However, these changes do not completely explain the increase [92]. ASD is highly heritable, but approximately 20% of the risk is environmental [93]. In a previous analysis of the effect of phthalate exposure on the neurodevelopmental outcome, phthalate exposure was not associated with ASD but some phthalate metabolites were associated with Non-TD [27]. We aimed to determine if differences in maternal serum metabolism were associated with neurodevelopmental outcomes and did not observe any associations.
Since the MARBLES cohort consists of younger siblings of children with ASD, the risk of ASD in this cohort is enriched [94]. This might partly explain why we did not observe significant associations between the maternal serum metabolome (or the individual metabolites that constitute it) and neurodevelopmental outcomes. Previous analyses of the maternal serum or plasma polar metabolome have identified metabolites and metabolic pathways associated with their child’s neurodevelopment, but these studies were conducted within the general population or explicitly included mothers with lower risk. In a study comparing plasma from pregnant women in the MARBLES cohort (at high risk of having a subsequent child with ASD) and women at low risk of having a subsequent child with ASD, metabolites related to 1-carbon metabolism and the transsulfuration pathways were useful in predicting the maternal risk group [42]. However, within the high-risk group from the MARBLES cohort, the maternal plasma metabolome was not predictive of ASD development in the child [42]. Similarly, metabolites related to one-carbon metabolism and the transsulfuration pathway differed in blood collected from non-pregnant mothers of children with and without ASD [95]. In an analysis of maternal serum samples collected around the 16th week of gestation as part of the California Prenatal Screening Program, metabolites related to glycosphingolipid, pyrimidine, bile acid, and steroid hormone metabolism differed between mothers of children with ASD and controls [43]. A more recent analysis of maternal third trimester plasma within the MARBLES cohort identified several lipid metabolites—related to fatty acid biosynthesis and essential fatty acid metabolism—negatively associated with Non-TD risk after FDR correction [44]. However, associations between polar metabolites positively associated with ASD risk and negatively associated with Non-TD risk did not survive FDR correction [44]. We speculate that in this high-familial risk cohort, maternal metabolic phenotype might be more similar than among mothers with low risk or that genetic contributions to ASD risk may outweigh metabolic contributions.
We also investigated associations between the placental metabolome and neurodevelopmental outcomes using PERMANOVA. The placenta, through its roles from nutrient transport to neurotransmitter biosynthesis, is known to influence fetal brain development [96]. We observed that the placental metabolome and neurodevelopmental outcome were associated in male children when we stratified by sex, though individual metabolite concentrations did not differ significantly between groups after FDR correction. However, males are at increased risk of neurodevelopmental disorders, including ASD [97]. Several theories to explain this sex bias have been put forward. One of these is the female protective effect theory, which holds that females have a higher threshold for developing a neurodevelopmental disorder and thus require more risks or exposures than males [98]. The biological underpinnings of a protective effect are yet unclear, but roles for the sex chromosomes and sex hormones have been proposed [98]. Fetal sex also plays a role in the development and function of the placenta, which shares its sex with the fetus as fetal tissue [99]. Males are at higher risk of a variety of adverse pregnancy and birth outcomes, including intrauterine growth restriction, preterm birth, and gestational diabetes [100,101]. Furthermore, males tend to experience more adverse effects during pregnancy than do similarly exposed females [101]. Sex differences in placental development have been implicated in sex differences in neurodevelopment through the placenta-brain axis [96]. Thus, the association between neurodevelopmental outcome and the placental metabolome only in males could be tied to sex differences in how the placenta processes exposures.

4.3. Limitations

These findings are from an enriched risk cohort for ASD, so they may not be generalizable to the broader population without a family history of ASD. Given that ASD is highly heritable, and this cohort is enriched for this heritable risk by using younger siblings, differences in maternal serum metabolism related to neurodevelopment might not be observed in this dataset. Additionally, while we considered a variety of potential confounding factors in this study, we could not rule out potential confounding by unmeasured factors, including other environmental toxicants. Indeed, pregnant women in the U.S. are exposed to a cocktail of chemical analytes from multiple chemical classes, including phthalates, pesticides, polybrominated diphenyl ethers, and toxic metals [10,12]. It is also possible that the phthalate measurement did not fully capture phthalate exposure over the course of pregnancy. In the MARBLES cohort, within-subject variability for urinary phthalate metabolites collected across the second and third trimester varied by subject and metabolite [102]. Furthermore, phthalate exposures are known to change over the course of pregnancy [103,104]. Finally, the limitations related to sample processing, storage, and other potential sources of variability might obscure potential associations between the serum and placental metabolomes and neurodevelopmental outcomes. In combination with the relatively small sample size, these sources of variability may have reduced our ability to detect subtle differences in metabolite concentrations that may be detected in larger studies.

5. Conclusions

Prenatal exposure to phthalate mixtures is related to changes in placental metabolism. Larger studies will be necessary to understand the metabolic and developmental effects of phthalate mixtures as sources and patterns of exposure are poised to change. We observed an association between placental metabolites and phthalate metabolite concentrations, despite the important effect of the birth year on the placental metabolome. Further, we observed that the placental metabolome was associated with neurodevelopmental outcomes in males. Additional studies investigating the associations between placental metabolism, prenatal phthalate exposures, and neurodevelopmental outcomes are warranted.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo12090829/s1, Table S1: Distribution of phthalate metabolite concentrations, Table S2: Estimates of the total and direct effects of each phthalate metabolite on each serum metabolite, Table S3: Estimates of the total and direct effects the phthalate mixture on each serum metabolite, Table S4: Associations between the serum metabolome and neurodevelopmental outcome stratified by sex, Table S5: Associations between serum metabolites and neurodevelopmental outcome, Table S6: Estimates of the total and direct effects of each phthalate metabolite on each placenta metabolite, Table S7: Associations between placenta metabolites and birth year, Table S8: Associations between placenta metabolites and neurodevelopmental outcome overall and stratified by sex, Figure S1: Correlation heatmap between placental and serum metabolites, Figure S2: Directed acyclic graphs (DAGs) for the association between phthalate exposure, serum metabolome, and neurodevelopmental outcome, Figure S3: Directed acyclic graph (DAG) for the association between phthalate exposure, placenta metabolome, and neurodevelopmental outcome.

Author Contributions

Conceptualization, C.K.W., D.J.T. and C.M.S.; Methodology, C.M.S.; Software, C.M.S.; Formal Analysis, M.P. and H.-M.S.; Investigation, M.P.; Resources, C.M.S.; Data Curation, R.J.S., S.O., P.K., I.H.-P. and M.P.; Writing—Original Draft Preparation, M.P.; Writing—Review and Editing, R.J.S., S.O., D.J.T., I.H.-P., C.K.W., C.M.S., H.-M.S. and P.K.; Visualization, M.P.; Supervision, C.M.S.; Project Administration, R.J.S., I.H.-P., C.K.W. and C.M.S.; Funding Acquisition, C.K.W., C.M.S., R.J.S., I.H.-P. and H.-M.S. All authors have read and agreed to the published version of the manuscript.

Funding

The project described was supported by Grant Number T32 ES007059, R21ES028129, R21ES025551, R01ES020392, R01ES028089, R24ES028533, P01ES11269, and R01ES015359 from NIH NIEHS. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS, NICHD, or NIH. This work was also supported by US EPA STAR grants R829388 and R833292 and also supported by the USDA National Institute of Food and Agriculture Hatch project 1021411. This work was also supported by Simons Foundation Autism Research Initiative (SFARI) 863967 and University of California Davis Medical Investigation of Neurodevelopmental Disorder (MIND) Institute Intellectual and Developmental Disabilities Research Center (IDDRC) U54-HD079125 from the NIH NICHD.

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of the University of California, Davis (protocol code 225645; date of initial approval 6/8/2006, with annual reapproval since).

Informed Consent Statement

Informed consent was obtained from all maternal subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the need to protect participant privacy.

Acknowledgments

The authors express their sincere appreciation for the families participating in the MARBLES project. The authors would also like to acknowledge Antonia Calafat and the late Xiaoyun Ye (Centers for Disease Control and Prevention, CDC) for the measurement of urinary phthalate metabolites.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Rowdhwal, S.S.S.; Chen, J. Toxic Effects of Di-2-Ethylhexyl Phthalate: An Overview. Bio. Med. Res. Int. 2018, 2018, 1750368. [Google Scholar] [CrossRef] [PubMed]
  2. Martino-Andrade, A.J.; Liu, F.; Sathyanarayana, S.; Barrett, E.S.; Redmon, J.B.; Nguyen, R.H.N.; Levine, H.; Swan, S.H.; the TIDES Study Team. Timing of Prenatal Phthalate Exposure in Relation to Genital Endpoints in Male Newborns. Andrology 2016, 4, 585–593. [Google Scholar] [CrossRef] [PubMed]
  3. Dales, R.E.; Kauri, L.M.; Cakmak, S. The Associations between Phthalate Exposure and Insulin Resistance, β-Cell Function and Blood Glucose Control in a Population-Based Sample. Sci. Total Environ. 2018, 612, 1287–1292. [Google Scholar] [CrossRef] [PubMed]
  4. Høyer, B.B.; Lenters, V.; Giwercman, A.; Jönsson, B.A.G.; Toft, G.; Hougaard, K.S.; Bonde, J.P.E.; Specht, I.O. Impact of Di-2-Ethylhexyl Phthalate Metabolites on Male Reproductive Function: A Systematic Review of Human Evidence. Curr. Environ. Health Rep. 2018, 5, 20–33. [Google Scholar] [CrossRef]
  5. Kim, M.J.; Moon, S.; Oh, B.-C.; Jung, D.; Choi, K.; Park, Y.J. Association Between Diethylhexyl Phthalate Exposure and Thyroid Function: A Meta-Analysis. Thyroid 2019, 29, 183–192. [Google Scholar] [CrossRef]
  6. Shoshtari-Yeganeh, B.; Zarean, M.; Mansourian, M.; Riahi, R.; Poursafa, P.; Teiri, H.; Rafiei, N.; Dehdashti, B.; Kelishadi, R. Systematic Review and Meta-Analysis on the Association between Phthalates Exposure and Insulin Resistance. Environ. Sci. Pollut. Res. 2019, 26, 9435–9442. [Google Scholar] [CrossRef]
  7. Lee, D.-W.; Kim, M.-S.; Lim, Y.-H.; Lee, N.; Hong, Y.-C. Prenatal and Postnatal Exposure to Di-(2-Ethylhexyl) Phthalate and Neurodevelopmental Outcomes: A Systematic Review and Meta-Analysis. Environ. Res. 2018, 167, 558–566. [Google Scholar] [CrossRef]
  8. Koch, H.M.; Preuss, R.; Angerer, J. Di(2-Ethylhexyl)Phthalate (DEHP): Human Metabolism and Internal Exposure—An Update and Latest Results1. Int. J. Androl. 2006, 29, 155–165. [Google Scholar] [CrossRef]
  9. Zota, A.R.; Calafat, A.M.; Woodruff, T.J. Temporal Trends in Phthalate Exposures: Findings from the National Health and Nutrition Examination Survey, 2001–2010. Environ. Health Perspect. 2014, 122, 235–241. [Google Scholar] [CrossRef]
  10. Woodruff, T.J.; Zota, A.R.; Schwartz, J.M. Environmental Chemicals in Pregnant Women in the United States: NHANES 2003–2004. Environ. Health Perspect. 2011, 119, 878–885. [Google Scholar] [CrossRef] [Green Version]
  11. Shin, H.-M.; Dhar, U.; Calafat, A.M.; Nguyen, V.; Schmidt, R.J.; Hertz-Picciotto, I. Temporal Trends of Exposure to Phthalates and Phthalate Alternatives in California Pregnant Women during 2007–2013: Comparison with Other Populations. Environ. Sci. Technol. 2020, 54, 13157–13166. [Google Scholar] [CrossRef] [PubMed]
  12. Buckley, J.P.; Kuiper, J.R.; Bennett, D.H.; Barrett, E.S.; Bastain, T.; Breton, C.V.; Chinthakindi, S.; Dunlop, A.L.; Farzan, S.F.; Herbstman, J.B.; et al. Exposure to Contemporary and Emerging Chemicals in Commerce among Pregnant Women in the United States: The Environmental Influences on Child Health Outcome (ECHO) Program. Environ. Sci. Technol. 2022, 56, 6560–6573. [Google Scholar] [CrossRef] [PubMed]
  13. Wenzel, A.G.; Brock, J.W.; Cruze, L.; Newman, R.B.; Unal, E.R.; Wolf, B.J.; Somerville, S.E.; Kucklick, J.R. Prevalence and Predictors of Phthalate Exposure in Pregnant Women in Charleston, SC. Chemosphere 2018, 193, 394–402. [Google Scholar] [CrossRef] [PubMed]
  14. Chan, M.; Mita, C.; Bellavia, A.; Parker, M.; James-Todd, T. Racial/Ethnic Disparities in Pregnancy and Prenatal Exposure to Endocrine-Disrupting Chemicals Commonly Used in Personal Care Products. Curr. Environ. Health Rep. 2021, 8, 98–112. [Google Scholar] [CrossRef] [PubMed]
  15. Latini, G.; De Felice, C.; Presta, G.; Del Vecchio, A.; Paris, I.; Ruggieri, F.; Mazzeo, P. In Utero Exposure to Di-(2-Ethylhexyl)Phthalate and Duration of Human Pregnancy. Environ. Health Perspect. 2003, 111, 1783–1785. [Google Scholar] [CrossRef]
  16. Fromme, H.; Gruber, L.; Seckin, E.; Raab, U.; Zimmermann, S.; Kiranoglu, M.; Schlummer, M.; Schwegler, U.; Smolic, S.; Völkel, W. Phthalates and Their Metabolites in Breast Milk—Results from the Bavarian Monitoring of Breast Milk (BAMBI). Environ. Int. 2011, 37, 715–722. [Google Scholar] [CrossRef]
  17. Huang, P.-C.; Tsai, C.-H.; Liang, W.-Y.; Li, S.-S.; Huang, H.-B.; Kuo, P.-L. Early Phthalates Exposure in Pregnant Women Is Associated with Alteration of Thyroid Hormones. PLoS ONE 2016, 11, e0159398. [Google Scholar] [CrossRef]
  18. Kim, J.H.; Park, H.Y.; Bae, S.; Lim, Y.-H.; Hong, Y.-C. Diethylhexyl Phthalates Is Associated with Insulin Resistance via Oxidative Stress in the Elderly: A Panel Study. PLoS ONE 2013, 8, e71392. [Google Scholar] [CrossRef] [PubMed]
  19. Ferguson, K.K.; Cantonwine, D.E.; Rivera-González, L.O.; Loch-Caruso, R.; Mukherjee, B.; AnzalotaDel Toro, L.V.; Jiménez-Vélez, B.; Calafat, A.M.; Ye, X.; Alshawabkeh, A.N.; et al. Urinary Phthalate Metabolite Associations with Biomarkers of Inflammation and Oxidative Stress Across Pregnancy in Puerto Rico. Environ. Sci. Technol. 2014, 48, 7018–7025. [Google Scholar] [CrossRef]
  20. Duan, Y.; Wang, L.; Han, L.; Wang, B.; Sun, H.; Chen, L.; Zhu, L.; Luo, Y. Exposure to Phthalates in Patients with Diabetes and Its Association with Oxidative Stress, Adiponectin, and Inflammatory Cytokines. Environ. Int. 2017, 109, 53–63. [Google Scholar] [CrossRef]
  21. Ferguson, K.K.; Chen, Y.-H.; VanderWeele, T.J.; McElrath, T.F.; Meeker, J.D.; Mukherjee, B. Mediation of the Relationship between Maternal Phthalate Exposure and Preterm Birth by Oxidative Stress with Repeated Measurements across Pregnancy. Environ. Health Perspect. 2017, 125, 488–494. [Google Scholar] [CrossRef] [PubMed]
  22. Hansen, J.M.; Jones, D.P.; Harris, C. The Redox Theory of Development. Antioxid. Redox Signal. 2019, 32, 715–740. [Google Scholar] [CrossRef] [PubMed]
  23. Manivasagam, T.; Arunadevi, S.; Essa, M.M.; SaravanaBabu, C.; Borah, A.; Thenmozhi, A.J.; Qoronfleh, M.W. Role of Oxidative Stress and Antioxidants in Autism. In Personalized Food Intervention and Therapy for Autism Spectrum Disorder Management, Advances in Neurobiology; Essa, M.M., Qoronfleh, M.W., Eds.; Springer International Publishing: Cham, Switzerland, 2020; pp. 193–206. ISBN 978-3-030-30402-7. [Google Scholar]
  24. Oulhote, Y.; Lanphear, B.; Braun, J.M.; Webster, G.M.; Arbuckle, T.E.; Etzel, T.; Forget-Dubois, N.; Seguin, J.R.; Bouchard, M.F.; MacFarlane, A.; et al. Gestational Exposures to Phthalates and Folic Acid, and Autistic Traits in Canadian Children. Environ. Health Perspect. 2020, 128, 027004. [Google Scholar] [CrossRef] [PubMed]
  25. Day, D.B.; Collett, B.R.; Barrett, E.S.; Bush, N.R.; Swan, S.H.; Nguyen, R.H.N.; Szpiro, A.A.; Sathyanarayana, S. Phthalate Mixtures in Pregnancy, Autistic Traits, and Adverse Childhood Behavioral Outcomes. Environ. Int. 2021, 147, 106330. [Google Scholar] [CrossRef] [PubMed]
  26. Kim, J.I.; Lee, J.; Lee, K.-S.; Lee, Y.A.; Shin, C.H.; Hong, Y.-C.; Kim, B.-N.; Lim, Y.-H. Association of Phthalate Exposure with Autistic Traits in Children. Environ. Int. 2021, 157, 106775. [Google Scholar] [CrossRef]
  27. Shin, H.-M.; Schmidt, R.J.; Tancredi, D.; Barkoski, J.; Ozonoff, S.; Bennett, D.H.; Hertz-Picciotto, I. Prenatal Exposure to Phthalates and Autism Spectrum Disorder in the MARBLES Study. Environ. Health 2018, 17, 85. [Google Scholar] [CrossRef]
  28. Haggerty, D.K.; Strakovsky, R.S.; Talge, N.M.; Carignan, C.C.; Glazier-Essalmi, A.N.; Ingersoll, B.R.; Karthikraj, R.; Kannan, K.; Paneth, N.S.; Ruden, D.M. Prenatal Phthalate Exposures and Autism Spectrum Disorder Symptoms in Low-Risk Children. Neurotoxicol. Teratol. 2021, 83, 106947. [Google Scholar] [CrossRef]
  29. Lord, C.; Brugha, T.S.; Charman, T.; Cusack, J.; Dumas, G.; Frazier, T.; Jones, E.J.H.; Jones, R.M.; Pickles, A.; State, M.W.; et al. Autism Spectrum Disorder. Nat. Rev. Dis. Primer 2020, 6, 5. [Google Scholar] [CrossRef]
  30. Rosen, N.E.; Lord, C.; Volkmar, F.R. The Diagnosis of Autism: From Kanner to DSM-III to DSM-5 and Beyond. J. Autism Dev. Disord. 2021, 51, 4253–4270. [Google Scholar] [CrossRef]
  31. Ozonoff, S.; Heung, K.; Byrd, R.; Hansen, R.; Hertz-Picciotto, I. The Onset of Autism: Patterns of Symptom Emergence in the First Years of Life. Autism Res. Off. J. Int. Soc. Autism Res. 2008, 1, 320–328. [Google Scholar] [CrossRef] [Green Version]
  32. Ozonoff, S.; Young, G.S.; Brian, J.; Charman, T.; Shephard, E.; Solish, A.; Zwaigenbaum, L. Diagnosis of Autism Spectrum Disorder After Age 5 in Children Evaluated Longitudinally Since Infancy. J. Am. Acad. Child Adolesc. Psychiatry 2018, 57, 849–857.e2. [Google Scholar] [CrossRef] [PubMed]
  33. Bennett, D.H.; Busgang, S.A.; Kannan, K.; Parsons, P.J.; Takazawa, M.; Palmer, C.D.; Schmidt, R.J.; Doucette, J.T.; Schweitzer, J.B.; Gennings, C.; et al. Environmental Exposures to Pesticides, Phthalates, Phenols and Trace Elements Are Associated with Neurodevelopment in the CHARGE Study. Environ. Int. 2022, 161, 107075. [Google Scholar] [CrossRef] [PubMed]
  34. Pacyga, D.C.; Sathyanarayana, S.; Strakovsky, R.S. Dietary Predictors of Phthalate and Bisphenol Exposures in Pregnant Women. Adv. Nutr. 2019, 10, 803–815. [Google Scholar] [CrossRef] [PubMed]
  35. Hertz-Picciotto, I.; Schmidt, R.J.; Krakowiak, P. Understanding Environmental Contributions to Autism: Causal Concepts and the State of Science. Autism Res. 2018, 11, 554–586. [Google Scholar] [CrossRef]
  36. Tanner, S.; Thomson, S.; Drummond, K.; O’Hely, M.; Symeonides, C.; Mansell, T.; Saffery, R.; Sly, P.D.; Collier, F.; Burgner, D.; et al. A Pathway-Based Genetic Score for Oxidative Stress: An Indicator of Host Vulnerability to Phthalate-Associated Adverse Neurodevelopment. Antioxidants 2022, 11, 659. [Google Scholar] [CrossRef]
  37. Zhang, X.; Barr, D.B.; Dunlop, A.L.; Panuwet, P.; Sarnat, J.A.; Lee, G.E.; Tan, Y.; Corwin, E.J.; Jones, D.P.; Ryan, P.B.; et al. Assessment of Metabolic Perturbations Associated with Exposure to Phthalates among Pregnant African American Women. Sci. Total Environ. 2021, 818, 151689. [Google Scholar] [CrossRef]
  38. Grindler, N.M.; Vanderlinden, L.; Karthikraj, R.; Kannan, K.; Teal, S.; Polotsky, A.J.; Powell, T.L.; Yang, I.V.; Jansson, T. Exposure to Phthalate, an Endocrine Disrupting Chemical, Alters the First Trimester Placental Methylome and Transcriptome in Women. Sci. Rep. 2018, 8, 6086. [Google Scholar] [CrossRef]
  39. Paquette, A.G.; MacDonald, J.; Lapehn, S.; Bammler, T.; Kruger, L.; Day, D.B.; Price, N.D.; Loftus, C.; Kannan, K.; Marsit, C.; et al. A Comprehensive Assessment of Associations between Prenatal Phthalate Exposure and the Placental Transcriptomic Landscape. Environ. Health Perspect. 2021, 129, 97003. [Google Scholar] [CrossRef]
  40. Orozco, J.S.; Hertz-Picciotto, I.; Abbeduto, L.; Slupsky, C.M. Metabolomics Analysis of Children with Autism, Idiopathic-Developmental Delays, and Down Syndrome. Transl. Psychiatry 2019, 9, 1–15. [Google Scholar] [CrossRef]
  41. Smith, A.M.; Natowicz, M.R.; Braas, D.; Ludwig, M.A.; Ney, D.M.; Donley, E.L.R.; Burrier, R.E.; Amaral, D.G. A Metabolomics Approach to Screening for Autism Risk in the Children’s Autism Metabolome Project. Autism Res. Off. J. Int. Soc. Autism Res. 2020, 13, 1270–1285. [Google Scholar] [CrossRef]
  42. Hollowood, K.; Melnyk, S.; Pavliv, O.; Evans, T.; Sides, A.; Schmidt, R.J.; Hertz-Picciotto, I.; Elms, W.; Guerrero, E.; Kruger, U.; et al. Maternal Metabolic Profile Predicts High or Low Risk of an Autism Pregnancy Outcome. Res. Autism Spectr. Disord. 2018, 56, 72–82. [Google Scholar] [CrossRef] [PubMed]
  43. Ritz, B.; Yan, Q.; Uppal, K.; Liew, Z.; Cui, X.; Ling, C.; Inoue, K.; von Ehrenstein, O.; Walker, D.I.; Jones, D.P. Untargeted Metabolomics Screen of Mid-Pregnancy Maternal Serum and Autism in Offspring. Autism Res. 2020, 13, 1258–1269. [Google Scholar] [CrossRef] [PubMed]
  44. Schmidt, R.J.; Liang, D.; Busgang, S.A.; Curtin, P.; Giulivi, C. Maternal Plasma Metabolic Profile Demarcates a Role for Neuroinflammation in Non-Typical Development of Children. Metabolites 2021, 11, 545. [Google Scholar] [CrossRef]
  45. Zhu, Y.; Gao, H.; Huang, K.; Zhang, Y.; Cai, X.; Yao, H.; Mao, L.; Ge, X.; Zhou, S.; Xu, Y.; et al. Prenatal Phthalate Exposure and Placental Size and Shape at Birth: A Birth Cohort Study. Environ. Res. 2018, 160, 239–246. [Google Scholar] [CrossRef]
  46. Warner, G.R.; Dettogni, R.S.; Bagchi, I.C.; Flaws, J.A.; Graceli, J.B. Placental Outcomes of Phthalate Exposure. Reprod. Toxicol. 2021, 103, 1–17. [Google Scholar] [CrossRef]
  47. Hertz-Picciotto, I.; Schmidt, R.; Walker, C.K.; Bennett, D.H.; Oliver, M.; Shedd-Wise, K.M.; LaSalle, J.M.; Giulivi, C.; Puschner, B.; Thomas, J.; et al. A Prospective Study of Environmental Exposures and Early Biomarkers in Autism Spectrum Disorder: Design, Protocols, and Preliminary Data from the MARBLES Study. Environ. Health Perspect. 2018, 126, 117004. [Google Scholar] [CrossRef]
  48. Lord, C.; Risi, S.; Lambrecht, L.; Cook, E.H.; Leventhal, B.L.; DiLavore, P.C.; Pickles, A.; Rutter, M. The Autism Diagnostic Observation Schedule-Generic: A Standard Measure of Social and Communication Deficits Associated with the Spectrum of Autism. J. Autism Dev. Disord. 2000, 30, 205–223. [Google Scholar] [CrossRef] [PubMed]
  49. Mullen, E.M.; American Guidance Service. Mullen Scales of Early Learning; AGS: Circle Pines, MN, USA, 1995. [Google Scholar]
  50. Ozonoff, S.; Young, G.S.; Belding, A.; Hill, M.; Hill, A.; Hutman, T.; Johnson, S.; Miller, M.; Rogers, S.J.; Schwichtenberg, A.J.; et al. The Broader Autism Phenotype in Infancy: When Does It Emerge? J. Am. Acad. Child Adolesc. Psychiatry 2014, 53, 398–407.e2. [Google Scholar] [CrossRef]
  51. Hasegawa, Y.; Otoki, Y.; McClorry, S.; Coates, L.C.; Lombardi, R.L.; Taha, A.Y.; Slupsky, C.M. Optimization of a Method for the Simultaneous Extraction of Polar and Non-Polar Oxylipin Metabolites, DNA, RNA, Small RNA, and Protein from a Single Small Tissue Sample. Methods Protoc. 2020, 3, 61. [Google Scholar] [CrossRef]
  52. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  53. Keil, A.P.; Buckley, J.P.; O’Brien, K.M.; Ferguson, K.K.; Zhao, S.; White, A.J. A Quantile-Based g-Computation Approach to Addressing the Effects of Exposure Mixtures. Environ. Health Perspect. 2020, 128, 47004. [Google Scholar] [CrossRef] [PubMed]
  54. Anderson, M.J. Permutational Multivariate Analysis of Variance (PERMANOVA). In Wiley StatsRef: Statistics Reference Online; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2017; pp. 1–15. ISBN 978-1-118-44511-2. [Google Scholar]
  55. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package 2020. Available online: https://CRAN.R-project.org/package=vegan (accessed on 21 August 2021).
  56. Anderson, M.J. Permutation Tests for Univariate or Multivariate Analysis of Variance and Regression. Can. J. Fish. Aquat. Sci. 2001, 58, 626–639. [Google Scholar] [CrossRef]
  57. Textor, J.; van der Zander, B.; Gilthorpe, M.S.; Liśkiewicz, M.; Ellison, G.T. Robust Causal Inference Using Directed Acyclic Graphs: The R Package ‘Dagitty’. Int. J. Epidemiol. 2016, 45, 1887–1894. [Google Scholar] [CrossRef] [PubMed]
  58. Schmidt, R.J.; Iosif, A.-M.; Angel, E.G.; Ozonoff, S. Association of Maternal Prenatal Vitamin Use With Risk for Autism Spectrum Disorder Recurrence in Young Siblings. JAMA Psychiatry 2019, 76, 391–398. [Google Scholar] [CrossRef]
  59. Krakowiak, P.; Walker, C.K.; Bremer, A.A.; Baker, A.S.; Ozonoff, S.; Hansen, R.L.; Hertz-Picciotto, I. Maternal Metabolic Conditions and Risk for Autism and Other Neurodevelopmental Disorders. Pediatrics 2012, 129, e1121–e1128. [Google Scholar] [CrossRef]
  60. Carpenter, D.O.; Arcaro, K.; Spink, D.C. Understanding the Human Health Effects of Chemical Mixtures. Environ. Health Perspect. 2002, 110, 25–42. [Google Scholar] [CrossRef]
  61. Yaghjyan, L.; Sites, S.; Ruan, Y.; Chang, S.-H. Associations of Urinary Phthalates with Body Mass Index, Waist Circumference, and Serum Lipids Among Females: National Health and Nutrition Examination Survey 1999–2004. Int. J. Obes. 2015, 39, 994–1000. [Google Scholar] [CrossRef]
  62. Adeva-Andany, M.M.; Calvo-Castro, I.; Fernández-Fernández, C.; Donapetry-García, C.; Pedre-Piñeiro, A.M. Significance of L-Carnitine for Human Health. IUBMB Life 2017, 69, 578–594. [Google Scholar] [CrossRef]
  63. Xu, Y.; Knipp, G.T.; Cook, T.J. Effects of Di-(2-Ethylhexyl)-Phthalate and Its Metabolites on the Lipid Profiling in Rat HRP-1 Trophoblast Cells. Arch. Toxicol. 2006, 80, 293–298. [Google Scholar] [CrossRef]
  64. Xu, Y.; Agrawal, S.; Cook, T.J.; Knipp, G.T. Maternal Di-(2-Ethylhexyl)-Phthalate Exposure Influences Essential Fatty Acid Homeostasis in Rat Placenta. Placenta 2008, 29, 962–969. [Google Scholar] [CrossRef] [Green Version]
  65. Zhang, S.; Wu, Z.; Heng, J.; Tian, M.; Chen, J.; Chen, F.; Guan, W. L-Carnitine Increases Cell Proliferation and Amino Acid Transporter Expression via the Activation of Insulin-like Growth Factor I Signaling Pathway in Rat Trophoblast Cells. Food Sci. Nutr. 2020, 8, 3298–3307. [Google Scholar] [CrossRef] [PubMed]
  66. Arenas, J.; Rubio, J.C.; Martín, M.A.; Campos, Y. Biological Roles of L-Carnitine in Perinatal Metabolism. Early Hum. Dev. 1998, 53, S43–S50. [Google Scholar] [CrossRef]
  67. Bai, M.; Zeng, Q.; Chen, Y.; Chen, M.; Li, P.; Ma, Z.; Sun, D.; Zhou, H.; Zheng, C.; Zeng, S.; et al. Maternal Plasma L-Carnitine Reduction During Pregnancy Is Mainly Attributed to OCTN2-Mediated Placental Uptake and Does Not Result in Maternal Hepatic Fatty Acid β-Oxidation Decline. Drug Metab. Dispos. Biol. Fate Chem. 2019, 47, 582–591. [Google Scholar] [CrossRef] [PubMed]
  68. Fagerberg, L.; Hallström, B.M.; Oksvold, P.; Kampf, C.; Djureinovic, D.; Odeberg, J.; Habuka, M.; Tahmasebpoor, S.; Danielsson, A.; Edlund, K.; et al. Analysis of the Human Tissue-Specific Expression by Genome-Wide Integration of Transcriptomics and Antibody-Based Proteomics. Mol. Cell. Proteom. MCP 2014, 13, 397–406. [Google Scholar] [CrossRef]
  69. Chang, T.-T.; Shyu, M.-K.; Huang, M.-C.; Hsu, C.-C.; Yeh, S.-Y.; Chen, M.-R.; Lin, C.-J. Hypoxia-Mediated Down-Regulation of OCTN2 and PPARα Expression in Human Placentas and in BeWo Cells. Mol. Pharm. 2011, 8, 117–125. [Google Scholar] [CrossRef] [PubMed]
  70. Cantonwine, D.E.; Meeker, J.D.; Ferguson, K.K.; Mukherjee, B.; Hauser, R.; McElrath, T.F. Urinary Concentrations of Bisphenol A and Phthalate Metabolites Measured during Pregnancy and Risk of Preeclampsia. Environ. Health Perspect. 2016, 124, 1651–1655. [Google Scholar] [CrossRef]
  71. Bedell, S.M.; Lyden, G.R.; Sathyanarayana, S.; Barrett, E.S.; Ferguson, K.K.; Santilli, A.; Bush, N.R.; Swan, S.H.; McElrath, T.F.; Nguyen, R.H.N. First- and Third-Trimester Urinary Phthalate Metabolites in the Development of Hypertensive Diseases of Pregnancy. Int. J. Environ. Res. Public. Health 2021, 18, 10627. [Google Scholar] [CrossRef]
  72. Peluso, G.; Barbarisi, A.; Savica, V.; Reda, E.; Nicolai, R.; Benatti, P.; Calvani, M. Carnitine: An Osmolyte That Plays a Metabolic Role. J. Cell. Biochem. 2000, 80, 1–10. [Google Scholar] [CrossRef]
  73. Brusati, V.; Józwik, M.; Józwik, M.; Teng, C.; Paolini, C.; Marconi, A.M.; Battaglia, F.C. Fetal and Maternal Non-Glucose Carbohydrates and Polyols Concentrations in Normal Human Pregnancies at Term. Pediatr. Res. 2005, 58, 700–704. [Google Scholar] [CrossRef]
  74. Jauniaux, E.; Hempstock, J.; Teng, C.; Battaglia, F.C.; Burton, G.J. Polyol Concentrations in the Fluid Compartments of the Human Conceptus during the First Trimester of Pregnancy: Maintenance of Redox Potential in a Low Oxygen Environment. J. Clin. Endocrinol. Metab. 2005, 90, 1171–1175. [Google Scholar] [CrossRef] [Green Version]
  75. Arroyo, J.A.; Garcia-Jones, P.; Graham, A.; Teng, C.C.; Battaglia, F.C.; Galan, H.L. Placental TonEBP/NFAT5 Osmolyte Regulation in an Ovine Model of Intrauterine Growth Restriction. Biol. Reprod. 2012, 86, 94. [Google Scholar] [CrossRef] [PubMed]
  76. Wang, B. Molecular Mechanism Underlying Sialic Acid as an Essential Nutrient for Brain Development and Cognition. Adv. Nutr. 2012, 3, 465S–472S. [Google Scholar] [CrossRef] [PubMed]
  77. Bian, D.; Wang, X.; Huang, J.; Chen, X.; Li, H. Maternal Neu5Ac Supplementation During Pregnancy Improves Offspring Learning and Memory Ability in Rats. Front. Nutr. 2021, 8. [Google Scholar] [CrossRef] [PubMed]
  78. Abeln, M.; Albers, I.; Peters-Bernard, U.; Flächsig-Schulz, K.; Kats, E.; Kispert, A.; Tomlinson, S.; Gerardy-Schahn, R.; Münster-Kühnel, A.; Weinhold, B. Sialic Acid Is a Critical Fetal Defense against Maternal Complement Attack. J. Clin. Investig. 2019, 129, 422–436. [Google Scholar] [CrossRef]
  79. Ma, Z.; Yang, H.; Peng, L.; Kuhn, C.; Chelariu-Raicu, A.; Mahner, S.; Jeschke, U.; von Schönfeldt, V. Expression of the Carbohydrate Lewis Antigen, Sialyl Lewis A, Sialyl Lewis X, Lewis X, and Lewis Y in the Placental Villi of Patients with Unexplained Miscarriages. Front. Immunol. 2021, 12, 679424. [Google Scholar] [CrossRef] [PubMed]
  80. Marini, M.; Bonaccini, L.; Thyrion, G.D.Z.; Vichi, D.; Parretti, E.; Sgambati, E. Distribution of Sugar Residues in Human Placentas from Pregnancies Complicated by Hypertensive Disorders. Acta Histochem. 2011, 113, 815–825. [Google Scholar] [CrossRef]
  81. Chiodelli, P.; Rezzola, S.; Urbinati, C.; Federici Signori, F.; Monti, E.; Ronca, R.; Presta, M.; Rusnati, M. Contribution of Vascular Endothelial Growth Factor Receptor-2 Sialylation to the Process of Angiogenesis. Oncogene 2017, 36, 6531–6541. [Google Scholar] [CrossRef]
  82. Chen, D.; Zheng, J. Regulation of Placental Angiogenesis. Microcirculation 2014, 21, 15–25. [Google Scholar] [CrossRef]
  83. Seymore, T.N.; Rivera-Núñez, Z.; Stapleton, P.A.; Adibi, J.J.; Barrett, E.S. Phthalate Exposures and Placental Health in Animal Models and Humans: A Systematic Review. Toxicol. Sci. 2022, 188, kfac060. [Google Scholar] [CrossRef]
  84. Landaas, S. The Formation of 2-Hydroxybutyric Acid in Experimental Animals. Clin. Chim. Acta Int. J. Clin. Chem. 1975, 58, 23–32. [Google Scholar] [CrossRef]
  85. Lord, R.S.; Bralley, J.A. Clinical Applications of Urinary Organic Acids. Part I: Detoxification Markers. Altern. Med. Rev. J. Clin. Ther. 2008, 13, 205–215. [Google Scholar]
  86. Amara, I.; Salah, A.; Timoumi, R.; Annabi, E.; Scuto, M.; Trovato, A.; Neffati, F.; Calabrese, V.; Abid-Essefi, S. Effect of Di(2-Ethylhexyl) Phthalate on Nrf2-Regulated Glutathione Homeostasis in Mouse Kidney. Cell Stress Chaperones 2020, 25, 919–928. [Google Scholar] [CrossRef] [PubMed]
  87. Mohammadi, H.; Ashari, S. Mechanistic Insight into Toxicity of Phthalates, the Involved Receptors, and the Role of Nrf2, NF-ΚB, and PI3K/AKT Signaling Pathways. Environ. Sci. Pollut. Res. 2021, 28, 35488–35527. [Google Scholar] [CrossRef]
  88. Wagner-Golbs, A.; Neuber, S.; Kamlage, B.; Christiansen, N.; Bethan, B.; Rennefahrt, U.; Schatz, P.; Lind, L. Effects of Long-Term Storage at −80 °C on the Human Plasma Metabolome. Metabolites 2019, 9, 99. [Google Scholar] [CrossRef] [PubMed]
  89. Serkova, N.; Bendrick-Peart, J.; Alexander, B.; Tissot van Patot, M.C. Metabolite Concentrations in Human Term Placentae and Their Changes Due to Delayed Collection After Delivery. Placenta 2003, 24, 227–235. [Google Scholar] [CrossRef]
  90. Walejko, J.M.; Chelliah, A.; Keller-Wood, M.; Gregg, A.; Edison, A.S. Global Metabolomics of the Placenta Reveals Distinct Metabolic Profiles between Maternal and Fetal Placental Tissues Following Delivery in Non-Labored Women. Metabolites 2018, 8, 10. [Google Scholar] [CrossRef]
  91. Mitro, S.D.; Wu, J.; Rahman, M.L.; Cao, Y.; Zhu, Y.; Chen, Z.; Chen, L.; Li, M.; Hinkle, S.N.; Bremer, A.A.; et al. Longitudinal Plasma Metabolomics Profile in Pregnancy—A Study in an Ethnically Diverse U.S. Pregnancy Cohort. Nutrients 2021, 13, 3080. [Google Scholar] [CrossRef]
  92. Hertz-Picciotto, I.; Delwiche, L. The Rise in Autism and the Role of Age at Diagnosis. Epidemiol. Camb. Mass 2009, 20, 84–90. [Google Scholar] [CrossRef]
  93. Bai, D.; Yip, B.H.K.; Windham, G.C.; Sourander, A.; Francis, R.; Yoffe, R.; Glasson, E.; Mahjani, B.; Suominen, A.; Leonard, H.; et al. Association of Genetic and Environmental Factors with Autism in a 5-Country Cohort. JAMA Psychiatry 2019, 76, 1035–1043. [Google Scholar] [CrossRef]
  94. Newschaffer, C.J.; Croen, L.A.; Fallin, M.D.; Hertz-Picciotto, I.; Nguyen, D.V.; Lee, N.L.; Berry, C.A.; Farzadegan, H.; Hess, H.N.; Landa, R.J.; et al. Infant Siblings and the Investigation of Autism Risk Factors. J. Neurodev. Disord. 2012, 4, 7. [Google Scholar] [CrossRef]
  95. Hollowood-Jones, K.; Adams, J.B.; Coleman, D.M.; Ramamoorthy, S.; Melnyk, S.; James, S.J.; Woodruff, B.K.; Pollard, E.L.; Snozek, C.L.; Kruger, U.; et al. Altered Metabolism of Mothers of Young Children with Autism Spectrum Disorder: A Case Control Study. BMC Pediatr. 2020, 20, 557. [Google Scholar] [CrossRef] [PubMed]
  96. Rosenfeld, C.S. The Placenta-Brain-Axis. J. Neurosci. Res. 2021, 99, 271–283. [Google Scholar] [CrossRef] [PubMed]
  97. Loomes, R.; Hull, L.; Mandy, W.P.L. What Is the Male-to-Female Ratio in Autism Spectrum Disorder? A Systematic Review and Meta-Analysis. J. Am. Acad. Child Adolesc. Psychiatry 2017, 56, 466–474. [Google Scholar] [CrossRef]
  98. Ferri, S.L.; Abel, T.; Brodkin, E.S. Sex Differences in Autism Spectrum Disorder: A Review. Curr. Psychiatry Rep. 2018, 20, 9. [Google Scholar] [CrossRef] [PubMed]
  99. Meakin, A.S.; Cuffe, J.S.M.; Darby, J.R.T.; Morrison, J.L.; Clifton, V.L. Let’s Talk about Placental Sex, Baby: Understanding Mechanisms That Drive Female- and Male-Specific Fetal Growth and Developmental Outcomes. Int. J. Mol. Sci. 2021, 22, 6386. [Google Scholar] [CrossRef]
  100. Inkster, A.M.; Fernández-Boyano, I.; Robinson, W.P. Sex Differences Are Here to Stay: Relevance to Prenatal Care. J. Clin. Med. 2021, 10, 3000. [Google Scholar] [CrossRef]
  101. DiPietro, J.A.; Voegtline, K.M. The Gestational Foundation of Sex Differences in Development and Vulnerability. Neuroscience 2017, 342, 4–20. [Google Scholar] [CrossRef]
  102. Shin, H.-M.; Bennett, D.H.; Barkoski, J.; Ye, X.; Calafat, A.M.; Tancredi, D.; Hertz-Picciotto, I. Variability of Urinary Concentrations of Phthalate Metabolites during Pregnancy in First Morning Voids and Pooled Samples. Environ. Int. 2019, 122, 222–230. [Google Scholar] [CrossRef]
  103. James-Todd, T.; Meeker, J.; Huang, T.; Hauser, R.; Seely, E.; Ferguson, K.; Rich-Edwards, J.; McElrath, T. Racial and Ethnic Variations in Phthalate Metabolite Concentrations across Pregnancy. J. Expo. Sci. Environ. Epidemiol. 2017, 27, 160–166. [Google Scholar] [CrossRef] [PubMed]
  104. Wu, H.; Kupsco, A.J.; Deierlein, A.L.; Just, A.C.; Calafat, A.M.; Oken, E.; Braun, J.M.; Mercado-Garcia, A.; Cantoral, A.; Téllez-Rojo, M.M.; et al. Trends and Patterns of Phthalates and Phthalate Alternatives Exposure in Pregnant Women from Mexico City during 2007–2010. Environ. Sci. Technol. 2020, 54, 1740–1749. [Google Scholar] [CrossRef]
Figure 1. Spearman correlations between phthalate metabolites and serum metabolites (left) and placenta metabolites (right). Positive associations (red) identify metabolites where higher urinary phthalate metabolite concentrations correlated with higher metabolite concentrations, while negative associations (blue) identify metabolites where higher urinary phthalate metabolite concentrations correlated with lower metabolite concentrations. Metabolites that remained significant after FDR correction are marked with stars. For FDR-corrected p < 0.001, ***; p < 0.01, **; and p < 0.05, *.
Figure 1. Spearman correlations between phthalate metabolites and serum metabolites (left) and placenta metabolites (right). Positive associations (red) identify metabolites where higher urinary phthalate metabolite concentrations correlated with higher metabolite concentrations, while negative associations (blue) identify metabolites where higher urinary phthalate metabolite concentrations correlated with lower metabolite concentrations. Metabolites that remained significant after FDR correction are marked with stars. For FDR-corrected p < 0.001, ***; p < 0.01, **; and p < 0.05, *.
Metabolites 12 00829 g001
Table 1. Demographic characteristics of study subjects by the child of interest’s neurodevelopmental outcome: typically developing (TD), ASD, and non-typically developing (Non-TD). For each neurodevelopmental outcome, the mean ± standard deviation for continuous variables and proportion (count) for categorical variables are reported. Differences in maternal and child factors between neurodevelopmental groups were assessed using chi-squared tests for categorical variables (or Fisher’s exact test, indicated by *) and Kruskal–Wallis tests for continuous variables.
Table 1. Demographic characteristics of study subjects by the child of interest’s neurodevelopmental outcome: typically developing (TD), ASD, and non-typically developing (Non-TD). For each neurodevelopmental outcome, the mean ± standard deviation for continuous variables and proportion (count) for categorical variables are reported. Differences in maternal and child factors between neurodevelopmental groups were assessed using chi-squared tests for categorical variables (or Fisher’s exact test, indicated by *) and Kruskal–Wallis tests for continuous variables.
SerumPlacenta
TD
(n = 66)
ASD
(n = 26)
Non-TD
(n = 14)
pTD
(n = 83)
ASD
(n = 32)
Non-TD
(n = 17)
p
Maternal age at birth, years35 ± 536 ± 534 ± 40.3835 ± 535 ± 533 ± 40.40
Birth weight, g3412 ± 4483465 ± 4283287 ± 2680.203435 ± 4863445 ± 4433344 ± 3780.81
Birth year2011 ± 22011 ± 22012 ± 20.532010 ± 22010 ± 22010 ± 20.15
Pre-pregnancy BMI25.7 ± 5.027.1 ± 7.427.6 ± 10.50.8526.3 ± 6.028.1 ± 8.229.2 ± 9.80.45
Time since last meal or snack, minutes102 ± 75101 ± 6694 ± 680.95----
Gestational age at birth, weeks38.8 ± 1.539.4 ± 0.839.1 ± 0.90.3538.9 ± 1.339.6 ± 1.038.9 ± 1.30.04
Gestational age at serum collection, days233 ± 19230 ± 18238 ± 190.50----
Prenatal vitamin use in the first month of pregnancy36 (54.5%)11 (42.3%)8 (57.1%)0.5250 (60.2%)10 (31.2%)8 (47.1%)0.02
Male child34 (51.5%)18 (69.2%)8 (57.1%)0.3043 (51.8%)23 (71.9%)10 (58.8%)0.15
Cesarean delivery27 (40.9%)9 (34.6%)4 (28.6%)0.6448 (57.8%)24 (75%)11 (64.7%)0.23
Delivery payer: Government16 (24.2%)6 (23.1%)4 (28.6%)0.8717 (20.5%)7 (21.9%)6 (35.3%)0.34
Homeowner42 (63.6%)15 (57.7%)8 (57.1%)0.8247 (56.6%)17 (53.1%)11 (64.7%)0.74
At least a bachelor’s degree37 (56.1%)13 (50%)5 (35.7%)0.3749 (59%)14 (43.8%)5 (29.4%)0.05
Maternal race/ethnicity 0.43 0.44
Non-Hispanic, White36 (54.5%)14 (53.8%)5 (35.7%) 50 (60.2%)16 (50%)8 (47.1%)
Historically marginalized groups
Black/African-American1 (1.5%)1 (3.8%)2 (14.3%) 3 (3.6%)3 (9.4%)2 (11.8%)
Asian12 (18.2%)5 (19.2%)3 (21.4%) 11 (13.3%)4 (12.5%)3 (17.6%)
Hispanic, white14 (21.2%)4 (15.4%)3 (21.4%) 17 (20.5%)7 (21.9%)3 (17.6%)
Hispanic, non-white1 (1.5%)1 (3.8%)1 (7.1%) 0 (0%)1 (3.1%)1 (5.9%)
Multi-racial2 (3.0%)1 (3.8%)0 (0%) 2 (2.4%)1 (3.1%)0 (0%)
Type 2 diabetes *2 (3%)0 (0%)0 (0%)>0.92 (2.4%)0 (0%)0 (0%)>0.9
Gestational diabetes *15 (22.7%)5 (19.2%)1 (7.1%)0.4315 (18.1%)8 (25%)1 (5.9%)0.27
Hypertension *6 (9.1%)3 (11.5%)2 (14.3%)0.725 (6%)3 (9.4%)2 (11.8%)0.57
Preeclampsia *3 (4.5%)1 (3.8%)1 (7.1%)0.814 (4.8%)2 (6.2%)2 (11.8%)0.35
Maternal metabolic condition * 0.37 0.23
BMI ≤ 25 and no metabolic conditions25 (37.9%)8 (30.8%)7 (50.0%) 34 (41.0%)8 (25.0%)7 (41.2%)
25 < BMI < 30 and no metabolic conditions14 (21.2%)6 (23.1%)5 (35.7%) 16 (19.3%)7 (21.9%)5 (29.4%)
BMI ≥ 30 and no other metabolic conditions8 (12.1%)4 (15.4%)0 (0%) 14 (16.9%)5 (15.6%)2 (11.8%)
Any hypertensive disorder (without any diabetes) at any BMI2 (3.0%)3 (11.5%)1 (7.1%) 2 (2.4%)4 (12.5%)2 (11.8%)
Diabetes at any BMI17 (25.8%)5 (19.2%)1 (7.1%) 17 (20.5%)8 (25.0%)1 (5.9%)
Table 2. Estimates of the total and direct effects of a 1 quartile increase of each phthalate metabolite in the mixture of each placental metabolite. The estimate and 95% confidence interval (CI) for each placental metabolite were bootstrapped independently using quantile g-computation models. The total effect was adjusted for birth year, homeownership, maternal education, and maternal race/ethnicity. The direct effect models were adjusted for birth year, delivery mode, gestational age at birth, homeownership, maternal education, maternal metabolic condition, and maternal race/ethnicity. All p-values were adjusted for false discovery rate (FDR) and both the original and FDR p-values are reported.
Table 2. Estimates of the total and direct effects of a 1 quartile increase of each phthalate metabolite in the mixture of each placental metabolite. The estimate and 95% confidence interval (CI) for each placental metabolite were bootstrapped independently using quantile g-computation models. The total effect was adjusted for birth year, homeownership, maternal education, and maternal race/ethnicity. The direct effect models were adjusted for birth year, delivery mode, gestational age at birth, homeownership, maternal education, maternal metabolic condition, and maternal race/ethnicity. All p-values were adjusted for false discovery rate (FDR) and both the original and FDR p-values are reported.
Total EffectDirect Effect
MetaboliteEstimate (95% CI)pFDR pEstimate (95% CI)pFDR p
1,3-Dihydroxyacetone0.0034 (−0.1237, 0.1306)0.95790.9579−0.0019 (−0.1452, 0.1415)0.97950.9983
2-Hydroxybutyrate−0.0804 (−0.1378, −0.0230)0.00710.0896−0.0925 (−0.1562, −0.0287)0.00530.0907
3-Hydroxybutyrate−0.0600 (−0.1955, 0.0754)0.38680.8106−0.0655 (−0.2129, 0.0820)0.38600.9046
4-Aminobutyrate0.0081 (−0.0747, 0.0910)0.84760.95790.0223 (−0.0632, 0.1079)0.61030.9679
4-Hydroxybutyrate−0.0103 (−0.1715, 0.1508)0.90020.9579−0.0166 (−0.1975, 0.1643)0.85740.9781
Acetate−0.0231 (−0.1363, 0.0901)0.68990.9292−0.0240 (−0.1405, 0.0925)0.68710.9764
Alanine−0.0146 (−0.0659, 0.0367)0.57800.9292−0.0122 (−0.0669, 0.0425)0.66320.9679
Arginine−0.0250 (−0.0814, 0.0313)0.38530.8106−0.0285 (−0.0954, 0.0384)0.40580.9046
Asparagine−0.0175 (−0.0706, 0.0356)0.51920.8869−0.0156 (−0.0729, 0.0416)0.59350.9679
Aspartate−0.0323 (−0.0800, 0.0153)0.18630.6288−0.0372 (−0.0891, 0.0147)0.16340.5882
Betaine−0.0473 (−0.0983, 0.0036)0.07140.3613−0.0445 (−0.1051, 0.0162)0.15370.5882
Carnitine−0.0738 (−0.1207, −0.0270)0.00250.0783−0.0785 (−0.1306, −0.0265)0.00380.0907
Choline−0.0231 (−0.0756, 0.0294)0.39030.8106−0.0237 (−0.0791, 0.0317)0.40400.9046
Creatine−0.0132 (−0.0538, 0.0274)0.52560.8869−0.0047 (−0.0546, 0.0452)0.85370.9781
Cystine 0.0588 (−0.0470, 0.1646)0.27860.8106 0.0647 (−0.0471, 0.1765)0.25930.7779
Ethanolamine−0.0173 (−0.0679, 0.0333)0.50390.8869−0.0214 (−0.0773, 0.0344)0.45420.9084
Formate−0.0382 (−0.1018, 0.0255)0.24250.7703−0.0336 (−0.0955, 0.0283)0.28940.8225
Fumarate−0.0628 (−0.1519, 0.0263)0.16970.6185−0.0622 (−0.1668, 0.0424)0.24630.7779
Glucitol−0.0912 (−0.1498, −0.0326)0.00290.0783−0.0956 (−0.1596, −0.0315)0.00420.0907
Glutamate−0.0160 (−0.0608, 0.0288)0.48550.8869−0.0191 (−0.0684, 0.0302)0.44860.9084
Glutamine−0.0185 (−0.0561, 0.0191)0.33640.8106−0.0257 (−0.0637, 0.0124)0.18930.6389
Glutathione−0.1151 (−0.2481, 0.0180)0.09280.4176−0.1437 (−0.2729, −0.0144)0.03150.2430
Glycerol−0.0149 (−0.0791, 0.0493)0.65030.9292−0.0160 (−0.0837, 0.0516)0.64320.9679
Glycine−0.0105 (−0.0579, 0.0369)0.66550.9292−0.0045 (−0.0585, 0.0494)0.86940.9781
Hypoxanthine−0.0167 (−0.0657, 0.0324)0.50700.8869−0.0173 (−0.0712, 0.0366)0.53100.9679
Inosine−0.0733 (−0.1778, 0.0312)0.17180.6185−0.0995 (−0.2103, 0.0112)0.08100.3365
Isoleucine−0.0099 (−0.0652, 0.0454)0.72610.9294−0.0032 (−0.0637, 0.0572)0.91670.9900
Kynurenine−0.0718 (−0.1461, 0.0025)0.06080.3613−0.0776 (−0.1594, 0.0042)0.06560.2952
Lactate 0.0089 (−0.0317, 0.0495)0.66870.9292 0.0082 (−0.0360, 0.0525)0.71640.9781
Leucine−0.0118 (−0.0693, 0.0458)0.68940.9292−0.0037 (−0.0666, 0.0591)0.90760.9900
Lysine−0.0307 (−0.0935, 0.0320)0.33950.8106−0.0192 (−0.0866, 0.0482)0.57790.9679
Methionine 0.0147 (−0.0464, 0.0758)0.63790.9292 0.0202 (−0.0459, 0.0864)0.54970.9679
myo-Inositol 0.0223 (−0.0222, 0.0668)0.32860.8106 0.0176 (−0.0357, 0.0709)0.51870.9679
N-Acetylneuraminate−0.0844 (−0.1460, −0.0228)0.00830.0896−0.0901 (−0.1559, −0.0244)0.00840.0907
NAD+−0.1044 (−0.2397, 0.0309)0.13330.5537−0.1315 (−0.2700, 0.0070)0.06540.2952
Niacinamide 0.0255 (−0.2252, 0.2762)0.84220.9579−0.0034 (−0.2684, 0.2616)0.97980.9983
O-Acetylcarnitine−0.0665 (−0.1150, −0.0181)0.00820.0896−0.0733 (−0.1256, −0.0209)0.00710.0907
O-Phosphocholine−0.1516 (−0.2954, −0.0077)0.04120.3178−0.1647 (−0.3153, −0.0140)0.03440.2430
O-Phosphoethanolamine−0.1191 (−0.2315, −0.0067)0.04010.3178−0.1212 (−0.2405, −0.0019)0.04900.2940
Ornithine−0.0055 (−0.1014, 0.0903)0.90980.9579 0.0129 (−0.0879, 0.1137)0.80240.9781
Pantothenate−0.0384 (−0.1141, 0.0373)0.32260.8106−0.0406 (−0.1219, 0.0407)0.32970.8902
Phenylalanine−0.0141 (−0.0703, 0.0421)0.62320.9292−0.0085 (−0.0702, 0.0532)0.78780.9781
Proline 0.0019 (−0.0495, 0.0533)0.94290.9579 0.0078 (−0.0470, 0.0627)0.78030.9781
Pyroglutamate−0.0272 (−0.0888, 0.0344)0.38900.8106−0.0281 (−0.0961, 0.0398)0.41880.9046
Serine−0.0044 (−0.0643, 0.0555)0.88620.9579 0.0001 (−0.0665, 0.0667)0.99830.9983
sn-Glycero−3-phosphocholine−0.0390 (−0.1353, 0.0574)0.42950.8590−0.0388 (−0.1299, 0.0523)0.40540.9046
Succinate−0.0096 (−0.1030, 0.0838)0.84030.9579−0.0243 (−0.1149, 0.0662)0.59950.9679
Taurine−0.0366 (−0.0761, 0.0028)0.07140.3613−0.0406 (−0.0819, 0.0008)0.05700.2952
Threonine 0.0048 (−0.0500, 0.0596)0.86350.9579 0.0103 (−0.0495, 0.0701)0.73630.9781
Tryptophan 0.0877 (−0.0075, 0.1829)0.07360.3613 0.1052 (0.0081, 0.2022)0.03600.2430
Tyrosine 0.0027 (−0.0511, 0.0565)0.92250.9579 0.0080 (−0.0508, 0.0668)0.79010.9781
Uracil 0.0090 (−0.0656, 0.0835)0.81410.9579 0.0070 (−0.0744, 0.0884)0.86710.9781
Uridine−0.0154 (−0.0953, 0.0644)0.70550.9292−0.0199 (−0.1088, 0.0690)0.66190.9679
Valine−0.0087 (−0.0597, 0.0424)0.74010.9294−0.0013 (−0.0575, 0.0549)0.96500.9983
Table 3. Associations between placenta metabolites and birth year. In samples with a molar sum of urinary DEHP metabolites below the median (n = 66), we conducted multiple linear regressions to measure the association between each log10-transformed placenta metabolite concentration and birth year after adjustment for gestational age at delivery, maternal metabolic condition, homeownership status, maternal race/ethnicity, and maternal education. We conducted permutation testing under the reduced model using 4999 permutations for each metabolite and corrected for false discovery rate (FDR). Only metabolites with FDR p-values < 0.1 are reported.
Table 3. Associations between placenta metabolites and birth year. In samples with a molar sum of urinary DEHP metabolites below the median (n = 66), we conducted multiple linear regressions to measure the association between each log10-transformed placenta metabolite concentration and birth year after adjustment for gestational age at delivery, maternal metabolic condition, homeownership status, maternal race/ethnicity, and maternal education. We conducted permutation testing under the reduced model using 4999 permutations for each metabolite and corrected for false discovery rate (FDR). Only metabolites with FDR p-values < 0.1 are reported.
MetaboliteEstimate (95% CI)r2pFDR p
Arginine−0.0502 (−0.0965, −0.0038)0.06810.03440.0978
Creatine0.0405 (0.0040, 0.0771)0.07120.03380.0978
Cystine−0.1050 (−0.1792, −0.0308)0.11100.00640.0432
Glutathione0.1238 (0.0390, 0.2087)0.11720.00520.0401
Inosine0.1328 (0.0678, 0.1978)0.20630.00000.0000
Isoleucine−0.0407 (−0.0765, −0.0049)0.07450.02780.0945
Leucine−0.0478 (−0.0842, −0.0114)0.09720.01220.0549
Lysine−0.0483 (−0.0912, −0.0054)0.07340.02800.0945
Methionine−0.0539 (−0.0941, −0.0137)0.10070.00740.0444
NAD+0.2018 (0.1186, 0.2849)0.26860.00000.0000
O-Phosphocholine0.1244 (0.0402, 0.2086)0.11980.00400.0360
O-Phosphoethanolamine0.0935 (0.0324, 0.1545)0.12750.00380.0360
Ornithine−0.0536 (−0.1026, −0.0045)0.06930.03200.0978
Phenylalanine−0.0461 (−0.0815, −0.0108)0.09600.01100.0549
Pyroglutamate−0.0974 (−0.1457, −0.0490)0.20180.00000.0000
Serine−0.0527 (−0.0933, −0.0121)0.09510.01160.0549
Tryptophan−0.0934 (−0.1659, −0.0208)0.09370.01380.0573
Tyrosine−0.0420 (−0.0773, −0.0067)0.08120.02180.0841
Uracil−0.0802 (−0.1318, −0.0286)0.13090.00280.0360
Table 4. The permutational multivariate analysis of variance (PERMANOVA) to assess the association between the placenta metabolome and neurodevelopmental outcome stratified by the child of interest’s sex. Placental metabolites were log10-transformed. The marginal effects of terms were tested under 9999 permutations.
Table 4. The permutational multivariate analysis of variance (PERMANOVA) to assess the association between the placenta metabolome and neurodevelopmental outcome stratified by the child of interest’s sex. Placental metabolites were log10-transformed. The marginal effects of terms were tested under 9999 permutations.
All ChildrenFemalesMales
TermdfFR2pdfFR2pdfFR2p
Birth year111.220.07770.000114.190.07130.001317.120.08520.0001
Delivery mode11.330.00920.206410.680.01160.710511.230.01470.2595
Gestational age at delivery, weeks11.150.00790.289510.950.01620.433311.110.01320.3224
Homeownership10.560.00390.815211.130.01920.308711.590.01900.1361
Maternal age, years10.740.00510.622710.910.01560.472710.470.00570.8711
Maternal race/ethnicity10.860.00590.509110.570.00970.824110.90.01080.4689
Maternal metabolic condition40.60.01670.959640.850.05780.678740.530.02510.9843
Birth weight, grams10.380.00260.957910.80.01360.558110.630.00760.7155
Prenatal vitamin use in the first month of pregnancy11.810.01260.079310.640.01100.747312.340.02800.0337
Neurodevelopmental outcome21.450.02010.127210.03420.416421.840.04400.0451
Residual117 0.8099 41 0.6980 61 0.7296
Total131 1.0000 55 1.0000 75 1.0000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Parenti, M.; Schmidt, R.J.; Ozonoff, S.; Shin, H.-M.; Tancredi, D.J.; Krakowiak, P.; Hertz-Picciotto, I.; Walker, C.K.; Slupsky, C.M. Maternal Serum and Placental Metabolomes in Association with Prenatal Phthalate Exposure and Neurodevelopmental Outcomes in the MARBLES Cohort. Metabolites 2022, 12, 829. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12090829

AMA Style

Parenti M, Schmidt RJ, Ozonoff S, Shin H-M, Tancredi DJ, Krakowiak P, Hertz-Picciotto I, Walker CK, Slupsky CM. Maternal Serum and Placental Metabolomes in Association with Prenatal Phthalate Exposure and Neurodevelopmental Outcomes in the MARBLES Cohort. Metabolites. 2022; 12(9):829. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12090829

Chicago/Turabian Style

Parenti, Mariana, Rebecca J. Schmidt, Sally Ozonoff, Hyeong-Moo Shin, Daniel J. Tancredi, Paula Krakowiak, Irva Hertz-Picciotto, Cheryl K. Walker, and Carolyn M. Slupsky. 2022. "Maternal Serum and Placental Metabolomes in Association with Prenatal Phthalate Exposure and Neurodevelopmental Outcomes in the MARBLES Cohort" Metabolites 12, no. 9: 829. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo12090829

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop