Racial and ethnic disparities in preterm birth: a mediation analysis incorporating mixtures of polybrominated diphenyl ethers

Background: Racial and ethnic disparities persist in preterm birth (PTB) and gestational age (GA) at delivery in the United States. It remains unclear whether exposure to environmental chemicals contributes to these disparities.

Objectives: We applied recent methodologies incorporating environmental mixtures as mediators in causal mediation analysis to examine whether racial and ethnic disparities in GA at delivery and PTB may be partially explained by exposures to polybrominated diphenyl ethers (PBDEs), a class of chemicals used as flame retardants in the United States.

Methods: Data from a multiracial/ethnic US cohort of 2008 individuals with low-risk singleton pregnancies were utilized, with plasma PBDE concentrations measured during early pregnancy. We performed mediation analyses incorporating three forms of mediators: (1) reducing all PBDEs to a weighted index, (2) selecting a PBDE congener, or (3) including all congeners simultaneously as multiple mediators, to evaluate whether PBDEs may contribute to the racial and ethnic disparities in PTB and GA at delivery, adjusted for potential confounders.

Results: Among the 2008 participants, 552 self-identified as non-Hispanic White, 504 self-identified as non-Hispanic Black, 568 self-identified as Hispanic, and 384 self-identified as Asian/Pacific Islander. The non-Hispanic Black individuals had the highest mean ∑PBDEs, the shortest mean GA at delivery, and the highest rate of PTB. Overall, the difference in GA at delivery comparing non-Hispanic Black to non-Hispanic White women was −0.30 (95% CI: −0.54, −0.05) weeks. This disparity reduced to −0.23 (95% CI: −0.49, 0.02) and −0.18 (95% CI: −0.46, 0.10) weeks if fixing everyone's weighted index of PBDEs to the median and the 25th percentile levels, respectively. The proportion of disparity mediated by the weighted index of PBDEs was 11.8%. No statistically significant mediation was found for PTB, other forms of mediator(s), or other racial and ethnic groups.

Conclusion: PBDE mixtures may partially mediate the Black vs. White disparity in GA at delivery. While further validations are needed, lowering the PBDEs at the population level might help reduce this disparity.

1 Introduction

Preterm birth (PTB) affects 9%–10% of pregnancies in the United States, and is associated with increased risk of maternal and neonatal morbidity and mortality (1, 2). There are pronounced racial and ethnic disparities in PTB in the United States, with rates disproportionately higher in non-Hispanic Black women than non-Hispanic White women (14% vs. 9%) (3, 4). These disparities may further contribute to higher infant mortality (4) among non-Hispanic Black relative to non-Hispanic White infants. For other groups, studies showed no significant difference in PTB rate comparing Asian or Hispanic women to White women, although the risk appeared higher in certain Asian subgroups (5). Therefore, identifying the potentially modifiable risk factors of PTB, especially those that are unevenly distributed across racial and ethnic groups, is important, to help understand and reduce the disparities in PTB.

The existing literature suggests that disparities in PTB are largely attributable to environmental factors rather than genetic variation (3, 6, 7). These include social stressors, physical stressors (such as environmental chemicals and pollutants), neighborhood variation, healthcare access/quality, and individual cultural practices (8). One study suggested that certain sociodemographic and perinatal health factors contributed to the Black vs. White disparity in PTB, although they reported that more than 60% of the disparities in PTB remained unexplained (9). Other studies also showed that the Black vs. White disparities in PTB persisted after accounting for socioeconomic status, access to care, or medical interventions (10–12). For environmental pollutants, multiple studies revealed associations of air pollution, lead, phthalates, and other chemicals with increased risk of PTB, and found higher exposure levels among non-Hispanic Black women compared with the non-Hispanic White women (13–20). However, it remains unclear whether and what proportion of racial and ethnic disparities in PTB is attributable to different exposures to these environmental factors. A causal mediation analysis (21) is needed to further explore the role of multiple environmental factors on the racial and ethnic disparities in PTB.

There have been recent calls for and developments in the methodology of evaluating environmental factors as potential mediators of health disparities (22–24). Furthermore, given that people are often simultaneously exposed to multiple environmental factors (25, 26), a growing body of conceptual models and statistical methods integrating the joint effects of multiple pollutants into a mediation analysis framework has been proposed, especially in the field of environmental health disparities (27–30). These methods can help quantify the proportion of disparity due to environmental factors, as well as the proportion of disparity that would remain if interventions were made to reduce the levels of these environmental factors. Despite the discussions on this framework and the related methods, a real-world, population-based application of these methods in evaluating the contribution of environmental chemicals/pollutants as a mixture to a health disparity question remains lacking.

One class of environmental chemicals, known as polybrominated diphenyl ethers (PBDEs), has been used as a flame retardant since the 1970s and remains to be detected in the US population even a decade after the voluntary phase out that began in 2004 (31–33). PBDEs have the potential to shed or volatilize into the environment (34). Human beings are exposed to PBDEs via inhalation of contaminated air, ingestion of contaminated food, and contact with indoor dust. We hypothesize that PBDEs might be potential mediators for the racial and ethnic disparities in PTB given the following evidence: (1) multiple studies showed higher exposure levels to PBDEs among non-Hispanic Black women compared with non-Hispanic White women (35–37); (2) studies have found associations between certain PBDE congeners and elevated risk of PTB (38–42). In this study, we aimed to use real-world data from a large, multicenter, multiracial/ethnic cohort of singleton pregnancies in the United States to evaluate whether and the extent to which exposure to PBDEs may contribute to the racial and ethnic disparity in PTB and gestational age at delivery, through applying causal mediation analyses incorporating these chemicals (individually and as mixtures) as potential mediators. Race and ethnicity are socially constructed, and racial/ethnic health disparities are driven by the root cause of structural/institutional racism. With that in mind, we present a causal diagram (43) of our research questions in Figure 1.

Figure 1

Figure 1. Causal diagram for studying racial and ethnic disparities in preterm birth and gestatonal age at delivery, mediated by polybrominated diphenyl ethers. The boxes in yellow correspond to the variables that we were able to measure with available data in this study. The boxes in gray correspond to the variables on the causal pathway that we were not able to measure in this study.

2 Methods

2.1 Study population

The study used data from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) Fetal Growth Studies—Singleton Cohort, a multicenter, multiracial/ethnic prospective study of 2,802 pregnant women recruited during 2009–2013 from 12 US clinical sites (44). Women aged 18–40 years with a singleton pregnancy were enrolled during 8–13 weeks of gestation and followed through delivery. Further details of the cohort can be found elsewhere (44, 45). For this study, we restricted to a subcohort of 2008 eligible women with a low-risk pregnancy (i.e., those with certain pre-existing medical conditions such as systemic diseases or past pregnancy complications were excluded from enrolling in the study) (44, 45) and without obesity [i.e., individuals whose body mass index (BMI) < 30 kg/m²], who had available data on gestational age at delivery and measurements of PBDEs from blood specimens. The rationale of these criteria and the numbers excluded are summarized in the Supplementary Material (Supplemental eMethod). Approval for human subjects’ research was obtained from the institutional review boards at all participating sites, and all participants provided informed consent.

2.2 Race and ethnicity

