Analysis of real-time mixture cytotoxicity data following repeated exposure using BK/TD models.

Cosmetic products generally consist of multiple ingredients. Thus, cosmetic risk assessment has to deal with mixture toxicity on a long-term scale which means it has to be assessed in the context of repeated exposure. Given that animal testing has been banned for cosmetics risk assessment, in vitro assays allowing long-term repeated exposure and adapted for in vitro - in vivo extrapolation need to be developed. However, most in vitro tests only assess short-term effects and consider static endpoints which hinder extrapolation to realistic human exposure scenarios where concentration in target organs is varies over time. Thanks to impedance metrics, real-time cell viability monitoring for repeated exposure has become possible. We recently constructed biokinetic/toxicodynamic models (BK/TD) to analyze such data (Teng et al., 2015) for three hepatotoxic cosmetic ingredients: coumarin, isoeugenol and benzophenone-2. In the present study, we aim to apply these models to analyze the dynamics of mixture impedance data using the concepts of concentration addition and independent action. Metabolic interactions between the mixture components were investigated, characterized and implemented in the models, as they impacted the actual cellular exposure. Indeed, cellular metabolism following mixture exposure induced a quick disappearance of the compounds from the exposure system. We showed that isoeugenol substantially decreased the metabolism of benzophenone-2, reducing the disappearance of this compound and enhancing its in vitro toxicity. Apart from this metabolic interaction, no mixtures showed any interaction, and all binary mixtures were successfully modeled by at least one model based on exposure to the individual compounds.

Perspectives for integrating human and environmental exposure assessments.

Integrated Risk Assessment (IRA) has been defined by the EU FP7 HEROIC Coordination action as "the mutual exploitation of Environmental Risk Assessment for Human Health Risk Assessment and vice versa in order to coherently and more efficiently characterize an overall risk to humans and the environment for better informing the risk analysis process" (Wilks et al., 2015). Since exposure assessment and hazard characterization are the pillars of risk assessment, integrating Environmental Exposure assessment (EEA) and Human Exposure assessment (HEA) is a major component of an IRA framework. EEA and HEA typically pursue different targets, protection goals and timeframe. However, human and wildlife species also share the same environment and they similarly inhale air and ingest water and food through often similar overlapping pathways of exposure. Fate models used in EEA and HEA to predict the chemicals distribution among physical and biological media are essentially based on common properties of chemicals, and internal concentration estimations are largely based on inter-species (i.e. biota-to-human) extrapolations. Also, both EEA and HEA are challenged by increasing scientific complexity and resources constraints. Altogether, these points create the need for a better exploitation of all currently existing data, experimental approaches and modeling tools and it is assumed that a more integrated approach of both EEA and HEA may be part of the solution. Based on the outcome of an Expert Workshop on Extrapolations in Integrated Exposure Assessment organized by the HEROIC project in January 2014, this paper identifies perspectives and recommendations to better harmonize and extrapolate exposure assessment data, models and methods between Human Health and Environmental Risk Assessments to support the further development and promotion of the concept of IRA. Ultimately, these recommendations may feed into guidance showing when and how to apply IRA in the regulatory decision-making process for chemicals.

BK/TD models for analyzing in vitro impedance data on cytotoxicity.

The ban of animal testing has enhanced the development of new in vitro technologies for cosmetics safety assessment. Impedance metrics is one such technology which enables monitoring of cell viability in real time. However, analyzing real time data requires moving from static to dynamic toxicity assessment. In the present study, we built mechanistic biokinetic/toxicodynamic (BK/TD) models to analyze the time course of cell viability in cytotoxicity assay using impedance. These models account for the fate of the tested compounds during the assay. BK/TD models were applied to analyze HepaRG cell viability, after single (48 h) and repeated (4 weeks) exposures to three hepatotoxic compounds (coumarin, isoeugenol and benzophenone-2). The BK/TD models properly fit the data used for their calibration that was obtained for single or repeated exposure. Only for one out of the three compounds, the models calibrated with a single exposure were able to predict repeated exposure data. We therefore recommend the use of long-term exposure in vitro data in order to adequately account for chronic hepatotoxic effects. The models we propose here are capable of being coupled with human biokinetic models in order to relate dose exposure and human hepatotoxicity.

Modelling the binding affinity of steroids to zebrafish sex hormone-binding globulin.

