- Open Access
A Cohort study evaluation of maternal PCB exposure related to time to pregnancy in daughters
Environmental Health volume 12, Article number: 66 (2013)
Polychlorinated biphenyls (PCBs) remain ubiquitous environmental contaminants. Developmental exposures are suspected to impact reproduction. Analysis of mixtures of PCBs may be problematic as components have a complex correlation structure, and along with limited sample sizes, standard regression strategies are problematic. We compared the results of a novel, empirical method to those based on categorization of PCB compounds by (1) hypothesized biological activity previously proposed and widely applied, and (2) degree of ortho- substitution (mono, di, tri), in a study of the relation of maternal serum PCBs and daughter’s time to pregnancy.
We measured PCBs in maternal serum samples collected in the early postpartum in 289 daughters in the Child Health and Development Studies birth cohort. We queried time to pregnancy in these daughters 28–31 years later. We applied a novel weighted quantile sum approach to find the bad-actor compounds in the PCB mixture found in maternal serum. The approach includes empirical estimation of the weights through a bootstrap step which accounts for the variation in the estimated weights.
Bootstrap analyses indicated the dominant functionality groups associated with longer TTP were the dioxin-like, anti-estrogenic group (average weight, 22%) and PCBs not previously classified by biological activity (54%). In contrast, the unclassified PCBs were not important in the association with shorter TTP, where the anti-estrogenic groups and the PB-inducers group played a more important role (60% and 23%, respectively). The highly chlorinated PCBs (average weight, 89%) were mostly associated with longer TTP; in contrast, the degree of chlorination was less discriminating for shorter TTP. Finally, PCB 56 was associated with the strongest relationship with TTP with a weight of 47%.
Our empirical approach found some associations previously identified by two classification schemes, but also identified other bad actors. This empirical method can generate hypotheses about mixture effects and mechanisms and overcomes some of the limitations of standard regression techniques.
Despite their ban in the 1970s, polychlorinated biphenyls (PCBs) remain ubiquitous environmental contaminants due to their chemical stability and lipophilic properties . In fact, one or more PCBs are detectable in over 80% of blood samples in contemporary participants in the National Health and Nutrition Examination Survey (NHANES; e.g.,  studying 2001–02 NHANES data for women of child bearing age). PCBs are hormonally active, and there is mounting evidence that in utero PCB exposure shortens length of gestation in humans  and is inversely associated with fetal growth . Our focus herein from a prospective study, prenatal exposure to PCBs may disrupt pregnancy in humans .
Fecundability (probability of pregnancy in each cycle) has been used as a sensitive marker to identify reproductive hazards in humans . Increased time to pregnancy (TTP) may indicate a problem at one or more of several stages of human reproduction, including gametogenesis, transport of gametes in both male and female reproductive tracts, fertilization, migration of the zygote to the uterus, implantation and early survival of the conceptus [4, 5]. Increased TTP is most likely multifactorial with risk factors in utero as demonstrated by effects of prenatal diethylstilbestrol exposure , childhood and puberty as demonstrated by effects of treatment for childhood cancer , and adult life where smoking has been extensively investigated .
Cohn et al.  tested the hypothesis that in utero exposure to PCBs alters TTP in humans using data from the Child Health Development Studies (CHDS). Their data included measures of maternal exposures in blood samples collected near delivery (days 1–3 after daughter’s birth during 1960–63) in relation to her daughter’s TTP decades later. Their analysis was based on a previous categorization of PCB congeners by Wolff et al.  based on functionality and inclusion in CHDS maternal serum samples:
● Group 1: potentially estrogenic and persistent: PCB 101, 187, 201;
● Group 2: potentially anti-estrogenic, immunotoxic, dioxin-like
o Group 2A: non-ortho or mono-ortho in structure: PCB 66, 74, 105, 118, 156, 167
o Group 2B: di-ortho and limited dioxin-like activity: PCB 138, 170
● Group 3: phenobarbital, CYP1A and CYP2B inducers: PCB 99, 153, 180, 203, 183.
Serum concentrations for two uncategorized PCBs (PCB56 and PCB146) were also included in the analysis herein.
The model building strategy conducted by Cohn et al.  included two steps. Step 1 entered congeners in a single model with potentially opposing effects (e.g., estrogenic, Group1, and anti-estrogenic, Group 2) to avoid “overlooking any associations for estrogenic compounds that might be confounded by the presence of correlated anti-estrogenic compounds.” The second step included remaining congeners from Group 3. Non-significant congeners were removed from the model when their p values exceeded 0.10. The initial model was adjusted for whether the daughter was breastfed (yes or no) and race (African American vs all other); the final model was further adjusted for the following maternal variables: age, body mass index, and lipids (triglycerides, cholesterol) to determine whether these altered associations observed. Following this strategy, PCB congeners 187, 156, and 99 in mother’s serum were associated with longer TTP in their daughters; PCB congeners 105, 138 and 183 were associated with shorter TTP.
Using the same data, our objectives are to use a holistic empirical approach to investigate the association between PCB functionality and TTP in humans. Recognizing the complex correlation structure among the PCB congeners (Additional file 1: Table S1) and with a limited sample size, standard regression strategies are problematic (e.g., ). Instead, we use a nonlinear weighted quartile sum (WQS) approach  nested within a parametric proportional hazards model. The method is nonlinear due to the empirical estimation of a weighted sum of quartiles with the typical regression coefficients. The index is based on quartiles of mixture components to reduce the impact of extreme values in typically right-skewed concentration distributions. Bootstrap distributions of the estimated weights are used to identify the strength of association between the Wolff et al. categories and TTP, the degree of chlorine substitution and TTP, and to define an ‘average weighted quartile sum index’ that is associated with TTP. The weighted index stabilizes the ill conditioning in the analysis due to multicollinearity .
Here we contribute to the literature in two significant fashions. First, we consider maternal prenatal exposures and TTP in offspring from a well- documented birth cohort. Previous work (e.g. ) has considered exposures and TTP only in the parental generation. Second, we analyze PCBs as mixtures using a novel empirical method.
Description of study sample
As described in Cohn et al. , the study sample consists of women born between 1960 and 1963 in the Oakland, California area, whose mothers were enrolled in the Kaiser Permanente Health Plan. These women and their mothers were among participants in the Child Health and Development Studies (CHDS), an investigation of prenatal determinants of infant, child, adolescent and adult health. Subjects were respondents to a follow-up study conducted in 1990–1991 which investigated prenatal determinants of TTP. Women were ages 28–31 years (65% participation rate after 30 years since their birth) when they completed a questionnaire on the number of non-contracepting menstrual cycles that preceded their most recent pregnancy; or the number of non-contracepting cycles during which they were at risk for pregnancy even when pregnancy was not achieved. Women who did not plan their pregnancies as well as women who did plan their pregnancies were included in the study as long as they could report their TTP, reducing planning bias. Details of the study sample are given below.
CHDS daughters born during 1960–1963 who resided in the San Francisco Bay area through ages 15–17 (N = 1,003) were eligible for this study if their mothers had reported information on smoking during pregnancy (N = 991). TTP data was originally collected for a study of prenatal tobacco exposure conducted by Dr. Pam Schwingl for her dissertation research. Of these 991 eligible women, 345 did not participate for the following reasons: no address, N = 117, no response to mail, N = 197, deceased or institutionalized, N = 7, refused, N = 24. There were 646 participants in the study (65% of 991 eligible), of which 357 were excluded from analysis of TTP for the following reasons: respondent was never at risk for pregnancy (N = 163), respondent could not estimate TTP due to sporadic birth control use or birth control failure (N = 119), incomplete questionnaire (N = 47), insufficient prenatal serum for assays (N = 18), younger sister (N = 5), under age 16 at most recent pregnancy “attempt” (N = 5).
Women in the following categories were included in the present study whether or not they had planned their pregnancy, as long as they could report TTP: women who omitted birth control just once, women who became pregnant during an interval between birth control methods or after stopping birth control, women who became pregnant when breast feeding, women who never used birth control, women who were attempting pregnancy after a prior pregnancy, and women who never became pregnant, but had been at risk for pregnancy for a known interval, whether or not they desired to become pregnant. Women were asked to report their TTP for their most recent pregnancy in order to minimize recall errors. Among the 289 women included in this report, 43% of women reported TTP for their first pregnancy, 26% for their second pregnancy, 15% for their third pregnancy, 8% for their fourth pregnancy, and 6% for their fifth or more, with 2% where gravidity before attempt was unknown. Women reporting on time to first pregnancy include women who were at risk for pregnancy, but failed to conceive.
In human observational studies where being at risk for pregnancy is not an experimental condition, there are expected differences between women who can report TTP and women who cannot. The two primary groups of women who are deleted from the TTP analysis are women never at risk for pregnancy (N = 163) and women who became pregnant and could not answer about TTP (N = 119). Compared to the analysis sample, women never at risk for pregnancy had higher personal incomes, were much more likely to have never been married, and were more highly educated and more likely to be Asian. Compared to the analysis sample, women who used birth control sporadically and could not report their TTP were more likely to be African American or Asian, slightly more likely to be never married, somewhat better educated, but were similar to the analysis sample for personal income. These differences result from personal, social and economic factors related to reproductive behavior, and are not surprising. Our results apply to those women who were at risk for pregnancy and could provide a TTP, including women who planned their pregnancies and women who did not plan their pregnancies. This research was approved by the Institutional Review Board of the Public Health Institute, and study subjects gave informed consent.
In utero exposure to PCBs among CHDS daughters was estimated by assays of serum samples drawn during the early postpartum period (days 1–3) from their mothers. PCBs were assayed in the laboratory of Dr. Mary Wolff. Details are provided in Cohn et al. .
Time to pregnancy was modeled using a two-parameter Weibull baseline hazard model, allowing for the “hazard” of pregnancy to be monotonic (increasing, decreasing, or constant) in a Weibull proportional hazards model. The model was parameterized using a nonlinear weighted quartile score regression strategy (described below). Following Cohn et al. , the initial model included two covariates: race (African American vs all other) and whether the daughter was breastfed (yes or no). The final model was further adjusted for maternal variables: age, body mass index, and maternal lipids (triglycerides, cholesterol) to determine whether these altered associations observed. The assumption of a Weibull hazard model with proportionality across covariate categories was verified using diagnostic plots of log(−log(S(t)) vs log(t), where S(t) is the survival function.
Construction of weighted quartile score Weibull proportional hazards model:
A quartile score (robust to extreme values) was calculated for each chemical: values in the 1st quartile, q = 0; values in the 2nd quartile, q = 1; values in the 3rd quartile, q = 2; values in the 4th quartile, q = 3. The body burden from the set of PCBs was calculated by summing the quartile scores across chemicals. However, the potency of the chemicals vary substantially; the inclusion of less potent chemicals may add “noise” in the score and mask the signal from more potent chemicals, and this problem may be exacerbated when the less potent chemicals are more prevalent than more toxic chemicals. Therefore, we used weighted quartile score regression within a Weibull proportional hazards model. The general approach was characterized by Carrico ; who extended the work of Gennings et al.  and Christensen et al.  with the addition of a bootstrap step for determining the weights.
Define t as the time (i.e., cycle) of pregnancy or last cycle of not pregnant. For each daughter, data are available for serum PCB congener concentrations from her mother within 3 days of her delivery along with covariates. A nonlinear Weibull proportional hazards model was parameterized as:
where t is the cycle number; γ and λ are Weibull parameters for shape and scale, respectively; qj is the quartile score for the jth chemical; the w j ∈ [0, 1] are estimated using a nonlinear programming optimization algorithm subject to the constraint that ∑ j w j = 1. The z’s represent covariates (e.g., race, breast fed indicator) with the corresponding unknown parameters estimated using maximum likelihood estimation. Initial estimates for the covariates were based on estimates from a proportional hazards model without adjusting for PCBs (from Proc LifeReg in SAS version 9.2); the weights were initially set to be equal and constrained to sum to 1.0. The trust region method was used for the primary analysis as it is found to be stable for small to moderate sample sizes . To ensure that the results were not dependent on the choice of the trust region method, the conjugate gradient method, the Newton–Raphson, the Nelder-Mead Simplex, and the Newton–Raphson with ridging as available in Proc NLP in SAS (version 9.2) were also assessed and the results were consistent using all algorithms. Following the approach of Carrico , the distributions of the weights were accommodated using 1000 bootstrap samples as follows:
● Bootstrap samples of size 289 (as in the observed data) were randomly generated by sampling with replacement and a set of weights and regression coefficients were estimated for each;
● In some datasets, without constraining the algebraic sign of β 1, the slope estimate for the weighted PCB score was estimated to be <0; in other datasets, the slope estimate for another PCB score was estimated to be >0.
● Histograms for the weights were constructed for the first 1000 bootstrap analyses with positive and significant slope estimates (p < 0.05), and for the first 1000 bootstrap analyses with negative and significant slope estimates.
● Average weighted quartile sum indices were constructed based on the mean weights from 1000 bootstrap samples with positive and significant slope estimates and from 1000 bootstrap samples with negative and significant slope estimates. These indices were tested for significance in the original dataset.
The weights to be averaged were selected based on an indication of significance in the bootstrap sample such that the weights are interpretable; i.e., when the slope estimate is near zero, the weights are not informative. Conditional on the average bootstrap weights, when β 1 in equation (1) is positive and significant, we claim the weighted quartile score is associated with an increase in the “hazard” of pregnancy, i.e., shorter TTP; when β 1 is negative and significant, we claim the weighted quartile score is associated with a decrease in the hazard of pregnancy, i.e., longer TTP. Both a positive case and a negative case were estimated separately from the observed data. From the model in (1) the relative hazard associated with an increase of 1 unit of the PCB scores was exp(β1).
Sums of estimated average weights were tabulated within the Wolff et al. PCB categorization considered by Cohn et al.  and the degree of chlorination considered by Kezios et al. . These sums may be interpreted as the strength (proportion) of the association due to the PCB category.
Ideally, when the sample size is large enough, the weighted quartile score regression strategy includes randomly splitting the data into a ‘test’ dataset and a ‘validation’ dataset. The test dataset is used to estimate the weights and the significance of slope parameters is based on the validation data where the weighted score is fixed. However, in the current case, the sample size is not considered large enough to (1) provide stable estimation of the weights in the test dataset; and (2) maintain adequate power in the validation dataset to test for the significance of β 1.
Of the 289 women included in the study sample, 239 women reported the number of non-contracepting cycles resulting in pregnancy. Fifty women reported either 12 or 13 non-contracepting cycles from which pregnancy was not achieved – a 17% censoring rate. The estimated median time to pregnancy was 2 cycles with the middle 50% (25th percentile to 75th percentile) between 1 and 6 cycles. The young age of the respondents (ages 27–31 at interview) resulted in short time interval between time of interview and their most recent pregnancy, or most recent at risk period (median interval of 3 years, with 25% reporting on a pregnancy or attempt time of 1 year of less, and 75% within the last 5 years.
Parameter estimates were found for the nonlinear Weibull proportional hazards model parameterized with the weights for the 18 PCBs and two covariates (race and breastfed indicators; equation (1)) using the trust region method. A bootstrap analysis was conducted to account for the variation/distribution of the weights. Both positive and negative slope estimates were evaluated using constraints on the coefficient when necessary. Supplementary figures include histograms for the bootstrap distributions for each PCB in a weighted quartile score associated with a positive effect on hazard of pregnancy (Additional file 1: Figure S1) and with a negative effect on the hazard (Additional file 1: Figure S2). These were used to construct histograms for the sum of weights within each functional category (Figure 1). In the analysis associated with longer TTP,
● the average weight from 1000 bootstrap distributions was 16% from estrogenic PCBs;
● the average weight was 22% from anti-estrogenic and dioxin-like PCBs;
● the average weight was 54% from the ‘other’ category; and
● the average weight for the other two groups were less than 10%.
This is in contrast to the weight distributions in the analysis associated with shorter TTP, a more complex scenario:
● the average weight associated with estrogenic PCBs was 17%;
● the average weight associated with anti-estrogenic PCBs that are dioxin-like was 45%;
● the average weight associated with anti-estrogenic PCBs that have limited dioxin activity was 14%; and
● the average weight for the PB-inducers was 23%.
For the chlorination categorization (Figure 2) in the analysis of longer TTP, the highly chlorinated PCBs were mostly related to longer TTP with average weight of 89%. In contrast, the degree of chlorine substitution was not as discriminating for short TTP with the average weights mostly similar (25%, 26% and 32% for mono, di, and higher) with the exception of lower weights for the tri-ortho group (14%).
A final model was estimated including the two initial covariates, the positive PCB score (using average bootstrap weights, Table 1), the negative PCB score (using average bootstrap weights, Table 1), and the maternal variables of age, BMI, cholesterol and triglycerides. None of the covariates were significant. Both PCB scores had p values less than 0.001; these p values should not be strictly interpreted, as the weighted scores were not a priori defined.
Bootstrap distributions of the relative hazard associated with a positive slope estimate and a negative slope estimate (based on different sets of PCBs) from the final model are provided in Figure 3. On average, the relative hazard of pregnancy associated with longer TTP is 0.70 per unit increase using the weights in Table 1, column 3. On average, the hazard of pregnancy increases 64% per unit increase in the PCB score (i.e., a relative hazard of 1.64) associated with shorter TTP using the weights in Table 1, column 5.
Finally, the cumulative probability of pregnancy was calculated for the four cases defined by the two scores (PCB + and PCB-) cut at their medians (Figure 4). The combination with the shortest TTP (high concentrations of PCB + and low concentrations of PCB-) has roughly a 10-20% increase in the probability of pregnancy compared to the combination with the longest TTP (low concentrations of PCB + and high concentrations of PCB-). In addition, the plot indicates the cumulative probability of pregnancy increases most rapidly roughly within the first three months.
A previous paper  demonstrated that in utero exposure to some PCBs may impact human reproduction, either by increasing or decreasing the time to pregnancy. In the studied population, women were exposed to a mix of PCB congeners before birth that predicted longer and shorter TTP and a reduced probability of contraception 30 years later.
The analysis of these same data using an alternative approach provides further insight into the possible link between different categories of PCB functionality and TTP. Using a weighted quartile score approach within a parametric proportional hazards model, empirical estimates of weights per PCB congener were estimated in a weighted quartile score that was positively associated with the hazard of pregnancy (i.e., shorter TTP) and a score that was negatively associated with the hazard of pregnancy (i.e., longer TTP). Generally, different congeners were identified in the two scores. In particular, using a bootstrap analysis to accommodate the variation in the estimated weights, the dominant functionality groups associated with longer TTP were the dioxin-like, anti-estrogenic group (average, 22%) and the ‘uncategorized’ group (PCBs not characterized by Wolff et al.; 54%). In contrast, the uncategorized group was not important in the association with shorter TTP, and the anti-estrogenic groups and the inducers group played an increased role (60% vs 22%; and 23% vs 9%; respectively). Using the second categorization, the highly chlorinated PCBs were mostly associated with longer TTP (89% of the weight), which was not so for short TTP.
Interpretation of the weights from weighted quartile sum regression was instructed through simulation studies conducted by Carrico  (detailed below). She demonstrated that clusters of high pairwise correlations (as in the PCBs) have individual weights that tend to be smaller compared to components that are less correlated. Thus, the weights should not be interpreted as indicators of relative importance in the association of interest, rather the correlation pattern should be considered. Here, since the pairwise correlations for PCB 56 were generally low (mean = 0.18, standard deviation =0.12; maximum = 0.39) and the weight was considerably larger than any other weight, we conclude there is an interpretable association with longer TTP.
Carrico  conducted extensive simulation studies to characterize WQS regression using environmentally relevant correlation patterns among mixture components. The strength of association between subsets of components and a continuous response variable were varied between correlations of 0.1 and 0.3 using observed correlation patterns from urinary concentrations of environmental chemicals from NHANES data. The method was evaluated in terms of the number of correctly identified bad actors and the number of incorrectly identified components. With a sample size of 250 in 1000 simulated studies where 8 of 11 components were associated with Y with a correlation of 0.3, WQS correctly identified at least 7 of the 8 components 75% of the time with none incorrectly identified more than 95% of the time. When the correlation decreased to 0.1 (a case of poor power), WQS correctly identified 4 or more of the 8 components 80% of the time with none incorrectly identified 35% of the time and in 60% of the cases only selected one component incorrectly. Carrico also conducted simulation studies using LASSO to detect the bad actors and found it tends to identify the correct bad actors but generally identifies incorrect components with high probability.
In our present study we found increased TTP with PCBs thought to be antiestrogenic (PCBs 66, 74, 105, 118, 156 and 167) as well as those not classified in any of the Wolf et al.  groupings (PCBs 56, 146). In contrast, shorter TTP was associated with PCBs thought to be anti-estrogenic (PCBs 66,74,105, 118, 156 and 167) and those thought to be PB-inducers (PCBs 99, 153, 180, 203 and 183). Buck Louis et al.  in the LIFE study also found reduced fecundability with maternal PCB congeners 118 and 167. They also found reduced fecundability with paternal PCB congeners 138, 167 and 170; for the latter two we found shorter TTP. Because the exact mechanisms by which these associations are likely, it is hard to discern why the results differed for paternal exposure patterns. One possible explanation is that we considered empirically derived mixtures of congeners while Buck Louis evaluated congeners individually.
Recent literature suggests that the developing reproductive system may be influenced by maternal exposures and that such outcomes may differ by sex. Jensen, et al.  in a study of Danish twins reported reductions in female fecundity and no relationship with male fecundity related to maternal smoking during pregnancy but not to exposure during childhood. These results support some earlier data ([18, 19]) but not other data ([20, 21]). Nevertheless, the data support animal studies suggesting that male and female mouse fetuses exposed to polycyclic aromatic hydrocarbons, the major toxic component of cigarette smoke, have reduced fertility, although the precise mechanisms of this association are unknown. More recent data followed a cohort of children exposed to PCBs as well as polychlorinated dibenzo-p-dioxins (PCDD/Fs) in utero and assessed reproductive development at age 8 . Girls exposed to the highest levels of PCBs (either individual congeners or total) had lower genital developmental ages and shorter fundi and uteri lengths and both boys and girls exposed to high levels of total toxic equivalents (TEQ; a combination of PCDD/Fs and specific dioxin like PCBs) had reduced concentrations of estradiol. These results are consistent with previous data that find reduced E2 levels among men exposed to TCDD during infancy  and with data that find delayed breast development in girls exposed to PCDD/Fs in utero . Although the precise biological mechanism for these associations is not known, it is thought that these endocrine disrupting chemicals may interfere with the normal action or regulation of reproductive endocrine hormones , with estrogen catabolism  or with alterations in the aryl hydrocarbon receptor signaling system . Delayed development, particularly in female 8-year olds support our results of reductions in fecundity among women exposed in utero to specific PCBs.
We are aware of no other prior studies of the impact of maternal perinatal serum levels of organochlorine compounds on time to pregnancy in daughters decades later, other than our work in the Child Health and Development Studies cohort [4, 28]. Previous work (e.g. ) has considered exposures during pregnancy and TTP only in the same (mother’s) generation. Data from the National Collaborative Perinatal Project find inconclusive evidence regarding maternal concurrent exposure to organochlorines and TTP again in the maternal generation . Results to date suggest that developmental exposures, like those that we report here, are more likely to affect time to pregnancy.
When we published our first paper on in maternal PCBs and daughter’s TTP , there was active work to determine whether the classification proposed by Wolf and colleagues  demonstrated predictive validity for a number of outcomes. By now it appears not to hold up well as we noted in our 2011 TTP paper  and in a recent paper on maternal breast cancer  and in other papers on infant health outcomes in the CHDS . It is for this very reason that we attempted an alternative approach that did not require a priori classification. In reporting our results, we did consider the Wolff et al. classification in order to provide information in this paper that could be compared to other papers based on this widely used classification. We do not mean to suggest that this classification is better than any other, but we are suggesting that an empirical approach, independent of a priori classification, can be very useful in light of a gap in understanding the biological effects of simultaneous exposure to multiple congeners in different proportions. The chlorination grouping demonstrated the benefit of combining a hypothesized categorization with an empirical approach in the finding that the higher chlorinated PCBs are mostly associated with longer TTP and not shorter TTP.
We were unable to stratify analysis by gravidity, as sample sizes would have been too small for analysis by the methods we used here. However, we previously compared results for daughters who were nulligravidas (N = 54), primigravidas (N = 70) and multigravidas (N = 159) in our prior publication that is the basis for this comparative methods study (Table 3 in ). In that comparison we showed, as expected, that TTPs were shorter for parous women, but that overall, PCB associations observed were consistent for all three gravidity groups. Current levels of PCBs were not assessed, but maternal and adult daughter exposures to tobacco (including passive smoking), alcohol and caffeine, daughter’s health history, including sexually transmitted diseases, race, education, income (mother and daughter), daughter’s history of oral contraceptive use and other contraceptive use, marital status, age at menarche, abortion history , age of mother at conception, and age of daughter at pregnancy or attempt were assessed. There was no evidence that these risk factors explained PCB associations  and so were not specifically investigated here.
It would have been ideal to measure maternal PCBs at multiple points in gestation in relation to daughter’s TTP 30 years later. However, it was our objective to conserve valuable maternal prenatal serum samples for other research questions that require precisely timed sampling, such as steroid hormones and thyroid hormones. We thus chose to use maternal postpartum samples as a proxy for daughter’s in utero exposure to PCBs. While it is possible that some PCB congeners are more subject to change during pregnancy , the postpartum blood samples in our study were drawn within 1–3 days of delivery, minimizing this source of misclassification. Based on rates of daily change in PCB levels reported by Bloom and colleagues , variability due to changes in PCB levels over three days after delivery would be small and considerably less than the substantial inter-individual variability of maternal PCB levels we observed in this sample . For example, based on mean daily rate of change estimates presented in Bloom et al., (, Table 4), these three days past delivery in our study would likely result in no or in extremely small differences between in utero PCB levels and measured PCB levels in postpartum samples collected so soon after delivery. Consistent with this argument, Longnecker et. al suggested that using postpartum maternal samples does not introduce serious misclassification in ranking daughters according to their in utero PCB exposures, since inter-individual variation in PCB levels greatly exceeds intra-individual variation due to timing of sampling .
One potential source of bias in this and our previous study is the retrospective assessment of time to pregnancy. As is typical, we chose a menstrual cycle as the time unit as each cycle offers a single ovulatory opportunity for pregnancy to occur . Numerous problems with this and other (i.e. prospective assessments) are described in the literature ([33–38]). As described in the methods, we asked each woman to report on her most recent pregnancy to avoid recall bias; nevertheless, as reported by Cooney et al.  and commented upon by Joffe , retrospective recall of TTP may result in misclassified true values, with the degree of measurement error proportional to gravidity, the time since pregnancy and duration of TTP. In the study of Cooney, et al. , where women who participated in a prospective study of TTP reported the same TTP retrospectively 10 years later, the distribution of reporting errors appears skewed to the left; however the sample size was small. In our study, where women were still reproductive age at interview (ages 27–31 years), the median time from the reported pregnancy to questionnaire was 3 years, with 25% reporting within 1 year and 75% reporting within 5 years of their most recent pregnancy or pregnancy risk period, suggesting that the bias may be smaller. As also indicated above in the methods, we included women who omitted birth control just once, women who became pregnant during an interval between birth control methods or after stopping birth control, women who became pregnant when breast feeding, women who never used birth control, women who were attempting pregnancy after a prior pregnancy, and women who never became pregnant, but had been at risk for pregnancy for a known interval, whether or not they desired to become pregnant. Many of these inclusions are an attempt to minimize other biases in the retrospective ascertainment of TTP as discussed by Joffe et al. , Key et al.  and Scheike and Keiding .
There are several studies which evaluate the associations between exposure to PCBs and other chlorinated compounds with regard to pregnancy. With one exception, previous studies base their exposure assessment on either fish consumption (paternal - ; maternal - ) or serum measures ([16, 29, 41–44]) proximal to the study. The notable exception is the analysis of exposure to contaminated cooking oil in the Yucheng cohort ; analyses here performed for women exposed during childhood and adolescence. Thus, our study ( and reported here) is unique in that it focuses on developmental exposure (i.e. during the in utero period) and fecundability outcomes in the female adult offspring. In Cohn et al. , we found longer TTP associated with prenatal exposure to PCBs 187, 156 and 99 and shorter TTP associated with prenatal exposure to PCBs 105, 138, and 183.
Women exposed to PCB/dibenzofuran/dibenzodioxins contaminated cooking oil  experienced prolonged TTP and reduced fertility, reported 24 years after the accident. Although TTP was queried for all pregnancies, only the first pregnancy was used in the analysis to minimize dependencies on subsequent TTP. TTP was defined as the duration between discontinuing contraception and the beginning of the last menstrual period prior to conception. TTP in exposed women was reduced by approximately 1 month compared to controls not exposed. We note that the Yucheng cohort has documented high exposure levels to this mixture of contaminants.
The New York State Angler’s Cohort is a prospective cohort study, which was initiated in 1991 as a result of concern regarding the role of environmental contaminants and reproductive outcomes. Prospective analyses using a subset of 102 women suggest that the weakly antiestrogenic PCB congener, PCB 205 and weakly estrogenic congener, PCB 206, as well as the antibacterial hexaclorobenzene were associated with increased TTP among women who became pregnant in the first three cycles following cessation of contraception . Further analyses found reduced fecundity for an a priori defined group of estrogenic PCB congeners (4_10, 5_8, 15_17, 18, 44, 47, 48, 52, 70, 99, 101, 136, 153 and 188) and a group of antiestrogenic PCB congeners (77_110, 105, 114, 126, 171_156, and 169), although the confidence limits surrounding the point estimates were wide and include the null .
An analysis of Swedish fishermen’s wives and sisters  found decreased TTP among those with PCB153 concentrations in the middle and high tertiles, compared to the lowest tertile, leading the investigators to posit that concentrations may not have been high enough to find an adverse association. All data collected in this study were retrospective and TTP was based on the first planned pregnancy. Moreover, PCB153 concentrations at the time of the first planned pregnancy were estimated based on a toxicokinetic model relating serum concentrations measured at the time of study, number of children, and breast feeding history. PCB153 is included in the group of estrogenic congeners but may not be the best for estimating exposure.
Data from the National Collaborative Perinatal Project  found a small association between exposure to the sum of 11 PCB congeners and TTP, such that women with concentrations in the highest quartile had longer TTP. TTP was assessed by questionnaire to the group of pregnant women. More recently, results from the LIFE study  reduced fecundability was associated with concentrations of individual PCB congeners (PCBs 118, 167, 209). Notably, the LIFE study prospectively recruited women attempting pregnancy following discontinuation of contraception. A final study considered the associations between PCBs and TTP among a subcohort of 332 pregnant women in the PELAGIE cohort . TTP was assessed via questionnaire. Reduced fecundability was found for exposure to the 10 PCB congeners (138, 153, 180, 118, 170, 187, 194, 183, and 203) and for the sum of the congeners. Many of the estimated associations were reduced with control for shellfish consumption.
In sum, there is some evidence of reduced fecundability and increased time to pregnancy in previous studies which consider both exposure and outcomes measured in adult women. Associations, however, are modest and conflicting in regard to which PCB congener or groups of congeners are important. Our study of in utero PCB exposure in relation to daughter’s fecundability is the only study to date where exposure was measured at this potentially critical window of susceptibility.
Results of the present study suggest that methods based on a priori classification of compounds could miss the contribution of compounds that have not been or cannot be classified. We suggest that empirical approaches, such as the one we report here, could supplement approaches based on expected biological activity or structure of mixture components. These supplemental empirical methods can generate new hypotheses and may stimulate the interest of toxicologists and basic scientists to investigate and explain the patterns reported. We suggest that the strategy of moving from observation in a human cohort to experimental science in the laboratory could accelerate discovery of new or different mechanisms.
Participants gave informed consent for a study on how early life factors influence reproduction with publications in aggregate form.
Child Health and Development Studies
National Health and Nutrition Examination Survey
A French birth cohort described in Chevrier et al. 
Toxic equivalent quotients
Time to pregnancy
Weighted quartile sum.
Kezios KL, Liu X, Cirillio PM, Kalantzi OI, Wang Y, Petreas MX, Park JS, Bradwin G, Cohn BA, Factor-Litvak P: Prenatal polychlorinated biphenyl exposure is associated with decreased gestational length but not birth weight: archived samples from the Child Health and Development Studies pregnancy cohort. Environ Health. 2012, 11: 49-10.1186/1476-069X-11-49.
Axelrad DA, Goodman S, Woodruff TJ: PCB body burdens in US women of childbearing age 2001–2002: an evaluation of alternate summary metrics of NHANES data. Environ Res. 2009, 109: 368-378. 10.1016/j.envres.2009.01.003.
Govarts E, Nieuwenhuijsen M, Schoeters G, Ballester F, Bloemen K, de Boer M, Chevrier C, Eggesbo M, Guxens M, Kramer U, Legler J, Martinez D, Palkovicova L, Patelarou E, Ranft U, Rautio A, Petersen MS, Slama R, Stigum H, Toft G, Trnovec T, Vandentorren S, Weihe P, Kuperus NW, Wilhelm M, Wittsiepe J, Bonde JP: Birth weight and prenatal exposure to polychlorinated biphenyls (PCBs) and dichlorodiphenyldichloroethylene (DDE): a meta-analysis within 12 European Birth Cohorts. Environ Health Perspect. 2012, 120: 162-170.
Cohn BA, Cirillo PM, Sholtz RI, Ferrara A, Park JS, Schwingl PJ: Polychlorinated biphenyl (PCB) exposure in mothers and time to pregnancy in daughters. Reprod Toxicol. 2011, 31: 290-296. 10.1016/j.reprotox.2011.01.004.
Baird DD, Wilcox AJ, Weinberg CR: Use of time to pregnancy to study environmental exposures. Am J Epidemiol. 1986, 124: 470-480.
Giusti RM, Iwamoto K, Hatch EE: Diethylstilbestrol revisited: a review of the long-term health effects. Ann Intern Med. 1995, 122: 778-788. 10.7326/0003-4819-122-10-199505150-00008.
Green DM, Kawashima T, Stovall M, Leisenring W, Sklar CA, Mertens AC, Donaldson SS, Byrne J, Robison LL: Fertility of female survivors of childhood cancer: a report from the childhood cancer survivor study. J Clin Oncol. 2009, 27: 2677-2685. 10.1200/JCO.2008.20.1541.
Dechanet C, Anahory T, Mathieu Daude JC, Quantin X, Reyftmann L, Hamamah S, Hedon B, Dechaud H: Effects of cigarette smoking on reproduction. Hum Reprod Update. 2011, 17: 76-95. 10.1093/humupd/dmq033.
Wolff MS, Camann D, Gammon M, Stellman SD: Proposed PCB congener groupings for epidemiological studies. Environ Health Perspect. 1997, 105: 13-14. 10.1289/ehp.9710513.
Myers RH: Classical and Modern Regression with Applications. 1990, Boston: PWS-KENT Publishing Co, 2
Carrico C: Characterization of a Weighted Quantile Score Approach for Highly Correlated Data in Risk Analysis Scenarios. 2013, PhD dissertation: Department of Biostatistics, Virginia Commonwealth University, Richmond, Virginia
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. 2008, 50: 1335-1342. 10.1097/JOM.0b013e31818f684d.
Gennings C, Sabo R, Carney E: Identifying subsets of complex mixtures most associated with complex diseases: polychlorinated biphenyls and endometriosis as a case study. Epidemiology. 2010, 21 (Suppl 4): S77-S84.
Christensen K, Carr C, Sanyal AJ, Gennings C: Mulitple classes of environmental chemicals are associated with liver disease: NHANES 2003–04. Int J Hyg Environ Health. 2013, Available online Mar 8 2013. [http://dx.doi.org/10.1016/j.ijheh.2013.01.005]
SAS: SAS.OR 9.22 User's Guide: Mathematical Programming. 2010, Cary NC, [http://support.sas.com/documentation/cdl/en/ormpug/63352/PDF/default/ormpug.pdf, page 367]
Buck Louis GM, Sundaram R, Schisterman EF, Sweeney AM, Lynch CD, Gore-Langton RE, Maisog J, Kim S, Chen Z, Barr DB: Persistent environmental pollutants and couple fecundity: the LIFE study. Environ Health Perspect. 2013, 121: 231-236.
Jensen TK, Joffe M, Scheike T, Skytthe A, Gaist D, Petersen I, Christensen K: Early exposure to smoking and future fecundity among Danish twins. Int J Androl. 2006, 29: 603-613. 10.1111/j.1365-2605.2006.00701.x.
Weinberg CR, Wilcox AJ, Baird DD: Reduced fecundability in women with prenatal exposure to cigarette smoking. Am J Epidemiol. 1989, 129: 1072-1078.
Jensen TK, Henriksen TB, Hjollund NH, Scheike T, Kolstad H, Giwercman A, Ernst E, Bonde JP, Skakkebaek NE, Olsen J: Adult and prenatal exposures to tobacco smoke as risk indicators of fertility among 430 Danish couples. Am J Epidemiol. 1998, 148: 992-997. 10.1093/oxfordjournals.aje.a009576.
Wilcox AJ, Baird DD, Weinberg CR: Do women with childhood exposure to cigarette smoking have increased fecundability?. Am J Epidemiol. 1989, 129: 1079-1083.
Joffe M, Barnes I: Do parental factors affect male and female fertility?. Epidemiology. 2000, 11: 700-705. 10.1097/00001648-200011000-00015.
Su PH, Huang PC, Lin CY, Ying TH, Chen JY, Wang SL: The effect of in utero exposure to dioxins and polychlorinated biphenyls on reproductive development in eight year-old children. Environ Int. 2012, 39: 181-187. 10.1016/j.envint.2011.09.009.
Mocarelli P, Gerthoux PM, Patterson DG, Milani S, Limonta G, Bertona M, Signorini S, Tramacere P, Colombo L, Crespi C, Brambilla P, Sarto C, Carreri V, Sampson EJ, Turner WE, Needham LL: Dioxin exposure, from infancy through puberty, produces endocrine disruption and affects human semen quality. Environ Health Perspect. 2008, 116: 70-77. 10.1289/ehp.116-a70.
Leijs MM, Koppe JG, Olie K, van Aalderen WM, Voogt P, Vulsma T, Westra M, ten Tusscher GW: Delayed initiation of breast development in girls with higher prenatal dioxin exposure; a longitudinal cohort study. Chemosphere. 2008, 73: 999-1004. 10.1016/j.chemosphere.2008.05.053.
Mouritsen A, Aksglaede L, Sorensen K, Mogensen SS, Leffers H, Main KM, Frederiksen H, Andersson AM, Skakkebaek NE, Juul A: Hypothesis: exposure to endocrine-disrupting chemicals may interfere with timing of puberty. Int J Androl. 2010, 33: 346-359. 10.1111/j.1365-2605.2010.01051.x.
Wang SL, Chang YC, Chao HR, Li CM, Li LA, Lin LY, Papke O: Body burdens of polychlorinated dibenzo-p-dioxins, dibenzofurans, and biphenyls and their relations to estrogen metabolism in pregnant women. Environ Health Perspect. 2006, 114: 740-745. 10.1289/ehp.8809.
Safe S, Wang F, Porter W, Duan R, McDougal A: Ah receptor agonists as endocrine disruptors: antiestrogenic activity and mechanisms. Toxicol Lett. 1998, 102–103: 343-347.
Cohn BA, Cirillo PM, Wolff MS, Schwingl PJ, Cohen RD, Sholtz RI, Ferrara A, Christianson RE, van den Berg BJ, Siiteri PK: DDT and DDE exposure in mothers and time to pregnancy in daughters. Lancet. 2003, 361: 2205-2206. 10.1016/S0140-6736(03)13776-2.
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: 523-532. 10.1093/aje/kwi240.
Cohn BA, Terry MB, Plumb M, Cirillo PM: Exposure to polychlorinated biphenyl (PCB) congeners measured shortly after giving birth and subsequent risk of maternal breast cancer before age 50. Breast Cancer Res Treat. 2012, 136: 267-275. 10.1007/s10549-012-2257-4.
Bloom MS, Buck Louis GM, Schisterman EF, Liu A, Kostyniak PJ: Maternal serum polychlorinated biphenyl concentrations across critical windows of human development. Environ Health Perspect. 2007, 115: 1320-1324. 10.1289/ehp.10086.
Longnecker MP, Klebanoff MA, Gladen BC, Berendes HW: Serial levels of serum organochlorines during pregnancy and postpartum. Arch Environ Health. 1999, 54: 110-114. 10.1080/00039899909602244.
Weinberg CR, Baird DD, Wilcox AJ: Sources of bias in studies of time to pregnancy. Stat Med. 1994, 13: 671-681. 10.1002/sim.4780130528.
Keiding N, Kvist K, Hartvig H, Tvede M, Juul S: Estimating time to pregnancy from current durations in a cross-sectional sample. Biostatistics. 2002, 3: 565-578. 10.1093/biostatistics/3.4.565.
Joffe M, Key J, Best N, Keiding N, Scheike T, Jensen TK: Studying time to pregnancy by use of a retrospective design. Am J Epidemiol. 2005, 162: 115-124. 10.1093/aje/kwi172.
Scheike TH, Keiding N: Design and analysis of time-to-pregnancy. Stat Methods Med Res. 2006, 15: 127-140. 10.1191/0962280206sm435oa.
Key J, Best N, Joffe M, Jensen TK, Keiding N: Methodological issues in analyzing time trends in biologic fertility: protection bias. Am J Epidemiol. 2009, 169: 285-293.
Cooney MA, Buck Louis GM, Sundaram R, McGuiness BM, Lynch CD: Validity of self-reported time to pregnancy. Epidemiology. 2009, 20: 56-59. 10.1097/EDE.0b013e31818ef47e.
Joffe M: Validity of self-reported time to pregnancy. Epidemiology. 2010, 21: 160-161. author reply 161
Buck GM, Mendola P, Vena JE, Sever LE, Kostyniak P, Greizerstein H, Olson J, Stephen FD: Paternal Lake Ontario fish consumption and risk of conception delay, New York State Angler Cohort. Environ Res. 1999, 80: S13-S18. 10.1006/enrs.1998.3926.
Chevrier C, Warembourg C, Gaudreau E, Monfort C, Le Blanc A, Guldner L, Cordier S: Organochlorine pesticides, polychlorinated biphenyls, seafood consumption, and time-to-pregnancy. Epidemiology. 2013, 24: 251-260. 10.1097/EDE.0b013e31827f53ec.
Buck GM, Vena JE, Greizerstein HB, Weiner JM, McGuinness B, Mendola P, Kostyniak PJ, Swanson M, Bloom MS, Olson JR: PCB congeners and pesticides and female fecundity, New York State Angler Prospective Pregnancy Study. Environ Toxicol Pharmacol. 2002, 12: 83-92. 10.1016/S1382-6689(02)00026-1.
Buck Louis GM, Dmochowski J, Lynch C, Kostyniak P, McGuinness BM, Vena JE: Polychlorinated biphenyl serum concentrations, lifestyle and time-to-pregnancy. Hum Reprod. 2009, 24: 451-458. 10.1093/humrep/den373.
Axmon ARL, Stromberg U, Jonsson B, Nilsson-Ehle P, Hagmar L: Polychlorinated biphenyls in serum and time to pregnancy. Environ Res. 2003, 96: 186-195.
Yang CY, Wang YJ, Chen PC, Tsai SJ, Guo YL: Exposure to a mixture of polychlorinated biphenyls and polychlorinated dibenzofurans resulted in a prolonged time to pregnancy in women. Environ Health Perspect. 2008, 116: 599-604. 10.1289/ehp.10715.
The authors gratefully acknowledge support from #T32 ES0007334, #UL1TR000058, #HHSN275201100020C, and #R01 ES0834.
The authors declare that they have no competing interests.
CG contributed to the statistical methods development, performed the statistical analysis, and wrote the first draft of the paper. CC contributed to the statistical methods development. PF-L participated in the design and coordination of the study and helped draft the manuscript. NK assisted in data management. PC and BC participated in the design and coordination of the study and helped draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Histograms for bootstrap distributions of weights for the PCB congeners in the weighted quartile score with a positive slope parameter in a Weibull proportional hazards model adjusted for race (African American vs all other) and whether the daughter was breast fed (yes or no). . Histograms for bootstrap distributions of weights for the PCB congeners in the weighted quartile score with a negative slope parameter in a Weibull proportional hazards model adjusted for race (African American vs all other) and whether the daughter was breast fed (yes or no). Table S1. Correlation estimates for the PCB congeners based on maternal serum concentrations (N = 289). (PDF 965 KB)
About this article
Cite this article
Gennings, C., Carrico, C., Factor-Litvak, P. et al. A Cohort study evaluation of maternal PCB exposure related to time to pregnancy in daughters. Environ Health 12, 66 (2013). https://doi.org/10.1186/1476-069X-12-66
- Endocrine disruptors
- Polychlorinated biphenyls
- Prenatal exposure