Self-identified race and ethnicity were collected at baseline in four categories: non-Hispanic White, non-Hispanic Black, Hispanic, and Asian/Pacific Islander. Further specifications such as self-reported Hispanic origin or Asian background were evaluated in secondary analyses. Non-Hispanic White was defined as the reference group. Too few Hispanic White (n = 4) and Hispanic Black (n = 4) participants were included to consider these groups separately. We use the self-identified race and ethnicity as the “predictor” parameter in the mediation analysis, while recognizing that race is a social construct (46) that may through racism impact differences in exposures to PBDEs and their sources, as well as differences in factors contributing to PTB or shorter gestational age at delivery, including pathophysiology and access to/quality of prenatal care (7) (Figure 1). As race and ethnicity are non-manipulable, the effect estimates from the mediation analysis should be interpreted as associations reflecting disparity-related (instead of causal/biological) information (23), but we maintained the usage of “effects” when describing these measures to be consistent with common causal mediation terminologies.

2.3 Outcomes

The primary outcomes of interest were: (1) gestational age at delivery (weeks), calculated as the difference between date of delivery (abstracted from medical records) and self-reported date of first day of last menstrual period (LMP) as validated by ultrasound (47); and (2) a binary outcome of PTB, defined as delivery prior to 37 weeks of gestation. As secondary outcomes, PTB was further categorized as very early or moderate (<34 weeks) and late (34 to <37 weeks) PTB.

2.4 Mediators

A set of potential mediators was determined based on prior knowledge (42), which included plasma concentrations of polybrominated biphenyl (PBB) 153 and 9 PBDEs (PBDE 28, 47, 85, 99, 100, 153, 154, 183, and 209) collected upon enrollment (median: 11 weeks of gestation). Details of the processing, measurement, and limits of quantification (LOQs) of these chemicals have been reported previously (48). All chemical concentrations were reported as ng/mL plasma. For this analysis, we restricted to six PBDEs with quantification rates >30% in this population, including PBDE 28, 47, 99, 100, 153, and 154. Machine-observed values were used for all chemicals in the analysis without substitution, including concentrations below the LOQ (49).

2.5 Covariates

The following covariates (collected from the baseline questionnaire unless otherwise specified) were incorporated into our mediation analyses, based on a priori knowledge of being potential confounders for the mediator-outcome associations: maternal age (years); prepregnancy BMI (kg/m²), calculated from self-recalled prepregnancy weight divided by measured height squared (50); parity (0, 1, 2+); education level (college degree, some college/undergraduate, graduate/postgraduate); marital status (married or living with partner, not married); family income during last year (<$30,000, $30,000–$49,999, $50,000–$99,999, ≥$100,000, not reported); plasma cotinine level (ng/mL), measured in specimens collected at enrollment (35); plasma total lipids (non-fasting) (ng/mL) at enrollment, quantified using commercially available enzymatic methods (51), and calculated as total cholesterol × 2.27 + triglycerides + 62.3 (52); total and sedentary activities [metabolic equivalent of task (MET) hours/week]; and acculturation status (US-born, recent immigrant, long-term immigrant) based on previous definitions (53). It is possible that race and ethnicity are associated with various downstream risk factors, which might violate the assumption of no mediator-outcome confounders affected by the exposure (54). To address this, we conducted sensitivity analyses using more generalized approaches (23), with details described in the statistical analysis.

2.6 Statistical analysis

2.6.1 Descriptive analysis

The characteristics of the study population were summarized with means ± standard deviations or numbers (percentages). Geometric means (GMs) and 95% confidence intervals (CIs) of lipid-adjusted PBDE congener concentrations and their molar sum (∑PBDEs) were calculated, stratified by race and ethnicity and by PTB status.

2.6.2 Mediation analysis

For mediation analysis, we natural log-transformed the machine-observed values of the chemical concentrations to account for skewedness of their distributions, and then performed standardization (subtracted the mean and divided by the standard deviation) to generate comparable scales. The total racial and ethnic (denoted by X) disparity in PTB or gestational age at delivery (denoted by Y) accounting for a set of covariates (denoted by C) was calculated using: $E[Y|X,C]={\alpha }_0+{\alpha }_{1}X+{\alpha }_{2}C$ (when Y represents continuous gestational age at delivery, in weeks), or $logit\{Pr[Y=1|X,C]\}={\alpha }_0+{\alpha }_{1}X+{\alpha }_{2}C$ (where $Y=1$ represents PTB and $Y=0$ represents non-PTB). The following forms of mediator(s) were then evaluated within a counterfactual framework using causal mediation models (for simplicity, we use a continuous variable Y as an illustration).

2.6.2.1 Reducing the PBDEs mixtures to a single mediator—weighted quantile sum

As the first approach, we reduced the dimensions of the PBDEs mixtures to a single summary index score via the weighted quantile sum (WQS) approach, which is a method that constructs a weighted index estimating the mixture effect associated with all predictor variables on an outcome (55). The weights for each PBDE were empirically determined using a 40%/60% split of training/validation sets from the data and 500 bootstrap samples for parameter estimation. Next, the WQS index was treated as a single summary measure of the PBDE congeners, and was included as a single mediator in the following models:

$E[Y|X,\mathrm{WQS},C]={\alpha }_0^{\mathrm{\prime }}+{\alpha }_1^{\mathrm{\prime }}X+{\alpha }_2^{\mathrm{\prime }}\mathrm{WQS}+{\alpha }_3^{\mathrm{\prime }}C$

$E[\mathrm{WQS}|X,C]={\beta }_0^{\mathrm{\prime }}+{\beta }_1^{\mathrm{\prime }}X+{\beta }_2^{\mathrm{\prime }}C$

The direct and indirect effects through this single mediator were estimated using standard regression-based methods (56).

2.6.2.2 Reducing the number of mediators—select specific mediator(s)

As the second approach, we reduced the number of mediators by selecting a single specific mediator based on the results of a previous study utilizing data from the same cohort of individuals, where multiple statistical approaches [including generalized linear models, principal component analysis, and Bayesian kernel machine regression (BKMR) (57)] have consistently demonstrated PBDE 153 being the main congener associated with shorter gestation and higher risk of PTB, after adjusting for race/ethnicity and other covariates (42). In this study, we further utilized a hierarchical BKMR variable selection approach based on correlation structures of PBDEs in this cohort (which address the potential bias introduced by highly correlated chemicals) to re-evaluate that PBDE 153 is the most important contributor that is associated with gestational age at delivery.

In this approach, we used a single mediator (PBDE 153 as an example) in the following models:

$E[Y|X,\mathrm{PBDE}153,C]={\alpha }_0^{\ast }+{\alpha }_1^{\ast }X+{\alpha }_2^{\ast }\mathrm{PBDE}153+{\alpha }_3^{\ast }C$

$E[\mathrm{PBDE}153|X,C]={\beta }_0^{\ast }+{\beta }_1^{\ast }X+{\beta }_2^{\ast }C$

The direct and indirect effects through this single mediator were estimated using regression-based methods (56).

2.6.2.3 Modeling all six PBDE congeners as multiple mediators—multiple regression

As the third approach, we included PBDE 28, 47, 99, 100, 153, and 154 simultaneously in the same model:

$\begin{array}{rl}& E\left[Y|X,\mathrm{PBDE}s28,47,99,100,153,154,C\right]\\ & ={\alpha }_0^{\mathrm{\prime }\mathrm{\prime }}+{\alpha }_1^{\mathrm{\prime }\mathrm{\prime }}X+{\alpha }_2^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}28+{\alpha }_3^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}47\\ & +{\alpha }_4^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}99+{\alpha }_5^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}100\\ & +{\alpha }_6^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}153+{\alpha }_7^{\mathrm{\prime }\mathrm{\prime }}\mathrm{PBDE}154+{\alpha }_8^{\mathrm{\prime }\mathrm{\prime }}C\end{array}$

