Toward a harmonized methodology to analyze field side effects of two pesticide products on earthworms at the EU level.

Before plant protection product (PPP) marketing authorization, a risk assessment for nontarget soil organisms (e.g., earthworms) is required as part of Regulation (EC) No. 1107/2009. Following a stepwise approach, higher tier earthworm field studies are needed if they cannot demonstrate low long-term risk based on laboratory studies. The European guidance for terrestrial ecotoxicology refers to ISO guideline 11268-3 as a standard to conduct earthworm field studies. Assessment of such studies may be challenging, as no European harmonized guidance is available to properly analyze the accuracy, representativeness, and appropriateness of experimental designs, as well as the statistical analysis robustness of results and their scientific reliability. Following the ISO guideline 11268-3, a field study was performed in 2016-2017 (Versailles, France). An assessment of the first year of this field study was performed in agreement with the quality criteria provided in 2006 in the guidance document published by de Jong and collaborators and recommendations by Kula and collaborators that allows describing the protocol and results of earthworm field studies. Not only did we underline the importance of a detailed analysis of raw data on the effects of pesticides on earthworms in field situations, but we also provided recommendations to harmonize protocols for assessing higher tier field studies devoted to earthworms to advance a better assessment of PPP fate and ecotoxicity. In particular, we provided practical field observations related to the study design, pesticide applications, and earthworm sampling. Concurrently, in addition to the conventional earthworm community study, we propose carrying out an assessment of soil function (i.e., organic matter decomposition, soil structuration, etc.) and calculating diversity indices to obtain information about earthworm community dynamics after the application of PPPs. Finally, through field observations, any relevant observation of external and/or internal recovery should be reported. Integr Environ Assess Manag 2023;19:254-271. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Soil dissipation and bioavailability to earthworms of two fungicides under laboratory and field conditions.

The representativeness of laboratory studies of the fate of pesticides in soil in field conditions is questionable. This study aimed at comparing the dissipation and bioavailability to earthworms of two fungicides, dimoxystrobin (DMX) and epoxiconazole (EPX), over 12 months under laboratory and field conditions. In both approaches, the fungicides were applied to the same loamy soil as a formulated mixture at several concentrations. We determined total DMX and EPX concentrations in the soil using exhaustive extraction, their environmental availability using mild extraction and their bioavailability through internal concentrations in exposed earthworms. The initial fungicide application appeared as much better controlled in terms of dose and homogeneity in the laboratory than in the field. One year after application, a similar dissipation rate was observed between the laboratory and field experiments (ca 80% and 60% for DMX and EPX, respectively). Similarly, the ratio of available/total concentrations in soil displayed the same trend whatever the duration and the conditions (field or lab), EPX being more available than DMX. Finally, the environmental bioavailability of the two fungicides to earthworms was heterogeneous in the field, but, in the laboratory, the bioaccumulation was evidenced to be dose-dependent only for DMX. Our findings suggest that the actual fate of the two considered fungicides in the environment is consistent with the one determined in the laboratory, although the conditions differed (e.g., presence of vegetation, endogeic earthworm species). This study allowed better understanding of the fate of the two considered active substances in the soil and underlined the need for more research dedicated to the link between environmental and toxicological bioavailability.

A two years field experiment to assess the impact of two fungicides on earthworm communities and their recovery.

