Perinatal urinary benzophenone-3 concentrations and glucose levels among women from a fertility clinic
Environmental Health volume 19, Article number: 45 (2020)
Subfertile women have higher risk of glucose intolerance during pregnancy. Studies suggest associations between several endocrine disrupting chemicals (EDCs) and pregnancy glucose levels. However, the association between benzophenone-3 (BP-3), an EDC widely found in sunscreen, and pregnancy glucose levels remains unclear. We aimed to assess the association between perinatal exposures to BP-3 and pregnancy glucose levels in subfertile women.
We evaluated 217 women from a prospective cohort based at a fertility clinic who had urinary BP-3 concentrations measured during 3-month preconception, first and/or second trimesters, and blood glucose measured at glucose load tests (GLTs) during late pregnancy. Multivariable linear and logistic regression models were used to assess associations between time-specific BP-3 in quartiles (Q1 – Q4) and mean glucose levels, as well as odds of abnormal GLT (glucose level ≥ 140 mg/dL), adjusting for potential confounders. Effect modification was assessed by age, season, BMI, infertility diagnosis, sex of fetus (es) and physical activity.
Women with higher first trimester BP-3 concentrations had lower mean glucose levels [mean glucose (95% CI) for Q4 vs Q1 = 103.4 (95.0, 112.5) vs. 114.6 (105.8, 124.2) mg/dL]. Women with higher second trimester BP-3 concentrations had lower odds of abnormal GLT [OR (95% CI) for Q3 vs. Q1 = 0.12 (0.01, 0.94)]. The associations between BP-3 and glucose levels were modified by several factors: women with female-factor infertility, urine collected during summer, older age, lower BMI, or carried female fetus (es) had the strongest inverse associations between BP-3 and glucose levels, while no associations were observed in the remaining subgroups.
Time-specific inverse associations between BP-3 and pregnancy glucose levels existed in subfertile women, and especially among certain subgroups of this high-risk-population.
Gestational diabetes mellitus (GDM) has increased over the past decades, currently affecting 7.6% of pregnant women in the United States . Certain subgroups of the population have even greater risk of GDM, including those who conceived through assisted reproductive techniques (ARTs) compared to women who conceived naturally . Subfertile women have almost double the risk of GDM conditions, such as polycystic ovarian syndrome (PCOS) or metabolic responses to ART (e.g. response to induced ovulation, underlying exacerbated metabolic factors during ART) that may lead to glucose dysregulation in pregnancy . With infertility affecting 10–20% of U.S. couples , the increased risk of GDM remains a concern, as it is associated with adverse outcomes including premature delivery, preeclampsia, cesarean delivery, or neonatal hypoglycemia . Interestingly, women with elevated glucose levels, even below the limit of clinical GDM, are at higher risk of these complications .
Aside from traditional lifestyle factors , there has been growing evidence suggesting that some environmental factors like endocrine disrupting chemicals (EDCs) may also contribute to elevating glucose levels during pregnancy . For example, persistent organic pollutants (POPs) were found to increase GDM risk in average-risk populations . In addition, we previously found associations between pregnancy glucose levels and certain non-persistent chemicals, including bisphenol A , certain parabens , and certain phthalates , in higher-risk (i.e. subfertile) populations from the same study cohort. However, benzophenone-3 (BP-3)—a highly ubiquitous EDC—has not been assessed in relation to glucose levels, particularly in the insulin-resistant period during pregnancy .
BP-3 is a derivative of benzophenone-type ultraviolet (UV) blockers. It was detected in over 96% of the U.S. general population between 2003 and 2012  and 100% of a diverse cohort of U.S. pregnant women between 2009 and 2010 . BP-3 is widely used in sunscreen products, as well as other consumer products including shampoos, fragrances, nail polishes, furniture, clothing, carpets, and in plastic products as UV stabilizers . While other EDCs have been found to increase the risk of glucose intolerance through regulating pancreatic β-cell function, inducing insulin resistance, or promoting oxidative stress [17,18,19], the role of BP-3 in glucose regulation remains understudied. A few epidemiological studies found positive associations between maternal benzophenone biomarkers and regulators for insulin metabolization (e.g. IGF-I and its binding protein IGFBP3)  or predictors for diabetes (e.g. markers of systematic inflammation and oxidative stress) [21, 22], yet another case-control study in Saudi Arabia suggested a potential inverse association between BP-3 concentrations and type 2 diabetes . However, the associations between perinatal BP-3 concentrations and pregnancy glucose levels have not been evaluated. It is also questionable how timing of exposure would have any differing effects on glucose levels, as BP-3 is non-persistent , and maternal glucose metabolism changes across pregnancy, with increasing insulin resistance during mid- to late-pregnancy .
Thus, we aimed to assess the association between time-specific (i.e. preconception, first trimester and second trimester) BP-3 urinary concentrations and blood glucose levels in late pregnancy among subfertile women seeking care at a fertility clinic. We evaluated whether season, maternal age, body mass index (BMI), infertility diagnosis, sex of fetus (es) and physical activity might be potential modifiers of the associations. Understanding the role of this widespread but understudied EDC and its relation to pregnancy glucose levels could provide implications for chemical exposures and GDM risk in a high-risk (i.e. subfertile) subgroup of women.
This study was conducted among a subset of women participating in the Environment and Reproductive Health (EARTH) Study, an ongoing prospective cohort assessing environmental impact on reproduction based on couples seeking treatment at the Massachusetts General Hospital (MGH) Fertility Center . In this study, we included women who: (1) provided at least one urine sample at preconception (defined as within 3 months prior to conception of the index pregnancy), first trimester (median: 8 weeks gestation), and/or second trimester (median: 22 weeks gestation); (2) had blood glucose levels measured from the 50-g glucose load test (GLT), as a part of the GDM screening test used universally in this cohort; (3) for those who were pregnant more than once during study participation (n = 15), only their first pregnancy data was included; (4) had a live birth. Women with a prior history of diabetes at baseline (including medical/self-reported diabetes, or reported to be on diabetes medications upon enrollment) were excluded (n = 6). In total, 217 women were included in the study, who provided 833 urine samples (469 samples from 178 women during 3-month preconception, 194 samples from 194 women during the first trimester, and 170 samples from 170 women during the second trimester). The study period ranged from 2009 to 2017. We obtained signed informed consent from all participants. This study was approved by the Partners IRB, Harvard IRB, and the Centers for Disease Control and Prevention (CDC).
Urine collection and BP-3 quantification
Spot urine samples were collected during preconception, first trimester, and/or second trimester in sterile polypropylene cups. Specific gravity (SG)  was measured with a handheld refractometer (National Instrument Company, Inc., Baltimore, MD, USA). Urine samples were divided, frozen (− 20 °C) and stored (− 80 °C) before being transported to the CDC (Atlanta, GA, USA) on dry ice for quantification of BP-3 concentrations. Online solid phase extraction along with high-performance liquid chromatography tandem mass spectrometry was utilized for quantification (detailed description elsewhere ). The limit of detection (LOD) was 0.4 ng/mL or 0.2 ng/mL depending on year of analysis.
Blood glucose levels during pregnancy were assessed using a 1-h non-fasting, 50-g glucose load test (GLT) at late-second/early-third trimester (median: 27 weeks gestation) as part of the two-step method with Carpenter-Coustan criteria for GDM screening in this study population. Data were abstracted from medical records. In this study, blood glucose levels were used as: (1) continuous glucose levels from GLT; (2) dichotomized glucose status, where abnormal GLT was defined as glucose level ≥ 140 mg/dL , in which women would have been referred for additional GDM screening based on the elevated glucose level.
Clinical data and covariates assessment
Sociodemographic and lifestyle factors were collected from questionnaires upon enrollment. Weight (kg) and height (m) were measured by trained study staff and body mass index (BMI) was calculated as weight/height2. Physician diagnosis of PCOS was abstracted from electronic medical records. Infertility was diagnosed by physicians based on the Society for Assisted Reproductive Technology definitions [i.e. female-factor (including diminished ovarian reserve, endometriosis, uterine, ovulatory, tubal, other), male-factor, or unexplained]. The study included pregnancies conceived by in vitro fertilization (IVF), intrauterine insemination (IUI), or natural conception. Number and sex of fetus (es) per birth were abstracted from maternal delivery records. Dates of urine sample collection and GLT were collected to determine seasonality (winter: Dec-Feb, spring: Mar-May, summer: Jun-Aug, fall: Sep-Nov).
Urinary dilution was accounted for by using SG-adjusted BP-3 concentration (ng/mL) using the formula Pc = P [(1.016–1) / (SG – 1)], where Pc is the SG-adjusted concentration, P is the measured concentration, and 1.016 is the mean SG concentration for all urine samples in this subset of the study population. Concentrations below the LOD were substituted by a value equal to the LOD divided by square root of 2 before SG-adjustment . During pregnancy, participants provided one urine sample at each trimester. During preconception, 82% (139/170) of the participants had multiple urine samples. Thus, preconception BP-3 concentrations were calculated as the geometric means of SG-adjusted BP-3 concentrations per participant. For each time point, SG-adjusted BP-3 concentrations were categorized into quartiles. The lowest quartile was used as referent category in all analyses. Due to a skewed distribution, blood glucose levels were log-transformed to achieve normality and then back transformed for better interpretability.
Baseline characteristics of the participants in this study were summarized across quartiles of SG-adjusted BP-3 concentrations in the form of mean ± standard deviation (SD) or number (percentages, %), with p-values obtained by Kruskal-Wallis test for the continuous variables or Fisher’s exact test for the categorical variables. Time-specific distribution of urinary BP-3 concentrations were summarized by calculating geometric means (95% CIs) and percentiles in the study. As a comparison, we calculated the weighted geometric means (95% CIs) and percentiles for BP-3 concentrations from 2009 to 2016 cycles of the National Health and Nutrition Examination Survey (NHANES) among pregnant women and among women aged 23–47 years (comparable to our study population). To identify potential seasonal patterns of SG-adjusted BP-3 concentrations, geometric means (95% CI) during each season were calculated and compared using linear regression (accounted for repeated measurements).
Multivariable linear regression models were used to identify the association between quartiles of SG-adjusted BP-3 concentrations in separate time periods and log-transformed blood glucose levels, adjusted for potential confounders. Population marginal means of glucose levels were used to reveal the population average glucose levels for each quartile of BP-3 concentrations. The median log SG-adjusted bezophenone-3 concentration in each quartile were applied to the models as a continuous variable to test for linear trend. We also used multivariable logistic regression models to identify the association between time-specific BP-3 concentrations and odds of abnormal GLT, defined as glucose level ≥ 140 mg/dL.
According to a priori knowledge with evidence of being risk factors for GDM [31, 32] and potentially associated with BP-3, we posited that maternal age, race/ethnicity, education level, physical activity, smoking status, family history of diabetes, seasonality, maternal BMI, PCOS status, number of fetus, sex of fetus (es) and infertility diagnosis and treatment were potential confounders. Models included the following covariates: maternal age (years), pre-pregnancy BMI (< 25 kg/m2, ≥25 kg/m2), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), infertility diagnosis (male-factor, female-factor, unexplained), physician-diagnosed PCOS (yes, no), infertility treatment (IVF, IUI, natural), fetus number (1, ≥2), sex of fetus (es) (male, female or unknown) and season of urine collection (spring, summer, fall, winter). We also evaluated potential effect modification by season (summer versus other), maternal age (< 37 versus ≥37 years), maternal BMI (< 25 versus ≥25 kg/m2), infertility diagnosis (male-factor, female-factor, versus unexplained), sex of fetus (es) (male versus female), and physical activity (< 5.5 versus ≥5.5 h/week) by adding interaction terms with BP-3 in multivariable models, as well as conducting stratified analyses.
To determine the robustness of results, we performed several sensitivity analyses. First, given the seasonal pattern of BP-3 concentrations in our study, we evaluated associations with and without seasonality in the model. Since outdoor activity may be associated with sunscreen use and glucose levels, we also further adjusted for moderate/heavy outdoor work. We then included time difference between the date of urine collection and GLT in the models. Since the collection of preconception urine samples preceded determination of the sex of fetus (es), to assess any potential mediation by sex, we also evaluated associations between preconception BP-3 and glucose without adjusting for sex of fetus (es). Moreover, we excluded urine samples collected on the same day of GLT (n = 10) and repeated the analyses to address potential issues of temporality. Lastly, we averaged BP-3 concentrations from the first and second trimesters and repeated the analyses. All statistical analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC). Two-sided p-values < 0.05 were considered statistically significant unless otherwise specified.
Baseline characteristics of the study population were summarized in Table 1. Participants had an average age of 34.8 years. Most were non-Hispanic white (82%), never-smokers (76%), and college graduates or higher (86%). 66% had BMI less than 25 kg/m2, 44% had unexplained infertility at baseline, and 60% underwent IVF treatment. The majority delivered a singleton (80%), and 52% delivered female fetus (es). The mean blood glucose level from the GLT was 116.3 mg/dL (SD: 28.3). In general, baseline characteristics were not significantly different comparing the highest quartile of BP-3 to the lowest quartile, except for infertility diagnosis and time difference between urine collection and GLT, where women in the highest quartile of first trimester BP-3 were less likely to be diagnosed with male-factor infertility, and those in the highest quartile of first/second trimester BP-3 tended to have longer time differences between urine collection and GLT.
The distribution of time-specific BP-3 concentrations was shown in Table 2. Detection frequencies were 99.6%, 99.5% and 98.8% for preconception, first trimester, and second trimester, respectively. Overall, the geometric means of SG-adjusted BP-3 concentrations were similar across time windows [geometric mean (95% CI) in ng/mL = 166.6 (130.1, 213.4), 142.4 (109.5, 185.2), 189.1 (139.4, 256.5), respectively]. Compared to the general U.S. population from the 2009–2016 NHANES database which included only pregnant women (BP-3 geometric mean = 42.9 ng/mL) and women with similar age to our study regardless of pregnancy status (BP-3 geometric mean = 39.4 ng/mL), the concentrations in our study population were much higher (See additional file: Table S1). Because sunscreen is the most common source of BP-3, whose use may be greatly dependent on season, the distributions of BP-3 concentrations across seasons were calculated (see Additional file 1: Figure S1). The geometric means (95% CI) in ng/mL of SG-adjusted BP-3 concentrations during summer [286.7 (203.2, 404.6)] were significantly higher than winter [120.3 (92.3, 156.8)] in all urine samples in this study population (p-value < 0.0001).
In Table 3, we present the association between time-specific BP-3 concentrations and pregnancy glucose levels. We did not observe any significant association between preconception or second trimester BP-3 concentrations and glucose levels (preconception: p-trend = 0.29; second trimester: p-trend = 0.28). However, we found that when assessing BP-3 concentrations during the first trimester and glucose levels, the adjusted average glucose levels among women in the highest BP-3 quartile were significantly lower than that among the lowest BP-3 quartile (103.4 versus 114.6 mg/dL, p-value = 0.02, p-trend = 0.07).
When assessing the association between quartiles of BP-3 and odds of abnormal GLT, there were no significant associations between preconception or first trimester BP-3 concentrations (p-trend = 0.49 and 0.28 respectively) and odds of abnormal GLT. When comparing higher quartile of second trimester BP-3 to the lowest quartile, Q3 vs. Q1 during second trimester was associated with lower odds of abnormal GLT [odds ratio (95% CI) = 0.12 (0.01, 0.94), p-value = 0.04), despite the overall non-significant trend (p-trend = 0.90) and no significant differences between Q4 vs. Q1 [odds ratio (95% CI) = 1.04 (0.25, 4.42)].
For potential heterogeneity of these associations, we observed effect modification by infertility diagnosis for first and second trimester BP-3 concentrations and glucose levels (p for interaction = 0.03 and 0.08, respectively). Specifically, we observed inverse associations between BP-3 concentrations in all time windows and glucose levels only among those with female-factor infertility (p-trend = 0.06, 0.01, and 0.005 during preconception, first and second trimester respectively, Fig. 1, see Additional file 1: Table S2). We also found effect modification by season of urine collection for the association between preconception BP-3 and blood glucose level (p for interaction = 0.005, Fig. 2, see Additional file 1: Table S3), where preconception BP-3 concentrations in urine samples collected during summer were inversely associated with blood glucose levels (p-trend < 0.0001) but not for other seasons. Associations were not significantly modified (p for interaction > 0.10) by age, BMI, sex of fetus (es) (see Additional file 1: Table S4 – S6), or physical activity (data not shown), but in stratified analyses there were inverse associations between glucose levels and: (1) preconception BP-3 among women aged ≥37 (p-trend = 0.04); (2) first trimester BP-3 among women with BMI < 25 kg/m2 (p-trend = 0.03); and (3) first trimester BP-3 among women with female fetus (es) (p-trend = 0.009).
The robustness of results was evaluated through sensitivity analyses (data not shown). First, excluding seasonality or sex of fetus (es), or including outdoor work, year of urine collection, and time difference from urine collection to GLT in the multivariable models had very little impact on all results. Secondly, when excluding urine samples collected on the same day of GLT to address for temporality, similar patterns of glucose levels across BP-3 quartiles remained. Lastly, when assessing average BP-3 concentrations from women with both first and second trimesters samples (n = 153), the adjusted average glucose levels among those in the highest BP-3 quartile (99.2 mg/dL) were still significantly lower than in the lowest quartile (110.1 mg/dL, p-value = 0.04), with a significant linear trend across quartiles (p-trend = 0.04).
In this study population of subfertile women, we found time-specific inverse associations between BP-3 urinary concentrations and blood glucose levels measured at late pregnancy. While higher first trimester BP-3 concentrations were associated with significantly lower mean glucose levels, preconception and second trimester BP-3 concentrations did not show this association with glucose levels. Interestingly, some negative associations were stronger among women with female-factor infertility, urine samples collected during summer, older age, lower BMI, or carried female fetus (es). Our results suggest that exposure to BP-3 might potentially decrease glucose levels during pregnancy among women seeking infertility treatment.
There have been few epidemiological studies evaluating exposures to BP-3 and diabetes-related outcomes. A case-control study in Saudi Arabia found that the urinary geometric mean concentrations of BP-3 were 4.2 times greater in the non-diabetic control group compared to the type 2 diabetes group , which is consistent with our findings of inverse associations between BP-3 and glucose levels. However, this study focused on type 2 diabetes in the general population instead of GDM or pregnancy glucose levels, and their sample size was relatively small (54 cases and 47 controls). One study from Denmark found positive associations between a benzophenone derivative and IGF-I and binding protein IGFBP3 among women pregnant with male fetuses , suggesting the potential for benzophenone to regulate insulin metabolism. However, they did not directly measure GDM or glucose levels as outcomes, and BP-3 was measured in serum. Another study from Puerto Rico showed an inverse association between BP-3 and c-reactive protein , a marker for systematic inflammation commonly perceived as a predictor for diabetes ; but again, their outcome was not glucose levels or diabetes, and the study was cross-sectional.
Several epidemiological studies examined the association between other EDCs and GDM or pregnancy glucose levels. Serum POPs concentrations were found to be positively associated with GDM risk . In addition, previous studies conducted among pregnant women from the EARTH Study showed that urinary bisphenol A was positively associated with glucose levels , urinary butylparaben was positively while propylparaben inversely associated with glucose levels , and urinary phthalates also showed associations with glucose levels in mixed directions where a positive association was seen between monoethyl phthalate and glucose but an inverse association was seen between mono-isobutyl phthalate and glucose . The association between BP-3 concentrations and pregnancy glucose levels, however, remains unclear. As such, additional work on this aspect is needed to increase the understanding of potential impact of another highly ubiquitous EDC on GDM risk factors.
Our results showed inverse associations between time-specific BP-3 and pregnancy glucose levels. Potential biological mechanisms for the effect of BP-3 on glucose remain unclear and inconsistent. Conflicting animal studies demonstrated that administration of benzophenone led to either decreased or increased glucose levels among rats [34, 35]. One possible mechanism may involve the reproductive-hormone-like potency of BP-3. BP-3 has been reported to possess estrogenic and anti-androgenic effects, but anti-estrogenic activities were also reported with inverse associations between mixtures of BP-3 and benzophenone-1 and estradiol levels among healthy premenopausal women [36, 37]. Moderate estrogens are known to increase insulin sensitivity, but estrogens at supraphysiological levels may induce insulin resistance . As such, it is still unclear in which direction BP-3 would act on the estrogen levels among pregnant women in late pregnancy. Future studies are needed to help identify the potential effect of BP-3 on estrogen and glucose levels among pregnant women. Another possible mechanism may involve thyroid hormone regulations. Triiodothyronine (T3) was found to be positively associated with risk factors for GDM . An animal study indicated that BP-3 could down-regulate genes related to thyroid stimulating hormones with significant decreases in T3 levels , and similar associations were found among pregnant women . Therefore, it is possible that exposure to BP-3 could lead to decreased level of T3 and decreased GDM risk, potentially via decreasing endogenous glucose production and plasma glucose . Again, future studies are needed to explore these proposed mechanisms. Studies have also shown borderline associations between BP-3 and enhanced glutathione peroxidase activity , which might lead to enhanced anti-inflammatory response and prevention of β-cell dysfunction [43, 44].
We found significant effect modification by infertility diagnosis, where BP-3 concentrations were inversely associated with blood glucose levels specifically among those diagnosed with female-factor infertility. While this finding could be due to chance given the multiple comparisons, it is also possible that some intrinsic mechanisms of BP-3 interacting with infertility might lead to reduced risk of GDM, including potential reductions in the levels of inflammatory markers (e.g. CRP or reactive oxygen species ) for pre-existing inflammation related to infertility. More research is needed to understand how BP-3 may act differently on glucose levels according to fertility conditions.
Although we found year-round high urinary concentrations of BP-3, they were especially high during summer, and interestingly the association between BP-3 concentrations and glucose levels were stronger during summer. In two other studies including women from the same study cohort, urinary BP-3 concentrations were positively associated with sunscreen use during the past 24 h , and sunscreen use was found to be inversely associated with glucose levels during pregnancy , which further strengthened our results as sunscreen is a major source of BP-3. The effect modification by season might suggest unmeasured confounding of sun exposure and vitamin D levels. It is possible that women with higher BP-3 concentrations tend to apply more sunscreens because of increased outdoor activity, therefore gaining more sun exposure that could impact levels of vitamin D in the body, while higher early pregnancy vitamin D levels were shown to be associated with decreased GDM risk . This may also be the reason why we saw inverse associations in women with lower BMI, as this subgroup could spend more time outdoors and have a healthier lifestyle, but we are uncertain why the inverse associations were stronger in women with older age. Due to limitations of our study, we were unable to evaluate this in greater detail and more studies are needed to understand the impact of time spent outdoors, vitamin D levels, and vitamin supplements on associations between BP-3 and glucose levels. In sex-stratified analyses, an inverse association between first trimester BP-3 concentrations and glucose levels was found specifically among women who carried female fetus (es). Although still subject to sample size limitation, this may still point to potential sex-specific effects of BP-3 during pregnancy, possibly through a joint effect of BP-3 and fetus sex on glucose regulation [49, 50]. Further research is needed to identify possible drivers of sex-specific effects of BP-3 during pregnancy.
Our study has several strengths. First, our study is the first to report the associations between BP-3 and glucose levels during pregnancy in a modest-sized subfertile cohort, with measures of BP-3 in multiple time periods before and during pregnancy to evaluate potential sensitive time windows of exposure. Secondly, for preconception, we averaged multiple urine samples for women who had more than one sample, which accounted for possible variability in the exposure measurement. Thirdly, we accounted for a diversity of potential confounders, including lifestyle factors, seasonality, and reproductive conditions, as well as evaluated possible effect modifiers. Fourthly, we conducted sensitivity analysis removing urine samples collected on the same day of GLT, to limit the possibility of reverse causation through a prospective design. Finally, we evaluated this association in a subfertile population at high risk of developing GDM.
Despite these strengths several limitations exist. First, urine samples were convenience samples during clinical visits. Therefore, we cannot rule out the possibility of non-differential exposure misclassification that may bias results toward the null. However, we found in sensitivity analysis that average pregnancy BP-3 concentrations (i.e. geometric mean from first and second trimesters) had similar inverse associations with glucose levels as seen at each individual time period. Second, we were unable to evaluate clinical diagnosis of GDM as an outcome since only 6 women had this condition. Third, due to limited sample size, the results in stratified analyses may be due to chance and should be interpreted with caution. Fourth, our study was based on a high-risk population of women seeking fertility treatment; therefore, our results may not be generalizable to healthier women without fertility problems. Also, BP-3 concentrations in our study population were much higher than pregnant women in a representative subset of the U.S. population during comparable time periods. Even the lowest quartile in our study would likely represent relatively higher concentrations in the general U.S. population. Therefore, we may not be able to capture the full dose-response relationship, specifically between very low BP-3 concentrations and glucose levels. One possible explanation for such big differences in urinary BP-3 concentrations may be differences in the proportion of non-Hispanic white women between both study populations (our study ~ 80% vs. NHANES ~ 40%) . However, demographic characteristics in our study are comparable to nationwide fertility clinics , providing insights to women with a higher baseline risk of GDM due to fertility issues. Lastly, residual (e.g. misclassification of infertility diagnosis) or unmeasured confounding (e.g. by dietary factors, outdoor activities, or vitamin D levels) is still possible. Further studies are needed to replicate, as well as account for more detailed nutritional/behavioral factors.
In conclusion, time-specific associations between higher urinary BP-3 concentrations and lower glucose levels were found among subfertile women. Associations were stronger among women with female-factor infertility, urine samples collected during summer, older age, lower BMI, or carried female fetus (es). While future studies are needed for further understanding of these associations, our results suggest that BP-3, an EDC with ubiquitous existence in consumer products, might affect glucose levels during pregnancy among a group of women at higher risk of GDM.
Availability of data and materials
The datasets in this study are not publicly available due to institutional restrictions prohibiting the sharing of data that contain information that could compromise research participant’s privacy upon which participant consent was contingent. However, the data are available from the corresponding author, Z.W., and the PI of the EARTH Study, R.H. upon reasonable request.
Assisted reproductive technique
Body mass index
The Centers for Disease Control and Prevention
The Environment and Reproductive Health Study
Endocrine disrupting chemical
Gestational diabetes mellitus
Glucose load test
Insulin-like growth factor 1
In vitro fertilization
Limit of detection
Massachusetts General Hospital
National Health and Nutrition Examination Survey
Polycystic ovarian syndrome
Persistent organic pollutant
Quartiles (1st to 4th)
Casagrande SS, Linder B, Cowie CC. Prevalence of gestational diabetes and subsequent type 2 diabetes among U.S. women. Diabetes Res Clin Pract. 2018;141:200–8.
Ashrafi M, Gosili R, Hosseini R, Arabipoor A, Ahmadi J, Chehrazi M. Risk of gestational diabetes mellitus in patients undergoing assisted reproductive techniques. Eur J Obstet Gynecol Reprod Biol. 2014;176:149–52.
Holst S, Kjær SK, Jørgensen ME, Damm P, Jensen A. Fertility problems and risk of gestational diabetes mellitus: a nationwide cohort study. Fertil Steril. 2016;106:427–34.
Thoma ME, McLain AC, Louis JF, King RB, Trumble AC, Sundaram R, Louis GM. Prevalence of infertility in the United States as estimated by the current duration approach and a traditional constructed approach. Fertil Steril. 2013;99:1324–31.
HAPO Study Cooperative Research Group. Hyperglycemia and adverse pregnancy outcomes. N Engl J Med. 2008;358:1991–2002.
Riskin-Mashiah S, Younes G, Damti A, Auslender R. First-trimester fasting hyperglycemia and adverse pregnancy outcomes. Diabetes Care. 2009;32:1639–43.
Zhang C, Rawal S, Chong YS. Risk factors for gestational diabetes: is prevention possible? Diabetologia. 2016;59:1385–90.
Ehrlich S, Lambers D, Baccarelli A, Khoury J, Macaluso M, Ho SM. Endocrine disruptors: a potential risk factor for gestational diabetes mellitus. Am J Perinatol. 2016;33:1313–8.
Eslami B, Naddafi K, Rastkari N, Rashidi BH, Djazayeri A, Malekafzali H. Association between serum concentrations of persistent organic pollutants and gestational diabetes mellitus in primiparous women. Environ Res. 2016;151:706–12.
Chiu YH, Mínguez-Alarcón L, Ford JB, Keller M, Seely EW, Messerlian C, Petrozza J, Williams PL, Ye X, Calafat AM, Hauser R. Trimester-specific urinary bisphenol a concentrations and blood glucose levels among pregnant women from a fertility clinic. J Clin Endocrinol Metab. 2017;102:1350–7.
Bellavia A, Chiu YH, Brown FM, Mínguez-Alarcón L, Ford JB, Keller M, Petrozza J, Williams PL, Ye X, Calafat AM, Hauser R. Urinary concentrations of parabens mixture and pregnancy glucose levels among women from a fertility clinic. Environ Res. 2019;168:389–96.
James-Todd TM, Chiu YH, Messerlian C, Mínguez-Alarcón L, Ford JB, Keller M, Petrozza J, Williams PL, Ye X, Calafat AM, Hauser R. Trimester-specific phthalate concentrations and glucose levels among women from a fertility clinic. Environ Health. 2018;17:55.
Tchernof A. Sex differences in energy balance, body composition, and body fat distribution. In: Tsatsoulis WJ, Brown FM, editors. Diabetes in women: Pathphysiology and therapy. New York: US: Humana press/springer science+business media; 2009.
Han C, Lim YH, Hong YC. Ten-year trends in urinary concentrations of triclosan and benzophenone-3 in the general U.S. population from 2003 to 2012. Environ Pollut. 2016;208:803–10.
Mortensen ME, Calafat AM, Ye X, Wong LY, Wright DJ, Pirkle JL, Merrill LS, Moye J. Urinary concentrations of environmental phenols in pregnant women in a pilot study of the National Children’s study. Environ Res. 2014;129:32–8.
Calafat AM, Wong LY, Ye X, Reidy JA, Needham LL. Concentrations of the sunscreen agent benzophenone-3 in residents of the United States: National Health and nutrition examination survey 2003-2004. Environ Health Perspect. 2008;116:893–7.
Alonso-Magdalena P, Morimoto S, Ripoll C, Fuentes E, Nadal A. The estrogenic effect of bisphenol a disrupts pancreatic β-cell function in vivo and induces insulin resistance. Environ Health Perspect. 2006;114:106–12.
Holland N, Huen K, Tran V, Street K, Nguyen B, Bradman A, Eskenazi B. Urinary phthalate metabolites and biomarkers of oxidative stress in a Mexican-American cohort: variability in early and late pregnancy. Toxics. 2016;4:7.
Artacho-Cordón F, Ríos-Arrabal S, León J, Frederiksen H, Sáenz JM, Martín-Olmedo P, Fernández MF, Olea N, Arrebola JP. Adipose tissue concentrations of non-persistent environmental phenols and local redox balance in adults from southern Spain. Environ Int. 2019;133:105118.
Krause M, Frederiksen H, Sundberg K, Jørgensen FS, Jensen LN, Nørgaard P, Jørgensen C, Ertberg P, Petersen JH, Feldt-Rasmussen U, Juul A. Maternal exposure to UV filters: associations with maternal thyroid hormones, IGF-I/IGFBP3 and birth outcomes. Endocr Connect. 2018;7:334–46.
Watkins DJ, Ferguson KK, Del Toro LV, Alshawabkeh AN, Cordero JF, Meeker JD. Associations between urinary phenol and paraben concentrations and markers of oxidative stress and inflammation among pregnant women in Puerto Rico. Int J Hyg Environ Health. 2015;218:212–9.
Ferguson KK, Lan Z, Yu Y, Mukherjee B, McElrath TF, Meeker JD. Urinary concentrations of phenols in association with biomarkers of oxidative stress in pregnancy: assessment of effects independent of phthalates. Environ Int. 2019;131:104903.
Li AJ, Xue J, Lin S, Al-Malki AL, Al-Ghamdi MA, Kumosani TA, Kannan K. Urinary concentrations of environmental phenols and their association with type 2 diabetes in a population in Jeddah, Saudi Arabia. Environ Res. 2018;166:544–52.
Philippat C, Nakiwala D, Calafat AM, Botton J, De Agostini M, Heude B, Slama R. EDEN mother–child study group. Prenatal exposure to nonpersistent endocrine disruptors and behavior in boys at 3 and 5 years. Environ Health Perspect. 2017;125:097014.
Freemark M. Regulation of maternal metabolism by pituitary and placental hormones: roles in fetal development and metabolic programming. Horm Res Paediatr. 2006;65:41–9.
Messerlian C, Williams PL, Ford JB, Chavarro JE, Mínguez-Alarcón L, Dadd R, Braun JM, Gaskins AJ, Meeker JD, James-Todd T, Chiu YH. The Environment and Reproductive Health (EARTH) Study: a prospective preconception cohort. Hum Reprod Open. 2018;2018:hoy001.
Pearson MA, Lu C, Schmotzer BJ, Waller LA, Riederer AM. Evaluation of physiological measures for correcting variation in urinary output: implications for assessing environmental chemical exposure in children. J Expo Sci Environ Epidemiol. 2009;19:336.
Silva MJ, Samandar E, Preau JL Jr, Reidy JA, Needham LL, Calafat AM. Quantification of 22 phthalate metabolites in human urine. J Chromatogr B. 2007;860:106–12.
Gupta Y, Kalra B, Baruah M, Singla R, Kalra S. Updated guidelines on screening for gestational diabetes. Int J Women's Health. 2015;7:539.
Hornung RW, Reed LD. Estimation of average concentrations in the presence of nondetectable values. Appl Occup Environ Hyg. 1990;5:46–51.
Mirghani Dirar A, Doupis J. Gestational diabetes from a to Z. World J Diabetes. 2017;8:489–511.
Chiefari E, Pastore I, Puccio L, Caroleo P, Oliverio R, Vero A. P Foti D, Vero R, Brunetti a. impact of seasonality on gestational diabetes mellitus. Endocr Metab Immune Disord Drug Targets. 2017;17:246–52.
Oda E. High-sensitivity C-reactive protein, but not white blood cell count, independently predicted incident diabetes in a Japanese health screening population. Acta Diabetol. 2015;52:983–90.
Yamasaki K, Okuda H, Kikkawa H. Endocrine-mediated effects of a benzophenone related chemical, 2,3,4,4′-Tetrahydroxybenzophenone, based on uterotrophic assay, hershberger assay, and subacute oral toxicity study. J Clin Toxicol. 2012;2:129.
Burdock GA, Ford RA, Cliffs E, Manufacturers E. Safety evaluation of benzophenone. Food Chem Toxicol. 1991;29:741–50.
Kim S, Choi K. Occurrences, toxicities, and ecological risks of BP-3, a common component of organic sunscreen products: a mini-review. Environ Int. 2014;70:143–57.
Pollack AZ, Mumford SL, Krall JR, Carmichael AE, Sjaarda LA, Perkins NJ, Kannan K, Schisterman EF. Exposure to bisphenol a, chlorophenols, benzophenones, and parabens in relation to reproductive hormones in healthy women: a chemical mixture approach. Environ Int. 2018;120:137–44.
Mauvais-Jarvis F, Clegg DJ, Hevener AL. The role of estrogens in control of energy balance and glucose homeostasis. Endocr Rev. 2013;34:309–38.
Rawal S, Tsai MY, Hinkle SN, Zhu Y, Bao W, Lin Y, Panuganti P, Albert PS, Ma RC, Zhang C. A longitudinal study of thyroid markers across pregnancy and the risk of gestational diabetes. J Clin Endocrinol Metab. 2018;103:2447–56.
Lee J, Kim S, Park YJ, Moon HB, Choi K. Thyroid hormone-disrupting potentials of major benzophenones in two cell lines (GH3 and FRTL-5) and embryo-larval zebrafish. Environ Sci Technol. 2018;52:8858–65.
Aker AM, Johns L, McElrath TF, Cantonwine DE, Mukherjee B, Meeker JD. Associations between maternal phenol and paraben urinary biomarkers and maternal hormones during pregnancy: a repeated measures study. Environ Int. 2018;113:341–9.
Klieverik LP, Janssen SF, van Riel A, Foppen E, Bisschop PH, Serlie MJ, Boelen A, Ackermans MT, Sauerwein HP, Fliers E, Kalsbeek A. Thyroid hormone modulates glucose production via a sympathetic pathway from the hypothalamic paraventricular nucleus to the liver. Proc Natl Acad Sci. 2009;106:5966–71.
Li C, Deng X, Xie X, Liu Y, Friedmann Angeli JP, Lai L. Activation of glutathione peroxidase 4 as a novel anti-inflammatory strategy. Front Pharmacol. 2018;9:1120.
Koulajian K, Ivovic A, Ye K, Desai T, Shah A, George Fantus I, Ran Q, Giacca A. Overexpression of glutathione peroxidase 4 prevents β-cell dysfunction induced by prolonged elevation of lipids in vivo. Am J Physiol Endocrinol Metab. 2013;305(2):E254–62.
Qiu C, Sorensen TK, Luthy DA, Williams MA. A prospective study of maternal serum C-reactive protein (CRP) concentrations and risk of gestational diabetes mellitus. Paediatr Perinat Epidemiol. 2004;18:377–84.
Mínguez-Alarcón L, Chiu YH, Nassan FL, Williams PL, Petrozza J, Ford JB, Calafat AM, Hauser R, Chavarro JE. Earth Study Team. Urinary concentrations of benzophenone-3 and reproductive outcomes among women undergoing infertility treatment with assisted reproductive technologies. Sci Total Environ. 2019;678:390–8.
Bellavia A, Mínguez-Alarcón L, Ford JB, Keller M, Petrozza J, Williams PL, Hauser R, James-Todd T. EARTH study team. Association of self-reported personal care product use with blood glucose levels measured during pregnancy among women from a fertility clinic. Sci Total Environ. 2019;695:133855.
Arnold DL, Enquobahrie DA, Qiu C, Huang J, Grote N, VanderStoep A, Williams MA. Early pregnancy maternal vitamin D concentrations and risk of gestational diabetes mellitus. Paediatr Perinat Epidemiol. 2015;29:200–10.
Retnakaran R, Kramer CK, Ye C, Kew S, Hanley AJ, Connelly PW, Sermer M, Zinman B. Fetal sex and maternal risk of gestational diabetes mellitus: the impact of having a boy. Diabetes Care. 2015;38:844–51.
Jaskolka D, Retnakaran R, Zinman B, Kramer CK. Sex of the baby and risk of gestational diabetes mellitus in the mother: a systematic review and meta-analysis. Diabetologia. 2015;58:2469–75.
Mortensen ME, Caudill SP, Caldwell KL, Ward CD, Jones RL. Total and methyl mercury in whole blood measured for the first time in the US population: NHANES 2011–2012. Environ Res. 2014;134:257–64.
We gratefully acknowledge Xiaoliu Zhou, Tao Jia, and the late Xiaoyun Ye (CDC, Atlanta, GA) for providing technical assistance with BP-3 quantification. We also thank research staff Ramace Dadd and Patricia Morey at Harvard T.H. Chan School of Public Health, physicians and staff at MGH Fertility Center, all members of the EARTH study personnel, and all the study participants.
This work was funded by the National Institute of Environmental Health Sciences R01ES026166, R01ES022955, R01ES009718, P30ES000002.
Ethics approval and consent to participate
The study was approved by the Partners IRB, Harvard IRB, and the Centers for Disease Control and Prevention (CDC). We obtained signed informed consent from all participants.
Consent for publication
None of the authors have conflicts of interest to declare. The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the CDC. Use of trade names is for identification only and does not imply endorsement by the CDC, the Public Health Service, or the U.S. Department of Health and Human Services.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Distribution of urinary BP-3 concentrations (ng/mL) among NHANES women during 2009–2016. Legend: Abbreviations: BP-3, benzophenone-3; LOD, limit of detection (0.4 ng/ml in NHANES 2009–2016). All estimates incorporated NHANES multistage sub-sample weights for combined survey cycles. BP-3 concentrations below LOD were assigned a value equal to the LOD divided by square root of 2 in NHANES.a Women in NHANES 2009–2016 who had positive lab pregnancy test results or self-reported pregnant at the time of urine sample collection, with available BP-3 concentrations. Due to disclosure risks, pregnancy status was only released for women 20–44 years of age since 2007.b Women in NHANES 2009–2016 who were 23–47 years old, regardless of pregnancy status. Table S2. Quartiles of time-specific BP-3 urinary concentrations and blood glucose levels among pregnant women in the EARTH Study: stratified by infertility diagnosis. Legend: Abbreviations: LOD, limit of detection (0.4 ng/ml for some years and 0.2 ng/ml for other years). *Adjusted for maternal age at pregnancy (years), pre-pregnancy BMI (< 25 kg/m2, ≥25 kg/m2), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), physician-diagnosed PCOS (yes, no), fetus number in a pregnancy (1, ≥2), sex of fetus (es) (male, female or missing) and season (spring, summer, fall, winter). † Test for linear trend were performed using the median SG-adjusted urinary BP-3 concentration in each quartile as a continuous variable in the model, adjusted for the above covariates. ‡ Adjusted p for interaction was obtained by adding interaction terms of infertility diagnosis*the median SG-adjusted urinary BP-3 concentration in each quartile to the models described in †. § p-value for comparison against Q1 is < 0.05. Table S3. Quartiles of time-specific BP-3 urinary concentrations and blood glucose levels among pregnant women in the EARTH Study: stratified by season of urine collection. Legend: Abbreviations: LOD, limit of detection (0.4 ng/ml for some years and 0.2 ng/ml for other years). *Adjusted for maternal age at pregnancy (years), pre-pregnancy BMI (< 25 kg/m2, ≥25 kg/m2), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), infertility diagnosis (male factor, female factor, unexplained), physician-diagnosed PCOS (yes, no), type of infertility treatment (IVF, IUI, natural), fetus number in a pregnancy (1, ≥2), and sex of fetus (es) (male, female or missing). † Test for linear trend were performed using the median SG-adjusted urinary BP-3 concentration in each quartile as a continuous variable in the model, adjusted for the above covariates. ‡ Adjusted p for interaction was obtained by adding interaction terms of season*the median SG-adjusted urinary BP-3 concentration in each quartile to the models described in †. § p-value for comparison against Q1 is < 0.05. Table S4. Quartiles of time-specific BP-3 urinary concentrations and blood glucose levels among pregnant women in the EARTH Study: stratified by maternal age. Legend: Abbreviations: LOD, limit of detection (0.4 ng/ml for some years and 0.2 ng/ml for other years). *Adjusted for pre-pregnancy BMI (< 25 kg/m2, ≥25 kg/m2), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), infertility diagnosis (male factor, female factor, unexplained), physician-diagnosed PCOS (yes, no), type of infertility treatment (IVF, IUI, natural), fetus number in a pregnancy (1, ≥2), sex of fetus (es) (male, female or missing), and season (spring, summer, fall, winter). † Test for linear trend were performed using the median SG-adjusted urinary BP-3 concentration in each quartile as a continuous variable in the model, adjusted for the above covariates. ‡ Adjusted p for interaction was obtained by adding interaction terms of maternal age*the median SG-adjusted urinary BP-3 concentration in each quartile to the models described in †. § p-value for comparison against Q1 is < 0.05. Table S5. Quartiles of time-specific BP-3 urinary concentrations and blood glucose levels among pregnant women in the EARTH Study: stratified by maternal BMI. Legend: Abbreviations: LOD, limit of detection (0.4 ng/ml for some years and 0.2 ng/ml for other years). * Adjusted for maternal age at pregnancy (years), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), infertility diagnosis (male factor, female factor, unexplained), physician-diagnosed PCOS (yes, no), type of infertility treatment (IVF, IUI, natural), fetus number in a pregnancy (1, ≥2), sex of fetus (es) (male, female or missing), and season (spring, summer, fall, winter). † Test for linear trend were performed using the median SG-adjusted urinary BP-3 concentration in each quartile as a continuous variable in the model, adjusted for the above covariates. ‡ Adjusted p for interaction was obtained by adding interaction terms of maternal BMI*the median SG-adjusted urinary BP-3 concentration in each quartile to the models described in †. § p-value for comparison against Q1 is < 0.05. Table S6. Quartiles of time-specific BP-3 urinary concentrations and blood glucose levels among pregnant women in the EARTH Study: stratified by sex of fetus (es). Legend: Abbreviations: LOD, limit of detection (0.4 ng/ml for some years and 0.2 ng/ml for other years). * Adjusted for maternal age at pregnancy (years), pre-pregnancy BMI (< 25 kg/m2, ≥25 kg/m2), family history of diabetes (yes, no), baseline smoking status (never, ever), total physical activity (hours/week), race (white, non-white), education level (college graduate or higher, other), infertility diagnosis (male factor, female factor, unexplained), physician-diagnosed PCOS (yes, no), type of infertility treatment (IVF, IUI, natural), fetus number in a pregnancy (1, ≥2), and season (spring, summer, fall, winter). † Test for linear trend were performed using the median SG-adjusted urinary BP-3 concentration in each quartile as a continuous variable in the model, adjusted for the above covariates. ‡ Adjusted p for interaction was obtained by adding interaction terms of sex of fetus (es)*the median SG-adjusted urinary BP-3 concentration in each quartile to the models described in †. § p-value for comparison against Q1 is < 0.05. Figure S1. SG-adjusted urinary BP-3 concentrations (in ng/mL) across seasons among 217 pregnant women in the EARTH Study. Legend: Overall p-value obtained from Type 3 analysis in linear regression, accounted for repeated measures of BP-3 concentrations within participants (number of participants = 217, number of samples = 833).
About this article
Cite this article
Wang, Z., Mínguez-Alarcón, L., Williams, P.L. et al. Perinatal urinary benzophenone-3 concentrations and glucose levels among women from a fertility clinic. Environ Health 19, 45 (2020). https://doi.org/10.1186/s12940-020-00598-7
- Glucose levels
- Endocrine disrupting chemicals
- Environmental epidemiology