Skip to main content

Persistent organic pollutants, pre-pregnancy use of combined oral contraceptives, age, and time-to-pregnancy in the SELMA cohort

Abstract

Background

We are exposed to several chemicals such as persistent organic pollutants (POPs) in our everyday lives. Prior evidence has suggested that POPs may have adverse effects on reproductive function by disrupting hormone synthesis and metabolism. While there is age-related decline of fertility, the use of hormonal combined oral contraceptives (COCs) and its association to return of fertility remains controversial. The goal of this study is to investigate the association between exposure to POPs, both individually and as a mixture, and fecundability measured as time-to-pregnancy (TTP) according to pre-pregnancy use of COCs and age.

Methods

Using the SELMA (Swedish Environmental Longitudinal Mother and Child, Allergy and Asthma) study, we have identified 818 pregnant women aged 18–43 years (mean 29 years) with data on how long they tried to get pregnant and what was their most recently used contraceptive method. These data were collected at enrollment to the study (median week 10 of pregnancy). Concentrations of 22 POPs and cotinine were analyzed in the blood samples collected at the same time as the questions on TTP and pre-pregnancy use of contraceptive. Analyses were done on the association between POPs exposure and TTP measured as continuous (months) and binary (infertile for those with TTP > 12 months). To study the chemicals individually, Cox regression and logistic regression were used to estimate fecundability ratios (FRs) and odds ratios (ORs), respectively. Weighted quantile sum (WQS) regression was used to investigate the chemicals as a mixture where chemicals of concern were identified above the 7.6% threshold of equal weights. To perform the subgroup analysis, we stratified the sample according to use of COCs as the most recent pre-pregnancy contraception method and age (< 29 years, and ≥ 29 years). The models were adjusted for parity, regularity of menses, maternal body mass index (BMI) and smoking status, and stratified as described above.

Results

Prior to stratification, none of the POPs were associated with fecundability while increased exposure to HCB, PCB 74 and 118 had higher odds of infertility. Upon stratification, POP exposure was significantly associated with longer TTP in women aged ≥29 years who did not use COC. Specifically, PCBs 156, 180, 183, and 187 were associated with reduced fecundability while PCBs 99, 153, 156, 180, 183, and 187 had higher odds of infertility. As a mixture, we identified the chemicals of concern for a longer TTP include PCBs 118, 156, 183, and 187. Moreover, chemicals of concern identified with increased odds of infertility were PCB 74, 156, 183, 187, and transnonachlor.

Conclusion

Serum concentrations of selected POPs, both as individual chemicals and as a mixture, were significantly associated with lower fecundability and increased odds of infertility in women aged 29 years and above not using COC as their most recent pre-pregnancy contraceptive. Our findings suggest that pre-pregnancy use of oral contraceptive and age may modify the link between POPs and fecundability. The differences of specific chemicals in the individual analysis and as a mixture support the need to study combination effects of chemicals when evaluating reproductive outcomes.

Peer Review reports

Background

Persistent organic pollutants (POPs) such as organochlorine pesticides (OCPs), polychlorinated biphenyls (PCBs), and polybrominated diphenyl ethers (PBDEs) are organic substances that have long half-lives, bioaccumulate in the fatty tissues with increasing concentration towards the top of the food chain, and travel long distances in the atmosphere due to their semi-volatility [1]. They are widely used in agriculture, consumer products, and industrial products as well as unintentionally released as by-products from industrial processes and incineration. Exposure to these chemicals occurs through ingestion, inhalation, and absorption. We are exposed not only to a single chemical but to several chemicals simultaneously, resulting to mixtures that are associated with adverse outcomes [2,3,4]. Although many of these compounds have been strictly regulated in Europe including Sweden through the Stockholm Convention in 2001 [5], they still exist in the environment due to their high resistance to degradation. This has become a global health concern as studies of humans, wildlife populations, and epidemiological cohorts show associations between POP exposure and adverse effects on reproductive and endocrine functions [1, 6].

Since reproduction is regulated by hormones, infertility may occur when the endocrine system is disrupted. Infertility is the inability to achieve clinical pregnancy after 12 months of regular unprotected sexual intercourse [7]. By contrast, fertility is the actual production of a live offspring while fecundity is the ability to conceive given unprotected intercourse. Fecundability is the probability of becoming pregnant and is measured by time-to-pregnancy (TTP), which is the number of menstrual cycles or months without the use of contraceptives until a clinically detectable pregnancy is achieved [8, 9]. TTP has been used as a marker for fecundability of both parents to identify environmental exposures as potential hazards to human reproduction, encompassing series of biological events such as gametogenesis, fertilization, and implantation. Increased TTP may then indicate a problem in any or several of these stages [10].

There are several factors that affect female fertility. While biological factors such as age and BMI are inversely associated with fertility [11, 12], use of combined oral contraceptives (COCs), which is composed of both estrogen and progesterone, have contradicting studies about its association to return to fertility [13,14,15]. In addition, environmental risk factors such as exposure to POPs may also reduce fertility. Detection of environmental chemicals in human follicular fluid, serum, and seminal plasma led to the hypothesis of susceptibility of the reproductive function to disruption of biosynthesis and metabolism of hormones such as estrogen and progesterone [16, 17]. While POPs have been found to have effects on conception rates, menstrual cycles, and birth outcomes in animals such as rhesus monkeys and rodents [18, 19], previous human epidemiological studies have contradicting results. Some studies reported longer TTP for women exposed to POPs [20,21,22]. However, others did not find this association [23, 24]. These inconsistencies warrant further studies to explore the hazards of exposure to POPs and COC use.

Aim of the study

The goal of this study is to investigate the association between exposure to POPs, both individually and as a mixture, during early pregnancy and fecundability measured as time-to-pregnancy (TTP) in regard to pre-pregnancy use of COCs and age.

Methods

Study population