Recent EFSA (European Food Safety Authority) reports highlighted that the ecological risk assessment of pesticides needed to go further by taking more into account the impacts of chemicals on biodiversity under field conditions. We assessed the effects of two commercial formulations of fungicides separately and in mixture, i.e., Cuprafor Micro® (containing 500 g kg-1 copper oxychloride) at 4 (C1, corresponding to 3.1 mg kg-1 dry soil of copper) and 40 kg ha-1 (C10), and Swing® Gold (50 g L-1 epoxiconazole EPX and 133 g L-1 dimoxystrobin DMX) at one (D1, 5.81 10-2 and 1.55 10-1 mg kg-1 dry soil of EPX and DMX, respectively) and ten times (D10) the recommended field rate, on earthworms at 1, 6, 12, 18 and 24 months after the application following the international ISO standard no. 11268-3 to determine the effects on earthworms in field situations. The D10 treatment significantly reduced the species diversity (Shannon diversity index, 54% of the control), anecic abundance (29% of the control), and total biomass (49% of the control) over the first 18 months of experiment. The Shannon diversity index also decreased in the mixture treatment (both fungicides at the recommended dose) at 1 and 6 months after the first application (68% of the control at both sampling dates), and in C10 (78% of the control) at 18 months compared with the control. Lumbricus terrestris, Aporrectodea caliginosa, Aporrectodea giardi, Aporrectodea longa, and Allolobophora chlorotica were (in decreasing order) the most sensitive species to the tested fungicides. This study not only addressed field ecotoxicological effects of fungicides at the community level and ecological recovery, but it also pinpointed some methodological weaknesses (e.g., regarding fungicide concentrations in soil and statistics) of the guideline to determine the effects on earthworms in field situations.

An energy-based model to analyze growth data of earthworms exposed to two fungicides.

The pesticide risk assessment for earthworms is currently performed using standardized tests, the model species Eisenia fetida, and the analyses of the data obtained are performed with ad hoc statistical tools. We assessed the impact of two fungicides on the entire growth pattern of the earthworm species Aporrectodea caliginosa, which is highly representative of agricultural fields. Individuals of three different ages (from hatching to 56 days old) were exposed to Cuprafor micro® (copper oxychloride) and Swing® Gold (dimoxystrobin and epoxiconazole). Data were analyzed with an energy-based toxicodynamic model coupled with a toxicokinetic model. The copper fungicide caused a drastic growth inhibition once the no effect concentration (NEC), estimated at 65 mg kg-1 of copper, was exceeded. The Swing® Gold negatively affected the growth with NEC values estimated at 0.387 mg kg-1 and 0.128 mg kg-1 for the dimoxystrobin and the epoxiconazole in this fungicide formulation, respectively. The time-profile of the effects on A. caliginosa individuals was fully accounted for by the model, whatever their age of exposure. Furthermore, toxicity data analyses, supported by measurements of fungicide concentrations in earthworm at the end of the experiment, allowed bettering understanding of the mechanisms of action of the fungicides towards earthworm growth.

Earthworms Mitigate Pesticide Effects on Soil Microbial Activities.

Earthworms act synergistically with microorganisms in soils. They are ecosystem engineers involved in soil organic matter degradation and nutrient cycling, leading to the modulation of resource availability for all soil organisms. Using a soil microcosm approach, we aimed to assess the influence of the earthworm Aporrectodea caliginosa on the response of soil microbial activities against two fungicides, i.e., Cuprafor Micro® (copper oxychloride, a metal) and Swing® Gold (epoxiconazole and dimoxystrobin, synthetic organic compounds). The potential nitrification activity (PNA) and soil enzyme activities (glucosidase, phosphatase, arylamidase, and urease) involved in biogeochemical cycling were measured at the end of the incubation period, together with earthworm biomass. Two common indices of the soil biochemistry were used to aggregate the response of the soil microbial functioning: the geometric mean (Gmean) and the Soil Quality Index (SQI). At the end of the experiment, the earthworm biomass was not impacted by the fungicide treatments. Overall, in the earthworm-free soil microcosms, the two fungicides significantly increased several soil enzyme and nitrification activities, leading to a higher GMean index as compared to the non-treated control soils. The microbial activity responses depended on the type of activity (nitrification was the most sensitive one), on the fungicide (Swing® Gold or Cuprafor Micro®), and on the doses. The SQI indices revealed higher effects of both fungicides on the soil microbial activity in the absence of earthworms. The presence of earthworms enhanced all soil microbial activities in both the control and fungicide-contaminated soils. Moreover, the magnitude of the fungicide impact, integrated through the SQI index, was mitigated by the presence of earthworms, conferring a higher stability of microbial functional diversity. Our results highlight the importance of biotic interactions in the response of indicators of soil functioning (i.e., microbial activity) to pesticides.

