Abstract
Dichlorodiphenyldichloroethane (DDE) and polychlorinated biphenyls (PCB) are widespread environmental contaminants that have been postulated to increase the risk of diseases such as non–Hodgkin's lymphoma, breast cancer, as well as lead to early menopause. Studies assessing the effect of organochlorine exposure often can only measure organochlorine levels once, such as at study enrollment, which may not be an etiologically relevant time period. We assessed the temporal changes in DDE and PCBs and the predictors of those changes using interview data and DDE and PCB measures collected from 123 women who were enrolled in a baseline study from 1978 to 1982 and followed up in 2003 to 2004. Baseline and follow-up organochlorine levels were compared using Spearman correlations (rs), and predictors of the rate of change in log concentration were evaluated using linear regression models. Although serum concentrations dramatically declined (median follow-up to baseline concentration ratio was 16% for DDE and 45% for PCB), baseline and follow-up measures were strongly correlated for DDE (rs = 0.72) and moderately correlated for PCBs (rs = 0.43). Prediction of follow-up PCB levels was substantially improved (rs = 0.75) with data on initial concentration, length of lactation, baseline body mass index, and percent change in body fat, whereas DDE prediction improved slightly (rs = 0.83) with data on lactation and baseline body mass index. These findings suggest that a single organochlorine measure provides considerable information on relative ranking at distant times and that the predictive power can be improved, particularly for PCBs, with information on a few predictors. (Cancer Epidemiol Biomarkers Prev 2008;17(10):2729–36)
Introduction
Extensive use of the pesticide, dichlorodiphenyltrichloroethane (DDT), beginning in the 1940s and the industrial chemical, polychlorinated biphenyls (PCB), beginning in 1929, have led to their widespread distribution in the environment. Despite their ban since the 1970s, PCBs and the DDT metabolite, dichlorodiphenyldichloroethane (DDE), remain ubiquitous, raising concerns over their environmental and public health effect. A number of studies have evaluated associations of these compounds with human health effects, including non–Hodgkin's lymphoma (1), breast cancer (2, 3), and earlier menopause (4, 5), although clear associations have not been established.
Studies assessing the effect of these organochlorines on risk of late onset diseases or diseases with long latency periods often can only measure levels once; for example, they may be measured at enrollment into a prospective study of disease incidence. Although a single measure reflects exposure accumulation and excretion of a compound up to the time of sampling, it may not provide an accurate picture of the exposure level at a time of etiologic relevance, which may be years or decades apart from the time of sampling (6, 7).
The objective of this paper is to examine the ability to predict changes in organochlorine levels from baseline forward in time over a period of decades. Few studies to date have assessed long-term intraindividual changes in DDE or PCB concentrations and the predictors of those changes. These longitudinal studies (8-13) span ∼2 to 10 years in duration, and most have limited covariate data. Our study uses information collected from women who had enrolled in a baseline study of organochlorines (from 1978-1982) and were followed for ∼25 years.
Materials and Methods
Study Population
From 1978 to 1982, pregnant women (n = 877) from North Carolina were recruited into the baseline study of organochlorines and child development (14). Around the time of delivery, a questionnaire including information on demographics, prepregnancy weight, and reproductive and medical histories was administered; maternal and cord bloods, placenta, and breast milk were also collected for organochlorine determination. Additional samples were collected at subsequent study visits, including blood at the 6-wk visit and milk periodically for as long as the child was breastfed.
In 2003 to 2004, women were traced and invited into a follow-up study on menopause. Of 726 traceable women, 513 (70.7%) participated. Additional data including reproductive and lactation history (covering all pregnancies), height, and current weight were obtained via a telephone-administered questionnaire.
A sample of 126 individuals, constituting 82% with available blood sample from a 30% random sample of follow-up participants who were selected for a blood draw, had DDE and PCBs assessed in their serum samples. Of these, 123 had a serum organochlorine measure at delivery from the baseline study and will be the subjects for the current study. One individual had missing pregnancy data on the follow-up questionnaire and was excluded from our predictive modeling. Informed consent and Institutional Review Board approval were obtained for both baseline and follow-up studies.
Laboratory Assays
At baseline, p,p'-DDE and PCBs were assayed using packed column-gas chromatography with electron capture detection. PCBs were quantified as the sum of two peaks on a chromatogram (15). The coefficients of variation for serum and milk samples ranged from 3.7% to 11.2% for DDE and 3.1% to 9.2% for total PCBs (15). Due to the highly variable fat content in breast milk, milk concentrations were expressed as grams of organochlorine per gram of lipid, whereas concentrations in other tissues were reported per gram of tissue.
At follow-up, organochlorines were assayed using gas chromatography coupled with mass spectrometry. Fourteen congeners were quantitated; for this study, we focused on the total of 8 congeners (PCB# 99, 118, 138, 153, 156, 170, 180, and 187) having the lowest proportion of nondetects (<18%, compared with 52-99% for the remaining 6 congeners: PCB-28, 52, 101, 105, 128, and 183) and the highest relative abundance (together comprising ∼94% of all PCBs measured). The coefficients of variation were 5.1% for DDE and 5.6%, 4.1%, 3.5%, 3.8%, 11.4%, 4.7%, 3.5%, and 4.0%, respectively, for PCB-99, 118, 138, 153, 156, 170, 180, and 187. Total lipids were also estimated from total cholesterol, phospholipids, and triglycerides using the Phillips et al. approach (16), thereby allowing concentrations to be expressed per gram of lipid. A validation study shows that this approach gives values very similar to those obtained using the standard gravimetric lipid analysis (17).
Accounting for Change in PCB Quantitation Methods
Although the assays have changed since the baseline study, DDE quantitation (which involves measurement of a single peak) remained the same. However, PCB quantitation methods used at baseline differed from those used at follow-up and do not produce directly comparable results. To convert baseline PCB concentrations into their follow-up equivalents, we used data from 10 baseline milk samples reanalyzed by chromatography coupled with mass spectrometry (18). Baseline levels were multiplied by the median ratio (0.236127) of concentration using the follow-up method to concentration using the baseline method.
Handling Nondetects
One individual at baseline had nondetectable serum DDE and 18 (14.6%) had nondetectable serum PCBs. Serum organochlorine concentrations below the limit of detection at baseline were assigned estimated values using measurements from all available specimens at baseline. The concentrations from each specimen were previously scaled to reflect the concentration in milk lipids at birth of the index child and then averaged to create a composite measure for each woman (14). For this analysis, these composite measures were converted into their approximate serum equivalents by multiplying by the median ratio of serum to composite concentrations (5.0071367 for DDE and 5.4120890 for PCBs) obtained from members of the baseline population with detectable serum values. To evaluate the validity of this imputation, we compared estimated and actual levels among those with detectable concentrations in serum at delivery: mean % difference and SD were 6.1 + 36.9 for DDE and 3.8 + 24.3 for PCBs. Pearson and Spearman correlations were high between estimated and observed values for both DDE (0.90 and 0.91, respectively) and PCBs (0.91 and 0.84).
At follow-up, 1 individual had nondetectable serum DDE, whereas 28 (22.8%) had at least 1 nondetectable PCB congener. Nondetects at follow-up were assigned a value of the limit of detection divided by the square root of 2 (19). Because the nondetectable congeners were usually the minor ones, there were only 6 individuals in which the imputation method (using extremes of limit of detection or 0) made a >10% difference in their estimated total PCBs.
Accounting for Pregnancy-Related Increases at Baseline
As serum lipids were not available at baseline, our primary goal was to explore changes in wet weight serum organochlorine measures over time. Measures at baseline were acquired at or around the time of delivery; however, short-term changes due to pregnancy are not the changes of interest. Given the characteristic increase in lipid content across pregnancy (20, 21) and the resulting elevation in organochlorines, wet weight serum measures were adjusted to reflect estimated nonpregnant levels. Based on mean serum total lipids of 5.9 g/L among >2,000 nonpregnant Caucasian women, ages 20 to 34 y, participating in the 1976 to 1980 National Health and Nutrition Examination Survey (22, 23) and 8.5 g/L in a 1977 to 79 study of 553 pregnant Caucasian women, ages 20 to 41, at 36 wk of gestation (21), we estimated a 30% increase in total lipids across pregnancy. To account for this increase, wet weight serum concentrations at baseline were reduced by an equivalent amount.
Secondary Analysis: Lipid-Adjusted Serum Measures
Although our primary analysis was based on wet weight organochlorines, we conducted a secondary analysis on estimated lipid-adjusted measures. Lipid-adjusted serum concentrations at baseline were estimated by dividing the aforementioned “composite” milk values by a correction factor (1.5) that accounts for the difference in lipid concentrations between milk and serum. This factor is the approximate ratio of concentrations in milk lipids to those in serum lipids for both DDE (24, 25) and PCBs (24, 26).
Predictive Modeling
Although changes in serum DDE and PCB levels over time are driven by complex processes, we posit that the time course of such changes during the period under study can be reasonably approximated by an exponential decay. Thus, our basic predictive model was:
which can be rewritten as
where [Follow-up] and [Baseline] are the concentrations, respectively, at follow-up and baseline, and “follow-up time” represents years from baseline to follow-up. Slope represents change in log concentrations per unit of follow-up time. The slope can be affected by various determinants, so we modeled it using linear regression; the full model is then of the form:
where α is the intercept and βk is the regression coefficient for the k-th predictor, Xk.
Lactational transfer from mother to infant is an important route of organochlorine excretion (14, 27). In this study, lactation duration, from the follow-up questionnaire, is defined as the number of weeks she breast fed twice or more per day. Given the diminishing role of exclusive breastfeeding as lactation proceeds, the effect of breastfeeding on the rate of organochlorine excretion should lessen later in lactation. Thus, lactation duration was treated as a continuous variable and modeled in a piecewise linear fashion with the points where the slope changes prespecified at 26 and 39 wk of lactation. This was accomplished by defining 3 variables representing the number of weeks of lactation during the 3 intervals: 0 to 26th, 27th to 39th, and 40th or more week of lactation. The number of weeks of lactation for each pregnancy from baseline to follow-up was partitioned among these three variables. For women with multiple lactations between baseline and follow-up, the contribution from each pregnancy was summed. To illustrate, for a woman with one child whose lactation lasted 60 wk, a value of 26, 13, and 21 wk was assigned to the first, second, and third lactation variables, respectively. If she also had a second child who was breastfed for 30 wk, then an additional 26 wk was assigned to the first lactation variable and 4 wk to the second, giving final values of 52, 17, and 21 for the three lactation variables. To operationalize, the 3 pieces contribute to the overall slope of organochlorine decline based on their regression coefficients: slope contribution = β1(piece 1) + β2(piece2) + β3(piece3).
Previous studies on PCBs and dioxins show that high initial concentration was associated with faster decline over time (28, 29). We included the log of initial (baseline) concentration of serum DDE or PCBs as a linear variable in the model.
Body mass index (BMI) may be associated with an individual's ability to metabolize organochlorines and other xenobiotics, with heavier individuals having slower metabolism (30). Baseline BMI was included in the model as a categorical predictor, with cutpoints at 20 and 23 kg/m2, due to its narrow distribution.
Because fat is where organochlorines are sequestered, an increase in body fat without an increase in organochlorines might dilute organochlorine concentrations. Thus, we included a measure of percent change in body fat (kg) based on a validated formula for percent body fat, which we then converted to body fat by multiplying by body weight (31):
Usual weight before baseline study pregnancy and current weight were obtained, respectively, from the baseline and follow-up questionnaires. Percentage change in body fat, defined as the difference in body fat from baseline to follow-up divided by baseline body fat, was included as a linear variable.
Mother's date of birth reflects the secular trend in environmental DDE and PCB levels and, hence, the exposure potential. We included calendar date of birth (measured in days with a SAS reference date of January 1, 1960) as a linear variable in the model. Another possible predictor is baseline maternal age; organochlorine metabolism may change with age. However, because age at baseline is highly correlated with date of birth (r = 0.97), it was not included.
Statistical Analysis
Linear regression models were fit using SAS version 9.1. Each predictor was examined to obtain the fewest categories needed to adequately characterize the relationship between that variable and the slope. The predictors were added to the model with initial concentration first followed by other predictors in successive order of their assumed strength of effect. The statistical significance of each potential predictor was assessed using the F test. Spearman correlations comparing actual versus model-predicted concentrations at follow-up were used to assess the predictiveness of each model. To illustrate the model fits, model-predicted changes in organochlorine levels were plotted along with actual changes for the 122 individuals with nonmissing questionnaire data.
Results
The mean age of our study sample at baseline was 29 years. Other descriptive information is shown in Table 1. Among the 122 women with pregnancy data, there were 207 births during the study period (i.e., from baseline to follow-up, including baseline pregnancy) and 39% of these are the only birth during the study period. Ninety percent of the babies were breastfed, and 52% were breastfed for >26 weeks. Approximately 90% of women had a prepregnancy BMI of <25 kg/m2 at baseline, and all but two had increased their body fat mass (65% with at least a 50% increase) by the time of follow-up.
Variable . | n . | % . | ||
---|---|---|---|---|
Age at baseline (y) | ||||
19-24 | 24 | 19 | ||
25-29 | 54 | 44 | ||
30-34 | 39 | 32 | ||
35-37 | 6 | 5 | ||
No. of births (baseline to follow-up) | ||||
1 | 48 | 39 | ||
2 | 56 | 46 | ||
3 | 17 | 14 | ||
6 | 1 | 1 | ||
No. of babies breast-fed (baseline to follow-up)* | ||||
0 | 12 | 10 | ||
1 | 45 | 37 | ||
2 | 48 | 39 | ||
3 | 16 | 13 | ||
6 | 1 | 1 | ||
No. weeks of lactation for each birth (baseline to follow-up)* | ||||
0-13 | 49 | 24 | ||
14-26 | 50 | 24 | ||
27-39 | 31 | 15 | ||
40-52 | 54 | 26 | ||
53-165 | 23 | 11 | ||
Baseline prepregnancy BMI (kg/m2) | ||||
16.6-19.9 | 47 | 38 | ||
20-22.9 | 52 | 42 | ||
23-32.9 | 24 | 20 | ||
Percent change in body fat | ||||
−6 to 0 | 2 | 2 | ||
1-24 | 13 | 11 | ||
25-49 | 27 | 22 | ||
50-99 | 49 | 40 | ||
100-149 | 18 | 14 | ||
150-294 | 14 | 11 |
Variable . | n . | % . | ||
---|---|---|---|---|
Age at baseline (y) | ||||
19-24 | 24 | 19 | ||
25-29 | 54 | 44 | ||
30-34 | 39 | 32 | ||
35-37 | 6 | 5 | ||
No. of births (baseline to follow-up) | ||||
1 | 48 | 39 | ||
2 | 56 | 46 | ||
3 | 17 | 14 | ||
6 | 1 | 1 | ||
No. of babies breast-fed (baseline to follow-up)* | ||||
0 | 12 | 10 | ||
1 | 45 | 37 | ||
2 | 48 | 39 | ||
3 | 16 | 13 | ||
6 | 1 | 1 | ||
No. weeks of lactation for each birth (baseline to follow-up)* | ||||
0-13 | 49 | 24 | ||
14-26 | 50 | 24 | ||
27-39 | 31 | 15 | ||
40-52 | 54 | 26 | ||
53-165 | 23 | 11 | ||
Baseline prepregnancy BMI (kg/m2) | ||||
16.6-19.9 | 47 | 38 | ||
20-22.9 | 52 | 42 | ||
23-32.9 | 24 | 20 | ||
Percent change in body fat | ||||
−6 to 0 | 2 | 2 | ||
1-24 | 13 | 11 | ||
25-49 | 27 | 22 | ||
50-99 | 49 | 40 | ||
100-149 | 18 | 14 | ||
150-294 | 14 | 11 |
NOTE: One individual had missing pregnancy and lactation data.
Among the 207 births delivered to 122 mothers from baseline to follow-up.
At baseline, the median wet weight concentrations were 8.5 μg/L for DDE and 1.5 μg/L for PCBs (Table 2). At follow-up, medians fell to 1.2 and 0.7 μg/L, respectively; expressing each follow-up concentration as a percent of her baseline, medians were 16% for DDE and 45% for PCBs. There was a dramatic decline in DDE levels, with 90% of women having follow-up concentrations <30% of baseline, whereas the variability in individual changes was greater for PCBs. Changes over time in lipid-adjusted measures were of similar magnitude.
Organochlorine . | Min . | Percentile . | . | . | . | . | Max . | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | 10 . | 25 . | 50 . | 75 . | 90 . | . | |||||||
Wet weight measures | ||||||||||||||
DDE (μg/L serum) | ||||||||||||||
Baseline | 1.3 | 3.8 | 5.8 | 8.5 | 12.5 | 19.8 | 41.0 | |||||||
Follow-up | 0.01 | 0.5 | 0.7 | 1.2 | 2.3 | 3.9 | 11.0 | |||||||
PCBs (μg/L serum) | ||||||||||||||
Baseline | 0.6 | 0.9 | 1.2 | 1.5 | 2.1 | 2.7 | 5.2 | |||||||
Follow-up | 0.2 | 0.4 | 0.5 | 0.7 | 1.0 | 1.2 | 2.2 | |||||||
Lipid-adjusted measures | ||||||||||||||
DDE (μg/L serum) | ||||||||||||||
Baseline | 0.3 | 0.7 | 1.2 | 1.7 | 2.5 | 3.7 | 7.4 | |||||||
Follow-up | 0.002 | 0.1 | 0.1 | 0.2 | 0.4 | 0.5 | 1.9 | |||||||
PCBs (μg/L serum) | ||||||||||||||
Baseline | 0.1 | 0.2 | 0.2 | 0.3 | 0.4 | 0.5 | 0.7 | |||||||
Follow-up | 0.03 | 0.1 | 0.1 | 0.1 | 0.1 | 0.2 | 0.3 |
Organochlorine . | Min . | Percentile . | . | . | . | . | Max . | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | 10 . | 25 . | 50 . | 75 . | 90 . | . | |||||||
Wet weight measures | ||||||||||||||
DDE (μg/L serum) | ||||||||||||||
Baseline | 1.3 | 3.8 | 5.8 | 8.5 | 12.5 | 19.8 | 41.0 | |||||||
Follow-up | 0.01 | 0.5 | 0.7 | 1.2 | 2.3 | 3.9 | 11.0 | |||||||
PCBs (μg/L serum) | ||||||||||||||
Baseline | 0.6 | 0.9 | 1.2 | 1.5 | 2.1 | 2.7 | 5.2 | |||||||
Follow-up | 0.2 | 0.4 | 0.5 | 0.7 | 1.0 | 1.2 | 2.2 | |||||||
Lipid-adjusted measures | ||||||||||||||
DDE (μg/L serum) | ||||||||||||||
Baseline | 0.3 | 0.7 | 1.2 | 1.7 | 2.5 | 3.7 | 7.4 | |||||||
Follow-up | 0.002 | 0.1 | 0.1 | 0.2 | 0.4 | 0.5 | 1.9 | |||||||
PCBs (μg/L serum) | ||||||||||||||
Baseline | 0.1 | 0.2 | 0.2 | 0.3 | 0.4 | 0.5 | 0.7 | |||||||
Follow-up | 0.03 | 0.1 | 0.1 | 0.1 | 0.1 | 0.2 | 0.3 |
NOTE: PCB concentration comprises PCB-99, 118, 138, 153, 156, 170, 180, and 187.
The correlations between the two wet weight measures taken >20 years apart were relatively strong (Fig. 1). DDE at baseline was highly correlated with DDE at follow-up (Spearman rs = 0.72), whereas the correlation for PCBs was lower but still sizable (rs = 0.43). Correlations were similar for lipid-adjusted measures (0.73 and 0.50, respectively).
Higher initial concentration, lactation before 39 weeks, and an increase fat mass were associated with an increase in the rate of elimination of wet weight DDE and PCBs, whereas higher baseline BMI tended to decrease the rate of decay (Table 3). Similar results were found for lipid-adjusted values (data not shown). Mother's date of birth is not presented, as its effect was not significant for either DDE or PCBs. Table 4 shows the statistical significance of successive additional terms in the model for the slope and the correlation of the resulting predicted follow-up concentration with observed concentration. Lactation duration and baseline BMI were the only statistically significant predictors of DDE slope, whereas for PCBs, initial concentration and percent change in body fat were also significant. Predicted follow-up concentrations from the full model were highly correlated with actual follow-up concentrations for both DDE (rs = 0.83) and PCBs (rs = 0.75; Fig. 2; Table 4).
Predictors . | DDE . | . | PCB . | . | ||||
---|---|---|---|---|---|---|---|---|
. | β . | SE . | β . | SE . | ||||
Intercept | −0.061983 | 0.0100 | −0.012182 | 0.0036 | ||||
Log of initial concentration (μg/L) | −0.004270 | 0.0036 | −0.016655 | 0.0031 | ||||
Total weeks of breastfeeding in each lactation period | ||||||||
0-26th | −0.000202 | 0.0001 | −0.000201 | 0.0001 | ||||
27th-39th | −0.000295 | 0.0003 | −0.000427 | 0.0001 | ||||
>39th | 0.000052 | 0.0001 | −0.000086 | 0.0001 | ||||
Baseline BMI (kg/m2) | ||||||||
16.6-19.9 | Reference | Reference | Reference | Reference | ||||
20−22.9 | 0.005467 | 0.0053 | 0.002533 | 0.0029 | ||||
≥23 | 0.031238 | 0.0063 | 0.011727 | 0.0036 | ||||
Percent change in body fat | −0.008504 | 0.0046 | −0.006323 | 0.0025 |
Predictors . | DDE . | . | PCB . | . | ||||
---|---|---|---|---|---|---|---|---|
. | β . | SE . | β . | SE . | ||||
Intercept | −0.061983 | 0.0100 | −0.012182 | 0.0036 | ||||
Log of initial concentration (μg/L) | −0.004270 | 0.0036 | −0.016655 | 0.0031 | ||||
Total weeks of breastfeeding in each lactation period | ||||||||
0-26th | −0.000202 | 0.0001 | −0.000201 | 0.0001 | ||||
27th-39th | −0.000295 | 0.0003 | −0.000427 | 0.0001 | ||||
>39th | 0.000052 | 0.0001 | −0.000086 | 0.0001 | ||||
Baseline BMI (kg/m2) | ||||||||
16.6-19.9 | Reference | Reference | Reference | Reference | ||||
20−22.9 | 0.005467 | 0.0053 | 0.002533 | 0.0029 | ||||
≥23 | 0.031238 | 0.0063 | 0.011727 | 0.0036 | ||||
Percent change in body fat | −0.008504 | 0.0046 | −0.006323 | 0.0025 |
Model* . | Predictors . | Statistical significance (P) of successive terms in model for slope . | . | Spearman correlation for observed and predicted follow-up concentration . | . | ||
---|---|---|---|---|---|---|---|
. | . | DDE . | PCBs . | DDE . | PCBs . | ||
I | Intercept | N/a | N/a | 0.72 | 0.43 | ||
II | Log of initial concentration (μg/L) | 0.75 | <0.01 | 0.73 | 0.46 | ||
III | Weeks of lactation | 0.05 | <0.01 | 0.75 | 0.67 | ||
IV | Baseline BMI | <0.01 | 0.01 | 0.82 | 0.72 | ||
V | Percent change in body fat | 0.07 | 0.01 | 0.83 | 0.75 |
Model* . | Predictors . | Statistical significance (P) of successive terms in model for slope . | . | Spearman correlation for observed and predicted follow-up concentration . | . | ||
---|---|---|---|---|---|---|---|
. | . | DDE . | PCBs . | DDE . | PCBs . | ||
I | Intercept | N/a | N/a | 0.72 | 0.43 | ||
II | Log of initial concentration (μg/L) | 0.75 | <0.01 | 0.73 | 0.46 | ||
III | Weeks of lactation | 0.05 | <0.01 | 0.75 | 0.67 | ||
IV | Baseline BMI | <0.01 | 0.01 | 0.82 | 0.72 | ||
V | Percent change in body fat | 0.07 | 0.01 | 0.83 | 0.75 |
Abbreviation: N/a, not applicable.
Each model includes all variables from the preceding model.
The Supplementary figure illustrates the predicted changes in DDE and PCB levels, and the effect of lactation and baseline BMI. Overall, there is a rapid decline in DDE and shallower reduction in PCB. The effect of the heaviest versus lightest category of baseline BMI was considerable for both DDE and PCBs. The drop in concentrations as a result of breastfeeding 1 child for 1 year versus no breastfeeding was evident for both compounds, although weaker for DDE. Although not shown, the effect of initial concentration (75th versus 25th percentile) on PCB decline was also notable, but the effect of an increase in fat mass was comparatively small.
Discussion
We found a substantial decline in intraindividual DDE levels over an approximate 25-year time span (1978-1982 to 2003-2004), with a smaller drop in PCB concentrations. These declines correspond with the reduction in environmental levels because of the banning of DDT in 1972 and PCBs in 1977 in the United States (32-35). The larger decline in DDE compared with PCBs is similar to those of previous studies conducted among Swedish men from 1991 to 2001 (12) and Great Lakes fisheaters from 1982 to 1989 (9). In shorter studies, with median follow-up duration of 25.4 months (11) and 5 years (8), the differences in the decline between DDE and PCBs were less obvious. Peak production and the restrictions on usage occurred later for PCBs than for DDT and may partly explain the smaller decline in PCBs. Differences between regions and sampling time because of organochlorine restrictions may also contribute to the differences between studies.
Our baseline levels of DDE were highly correlated with those at follow-up (Spearman r = 0.72), whereas correlations were lower for PCBs. Higher correlations for DDE than for PCBs were also found in the Wolff et al. (11) study with a median of 25.4 months of follow-up (0.95 versus 0.83, respectively) and the Hoyer et al. (8) study with a maximum of 7 years follow-up (0.79 versus 0.64, respectively). These differences became less evident when the specific congener PCB-153 was assessed; DDE and PCB-153 correlations were respectively 0.92 and 0.90 in a 1991 to 2001 study by Hagmar et al. (12) and 0.79 and 0.68 in the Hoyer et al. (8) study. The latter study also showed that correlations varied by PCB congener. Individual congeners could not be evaluated in our study because they were not measured at baseline, and the fact that total PCBs is a mix of congeners with different rates of decomposition may partly explain the lower correlations in PCBs.
A number of previous studies have used occupational cohorts or individuals exposed to very high, acute doses to estimate half-lives of organochlorines, with estimates being ∼6 to 10 years for DDE (36, 37) and a few weeks to over 10 years for individual PCB congeners (38). We examined a group of women with protracted, albeit declining, environmental exposures to organochlorines, with attention to the influence of individual factors such as lactation and weight fluctuations on DDE and PCB levels. Based on the predicted curves (Supplementary Figure), the apparent half-lives we predict during this time period are ∼10 years for DDE and 25 years for PCBs, although they vary for women with different levels of the predictors. Because exposure is ongoing, although declining, in our study population, half-life estimation is expected to be longer compared with the aforementioned studies of occupational cohorts.
Breastfeeding was important in the elimination of organochlorines in previous studies of ∼2 years (14) and 4 years of follow-up (39), and our study indicates that the effect of lactation still exists ∼25 years later. In contrast to previous studies, our study attempted to differentiate the elimination potential of different intervals of lactation. In our study, lactation beyond 9 months had less of an effect on DDE and PCB levels than the earlier months. This change in effect across intervals of lactation is expected, given that breastfeeding is likely to diminish in the later months as solid foods become an important component of diet. The effect of one year of lactation was more dramatic for PCBs than for DDE (Supplementary Figure).
We found an association between higher initial concentration and faster elimination of PCBs, consistent with studies among those occupationally exposed (29, 40). A study of 701 women environmentally exposed to PCBs (median, 31 months follow-up) also reported similar findings; those whose PCBs had decreased were over thrice more likely to have higher initial concentrations than those who had maintained their levels (10). A possible explanation for this phenomenon that is most likely in our population is higher induction of enzymes at higher doses (28). Although faster elimination of DDT from adipose tissue was associated with higher initial concentration in an experimental study of three individuals taking varying high level doses of DDT (41), similar data for DDE is not available.
Baseline BMI had a strong effect on the rate of change of both DDE and PCBs. This is consistent with BMI-dependent variation in DDE levels in previous studies (11, 13), showing a positive association or correlation between a single BMI measure and DDE half-life. Sweeney et al. (10) also found a small positive association between baseline BMI and an increase in PCB levels, although another study did not (11). It is possible that the effect of higher BMI on the rate of decline may be associated with reduced metabolism by xenobiotic-metabolizing enzymes, such as cytochrome P450 (CYP). CYPs of the 1A, 2B, and 3A family are most likely involved in the metabolism of various PCB congeners (33, 42, 43), and obesity has been associated with decreased activity in CYP 3A4 in several studies (30).
It should be noted that despite the narrow range of BMI at baseline, the effect of a single measure of BMI was impressive. The relationship between baseline BMI and organochlorine changes over time could not be explained by variations in baseline age or initial concentration. Whether the effect of baseline BMI would be stronger at higher BMI remains unclear.
Studies assessing the influence of weight change on organochlorine change have been few, particularly for long term trends (8, 12, 44-46). Increased fat as a result of weight gain may lead to the dilution of organochlorine concentrations, thereby resulting in a negative association between body fat and serum organochlorine levels (47). Our finding of a significant effect of fat mass change for PCBs agrees with other studies (with 5 or more years of follow-up) that use relative change in BMI (12) or absolute weight change (8) as the predictors of interest.
The results for BMI and fat mass change suggest that different aspects of weight function differently in the elimination of serum organochlorines and, therefore, should be distinguished. Conceptually, BMI may be a marker of metabolism, whereas fat mass change may be a reflection of passive dilution of organochlorines.
Mother's date of birth, and consequently, age, was not a significant predictor of change for either DDE or PCBs, which mirrors the lack of association found in other studies between age at study entry and change in DDE (12) or PCB levels (10). As postulated by several theoretical models, the exposure patterns over time are a result of complex interactions between cohort-related exposures and breastfeeding patterns as well as age-dependent growth and weight/fat content, dietary intake and composition, and metabolism (47, 48). However, given the small variation in birth year and age in our study population and the fact that everyone was born before the ban, we would need a much larger sample size to more fully assess the effect of birth cohort and age.
Although the factors considered in our study help to explain a large portion of the variation in DDE and PCB levels over time, other potentially important covariates such as changes in diet, particularly fish consumption, were not measured. A longitudinal study (35) conducted from 1980 to 1995 and a serial cross-sectional study from 1973 to 1993 (49) of Great Lakes fisheaters and “noneaters” showed fish consumption in the previous year to be predictive of PCB body burden. Likewise, whale blubber consumption in the Faroes Islands have been associated with increased levels of highly chlorinated PCBs (50).
A limitation of this study is that baseline and follow-up samples were analyzed by different labs at different times using different assays. The most important difference, however, is the change in PCB quantitation methods. To obtain comparable PCB values from the old and new techniques, we used conversion factors derived from a set of 10 breast milk samples spanning the range of organochlorine values. Although derived from a small sample, our conversion factor seemed reasonable, given the similarity in values and high Spearman and Pearson correlations (0.83 and 0.96, respectively) between the original and converted values.
We analyzed both estimated wet weight and lipid-adjusted measures. These measures involved various assumptions, including those regarding methodology, lipid concentration, and nondetects. Similar results for both analyses, nevertheless, provide further assurance of our outcome and conclusions.
“Weeks of lactation” was modeled as 3 periods of lactation duration, and additional cutpoints (at 13 and 52 weeks) did not significantly alter the results. However, “weeks of lactation” alone does not fully capture the extent to which organochlorines are being transferred, as there are no data on the amount of milk expressed per breastfeeding session or the duration of weaning for all relevant pregnancies. Additionally, for women who fed multiple children, we assumed the effect of lactation would not vary across pregnancies; however, it is not possible to evaluate this assumption given the small sample size and the lack of organochlorine measurement at each pregnancy.
Although our findings agree with those of previous studies, there are limits to the generalizability. Our study population consists of women who were highly educated, mostly Caucasian (95%) and pregnant at baseline. We accounted for the changes in lipid profile due to pregnancy by reducing the baseline serum concentrations by 30%. We did a sensitivity analysis using, instead, ranges of 20% to 50%. The magnitude of the reduction in levels over the follow-up period was affected (i.e., median relative change in levels using a 20%, 30%, and 50% weight were, respectively, 14%, 16%, and 23% for DDE and 39%, 45%, and 63% for PCBs); however, there was virtually no effect on the model predictors.
Our participants were also relatively young and may differ from those not studied in terms of their age- and cohort-related exposure potential, dietary pattern, BMI, and metabolism. More importantly, our results may not be generalizable to other regions of the world with vastly different exposures, and the findings can be extrapolated beyond the time frame of our study (1978-1982 to 2003-2004) only with caution. Although our assumption of an exponential decay is plausible, the shape of the concentration curve between the two measurements cannot be determined given only two measures. Multiple measures would be required to more fully evaluate the shape. Nevertheless, our data may be useful to those conducting studies on the role of these organochlorines in the development of chronic diseases and have only one DDE or PCB measure, at perhaps, a time period that is not optimal in terms of etiologic relevance.
We have focused on predicting future levels based on past measurements. Others have attempted to predict past exposure from current concentrations (51, 52). Although our model coefficients cannot be directly used for backward estimations because initial (baseline) concentration, itself, is a predictor in the model, our study is useful in identifying a handful of important predictors to consider in such estimation models. A single measure of DDE is highly predictive of a woman's relative exposure over a time span of ∼25 years, and can be further improved with data on lactation and baseline BMI. A single measure of PCBs is also predictive of future levels. Although initial concentration and percent body fat change contribute to the prediction of PCB, but not DDE values, lactation duration and baseline BMI are the predominant contributors to both DDE and PCBs and information on these predictors can be obtained retrospectively.
This study addressing intraindividual variations and influences on body burden of DDE and PCBs over time in a general population with relatively low level environmental exposure includes measures spanning ∼25 years apart. Long term trends and predictions are of interest for public health assessment as well as studies of cancer and other diseases with long latency periods, where the levels of etiologic relevance may be decades apart from the time of measurement. Future studies to replicate or confirm our findings (particular in other periods and populations) are warranted and would benefit from multiple organochlorine measures and additional information on other potential predictors.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Grant support: Intramural Research Program of the NIH, National Institute of Environmental Health Sciences.
Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).
Acknowledgments
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.