This study utilized the Swedish Environmental Longitudinal, Mother and child, Asthma and allergy (SELMA) study, a pregnancy cohort with the primary aim to investigate the importance of early life exposure to environmental toxicants with focus on endocrine disrupting chemicals (EDCs) during the pregnancy and infancy period for health and development in the children [25]. Pregnant women were recruited from September 2007 to March 2010 during their first visit (median = 10 weeks age of gestation) at antenatal care centers in Värmland, Sweden. Out of 8394 women, 7119 were invited and 2582 (39%) consented to participate. Maternal blood was collected, and questionnaires were used to gather information regarding their medical history, lifestyles, and socioeconomic status. Questions on TTP were asked by midwife at enrollment. From the original cohort, only women with complete data on TTP, serum POP concentrations, age, parity, preventive method, regularity of menses, BMI, and lifestyle factors were included in this study (n = 818) (Supp Fig. 1).

Time-to-pregnancy (TTP)

Fecundability or the probability of pregnancy in each cycle was measured through TTP, a sensitive and convenient marker used to study environmental exposures. Women were asked “How long have you been trying to get pregnant?” with answers in years and months. Additional questions on pre-pregnancy contraceptive use and reproductive history were also inquired. They were asked “What was the preventive method you used before you became pregnant?” and when they stopped using it. They were also asked if they have regular or irregular menses.

Chemical analyses of POPs in serum

Concentrations of 22 POPs (nine OCPs namely pentachlorobenzene (PeCB), hexachlorobenzene (HCB), three isomers of hexachlorocyclohexane (α–HCH, β–HCH, γ–HCH), oxychlordane, transnonachlor, dichlorodiphenyltrichloroethane (p,p´-DDT), dichlorodiphenyldichloroethylene (p,p´-DDE), ten PCB congeners namely, PCB 74, 99, 118, 138, 153, 156, 170, 180, 183, and 187 as well as three PBDEs namely PBDE 47, 99, and 153) were analyzed in the blood samples collected during the first healthcare visit of current pregnancy, median 10 weeks age of gestation, as described previously [26]. Briefly, 200 ul of serum sample were enriched with 400 pg 13C-labelled internal standards of each compound. POPs were extracted with 2 ml dichloromethane-hexane (1:4). Extracts were cleaned with multilayer silica columns and the eluate was concentrated for gas chromatography – tandem mass spectrometry (GC-MS/MS) analysis (Agilent 7010 GC-MS/MS system, Wilmington, DE, USA). Control serum sample (NIST SRM 1958) and an in-house low-concentration control sample (1 to 9 dilution of SRM 1958 with new born calf serum) were included. Concentrations were reported as wet weight (pg/ml) and the limits of quantification (LOQ) ranged from 5 to 40 pg/ml. LOQ was defined as the concentration corresponding to ten times the standard deviation of the signal-to-noise ratio. Of the 22 POPs analyzed, 13 were detected in more than 70% of the included women.

Analysis of cotinine in serum samples

Cotinine was measured in serum samples and used as a biomarker for tobacco smoke exposure. Isotopically labelled internal standards were added to 100 μl serum. Samples were digested with glucuronidase and proteins were precipitated using acetonitrile. Analyses were performed using triple quadrupole mass spectrometry (QTRAP 5500; AB Sciex, Foster City, USA) coupled to a liquid chromatography system (UFLCXR, Shimadzu Corporation, Kyoto, Japan) (LC-MS/MS). Subjects were categorized as non-smokers for cotinine concentrations below 0.2 ng/ml, passive smokers with cotinine concentrations 0.2–15 ng/ml and active smokers for cotinine above 15 ng/ml [27].

Statistical analyses

POPs detected in less than 70% of the samples were not included in the analyses. Among those chemicals detected in more than 70% of the mothers, concentrations below the LOQ were replaced with LOQ/2. Characteristics of the cohort were presented as n (%), unless otherwise specified. Chi-square test, Wilcoxon rank-sum test and t-test test were used to analyze demographics and chemicals among groups. Correlation among chemicals were determined using Spearman’s rank correlation. To compare exposure patterns between Sweden and the United States, the geometric mean of serum concentrations of the 13 POPs were used to determine cumulative burden and its relative proportion in the SELMA cohort as well as in women aged 20–39 years in the National Health and Nutrition Examination Survey (NHANES) in 2007–2010 conducted by Centers for Disease Control and Prevention in the United States [28, 29].

Analyses were done on the association of POPs on the outcome variable TTP, both as continuous (measured in months) and binary (infertile for those with TTP > 12 months). Chemicals were analyzed both as continuous (log10 transformed) and quartiles. Discrete-time Cox regression models were used to estimate hazard ratios of fecundability (FR). The proportional hazards assumption was met based on the statistical tests on the scale Schoenfeld residuals. The FR represents the per-cycle probability of conception in subgroups with higher concentration POP exposure relative to a reference group with lower concentration POP exposure [30]. We censored the data when TTP exceeded 12 months. Sensitivity analyses were done with censored data at 10 and 14 months. An FR > 1 denotes higher fecundability, which means shorter TTP and FR < 1 denotes lower fecundability, which means longer TTP. Moreover, logistic regression was used to estimate odds ratio (OR) for infertility. An OR > 1 denotes higher odds for infertility and OR < 1 denotes lower odds for infertility. OR represents the odds that an outcome happens given a particular exposure compared to the odds of an outcome happening without that exposure [31]. To test the linear trend of the FR and OR, P-trend was calculated using the quartiles in a continuous form.

Aside from analyzing the individual compounds, a mixture approach was also performed. Due to the high correlation among chemicals, survival analysis was not used to explore mixture effects. Instead, weighted quantile sum (WQS) regression, a strategy for estimating empirical weights for a weighted sum of quantiled concentrations most associated with the health outcome (TTP and binary infertility), was performed [32]. Chemicals of concern were identified with a threshold of 7.6% (100/number of chemicals). Bootstrapping was set to 100. The dataset was split into a 35% training set and a 65% holdout validation set.