Effects of two common fungicides on the reproduction of Aporrectodea caliginosa in natural soil.

The use of pesticides in agroecosystems can have negative effects on earthworms, which play key roles in soil functioning such as organic matter decomposition. The aim of this study was to assess the effects of two fungicides (Cuprafor micro®, composed of copper oxychloride, and Swing Gold®, composed of epoxiconazole (EPX) and dimoxystrobin (DMX)) on earthworm reproduction by exposing adults and cocoons. First, adult Aporrectodea caliginosa individuals were exposed for 28 days to 3.33, 10 and 30 times the recommended dose (RD) of Cuprafor micro® corresponding to 25.8, 77.5 and 232.5 mg kg-1 dry soil of copper, respectively, and 0.33, 1 and 3 times the RD of Swing Gold® (corresponding to 5.2 × 10-2 mg DMX kg-1 + 1.94 × 10-2 mg EPX kg-1, 1.55 × 10-1 mg DMX kg-1 + 5.81 × 10-2 mg EPX kg-1 and 4.62 × 10-1 mg DMX kg-1 + 1.74 × 10-1 mg EPX kg-1 respectively), in addition to a control soil with no fungicide treatment. Cocoon variables (production, weight, hatching success, hatching time) were monitored. Second, "naïve" cocoons produced by uncontaminated earthworms were exposed to soils contaminated by the same concentrations of the two fungicides, and we assessed hatching success and hatching time. In the first experiment, cocoon production was halved at the highest copper concentration (232.5 mg Cu kg-1 of dry soil) as compared to the control. Cocoons took 5 more days to hatch, and the hatching success decreased by 35% as compared to the control. In the Swing Gold® treatments, cocoon production was reduced by 63% at 3 times the RD, and the hatching success significantly decreased by 16% at the RD. In the second experiment, only the hatching success of cocoons was impacted by Swing Gold® at 3 times the RD (30% less hatching). It is concluded that the cocoon stock in the soil is crucial for the renewal of populations in the field. The most sensitive endpoint was the hatching success of the cocoons produced by exposed adults. This endpoint and the effects observed on the "naïve" cocoons could be taken into account in pesticide risk assessment.

Aporrectodea caliginosa, a relevant earthworm species for a posteriori pesticide risk assessment: current knowledge and recommendations for culture and experimental design.

Ecotoxicological tests with earthworms are widely used and are mandatory for the risk assessment of pesticides prior to registration and commercial use. The current model species for standardized tests is Eisenia fetida or Eisenia andrei. However, these species are absent from agricultural soils and often less sensitive to pesticides than other earthworm species found in mineral soils. To move towards a better assessment of pesticide effects on non-target organisms, there is a need to perform a posteriori tests using relevant species. The endogeic species Aporrectodea caliginosa (Savigny, 1826) is representative of cultivated fields in temperate regions and is suggested as a relevant model test species. After providing information on its taxonomy, biology, and ecology, we reviewed current knowledge concerning its sensitivity towards pesticides. Moreover, we highlighted research gaps and promising perspectives. Finally, advice and recommendations are given for the establishment of laboratory cultures and experiments using this soil-dwelling earthworm species.

How to assess the feeding activity in ecotoxicological laboratory tests using enchytraeids?

The risk assessment of pesticides on soil fauna is an issue to protect agroecosystem sustainability. Enchytraeids are recognized as relevant soil bioindicators of chemical stress in agroecosystems. In laboratory, the reproduction test was found to be sensitive to reveal chemical impacts on enchytraeids. However, it does not allow to assess the impacts on ecological functions in which enchytraeids are involved. The objectives of this study were (i) to explore the feasibility of the bait-lamina test with enchytraeids under laboratory conditions and (ii) to compare its sensitivity with the Enchytraeid Reproduction Test. We exposed individuals of Enchytraeus albidus to two commercial formulations of fungicides widely used in Europe. The EC50 reproduction for the Swing® Gold (50 g L-1 epoxiconazole and 133 g L-1 dimoxystrobin) and the Cuprafor micro® (50% copper oxychloride) were respectively estimated at 1.66 ± 0.93 times the recommended dose and > 496 mg kg-1 of copper. However, no impact was found on the feeding activity of enchytraeids. The bait-lamina test thus appeared less sensitive than the Enchytraeid Reproduction Test to the tested fungicides. Despite that, this test which is achievable under laboratory conditions and allows assessing indirect effects of pesticides is quick, cheap, and easy to perform. It would deserve to be used to explore longer-exposure effects through the repeated addition of bait-lamina sticks.