The circulating endogenous steroids are transported in the bloodstream. These are bound to a highly specific sex hormone-binding globulin (SHBG) and in lower affinity to proteins such as the corticosteroid-binding protein and albumin in vertebrates, including fish. It is generally believed that the glycoprotein SHBG protects these steroids from rapid metabolic degradation and thus intervenes in its availability at the target tissues. Endocrine disrupters binding to SHBG affect the normal activity of natural steroids. Since xenobiotics are primarily released in the aquatic environment, there is a need to evaluate the binding affinity of xenosteroid mimics on fish SHBG, especially in zebrafish (Danio rerio), a small freshwater fish originating in India and widely employed in ecotoxicology, toxicology, and genetics. In this context, a zebrafish SHBG (zfSHBG) homology model was developed using the human SHBG (hSHBG) receptor structure as template. It was shown that interactions with amino acids Ser-36, Asp-59 and Thr-54 were important for binding affinity. A ligand-based pharmacophore model was also developed for both zfSHBG and hSHBG inhibitors that differentiated binders from non-binders, but also demonstrated structural requirements for zfSHBG and hSHBG ligands. The study provides insights into the mechanism of action of endocrine disruptors in zebrafish as well as providing a useful tool for identifying anthropogenic compounds inhibiting zfSHBG.

Perspectives for integrating human and environmental risk assessment and synergies with socio-economic analysis.

For more than a decade, the integration of human and environmental risk assessment (RA) has become an attractive vision. At the same time, existing European regulations of chemical substances such as REACH (EC Regulation No. 1907/2006), the Plant Protection Products Regulation (EC regulation 1107/2009) and Biocide Regulation (EC Regulation 528/2012) continue to ask for sector-specific RAs, each of which have their individual information requirements regarding exposure and hazard data, and also use different methodologies for the ultimate risk quantification. In response to this difference between the vision for integration and the current scientific and regulatory practice, the present paper outlines five medium-term opportunities for integrating human and environmental RA, followed by detailed discussions of the associated major components and their state of the art. Current hazard assessment approaches are analyzed in terms of data availability and quality, and covering non-test tools, the integrated testing strategy (ITS) approach, the adverse outcome pathway (AOP) concept, methods for assessing uncertainty, and the issue of explicitly treating mixture toxicity. With respect to exposure, opportunities for integrating exposure assessment are discussed, taking into account the uncertainty, standardization and validation of exposure modeling as well as the availability of exposure data. A further focus is on ways to complement RA by a socio-economic assessment (SEA) in order to better inform about risk management options. In this way, the present analysis, developed as part of the EU FP7 project HEROIC, may contribute to paving the way for integrating, where useful and possible, human and environmental RA in a manner suitable for its coupling with SEA.

Exposure assessment of phthalates in French pregnant women: results of the ELFE pilot study.

The ubiquitous use of phthalate esters in plastics, building material, medical devices, personal care products and food packaging materials results in a widespread exposure of general population. This study reports measurement of urinary concentration of phthalate metabolites in France and provides a first assessment of the exposure of French pregnant women to this chemical class. For the majority of the phthalate metabolites, concentrations measured in urine were similar to those reported in previous studies except for two phthalates that were characterized by high concentrations of metabolites if compared to previous European and American studies: DiNP (Di-iso-nonylphthalate) and DEHP (Di(2-ethylhexyl)phthalate). In a second part of the study, a pharmacokinetic model was used in order to gain understanding on exposure to DEHP. A high concentration of the primary metabolite of DEHP, MEHP (Mono(2-ethylhexyl)phthalate), was thus identified probably because of a very recent exposure to perfusion materials at the hospital. Pharmacokinetics modelling highlighted that gathering data on the time gap between exposure and biomonitoring is an essential information requirement for reconstructing the dose of non persistent pollutants. Information about exposure pathway is also crucial for conducting effective reverse dosimetry.

Fluoxetine effects assessment on the life cycle of aquatic invertebrates.