Because interaction tests between chemicals and covariates showed significant results for age and use of COCs, we stratified the sample based on these two factors, namely age (with the cohort mean 29 years old as cut off) and use of COC (yes, no) as the most recent pre-pregnancy contraception method. Quartiles were specific for each age group to ensure equal n in each quartile. All models were adjusted for well-established risk factors such as parity, education, regularity of menses (regular, irregular), maternal body mass index (BMI), and smoking status based on serum cotinine concentrations as described above. We also tested other covariates such as education (lower secondary, upper secondary/vocational studies, university/college), physical activity (< 1 h/week, 1 h/week, 1–2 h/week, ≥3 h/week), and alcohol intake (never, seldom, once every other week, at least once a week). The only paternal covariate that was tested was BMI. However, none of these other covariates modified the FR and OR by 10% or more. Hence, they were not added to the model.

These analyses were performed using IBM SPSS Statistics, version 22.0 (IBM Corp., Armonk, NY, USA) and the R programming language [33] through RStudio [34] with the packages gWQS [35], ggplot2 [36], ggpubr [37], and reshape2 [38].

Results

Of the 818 women included in this study, the median TTP was 3 months with the interquartile range (25th to 75th percentile) between 1 and 6 months. Age was positively correlated with TTP (rho = 0.11, p = 0.002). Women who did not use COC had shorter TTP compared to women who used COC, both for those aged < 29 years (β = − 1, 95% CI = − 6.4, − 1) and aged ≥29 years (β = − 1, 95% CI = − 2, − 1) (Table 1). Seventy-nine women reported more than 12 months until pregnancy was achieved, resulting in a 9.7% censoring rate. Parity was higher among non-COC user women compared to COC user women. (Table 1).

Table 1 Cohort characteristics. Values are mean (SD), median (min-max) or n (%)

Out of the 22 POPs studied, 13 chemicals, consisting of three OCPs and ten congeners of PCBs, were detected in more than 70% of the 818 women. Concentrations of the 13 chemicals were significantly positively correlated with each other (rho = 0.48–0.98, p < 0.001) (Fig. 1). The cumulative burden of serum POPs in the SELMA cohort varied from 500 to 650 pg/mL in the subgroups and was half of that of the women in the United States (Fig. 2). Similar exposure profiles were seen among all groups in the SELMA cohort, mainly consisting of p,p’-DDE (31%) followed by PCB 153 (18%), PCB 180 (13%), PCB 138 (12%), HCB (8%) and PCB 170 (7%). However, exposure profile in women in NHANES was mainly p,p’-DDE (80%) and HCB (4%) (Fig. 2) [28, 29]. Compared to exposure in women aged 20–39 years in the United States [28, 29], the SELMA cohort had lower concentrations of transnonachlor (p < 0.001), p,p’-DDE (p < 0.001), PCB 74 (p < 0.001), and 99 (p < 0.05) but higher concentrations of PCB 138 (p < 0.001), 153 (p < 0.001), 156 (p < 0.001), 170 (p < 0.001), 180 (p < 0.001), 183 (p < 0.001), and 187 (p < 0.001).

Fig. 1
figure 1

Correlation of POPs in maternal serum. Chemicals included are those detected in more than 70% of the samples. Data are presented as Spearman’s rank correlation rho ranging from −1.0 (dark blue) to 1.0 (dark red). All correlations are significant (p < 0.001). Number of observations can be found in Table 2. Abbreviations: HCB - hexachlorobenzene; p,p´-DDE - dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyl

Fig. 2
figure 2

(a) Cumulative concentration and (b) relative proportion of POPs in NHANES 2007–2010 and SELMA cohort. Geometric means of serum concentrations of the 13 analyzed compounds were used to determine cumulative burden and relative proportion in women aged 20–39 years in NHANES conducted in 2007–2010 in the United States and Swedish SELMA cohort both unstratified and stratified. Abbreviations: POPs – persistent organic pollutants; COC – combined oral contraceptives; HCB - hexachlorobenzene; p,p´-DDE – dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyls; NHANES – National Health and Nutrition Examination Survey

The serum concentrations of POPs between the groups were compared (Table 2). Among women < 29 years old, serum concentrations of HCB, transnonachlor, p,p’-DDE, PCBs 74, 99, 118, 138, 153, 156, 183, and PBDE 153 were lower in non-COC users compared to those of COC users. Among women ≥29 years old, there was no difference in POPs concentrations between non-COC and COC users (Table 2).

Table 2 Persistent organic pollutants (POPs) in maternal serum according to pre-pregnancy use of combined oral contraceptives (COCs) and age. Data are presented as GM (95% CI), pg/ml. Levels below LOQ replaced with LOQ/2

Without age and COC stratification, none of the chemicals (as continuous or quartiles) had significant FR after adjustment with age, pre-pregnancy use of COCs, parity, maternal BMI and smoking status (Suppl Table 1). However, upon stratification, significant associations were found. Among women aged < 29 years who used COC, those exposed within the third quartiles of HCB and PCB 138 had 37% reduced fecundability (i.e. longer TTP) compared to those exposed to their respective first quartiles (FR 0.63; 95% CI 0.4–0.99 and FR 0.63; 95% CI 0.4–0.98, respectively) (Fig. 3). Among women aged ≥29 years who did not use COC, those exposed in the third quartile of PCB 156 (FR 0.67; 95% CI 0.46–0.95), and fourth quartiles of PCB 180 (FR 0.69; 95% CI 0.48–0.98), 183 (FR 0.69; 95% CI 0.49–0.98) and 187 (FR 0.66; 95% CI 0.46–0.95) had longer TTP compared their respective first quartiles (Fig. 3). There was a 47% (FR 0.53; 95% CI 0.28–0.99) and 45% (FR 0.55; 95% CI 0.31–0.98) decrease in fecundability for every 10-fold increase of PCB 180 and 187, respectively (Fig. 3). On the other hand, older women who used COCs had a reduction in fecundability for those in the second and fourth quartiles of p,p’-DDE (FR 0.56; 95% CI 0.36–0.87 and FR 0.65; 95% CI 0.42–0.99, respectively), and second quartile of PCB 183 (FR 0.60; 95% CI 0.39–0.93) compared to their quartiles (Fig. 3). Sensitivity analyses with censored data at 10 and 14 months showed similar results (Supp Figs. 2 and 3).