Short-term effects of two fungicides on enchytraeid and earthworm communities under field conditions.

Because of the wide use of pesticides in agriculture, there is still a need of higher-tier field studies to assess ecotoxicological effects of commercial formulations on a wider variety of non-target soil organisms such as soil annelids. We here tested the effects of different concentrations of two fungicide formulations, i.e., Cuprafor Micro® (composed of 500 g kg-1 copper oxychloride) and Swing Gold® (composed of 50 g l-1 epoxiconazole and 133 g l-1 dimoxystrobin) on two families of terrestrial oligochaetes (Lumbricidae and Enchytraeidae) after 1 month of exposure. We also assessed the feeding activity of soil organisms using the bait lamina method. Along with the feeding activity, the enchytraeid density, diversity and communities were not different in the control and the contaminated plots. By contrast, epigeic earthworms were absent and earthworm diversity and densities of anecic species decreased significantly in the plots contaminated at ten times the recommended dose of the Swing Gold® formulation. The copper fungicide (at 0.75 and 7.5 kg Cu ha-1) and the treatment with the pesticide mixture (Cuprafor Micro® at 0.75 kg Cu ha-1 and Swing Gold® at the recommended dose) did not affect Oligochaeta communities compared with the control, except the Shannon index for earthworms in the mixture of both fungicides. Responses of the two annelid families to the tested pesticides were different with higher effects observed on the diversity and the community structure of earthworms compared with enchytraeids. This study allowed detecting early changes on oligochaete populations after pesticide application.

Refining uptake and depuration constants for fluoroalkyl chemicals in Chironomus riparius larvae on the basis of experimental results and modelling.

The aims of this study were to determine depuration rates for a range of per- and polyfluoroalkyl substances (PFASs) using Chironomus riparius, and to test a concentration-dependency hypothesis for the long-chain perfluorotridecanoic acid (PFTrDA) for this species. Midge larvae were exposed to field sediments collected downstream of a fluorotelomer plant, and to the same sediment spiked with PFTrDA. Elimination kinetics results indicated complete elimination of all PFASs by chironomids after 42h. These data were used to develop two PFTrDA bioaccumulation models accounting for chironomid growth and for compound concentration dependency or not. There was much better agreement between observed and simulated data under the concentration-dependency hypothesis than under the alternative one (passive diffusion). The PFTrDA uptake rate derived from the concentration-dependency model equaled 0.013 ± 0.008gocgwwh-1, and the depuration rate 0.032 ± 0.009h-1.

Differences in sensitivity between earthworms and enchytraeids exposed to two commercial fungicides.

