Study population
The present study used data from the Mothers and Children's Environmental Health (MOCEH) study. The MOCEH study was a prospective birth cohort study to investigate children’s and adolescents' growth, development, and cognitive function before and after birth with environmental exposures based on hospitals and communities in South Korea. Details of the MOCEH study have been previously described [36, 37]. Briefly, from 2006 to 2010, a total of 1,751 pregnant women in their first trimester (less than 20 weeks of gestation) was enrolled from three university hospitals located in three cities in Korea: Seoul (metropolitan area), Cheonan (urban area), and Ulsan (metropolitan and industrial area). Using a schedule of planned visits, the mother-child pairs were monitored based on the protocols at the enrolled hospitals. Each participant provided written informed consent at the first visit. The Institutional Review Boards approved this study at Ewha Woman’s University (IRB No. 12-07B-15), Dankook University Hospital (IRB No. 2011-09-0340), and Ulsan University Hospital (IRB No. 06–29).
The flow chart of inclusion is presented in Fig. 1. Among the enrolled participants, 567 remained after excluding those who did not participate in the IQ test (n = 1,184). Of these, 181 children had measured metabolite concentrations of EDCs at ages of 24 and 36 months. After further excluding participants who had missing information on covariates (n = 28), 152 children were included in the analyses.
Measurement Of BPA And Phthalate Exposure
We measured children’s urinary concentrations of BPA, MEOHP, MEHHP, and MnBP at ages of 24 and 36 months. Samples of spot urine in 100-mL sterile cups were kept in freezers at − 20°C until they were delivered to the laboratory. The measurement method used to achieve accuracy and reliability of metabolite measurements is described in the Centers for Disease Control and Prevention Manual of Laboratory Procedures [38, 39]. The limits of detection (LODs) of BPA, MEOHP, MEHHP, and MnBP were 0.275, 0.65, 0.53, and 0.87 µg/L. None of the 4 chemicals samples contained levels below the LOD. Quality control analysis was conducted using functional materials from a German External Quality Assessment Scheme for phthalates and BPA in the laboratory. Urine creatinine concentrations were used to convert BPA and phthalate concentrations [36, 40].
Cognitive Assessment
Trained physicians measured the children’s IQ with the Korean Wechsler Preschool and Primary Scale of Intelligence, fourth edition (K-WPPSI-IV) at age 5 years. The K-WPPSI-IV is the most common test for comprehensive assessment of cognitive function of children in Korea. The K-WPPSI-IV consists of 15 subdomains, five basic indicators, and four additional indicators [41]. The MOCEH database management system validated answers and numerical data entry with a web-based database program using a mechanized screening system.
Potential Confounding Factors
Potential confounding factors were selected a priori from the literature, and a directed acyclic graph was used to select covariates included in the models. Information on child sex was collected from delivery room records. Maternal age, education level, household income level, and child’s exposure to secondhand smoke were collected through a questionnaire at the planned visits. Maternal IQ was obtained with the Korean Wechsler Adult Intelligence Scale, fourth edition, at the pre- or postnatal period for those who agreed to take the test, and there were no significant differences in characteristics of those who did or did not take the test (Table S1). The following variables were considered as categorical variables: maternal education level (below bachelor’s degree, bachelor’s degree, higher than bachelor’s degree), household monthly income level (< KRW 2,000,000, KRW 2,000,000–4,000,000, > KRW4,000,000), child sex (boy or girl), and exposure to secondhand smoke (yes or no). The variables considered as continuous were maternal age and IQ and child’s age. Smoking during pregnancy was not included because none of the mothers were smokers.
Statistical analysis
For categorical variables, data were presented as frequency and percentage, and a chi-square test was used for comparison. For continuous variables, data were presented as mean and standard deviation, and the t-test and Wilcoxon rank sum test were used for comparison. Pearson’s correlation test was conducted for each of the four phthalate metabolites at each visit. All concentrations of BPA and phthalate metabolites were log-transformed in statistical analyses. To simplify comparison of the magnitude of the effects of different metabolites, we took the concentrations of each BPA and phthalate metabolite on a log scale and normalized them to mean and standard deviation.
First, multiple linear regression was conducted to examine the association of children’s exposure to each chemical at 24 and 36 months of age with children’s IQ. In model 1, we adjusted for potentially confounding factors of maternal age and education, household income, child’s secondhand smoke exposure, child’s sex, and child’s age. In model 2, we additionally adjusted for maternal IQ.
Second, we used a weighted quantile sum (WQS) regression [42] and quantile-based g-computation [43] to assess the associations between exposure to mixtures of BPA and phthalates and children’s IQ. Both statistical models combined each estimated associations into weighted exponents and evaluated specific weights for individual chemicals in the mixture. The estimated weights described the relative contribution, with higher weights indicating a larger contribution to the overall association. A bootstrap procedure with 100 repetitions was conducted to calculate the statistical significance of the resulting WQS index. For the WQS regression analyses, we used 40% of the data as a training set and the remaining 60% as the validation set. The training data set was used to estimate chemical weights, and the validation set tested the WQS index. The analysis was based on the following formula:
\(Y={\beta }_{0}+\beta \left(\sum _{i}^{c}{w}_{i}{q}_{i}\right)+\) covariates
where Y is the outcome (child’s IQ), β0 is the intercept, and β is the WQS index estimate. qi is the quantile of the ith element, wi is the weight (estimated) associated with the ith element. The \(\left(\sum _{i}^{c}{w}_{i}{q}_{i}\right)\) term represents the sum of the weights of the components included. The β of the WQS index represents the decrease in child's IQ for a one-quartile increase in the WQS index. The β of the quantile-based g-computation index represents the decrease in child’s IQ with an increase in the first quartile for all exposures included in the analysis. WQS estimates unidirectionally (positive or negative) for all individual associations within a mixture, but in quantile-based g-computations individual contributions to mixture association can be bidirectional. Therefore, we performed WQS analyses assuming negative associations, assuming the increased exposure to EDCs mixture would lower the IQ. The WQS and quantile-based g-computation analyses were adjusted for the same potential confounders of the multiple regression. The “qWQS” and” qgcomp” packages (version 0.2.0) in R version 4.1.2 were used for analysis [44].
We conducted sex-stratified analyses to examine sex differences using multiple regression, WQS, and quantile-based g-computation analyses. Due to the small sample size, we could not adjust mother’s IQ in the sex-stratified analysis.