Fig. 3
figure 3

Associations of POPs exposure with TTP using Cox regression. Dot shows fecundability ratio (FR) with 95% confidence interval error bars of chemicals as continuous (log) and quartiles (Q1 as reference). FR > 1 denotes higher fecundability and shorter TTP while FR < 1 denotes lower fecundability and longer TTP. FR indicated with * for p value < 0.05. P trend shown for significant linear trend of FR. Abbreviations: COC – combined oral contraceptives; HCB - hexachlorobenzene; p,p´-DDE – dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyls; Q1 – first quartile, Q2 – second quartile, Q3 – third quartile, Q4 – fourth quartile

Prior to age and COC stratification, every 10-fold increase in exposure to PCB 74 and 118 had a 3-fold (OR 3; 95% CI 1.33, 6.79) and 4-fold (OR 4.06; 95% CI 1.28, 12.88) odds to be infertile, respectively (Suppl Table 1). The third quartile of HCB had 125% (OR 2.25; 95% CI 1.06, 4.78) higher odds to be infertile compared to the first quartile (Suppl Table 1). After stratification, every 10-fold increase in exposure to PCB 187 had an 8-fold higher odds for infertility in women aged ≥29 years old who did not use COC (OR 7.99; 95% CI 1.02, 62.83) (Fig. 4). There were also higher odds for infertility for those in the second quartile of PCB 99 (OR 5.62; 95% CI 1.15, 27.5), fourth quartile of PCB 153 (OR 4.53; 95% CI 1.23, 16.64), third quartile of PCB 156 (OR 5.04; 95% CI 1.47, 17.22), fourth quartile of PCB 180 (OR 6.08; 95% CI 1.42, 26.07), fourth quartile of PCB 183 (OR 4.38; 95% CI 1.11, 17.21), and third and fourth quartiles of PCB 187 (OR 3.87; 95% CI 1.11, 13.5 and OR 4.26; 95% CI 1.17, 15.48) (Fig. 4). However, none of the chemicals was significantly associated to infertility in the three other groups of women (Fig. 4).

Fig. 4
figure 4

Associations of POPs exposure with infertility using logistic regression. Dot shows odds ratio (OR) with 95% confidence interval (CI) error bars of chemicals as continuous (log) and quartiles (Q1 as reference). OR > 1 denotes higher odds for infertility while OR < 1 denotes lower odds for infertility. OR indicated with * for p value < 0.05. P trend shown for significant linear trend of OR. Abbreviations: COC – combined oral contraceptives; HCB - hexachlorobenzene; p,p´-DDE – dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyls; Q1 – first quartile, Q2 – second quartile, Q3 – third quartile, Q4 – fourth quartile

Because the chemicals were highly correlated (Fig. 1), WQS regression was performed to study the exposure to POPs as a mixture. The WQS-index was not significantly associated with TTP nor infertility prior to stratification. However, after stratification, there was a 9% increase of the log transformed TTP for every one quartile increase of the WQS-index in women aged ≥29 years old not using COCs (β = 0.09; SE = 0.04; p < 0.05) (Fig. 5 and Table 3). Identified chemicals of concern (i.e. weights above 7.6% cutoff) were PCB 183 (35.6%), PCB 156 (23.4%), PCB 187 (14.7%) and PCB118 (10.9%). There was also a 79% increase in odds for infertility (OR 1.79; 95% CI 1.03, 3.11) (Fig. 5 and Table 3). Chemicals of concern were PCB 187 (42.2%), PCB 156 (22.1%), PCB 74 (9.8%), PCB 183 (9.2%), and transnonachlor (8.0%). None of the WQS-indexes were significantly associated with TTP nor infertility in the three other groups of women.

Fig. 5
figure 5

Association between POP mixture and fertility outcomes in women ≥29 years old not using COCs. Graph shows shape and direction of association between weighted quantile sum (WQS) index and (a) time-to-pregnancy and (b) infertility. Logistic regression fit was for women with parity = 1, non-smoker, regular menses, and BMI = 24 kg/m2

Table 3 Weighted quantile sum index stratified according to age and use of combined oral contraceptives (COCs). Data presented are estimate (SE) and OR (95% CI)

Discussion

The main finding of our study showed that serum concentrations of certain POPs were significantly associated with lower fecundability and higher odds of infertility in women aged ≥29 years who did not use COC as their most recent pre-pregnancy contraceptive even though none of the POPs were associated with fecundability prior to stratification. Our findings suggest that pre-pregnancy use of combined oral contraceptives and age may modify the link between POPs and fecundability. To our knowledge, this is the first study to investigate the chemicals both individually and as a mixture as well as focus on POPs exposure and the effect modification by pre-pregnancy COC use and age in relation to fecundability.

In the individual chemical stratified analyses, higher concentrations of PCBs 153, 180, 183 and 187 were associated either with lower fecundability and/or higher odds for infertility in at least one stratum. However, this monotonic dose response was not observed with HCB, p,p’-DDE, PCBs 99, 138, 156, and 183, where only either the second or third quartiles were significant, but not the fourth quartile. Nonetheless, these chemicals were still associated with reduced fecundability, albeit challenging conventional toxicological dogmas such as monotonicity. There was a positive correlation among the chemicals analyzed, where higher exposure to one POP was linked to increased exposure to other POPs. This demonstrates the need for taking mixture exposures into account since we are exposed to several compounds at the same time. In the current study, we used the WQS regression in determining the chemicals of concern in mixtures that were associated with longer TTP and infertility, which was significant among women aged ≥29 years not using COCs. In this stratum, while the chemicals in the individual analyses were in part similar to those chemicals of concern in the mixture (i.e. PCB 156, 183 and 187 for both TTPand infertility), there were still chemicals that were significant only for the individual analyses (i.e. PCB 180 for TTP; PCB 99, 153, and 180 for infertility) but not for the mixtures and vice versa (i.e. PCB 118 for TTP; PCB 74 and transnonachlor for infertility). Hence, it is important to study chemicals as a mixture because possible antagonistic, synergistic or additive properties between chemicals may be overlooked when chemicals are assessed individually, which may lead to underestimation of hazards and risks [39]. As of yet, there is no single overarching systematic, comprehensive, and integrated approach in understanding mixtures but there is continuous effort globally to develop strategies in studying mixtures. To quantify the effects of mixtures, several statistical approaches have been proposed from simple linear regression to advanced machine learning techniques but there is no one particular method that outperforms the other [40]. This study supports the need to design and optimize approaches in addressing mixtures in epidemiological studies.

