- Open Access
- Open Peer Review
Urinary trace metals, maternal circulating angiogenic biomarkers, and preeclampsia: a single-contaminant and mixture-based approach
Environmental Health volume 18, Article number: 63 (2019)
Exposures to toxic metals and deficiencies in essential metals disrupt placentation and may contribute to preeclampsia. However, effects of exposure to combinations of metals remain unknown.
We investigated the relationship between urinary trace metals, circulating angiogenic biomarkers, and preeclampsia using the LIFECODES birth cohort.
Urine samples collected during pregnancy were analyzed for 17 trace metals and plasma samples were analyzed for soluble fms-like tyrosine-1 (sFlt-1) and placental growth factor (PlGF). Cox proportional hazard models were used to estimate the hazard ratios (HR) of preeclampsia associated with urinary trace metals. Linear regression models were used to estimate the relationship between urinary trace metals and angiogenic biomarkers. Principal components analysis (PCA) was used to identify groups of metals and interactions between principal components (PCs) loaded by toxic and essential metals were examined.
In single-contaminant models, several toxic and essential metals were associated with lower PlGF and higher sFlt-1/PlGF ratio. Detection of urinary chromium was associated with preeclampsia: HR (95% Confidence Interval [CI]) = 3.48 (1.02, 11.8) and an IQR-increase in urinary selenium was associated with reduced risk of preeclampsia (HR: 0.28, 95% CI: 0.08, 0.94). Using PCA, 3 PCs were identified, characterized by essential metals (PC1), toxic metals (PC2), and seafood-associated metals (PC3). PC1 and PC2 were associated with lower PlGF levels, but not preeclampsia risk in the overall cohort.
Trace urinary metals may be associated with adverse profiles of angiogenic biomarkers and preeclampsia.
Preeclampsia is a reproductive disorder that is diagnosed by de novo hypertension and proteinuria after 20 weeks gestation . In the United States, the incidence of preeclampsia has increased over the past several decades and is estimated to affect approximately 5% of all pregnancies [2, 3]. Importantly, the underlying etiology of preeclampsia remains unknown. The origins of preeclampsia are hypothesized to lie within the placenta. Preeclamptic pregnancies are characterized by a placenta with shallow invasion into the maternal decidua and incomplete remodeling of maternal spiral arteries . As a consequence, the preeclamptic placenta is poorly perfused and hypoxic, giving rise to the systemic maternal endothelial dysfunction that characterizes the disorder [3, 4].
Given that preeclampsia is characterized by impaired spiral artery remodeling, angiogenic biomarkers may serve as predictors of preeclampsia risk [5, 6]. These biomarkers include soluble fms-like tyrosine kinase-1 (sFlt-1) and placental growth factor (PlGF). Both sFlt-1 and PlGF are members of the vascular endothelial growth factor (VEGF) family and play important roles in adapting maternal spiral arteries during placentation. More specifically, higher levels of maternal circulating sFlt-1 indicate anti-angiogenic activity, while lower levels of circulating PlGF reflect reduced placental vascularization . Both biomarkers have been shown to predict preeclampsia . In addition, the sFlt-1/PlGF ratio (i.e., a ratio of anti-angiogenic to angiogenic activity) may more accurately predict preeclampsia than either biomarker alone .
Prenatal environmental metals exposure may contribute to preeclampsia risk. Previous studies have reported associations between toxic metals, such as cadmium (Cd) and lead (Pb), and an increased risk of preeclampsia [10,11,12,13,14,15]. Alternatively, higher exposure to essential metals, such as selenium (Se) and zinc (Zn), within a recommended range, have been routinely associated with a decreased risk of preeclampsia [16, 17]. Possible mechanisms for these associations include that toxic metals may impair trophoblast invasiveness [18,19,20,21], generate placental oxidative stress , or lead to maternal immunologic abnormalities , whereas essential metals have antioxidant properties that may reduce the effect of toxic exposures and promote normal placentation . In addition, there are known toxicokinetic and toxicodynamic interactions between toxic and essential metals [25,26,27]. At an epidemiologic level, this has been demonstrated with respect to both placental functioning and preeclampsia, suggesting that essential metals may mitigate the relationship between toxic metals and preeclampsia [10, 28, 29].
While there is evidence for the role of metals in the development of preeclampsia, many of these studies come from populations with exceptionally high levels of toxic metals exposure [11, 12, 14]. Evidence is more limited in populations with lower exposures to toxic metals. In addition, while studies have examined relationships between individual trace metals and preeclampsia, very few studies have investigated metals in combination with one another . In the present study, we performed an exploratory analysis of the relationship between 17 individual trace metals, maternal circulating levels of angiogenic biomarkers, and preeclampsia in the LIFECODES birth cohort. In addition, we investigated the impact of metals mixtures using principal components analysis (PCA) to examine the relationship between groups of urinary metals, preeclampsia, and levels of circulating angiogenic biomarkers. We selected PCA, rather than other statistical methods for analyzing chemical mixtures, as a mixtures-based approach as a method of dimension reduction and to examine groups of metals that may share important exposure sources or toxicities. Lastly, potential interactions between principal components characterized by toxic and essential metals were also investigated.
LIFECODES is an ongoing prospective birth cohort that began recruitment in 2006 at Brigham and Women’s Hospital (BWH) in Boston, MA. The study has few exclusion criteria. Specifically, women are eligible for enrollment into the LIFECODES cohort if they are at least 18 years of age, are seeking prenatal care before 15 weeks gestation, and intend on delivery at BWH . At the first visit (median, 10 weeks gestation), women completed detailed questionnaires of demographic information and medical history and provided informed consent. Gestational age was defined according to the guidelines of the American College of Obstetricians and Gynecologists (ACOG), with the last menstrual period verified by first trimester ultrasound scanning . Additional questionnaire data and urine samples are collected at three subsequent study visits (median, 18, 26, and 35 weeks gestation) and stored at − 80 °C until analysis . All research was approved by the Institutional Review Board at BWH and was deemed exempt by the University of Michigan and the National Institute of Environmental Health Sciences (NIEHS).
The present analysis used subjects from a case-control study of preterm birth recruited from 2006 to 2008, which is nested in the parent LIFECODES cohort (Additional file 1: Figure S1). This case-control study included almost all cases of singleton preterm birth during this period and unmatched singleton controls were included in a 1:3 ratio . The primary purpose of this case-control study was to examine the relationship between prenatal environmental exposures and preterm birth . Participants in the case-control study were included in this analysis if they had urine samples from the 3rd study visit available for metals analysis (n = 390). Importantly, demographic characteristics in the subset of women with urine samples at the 3rd study visit is highly similar to both the larger case-control study of preterm birth  and the overall LIFECODES cohort . To account for the disproportionate number of preterm births in our analysis, we applied inverse probability weighting based on preterm birth case status to all statistical analyses . Thus, results are generalizable to what would have been observed in the parent LIFECODES birth cohort.
After delivery, preeclampsia diagnosis and diagnosis date were abstracted from medical records. Preeclampsia was defined according to ACOG guidelines between 2006 and 2008: elevated maternal blood pressure (> 140 mmHg systolic and/or > 90 mmHg diastolic) and proteinuria (> 300 mg/24 h or a protein/creatinine ratio > 0.20) after 20 weeks gestation . All cases of preeclampsia were verified by two maternal-fetal medicine specialists and diagnosis date was used to determine the gestational age of preeclampsia onset. In the case of conflict, a third specialist reviewed the case. Among the 390 participants, 28 women developed preeclampsia. Among these 28 preeclampsia cases, one woman developed preeclampsia postpartum. For the present analysis, we excluded non-preeclamptics with gestational hypertension (n = 7) to provide a cleaner comparison of women who developed preeclampsia against women who did not develop any form of pregnancy-induced hypertension, resulting in a final sample size of 383.
Urinary trace metals analysis
Urine samples from the 3rd study visit were sent to NSF International (Ann Arbor, MI, USA) for trace metals analysis, in collaboration with the Children’s Health Exposure Analysis Research (CHEAR) Program. Trace metals were analyzed using inductively coupled plasma-mass spectrometry (ICP-MS), with methods described in further detail elsewhere . The metals analyzed included arsenic (As), barium (Ba), beryllium (Be), Cd, copper (Cu), chromium (Cr), mercury (Hg), manganese (Mn), molybdenum (Mo), nickel (Ni), Pb, Se, tin (Sn), thallium (Tl), uranium (U), tungsten (W), and Zn.
When metal concentrations were measured below the limit of detection (LOD) by ICP-MS, the machine-read values were used in analysis . If the reported value was below zero or blank, the value was replaced with LOD/√2 [36, 37]. Metals with very low detection rates (< 30%) were represented as detect versus non-detect in all subsequent analyses.
Plasma biomarkers of angiogenesis
Measurement of circulating maternal sFlt-1 and PlGF were measured in plasma samples collected at the 3rd study visit [5, 38]. Both sFlt-1 and PlGF were measured using ARCHITECT immunoassays (Abbot Laboratory, Abbott Park, IL, USA). Specifically, unbound PlGF concentrations from 1 to 1500 pg/mL were measured and total (unbound and bound) sFlt-1 concentrations from 0.10–150 ng/mL were measured. The ratio of sFlt-1 to PlGF was also calculated for analysis . The combined intra- and interassay coefficients of variation were < 7% for both PlGF and sFlt-1 .
First, differences in demographic and clinical characteristics based on preeclampsia status were examined using t-tests, chi-square tests, or Fisher’s exact tests, where appropriate. Second, we analyzed trace metal distributions by calculating median and interquartile ranges (IQR) by preeclampsia status.
Third, the relationship between urinary trace metals and the hazard ratio (HR) of preeclampsia were estimated using Cox proportional hazard models. Both unadjusted and adjusted Cox models were used to estimate the relationship between an IQR increase in urinary trace metals and the HR (95% Confidence Interval [CI]) of preeclampsia. The proportional hazards assumption was assessed by creating time-dependent variables for each predictor in the model. Tests of these time-dependent variables indicated that the proportional hazards assumption was satisfied. Trace urinary metals were log-transformed in order to improve model fit (i.e. lower AIC). Crude models were created with adjustment for gestational age at sample collection and, for metals that were modeled continuously, urinary specific gravity (SG). Potential confounders for adjusted analyses were identified using a directed acyclic graph (DAG) and included age (years), body mass index (BMI; kg/m2), maternal race (White/African American/Other), maternal education (< high school/technical college/ junior or some college/ > college), maternal health insurance status (private or HMO or self-pay/public), maternal tobacco use during pregnancy (yes/no), alcohol use during pregnancy (yes/no), parity (nulliparous/parous), previous preeclampsia diagnosis (yes/no), gestational diabetes (yes/no), use of assisted reproductive technology (ART; yes/no), chronic hypertension (yes/no), self-reported use of multivitamins (yes/no), calcium supplements (yes/no) and iron supplements (yes/no) during pregnancy, and infant sex (male/female) (Additional file 1: Figure S2). Factors were retained in adjusted analyses if their inclusion influenced effect estimates by > 10%. As a secondary analysis, we recreated the same statistical models excluding BMI as a covariate because of its potential role as a collider .
Fourth, the relationships between urinary trace metals and maternal angiogenic biomarkers were assessed using linear regression models. Trace metal species were log-transformed to improve model fit (i.e. lower AIC) and maternal angiogenic biomarkers were log-transformed to approximate a normal distribution. To ensure comparability between analyses, crude and adjusted models were constructed using the same confounders as the Cox models. Given that both the exposure and the outcome variables were log-transformed, beta estimates and standard errors from these regression models were converted to percent change (95% CI) associated with an IQR increase in trace metals concentration for interpretability.
Lastly, in addition to using single-contaminant models, urinary trace metals were investigated as a mixture using PCA. PCA was applied to the set of 13 urinary trace metals with > 30% detection to produce principal components (PCs). Urine metals measures were log-transformed and adjusted for SG prior to performing PCA to ensure normality and take the effect of hydration on the correlation of urinary metals into account. SG-correction was used to account for urine dilution using the following formula: MetalSG = Metal[(1.015–1)/(SG-1)], where MetalSG is the specific gravity-corrected metal concentration, Metal is the raw urinary trace metal concentration, 1.015 is the median specific gravity of all participants, and SG is the specific gravity level for the sample . A combination of the “Eigenvalue-one” method and the Scree test was used to determine a meaningful set of components. Additionally, PCs were only retained if at least three trace metals loaded onto them with a factor loading > 0.40. Varimax rotation was used to maximize the variance in the factor pattern matrix . While 13 urinary metals were initially included in PCA, Ba and Mo loaded onto multiple components, indicating that they were complex items . Therefore, PCA was re-created excluding these metals. These PCA scores were then fit to  unadjusted and adjusted Cox proportional hazard models to estimate the relationship between a 1-unit change in PC score and the HR of preeclampsia, and  unadjusted and adjusted linear regression models to estimate the relationship between a 1-unit change in PC score and the percent change in circulating angiogenic biomarkers. Given that PCA revealed components loaded by toxic and essential metals, we examined interactions between groups of toxic and essential metals using nested interaction and main effects models. Significant interactions were identified using a likelihood ratio test (LRT). Interactions were considered significant if the LRT p-value < 0.10. For these models, the PC loaded by essential metals was dichotomized at the median (< median vs. > median) in order to report the association between toxic metals and preeclampsia among individuals with low vs. high levels of essential metals.
All analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, USA). Unless otherwise stated, significance was defined as p-value < 0.05.
Subject demographics and urinary metals
Study sample demographics are summarized in Table 1. Briefly, the weighted study population had a median age of 32.7 years, a pre-pregnancy BMI of 23.9 kg/m2, was predominantly White (60.3%), and well-educated (86.5% postsecondary education). Women who developed preeclampsia were more likely to have a higher pre-pregnancy BMI (29.3 vs 23.8 kg/m2; p < 0.01), have been diagnosed with preeclampsia in a previous pregnancy (8.3% vs. 2.3%; p = 0.03), have used ART to become pregnant (24.3% vs. 8.3%; p = 0.02), have chronic hypertension (12.5% vs. 2.3%; p < 0.01), and be carrying a female fetus (77.2% vs. 42.8% male; p < 0.01). Women who developed preeclampsia were also more likely to report using calcium supplements during pregnancy (27.1% vs. 13.3%; p = 0.07) compared to women who did not. With respect to angiogenic biomarkers, circulating levels of PlGF were lower (275 vs. 456 pg/mL; p < 0.01) and the sFlt-1/PlGF ratio was higher (23.2 vs. 13.3; p = 0.04) among women who developed preeclampsia compared to those that did not. Interestingly, although previous analyses in the LIFECODES cohort have demonstrated differences in the trajectory of sFlt-1 levels between preeclamptics and non-preeclamptics throughout pregnancy , circulating sFlt-1 levels at the 3rd study visit did not differ by preeclampsia status.
Medians (IQRs) of SG-adjusted trace metals in 3rd study visit urine samples by preeclampsia status are shown in Table 2. Most trace metals had high detection rates. However, Be, Cr, U, and W were detected in < 30% of the samples. These metals were treated as detect vs. non-detect in all subsequent analyses. Overall, the distribution of metals in the urine was similar in women who developed preeclampsia and women who did not.
Urinary trace metals and preeclampsia risk
After adjusting for smoking during pregnancy, maternal race, maternal educational attainment, insurance status, infant sex, ART use, self-reported calcium supplementation, and pre-pregnancy BMI, the detection of Cr in urine was positively associated with preeclampsia risk (HR: 3.48, 95% CI: 1.02, 11.8) (Table 3). This association was similar in models excluding BMI (see Additional file 1: Table S1). However, it should be noted that Cr was only detected in 15% (N = 57) of all participants and 22% (N = 7) of the women who developed preeclampsia and this association is based on a limited number of observations. In addition, an IQR-increase in urinary Se was associated with a reduction in the risk of preeclampsia (HR: 0.28, 95% CI: 0.08, 0.94).
Urinary trace metals and circulating maternal angiogenic biomarkers
After adjusting for potential confounders, urinary Cu was associated with higher circulating sFlt-1 levels (% Change: 11.5, 95% CI: 0.18, 24.1) (Table 4). No other statistically significant associations were observed between urinary metals and sFlt-1 levels. Several metals were associated with lower levels of circulating PlGF. Namely, an IQR-increase in urinary Cd (% Change: -6.99, 95% CI: − 13.1, − 0.47), urinary Pb (% Change: -7.20, 95% CI: − 11.8, − 2.33), and detection of urinary Cr (% Change: -24.5, 95% CI: − 38.2, − 7.77) was associated with lower levels of circulating PlGF. In addition, urinary Cu and Se were also associated with lower circulating PlGF levels. Urinary Cu was also associated with a 23.7% (95% CI: 6.44, 43.8) higher sFlt-1/PlGF ratio (Table 4). Results from crude models and models excluding BMI remained highly similar (see Additional file 1: Table S2).
Trace metals mixtures, preeclampsia risk, and Angiogenic biomarkers
Using PCA, we identified three primary PCs, which explained 46.0% of the variance in urinary metals. PC1 was characterized by higher loading of essential metals (Cu, Se, and Zn); PC2 was characterized by higher loading of toxic metals (As, Mn, and Pb); and PC3 was characterized by higher loading of seafood-associated metals (As, Hg, and Sn) [43,44,45]. These factor loading patterns have been previously observed in our examination of urinary trace metals and preterm birth in this study population . Factor loadings and communalities for urinary metals are displayed in (see Additional file 1: Table S3).
With respect to preeclampsia, neither PC1 (Cu, Se, and Zn) (HR: 0.89, 95% CI: 0.35, 2.29), PC2 (Cd, Mn, and Pb) (HR: 1.63, 95% CI: 0.74, 3.61), nor PC3 (As, Hg, and Sn) (HR: 0.75, 95% CI: 0.39, 1.46) were statistically significantly associated with preeclampsia (Table 3). Interactions between PCs loaded by toxic/seafood-associated metals (i.e. PC2 and PC3) and PCs loaded by essential metals (PC1) were further investigated (Fig. 1). While there was no statistically significant interaction between PC1 (Cu, Se, and Zn) and PC2 (Cd, Mn, and Pb) (pLRT = 0.12), the association between PC2 (Cd, Mn, and Pb) and preeclampsia (HR: 3.54, 95% CI: 1.09, 11.5) was significant among individuals with low levels of PC1 (Cu, Se, and Zn). On the other hand, this association was null among individuals with higher levels of PC1(Cu, Se, and Zn).
After adjusting for potential confounders, PC1 (Cu, Se, and Zn) (% Change: -13.5, 95% CI: − 22.6, − 3.41) and PC2 (Cd, Mn, and Pb) (% Change: -11.9, 95% CI: -20.7, -2.12) were associated with lower circulating PlGF levels (Table 4). The associations between PCs and angiogenic biomarkers were similar in both crude and adjusted models with BMI excluded (see Additional file 1: Table S2).
In this study, the relationships between 17 trace urinary metals, preeclampsia risk, and levels of maternal circulating angiogenic biomarkers were investigated using both a single-contaminant and mixtures-based approach. In our single-contaminant approach, we observed an increased risk of preeclampsia associated with the detection of Cr in the urine and a reduced risk of preeclampsia associated with an IQR-increase in urinary Se. Urinary Cd, along with other toxic and essential metals, was also associated with reductions in circulating PlGF levels and urinary Cu was associated with higher sFlt-1 levels and a higher sFlt-1/PlGF ratio, suggesting an association with impaired placentation and preeclampsia risk. In our mixtures approach, we used PCA to further investigate how combinations of urinary metals were associated with the risk of preeclampsia and circulating angiogenic biomarker levels. Individual PCs were not associated with preeclampsia risk, but both PC1 (Cu, Se, and Zn) and PC2 (Cd, Mn, and Pb) were associated with lower levels of circulating PlGF. In addition, while we did not observe a statistically significant interaction between PC1 (Cu, Se, and Zn) and PC2 (Cd, Mn, and Pb), the association between PC2 (Cd, Mn, and Pb) and preeclampsia was significant among individuals with low levels of PC1 (Cu, Se, and Zn).
In single-contaminant models, detection of Cr in the urine was associated with an increased risk of preeclampsia and a reduction in circulating PlGF levels, suggesting that Cr is associated with a reduction in placentation and placental function. While the effect estimates were relatively imprecise, other studies have also noted higher levels of total urinary or hair Cr in preeclampsia cases compared to normotensive controls [11, 15]. Urinary Cr levels in this cohort were low, with most samples being below the LOD, but the primary source of exposure to Cr in the general population is from food, particularly from meat and fish products . Recent toxicologic studies have also demonstrated that hexavalent Cr exposure leads to increased expression of oxidative stress markers and apoptotic signaling in trophoblast cells and mouse placenta [47,48,49], providing further molecular-level evidence for a potential link between Cr exposure and preeclampsia. Despite the consistency of association between urinary Cr, angiogenic biomarkers, and preeclampsia, we did not feel that it was appropriate to apply mediation analysis methodologies to this cohort given the small number of samples with Cr above the LOD, the relatively high degree of imprecision in the estimates, and the lack of clear connection between Cr exposure and preeclampsia in current epidemiologic and toxicologic literature. Nevertheless, we feel that these results warrant further research.
Although we did not find associations between other toxic metals and the risk of preeclampsia, we identified associations between urinary Cd and Pb and lower circulating PlGF levels. These results may provide additional mechanistic evidence for previous studies that have reported associations between prenatal Cd and Pb exposure and preeclampsia [10, 13, 50]. Both Cd and Pb are ubiquitous environmental exposures. Both may be found in contaminated air, soil, or water. Cd levels tend to be higher in certain foods, such as rice and grains, and both Cd and Pb can be found in drinking water systems or housing, particularly in areas with old or failing infrastructure [51,52,53]. Similarly, the PC loaded by toxic metals (PC2: Cd, Pb, and Mn) was not associated with preeclampsia risk, but was associated with lower circulating PlGF levels.
With respect to essential metals, we observed a reduction, albeit imprecise, in the risk of preeclampsia associated with an increase in urinary Se. This observation is consistent with the literature, which suggests that essential metals, including Se, may reduce the risk of preeclampsia [16, 17]. In contrast, urinary Cu and Se were associated with lower circulating PlGF levels. Urinary Cu was also associated with higher sFlt-1 and the sFlt-1/PlGF ratio. Likewise, in our mixtures-based approach, the PC for essential metals (PC1: Cu, Se, and Zn) was associated with lower circulating PlGF levels. These findings may suggest that, instead of having a beneficial effect, levels of essential metals are associated with impaired placental angiogenesis and placentation. This is consistent with previously observed associations between urinary Cu and the risk of preterm birth, specifically placentally-mediated preterm birth, in this study population . In addition, recent meta-analyses have reported consistent observations that higher plasma or serum Cu is associated with preeclampsia risk [54, 55]. While Cu is an integral component of the antioxidant system, modest increases in Cu may lead to the production of reactive oxygen species and oxidative stress, which has been linked to the development of preeclampsia [9, 56]. Alternatively, it has been suggested that higher levels of essential metals, including Cu, may be the result of preeclampsia, where women developing preeclampsia may increase the uptake of essential metals in attempt to buffer the higher levels oxidative stress associated with the disease .
In addition to PC1 (Cu, Se, and Zn) and PC2 (Cd, Mn, and Pb), we identified a third PC that was characterized by loading from seafood-associated metals (PC3: As, Hg, and Sn). Notably, while inorganic As is considered a toxic metal, non-toxic organic forms of As, such as arsenobetaine, are widely found in seafood . In populations with frequent seafood consumption, including those within the United States, much of the total urinary As may be organic arsenicals from such foods [58,59,60]. In addition, although seafood-associated metals (PC3: As, Hg, and Sn) grouped independently of other toxic metals (PC2: Cd, Mn, and Pb), they should not be considered non-toxic. For instance, pregnant women are advised against consuming certain types of seafood on the basis that it may contain high levels of methylmercury . Despite this, we did not identify associations between PC3 (As, Hg, and Sn) and either preeclampsia or angiogenic biomarkers. Disentangling the harmful effects of these trace metals from the beneficial effects of seafood consumption (i.e. fatty acids) is complex and beyond the scope of this analysis.
Given known toxicologic interactions between toxic and essential metals, we investigated statistical interactions between toxic and essential metals with respect to both preeclampsia risk and levels of circulating angiogenic biomarkers. While there was no significant interaction between PC1 (Cu, Se, and Zn) and PC2 (Cd, Mn, and Pb), the association between PC2 (Cd, Mn, and Pb) and preeclampsia was significant among individuals with lower levels of PC1 (Cu, Se, and Zn). Previously, significant interactions between placental Cd and Se levels were reported in a case-control study of preeclampsia . At a molecular level, these findings are supported by evidence of toxic and essential metals interactions within the placenta. For example, antagonistic interactions between maternal toenail Cd and Se levels have been noted with respect to placental apoptotic gene expression . Taken together, these results may suggest that individuals with low levels of essential metals may be more susceptible to the effects of toxic exposures on placental development and/or functioning.
This study has several strengths. First, the LIFECODES birth cohort has a prospective study design. Urine samples were collected at the 3rd study visit and no patients had developed preeclampsia at the time of sample collection, reducing the likelihood that these results were influenced by reverse causality. Nevertheless, it is possible that preeclampsia was subclinical at the time of sample collection (median, 26 weeks gestation) because the origins of preeclampsia may lie in the early stages of pregnancy. Second, this study is among the first to employ a mixtures approach in the context of preeclampsia. This is of critical importance given that populations are exposed to complex mixtures of toxicants, rather than single-contaminants, and represents a major public health research objective .
Despite these strengths, this study also has several limitations. First, this study had a small number of preeclampsia cases. Thus, it was not possible to examine the relationship between urinary metals and individual subtypes of preeclampsia (i.e., early- vs. late-onset). There is evidence that these subtypes may have different underlying etiologies . It should also be noted that the sampling design of the case-control study that this analysis is based in may have resulted in an oversampling of early-onset preeclampsia. However, this oversampling should be partially accounted for by the inverse weighting approach used to account for selection into the study on preterm birth status. In addition, the limited sample size precluded an investigation of different dose-response shapes for trace metals and preeclampsia. Some, such as Mn, may feature non-monotonic dose-response curves that are not well captured by linear modeling. Second, with respect to exposure assessment, biomarkers of exposure were only available as urinary measures and at a single timepoint. Exposure to some metals may be better captured using other matrices, such as blood biomarkers for Se, Mn, Zn or Pb exposure [64,65,66,67]. Moreover, multiple metals measures taken throughout pregnancy would have allowed for a better characterization of metals exposure throughout pregnancy. In particular, metals measures taken earlier during pregnancy may better capture exposure at the time of early preeclampsia development. Third, we used PCA as a mixtures-based approach to examine urinary trace metals and preeclampsia. Like all statistical methods for examining chemical mixtures, PCA has limitations . Unsupervised PCA reduces dimensionality by selecting linear combinations of variables without considering how they relate to the outcome of interest. Thus, while it identifies chemicals with similar exposure patterns, it may not have selected the most relevant groups of trace metals for preeclampsia or angiogenic biomarkers. Moreover, because exposure patterns vary across populations, the PCs identified in this study may not be generalizable to other studies. In addition, by assuming linear combinations of variables within a PC, it assumes no interactions between metals within each cluster. Lastly, the generalizability of this study may be limited given that participants were recruited from BWH, a tertiary care center, which receives higher-risk pregnancies than other hospitals. However, this population does exhibit similar concentrations of urinary trace metals to those observed in the general United States population .
This study investigated 17 urinary metals, individually and in mixtures, in association with circulating maternal angiogenic biomarkers and preeclampsia. This study is among the first to examine the relationship between prenatal trace metals exposure and maternal angiogenic biomarkers, as well as the relationship between metals mixtures and preeclampsia. The results from both our single-contaminant and mixtures-based approach suggests that both toxic and essential urinary metals are associated with adverse maternal profiles of angiogenic biomarkers. In addition, there may be important interactions between toxic and essential metals in the context of preeclampsia.
Availability of data and materials
The datasets generated and analyzed during the current study are not publically available due to the sensitive nature of human biological and environmental exposure data, but are available from the corresponding author upon reasonable request.
American College of Obstetricians and Gynecologists
Assisted reproductive technology
Body mass index
Brigham and Women’s Hospital
Children’s Health Exposure Analysis Research
Directed acyclic graph
Inductively coupled plasma-mass spectrometry
Limit of detection
Likelihood ratio test
National Institue of Environmental Health Sciences
Principal components analysis
Placental growth factor
soluble fms-like tyrosine kinase-1
Vascular endothelial growth factor
ACOG. Hypertension in pregnancy. Report of the American College of Obstetricians and Gynecologists’ task force on hypertension in pregnancy. Obstet Gynecol. 2013;122(5):1122.
Wallis AB, Saftlas AF, Hsia J, Atrash HK. Secular trends in the rates of preeclampsia, eclampsia, and gestational hypertension, United States, 1987-2004. Am J Hypertens. 2008;21(5):521–6 PubMed PMID: 18437143. Epub 2008/04/26. eng.
Steegers EA, von Dadelszen P, Duvekot JJ, Pijnenborg R. Pre-eclampsia. Lancet. 2010;376(9741):631–44 PubMed PMID: 20598363. Epub 2010/07/06. eng.
Fisher SJ. Why is placentation abnormal in preeclampsia? Am J Obstet Gynecol. 2015;213(4 Suppl):S115–22 PubMed PMID: 26428489. Pubmed Central PMCID: PMC4592742. Epub 2015/10/03. eng.
McElrath TF, Lim KH, Pare E, Rich-Edwards J, Pucci D, Troisi R, et al. Longitudinal evaluation of predictive value for preeclampsia of circulating angiogenic factors through pregnancy. Am J Obstet Gynecol. 2012;207(5):407 e1–7 PubMed PMID: 22981320. Epub 2012/09/18. eng.
Levine RJ, Maynard SE, Qian C, Lim KH, England LJ, Yu KF, et al. Circulating angiogenic factors and the risk of preeclampsia. N Engl J Med. 2004;350(7):672–83 PubMed PMID: 14764923. Epub 2004/02/07. eng.
Coolman M, Timmermans S, de Groot CJ, Russcher H, Lindemans J, Hofman A, et al. Angiogenic and fibrinolytic factors in blood during the first half of pregnancy and adverse pregnancy outcomes. Obstet Gynecol. 2012;119(6):1190–200 PubMed PMID: 22617584. Epub 2012/05/24. eng.
Allen RE, Rogozinska E, Cleverly K, Aquilina J, Thangaratinam S. Abnormal blood biomarkers in early pregnancy are associated with preeclampsia: a meta-analysis. Eur J Obstet Gynecol Reprod Biol. 2014;2014(182):194–201.
Verlohren S, Herraiz I, Lapaire O, Schlembach D, Moertl M, Zeisler H, et al. The sFlt-1/PlGF ratio in different types of hypertensive pregnancy disorders and its prognostic potential in preeclamptic patients. Am J Obstet Gynecol. 2012;206(1):58 e1–8 PubMed PMID: 22000672. Epub 2011/10/18. eng.
Laine JE, Ray P, Bodnar W, Cable PH, Boggess K, Offenbacher S, et al. Placental cadmium levels are associated with increased preeclampsia risk. PloS One. 2015;10(9):e0139341 PubMed PMID: 26422011. Pubmed Central PMCID: PMC4589375. Epub 2015/10/01. eng.
Elongi Moyene JP, Scheers H, Tandu-Umba B, Haufroid V, Buassa-Bu-Tsumbu B, Verdonck F, et al. Preeclampsia and toxic metals: a case-control study in Kinshasa, DR Congo. Environ health. 2016;15:48 PubMed PMID: 27044488. Pubmed Central PMCID: PMC4820935. Epub 2016/04/06. eng.
Jameil NA. Maternal serum lead levels and risk of preeclampsia in pregnant women: a cohort study in a maternity hospital, Riyadh, Saudi Arabia. Int J Clin Exp Pathol. 2014;7(6):3182–9 PubMed PMID: 25031738. Pubmed Central PMCID: PMC4097270. Epub 2014/07/18. eng.
Wang F, Fan F, Wang L, Ye W, Zhang Q, Xie S. Maternal cadmium levels during pregnancy and the relationship with preeclampsia and fetal biometric parameters. Biol Trace Elem Res. 2018;186(2):322-9. PubMed PMID: 29651732. Epub 2018/04/14. eng.
Vigeh M, Yokoyama K, Ramezanzadeh F, Dahaghin M, Sakai T, Morita Y, et al. Lead and other trace metals in preeclampsia: a case-control study in Tehran, Iran. Environ Res. 2006;100(2):268–75 PubMed PMID: 16029873. Epub 2005/07/21. eng.
Maduray K, Moodley J, Soobramoney C, Moodley R, Naicker T. Elemental analysis of serum and hair from pre-eclamptic South African women. J Trace Elem Med Biol. 2017;43:180–6 PubMed PMID: 28325649. Epub 2017/03/23. eng.
Zhu Q, Zhang L, Chen X, Zhou J, Liu J, Chen J. Association between zinc level and the risk of preeclampsia: a meta-analysis. Arch Gynecol Obstet. 2016;293(2):377–82 PubMed PMID: 26386964. Epub 2015/09/21. eng.
Xu M, Guo D, Gu H, Zhang L, Lv S. Selenium and preeclampsia: a systematic review and meta-analysis. Biol Trace Elem Res. 2016;171(2):283–92 PubMed PMID: 26516080. Epub 2015/10/31. eng.
Brooks SA, Martin E, Smeester L, Grace MR, Boggess K, Fry RC. miRNAs as common regulators of the transforming growth factor (TGF)-beta pathway in the preeclamptic placenta and cadmium-treated trophoblasts: Links between the environment, the epigenome and preeclampsia. Food Chem Toxicol. 2016;98(Pt A):50–7 PubMed PMID: 27375191. Pubmed Central PMCID: PMC5156314. Epub 2016/10/30. eng.
Brooks SA, Fry RC. Cadmium inhibits placental trophoblast cell migration via miRNA regulation of the transforming growth factor beta (TGF-beta) pathway. Food Chem Toxicol. 2017;109(Pt 1):721–6 PubMed PMID: 28774740. Pubmed Central PMCID: PMC5656529. Epub 2017/08/05. eng.
Alvarez MM, Chakraborty C. Cadmium inhibits motility factor-dependent migration of human trophoblast cells. Toxicol In Vitro. 2011;25(8):1926–33 PubMed PMID: 21745561. Epub 2011/07/13. eng.
Erboga M, Kanter M. Effect of cadmium on trophoblast cell proliferation and apoptosis in different gestation periods of rat placenta. Biol Trace Elem Res. 2016;169(2):285–93 PubMed PMID: 26170172. Epub 2015/07/15. eng.
Zhang X, Xu Z, Lin F, Wang F, Ye D, Huang Y. Increased oxidative DNA damage in placenta contributes to cadmium-induced Preeclamptic conditions in rat. Biol Trace Elem Res. 2016;170(1):119–27 PubMed PMID: 26194818. Epub 2015/07/22. eng.
Zhang Q, Huang Y, Zhang K, Huang Y, Yan Y, Wang F, et al. Cadmium-induced immune abnormality is a key pathogenic event in human and rat models of preeclampsia. Environ Pollut. 2016;218:770–82 PubMed PMID: 27511439. Epub 2016/08/12. eng.
Ma Y, Shen X, Zhang D. The relationship between Serum Zinc level and Preeclampsia: a meta-analysis. Nutrients. 2015;7(9):7806–20 09/15 07/23/received 09/07/accepted. PubMed PMID: PMC4586561.
Goyer RA. Toxic and essential metal interactions. Annu Rev Nutr. 1997;17(1):37–50.
Brzoska MM, Moniuszko-Jakoniuk J. Interactions between cadmium and zinc in the organism. Food Chem Toxicol. 2001;39(10):967–80 PubMed PMID: 11524135. Epub 2001/08/29. eng.
Zwolak I, Zaporowska H. Selenium interactions and toxicity: a review. Selenium interactions and toxicity. Cell Biol Toxicol. 2012;28(1):31–46 PubMed PMID: 21913064. Epub 2011/09/14. eng.
Al-Saleh I, Al-Rouqi R, Obsum CA, Shinwari N, Mashhour A, Billedo G, et al. Interaction between cadmium (cd), selenium (se) and oxidative stress biomarkers in healthy mothers and its impact on birth anthropometric measures. Int J Hyg Environ Health. 2015;218(1):66–90 PubMed PMID: 25239374. Epub 2014/09/23. eng.
Everson TM, Kappil M, Hao K, Jackson BP, Punshon T, Karagas MR, et al. Maternal exposure to selenium and cadmium, fetal growth, and placental expression of steroidogenic and apoptotic genes. Environ Res. 2017;158:233–44 PubMed PMID: 28662449. Pubmed Central PMCID: PMC5554457. Epub 2017/07/01. eng.
Cantonwine DE, Meeker JD, Ferguson KK, Mukherjee B, Hauser R, McElrath TF. Urinary concentrations of bisphenol a and phthalate metabolites measured during pregnancy and risk of preeclampsia. Environ Health Perspect. 2016;124(10):1651–5.
ACOG. Committee opinion no 611: method for estimating due date. Obstet Gynecol. 2014;124(4):863–6.
Kim SS, Meeker JD, Carroll R, Zhao S, Mourgas MJ, Richards MJ, et al. Urinary trace metals individually and in mixtures in association with preterm birth. Environ Int. 2018;2018(121):582–90.
Cantonwine DE, Ferguson KK, Mukherjee B, Chen Y-H, Smith NA, Robinson JN, et al. Utilizing longitudinal measures of fetal growth to create a standard method to assess the impacts of maternal disease and environmental exposure. PLoS One. 2016;11(1):e0146532.
Ferguson KK, McElrath TF, Meeker JD. Environmental phthalate exposure and preterm birth. JAMA Pediatr. 2014;168(1):61–7 PubMed PMID: PMC4005250.
Jiang Y, Scott AJ, Wild CJ. Secondary analysis of case-control data. Stat Med. 2006;25(8):1323–39 PubMed PMID: 16220494. Epub 2005/10/13. eng.
Reed LD, Hornung RW. Estimation of average concentration in the presence of nondetectable values. Appl Occup Environ Hyg. 1990;5(1):46–51.
Hornung RW, Reed LD. Estimation of average concentration in the presence of nondetectable values. Appl Occup Environ Hyg. 1990;5(1):46–51.
Ferguson KK, McElrath TF, Cantonwine DE, Mukherjee B, Meeker JD. Phthalate metabolites and bisphenol-A in association with circulating angiogenic biomarkers across pregnancy. Placenta. 2015;36(6):699–703 PubMed PMID: PMC4441857.
Greenland S. Quantifying biases in causal models: classical confounding vs collider-stratification bias. Epidemiology. 2003;14(3):300–6 PubMed PMID: 12859030. Epub 2003/07/16. eng.
Ferguson KK, McElrath TF, Meeker JD. Environmental phthalate exposure and preterm birth. JAMA Pediatr. 2014;168(1):61–7 PubMed PMID: 24247736. Pubmed Central PMCID: PMC4005250. Epub 2013/11/20. eng.
O’Rourke N, Hatcher L. A step-by-step approach to using SAS for factor analysis and structural equation modeling. 2nd ed. Cary, NC: SAS Institute, Inc; 2013.
Honigberg MC, Cantonwine DE, Thomas AM, Lim KH, Parry SI, McElrath TF. Analysis of changes in maternal circulating angiogenic factors throughout pregnancy for the prediction of preeclampsia. J Perinatol. 2016;36(3):172–7 PubMed PMID: 26583938. Epub 2015/11/20. eng.
ATSDR. Toxicological profile for mercury. 1999.
ATSDR. Toxicological profile for tin and tin compounds. 2005.
ATSDR. Toxicological profile for arsenic. 2007.
ATSDR. Toxicological Profile for Chromium. 2012.
Banu SK, Stanley JA, Taylor RJ, Sivakumar KK, Arosh JA, Zeng L, et al. Sexually Dimorphic Impact of Chromium Accumulation on Human Placental Oxidative Stress and Apoptosis. Toxicol Sci. 2018;161(2):375–87 PubMed PMID: 29069462. Pubmed Central PMCID: PMC5837378. Epub 2017/10/27. eng.
Banu SK, Stanley JA, Sivakumar KK, Arosh JA, Taylor RJ, Burghardt RC. Chromium VI - Induced developmental toxicity of placenta is mediated through spatiotemporal dysregulation of cell survival and apoptotic proteins. Reprod Toxicol. 2017;68:171–90 PubMed PMID: 27443218. Pubmed Central PMCID: PMC5830189. Epub 2016/07/23. eng.
Lee H, Chun JH, Moon DH, Lee CU, Kang SG, Son BC, et al. Effects of chromium (VI) exposure on the placental function and reproduction in rats. J Prev Med Public Health. 2004;37(2):157–65 PubMed PMID: 25178447. Epub 2004/05/01. eng.
Kennedy DA, Woodland C, Koren G. Lead exposure, gestational hypertension and pre-eclampsia: a systematic review of cause and effect. J Obstet Gynaecol. 2012;32(6):512–7 PubMed PMID: 22779950. Epub 2012/07/12. eng.
Shiue I, Bramley G. Environmental chemicals mediated the effect of old housing on adult health problems: US NHANES, 2009-2010. Environ Sci Pollut Res Int. 2015;22(2):1299–308 PubMed PMID: 25138559. Epub 2014/08/21. eng.
Tyrrell J, Melzer D, Henley W, Galloway TS, Osborne NJ. Associations between socioeconomic status and environmental toxicant concentrations in adults in the USA: NHANES 2001-2010. Environ Int. 2013;59:328–35 PubMed PMID: 23892225. Epub 2013/07/31. eng.
Kim K, Melough MM, Vance TM, Noh H, Koo SI, Chun OK. Dietary cadmium intake and sources in the US. Nutrients. 2018;11(1):2. PubMed PMID: 30577418. Pubmed Central PMCID: PMC6356330. Epub 2018/12/24. eng. https://www.mdpi.com/about/announcements/784.
Fan Y, Kang Y, Zhang M. A meta-analysis of copper level and risk of preeclampsia: evidence from 12 publications. Biosci Rep. 2016;36(4):e00370. PubMed PMID: 27407173. Pubmed Central PMCID: PMC4995498. Epub 2016/07/14. eng.
Song X, Li B, Li Z, Wang J, Zhang D. High serum copper level is associated with an increased risk of preeclampsia in Asians: a meta-analysis. Nutr Res. 2017;39:14–24 PubMed PMID: 28385285. Epub 2017/04/08. eng.
Mistry HD, Gill CA, Kurlak LO, Seed PT, Hesketh JE, Meplan C, et al. Association between maternal micronutrient status, oxidative stress, and common genetic variants in antioxidant enzymes at 15 weeks gestation in nulliparous women who subsequently develop preeclampsia. Free Radical Biol Med. 2015;78:147–55 PubMed PMID: 25463281. Pubmed Central PMCID: PMC4291148. Epub 2014/12/03. eng.
Cullen WR, Reimer KJ. Arsenic speciation in the environment. Chem Rev. 1989;89(4):713–64.
Aylward LL, Ramasamy S, Hays SM, Schoeny R, Kirman CR. Evaluation of urinary speciated arsenic in NHANES: issues in interpretation in the context of potential inorganic arsenic exposure. Regul Toxicol Pharmacol. 2014;69(1):49–54 PubMed PMID: 24582650.
Jones MR, Tellez-Plaza M, Vaidya D, Grau M, Francesconi KA, Goessler W, et al. Estimation of inorganic arsenic exposure in populations with frequent seafood intake: evidence from MESA and NHANES. Am J Epidemiol. 2016;184(8):590–602 PubMed PMID: 27702745. Pubmed Central PMCID: 5065621.
Navas-Acien A, Francesconi KA, Silbergeld EK, Guallar E. Seafood intake and urine concentrations of total arsenic, dimethylarsinate and arsenobetaine in the US population. Environ Res. 2011;111(1):110–8 PubMed PMID: 21093857. Epub 11/19. eng.
Regist F. Advice about eating fish, from the environmental protection agency and the food and drug administration; revised fish advice; availability; 2017. p. 6571–4.
Carlin DJ, Rider CV, Woychik R, Birnbaum LS. Unraveling the health effects of environmental mixtures: an NIEHS priority. Environ Health Perspect. 2013;121(1):A6–8 PubMed PMID: 23409283. Pubmed Central PMCID: PMC3553446. Epub 2013/02/15. eng.
von Dadelszen P, Magee LA, Roberts JM. Subclassification of preeclampsia. Hypertens Pregnancy. 2003;22(2):143–8 PubMed PMID: 12908998. Epub 2003/08/12. eng.
Barbosa F, Tanus-Santos JE, Gerlach RF, Parsons PJ. A critical review of biomarkers used for monitoring human exposure to lead: advantages, limitations, and future needs. Environ Health Perspect. 2005;113(12):1669–74 08/11 01/11/received 08/10/accepted. PubMed PMID: PMC1314903.
Lowe NM, Fekete K, Decsi T. Methods of assessment of zinc status in humans: a systematic review. Am J Clin Nutr. 2009;89(6):2040S–51S.
Ashton K, Hooper L, Harvey LJ, Hurst R, Casgrain A, Fairweather-Tait SJ. Methods of assessment of selenium status in humans: a systematic review. Am J Clin Nutr. 2009;89(6):2025S–39S PubMed PMID: 19420095. Epub 2009/05/08. eng.
O’Neal SL, Zheng W. Manganese toxicity upon overexposure: a decade in review. Curr Environ Health Rep. 2015;2(3):315–28 PubMed PMID: 26231508. Pubmed Central PMCID: PMC4545267. Epub 2015/08/02. eng.
Braun JM, Gennings C, Hauser R, Webster TF. What Can Epidemiological Studies Tell Us about the Impact of Chemical Mixtures on Human Health? Environ Health Perspect. 2016;124(1):A6–9 PubMed PMID: 26720830. Pubmed Central PMCID: PMC4710611. Epub 2016/01/01. eng.
Ethics approval and consent to participate
All participants provided written informed consent and research was approved by the Institutional Review Board at Brigham and Women’s Hospital. This research was deemed exempt by the Institutional Review Boards at the University of Michigan and the National Institute of Environmental Health Sciences.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Flow-through diagram displaying study design. Figure S2. Simplified directed acyclic graph (DAG) used to identify potential confounders and/or colliders. Table S1. Unadjusted and adjusted relationship between urinary metals and the HR (95% CI) of preeclampsia. Table S2. Unadjusted and adjusted relationship between urinary trace metals and the percent change (95% CI) in circulating maternal angiogenic biomarkers. Table S3. Standardized and rotated loading factors and communalities for each variable within each principal component. (DOCX 208 kb)