The use of pesticides in crop fields may have negative effects on soil Oligochaeta Annelida, i.e., earthworms and enchytraeids, and thus affect soil quality. The aim of this study was to assess the effects of two commercial fungicide formulations on the earthworm Aporrectodea caliginosa and the enchytraeid Enchytraeus albidus in a natural soil. The fungicides were Cuprafor micro® (copper oxychloride), commonly used in organic farming, and Swing Gold® (epoxiconazole and dimoxystrobin), a synthetic fungicide widely used in conventional farming to protect cereal crops. Laboratory experiments were used to assess the survival, biomass loss and avoidance behaviour. No lethal effect was observed following exposure to the copper fungicide for 14 days, even at 5000mgkg-1 of copper, i.e. 650 times the recommended dose (RD). However, a significant decrease in biomass was observed from 50mgkg-1 of copper (6.5 times the RD) for A. caliginosa and at 5000mgkg-1 of copper (650 times the RD) for E. albidus. These sublethal effects suggest that a longer period of exposure would probably have led to lethal effects. The EC50 avoidance for the copper fungicide was estimated to be 51.2mgkg-1 of copper (6.7 times the RD) for A. caliginosa, and 393mgkg-1 of copper (51 times the RD) for E. albidus. For the Swing Gold® fungicide, the estimated LC50 was 7.0 10-3mLkg-1 (6.3 times the RD) for A. caliginosa and 12.7 10-3mLkg-1 (11.0 times the RD) for E. albidus. No effect on biomass or avoidance was observed at sublethal concentrations of this synthetic fungicide. It was concluded that enchytraeids were less sensitive than earthworms to the two commercial fungicides in terms of mortality, biomass loss and avoidance behaviour. Therefore we discuss the different strategies possibly used by the two Oligochaeta species to cope with the presence of the pesticides were discussed, along with the potential consequences on the soil functions.

Multiscale modelling approaches for assessing cosmetic ingredients safety.

The European Union's ban on animal testing for cosmetic ingredients and products has generated a strong momentum for the development of in silico and in vitro alternative methods. One of the focus of the COSMOS project was ab initio prediction of kinetics and toxic effects through multiscale pharmacokinetic modeling and in vitro data integration. In our experience, mathematical or computer modeling and in vitro experiments are complementary. We present here a summary of the main models and results obtained within the framework of the project on these topics. A first section presents our work at the organelle and cellular level. We then go toward modeling cell levels effects (monitored continuously), multiscale physiologically based pharmacokinetic and effect models, and route to route extrapolation. We follow with a short presentation of the automated KNIME workflows developed for dissemination and easy use of the models. We end with a discussion of two challenges to the field: our limited ability to deal with massive data and complex computations.

Zebrafish-based reporter gene assays reveal different estrogenic activities in river waters compared to a conventional human-derived assay.

Endocrine disrupting chemicals (EDCs) act on the endocrine system through multiple mechanisms of action, among them interaction with estrogen receptors (ERs) is a well-identified key event in the initiation of adverse outcomes. As the most commonly used estrogen screening assays are either yeast- or human-cell based systems, the question of their (eco)toxicological relevance when assessing risks for aquatic species can be raised. The present study addresses the use of zebrafish (zf) derived reporter gene assays, both in vitro (i.e. zf liver cell lines stably expressing zfERα, zfERß1 and zfERß2 subtypes) and in vivo (i.e. transgenic cyp19a1b-GFP zf embryos), to assess estrogenic contaminants in river waters. By investigating 20 French river sites using passive sampling, high frequencies of in vitro zfER-mediated activities in water extracts were measured. Among the different in vitro assays, zfERß2 assay was the most sensitive and responsive one, enabling the detection of active compounds at all investigated sites. In addition, comparison with a conventional human-based in vitro assay highlighted sites that were able to active zfERs but not human ER, suggesting the occurrence of zf-specific ER ligands. Furthermore, a significant in vivo estrogenic activity was detected at the most active sites in vitro, with a good accordance between estradiol equivalent (E2-EQ) concentrations derived from both in vitro and in vivo assays. Overall, this study shows the relevance and usefulness of such novel zebrafish-based assays as screening tools to monitor estrogenic activities in complex mixtures such as water extracts. It also supports their preferred use compared to human-based assays to assess the potential risks caused by endocrine disruptive chemicals for aquatic species such as fish.

Prediction of the metabolic clearance of benzophenone-2, and its interaction with isoeugenol and coumarin using cryopreserved human hepatocytes in primary culture.