Comparison of our results to previous findings is challenging. We have found that most previous studies only utilized pre-pregnancy use of COCs to adjust their model but did not examine it further through stratification. Nonetheless, these studies have shown associations of POPs to longer TTP [23, 24] while other studies failed to replicate these results [20, 22]. Similar to our results, it has also been shown that PCB mixtures were associated with longer TTP [41]. In our study, stratification allowed studying the association of the chemicals in four homogenous groups taking into account two important factors: age and COC use. While this approach prevented direct comparison between groups, we were able to draw specific inferences on associations of chemicals on each stratum.

Previous studies show inconclusive evidence on return to fertility after cessation of COCs. While it was seen that TTP in former COC users were comparable with other contraceptive methods [13], another study showed reduction in TTP after COC use [14]. On the other hand, duration of use may play a role in this association as it was also shown that short-term use of COC may cause delay in return-to-fertility or longer TTP while long-term use of COC has no deleterious effects on fecundability [15]. However, other studies have shown that external factors such as contraceptive use and lifestyle were not as important predictors of TTP as biological factors such as age of conception and menstrual cycle length [42]. Nonetheless, our findings show that both COC use and age play a role on TTP and its relationship with POPs.

The proportions of chemicals differed between SELMA and NHANES, showing that exposure patterns may vary between countries. While both cohorts were of similar reproductive age and time period, other factors such as occupation, race and lifestyle factors may contribute to these differences. This further emphasizes the need to focus further studies on mixture approach in analyzing chemicals exposure. Importantly, POPs are not the only chemicals women of reproductive age are exposed to. Recent biomonitoring studies consistently find more than 300 environmental chemicals or their metabolites in the serum of nearly all analyzed individuals [43]. It is of high priority to study how extensive mixtures of chemicals collectively modulate reproductive functions in women.

Possible biological mechanism

One of the molecular mechanisms that possibly could explain our results is the well-known impact of COCs on adenosine triphosphate (ATP)-binding cassette (ABC) transporters [44]. ABC transporters are highly conserved transmembrane proteins that are responsible for the efflux of substances towards the extracellular space. With this efflux, substances are eliminated and prevented from accumulating inside the cell [45]. While it aids in many physiological processes such as steroid production and immunological processes, it also provides support by efflux of environmental toxins out of the cells [46]. A number of ABC proteins, such as P-glycoprotein encoded by the ABCB1 gene, were shown to be upregulated by estrogen and progesterone (i.e. hormones that are present in COCs) [44, 47]. With the increased production of the ABC transporters, environmental chemicals such as POPs could possibly be eliminated from the cell [46], suggesting to less toxic effects on reproductive functions. It has also been shown that there is an age-related downregulation of ABC transporters [48], which possibly explains observing the association mostly in women older than 29 years old. Further studies are warranted to confirm this hypothesis.

Limitations of the study

Our study has limitations. The chemical concentrations were not lipid adjusted. In addition, our study focused on combined oral contraceptives, excluding other hormonal contraception methods such as progesterone only pills and intrauterine device with hormones. There was neither data on the duration and type of the COC use prior to discontinuation. Moreover, while we recognize that paternal factors may also affect TTP, the only paternal data in the cohort was BMI, which was eventually not included in the model because it did not modify the FR more than 10%. Unplanned pregnancies were also not available in the cohort. Lack of these data may lead to overestimation of findings, hence, further studies are needed to confirm these results. Because this study utilized a pregnancy cohort, the least fertile couples were less represented. In addition, different statistical analyses were performed when studying the chemicals as individual compounds and as a mixture. While Cox proportional hazard was not a problem for possible presence of non-monotonic dose response, it cannot handle highly correlated predictors and cannot be used for mixture analysis. On the other hand, WQS was appropriate for handling intracollinearity but limited to monotonic dose response. Hence, comparison of results between the analyses must be taken with caution.

Conclusions

Although the studied 13 chemicals were banned several decades ago, many of them were still detected in majority of the serum samples collected in 818 pregnant women during the period 2007–2010 living in Sweden. Even though none of the POPs were associated with fecundability prior to stratification, our results indicate that higher POP exposure, both as individual chemicals and as a mixture, may be associated with lower fecundability and increased infertility in women ≥29 years old not using COC as their most recent pre-pregnancy contraception method. This is important as more and more women postpone childbearing, which could make the chemical effects more profound. This observation could also help explain some of the contradictory findings in the literature. Use of COCs may modify the relationship between POPs and TTP. This study also supports studying mixtures in epidemiological studies to further understand combination effects of chemicals.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

Abbreviations

ABC:

ATP-binding cassette

ATP:

Adenosine triphosphate

BMI:

Body mass index

CI:

Confidence interval

COCs:

Combined oral contraceptives

FR:

Fecundability ratio

GM:

Geometric mean

HCB:

Hexachlorobenzene

HCH:

Hexachlorocyclohexane

GC-MS/MS:

Gas chromatography – tandem mass spectrometry

LC-MS/MS:

Liquid chromatography/ triple quadrupole mass spectrometry

LOQ:

Limit of quantification

NHANES:

National Health and Nutrition Examination Survey

OCP:

Organochlorine pesticide

OR:

Odds ratio

PBDE:

Polybrominated diphenyl ether

PCB:

Polychlorinated biphenyl

PeCB:

Pentachlorobenzene

POPs:

Persistent organic pollutants

p,p´-DDT:

Dichlorodiphenyltrichloroethane

p,p´-DDE:

Dichlorodiphenyldichloroethylene

SELMA:

