Our findings suggest that maternal hardships that co-occur with gestational exposure to some OCs and metals, may interact and produce greater detrimental effects on fetal growth than either exposure alone. For instance, trans-nonachlor and oxychlordane were more strongly associated with decreased infant birth weight among females with lower income compared to those with higher income. This finding is not unexpected as it has long been recognized that females with lower income have disproportionately higher exposure for environmental stressors as well as higher risk for impaired fetal growth [47]. For instance, Borders et al. (2007) found, among their population of 1,363 pregnant American females with low income, that maternal social hardships and low birth weight were strongly related. They also reported that pregnant females with lower income are more likely to face food insecurity and consume inadequate supplemental folic acid compared to those with higher income [48]. This relationship among poverty, food insecurity, and folic acid supplementation is well supported [27, 28, 49]. In our study, 104 (5%) participants had low supplemental folic acid intake; only 10 (1%) experienced both low income and low supplemental folic acid intake. Among those infants born to females with increasing PCB 118 concentrations and low supplemental folic acid intake, we observed that every 2-fold increase in PCB 118 concentration was associated with an 87 g reduction in birth weight among females with low supplemental folic acid intake, compared to a 7 g reduction among females with the recommended supplemental folic acid intake level.
Beside maternal income, maternal education also had a strong and positive association with birth weight [50]. In our study, Pb x low education was one of the selected interaction terms by elastic net. Pb exposure is unequally distributed across populations where higher exposure to Pb is typically found in communities of lower socioeconomic status and among individuals with less access to resources including financial, educational, social, and health [51]. Furthermore, both Pb exposure and low maternal education can independently contribute to lower birth weight [6, 50, 52–54]. In combination, we found that Pb was associated with a greater reduction in birth weight among females with low educational status (\(\beta\) = -100 g; 95% CI: -215, 16) compared to those who had higher education (\(\beta\) = -34 g; 95% CI: -64, -3) (Table 4).
Another hardship that has been known to affect birth weight is immigration status. When we considered the females’ immigration status along with their exposure to OCs (i.e., trans-nonachlor and PCB 180), immigrants showed a lower mean birth weight compared with Canadian-born participants (Fig. 2A, 2C and Table 4). This birth weight reduction is supported by an increasing number of studies that found that immigration status affects fetal growth [55] and that the birth weights of babies born to immigrants are generally lower than those of babies born to Canadian-born females [56, 57]. For example, South Asian born females tend to give birth to smaller babies than non-migrant females [58]. However, in a systematic review where Gagnon et al. (2009) explored whether immigrants have poorer infant outcomes compared to Canadian-born females, the authors found that being an immigrant was not a consistent marker for poorer infant outcomes [56]. Nevertheless, among our MIREC population, immigrants had lower mean birth weights and when considering the effects of environmental chemical exposures, they experienced a slight increase in vulnerability compared with Canadian-born females. Interestingly, we also observed a greater reduction in birth weight among infants born to students with increasing concentrations of OCs (i.e., Aroclor 1260 and PCB 180), compared to those of non-students. There is a dearth of studies on environmental chemical exposures in the context of pregnant students. Available studies typically examined behaviours in reducing exposure [59], knowledge and awareness of exposures to environmental chemicals [60], or exposure assessment on campus or in laboratories [61]. Future work, therefore, should include pregnant students as the intersection of their experiences and exposure status may generate differential impacts on their pregnancy outcomes.
The stronger relationships we observed for females with both higher environmental chemical concentrations and a maternal hardship is consistent with Knudson’s “two-hit” hypothesis [62] where two distinct gestational stressors (e.g., gestational exposures to environmental chemical and maternal hardship) showed different associations with birth weight when combined compared to either “hit” alone. To further explore the impact of cumulative effects, we summed the number of hardships each participant faced and found that as the number of hardships increased, the mean birth weight decreased (Table S1). Specifically, the presence of one hardship was associated with a 59 g (95% CI: -105, -12) lower birth weight while the presence of two or more hardships was associated with an 87 g (95% CI: -139, -34) lower birth weight, compared to no hardship (Table S2). While the findings of the present study support the hypothesis that exposure to two stressors can have different effects compared to exposure to a single stressor, these findings may depend on the nature of the stressor(s). For instance, infants born to MIREC participants with increasing Mn or Hg concentrations and reported maternal hardship had higher birth weights compared to those born to participants who did not report hardship.
Studies had shown that Mn exhibits a nonlinear inverse U-shaped relationship with birth weight where lower birth weight was observed at both low and high concentrations of blood Mn level [63, 64]. Despite that, most studies, including the present one, assumed linear functions. In our previous study, using a flexible Bayesian Kernel Machine Regression method that can accommodate non-linearity, we found that Mn showed a positive and linear relationship with infant birth weight among MIREC participants [6]. It is possible that we have captured the ascending segment of the inverted U-shaped Mn-birth weight relationship, which showed an association with higher birth weight in our adjusted regression analysis. As a result, when examined by maternal hardship status, we found evidence that infants born to MIREC participants with both higher Mn concentration and low income had a higher birth weight, while those with no reported hardship showed no association with birth weight. However, the pattern doesn't hold true for lone parenthood as infants born to single females experienced approximately 200 g greater birth weight reduction compared to those born to married parents. The limited research on lone parenthood (or single mother families) and environmental chemical exposures has focused on exposure to air pollutants. Two US-based environmental inequality studies reported that single-mother families are more likely to live in neighborhoods with higher pollution compared to other family types (e.g., married or single father families) [65, 66]. Another US air pollution study found that living in a poor neighbourhood may increase the risk of exposure to OCs and metals in the air [67]. We are not aware of any study that has examined the combined effects of lone parenthood and Mn exposure on pregnancy outcomes. Research that assesses the relationship between wider ranges in Mn concentrations and fetal growth is highly recommended to consider the shape of the Mn-birth weight function as well as the impact of lone parenthood.
Like Mn, total Hg also showed a positive relationship with birth weight among those with maternal hardships. Although Hg is a toxic chemical that can freely cross the placenta and potentially disrupt a range of important pregnancy processes [68], in a recent systematic review, Hg was reported to have a minimal to null association with birth weight [69]. However, one particular study reported that the lowest tertile hair Hg concentration was associated with a higher risk of having infants with low birth weight (adjusted odds ratio (aOR) = 7.2; 95% CI: 1.5, 35.6) compared to those with higher Hg concentration (aOR = 0.52; 95% CI: 0.17, 1.55) [70]. As Hg exposure generally results from maternal fish consumption, the authors hypothesized that the consumption of contaminated fish with high levels of Hg may have coincided with an increased intake of selenium (Se), an essential element found in fish. Se, in this case, may have moderated the toxic effects of Hg [71, 72] and thus resulted in a decreased OR for low-birth-weight infants among those with higher Hg exposure. Other studies have also reported that Hg levels have a positive association with socioeconomic status among females of childbearing age [73, 74]. Specifically, pregnant participants with low income and low education consumed more fish per week compared to pregnant participants with higher income and higher educational status, but they were found to have lower blood Hg levels [74]. This is because the type of fish consumed by different demographic groups can also affect Hg exposure levels [74]. Therefore, future studies should include Se, the types of fish consumed, and if possible, the whole diet when assessing the potential toxic effects of Hg.
Among the MIREC study participants, we observed that maternal hardships can modify the strength of the relationships (i.e., steeper negative slope and lower birth weight) as well as the direction of the relationships (i.e., positive slope and higher birth weight) between biomarkers of exposures to OCs and metals during pregnancy and infant birth weight. To our knowledge, this study is the largest single cohort study conducted on the modifying effect of maternal hardship on the relationship between gestational environmental chemical concentrations and birth weight. With the benefits of a large sample size and the use of elastic net regularization technique, we had higher statistical power to detect associations and identify key variables. However, our findings should be interpreted with caution due to the following limitations. First, we used multivariable linear regression to assess the independent effects of each maternal hardship and environmental chemical on birth weight and did not account for non-linear relationships. Second, we performed post-selection inference by first using elastic net for variable selection, then using OLS to derive statistical inference from the best-fitting model. A limitation of variable selection methods such as elastic net is their instability where any change in observations may change the model selected [75]. As OLS assumes that we have obtained a perfect model, which may not be the case, the resulting 95% CIs may be inaccurate and too narrow. Accounting for the effects of variable selection is a challenging task and work is underway to develop tools for selective inference so that we can properly assess the strength of the relationships [75]. Third, we used Lubin’s single imputation approach for measurements below the LOD and the standard errors may be biased. However, as the majority of our environmental chemicals were detected in over 90% of the participants, the impact of bias, if any exists, should be minimal [19]. Fourth, maternal hardships were based on a positive-negative dichotomy (i.e., yes-no hardship), which produced easy to interpret findings but essentially weighted all hardships equally. This assumption may have overlooked important and differential aspects of the potential effect of each hardship on the outcome. Furthermore, low maternal hardship proportions of 2 to 25% and missing maternal hardship values may bias the effect estimates and reduce our study power to detect a statistically significant association. Regardless, our use of dichotomized hardships served as an important starting point for the conceptualization of maternal hardships and highlighted the importance for future research to include parental hardships when examining the effects of environmental chemical exposure. Additionally, our use of objective maternal hardships may not fully represent a woman’s experience and well-being [76]. Obtaining qualitative data and additional quantitative data such as food security, stressful life events, and resilience will enhance our understanding of maternal hardships and enrich our current findings. Fifth, we examined only OCs and metals and did not consider co-exposures among the environmental chemicals or hardships. Therefore, we may have missed other environmental chemicals-birth weight associations modified by maternal hardships or any potential 2-way interactions between environmental chemicals or between hardships. Lastly, MIREC is not a representative sample of all Canadian females because it is not population-based [11]. As a result, the generalizability of our results may be limited.