Abstract
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
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/m2], 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/m2), 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: (when Y represents continuous gestational age at delivery, in weeks), or (where represents PTB and 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: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: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:along with six separate regression models estimating each mediator as a function of the exposure: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
| Overall | Non-Hispanic White | Non-Hispanic Black | Hispanic | Asian/Pacific Islander | |
|---|---|---|---|---|---|
| (n = 2,008) | (n = 552) | (n = 504) | (n = 568) | (n = 384) | |
| Age (years) | 28.3 ± 5.4 | 30.3 ± 4.4 | 25.6 ± 5.5 | 27.1 ± 5.5 | 30.6 ± 4.5 |
| Prepregnancy BMI (kg/m2) | 23.6 ± 3.0 | 23.3 ± 2.8 | 24.2 ± 3.1 | 24.4 ± 2.8 | 22.2 ± 2.6 |
| Parity, n (%) | |||||
| 0 | 979 (48.8) | 299 (54.2) | 253 (50.2) | 223 (39.3) | 204 (53.1) |
| 1 | 689 (34.3) | 184 (33.3) | 157 (31.2) | 204 (35.9) | 144 (37.5) |
| 2+ | 340 (16.9) | 69 (12.5) | 94 (18.7) | 141 (24.8) | 36 (9.4) |
| Education level | |||||
| Less than college degree | 554 (27.6) | 30 (5.4) | 191 (37.9) | 267 (47.0) | 66 (17.2) |
| Some college or undergraduate | 1,086 (54.1) | 336 (60.9) | 271 (53.8) | 282 (49.6) | 197 (51.3) |
| Graduate or postgraduate | 368 (18.3) | 186 (33.7) | 42 (8.3) | 19 (3.3) | 121 (31.5) |
| Marital statusa, n (%) | |||||
| Married or living with partner | 1,539 (76.6) | 518 (93.8) | 251 (49.8) | 417 (73.4) | 353 (91.9) |
| Not married | 467 (23.3) | 33 (6.0) | 252 (50.0) | 151 (26.6) | 31 (8.1) |
| Family income during last year, n (%) | |||||
| Less than $30,000 | 470 (23.4) | 21 (3.8) | 205 (40.7) | 196 (34.5) | 48 (12.5) |
| $30,000–$49,999 | 288 (14.3) | 40 (7.2) | 90 (17.9) | 125 (22.0) | 33 (8.6) |
| $50,000–$99,999 | 451 (22.5) | 166 (30.1) | 86 (17.1) | 95 (16.7) | 104 (27.1) |
| $100,000 or more | 523 (26.0) | 305 (55.3) | 57 (11.3) | 54 (9.5) | 107 (27.9) |
| Unknown | 276 (13.7) | 20 (3.6) | 66 (13.1) | 98 (17.3) | 92 (24.0) |
| Plasma cotinine (ng/mL)a | 1.1 ± 12.7 | 1.1 ± 13.2 | 2.7 ± 20.8 | 0.3 ± 4.4 | 0.0 ± 0.2 |
| Plasma total lipids (non-fasting) (mg/dL)a,b | 610.5 ± 98.7 | 613.3 ± 95.9 | 580.8 ± 99.0 | 628.3 ± 100.4 | 619.2 ± 91.2 |
| Total activity (MET hours per week)a | 323.1 ± 167.8 | 326.0 ± 147.8 | 354.8 ± 200.6 | 307.3 ± 158.9 | 300.3 ± 153.9 |
| Sedentary activity (MET hours per week)a | 26.1 ± 18.3 | 20.7 ± 12.2 | 37.1 ± 22.2 | 22.6 ± 17.0 | 24.4 ± 15.6 |
| Acculturationa, n (%) | |||||
| US-born | 1,322 (65.8) | 514 (93.1) | 462 (91.7) | 242 (42.6) | 104 (27.1) |
| Recent immigrant (<10 years) | 304 (15.1) | 16 (2.9) | 16 (3.2) | 117 (20.6) | 155 (40.4) |
| Long-term immigrant (≥10 years) | 379 (18.9) | 20 (3.6) | 26 (5.2) | 208 (36.6) | 125 (32.6) |
| Gestational age at delivery (weeks) | 39.2 ± 1.7 | 39.3 ± 1.5 | 39.0 ± 2.1 | 39.3 ± 1.5 | 39.3 ± 1.3 |
| Preterm birth, n (%) | 118 (5.9) | 28 (5.1) | 46 (9.1) | 26 (4.6) | 18 (4.7) |
Characteristics of the study population by race and ethnicity, NICHD Fetal Growth Study–Singleton Cohort (n = 2,008).
Means ± SD for continuous variables. N (%) for categorical variables.
Numbers may not add up to total numbers owing to missing values. Variables with missing values (missing rate) included: marital status (0.1%), plasma cotinine (1.6%), plasma total lipids (1.1%), AHEI 2010 score (37.4%), total activity (0.2%), sedentary activity (0.2%), and acculturation (0.1%).
Total lipids = total cholesterol × 2.27 + triglycerides + 62.3.
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
| Chemicals, ng/g lipid | Race and ethnicity | PTB status | ||||
|---|---|---|---|---|---|---|
| Non-Hispanic White | Non-Hispanic Black | Hispanic | Asian/Pacific Islander | Preterm | Non-preterm | |
| (n = 552) | (n = 504) | (n = 568) | (n = 384) | (n = 118) | (n = 1,890) | |
| PBDE 28 | 0.20 (0.18, 0.22) | 0.29 (0.26, 0.33) | 0.30 (0.26, 0.33) | 0.25 (0.22, 0.29) | 0.28 (0.22, 0.37) | 0.25 (0.24, 0.27) |
| PBDE 47 | 4.48 (3.77, 5.33) | 8.89 (7.49, 10.54) | 5.53 (4.55, 6.72) | 3.62 (2.88, 4.55) | 6.10 (4.12, 9.03) | 5.38 (4.87, 5.94) |
| PBDE 99 | 0.18 (0.13, 0.25) | 0.52 (0.38, 0.73) | 0.43 (0.32, 0.58) | 0.24 (0.16, 0.34) | 0.29 (0.14, 0.58) | 0.32 (0.27, 0.38) |
| PBDE 100 | 0.55 (0.46, 0.67) | 1.44 (1.19, 1.75) | 1.01 (0.85, 1.21) | 0.52 (0.41, 0.65) | 0.84 (0.54, 1.30) | 0.82 (0.74, 0.91) |
| PBDE 153 | 1.37 (1.16, 1.61) | 1.90 (1.61, 2.25) | 0.83 (0.73, 0.94) | 0.83 (0.71, 0.98) | 1.35 (0.95, 1.93) | 1.16 (1.07, 1.26) |
| PBDE 154 | 0.05 (0.04, 0.07) | 0.11 (0.08, 0.15) | 0.07 (0.05, 0.09) | 0.12 (0.09, 0.17) | 0.07 (0.04, 0.14) | 0.08 (0.07, 0.09) |
| ∑PBDEs, pmol/g lipida | 27.8 (25.2, 30.8) | 47.4 (42.8, 52.5) | 32.7 (29.6, 36.1) | 23.4 (20.6, 26.5) | 36.6 (29.4, 45.6) | 31.9 (30.2, 33.7) |
Geometric means (95% confidence intervals) of lipid-adjusted PBDE concentrations, stratified by race and ethnicity and by preterm birth status.
Geometric means (95% CI) calculated using all machine-observed values including those below the LOQs, where zero and negative values were assigned the value of (lowest positive value)/2. PTB, preterm birth.
∑PBDEs (in pmol/g lipid) refers to the molar sum of PBDEs, which was calculated by dividing each lipid-adjusted chemical concentration by its molecular weight and summing all detectable concentrations.
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 ORTE (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
Table 3
| Race and ethnicity | Adjusteda β (95% CIb) for gestational age at delivery, weeks | ||||||
|---|---|---|---|---|---|---|---|
| Natural direct effect (βNDE) | Natural indirect effect (βNIE) | Controlled direct effects (CDEs), fixing the WQS index at the 25th percentile, median, and 75th percentile | Total effect(βTE) | Proportion mediated (PM), % | |||
| βCDE(25th) | βCDE(median) | βCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | −0.26 (−0.51, −0.01) | −0.04 (−0.07, 0.00) | −0.18 (−0.46, 0.10) | −0.23 (−0.49, 0.02) | −0.32 (−0.57, −0.07) | −0.30 (−0.54, −0.05) | 11.8% (−3.8%, 27.5%) |
| Hispanic | 0.07 (−0.18, 0.32) | −0.01 (−0.03, 0.01) | 0.12 (−0.16, 0.40) | 0.09 (−0.17, 0.34) | 0.03 (−0.23, 0.29) | 0.06 (−0.19, 0.31) | −20.7% (−123.7%, 82.2%) |
| Asian/Pacific Islander | −0.05 (−0.32, 0.21) | −0.01 (−0.03, 0.02) | −0.02 (−0.32, 0.28) | −0.04 (−0.31, 0.23) | −0.08 (−0.35, 0.20) | −0.06 (−0.32, 0.20) | 13.1% (−50.6%, 76.7%) |
| Adjusteda OR (95% CIb) for preterm birth | |||||||
| Natural direct effect (ORNDE) | Natural indirect effect (ORNIE) | Controlled direct effects (CDEs), fixing the WQS index at the 25th percentile, median, and 75th percentile | Total effect (ORTE) | Proportion mediated (PM), % | |||
| ORCDE(25th) | ORCDE(median) | ORCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | 1.79 (0.98, 3.25) | 1.02 (0.98, 1.06) | 1.77 (0.89, 3.52) | 1.78 (0.95, 3.32) | 1.79 (0.98, 3.25) | 1.82 (1.00, 3.31) | 4.0% (−5.1%, 13.1%) |
| Hispanic | 0.89 (0.45, 1.75) | 1.01 (0.99, 1.03) | 0.83 (0.38, 1.83) | 0.85 (0.42, 1.73) | 0.89 (0.45, 1.75) | 0.90 (0.46, 1.77) | −7.4% (−60.1%, 45.2%) |
| Asian/Pacific Islander | 0.97 (0.47, 2.01) | 1.01 (0.98, 1.03) | 0.81 (0.35, 1.87) | 0.88 (0.41, 1.87) | 0.96 (0.47, 1.99) | 0.97 (0.47, 2.01) | −29.6% (−807.4%, 748.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.
Adjusted for maternal age (years), prepregnancy BMI (kg/m2), 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 or more, not reported), plasma cotinine level (ng/mL), plasma total lipids (ng/mL), total activity (MET hours per week), sedentary activity (MET hours per week), and acculturation (US born, recent immigrant, long-term immigrant). Observations with missing covariates were excluded from the adjusted models.
Standard errors for calculating the 95% CIs obtained using the delta method, based on point estimates obtained using closed-form parameter function estimation.
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
| Race and ethnicity | Adjusteda β (95% CIb) for gestational age at delivery, weeks | ||||||
|---|---|---|---|---|---|---|---|
| Natural direct effect (βNDE) | Natural indirect effect (βNIE) | Controlled direct effects (CDEs), fixing PBDE 153 concentration at the 25th percentile, median, and 75th percentile | Total effect (βTE) | Proportion mediated (PM), % | |||
| βCDE(25th) | βCDE(median) | βCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | −0.27 (−0.52, −0.03) | −0.02 (−0.06, 0.01) | −0.18 (−0.45, 0.09) | −0.18 (−0.45, 0.09) | −0.29 (−0.54, −0.05) | −0.30 (−0.54, −0.05) | 7.9% (−5.4%, 21.1%) |
| Hispanic | 0.05 (−0.20, 0.30) | 0.01 (−0.01, 0.03) | 0.07 (−0.21, 0.34) | 0.07 (−0.21, 0.34) | 0.05 (−0.21, 0.30) | 0.06 (−0.19, 0.31) | 18.4% (−52.2%, 77.6%) |
| Asian/Pacific Islander | −0.07 (−0.33, 0.19) | 0.01 (−0.01, 0.02) | −0.06 (−0.35, 0.23) | −0.06 (−0.35, 0.23) | −0.07 (−0.33, 0.19) | −0.06 (−0.32, 0.20) | −9.0% (−57.8%, 39.8%) |
| Adjusteda OR (95% CIb) for preterm birth (cutoff: 37 weeks) | |||||||
| Natural direct effect (ORNDE) | Natural indirect effect (ORNIE) | Controlled direct effects (CDEs), fixing PBDE 153 concentration at the 25th percentile, median, and 75th percentile | Total effect (ORTE) | Proportion mediated (PM), % | |||
| ORCDE(25th) | ORCDE(median) | ORCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | 1.81 (0.99, 3.32) | 1.01 (0.98, 1.04) | 1.92 (0.97, 3.79) | 1.92 (0.97, 3.79) | 1.81 (0.99, 3.31) | 1.82 (1.00, 3.33) | 1.3% (−5.0%, 7.6%) |
| Hispanic | 0.89 (0.45, 1.78) | 0.99 (0.94, 1.04) | 0.91 (0.42, 1.97) | 0.91 (0.42, 1.97) | 0.89 (0.45, 1.78) | 0.88 (0.45, 1.76) | 7.3% (−50.4%, 64.9%) |
| Asian/Pacific Islander | 0.98 (0.48, 2.03) | 0.99 (0.97, 1.01) | 0.90 (0.39, 2.07) | 0.90 (0.39, 2.07) | 0.98 (0.48, 2.03) | 0.97 (0.47, 2.01) | 33.9% (−911.4%, 979.2%) |
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.
Adjusted for maternal age (years), prepregnancy BMI (kg/m2), 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, or more, not reported), plasma cotinine level (ng/mL), plasma total lipids (ng/mL), total activity (MET hours per week), sedentary activity (MET hours per week), and acculturation (US born, recent immigrant, long-term immigrant). Observations with missing covariates were excluded from the adjusted models.
Standard errors for calculating the 95% CIs obtained using the delta method, based on point estimates obtained using closed-form parameter function estimation.
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
| Race and ethnicity | Adjusteda β (95% CIb) for gestational age at delivery, weeks | ||||||
|---|---|---|---|---|---|---|---|
| Natural direct effect (βNDE) | Natural indirect effect (βNIE) | Controlled direct effects (CDEs), fixing all chemical concentrations at the 25th percentile, median, and 75th percentile | Total effectc (βTE) | Proportion mediated (PM), % | |||
| βCDE(25th) | βCDE(median) | βCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | −0.25 (−0.54, −0.02) | −0.05 (−0.16, 0.05) | −0.25 (−0.53, 0.02) | −0.14 (−0.47, 0.17) | −0.27 (−0.61, −0.03) | −0.29 (−0.57, −0.11) | 16.3% (−16.2%, 81.1%) |
| Hispanic | 0.05 (−0.19, 0.27) | 0.01 (−0.05, 0.08) | 0.20 (−0.07, 0.46) | 0.23 (−0.02, 0.52) | −0.05 (−0.33, 0.21) | 0.06 (−0.17, 0.26) | 10.4% (−233.6%, 335.5%) |
| Asian/Pacific Islander | −0.04 (−0.29, 0.20) | −0.03 (−0.14, 0.05) | 0.06 (−0.28, 0.36) | 0.12 (−0.22, 0.43) | −0.13 (−0.43, 0.15) | −0.07 (−0.39, 0.18) | 44.8% (−203.2%, 212.7%) |
| Adjusteda OR (95% CIb) for preterm birth | |||||||
| Natural direct effect (ORNDE) | Natural indirect effect (ORNIE) | Controlled direct effects (CDEs), Fixing all chemical concentrations at the 25th percentile, median, and 75th percentile | Total effectc (ORTE) | Proportion mediated (PM), % | |||
| ORCDE(25th) | ORCDE(median) | ORCDE(75th) | |||||
| Non-Hispanic White | REF | REF | REF | REF | REF | REF | REF |
| Non-Hispanic Black | 1.66 (0.93, 3.16) | 1.04 (0.93, 1.15) | 1.66 (0.93, 3.21) | 1.66 (0.93, 3.20) | 1.67 (0.93, 3.19) | 1.73 (1.03, 3.04) | 9.7% (−27.8%, 47.3%) |
| Hispanic | 0.89 (0.51, 1.75) | 1.01 (0.52, 1.74) | 0.89 (0.51, 1.76) | 0.89 (0.51, 1.76) | 0.89 (0.51, 1.75) | 0.90 (0.52, 1.74) | −8.7% (−150.4%, 170.8%) |
| Asian/Pacific Islander | 0.97 (0.41, 1.93) | 1.00 (0.92, 1.09) | 0.97 (0.41, 1.95) | 0.97 (0.41, 1.95) | 0.97 (0.40, 1.94) | 0.97 (0.41, 2.02) | 4.8% (−100.2%, 103.3%) |
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.
Adjusted for maternal age (years), pre-pregnancy BMI (kg/m2), 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, or more, not reported), plasma cotinine level (ng/ml), plasma total lipids (ng/mL), total activity (MET hours per week), sedentary activity (MET hours per week), and acculturation (US born, recent immigrant, long-term immigrant). Observations with missing covariates were excluded from the adjusted models.
95% CIs obtained using the bootstrapping method, based on point estimates obtained using direct counterfactual imputation estimation.
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.
Statements
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.
MartinJAOstermanMJK. Describing the increase in preterm births in the United States, 2014–2016. NCHS Data Brief. (2018) (312):1–8. PMID: 30044213
2.
Institute of Medicine (US) Committee on Understanding Premature Birth and Assuring Healthy Outcomes. Preterm birth: causes, consequences, and prevention. In: BehrmanREButlerAS, 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.
BurrisHHCollinsJWWrightRO. Racial/ethnic disparities in preterm birth: clues from environmental exposures. Curr Opin Pediatr. (2011) 23(2):227–32. 10.1097/MOP.0b013e328344568f
4.
BurrisHHLorchSAKirpalaniHPursleyDMElovitzMACloughertyJE. Racial disparities in preterm birth in USA: a biosensor of physical and social environmental exposures. Arch Dis Child. (2019) 104(10):931–5. 10.1136/archdischild-2018-316486
5.
SchaafJLiemSMolBAbu-HannaARavelliA. Ethnic and racial disparities in the risk of preterm birth: a systematic review and meta-analysis. Amer J Perinatol. (2012) 30(6):433–50. 10.1055/s-0032-1326988
6.
AnumEASpringelEHShriverMDStraussJF. Genetic contributions to disparities in preterm birth. Pediatr Res. (2009) 65(1):1–9. 10.1203/PDR.0b013e31818912e7
7.
BravemanPDominguezTPBurkeWDolanSMStevensonDKJacksonFMet alExplaining 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. 10.3389/frph.2021.684207
8.
BurrisHHWrightCJKirpalaniHCollinsJWJr.LorchSAElovitzMAet alThe promise and pitfalls of precision medicine to resolve black–white racial disparities in preterm birth. Pediatr Res. (2020) 87(2):221–6. 10.1038/s41390-019-0528-z
9.
ThomaMEDrewLBHiraiAHKimTYFenelonAShenassaED. Black–white disparities in preterm birth: geographic, social, and health determinants. Am J Prev Med. (2019) 57(5):675–86. 10.1016/j.amepre.2019.07.007
10.
JohnsonJDGreenCAVladutiuCJManuckTA. Racial disparities in prematurity persist among women of high socioeconomic Status. Am J Obstet Gynecol MFM. (2020) 2(3):100104. 10.1016/j.ajogmf.2020.100104
11.
ThurstonHFieldsBEWhiteJ. Does increasing access to prenatal care reduce racial disparities in birth outcomes?J Pediatr Nurs. (2021) 59:96–102. 10.1016/j.pedn.2021.01.012
12.
TimofeevJSinghJIstwanNRheaDDriggersRW. Spontaneous preterm birth in African-American and Caucasian women receiving 17α-hydroxyprogesterone caproate. Am J Perinatol. (2014) 31(1):55–60. 10.1055/s-0033-1334452
13.
LiuWYYuZBQiuHYWangJBChenXYChenK. Association between ambient air pollutants and preterm birth in Ningbo, China: a time-series study. BMC Pediatr. (2018) 18(1):305. 10.1186/s12887-018-1282-9
14.
AshMBoyceJK. Racial disparities in pollution exposure and employment at US industrial facilities. Proc Natl Acad Sci U.S.A. (2018) 115(42):10636–41. 10.1073/pnas.1721640115
15.
MullerCSampsonRJWinterAS. Environmental inequality: the social causes and consequences of lead exposure. Annu Rev Sociol. (2018) 44(1):263–82. 10.1146/annurev-soc-073117-041222
16.
AndrewsKWSavitzDAHertz-PicciottoI. 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. 10.1002/ajim.4700260103
17.
FergusonKKMcElrathTFMeekerJD. Environmental phthalate exposure and preterm birth. JAMA Pediatr. (2014) 168(1):61–7. 10.1001/jamapediatrics.2013.3699
18.
SilvaMJBarrDBReidyJAMalekNAHodgeCCCaudillSPet alUrinary 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. 10.1289/ehp.6723
19.
DaoudaMHennemanLKioumourtzoglouMAGemmillAZiglerCCaseyJA. 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. 10.1088/1748-9326/abe4f7
20.
HarrisSMColacinoJBuxtonMCroxtonLNguyenVLoch-CarusoRet alA 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. 10.1007/s43032-022-00870-w
21.
VanderWeeleTJ. Mediation analysis: a practitioner’s guide. Annu Rev Public Health. (2016) 37:17–32. 10.1146/annurev-publhealth-032315-021402
22.
James-ToddTMChiuYHZotaAR. 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. 10.1007/s40471-016-0073-9
23.
NaimiAISchnitzerMEMoodieEEMBodnarLM. Mediation analysis for health disparities research. Am J Epidemiol. (2016) 184(4):315–24. 10.1093/aje/kwv329
24.
SchulzANorthridgeME. Social determinants of health: implications for environmental health promotion. Health Educ Behav. (2004) 31(4):455–71. 10.1177/1090198104265598
25.
AylwardLLKirmanCRSchoenyRPortierCJHaysSM. 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. 10.1289/ehp.1205740
26.
KortenkampAFaustMScholzeMBackhausT. Low-level exposure to multiple chemicals: reason for human health concerns?Environ Health Perspect. (2007) 115(Suppl 1):106–14. 10.1289/ehp.9358
27.
BellaviaAZotaARValeriLJames-ToddT. Multiple mediators approach to study environmental chemicals as determinants of health disparities. Environ Epidemiol. (2018) 2(2):e015. 10.1097/EE9.0000000000000015
28.
BellaviaAJames-ToddTWilliamsPL. Approaches for incorporating environmental mixtures as mediators in mediation analysis. Environ Int. (2019) 123:368–74. 10.1016/j.envint.2018.12.024
29.
BellaviaAValeriL. Decomposition of the total effect in the presence of multiple mediators and interactions. Am J Epidemiol. (2018) 187(6):1311–8. 10.1093/aje/kwx355
30.
TaylorKWJoubertBRBraunJMDilworthCGenningsCHauserRet alStatistical approaches for assessing health effects of environmental chemical mixtures in epidemiology: lessons from an innovative workshop. Environ Health Perspect. (2016) 124(12):A227–9. 10.1289/EHP547
31.
EskenaziBChevrierJRauchSAKogutKHarleyKGJohnsonCet alIn utero and childhood polybrominated diphenyl ether (PBDE) exposures and neurodevelopment in the CHAMACOS study. Environ Health Perspect. (2013) 121(2):257–62. 10.1289/ehp.1205597
32.
ChevrierJHarleyKGBradmanAGharbiMSjödinAEskenaziB. Polybrominated diphenyl ether (PBDE) flame retardants and thyroid hormone during pregnancy. Environ Health Perspect. (2010) 118(10):1444–9. 10.1289/ehp.1001905
33.
Johnson-RestrepoBKannanK. An assessment of sources and pathways of human exposure to polybrominated diphenyl ethers in the United States. Chemosphere. (2009) 76(4):542–8. 10.1016/j.chemosphere.2009.02.068
34.
AllenJGMcCleanMDStapletonHMWebsterTF. Linking PBDEs in house dust to consumer products using x-ray fluorescence. Environ Sci Technol. (2008) 42(11):4222–8. 10.1021/es702964a
35.
Buck LouisGMZhaiSSmarrMMGrewalJZhangCGrantzKLet alEndocrine disruptors and neonatal anthropometry, NICHD Fetal Growth Studies—singletons. Environ Int. (2018) 119:515–26. 10.1016/j.envint.2018.07.024
36.
NguyenVKKahanaAHeidtJPolemiKKvasnickaJJollietOet alA comprehensive analysis of racial disparities in chemical biomarker concentrations in United States women, 1999–2014. Environ Int. (2020) 137:105496. 10.1016/j.envint.2020.105496
37.
VarshavskyJRSenSRobinsonJFSmithSCFrankenfieldJWangYet alRacial/ethnic and geographic differences in polybrominated diphenyl ether (PBDE) levels across maternal, placental, and fetal tissues during mid-gestation. Sci Rep. (2020) 10:12247. 10.1038/s41598-020-69067-y
38.
HarleyKGChevrierJAguilar SchallRSjödinABradmanAEskenaziB. Association of prenatal exposure to polybrominated diphenyl ethers and infant birth weight. Am J Epidemiol. (2011) 174(8):885–92. 10.1093/aje/kwr212
39.
EickSMHom ThepaksornEKIzanoMACushingLJWangYSmithSCet alAssociations 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. 10.1186/s12940-020-00654-2
40.
BehniaFPeltierMRSaadeGRMenonR. Environmental pollutant polybrominated diphenyl ether, a flame retardant, induces primary amnion cell senescence. Am J Reprod Immunol. (2015) 74(5):398–406. 10.1111/aji.12414
41.
GaoYChenLWangCZhouYWangYZhangYet alExposure to polybrominated diphenyl ethers and female reproductive function: a study in the production area of Shandong, China. Sci Total Environ. (2016) 572:9–15. 10.1016/j.scitotenv.2016.07.181
42.
WangZZhangCWilliamsPLBellaviaAWylieBJHackerMRet alPolybrominated diphenyl ethers in early pregnancy and preterm birth: findings from the NICHD Fetal Growth Studies. Int J Hyg Environ Health. (2022) 243:113978. 10.1016/j.ijheh.2022.113978
43.
HoweCJBaileyZDRaifmanJRJacksonJW. Recommendations for using causal diagrams to study racial health disparities. Am J Epidemiol. (2022) 191(12):1981–9. 10.1093/aje/kwac140
44.
GrewalJGrantzKLZhangCSciscioneAWingDAGrobmanWAet alCohort profile: NICHD Fetal Growth Studies—singletons and twins. Int J Epidemiol. (2018) 47(1):25–25l. 10.1093/ije/dyx161
45.
ZhangCHedigerMLAlbertPSGrewalJSciscioneAGrobmanWAet alAssociation of maternal obesity with longitudinal ultrasonographic measures of fetal growth. JAMA Pediatr. (2018) 172(1):24–31. 10.1001/jamapediatrics.2017.3785
46.
KriegerN. Refiguring “race”: epidemiology, racialized biology, and biological expressions of race relations. Int J Health Serv. (2000) 30(1):211–6. 10.2190/672J-1PPF-K6QT-9N7U
47.
SkupskiDWOwenJKimSFuchsKMAlbertPSGrantzKL. Estimating gestational age from ultrasound fetal biometrics. Obstet Gynecol. (2017) 130(2):433–41. 10.1097/AOG.0000000000002137
48.
MaWLGaoCBellEMDruschelCMCagganaMAldousKMet alAnalysis 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. 10.1016/j.envres.2014.05.029
49.
SchistermanEFVexlerAWhitcombBWLiuA. The limitations due to exposure detection limits for regression models. Am J Epidemiol. (2006) 163(4):374–83. 10.1093/aje/kwj039
50.
PughSJAlbertPSKimSGrobmanWHinkleSNNewmanRBet alPatterns 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. 10.1016/j.ajog.2017.05.013
51.
AkinsJRWaldrepKBernertJT. The estimation of total serum lipids by a completely enzymatic “summation” method. Clin Chim Acta. (1989) 184(3):219–26. 10.1016/0009-8981(89)90054-5
52.
BernertJTTurnerWEPattersonDGNeedhamLL. Calculation of serum “total lipid” concentrations for the adjustment of persistent organohalogen toxicant measurements in human samples. Chemosphere. (2007) 68(5):824–31. 10.1016/j.chemosphere.2007.02.043
53.
MitroSDChuMTDodsonREAdamkiewiczGChieLBrownFMet alPhthalate metabolite exposures among immigrants living in the United States: findings from NHANES, 1999–2014. J Expo Sci Environ Epidemiol. (2019) 29(1):71–82. 10.1038/s41370-018-0029-x
54.
VanderweeleTJVansteelandtS. Conceptual issues concerning mediation, interventions and composition. Stat Interface. (2009) 2(4):457–68. 10.4310/SII.2009.v2.n4.a7
55.
CarricoCGenningsCWheelerDCFactor-LitvakP. 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. 10.1007/s13253-014-0180-3
56.
ValeriLVanderWeeleTJ. 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. 10.1037/a0031034
57.
BobbJFValeriLClaus HennBChristianiDCWrightROMazumdarMet alBayesian Kernel machine regression for estimating the health effects of multi-pollutant mixtures. Biostatistics. (2015) 16(3):493–508. 10.1093/biostatistics/kxu058
58.
VanderWeeleTJVansteelandtS. Mediation analysis with multiple mediators. Epidemiol Method. (2014) 2(1):95–115. 10.1515/em-2012-0010
59.
VanderWeeleTJ. A unification of mediation and interaction: a 4-way decomposition. Epidemiology. (2014) 25(5):749–61. 10.1097/EDE.0000000000000121
60.
NaimiAIColeSRKennedyEH. An introduction to g methods. Int J Epidemiol. (2017) 46(2):756–62. 10.1093/ije/dyw323
61.
RobinsJ. 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. 10.1016/0270-0255(86)90088-6
62.
ShiBChoiratCCoullBAVanderWeeleTJValeriL. CMAverse: a suite of functions for reproducible causal mediation analyses. Epidemiology. (2021) 32(5):e20. 10.1097/EDE.0000000000001378
63.
KeilAPBuckleyJPO ‘Brien KatieMFergusonKKZhaoSWhiteAJ. A quantile-based g-computation approach to addressing the effects of exposure mixtures. Environ Health Perspect. (2020) 128(4):047004. 10.1289/EHP5838
64.
YimGMinatoyaMKioumourtzoglouMABellaviaAWeisskopfMIkeda-ArakiAet alThe associations of prenatal exposure to dioxins and polychlorinated biphenyls with neurodevelopment at 6 months of age: multi-pollutant approaches. Environ Res. (2022) 209:112757. 10.1016/j.envres.2022.112757
65.
ValeriLLinXVanderWeeleTJ. 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. 10.1002/sim.6295
66.
RenzettiSGenningsCCurtinP. 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.
DongarwarDTahseenDWangLAliyuMHSalihuHM. Temporal trends in preterm birth phenotypes by plurality: Black–White disparity over half a century. J Perinatol. (2021) 41(2):204–11. 10.1038/s41372-020-00912-8
68.
ZotaARAdamkiewiczGMorello-FroschRA. Are PBDEs an environmental equity concern? Exposure disparities by socioeconomic status. Environ Sci Technol. (2010) 44(15):5691–2. 10.1021/es101723d
69.
LiuRNelsonDOHurleySPetreasMParkJSWangYet alAssociation between serum polybrominated diphenylether levels and residential proximity to solid waste facilities. Environ Sci Technol. (2016) 50(7):3945–53. 10.1021/acs.est.5b04715
70.
ZotaARRudelRAMorello-FroschRABrodyJG. Elevated house dust and serum concentrations of PBDEs in California: unintended consequences of furniture flammability standards?Environ Sci Technol. (2008) 42(21):8158–64. 10.1021/es801792z
71.
HarvilleEWKnoeppLRWallaceMEMillerKS. Cervical pathways for racial disparities in preterm births: the preterm prediction study. J Matern Fetal Neonatal Med. (2019) 32(23):4022–8. 10.1080/14767058.2018.1484091
72.
KramerMRCooperHLDrews-BotschCDWallerLAHogueCR. 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. 10.1016/j.socscimed.2010.09.011
73.
HeddersonMMXuFDayoOMLiuESridharSLeeCet alContribution of maternal cardiometabolic risk factors to racial-ethnicity disparities in preterm birth subtypes. Am J Obstet Gynecology MFM. (2022) 4(3):100608. 10.1016/j.ajogmf.2022.100608
74.
SinghGKYuSM. 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. 10.2105/AJPH.86.6.837
75.
CaseyJADaoudaMBabadiRSDoVFloresNMBerzanskyIet alMethods in public health environmental justice research: a scoping review from 2018 to 2021. Curr Environ Health Rep. (2023) 10(3):312–36. 10.1007/s40572-023-00406-7
Summary
Keywords
health disparities, race and ethnicity, preterm birth, gestational age, chemical stressors, polybrominated diphenyl ethers, mediation analysis, environmental mixtures
Citation
Wang Z, Zhang C, Williams PL, Bellavia A, Wylie BJ, Kannan K, Bloom MS, Hunt KJ and James-Todd T (2024) Racial and ethnic disparities in preterm birth: a mediation analysis incorporating mixtures of polybrominated diphenyl ethers. Front. Reprod. Health 5:1285444. doi: 10.3389/frph.2023.1285444
Received
30 August 2023
Accepted
12 December 2023
Published
08 January 2024
Volume
5 - 2023
Edited by
Darlene Dixon, National Institute of Environmental Health Sciences (NIH), United States
Reviewed by
Caitlin G. Howe, Dartmouth College, United States
Kaye Enid Basford, The University of Queensland, Australia
Updates
Copyright
© 2024 Wang, Zhang, Williams, Bellavia, Wylie, Kannan, Bloom, Hunt and James-Todd.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zifan Wang zwang@hsph.harvard.edu
Disclaimer
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.