along with six separate regression models estimating each mediator as a function of the exposure:

$E[\mathrm{PBDE}28|X,C]={\beta }_{0\mathrm{\_}28}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}28}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}28}^{\mathrm{\prime }\mathrm{\prime }}C$

$E[\mathrm{PBDE}47|X,C]={\beta }_{0\mathrm{\_}47}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}47}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}47}^{\mathrm{\prime }\mathrm{\prime }}C$

$E[\mathrm{PBDE}99|X,C]={\beta }_{0\mathrm{\_}99}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}99}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}99}^{\mathrm{\prime }\mathrm{\prime }}C$

$E[\mathrm{PBDE}100|X,C]={\beta }_{0\mathrm{\_}100}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}100}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}100}^{\mathrm{\prime }\mathrm{\prime }}C$

$E[\mathrm{PBDE}153|X,C]={\beta }_{0\mathrm{\_}153}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}153}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}153}^{\mathrm{\prime }\mathrm{\prime }}C$

$E[\mathrm{PBDE}154|X,C]={\beta }_{0\mathrm{\_}154}^{\mathrm{\prime }\mathrm{\prime }}+{\beta }_{1\mathrm{\_}154}^{\mathrm{\prime }\mathrm{\prime }}X+{\beta }_{2\mathrm{\_}154}^{\mathrm{\prime }\mathrm{\prime }}C$

The direct and indirect effects (specifically, the joint mediated effect through the set of mediators) were estimated using regression-based methods for multiple mediators (58).

In all the approaches, we estimated the following measures of the disparities in gestational age at delivery and PTB mediated by PBDEs, comparing each of the race and ethnicity groups to the non-Hispanic White group: the total effect (TE), the controlled direct effects (CDEs) while fixing the mediator(s) at various levels, the natural direct and indirect effects (NDE; NIE), and the overall percent mediated (PM) calculated as (NIE/TE) × 100%. All models used regression-based methods, and 95% CIs were obtained via the delta method (from closed-form parameter function estimation in single-mediator models) or bootstrapping (from direct counterfactual imputation estimation in multiple-mediator models). We further extended the models to allow for potential exposure–mediator or mediator–mediator interaction (29, 56, 59).

2.6.3 Secondary and sensitivity analysis

As secondary or sensitivity analyses, we evaluated the outcomes and mediator (WQS index) stratified by finer specifications of race and ethnicity including Hispanic origin or Asian background. We further conducted mediation analysis comparing selective subgroups to non-Hispanic White women. We also performed mediation analysis for PTB subcategories (very early/moderate PTB and late PTB). Furthermore, we evaluated mediation through the WQS index for the absolute risk difference (RD) of PTB using the g-formula approach (60).

Given that some of the proposed mediator-outcome confounders might be downstream factors of racism, hence potentially having an association with race and ethnicity, we conducted sensitivity analyses using the more generalized g-formula approach (23, 60–62), which allowed for a vector of the mediator-outcome confounders potentially affected by the exposure to be accounted for in the analysis.

We also performed the following analyses to evaluate the robustness of our main findings. First, we modeled the WQS index as a binary mediator (≥median vs. <median). Second, we evaluated potential non-linearity via categorizing the PBDEs into <LOQ and quartiles above LOQ, and the WQS index into quintiles, and we used these quantile measures as mediators. Given WQS regression's assumption of unidirectionality, we in addition explored the application of quantile g-computation (63), a flexible extension of WQS estimating the joint effects of a mixture while allowing for chemicals to act on both directions, although with the limitation of being subject to multicollinearity in the presence of highly correlated chemicals within a mixture (64). From the quantile g-computation results, we identified the PBDEs that contributed to the associations with shorter gestational age at delivery, and further created a weighted index of these chemicals as a mediator. Lastly, we conducted sensitivity analysis considering potential measurement errors of the mediator (65).

2.6.4 Statistical software