Benzophenone-2 (BP2) is widely used as a UV screen in both industrial products and cosmetic formulations, where it is frequently found associated with fragrance compounds, such as isoeugenol and coumarin. BP2 is now recognized as an endocrine disruptor, but to date, no information has been reported on its fate in humans. The intrinsic clearance (Clint) and metabolic interactions of BP2 were explored using cryopreserved human hepatocytes in primary cultures. In vitro kinetic experiments were performed to estimate the Michaelis-Menten parameters. The substrate depletion method demonstrated that isoeugenol was cleared more rapidly than BP2 or coumarin (Clint = 259, 94.7 and 0.40 µl/min/10(6) cells respectively). This vitro model was also used to study the metabolic interactions between BP2 and isoeugenol and coumarin. Coumarin exerted no effects on either isoeugenol or BP2 metabolism, because of its independent metabolic pathway (CYP2A6). Isoeugenol appeared to be a potent competitive substrate inhibitor of BP2 metabolism, equivalent to the specific UGT1A1 substrate: estradiol. Despite the fact that inhibition of UGT by xenobiotics is not usually considered to be a major concern, the involvement of UGT1A1 in BP2 metabolism may have pharmacokinetic and pharmacological consequences, due to the its polymorphisms in humans and its pure estrogenic effect.

Transgenerational Adaptation to Pollution Changes Energy Allocation in Populations of Nematodes.

Assessing the evolutionary responses of long-term exposed populations requires multigeneration ecotoxicity tests. However, the analysis of the data from these tests is not straightforward. Mechanistic models allow the in-depth analysis of the variation of physiological traits over many generations, by quantifying the trend of the physiological and toxicological parameters of the model. In the present study, a bioenergetic mechanistic model has been used to assess the evolution of two populations of the nematode Caenorhabditis elegans in control conditions or exposed to uranium. This evolutionary pressure resulted in a brood size reduction of 60%. We showed an adaptation of individuals of both populations to experimental conditions (increase of maximal length, decrease of growth rate, decrease of brood size, and decrease of the elimination rate). In addition, differential evolution was also highlighted between the two populations once the maternal effects had been diminished after several generations. Thus, individuals that were greater in maximal length, but with apparently a greater sensitivity to uranium were selected in the uranium population. In this study, we showed that this bioenergetics mechanistic modeling approach provided a precise, certain, and powerful analysis of the life strategy of C. elegans populations exposed to heavy metals resulting in an evolutionary pressure across successive generations.

An individual-based model of zebrafish population dynamics accounting for energy dynamics.

Developing population dynamics models for zebrafish is crucial in order to extrapolate from toxicity data measured at the organism level to biological levels relevant to support and enhance ecological risk assessment. To achieve this, a dynamic energy budget for individual zebrafish (DEB model) was coupled to an individual based model of zebrafish population dynamics (IBM model). Next, we fitted the DEB model to new experimental data on zebrafish growth and reproduction thus improving existing models. We further analysed the DEB-model and DEB-IBM using a sensitivity analysis. Finally, the predictions of the DEB-IBM were compared to existing observations on natural zebrafish populations and the predicted population dynamics are realistic. While our zebrafish DEB-IBM model can still be improved by acquiring new experimental data on the most uncertain processes (e.g. survival or feeding), it can already serve to predict the impact of compounds at the population level.

Biodistribution and Clearance of TiO2 Nanoparticles in Rats after Intravenous Injection.

Titanium dioxide (TiO2) nanoparticles are used in many applications. Due to their small size, easy body penetration and toxicological adverse effects have been suspected. Numerous studies have tried to characterize TiO2 translocation after oral, dermal or respiratory exposure. In this study, we focused on TiO2 nanoparticle biodistribution, clearance and toxicological effects after intravenous injection, considering TiO2 translocation in the blood occurs. Using ICP-OES, transmission electron microscopy, and histological methods, we found TiO2 accumulation in liver, lungs and spleen. We estimated TiO2 nanoparticles' half life in the body to about 10 days. Clinical biomarkers were also quantified for 56 days to identify potential toxicological impact on lungs, blood, liver, spleen and kidneys. Results showed absence of toxicological effects after TiO2 intravenous injection at concentrations of 7.7 to 9.4 mg/kg.

Microscale Ecotoxicity Testing of Moselle River Watershed (Lorraine Province, France) Sediments.