Fluoxetine is a serotonin re-uptake inhibitor, generally used as an antidepressant. It is suspected to provoke substantial effects in the aquatic environment. This study reports the effects of fluoxetine on the life cycle of four invertebrate species, Daphnia magna, Hyalella azteca and the snail Potamopyrgus antipodarum exposed to fluoxetine spiked-water and the midge Chironomus riparius exposed to fluoxetine-spiked sediments. For D. magna, a multi-generational study was performed with exposition of newborns from exposed organisms. Effects of fluoxetine could be found at low measured concentrations (around 10microgl(-1)), especially for parthenogenetic reproduction of D. magna and P. antipodarum. For daphnids, newborns length was impacted by fluoxetine and the second generation of exposed individuals showed much more pronounced effects than the first one, with a NOEC of 8.9microgl(-1). For P. antipodarum, significant decrease of reproduction was found for concentrations around 10microgl(-1). In contrast, we found no effect on the reproduction of H. azteca but a significant effect on growth, which resulted in a NOEC of 33microgl(-1), expressed in nominal concentration. No effect on C. riparius could be found for measured concentrations up to 59.5mgkg(-1). General mechanistic energy-based models showed poor relevance for data analysis, which suggests that fluoxetine targets specific mechanisms of reproduction.

Do differences between metal body residues reflect the differences between effects for Chironomus riparius exposed to different sediments?

Sediment characteristics are well known to interfere with toxicity, mainly through differences in terms of bioaccumulation. Here, with chironomids exposed to zinc in an artificial and a field sediment, we investigated the differences of zinc accumulation and of effects on the life cycle, at individual and population level. We used biology and energy-based modeling to analyze the data at all the levels of biological organization. This permits a reliable estimation of thresholds values for tissue residues. Differences in zinc tissue residues accounted for most of the differences between the results for both sediments (a factor of 11 for differences from 20 to 100 depending on the parameter which is considered). Taking into account accumulation and background variability, the differences relative to thresholds could be accounted for. However, it appeared that, once the threshold was passed, effects were much more pronounced for organisms exposed to artificial sediment compared to field sediment. We concluded that some sediment characteristics can enhance toxicity, in addition to their influence on the compound accumulation, even if the latter was the major source of differences in our study.

Energy-based modeling to study population growth rate and production for the midge Chironomus riparius in ecotoxicological risk assessment.

Effects of toxicants are commonly assessed at individual level, whereas the aim of ecotoxicology is to protect ecosystems. We recently built energy-based models to describe and predict growth, emergence and reproduction of the midge Chironomus riparius [Péry (2002) Environ. Toxicol. Chem. 21, 2507-13]. Here we use these models to derive effects at the population level from effects at individual levels. The first endpoint we consider is population growth rate, which is a commonly studied endpoint at the population level. This parameter is informative relative to the risk of disappearance of the population. We also examined the production of organisms per generation or per time unity at population equilibrium, the study of which is allowed by our models and accounts for energy transfers. Such a study is crucial to predict effects on ecosystems, for species of the Chironomidae are keystone species, being the main food source of many other species, including birds and fish. We show in our study that the disappearance of the population can only occur in cases of very severe toxicity (99% decrease of reproduction, more than 97% mortality during a 10 days survival test or a difference of mean emergence times between males and females of more than 10 days). Concerning production of organisms, we show that reproduction decrease has little effect on it, that mortality of young larvae has an impact that cannot be neglected and that mortality of old larvae and delay of emergence has a strong effect. Our study suggests that bioassays should focus on an EC50 for reproduction, a LC30 for young instars, and a No Effect Concentration for old instars (growth and mortality) to prevent effects at the population level.

Survival tests with Chironomus riparius exposed to spiked sediments can profit from DEBtox model.

DEBtox model is a biologically based model used to analyse aquatic toxicity data (The analysis of aquatic toxicity data, VU University Press, Amsterdam, The Netherlands, 1996, 149pp.). To date, it has not been used to analyse sediment survival tests, mainly because this would require daily counting of survivors, which is difficult for benthic organisms. In the present study, we adapted survival toxicity tests with the midge Chironomus riparius to permit survival data analysis with DEBtox. To validate the adaptation proposed, we exposed organisms to two chemicals, copper and methiocarb. We exposed larvae of second, third and fourth instar to different concentrations of the tested compounds and two different diets to assess the influence of diet and instar on DEBtox parameter estimates. Daily counting of organisms did not affect survival and did not lead to much more effort when compared to usual tests. Moreover, the analysis profited much from use of the DEBtox model. It was possible, with only survival data to estimate and validate the kinetics of the compounds, to predict survival during the recovery of exposed organisms and to assess the influence of food availability on toxicity. Food availability did not influence methiocarb toxicity, but copper was more toxic to food-limited organisms when a threshold concentration was exceeded. Comparison of parameters estimate also helped in understanding the differences in toxicity responses between instars. For the two compounds we studied, the difference was entirely explained by differences in threshold values.