Application of ensemble methods to analyse the decline of organochlorine pesticides in relation to the interactions between age, gender and time
Authors:
Aleysha Thomas aff001; Nicole M. White aff001; Leisa-Maree Leontjew Toms aff002; Kerrie Mengersen aff001
Authors place of work:
ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, Queensland, Australia
aff001; School of Public Health and Social Work, Queensland University of Technology, Brisbane, Queensland, Australia
aff002
Published in the journal:
PLoS ONE 14(11)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0223956
Summary
Organochlorine pesticides (OCPs) are toxic chemicals that persist in human tissue. Short and long term exposure to OCPs have been shown to have adverse effects on human health. This motivates studies into the concentrations of pesticides in humans. However these studies typically emphasise the analysis of the main effects of age group, gender and time of sample collection. The interactions between main effects can distinguish variation in OCP concentration such as the difference in concentrations between genders of the same age group as well as age groups over time. These are less studied but may be equally or more important in understanding effects of OCPs in a population. The aim of this study was to identify interactions relevant to understanding OCP concentrations and utilise them appropriately in models. We propose a two stage analysis comprising of boosted regression trees (BRTs) and hierarchical modelling to study OCP concentrations. BRTs are used to discover influential interactions between age group, gender and time of sampling. Hierarchical models are then employed to test and infer the effect of the interactions on OCP concentrations. Results of our analysis show that the best fitting model of an interaction effect varied between OCPs. The interaction between age group and gender was most influential for hexachlorobenzene (HCB) concentrations. There was strong evidence of an interaction effect between age group and time for β-hexachlorocyclohexane (β-HCH) concentrations in >60 year olds as well as an interaction effect between age group and gender for HCB concentrations for adults aged >45 years. This study highlights the need to consider appropriate interaction effects in the analysis of OCP concentrations and provides further insight into the interplay of main effects on OCP concentration trends.
Keywords:
Lipids – Age groups – Normal distribution – Water pollution – Decision trees – Pollutants – Pesticides
Introduction
Organochlorine pesticides (OCPs) are a class of lipophilic persistent organic pollutants (POPs) that exist in the environment and accumulate in human tissue [1–4]. POPs are subject to long range transport via air, water and dust [4, 5] and the concentrations of POPs are known to degrade slowly in humans [6, 7]. In Australia, pesticides including dichlorodiphenyltrichloroethane (DDT), hexachlorocyclohexanes (HCHs) and hexachlorobenzene (HCB) were routinely used for agricultural applications prior to their deregistration in 1985 [8, 9]. Short and long-term exposure to POPs, including OCPs, has been associated with negative health impacts such as cancers and lymphoma, as well as reproductive and neurological effects [10, 11]. These adverse effects have motivated epidemiological studies into POP concentrations in blood serum [12–14], other tissues [15, 16], and changes in concentrations over time [17–22].
The individual effects of age, gender and time of sample collection on POP concentrations are a widely studied subject of monitoring reports. In studies conducted on serum and breast milk samples collected in the late 1990’s and early 2000’s (after the production and use of the POPs) researchers around the world, including in Australia [17, 18], have identified that POP concentrations are higher in older adults [17, 23–25], vary between genders [17, 26, 27] and generally decrease over time [17–20]. This paper analyses five OCP concentrations sampled from children and adults in Australia between 2002 and 2012 [22]. Beyond individual effects, there are gaps in our understanding of how the interactions between age, gender and time affect the levels of POP concentrations in a population [21]. For example, an interaction between age group and gender would indicate differing concentrations between men and women of the same age group. Similarly, an interaction between age group and time of sample collection would suggest non-parallel age-specific trends in OCP concentrations over time and provide insight into these trends as the population ages. Studies in the field have not yet distinguished the effects of such interactions on the trend of OCP concentrations [21].
Interactions between age group, gender and time of sample collection on long term OCP concentrations have previously been statistically evaluated by Thomas et al. [22]. Similar data has been studied by Woodruff et al [28]. Using analysis of variance (ANOVA) and multiple comparisons, these authors identified significant interactions between age group and gender as well as time of sample collection and gender. Other studies have used similar methods to identify trends in OCP concentrations but did not investigate interaction effects [17, 18].
The paper proposes a two-stage statistical model for the discovery of influential interaction effects between the key covariates age group, gender and time. In the first stage of the analysis, boosted regression trees (BRTs) are applied to elucidate interactions for each OCP independently. Resulting interactions are incorporated into an hierarchical model for describing temporal trends in OCP concentrations as a function of the interactions. BRTs are a flexible, data-driven method to identify interactions that affect the response variable. Hierarchical models are a generalisation of linear models where parameters can vary at more than one level. The proposed hierarchical model in this study is formulated within the Bayesian framework so that a range of prior distributions can be used to incorporate the structure of the data into the model.
Materials and methods
This section gives a brief summary of the data analysed followed by a description of each model of the two stage analysis namely, BRTs and hierarchical models. The methods adopted to evaluate the performance of the models are briefly summarised here.
Data
The two stage analysis was applied to concentrations of HCB, β-HCH, trans-nonachlor, p,p’-DDE and p,p’-DDT measured from pooled samples of de-identified surplus human blood serum from males and females collected in Brisbane, Australia in 2002/03, 2006/07, 2008/09, 2010/11 and 2012/13 across the age groups 5-15, 16-30, 31-45, 46-60 and >60 years. The age group 0-4 years was available but excluded from the analysis due to missing data for the samples collected in 2002/03.
Each pool consisted of 100 individual serum samples for 2002/03, 2008/09, 2010/11 and 2012/13. Pooled samples collected in 2006/07 comprised up to 30 samples per pool. A total of 12,175 individuals made up of 183 pools. There are more pools for children aged 5-15 years in 2006/07 as these samples were originally collected for a different study as pools for the age groups 5-9, 9-12 ad 12-15 years [29].
OCP concentrations in serum are broadly representative of a person’s body burden of OCPs. There are factors which can affect an individual’s body burden such as exposure to OCP through occupation or lifestyle as well as factors affecting elimination such as breastfeeding or lipid mobilization as occurs during weight gain and loss [30]. An advantage of using pooled samples is that is has the effect of diluting any high or low concentrations to estimate the ‘average’ OCP concentration in the pool and eliminate variability.
Further details on data collection are provided in the Supporting Information (S1 File) and a summary of the data are provided in Thomas et al. [22]. The number of pools sampled at each time for each gender and age group is summarised in Table 1 below.
Ethics approval for this study was granted by The University of Queensland (UQ) Medical Research Ethics Committee and Queensland University of Technology (QUT) Ethics Committee (Ethics Approval Number: 2013000317). As per the ethics approval, the samples were de-identified and pooled therefore no consent was required from the participants.
Boosted regression trees
Boosted regression trees (BRTs) offer a non-parametric approach to predictive modelling by combining Classification and Regression trees (CART) [31] with a boosting algorithm [32] to refine model predictions. A CART is a decision tree that determines covariate-based rules to recursively split the data into subgroups of observations with a similar response [33]. Boosting is a stage wise process whereby multiple simple models are combined to enhance predictive performance [34]. In a BRT, the loss in predictive performance is represented by a loss function. For continuous responses, changes in deviance is a commonly chosen loss function [35]. The initial regression tree has the lowest possible loss function. Using the boosting algorithm, regression trees are sequentially applied to the residuals of the previous tree until there is no further improvement in predictive performance [36].
BRTs were adopted in this paper to determine the relative importance of age group, gender and time, and their interaction, for predicting OCP concentrations. For the data described in the earlier section, BRTs were applied separately to concentration measurements for HCB, β-HCH, trans-nonachlor, p,p’-DDT and p,p’-DDE. For the first inference, BRTs were used to determine the relative importance of predictor variables and possible interactions for each OCP.
The relative importance of a predictor variable in each BRT was calculated by weighting the number of times the variable was selected for a split by the squared improvement to the model caused by the split [36]. This improvement was determined by the least squares improvement criterion described in Friedman et al. [37]. The weighted value was averaged over all trees and scaled so that the sum of the relative importance of all predictor variables added up to 100 [36]. On this scale, larger values indicated stronger influence of the predictor variable on the response. The partial dependence plot of the categorical variables on the response illustrates the marginal effect of each fitted factor; a marginal effect greater than zero indicates higher concentrations than expected and values below zero indicate lower concentrations.
To evaluate the relative importance of interaction effects, each two-way interaction between predictor variables is related to BRT predictions using a linear model [36]. The residual variance of the linear model is then used as a measure of relative strength of the interaction fitted by the BRT [36].
All BRT models were fit in R v3.3.2 using the gbm package and brt extension [35, 36] with 50% of the data used for each training and test set. Results of the BRT analysis are detailed in the Results section. Important interactions identified by BRTs motivated the development of models for the second stage of the analysis, the main results of the BRT analysis are also referenced in the next section.
Hierarchical model
BRTs are well suited for modelling the influence of predictors and possible interactions but cannot be used to test differences between sub-groups or assess the effect sizes. Therefore, hierarchical models were used to infer the effects of predictors and interactions on OCP concentrations. Hierarchical linear models expand upon the linear modelling framework for the analysis of nested data structures [38].
Three hierarchical models are outlined in this section, labelled Models 1-3. Each model assumed a common response variable, yagt, defined as the observed OCP concentration for age group a (a = 1, …, 5), gender g (g = 1, 2) and time t (t = 1, …, 5). To reduce notation, each model is defined for a single OCP.
Model 1 was motivated by the interaction between age group and time, which was identified as the most influential interaction for β-HCH, Trans-nonachlor, p,p’-DDE and p,p’-DDT. This model took the following form for each yagt:
For HCB concentrations, the BRT analysis revealed a key interaction between age group and gender. This finding motivated a second hierarchical model to account for this interaction, labelled Model 2 and defined in Eq (2).
In Model 2, the random effect, αag was included to account for interactions between age group and gender. An intercept term for each age group, δa, was also included. Model 2 specifies that the changes in OCP concentrations can be described by the differences between age groups and genders and are not affected by the time.
Model 3 accounted for all interactions with age group, therefore combining BRT findings across all OCPs:
Note that in Model 2, αag was described by a normal distribution with a mean of δa and an overall variance, whereas in Model 3, it was described by a normal distribution with a mean of δg and an overall variance. As in Model 1, the priors for βa0, δa and δg were non-informative normal distributions with a mean of zero and large variance. All variances were also given non-informative priors.
The second inference focused on the interaction effects on OCP concentrations. The posterior means and 95% credible intervals of the OCP concentrations were calculated from the posterior distribution of the parameter estimates to compare interactions and visualise the temporal trends.
The third inference summarised the changes in OCP concentrations over time. The percent change and 95% credible interval of OCP concentrations between the years of sample collection was calculated from the posterior distribution of parameter estimates generated by the hierarchical model. For an overall view of the change, the percent change in OCP concentrations between 2002/03 and 2012/13 was calculated. The change between each year of sample collection was also calculated to distinguish any specific trend. To evaluate cohort effects, estimated OCP concentrations for each age group in 2002/03 were compared with their extrapolated age group in 2012/13.
The proposed hierarchical models were estimated by Markov Chain Monte Carlo (MCMC) using ‘JAGS’ in R version v3.3.2 [39–41]. Parameter estimates were based on two MCMC chains each comprising 10,000 iterations following a burn-in period of 10,000 iterations. The MCMC chains reached convergence, as assessed by the Gelman-Rubin diagnostic (S1 Fig).
Model assessment
The BRT model fit was graphically evaluated by comparing the observed OCP concentrations to the model’s predicted OCP concentrations. In a model with a good fit, the difference between the observed and predicted OCP concentrations would be small. On a graph with identical axes, the plotted values would be close to or on the diagonal line; larger differences would be plotted farther from the diagonal line.
The hierarchical models were assessed with two approaches, namely the predictive intervals (PI) [42] and deviance information criterion (DIC) [43].
For each hierarchical model, the (1 − α)% predictive intervals (PI) were determined for each observation, yi. In the study, three significance levels were considered, namely α = 0.2, 0.1 and 0.05. The proportion of observed values that fell into the respective PIs was calculated. For an adequate model, approximately (1 − α) of the observations should fall in their respective intervals [42].
The DIC [43] for each hierarchical model based on data y with parameters θ has the following form,
All models were also evaluated with root mean squared error (RMSE). RMSE is a measure of the differences between observed values and predicted values from a model. This measurement of goodness of fit was calculated by Eq 7.
Sensitivity assessment
Given the discrepancy in the number of samples for the 5-15 age group sampled in 2006/07 (Table 1), it was possible that the unbalanced nature of the data could distort the modelling results of the analysis. The sensitivity of the models with respect to the number of samples was therefore assessed as follows. Four pooled samples of the 5-15 age group in 2006/07 were randomly selected and analysed with the remaining age groups under Models 1, 2 and 3. Subsequent parameter estimates were compared with estimates obtained using the original number of samples.
Results
Model assessment
BRTs
Model assessment demonstrated larger RMSE values for β-HCH and p,p’-DDT concentrations compared to other OCPs, indicative of larger residuals for the former (Table 2). Graphical evaluation of the observed and predicted BRTs values indicated an acceptable goodness of fit for all OCPs (Fig 1). The predictions of β-HCH and p,p’-DDT concentrations from the BRT models were less accurate than those of HCB, Trans-nonachlor and p,p’-DDE concentrations.
Graphical comparison of observed (horizontal axis) and predicted (vertical axis) log-transformed OCP concentrations to illustrate accuracy of BRT models for HCB, β-HCH, Trans-nonachlor, p,p’-DDE and p,p’-DDT concentrations. BRT models for HCB and p,p’-DDE had more accurate predictions than models for β-HCH, Trans-nonachlor and p,p’-DDT. All points in these plots have partial transparency so that a darker shade is indicative of overlapping points. Unit of OCP concentrations were in ng/g lipid.
Hierarchical model
The best fitting model varied between OCPs (Table 3). The PIs obtained for Models 2 and 3 for Trans-nonachlor, p,p’-DDE and p,p’-DDT were similar. Based on this metric, Model 1 was best suited for β-HCH, Model 2 for HCB, Trans-nonachlor and p,p’-DDE, and Model 3 for p,p’-DDT. The DIC and RMSE values were smallest for Model 3 indicating that the most complex model had the best fit for all OCPs both with and without adjustment to the number of variables. The standard deviation of the posterior mean was also smallest for Model 3.
Variable importance
The relative importance of variables and factors was consistent for all OCPs while the relative contribution of interactions differed between OCPs.
Results of the BRT analysis (Fig 2) showed that the variable with the highest relative importance was age group followed by time and gender for all OCPs. The importance of age group in the analysis reflected the known persistence and accumulation of OCPs in the body [24, 44]. The relative importance scores for gender was lowest among the predictors. The tree analysis also described the relative importance of interactions between the variables (Fig 3). Age group x gender had the largest relative importance among interaction terms for HCB concentrations; age group x time was the most important interaction term for β-HCH, Trans-nonachlor, p,p’-DDE and p,p’ -DDT concentrations. These terms were later fed into the hierarchical model.
Relative importance of variables in predicting the concentrations of β-HCH, HCB, p,p’-DDE, p,p’-DDT and Trans-nonachlor based on a BRT analysis. In order of relative importance (left to right), the variables are age group, time and gender.
Relative importance of interactions in predicting the concentrations of β-HCH, HCB, p,p’-DDE, p,p’-DDT and Trans-nonachlor based on a BRT analysis. The interactions from left to right are between age group and time, age group and gender, and time and gender.
Factors that strongly influenced the concentration of HCB and p,p’-DDE included the older age groups 46-60 and >60 years, the earlier years 2002/03 and 2006/07 and gender (Fig 4). β-HCH and p,p’ -DDT concentrations were similarly impacted by children aged 0-4 years as well as adults aged 46-60 and >60 years, the years 2002/03, 2006/07 and 2008/09 and gender. The factors that influenced the concentration of Trans-nonachlor included adults aged 31-45, 46-60 and >60 years, the years 2002/03, 2006/07 and 2008/09 and gender.
Partial dependence plots of the marginal effects of age group (top), time (middle) and gender (bottom) on the concentrations of HCB, β-HCH, Trans-nonachlor, p,p’-DDE and p,p’-DDT where the zero-point of the marginal effect indicates the model expected OCP concentration. Values above zero indicate higher concentrations than expected and values below zero indicate lower concentrations.
Interaction effect on OCP concentrations
Hierarchical models demonstrated evidence of interaction effects between age group and time for β-HCH concentrations as well as age group and gender for HCB concentrations. Older age groups were also found to have higher HCB, Trans-nonchlor and p,p’-DDE concentrations than their younger counterparts.
Results of the hierarchical model are illustrated in Fig 5. There was strong evidence of an interaction effect between age group and time in β-HCH concentrations for adults aged >60 years; this age group had higher concentrations compared to children aged <17 years across all years of sample collection. The effect of an interaction between age group and gender was strongly supported in HCB concentrations for adults aged >45 years where women had higher concentrations than men; there was no evidence of the interaction for Trans-nonachlor and p.p’-DDE. The concentration of HCB in adults aged >60 years was higher than other age groups. The p,p’-DDE concentrations were higher in adults aged >60 years compared to 45-60 years and both age groups had higher concentrations than the younger age groups. The concentrations of Trans-nonachlor were higher in >60, 45-60 and 31-45 year olds compared to the younger age groups. There was no strong support for interaction effects between age group and time as well as age group and gender on p,p’-DDT concentrations.
Posterior mean log-transformed concentrations of HCB, β-HCH, Trans-nonachlor, p,p’-DDE and p,p’-DDT plotted over the years 2002/03 to 2012/13 for age groups 31-45, 46-60 and >60 year olds. HCB, Trans-nonachlor and p,p’-DDE were specified by Model 2, p,p’-DDT by Model 3 and β-HCH by Model 1. Line type indicates the gender, as described in the figure legend. The light grey bands around the posterior means indicate the 95% credible interval and the dark grey bands indicate an overlap of the 95% credible intervals where there is not enough evidence for the OCP concentration. Units of OCP concentrations were in ng/g lipid.
Results of the sensitivity assessment with respect to the age group 5-15 years sampled in 2006/07 showed minimal impact on the posterior mean estimates. Therefore the discrepancy in the number of pooled samples did not distort the modelling results of the analysis. Results can be found in the Supporting Information (S2 Fig).
Changes in OCP concentrations
All age groups for β-HCH and p,p’-DDT showed a decrease in concentrations between 2002/03 and 2012/13 (Fig 6). There was strong evidence of a decrease in β-HCH concentration between 2002/03 and 2012/13 for adults aged 45-60 years. There were no distinct changes in β-HCH and p,p’-DDT concentrations between each consecutive time of sample collection (Fig 7). There was no strong evidence of a change in β-HCH and p,p’-DDT concentrations between an age group in 2002/03 and its extrapolated age group in 2012/13 (Fig 8). As time was not specified in Model 2, HCB, Trans-nonachlor and p,p’-DDE were excluded from this analysis. The interpreted effect of age on the concentrations of HCB, Trans-nonachlor and p,p’-DDE could also be related to birth cohort effects, taking into account the historical use of OCPs prior to 1985.
The percent change (%) of β-HCH and p,p’-DDT concentrations (ng/g lipid) between 2002/03 and 2012/13 for age groups 5-15, 16-30, 31-45, 46-60 and >60. The error bars for each age group indicate the 95% credible interval of the percent change.
The percent change (%) of β-HCH and p,p’-DDT concentrations (ng/g lipid) between each consecutive sampling times for age groups 5-15, 16-30, 31-45, 46-60 and >60. For each age group, the sampling times listed from left to right indicate the differences between 2002/03 and 2006/07, 2006/07 and 2008/09, 2008/09 and 2010/11, and 2010/11 and 2012/13. The error bars for each age group indicate the 95% credible interval of the percent change.
The percent change (%) of β-HCH and p,p’-DDT concentrations (ng/g lipid) between an age group sampled in 2002/03 and its adjacent older age group sampled in 2012/13. The error bars for each age group indicate the 95% credible interval of the percent change. L-R: Difference 1: 5-15 year olds in 2002/03 and 16-30 year olds in 2012/13, Difference 2: 16-30 year olds in 2002/03 and 31-45 year olds in 2012/13, Difference 3: 31-45 year olds in 2002/03 and 45-60 year olds in 2012/13, Difference 4: 45-60 year olds in 2002/03 and >60 year olds in 2012/13.
Discussion
This paper has proposed a two-stage analysis for the exploration of age group, gender and temporal effects, and their interaction, on OCP concentrations. Using BRTs and hierarchical modelling, three statistical inferences were presented, using data on serum OCP concentrations. These inferences focused on the relative importance of variables and interactions, the effects of interactions on OCP concentrations and the change in OCP concentrations over time. The inferences highlighted strong evidence for an interaction effect between age group and time for β-HCH concentrations in adults aged >60 years as well as age group and gender for HCB concentrations in adults aged >45 years.
The application of BRTs for identifying key interactions was motivated by its underlying non-parametric approach that allowed the tree to model non-linear associations as well as the model’s robustness to outliers and missing data [32, 36, 45]. Information on the relative importance of all possible interactions, without needing to specify these terms prior to model implementation, is a key benefit of BRTs. Important interaction terms from the BRT were utilised in a hierarchical model which allowed multiple levels of the data to be included in a single model [46].
Previous analyses of OCP concentrations over multiple years have demonstrated a decline in OCP concentrations over time [18, 19, 47, 48]. A major cause of the declining trend has been attributed to the ban of the commercial or agricultural use of OCPs in various countries [18, 19].
In this paper, there was a similar negative trend for both β-HCH and p,p’-DDT concentrations, but the analysis showed no strong evidence that there was a substantive decrease in these OCP concentrations. This is not unexpected, since the study period spans years, long after major historical emissions [8] and the smaller ongoing emissions during the study period may result in lower and more diffuse exposure in humans, based on latent pathways into human serum [18, 49, 50].
In studies where OCP concentrations were sampled from various age groups, lower concentration of chemicals have been found in younger age groups [27, 51–53]. The same trend was observed in this study with the two-stage analysis providing strong evidence for larger OCP concentrations in older age groups. The difference in concentrations between age groups has been attributed to the implementation of bans on OCPs; older age groups had a longer exposure history to OCPs as they had lived through the agricultural use of the chemicals and hence had accumulated higher concentrations [54, 55] compared to their younger counterparts.
A number of studies had also found differences in OCP concentrations between genders; some studies found higher concentrations in women [26, 56, 57] while other studies found higher concentrations in men [27, 58, 59]. In the current study, the analysis showed higher concentrations of HCB, β-HCH, p,p’-DDE and p,p’-DDT in women, and strong evidence of this in HCB concentrations. Hypotheses on why there is a difference in concentrations between genders differ depending on age of the females. For parous (corresponding to our 15-30 and 31-45 year olds) females there is excretion due to placental transfer and breastfeeding [60]. In older females, this factor as discussed by Salihovic et al. (2012) is no longer relevant as the time since breastfeeding has long since past. While Salihovic et al. [59] suggest that females have higher body fat than males, Porta et al. [26] investigated OCPs and BMI and could not explain the higher concentrations in older females compared to males. Further work is required on gender differences in regard to OCP concentrations if there are likely to be any clinical repercussions.
Further to the main effects of age group, gender and time, most studies on OCP have not analysed the effect of interactions on the concentrations [21]. Thomas et al. [22] found that the interaction between age group and time had a significant effect on the HCB concentration and was best implemented for modelling β-HCH and p,p’-DDT. Thomas et al. [22] also found a significant interaction effect between age group and gender on p,p’-DDT concentrations. In this study, HCB, Trans-nonachlor and p,p’-DDE concentrations were best modelled by the interaction between age group and gender, without a temporal effect. This suggested that there was negligible change in OCP concentrations over the times of sample collection that were studied. This is possibly related to the indirect pathways of HCB, Trans-nonachlor and p,p’-DDE [18, 49, 50] into human serum long after its historical use [8]. Gender did not have a strong effect on the β-HCH concentrations which were modelled as a function of the interaction between age group and time. p,p’-DDT concentrations were modelled by interactions between age group and gender as well age group and time demonstrating the effect of age group, gender and time on the change in concentrations.
Limitations
Results did not provide a complete overview of temporal trends across all age groups as no data was collected for the years 2002/03 for 0-4 year olds, and 2004/05 for all age groups.
The analysis was unable to account for confounding factors such as metabolism rate, diet, place of residence and weight loss due to the pooled sampling method adopted to collect data.
Extensions
The hierarchical model utilised in this paper could be further extended to incorporate prior knowledge, elicited from experts or meta-analyses of published studies. It may be more appropriate for prior information to come from experts in the field as there are no standardised methods for laboratory analysis, detection limits, number of people sampled, age ranges, years sampled, number of years sampled and whether samples were pooled to lower laboratory costs [26, 61]. Prior information can be included as information on the coefficient values for the interactions between age group and gender, and age group and time.
In the current study, each OCP was modelled individually. This is useful as each chemical has a different history of production and use. However, OCPs could be modelled in a multivariate analysis that accounts for any positive correlations between OCPs such as that between p,p’-DDT and p,p’-DDE.
The model can also be applied to other POPs such as perfluorinated chemicals and polybrominated diphenyl ethers in order to investigate the effects of age group, gender and time as well as any interactions between these variables.
Supporting information
S1 Fig [tif]
Graphical summary of Gelman-Rubin diagnostic for Models 1, 2 and 3 of HCB, -HCH, Trans-nonachlor, ’-DDE and ’-DDT to demonstrate model convergence.
S2 Fig [tif]
Comparison of the posterior mean of Models 1, 2 and 3 for HCB, -HCH, Trans-nonachlor, ’-DDE and ’-DDT concentrations for males and females in the age group 5-15 sampled in 2006/07 using two datasets.
S1 File [pdf]
Details of sample collection and analysis.
Zdroje
1. Laug EP, Kunze FM, Prickett C. Occurrence of DDT in human fat and milk. AMA archives of industrial hygiene and occupational medicine. 1951;3(3):245–246. 14810251
2. Egan H, Goulding R, Roburn J, Tatton J. Organo-chlorine pesticide residues in human fat and human milk. BMJ. 1965;2(5453):66–69. doi: 10.1136/bmj.2.5453.66 14308202
3. Johnson-Restrepo B, Kannan K, Rapaport DP, Rodan BD. Polybrominated Diphenyl Ethers and Polychlorinated Biphenyls in Human Adipose Tissue from New York. Environmental Science & Technology. 2005;39(14):5177–5182. doi: 10.1021/es050399x
4. Kang JH, Son MH, Do Hur S, Hong S, Motoyama H, Fukui K, et al. Deposition of organochlorine pesticides into the surface snow of East Antarctica. Science of the total environment. 2012;433:290–295. doi: 10.1016/j.scitotenv.2012.06.037 22796727
5. Shen L, Wania F, Lei YD, Teixeira C, Muir DC, Bidleman TF. Atmospheric distribution and long-range transport behavior of organochlorine pesticides in North America. Environmental Science & Technology. 2005;39(2):409–420. doi: 10.1021/es049489c
6. Mrema EJ, Rubino FM, Brambilla G, Moretto A, Tsatsakis AM, Colosio C. Persistent organochlorinated pesticides and mechanisms of their toxicity. Toxicology. 2013;307:74–88. doi: 10.1016/j.tox.2012.11.015 23219589
7. Jung D, Becher H, Edler L, Flesch-Janys D, Gurn P, Konietzko J, et al. Elimination of β-hexachlorocyclohexane in occupationally exposed persons. Journal of toxicology and environmental health. 1997;51(1):23–34. doi: 10.1080/00984109708984009 9169059
8. Australian Academy of Science. The use of DDT in Australia. Reports of the Australian Academy of Science. 1972; No. 14. Canberra.
9. Australian Pesticides and Veterinary Medicines Authority. The history of ‘Organochlorine’ pesticides in Australia; 2013. http://www.apvma.gov.au/products/review/completed/organochlorines_history.php. Accessed 29/05/13.
10. Longnecker MP, Rogan WJ, Lucier G. The human health effects of DDT (dichlorodiphenyltrichloroethane) and PCBs (polychlorinated biphenyls) and an overview of organochlorines in public health 1. Annual review of public health. 1997;18(1):211–244. doi: 10.1146/annurev.publhealth.18.1.211 9143718
11. Beard J, Australian Rural Health Research Collaboration. DDT and human health. Science of the total environment. 2006;355(1):78–89. doi: 10.1016/j.scitotenv.2005.02.022 15894351
12. Arrebola J, Cuellar M, Claure E, Quevedo M, Antelo S, Mutch E, et al. Concentrations of organochlorine pesticides and polychlorinated biphenyls in human serum and adipose tissue from Bolivia. Environmental Research. 2012;112:40–47. doi: 10.1016/j.envres.2011.10.006 22078547
13. Waliszewski SM, Caba M, Herrero-Mercado M, Saldariaga-Noreña H, Meza E, Zepeda R, et al. Organochlorine pesticide residue levels in blood serum of inhabitants from Veracruz, Mexico. Environmental Monitoring and Assessment. 2012;184(9):5613–5621. doi: 10.1007/s10661-011-2366-2 21922174
14. Jakszyn P, Goñi F, Etxeandia A, Vives A, Millán E, López R, et al. Serum levels of organochlorine pesticides in healthy adults from five regions of Spain. Chemosphere. 2009;76(11):1518–1524. doi: 10.1016/j.chemosphere.2009.05.048 19586652
15. Jensen AA. Polychlorobiphenyls (PCBs), polychlorodibenzo-p-dioxins (PCDDs) and polychlorodibenzofurans (PCDFs) in human milk, blood and adipose tissue. Science of the total environment. 1987;64(3):259–293. doi: 10.1016/0048-9697(87)90250-6 3110947
16. Waliszewski S, Aguirre A, Infanzon R, Silva C, Siliceo J. Organochlorine pesticide levels in maternal adipose tissue, maternal blood serum, umbilical blood serum, and milk from inhabitants of Veracruz, Mexico. Archives of Environmental Contamination and Toxicology. 2001;40(3):432–438. doi: 10.1007/s002440010194 11443377
17. Toms LM, Harden F, Hobson P, Sjödin A, Mueller J. Temporal trend of organochlorine pesticides in Australia. Organohalogen Compounds. 2012;74: 775–778.
18. Mueller JF, Harden F, Toms LM, Symons R, Fürst P. Persistent organochlorine pesticides in human milk samples from Australia. Chemosphere. 2008;70(4):712–720. doi: 10.1016/j.chemosphere.2007.06.037 17675211
19. Solomon GM, Weiss PM. Chemical contaminants in breast milk: time trends and regional variability. Environmental Health Perspectives. 2002;110(6):A339–A347. doi: 10.1289/ehp.021100339 12055065
20. Polder A, Gabrielsen G, Odland J, Savinova T, Tkachev A, Løken K, et al. Spatial and temporal changes of chlorinated pesticides, PCBs, dioxins (PCDDs/PCDFs) and brominated flame retardants in human breast milk from Northern Russia. Science of the total environment. 2008;391(1):41–54. doi: 10.1016/j.scitotenv.2007.10.045 18063018
21. Porta M, Puigdomènech E, Ballester F, Selva J, Ribas-Fitó N, Llop S, et al. Monitoring concentrations of persistent organic pollutants in the general population: the international experience. Environment international. 2008;34(4):546–561. doi: 10.1016/j.envint.2007.10.004 18054079
22. Thomas A, Toms LML, Harden FA, Hobson P, White NM, Mengersen KL, et al. Concentrations of organochlorine pesticides in pooled human serum by age and gender. Environmental Research. 2017;154:10–18. doi: 10.1016/j.envres.2016.12.009 27992738
23. Wolff MS, Anderson HA, Britton JA, Rothman N. Pharmacokinetic variability and modern epidemiology-the example of dichlorodiphenyltrichloroethane, body mass index, and birth cohort. Cancer Epidemiology Biomarkers & Prevention. 2007;16(10):1925–1930. doi: 10.1158/1055-9965.EPI-07-0394
24. Laden F, Neas LM, Spiegelman D, Hankinson SE, Willett WC, Ireland K, et al. Predictors of plasma concentrations of DDE and PCBs in a group of US women. Environmental health perspectives. 1999;107(1):75–81. doi: 10.1289/ehp.9910775 9872720
25. James RA, Hertz-Picciotto I, Willman E, Keller JA, Charles MJ. Determinants of serum polychlorinated biphenyls and organochlorine pesticides measured in women from the child health and development study cohort, 1963-1967. Environmental health perspectives. 2002;110(7):617–624. doi: 10.1289/ehp.02110617 12117636
26. Porta M, López T, Gasull M, Rodríguez-Sanz M, Garí M, Pumarega J, et al. Distribution of blood concentrations of persistent organic pollutants in a representative sample of the population of Barcelona in 2006, and comparison with levels in 2002. Science of the total environment. 2012;423:151–161. doi: 10.1016/j.scitotenv.2012.02.001 22397902
27. Wang HS, Chen ZJ, Wei W, Man YB, Giesy JP, Du J, et al. Concentrations of organochlorine pesticides (OCPs) in human blood plasma from Hong Kong: Markers of exposure and sources from fish. Environment international. 2013;54:18–25. doi: 10.1016/j.envint.2013.01.003 23376599
28. Woodruff TJ, Zota AR, Schwartz JM. Environmental chemicals in pregnant women in the United States: NHANES 2003-2004. Environmental Health Perspectives. 2011;119(6):878–885. doi: 10.1289/ehp.1002727 21233055
29. Toms LM, Thompson J, Rotander A, Hobson P, Calafat AM, Kato K, et al. Decline in perfluorooctane sulfonate and perfluorooctanoate serum concentrations in an Australian population from 2002 to 2011. Environment international. 2014;71:74–80. doi: 10.1016/j.envint.2014.05.019 24980755
30. Frugé AD, Cases MG, Schildkraut JM, Demark-Wahnefried W. Associations between obesity, body fat distribution, weight loss and weight cycling on serum pesticide concentrations. Journal of food & nutritional disorders. 2016;5(3).
31. Breiman L, Friedman J, Stone CJ, Olshen RA. Classification and regression trees. Boca Raton, Florida: Chapman & Hall/CRC; 1984.
32. Schapire RE. In: Denison DD, Hansen MH, Holmes CC, Mallick B, Yu B, editors. The Boosting Approach to Machine Learning: An Overview. New York, NY: Springer; 2003. p. 149–171.
33. Zhang H, Singer BH. Recursive partitioning and applications. 2nd ed. New York, NY: Springer; 2010.
34. Friedman J, Hastie T, Tibshirani R. Additive logistic regression: a statistical view of boosting. The annals of statistics. 2000;28(2):337–407. doi: 10.1214/aos/1016218222
35. Greenwell B, Boehmke B, Cunningham J, GBM Developers. gbm: Generalized Boosted Regression Models; 2019. R package version 2.1.5. Available from: https://CRAN.R-project.org/package=gbm.
36. Elith J, Leathwick JR, Hastie T. A working guide to boosted regression trees. Journal of Animal Ecology. 2008;77(4):802–813. doi: 10.1111/j.1365-2656.2008.01390.x 18397250
37. Friedman JH. Greedy function approximation: a gradient boosting machine. Annals of statistics. 2001;29(5):1189–1232. doi: 10.1214/aos/1013203450
38. Walsh CD, Mengersen KL. Ordering of hierarchies in hierarchical models: Bone mineral density estimation. In: Alston C, Mengersen K, Pettitt A, editors. Case Studies in Bayesian Statistical Modelling and Analysis. Chichester, West Sussex: John Wiley & Sons, Ltd; 2012. p. 159–170.
39. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria; 2016. Available from: https://www.R-project.org/.
40. Plummer M, Best N, Cowles K, Vines K. CODA: Convergence Diagnosis and Output Analysis for MCMC. R News. 2006;6(1):7–11. Available from: https://journal.r-project.org/archive/.
41. Plummer M. rjags: Bayesian Graphical Models using MCMC; 2019. R package version 4-9. Available from: https://CRAN.R-project.org/package=rjags.
42. Gelman A, Meng XL, Stern H. Posterior predictive assessment of model fitness via realized discrepancies. Statistica sinica. 1996;6:733–807.
43. Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A. Bayesian measures of model complexity and fit. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2002;64(4):583–639. doi: 10.1111/1467-9868.00353
44. Wolff MS, Deych E, Ojo F, Berkowitz GS. Predictors of organochlorines in New York City pregnant women, 1998-2001. Environmental Research. 2005;97(2):170–177. doi: 10.1016/j.envres.2004.07.014 15533333
45. Leathwick J, Elith J, Francis M, Hastie T, Taylor P. Variation in demersal fish species richness in the oceans surrounding New Zealand: an analysis using boosted regression trees. Marine Ecology Progress Series. 2006;321:267–281. doi: 10.3354/meps321267
46. Gelman A, Hill J, Yajima M. Why we (usually) don’t have to worry about multiple comparisons. Journal of Research on Educational Effectiveness. 2012;5(2):189–211. doi: 10.1080/19345747.2011.618213
47. Norén K, Meironyté D. Certain organochlorine and organobromine contaminants in Swedish human milk in perspective of past 20–30 years. Chemosphere. 2000;40(9):1111–1123. doi: 10.1016/s0045-6535(99)00360-4 10739053
48. Mannetje A, Coakley J, Bridgen P, Brooks C, Harrad S, Smith AH, et al. Current concentrations, temporal trends and determinants of persistent organic pollutants in breast milk of New Zealand women. Science of the Total Environment. 2013;458:399–407. doi: 10.1016/j.scitotenv.2013.04.055 23685364
49. Dewan P, Jain V, Gupta P, Banerjee BD. Organochlorine Pesticide Residues In Maternal Blood, Cord Blood, Placenta, And Breastmilk And Their Relation To Birth Size. Chemosphere. 2013;90(5):1704–1710. doi: 10.1016/j.chemosphere.2012.09.083 23141556
50. Siddiqui M, Saxena M, Bhargava A, Seth T, Murti CK, Kutty D. Agrochemicals In The Maternal Blood, Milk, And Cord Blood: A Source Of Toxicants For Prenates And Neonates. Environmental Research. 1981;24(1):24–32. doi: 10.1016/0013-9351(81)90128-6 6163625
51. Meng G, Feng Y, Nie Z, Wu X, Wei H, Wu S, et al. Internal exposure levels of typical POPs and their associations with childhood asthma in Shanghai, China. Environmental research. 2016;146:125–135. doi: 10.1016/j.envres.2015.12.026 26748225
52. Schade G, Heinzow B. Organochlorine pesticides and polychlorinated biphenyls in human milk of mothers living in northern Germany: current extent of contamination, time trend from 1986 to 1997 and factors that influence the levels of contamination. Science of the total environment. 1998;215(1):31–39. doi: 10.1016/s0048-9697(98)00008-4 9599454
53. Bjermo H, Darnerud PO, Lignell S, Pearson M, Rantakokko P, Nälsén C, et al. Fish intake and breastfeeding time are associated with serum concentrations of organochlorines in a Swedish population. Environment International. 2013;51:88–96. doi: 10.1016/j.envint.2012.10.010 23201820
54. Saoudi A, Fréry N, Zeghnoun A, Bidondo ML, Deschamps V, Göen T, et al. Serum levels of organochlorine pesticides in the French adult population: the French National Nutrition and Health Study (ENNS), 2006–2007. Science of the Total Environment. 2014;472:1089–1099. doi: 10.1016/j.scitotenv.2013.11.044 24361744
55. Luo D, Zhou T, Tao Y, Feng Y, Shen X, Mei S. Exposure to organochlorine pesticides and non-Hodgkin lymphoma: a meta-analysis of observational studies. Scientific reports. 2016;6:25768. doi: 10.1038/srep25768 27185567
56. Malarvannan G, Kunisue T, Isobe T, Sudaryanto A, Takahashi S, Prudente M, et al. Organohalogen compounds in human breast milk from mothers living in Payatas and Malate, the Philippines: levels, accumulation kinetics and infant health risk. Environmental Pollution. 2009;157(6):1924–1932. doi: 10.1016/j.envpol.2009.01.010 19232797
57. Tue NM, Sudaryanto A, Minh TB, Nhat BH, Isobe T, Takahashi S, et al. Kinetic differences of legacy organochlorine pesticides and polychlorinated biphenyls in Vietnamese human breast milk. Chemosphere. 2010;81(8):1006–1011. doi: 10.1016/j.chemosphere.2010.09.013 20870263
58. Hardell E, Carlberg M, Nordström M, van Bavel B. Time trends of persistent organic pollutants in Sweden during 1993–2007 and relation to age, gender, body mass index, breast-feeding and parity. Science of the total environment. 2010;408(20):4412–4419. doi: 10.1016/j.scitotenv.2010.06.029 20643475
59. Salihovic S, Lampa E, Lindström G, Lind L, Lind PM, van Bavel B. Circulating levels of persistent organic pollutants (POPs) among elderly men and women from Sweden: results from the Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS). Environment International. 2012;44:59–67. doi: 10.1016/j.envint.2012.01.011 22361238
60. LaKind JS, Berlin CM Jr, Sjödin A, Turner W, Wang RY, Needham LL, et al. Do Human Milk Concentrations Of Persistent Organic Chemicals Really Decline During Lactation? Chemical Concentrations During Lactation And Milk/serum Partitioning. Environmental Health Perspectives. 2009;117(10):1625–1631. doi: 10.1289/ehp.0900876 20019916
61. Kang JH, Park H, Chang YS, Choi JW. Distribution of organochlorine pesticides (OCPs) and polychlorinated biphenyls (PCBs) in human serum from urban areas in Korea. Chemosphere. 2008;73(10):1625–1631. doi: 10.1016/j.chemosphere.2008.07.087 18829066
Článek vyšel v časopise
PLOS One
2019 Číslo 11
- Případ Bulovka: Jak postupují jiné nemocnice, aby podobným tragédiím zabránily?
- Jak často se u masových střelců vyskytují duševní choroby?
- FDA varuje před selfmonitoringem cukru pomocí chytrých hodinek. Jak je to v Česku?
- Metamizol jako analgetikum první volby: kdy, pro koho, jak a proč?
- Není statin jako statin aneb praktický přehled rozdílů jednotlivých molekul
Nejčtenější v tomto čísle
- A daily diary study on maladaptive daydreaming, mind wandering, and sleep disturbances: Examining within-person and between-persons relations
- Molecular validation of clinical Pantoea isolates identified by MALDI-TOF
- Pathways to conspiracy: The social and linguistic precursors of involvement in Reddit’s conspiracy theory forum
- Urticating setae of tarantulas (Araneae: Theraphosidae): Morphology, revision of typology and terminology and implications for taxonomy