All causal mediation analyses were conducted using the CMAverse (v.0.1.0) package in R (https://bs1125.github.io/CMAverse/) (62). The WQS analyses were conducted using the gWQS (v.3.0.0) package in R (https://cran.r-project.org/web/packages/gWQS) (66).

3 Results

Among the 2008 women included in the study, 552 (27.5%) self-identified as non-Hispanic White, 504 (25.1%) self-identified as non-Hispanic Black, 568 (28.3%) self-identified as Hispanic, and 384 (19.1%) self-identified as Asian/Pacific Islander (Table 1). There were several differences in characteristics across these groups (Table 1). On average, compared with non-Hispanic White women, non-Hispanic Black women were younger, had higher BMI, lower education level, and less family income, and were more likely to be unmarried. Non-Hispanic Black women also had the highest plasma cotinine level and total and sedentary activity levels compared with other groups. Hispanic women had the highest mean BMI and plasma total lipid level, the lowest percentage of being nulliparous, and the highest percentages of attaining less than a college degree or being long-term immigrants. Asian/Pacific Islander women had the highest mean age, the lowest mean BMI, plasma cotinine level, and total activity level, as well as the highest percentage of being recent immigrants. Non-Hispanic Black women had shorter mean gestational ages at delivery (39.0 vs. 39.3 weeks) and higher risks of PTB (9.1% vs. 5.1%) compared with non-Hispanic White women. The outcomes among Hispanic or Asian/Pacific Islander women were similar to those of the non-Hispanic White women.

Table 1

Table 1. Characteristics of the study population by race and ethnicity, NICHD Fetal Growth Study–Singleton Cohort (n = 2,008).

In Table 2, non-Hispanic Black women have higher GMs of all six PBDE congeners and ∑PBDEs than non-Hispanic White women. Across all groups, Hispanic and Asian/Pacific Islander women had the highest GMs of PBDE 28 and PBDE 154, respectively. When comparing those with vs. without PTB in the study population, four PBDEs (i.e., PBDE 28, 47, 100, and 153) and ∑PBDEs had higher GMs.

Table 2

Table 2. Geometric means (95% confidence intervals) of lipid-adjusted PBDE concentrations, stratified by race and ethnicity and by preterm birth status.

We used a WQS index to estimate the mixture effect of six PBDEs on gestational age at delivery (weights for each PBDE shown in Figure 2). The association of a 1-unit increase in the WQS index with gestational age at delivery was β (95% CI) = −0.20 (−0.35, −0.05) weeks, adjusted for race and ethnicity and other covariates. Table 3 provides results from the mediation analysis, where the WQS index was considered a potential mediator for the racial and ethnic disparity of gestational age at delivery or PTB. Comparing non-Hispanic Black women with non-Hispanic White women, the covariate-adjusted difference in gestational age at delivery was β_TE (95% CI) = −0.30 (−0.54, −0.05) weeks. The CDEs (95% CIs) when fixing everyone's WQS index levels at the 25th, 50th, and 75th percentiles were −0.18 (−0.46, 0.10), −0.23 (−0.49, 0.02), and −0.32 (−0.57, −0.07) weeks, respectively. Overall, a suggestive NIE of β_NIE (95% CI) = −0.04 (−0.07, 0.00) weeks were mediated through the WQS index (proportion mediated = 11.8%). The odds ratio of PTB comparing non-Hispanic Black with non-Hispanic White women was OR_TE (95% CI) = 1.82 (1.00, 3.31), yet no statistically significant NIE was found. In addition, no statistically significant disparity was found when comparing Hispanic or Asian/Pacific Islander women with non-Hispanic White women.

Figure 2

Figure 2. Weights for each PBDE congener from the weighted quantile sum index. WQS, weighted quantile sum.

Table 3

Table 3. Estimates of direct and indirect effects mediated through a weighted quantile sum exposure index of PBDEs for the associations of race and ethnicity with gestational age at delivery and preterm birth.

The correlation coefficients between PBDEs are shown in Supplementary Figure S1. PBDE 28, 47, 99, and 100 were moderately to highly correlated, and PBDE 153 and 154 were weakly correlated. Using BKMR with hierarchical variable selection (based on the correlation structure, PBDE 28, 47, 99, and 100 were assigned as Group 1, and PBDE 153 and 154 were assigned as Group 2), we found that Group 2 was of relatively greater importance, and PBDE 153 was the most important chemical within Group 2 [reflected by the posterior inclusion probabilities (PIPs) shown in Supplementary Table S1] that was associated with shorter gestational age at delivery (Supplementary Figure S2). There were no qualitative interactions between the PBDEs (Supplementary Figure S3). Thus, for the single-mediator model, we included PBDE 153 as the mediator. Table 4 provides results from the mediation analysis, where only PBDE 153 was considered as a potential mediator. The CDEs (95% CIs) when fixing everyone's PBDE 153 levels at the 25th, 50th, and 75th percentiles were −0.18 (−0.45, 0.09), −0.18 (−0.45, 0.09), and −0.29 (−0.54, −0.05) weeks, respectively, yet with a non-significant NIE mediated via PBDE 153 (proportion mediated = 7.9%). No statistically significant NIE was found for the non-Hispanic Black vs. non-Hispanic White disparity in PTB.

Table 4

Table 4. Estimates of direct and indirect effects mediated through PBDE 153 concentrations for the associations of race and ethnicity with gestational age at delivery and preterm birth.

Table 5 provides results when all six PBDEs were included simultaneously as multiple mediators. The CDEs (95% CIs) when fixing all PBDEs at the 25th, 50th, and 75th percentiles were −0.25 (−0.53, 0.02), −0.14 (−0.47, 0.17), and −0.27 (−0.61, −0.03) weeks, respectively, yet with a non-significant NIE jointly mediated via PBDE 28, 46, 99, 100, 153, and 154 (proportion mediated = 16.3%). No statistically significant NIE was found for the disparity in PTB. No exposure–mediator(s) or mediator–mediator interaction was found for any of the aforementioned analyses (p-values for interactions >0.05).

Table 5

Table 5. Estimates of direct and indirect effects mediated through concentrations of all six PBDE congeners (PBDE 28, 47, 99, 100, 153, and 154) for the associations of race and ethnicity with gestational age at delivery and preterm birth.

As secondary analysis, the outcomes stratified by further specified Hispanic origin or Asian background are provided in Supplementary Table S2. Several subgroups had shorter mean gestational ages at delivery than non-Hispanic White participants. Among them, those reporting Filipino background also had higher mean WQS index levels than non-Hispanic White women (Supplementary Figure S4). Mediation analysis results comparing Filipino to non-Hispanic White women are provided in Supplementary Table S3. The proportion mediated by WQS index for the shorter gestational age at delivery was 16% (95% CI: −11%, 61%) when comparing the n = 45 women with Filipino background with non-Hispanic White women.

The results from the sensitivity analyses to evaluate the robustness of our main results are shown in the Supplementary Materials. In summary, evaluating PTB in subcategories showed no significant findings (Supplementary Table S4), and evaluating the absolute risk of PTB on the risk difference scale showed similar findings (Supplementary Table S5). Using g-estimation yielded similar patterns of mediation (Supplementary Table S6). Modeling the WQS index as a binary mediator dichotomized at the median had very little impact on the indirect effect estimates (Supplementary Table S7). The associations between PBDE 153 or the WQS index and gestational age at delivery were linear (p-trend < 0.05, Supplementary Table S8), and modeling WQS or PBDE 153 as quantiles showed similar mediation effects (Supplementary Tables S9, S10). Using a weighted exposure index of 4 PBDEs (PBDE 28,99, 100, and 153) based on quantile g-computation analysis (selecting the PBDEs with weights toward an association with shortened gestational age at delivery, Supplementary Figure S5) yielded similar mediation patterns as the WQS index (Supplementary Table S11). Finally, Supplementary Figure S6 showed similar estimates across various magnitudes of potential mediator measurement error.

4 Discussion

In this multiracial/ethnic cohort of pregnant women in the United States, we demonstrated shorter gestational age at delivery, higher risk of PTB, and higher exposure levels to PBDEs among non-Hispanic Black than in non-Hispanic White women, and we evaluated potential mediation by PBDEs for the racial and ethnic disparities in gestational age at delivery and PTB utilizing several recently developed causal mediation approaches. In particular, we observed that a weighted index summarizing PBDEs as a mixture had a suggestive mediating role in the Black vs. White disparity in gestational age at delivery that accounted for 11.8% of the total disparity. No significant mediation was found for the disparity of PTB, or from evaluating other forms of PBDE mediators. We also revealed disparities in gestational age at delivery comparing the Filipino subgroup with non-Hispanic White women, although no significant mediation via PBDEs was found. While further validation using larger datasets are needed, our results point to the possibility of PBDE mixtures acting as mediators for the existing racial and ethnic disparities in gestational age at delivery.

Our observation of a higher risk of PTB and shorter mean gestational age at delivery among non-Hispanic Black women compared with non-Hispanic White women is consistent with previous reports (4, 9, 67). However, the PTB risks were lower than the general population since this study consisted of relatively healthy, non-obese individuals. We also observed higher average concentrations of PBDEs comparing non-Hispanic Black with non-Hispanic White women, which aligned with previous studies that reported similar disparity patterns among certain PBDEs (e.g., 28, 47, 99, and 100) (36, 68). This disparity may be explained by differences in social and contextual factors contributing to sources of PBDE exposures, such as differences in residential neighborhoods (69), housing (e.g., indoor dust) (70), and furniture PBDE exposures (68). Furthermore, we observed higher mean levels of four PBDE congeners and ∑PBDEs in those who delivered preterm compared with non-preterm, suggesting that PBDEs, either individually or as mixtures, might be potential mediator(s) accounting for part of the disparities in PTB or gestational age at delivery. In our causal mediation analyses, we found that a weighted index of all PBDEs (i.e., the WQS index) accounted for 11.8% of the total Black vs. White disparity in gestational age at delivery. Particularly, we found that the CDEs were closer to the null when fixing everyone's WQS index at lower levels, suggesting the potential benefit in reducing the existing disparity in gestational age at delivery by intervening on PBDE levels in the entire population. Conversely, the proportion mediated by the WQS index for the Black vs. White disparity in PTB was only 4% (and non-significant), which could be explained either by lower statistical power owing to a limited number of events, or that the magnitude of mediation for gestational age at delivery might be relatively small to make a noticeable impact on the risk of PTB in this healthier population. Given this is the first study that evaluated the potential mediation role of PBDEs for this disparity, future studies of larger sample sizes or conducted among a higher-risk population might be needed to validate our findings. Past studies have revealed other mediators (such as socioeconomic and health factors, and access to healthcare) for the racial and ethnic disparity in PTB, but a large proportion of the disparity remained (9–11, 71–73). If PBDEs truly mediate part of the racial and ethnic disparity in length of gestation or PTB, then this class of chemicals might be an additional modifiable factor to help further alleviate this disparity.

Similarly to previous literature (5), we did not observe significant differences in gestational age at delivery or PTB comparing Hispanic or Asian/Pacific Islander women with non-Hispanic women. However, we did find a 38 per 1,000 births higher risk of PTB and 0.5-week shorter mean gestational age at delivery comparing a subgroup of Filipino women with non-Hispanic White women, which was consistent with previous studies showing that Filipino women had higher relative risk of PTB (compared with non-Hispanic White) than other Asian subgroups (74). Despite these Filipino women also having higher exposure levels to PBDEs, the results from mediation analysis were non-significant. This might be due to the small number of participants with various Asian backgrounds, although we could not rule out the possibility that there might be unmeasured confounding such as cultural, psychosocial, or early life factors that are driving this disparity, especially when more than half of the Filipino women in this study were immigrants. Future studies with more specific focus on these racial/ethnic minority subgroups that collect acculturation-related variables are needed to further explore this mediation.

In this study, we compared three different approaches of incorporating PBDEs as potential mediators of racial and ethnic disparities in gestational age at delivery and PTB: reducing to a WQS index, selecting a single PBDE 153 congener according to prior knowledge and its relative importance from the hierarchical BKMR selection, and including six PBDE congeners as multiple mediators. Overall, the estimated proportion mediated via the single PBDE 153 congener was smaller than that via the WQS index. This is possibly owing to the limitation of selecting mediator(s) a priori based on the mediator-outcome association alone, which might leave out important mediator(s) weakly associated with the outcome that may also contribute to the indirect effect. The estimated proportion jointly mediated by multiple PBDEs was higher than the proportion mediated by the WQS index, but with much wider CIs due to potential overfitting or multicollinearity. Our example showed that the WQS approach carries the advantage of reducing the PBDEs to a single score to avoid overfitting or multicollinearity, while preserving the information from each congener, serving as a suitable approach to explore the overall contribution of a chemical mixture to a health disparity question (28).

We acknowledge several limitations of this study. First, this study consisted of women with low-risk of adverse health outcomes at baseline and without obesity, so our findings might not be fully generalizable to the overall US population. Second, unmeasured confounding was inevitable, such as other geographic, psychosocial, or lifestyle factors. Third, statistical power was limited when evaluating potential mediation within certain subgroups. Lastly, we were not able to directly measure historical or contemporary environmental racism or adverse environments in these data that are contributing to (or on the causal pathway for) the observed disparities where race and ethnicity act as a proxy for these complex processes (75). Further studies are needed to inform interventions on the policies and systems level.

This study has many unique strengths. First, this study utilized a prospective cohort design in a large, racially/ethnically diverse population with clinically validated outcomes and a comprehensive set of covariates. Second, we applied different statistical approaches to evaluate mediation(s) through individual as well as mixtures of PBDEs. Third, efforts were made to evaluate mediation for disparities in subcategories of PTB, or among other under-studied racial and ethnic subgroups (e.g., based on Asian backgrounds). Lastly, we conducted various sensitivity analyses to validate the robustness of our findings.

5 Conclusions

In conclusion, in this multiracial/ethnic cohort of pregnant women in the United States, we found that non-Hispanic Black women had shorter gestational ages at delivery, higher risk of PTB, and higher exposures to PBDEs compared with non-Hispanic White women. Our mediation analysis provided suggestive evidence that the Black vs. White disparity in gestational age at delivery might be partially mediated by disparities in exposures to PBDEs. Lowering the PBDE exposures at the population level may help reduce this disparity.

Data availability statement

The datasets presented in this article are not readily available because of data usage and confidentiality agreements. Requests to access the data should be directed to the corresponding author, the senior author, and the Eunice Kennedy Shriver National Institute of Child Health and Human Development Intramural Research Program. Requests to access the datasets should be directed to Tamarra James-Todd, tjtodd@hsph.harvard.edu.

Ethics statement

The studies involving humans were approved by Eunice Kennedy Shriver National Institute of Child Health and Human Development. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

ZW: Conceptualization, Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review and editing. CZ: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Resources, Supervision, Writing – review and editing. PW: Conceptualization, Investigation, Methodology, Supervision, Validation, Writing – review and editing. AB: Conceptualization, Investigation, Methodology, Supervision, Writing – review and editing. BW: Investigation, Methodology, Writing – review and editing. KK: Investigation, Methodology, Validation, Writing – review and editing. MB: Funding acquisition, Investigation, Methodology, Writing – review and editing. KH: Funding acquisition, Investigation, Methodology, Writing – review and editing. TJ-T: Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article.

This work is supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development Intramural Research Program (contract numbers HHSN275200800013C, HHSN275200800002I, HHSN27500006, HHSN275200800003IC, HHSN275200800014C, HHSN275200800012C, HHSN275200800028C, HHSN275201000009C, HHSN275291199991I, and HHSN275201000001Z), the National Institute of Environmental Health Sciences (contract number P30ES000002), and the Office of the Director, National Institutes of Health under Award Number UG3OD023316. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Health. Blair J. Wylie is supported by NIEHS R01 ES028688.

Acknowledgments

The authors acknowledge the research teams at all participating clinical centers in the NICHD Fetal Growth Studies, including Christina Care Health Systems, University of California, Irvine, Long Beach Memorial Medical Center, Northwestern University, Medical University of South Carolina, Columbia University, New York Presbyterian Queens, Queens, St. Peters' University Hospital, University of Alabama at Birmingham, Women and Infants Hospital of Rhode Island, Fountain Valley Regional Hospital and Medical Center, and Tufts University. The authors also acknowledge C-TASC and The EMMES Corporations in providing data and imaging support for this multisite study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frph.2023.1285444/full#supplementary-material

References

1. Martin JA, Osterman MJK. Describing the increase in preterm births in the United States, 2014–2016. NCHS Data Brief. (2018) (312):1–8. PMID: 30044213

Google Scholar

2. Institute of Medicine (US) Committee on Understanding Premature Birth and Assuring Healthy Outcomes. Preterm birth: causes, consequences, and prevention. In: Behrman RE, Butler AS, editors. The National Academies Collection: Reports Funded by National Institutes of Health. Washington (DC): National Academies Press (US) (2007). Available at: https://www.ncbi.nlm.nih.gov/books/NBK11362/ (accessed December 04, 2023).

3. Burris HH, Collins JW, Wright RO. Racial/ethnic disparities in preterm birth: clues from environmental exposures. Curr Opin Pediatr. (2011) 23(2):227–32. doi: 10.1097/MOP.0b013e328344568f

PubMed Abstract | Crossref Full Text | Google Scholar

4. Burris HH, Lorch SA, Kirpalani H, Pursley DM, Elovitz MA, Clougherty JE. Racial disparities in preterm birth in USA: a biosensor of physical and social environmental exposures. Arch Dis Child. (2019) 104(10):931–5. doi: 10.1136/archdischild-2018-316486

PubMed Abstract | Crossref Full Text | Google Scholar

5. Schaaf J, Liem S, Mol B, Abu-Hanna A, Ravelli A. Ethnic and racial disparities in the risk of preterm birth: a systematic review and meta-analysis. Amer J Perinatol. (2012) 30(6):433–50. doi: 10.1055/s-0032-1326988

Crossref Full Text | Google Scholar

6. Anum EA, Springel EH, Shriver MD, Strauss JF. Genetic contributions to disparities in preterm birth. Pediatr Res. (2009) 65(1):1–9. doi: 10.1203/PDR.0b013e31818912e7

PubMed Abstract | Crossref Full Text | Google Scholar

7. Braveman P, Dominguez TP, Burke W, Dolan SM, Stevenson DK, Jackson FM, et al. Explaining the black–white disparity in preterm birth: a consensus statement from a multi-disciplinary scientific work group convened by the march of dimes. Front Reprod Health. (2021) 3:684207. doi: 10.3389/frph.2021.684207

PubMed Abstract | Crossref Full Text | Google Scholar

8. Burris HH, Wright CJ, Kirpalani H, Collins JW Jr., Lorch SA, Elovitz MA, et al. The promise and pitfalls of precision medicine to resolve black–white racial disparities in preterm birth. Pediatr Res. (2020) 87(2):221–6. doi: 10.1038/s41390-019-0528-z

PubMed Abstract | Crossref Full Text | Google Scholar

9. Thoma ME, Drew LB, Hirai AH, Kim TY, Fenelon A, Shenassa ED. Black–white disparities in preterm birth: geographic, social, and health determinants. Am J Prev Med. (2019) 57(5):675–86. doi: 10.1016/j.amepre.2019.07.007

PubMed Abstract | Crossref Full Text | Google Scholar

10. Johnson JD, Green CA, Vladutiu CJ, Manuck TA. Racial disparities in prematurity persist among women of high socioeconomic Status. Am J Obstet Gynecol MFM. (2020) 2(3):100104. doi: 10.1016/j.ajogmf.2020.100104

PubMed Abstract | Crossref Full Text | Google Scholar

11. Thurston H, Fields BE, White J. Does increasing access to prenatal care reduce racial disparities in birth outcomes? J Pediatr Nurs. (2021) 59:96–102. doi: 10.1016/j.pedn.2021.01.012

PubMed Abstract | Crossref Full Text | Google Scholar

12. Timofeev J, Singh J, Istwan N, Rhea D, Driggers RW. Spontaneous preterm birth in African-American and Caucasian women receiving 17α-hydroxyprogesterone caproate. Am J Perinatol. (2014) 31(1):55–60. doi: 10.1055/s-0033-1334452

PubMed Abstract | Crossref Full Text | Google Scholar

13. Liu WY, Yu ZB, Qiu HY, Wang JB, Chen XY, Chen K. Association between ambient air pollutants and preterm birth in Ningbo, China: a time-series study. BMC Pediatr. (2018) 18(1):305. doi: 10.1186/s12887-018-1282-9

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ash M, Boyce JK. Racial disparities in pollution exposure and employment at US industrial facilities. Proc Natl Acad Sci U.S.A. (2018) 115(42):10636–41. doi: 10.1073/pnas.1721640115

PubMed Abstract | Crossref Full Text | Google Scholar

15. Muller C, Sampson RJ, Winter AS. Environmental inequality: the social causes and consequences of lead exposure. Annu Rev Sociol. (2018) 44(1):263–82. doi: 10.1146/annurev-soc-073117-041222

Crossref Full Text | Google Scholar

16. Andrews KW, Savitz DA, Hertz-Picciotto I. Prenatal lead exposure in relation to gestational age and birth weight: a review of epidemiologic studies. Am J Ind Med. (1994) 26(1):13–32. doi: 10.1002/ajim.4700260103

PubMed Abstract | Crossref Full Text | Google Scholar

17. Ferguson KK, McElrath TF, Meeker JD. Environmental phthalate exposure and preterm birth. JAMA Pediatr. (2014) 168(1):61–7. doi: 10.1001/jamapediatrics.2013.3699

PubMed Abstract | Crossref Full Text | Google Scholar

18. Silva MJ, Barr DB, Reidy JA, Malek NA, Hodge CC, Caudill SP, et al. Urinary levels of seven phthalate metabolites in the U.S. Population from the National ealth and Nutrition Examination Survey (NHANES) 1999–2000. Environ Health Perspect. (2004) 112(3):331–8. doi: 10.1289/ehp.6723

PubMed Abstract | Crossref Full Text | Google Scholar

19. Daouda M, Henneman L, Kioumourtzoglou MA, Gemmill A, Zigler C, Casey JA. Association between county-level coal-fired power plant pollution and racial disparities in preterm births from 2000 to 2018. Environ Res Lett. (2021) 16(3):034055. doi: 10.1088/1748-9326/abe4f7

PubMed Abstract | Crossref Full Text | Google Scholar

20. Harris SM, Colacino J, Buxton M, Croxton L, Nguyen V, Loch-Caruso R, et al. A data mining approach reveals chemicals detected at higher levels in non-Hispanic black women target preterm birth genes and pathways. Reprod Sci. (2022) 29(7):2001–12. doi: 10.1007/s43032-022-00870-w

PubMed Abstract | Crossref Full Text | Google Scholar

21. VanderWeele TJ. Mediation analysis: a practitioner’s guide. Annu Rev Public Health. (2016) 37:17–32. doi: 10.1146/annurev-publhealth-032315-021402

PubMed Abstract | Crossref Full Text | Google Scholar

22. James-Todd TM, Chiu YH, Zota AR. Racial/ethnic disparities in environmental endocrine disrupting chemicals and women’s reproductive health outcomes: epidemiological examples across the life course. Curr Epidemiol Rep. (2016) 3(2):161–80. doi: 10.1007/s40471-016-0073-9

PubMed Abstract | Crossref Full Text | Google Scholar

23. Naimi AI, Schnitzer ME, Moodie EEM, Bodnar LM. Mediation analysis for health disparities research. Am J Epidemiol. (2016) 184(4):315–24. doi: 10.1093/aje/kwv329

PubMed Abstract | Crossref Full Text | Google Scholar

24. Schulz A, Northridge ME. Social determinants of health: implications for environmental health promotion. Health Educ Behav. (2004) 31(4):455–71. doi: 10.1177/1090198104265598

PubMed Abstract | Crossref Full Text | Google Scholar

25. Aylward LL, Kirman CR, Schoeny R, Portier CJ, Hays SM. Evaluation of biomonitoring data from the CDC National Exposure Report in a risk assessment context: perspectives across chemicals. Environ Health Perspect. (2013) 121(3):287–94. doi: 10.1289/ehp.1205740

PubMed Abstract | Crossref Full Text | Google Scholar

26. Kortenkamp A, Faust M, Scholze M, Backhaus T. Low-level exposure to multiple chemicals: reason for human health concerns? Environ Health Perspect. (2007) 115(Suppl 1):106–14. doi: 10.1289/ehp.9358

PubMed Abstract | Crossref Full Text | Google Scholar

27. Bellavia A, Zota AR, Valeri L, James-Todd T. Multiple mediators approach to study environmental chemicals as determinants of health disparities. Environ Epidemiol. (2018) 2(2):e015. doi: 10.1097/EE9.0000000000000015

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bellavia A, James-Todd T, Williams PL. Approaches for incorporating environmental mixtures as mediators in mediation analysis. Environ Int. (2019) 123:368–74. doi: 10.1016/j.envint.2018.12.024

PubMed Abstract | Crossref Full Text | Google Scholar

29. Bellavia A, Valeri L. Decomposition of the total effect in the presence of multiple mediators and interactions. Am J Epidemiol. (2018) 187(6):1311–8. doi: 10.1093/aje/kwx355

PubMed Abstract | Crossref Full Text | Google Scholar

30. Taylor KW, Joubert BR, Braun JM, Dilworth C, Gennings C, Hauser R, et al. Statistical approaches for assessing health effects of environmental chemical mixtures in epidemiology: lessons from an innovative workshop. Environ Health Perspect. (2016) 124(12):A227–9. doi: 10.1289/EHP547

PubMed Abstract | Crossref Full Text | Google Scholar

31. Eskenazi B, Chevrier J, Rauch SA, Kogut K, Harley KG, Johnson C, et al. In utero and childhood polybrominated diphenyl ether (PBDE) exposures and neurodevelopment in the CHAMACOS study. Environ Health Perspect. (2013) 121(2):257–62. doi: 10.1289/ehp.1205597

PubMed Abstract | Crossref Full Text | Google Scholar

32. Chevrier J, Harley KG, Bradman A, Gharbi M, Sjödin A, Eskenazi B. Polybrominated diphenyl ether (PBDE) flame retardants and thyroid hormone during pregnancy. Environ Health Perspect. (2010) 118(10):1444–9. doi: 10.1289/ehp.1001905

PubMed Abstract | Crossref Full Text | Google Scholar

33. Johnson-Restrepo B, Kannan K. An assessment of sources and pathways of human exposure to polybrominated diphenyl ethers in the United States. Chemosphere. (2009) 76(4):542–8. doi: 10.1016/j.chemosphere.2009.02.068

PubMed Abstract | Crossref Full Text | Google Scholar

34. Allen JG, McClean MD, Stapleton HM, Webster TF. Linking PBDEs in house dust to consumer products using x-ray fluorescence. Environ Sci Technol. (2008) 42(11):4222–8. doi: 10.1021/es702964a

PubMed Abstract | Crossref Full Text | Google Scholar

35. Buck Louis GM, Zhai S, Smarr MM, Grewal J, Zhang C, Grantz KL, et al. Endocrine disruptors and neonatal anthropometry, NICHD Fetal Growth Studies—singletons. Environ Int. (2018) 119:515–26. doi: 10.1016/j.envint.2018.07.024

PubMed Abstract | Crossref Full Text | Google Scholar

36. Nguyen VK, Kahana A, Heidt J, Polemi K, Kvasnicka J, Jolliet O, et al. A comprehensive analysis of racial disparities in chemical biomarker concentrations in United States women, 1999–2014. Environ Int. (2020) 137:105496. doi: 10.1016/j.envint.2020.105496

PubMed Abstract | Crossref Full Text | Google Scholar

37. Varshavsky JR, Sen S, Robinson JF, Smith SC, Frankenfield J, Wang Y, et al. Racial/ethnic and geographic differences in polybrominated diphenyl ether (PBDE) levels across maternal, placental, and fetal tissues during mid-gestation. Sci Rep. (2020) 10:12247. doi: 10.1038/s41598-020-69067-y

PubMed Abstract | Crossref Full Text | Google Scholar

38. Harley KG, Chevrier J, Aguilar Schall R, Sjödin A, Bradman A, Eskenazi B. Association of prenatal exposure to polybrominated diphenyl ethers and infant birth weight. Am J Epidemiol. (2011) 174(8):885–92. doi: 10.1093/aje/kwr212

PubMed Abstract | Crossref Full Text | Google Scholar

39. Eick SM, Hom Thepaksorn EK, Izano MA, Cushing LJ, Wang Y, Smith SC, et al. Associations between prenatal maternal exposure to per- and polyfluoroalkyl substances (PFAS) and polybrominated diphenyl ethers (PBDEs) and birth outcomes among pregnant women in San Francisco. Environ Health. (2020) 19(1):100. doi: 10.1186/s12940-020-00654-2

PubMed Abstract | Crossref Full Text | Google Scholar

40. Behnia F, Peltier MR, Saade GR, Menon R. Environmental pollutant polybrominated diphenyl ether, a flame retardant, induces primary amnion cell senescence. Am J Reprod Immunol. (2015) 74(5):398–406. doi: 10.1111/aji.12414

PubMed Abstract | Crossref Full Text | Google Scholar

41. Gao Y, Chen L, Wang C, Zhou Y, Wang Y, Zhang Y, et al. Exposure to polybrominated diphenyl ethers and female reproductive function: a study in the production area of Shandong, China. Sci Total Environ. (2016) 572:9–15. doi: 10.1016/j.scitotenv.2016.07.181

PubMed Abstract | Crossref Full Text | Google Scholar

42. Wang Z, Zhang C, Williams PL, Bellavia A, Wylie BJ, Hacker MR, et al. Polybrominated diphenyl ethers in early pregnancy and preterm birth: findings from the NICHD Fetal Growth Studies. Int J Hyg Environ Health. (2022) 243:113978. doi: 10.1016/j.ijheh.2022.113978

PubMed Abstract | Crossref Full Text | Google Scholar

43. Howe CJ, Bailey ZD, Raifman JR, Jackson JW. Recommendations for using causal diagrams to study racial health disparities. Am J Epidemiol. (2022) 191(12):1981–9. doi: 10.1093/aje/kwac140

PubMed Abstract | Crossref Full Text | Google Scholar

44. Grewal J, Grantz KL, Zhang C, Sciscione A, Wing DA, Grobman WA, et al. Cohort profile: NICHD Fetal Growth Studies—singletons and twins. Int J Epidemiol. (2018) 47(1):25–25l. doi: 10.1093/ije/dyx161

PubMed Abstract | Crossref Full Text | Google Scholar

45. Zhang C, Hediger ML, Albert PS, Grewal J, Sciscione A, Grobman WA, et al. Association of maternal obesity with longitudinal ultrasonographic measures of fetal growth. JAMA Pediatr. (2018) 172(1):24–31. doi: 10.1001/jamapediatrics.2017.3785

PubMed Abstract | Crossref Full Text | Google Scholar

46. Krieger N. Refiguring “race”: epidemiology, racialized biology, and biological expressions of race relations. Int J Health Serv. (2000) 30(1):211–6. doi: 10.2190/672J-1PPF-K6QT-9N7U

PubMed Abstract | Crossref Full Text | Google Scholar

47. Skupski DW, Owen J, Kim S, Fuchs KM, Albert PS, Grantz KL. Estimating gestational age from ultrasound fetal biometrics. Obstet Gynecol. (2017) 130(2):433–41. doi: 10.1097/AOG.0000000000002137

PubMed Abstract | Crossref Full Text | Google Scholar

48. Ma WL, Gao C, Bell EM, Druschel CM, Caggana M, Aldous KM, et al. Analysis of polychlorinated biphenyls and organochlorine pesticides in archived dried blood spots and its application to track temporal trends of environmental chemicals in newborns. Environ Res. (2014) 133:204–10. doi: 10.1016/j.envres.2014.05.029

PubMed Abstract | Crossref Full Text | Google Scholar

49. Schisterman EF, Vexler A, Whitcomb BW, Liu A. The limitations due to exposure detection limits for regression models. Am J Epidemiol. (2006) 163(4):374–83. doi: 10.1093/aje/kwj039

PubMed Abstract | Crossref Full Text | Google Scholar

50. Pugh SJ, Albert PS, Kim S, Grobman W, Hinkle SN, Newman RB, et al. Patterns of gestational weight gain and birth weight outcomes in the NICHD Fetal Growth Study—singletons: a prospective study. Am J Obstet Gynecol. (2017) 217(3):346.e1–346.e11. doi: 10.1016/j.ajog.2017.05.013

PubMed Abstract | Crossref Full Text | Google Scholar

51. Akins JR, Waldrep K, Bernert JT. The estimation of total serum lipids by a completely enzymatic “summation” method. Clin Chim Acta. (1989) 184(3):219–26. doi: 10.1016/0009-8981(89)90054-5

PubMed Abstract | Crossref Full Text | Google Scholar

52. Bernert JT, Turner WE, Patterson DG, Needham LL. Calculation of serum “total lipid” concentrations for the adjustment of persistent organohalogen toxicant measurements in human samples. Chemosphere. (2007) 68(5):824–31. doi: 10.1016/j.chemosphere.2007.02.043

PubMed Abstract | Crossref Full Text | Google Scholar

53. Mitro SD, Chu MT, Dodson RE, Adamkiewicz G, Chie L, Brown FM, et al. Phthalate metabolite exposures among immigrants living in the United States: findings from NHANES, 1999–2014. J Expo Sci Environ Epidemiol. (2019) 29(1):71–82. doi: 10.1038/s41370-018-0029-x

PubMed Abstract | Crossref Full Text | Google Scholar

54. Vanderweele TJ, Vansteelandt S. Conceptual issues concerning mediation, interventions and composition. Stat Interface. (2009) 2(4):457–68. doi: 10.4310/SII.2009.v2.n4.a7

Crossref Full Text | Google Scholar

55. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of weighted quantile sum regression for highly correlated data in a risk analysis setting. J Agric Biol Environ Stat. (2015) 20(1):100–21. doi: 10.1007/s13253-014-0180-3

PubMed Abstract | Crossref Full Text | Google Scholar

56. Valeri L, VanderWeele TJ. Mediation analysis allowing for exposure–mediator interactions and causal interpretation: theoretical assumptions and implementation with SAS and SPSS macros. Psychol Methods. (2013) 18(2):137–50. doi: 10.1037/a0031034

PubMed Abstract | Crossref Full Text | Google Scholar

57. Bobb JF, Valeri L, Claus Henn B, Christiani DC, Wright RO, Mazumdar M, et al. Bayesian Kernel machine regression for estimating the health effects of multi-pollutant mixtures. Biostatistics. (2015) 16(3):493–508. doi: 10.1093/biostatistics/kxu058

PubMed Abstract | Crossref Full Text | Google Scholar

58. VanderWeele TJ, Vansteelandt S. Mediation analysis with multiple mediators. Epidemiol Method. (2014) 2(1):95–115. doi: 10.1515/em-2012-0010

PubMed Abstract | Crossref Full Text | Google Scholar

59. VanderWeele TJ. A unification of mediation and interaction: a 4-way decomposition. Epidemiology. (2014) 25(5):749–61. doi: 10.1097/EDE.0000000000000121

PubMed Abstract | Crossref Full Text | Google Scholar

60. Naimi AI, Cole SR, Kennedy EH. An introduction to g methods. Int J Epidemiol. (2017) 46(2):756–62. doi: 10.1093/ije/dyw323

PubMed Abstract | Crossref Full Text | Google Scholar

61. Robins J. A new approach to causal inference in mortality studies with a sustained exposure period—application to control of the healthy worker survivor effect. Mathemat Modell. (1986) 7(9):1393–512. doi: 10.1016/0270-0255(86)90088-6

Crossref Full Text | Google Scholar

62. Shi B, Choirat C, Coull BA, VanderWeele TJ, Valeri L. CMAverse: a suite of functions for reproducible causal mediation analyses. Epidemiology. (2021) 32(5):e20. doi: 10.1097/EDE.0000000000001378

PubMed Abstract | Crossref Full Text | Google Scholar

63. Keil AP, Buckley JP, O ‘Brien Katie M, Ferguson KK, Zhao S, White AJ. A quantile-based g-computation approach to addressing the effects of exposure mixtures. Environ Health Perspect. (2020) 128(4):047004. doi: 10.1289/EHP5838

PubMed Abstract | Crossref Full Text | Google Scholar

64. Yim G, Minatoya M, Kioumourtzoglou MA, Bellavia A, Weisskopf M, Ikeda-Araki A, et al. The associations of prenatal exposure to dioxins and polychlorinated biphenyls with neurodevelopment at 6 months of age: multi-pollutant approaches. Environ Res. (2022) 209:112757. doi: 10.1016/j.envres.2022.112757

PubMed Abstract | Crossref Full Text | Google Scholar

65. Valeri L, Lin X, VanderWeele TJ. Mediation analysis when a continuous mediator is measured with error and the outcome follows a generalized linear model. Stat Med. (2014) 33(28):4875–90. doi: 10.1002/sim.6295

PubMed Abstract | Crossref Full Text | Google Scholar

66. Renzetti S, Gennings C, Curtin P. gWQS: an R Package for Linear and Generalized Weighted Quantile Sum (WQS) Regression. (2020). Available at: https://www.semanticscholar.org/paper/gWQS%3A-An-R-Package-for-Linear-and-Generalized-Sum-Renzetti-Gennings/dade18b42723a8f9750fa5831ac06545942c02a3

67. Dongarwar D, Tahseen D, Wang L, Aliyu MH, Salihu HM. Temporal trends in preterm birth phenotypes by plurality: Black–White disparity over half a century. J Perinatol. (2021) 41(2):204–11. doi: 10.1038/s41372-020-00912-8

PubMed Abstract | Crossref Full Text | Google Scholar

68. Zota AR, Adamkiewicz G, Morello-Frosch RA. Are PBDEs an environmental equity concern? Exposure disparities by socioeconomic status. Environ Sci Technol. (2010) 44(15):5691–2. doi: 10.1021/es101723d

PubMed Abstract | Crossref Full Text | Google Scholar

69. Liu R, Nelson DO, Hurley S, Petreas M, Park JS, Wang Y, et al. Association between serum polybrominated diphenylether levels and residential proximity to solid waste facilities. Environ Sci Technol. (2016) 50(7):3945–53. doi: 10.1021/acs.est.5b04715

PubMed Abstract | Crossref Full Text | Google Scholar

70. Zota AR, Rudel RA, Morello-Frosch RA, Brody JG. Elevated house dust and serum concentrations of PBDEs in California: unintended consequences of furniture flammability standards? Environ Sci Technol. (2008) 42(21):8158–64. doi: 10.1021/es801792z

PubMed Abstract | Crossref Full Text | Google Scholar

71. Harville EW, Knoepp LR, Wallace ME, Miller KS. Cervical pathways for racial disparities in preterm births: the preterm prediction study. J Matern Fetal Neonatal Med. (2019) 32(23):4022–8. doi: 10.1080/14767058.2018.1484091

PubMed Abstract | Crossref Full Text | Google Scholar

72. Kramer MR, Cooper HL, Drews-Botsch CD, Waller LA, Hogue CR. Metropolitan isolation segregation and black-white disparities in very preterm birth: a test of mediating pathways and variance explained. Soc Sci Med. (2010) 71(12):2108–16. doi: 10.1016/j.socscimed.2010.09.011

PubMed Abstract | Crossref Full Text | Google Scholar

73. Hedderson MM, Xu F, Dayo OM, Liu E, Sridhar S, Lee C, et al. Contribution of maternal cardiometabolic risk factors to racial-ethnicity disparities in preterm birth subtypes. Am J Obstet Gynecology MFM. (2022) 4(3):100608. doi: 10.1016/j.ajogmf.2022.100608

Crossref Full Text | Google Scholar

74. Singh GK, Yu SM. Adverse pregnancy outcomes: differences between US- and foreign-born women in major US racial and ethnic groups. Am J Public Health. (1996) 86(6):837–43. doi: 10.2105/AJPH.86.6.837

PubMed Abstract | Crossref Full Text | Google Scholar

75. Casey JA, Daouda M, Babadi RS, Do V, Flores NM, Berzansky I, et al. Methods in public health environmental justice research: a scoping review from 2018 to 2021. Curr Environ Health Rep. (2023) 10(3):312–36. doi: 10.1007/s40572-023-00406-7

PubMed Abstract | Crossref Full Text | Google Scholar