Swedish Environmental Longitudinal, Mother and child, Asthma and allergy

SD:

Standard deviation

SE:

Standard error

TTP:

Time-to-pregnancy

References

  1. Bergman Å, Heindel J, Jobling S, Kidd K, Zoeller RT. State of the science of endocrine disrupting chemicals. United Nations Environment Programme and the World Health Organization; 2013.

  2. Kelley AS, Banker M, Goodrich JM, Dolinoy DC, Burant C, Domino SE, et al. Early pregnancy exposure to endocrine disrupting chemical mixtures are associated with inflammatory changes in maternal and neonatal circulation. Sci Rep. 2019;9(1):1–14.

    Article  CAS  Google Scholar 

  3. Tanner EM, Hallerbäck MU, Wikström S, Lindh C, Kiviranta H, Gennings C, et al. Early prenatal exposure to suspected endocrine disruptor mixtures is associated with lower IQ at age seven. Environ Int. 2020;134.

  4. Czarnota J, Gennings C, Wheeler DC. Assessment of weighted quantile sum regression for modeling chemical mixtures and cancer risk. Cancer Inform [Internet]. 2015;14(Suppl 2):159–171. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=4431483&tool=pmcentrez&rendertype=abstract.

  5. Secretariat of the Stockholm Convention. Stockholm Convention Overview [Internet]. 2008 [cited 2019 Feb 13]. Available from: http://www.pops.int/TheConvention/Overview/tabid/3351/Default.aspx.

  6. Gore AC, Chappell VA, Fenton SE, Flaws JA, Nadal A, Prins GS, et al. EDC-2: the Endocrine Society’s second scientific statement on endocrine-disrupting chemicals. Endocr Rev. 2015;36:E1–150.

    Article  CAS  Google Scholar 

  7. Infertility definitions and terminology [Internet]. WHO. World Health Organization; 2016. [cited 2017 Jul 21]. Available from: http://www.who.int/reproductivehealth/topics/infertility/definitions/en/.

  8. Baird DD, Wilcox AJ, Weinberg CR. Use of time to pregnancy to study environmental exposures. Am J Epidemiol. 1986;124(3):470–80.

    Article  CAS  Google Scholar 

  9. Habbema JD., Collins J, Evers JLH, Lunenfeld B, te Velde ER. Towards less confusing terminology in reproductive medicine:a proposal. Hum Reprod [Internet]. 2004;19(7):1497–501. [cited 2019 Aug 16]. Available from: https://watermark.silverchair.com/deh303.pdf?token=AQECAHi208BE49Ooan9kkhW_Ercy7Dm3ZL_9Cf3qfKAc485ysgAAAlkwggJVBgkqhkiG9w0BBwagggJGMIICQgIBADCCAjsGCSqGSIb3DQEHATAeBglghkgBZQMEAS4wEQQMgpze7PpsklTYaJ_LAgEQgIICDNFp3lU6GKuQ5tbJmn68x0DXnQNW0T3-WMdfh5dI64-ilei.

  10. Joffe M. Time to pregnancy: a measure of reproductive function in either sex. Asclepios project. Occup Environ Med [Internet]. 1997;54(5):289–95. Available from: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC1128774/.

  11. Gesink Law DC, Maclehose RF, Longnecker MP. Obesity and Time to Pregnancy. Hum Reprod [Internet]. 2007;22(2):414–20. [cited 2019 Apr 9]. Available from: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC1924918/pdf/nihms18978.pdf.

  12. The American College of Obstetricians and Gynecologists Committee on Gynecologic Practice, The Practice Committee of the American Society for Reproductive Medicine. Female age-related fertility decline. Fertil Steril [Internet]. 2014;101:633–4. [cited 2019 Apr 9]. Available from: http://fertstertforum.com/asrmpraccom-age-fertility-decline/.

  13. Barnhart KT, Schreiber CA. Return to fertility following discontinuation of oral contraceptives. Fertil Steril [Internet]. 2009;91(3):659–63. [cited 2019 Feb 13]. Available from: https://www.fertstert.org/article/S0015-0282(09)00049-1/pdf.

  14. Hassan MAM, Killick SR. Is previous use of hormonal contraception associated with a detrimental effect on subsequent fecundity? Hum Reprod [Internet]. 2004;19(2):344–51. [cited 2019 Feb 13]. Available from: https://0-academic-oup-com.brum.beds.ac.uk/humrep/article-abstract/19/2/344/2356296.

  15. Mikkelsen EM, Riis AH, Wise LA, Hatch EE, Rothman KJ, Toft Sørensen H. Pre-gravid oral contraceptive use and time to pregnancy: a Danish prospective cohort study. Hum Reprod [Internet]. 2013;28(5):1398–405. [cited 2018 Mar 26]. Available from: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC3627337/pdf/det023.pdf.

  16. Younglai E V, Foster WG, Hughes EG, Trim K, Jarrell JF. Levels of Environmental Contaminants in Human Follicular Fluid, Serum, and Seminal Plasma of Couples Undergoing In Vitro Fertilization. Arch Environ Contam Toxicol [Internet]. 2002;43:121–6. [cited 2017 Jul 19]. Available from: https://0-link-springer-com.brum.beds.ac.uk/content/pdf/10.1007%2Fs00244-001-0048-8.pdf.

  17. Björvang RD, Damdimopoulou P. Persistent environmental endocrine-disrupting chemicals in ovarian follicular fluid and in vitro fertilization treatment outcome in women Persistent environmental endocrine-disrupting chemicals in ovarian follicular fluid and in vitro fertilization treatm. Ups J Med Sci [Internet]. 2020;1–10. [cited 2020 Feb 25]. Available from: https://0-www-tandfonline-com.brum.beds.ac.uk/action/journalInformation?journalCode=iups20.

  18. Golub MS, Donald JM, Reyes JA. Reproductive Toxicity of Commercial PCB Mixtures: LOAELs and NOAELs from Animal Studies. Environ Health Perspect [Internet]. 1991;94:245–53. [cited 2019 Feb 13]. Available from: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC1567963/pdf/envhper00414-0235.pdf.

  19. Peterson RE, Michael Theobald H, Kimmel GL, Kimme GL. Developmental and Reproductive Toxicity of Dioxins and Related Compounds: Cross-Species Comparisons Developmental and Reproductive Toxicity of Dioxins and Related Compounds: Cross-Species Comparisons. Crit Rev Toxicol [Internet]. 1993;23(3):283–335. [cited 2019 Feb 13]. Available from: https://0-www-tandfonline-com.brum.beds.ac.uk/action/journalInformation?journalCode=itxc20.

  20. Harley KG, Marks AR, Chevrier J, Bradman A, Sjödin A, Eskenazi B. PBDE concentrations in Women’s serum and Fecundability. Environ Health Perspect. 2010;118(5):699–704.

    Article  CAS  Google Scholar 

  21. Harley KG, Marks AR, Bradman A, Barr DB, Eskenazi B. DDT exposure, work in agriculture, and time to pregnancy among farmworkers in California. J Occup Environ Med [Internet]. 2008;50(12):1335–1342. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2684791&tool=pmcentrez&rendertype=abstract.

  22. Cohn BA, Cirillo PM, Sholtz RI, Ferrara A, Park J, Schwingl PJ. Polychlorinated biphenyl (PCB) exposure in mothers and time to pregnancy in daughters. Reprod Toxicol. 2011;31(3):290–6.

    Article  CAS  Google Scholar 

  23. Gesink Law DC, Klebanoff MA, Brock JW, Dunson DB, Longnecker MP. Maternal serum levels of polychlorinated biphenyls and 1,1-dichloro-2,2- bis(p-chlorophenyl) ethylene (DDE) and time to pregnancy. Am J Epidemiol. 2005;162(6):523–32.

    Article  Google Scholar 

  24. Axmon A, Thulstrup AM, Rignell-Hydbom A, Pedersen HS, Zvyezday V, Ludwicki JK, et al. Time to pregnancy as a function of male and female serum concentrations of 2,2′4,4′5,5′-hexachlorobiphenyl (CB-153) and 1,1-dichloro-2,2-bis (p-chlorophenyl)-ethylene (p,p′-DDE). Hum Reprod. 2006;21(3):657–65.

    Article  CAS  Google Scholar 

  25. Bornehag CG, Moniruzzaman S, Larsson M, Lindström CB, Hasselgren M, Bodin A, et al. The SELMA study: A birth cohort study in sweden following more than 2000 mother-child pairs. Paediatr Perinat Epidemiol; 2012.

  26. Koponen J, Rantakokko P, Airaksinen R, Kiviranta H. Determination of selected perfluorinated alkyl acids and persistent organic pollutants from a small volume human serum sample relevant for epidemiological studies. J Chromatogr A [Internet]. 2013;1309:48–55. [cited 2017 Aug 1]. Available from: http://ac.els-cdn.com.proxy.kib.ki.se/S0021967313011242/1-s2.0-S0021967313011242-main.pdf?_tid=f6bdbf20-76c0-11e7-ba9e-00000aab0f01&acdnat=1501595861_6bae0125da14a452c7d3799313db2294.

  27. Shu H, Lindh CH, Wikströ S, Bornehag C-G. Temporal trends and predictors of perfluoroalkyl substances serum levels in Swedish pregnant women in the SELMA study. PLoS One [Internet]. 2018;13(12). [cited 2019 Feb 13]. Available from: https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0209255.

  28. Centers for Disease Control and Prevention (CDC), National Center for Health Statistics (NCHS). National Health and Nutrition Examination Survey Data [Internet]. Hyattsville, MD; 2008. Available from: https://wwwn.cdc.gov/nchs/nhanes/search/datapage.aspx?Component=Laboratory&CycleBeginYear=2007.

  29. Centers for Disease Control and Prevention (CDC), National Center for Health Statistics (NCHS). National Health and Nutrition Examination Survey Data [Internet]. Hyattsville, MD; 2010. Available from: https://wwwn.cdc.gov/nchs/nhanes/search/datapage.aspx?Component=Laboratory&CycleBeginYear=2007.

  30. Bello B, Heederik D, Kielkowski D, Wilson K. Increased time-to-pregnancy is associated with domestic work in South Africa. Reprod Health [Internet]. Reproductive Health. 2016;13(1):1–9. Available from: https://0-doi-org.brum.beds.ac.uk/10.1186/s12978-016-0224-y.

  31. Szumilas M. Explaining Odds Ratios. J Can Acad Child Adolesc Psychiatry [Internet]. 2010;19(3):227–9. [cited 2020 Feb 2]. Available from: http://www.csm-oxford.org.uk/.

  32. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of weighted Quantile sum regression for highly correlated data in a risk analysis setting. J Agric Biol Environ Stat. 2015;20(1):100–20.

    Article  Google Scholar 

  33. R Core Team. R: a language and environment for statistical computing [internet]. Vienna, Austria: R Foundation for statistical Computing; 2018. Available from: https://www.r-project.org/.

  34. RStudio Team. RStudio: integrated development environment for R [internet]. Boston, MA: RStudio, Inc.; 2016. Available from: http://www.rstudio.com/.

    Google Scholar 

  35. Renzetti S, Curtin P, Just AC, Bello G, Gennings C. gWQS: generalized weighted Quantile sum regression [internet]. 2018. Available from: https://cran.r-project.org/package=gWQS.

  36. Wickham H. ggplot2: elegant graphics for data analysis [internet]. New York: Springer-Verlag; 2009. Available from: http://ggplot2.org.

    Book  Google Scholar 

  37. Kassambara A. ggpubr: “ggplot2” Based Publication Ready Plots. R package version 0.1.6. 2017.

  38. Wickham H. Reshaping data with the reshape package. J Stat Softw [Internet]. 2007;21(12):1–20. Available from: http://www.jstatsoft.org/v21/i12/.

    Google Scholar 

  39. Bornehag CG, Kitraki E, Stamatakis A, Panagiotidou E, Rudén C, Shu H, et al. A novel approach to chemical mixture risk assessment—linking data from population-based epidemiology and experimental animal tests. Risk Anal. 2019;39(10):2259–71.

    Article  Google Scholar 

  40. Taylor KW, Joubert BR, Braun JM, Dilworth C, Gennings C, Hauser R, et al. Statistical approaches for assessing health effects of environmental chemical mixtures in epidemiology: lessons from an innovative workshop. Env Heal Perspect. 2016;124(12):A227–9.

    Article  Google Scholar 

  41. Gennings C, Carrico C, Factor-Litvak P, Krigbaum N, Cirillo PM, Cohn BA. A cohort study evaluation of maternal PCB exposure related to time to pregnancy in daughters. Environ Health [Internet]. 2013;12(1):66. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/23962309%5Cn, http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC3766643.

  42. Axmon A, Rylander L, Albin M, Hagmar L. Factors affecting time to pregnancy. Hum Reprod. 2006;21(5):1279–84.

    Article  CAS  Google Scholar 

  43. U.S. Department of Health and Human Services. Fourth National Report on Human Exposure to Environmental Chemicals [Internet]. 2018. Available from: https://www.cdc.gov/exposurereport/pdf/FourthReport_UpdatedTables_Volume1_Mar2018.pdf.

  44. Bloise E, Ortiga-Carvalho TM, Reis FM, Lye SJ, Gibb W, Matthews SG. ATP-binding cassette transporters in reproduction: a new frontier. Hum Reprod Update. 2016;22(2):164–181.

  45. Tarling EJ, De Aguiar Vallim TQ, Edwards PA. Role of ABC transporters in lipid transport and human disease. Trends Endocrinol Metab [Internet]. 2013;24(7):342–50. [cited 2019 Feb 13]. Available from: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC3659191/pdf/nihms-440487.pdf.

  46. Mazur CS, Marchitti SA, Dimova M, Kenneke JF, Lumen A, Fisher J. Human and Rat ABC Transporter Efflux of Bisphenol A and Bisphenol A Glucuronide: Interspecies Comparison and Implications for Pharmacokinetic Assessment. Toxicol Sci [Internet]. 2012;128(2):317–25. [cited 2019 Feb 13]. Available from: https://0-academic-oup-com.brum.beds.ac.uk/toxsci/article-abstract/128/2/317/1649887.

  47. Peng R, Zhang H, Zhang Y, Wei D-Y. Effects of the ABCB1 (1199G &gt; A) Polymorphism on Steroid Sex Hormone-Induced P-Glycoprotein Expression, ATPase Activity, and Hormone Efflux. Med Sci. 2015;3:124–37.

  48. Erdo F, Krajcsi P. Age-Related Functional and Expressional Changes in Efflux Pathways at the Blood-Brain Barrier. Front Ageing Neurosci [Internet]. 2019;11:1–16. [cited 2020 Feb 2]. Available from: https://www.frontiersin.org/articles/10.3389/fnagi.2019.00196/full.

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the Swedish Research Council FORMAS [project EDC2020], and Jane & Aatos Erkko Foundation [2015 “Keeping eggs healthy”]. The financial supporters had no role in the study design, collection and analysis of data, data interpretation or in writing the report. Open access funding provided by Karlstad University.