The ecotoxic potential of seven Moselle river watershed sediments was assessed with a battery of bioassays comprised of rapid phototrophic [LuminoTox solid phase (L-SPA) and elutriate (L-ELU) assays] and bacterial [Microtox solid phase assay (M-SPA)] exposure tests, as well as with two micro-invertebrate solid phase tests conducted with Hydra attenuata (lethal and sublethal effects solid phase assay, HL-SPA and HSL-SPA) and Chironomus riparius. Measured effects of sediments and their elutriates were varied and reflected responses that were ecotoxicity test-, endpoint- and site-dependent, suggesting some degree of risk toward benthic and water column organisms, respectively, at specific sites. Correlation analysis demonstrated that L-SPA and M-SPA ecotoxicity responses were significantly linked with the Hydra HSL-SPA assay, indicating their ability to predict ecotoxicity towards an invertebrate taxonomic group representing secondary consumers. While the L-SPA and M-SPA assays hold promise as rapid screens for sediment ecotoxicity, correlation analysis with grain size (L-SPA: r=-0.795, P=0.033; M-SPA: r=-0.73, P=0.07) points out that their responses can be influenced by the presence of fines (i.e., sediment particles ≤0.063 mm in size) and that this information is essential to properly interpret ecotoxicity data generated with these assays. Finally, notable differences observed in trophic level sensitivities once again recall the importance of employing a test battery to adequately appraise the ecotoxicity of sediments.

Energy-based modelling to assess effects of chemicals on Caenorhabditis elegans: a case study on uranium.

The ubiquitous free-living nematode Caenorhabditis elegans is a powerful animal model for measuring the evolutionary effects of pollutants which is increasingly used in (eco) toxicological studies. Indeed, toxicity tests with this nematode can provide in a few days data on the whole life cycle. These data can be analysed with mathematical tools such as toxicokinetic-toxicodynamic modelling approaches. In this study, we assessed how a chronic exposure to a radioactive heavy metal (uranium) affects the life-cycle of C. elegans using a mechanistic model. In order to achieve this, we exposed individuals to a range of seven concentrations of uranium. Growth and reproduction were followed daily. These data were analysed with a model for nematodes based on the Dynamic Energy Budget theory, able to handle a wide range of plausible biological parameters values. Parameter estimations were performed using a Bayesian framework. Our results showed that uranium affects the assimilation of energy from food with a no-effect concentration (NEC) of 0.42 mM U which would be the threshold for effects on both growth and reproduction. The sensitivity analysis showed that the main contributors to the model output were parameters linked to the feeding processes and the actual exposure concentration. This confirms that the real exposure concentration should be measured accurately and that the feeding parameters should not be fixed, but need to be reestimated during the parameter estimation process.

A physiologically based toxicokinetic model for the zebrafish Danio rerio.

Zebrafish (Danio rerio) is a widely used model for toxicological studies, in particular those related to investigations on endocrine disruption. The development and regulatory use of in vivo and in vitro tests based on this species can be enhanced by toxicokinetic modeling. For this reason, we propose a physiologically based toxicokinetic (PBTK) model for zebrafish describing the uptake and disposition of organic chemicals. The model is based on literature data on zebrafish, other cyprinidae and other fish families, new experimental physiological information (volumes, lipids and water contents) obtained from zebrafish, and chemical-specific parameters predicted by generic models. The relevance of available models predicting the latter parameters was evaluated with respect to gill uptake and partition coefficients in zebrafish. This evaluation benefited from the fact that the influence of confounding factors such as body weight and temperature on ventilation rate was included in our model. The predictions for six chemicals (65 data points) yielded by our PBTK model were compared to available toxicokinetics data for zebrafish and 88% of them were within a factor of 5 of the corresponding experimental values. Sensitivity analysis highlighted that the 1-octanol/water partition coefficient, the metabolism rate, and all the parameters that enable the prediction of assimilation efficiency and partitioning of chemicals need to be precisely determined in order to allow an effective toxicokinetic modeling.