Author information

Authors and Affiliations

Authors

Contributions

RDB wrote the paper, did the statistical analyses and data interpretation, and prepared the figures and Tables. CG, PIL and GH assisted in the statistical analyses and data interpretation. HK, PR, PR, and CL were responsible for the chemical analyses. PD was involved in interpreting data, and writing the paper. CGB assisted in the statistical analyses, interpretation of data, writing of the paper, and was responsible for the study design and PI of the SELMA study. All authors approved the last version of the manuscript.

Corresponding author

Correspondence to Carl-Gustaf Bornehag.

Ethics declarations

Ethics approval and consent to participate

The SELMA study has been approved by the Uppsala regional ethical committee (Dnr. 2007/062, 2015/177)). Participants provided their written consent in accordance with the Helsinki Declaration.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Estimated fecundability ratio (FR) for time-to-pregnancy (TTP) and odds ratio (OR) for infertility prior to stratification. Data are presented as FR (95% CI) or OR (95% CI). Supp Figure S1. Flowchart of recruitment process and analyses in the SELMA study. Abbreviations: COC – combined oral contraceptives; FR – fecundability ratio; OR – odds ratio; POPs – persistent organic pollutants; TTP – time-to-pregnancy; WQS – weighted quantile sum. Supp Figure S2. Sensitivity analyses with censoring at 10 months. Dot shows fecundability ratio (FR) with 95% confidence interval error bars of chemicals as continuous (log) and quartiles (Q1 as reference). FR > 1 denotes higher fecundability and shorter TTP while FR < 1 denotes lower fecundability and longer TTP. FR indicated with * for p value < 0.05. P trend shown for significant linear trend of FR. Abbreviations: COC – combined oral contraceptives; HCB - hexachlorobenzene; p,p´-DDE – dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyls; Q1 – first quartile, Q2 – second quartile, Q3 – third quartile, Q4 – fourth quartile. Supp Figure S3. Sensitivity analyses with censoring at 14 months. Dot shows fecundability ratio (FR) with 95% confidence interval error bars of chemicals as continuous (log) and quartiles (Q1 as reference). FR > 1 denotes higher fecundability and shorter TTP while FR < 1 denotes lower fecundability and longer TTP. FR indicated with * for p value < 0.05. P trend shown for significant linear trend of FR. Abbreviations: COC – combined oral contraceptives; HCB - hexachlorobenzene; p,p´-DDE – dichlorodiphenyldichloroethylene; PCB – polychlorinated biphenyls; Q1 – first quartile, Q2 – second quartile, Q3 – third quartile, Q4 – fourth quartile.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Björvang, R.D., Gennings, C., Lin, PI. et al. Persistent organic pollutants, pre-pregnancy use of combined oral contraceptives, age, and time-to-pregnancy in the SELMA cohort. Environ Health 19, 67 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-020-00608-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-020-00608-8

Keywords