A workflow to build PBTK models for novel species.

Physiology-based pharmacokinetic and toxicokinetic (PBPK/TK) models allow us to simulate the concentration of xenobiotica in the plasma and different tissues of an organism. PBPK/TK models are therefore routinely used in many fields of life sciences to simulate the physiological concentration of exogenous compounds in plasma and tissues. The application of PBTK models in ecotoxicology, however, is currently hampered by the limited availability of models for focal species. Here, we present a best practice workflow that describes how to build PBTK models for novel species. To this end, we extrapolated eight previously established rabbit models for several drugs to six additional mammalian species (human, beagle, rat, monkey, mouse, and minipig). We used established PBTK models for these species to account for the species-specific physiology. The parameter sensitivity in the resulting 56 PBTK models was systematically assessed to rank the relevance of the parameters on overall model performance. Interestingly, more than 80% of the 609 considered model parameters showed a negligible sensitivity throughout all models. Only approximately 5% of all parameters had a high sensitivity in at least one of the PBTK models. This approach allowed us to rank the relevance of the various parameters on overall model performance. We used this information to formulate a best practice guideline for the efficient development of PBTK models for novel animal species. We believe that the workflow proposed in this study will significantly support the development of PBTK models for new animal species in the future.

Demographic Toxicokinetic-Toxicodynamic Modeling of Lethal Effects.

The aquatic effect assessment of chemicals is largely based on standardized measures of toxicity determined in short-term laboratory tests which are designed to reduce variability. For this purpose, uniform individuals of a species are kept under environmental and chemical exposure conditions which are as constant as possible. In nature, exposure often appears to be pulsed, effects might last longer than a few days, sensitivity might vary among different sized organisms and populations are usually size or age structured and are subject to demographic processes. To overcome this discrepancy, we tested toxicokinetic-toxicodynamic models of different complexities, including body size scaling approaches, for their ability to represent lethal effects observed for Daphnia magna exposed to triphenyltin. The consequences of the different toxicokinetic and toxicodynamic assumptions for population level responses to pulsed exposure are tested by means of an individual based model and are evaluated by confronting model predictions with population data for various pulsed exposure scenarios. We provide an example where increased model complexity reduces the uncertainty in model outputs. Furthermore, our results emphasize the importance of considering population demography in toxicokinetics and toxicodynamics for understanding and predicting potential chemical impacts at higher levels of biological organization.

Cross-Species Extrapolation of Uptake and Disposition of Neutral Organic Chemicals in Fish Using a Multispecies Physiologically-Based Toxicokinetic Model Framework.

The potential to bioconcentrate is generally considered to be an unwanted property of a substance. Consequently, chemical legislation, including the European REACH regulations, requires the chemical industry to provide bioconcentration data for chemicals that are produced or imported at volumes exceeding 100 tons per annum or if there is a concern that a substance is persistent, bioaccumulative, and toxic. For the filling of the existing data gap for chemicals produced or imported at levels that are below this stipulated volume, without the need for additional animal experiments, physiologically-based toxicokinetic (PBTK) models can be used to predict whole-body and tissue concentrations of neutral organic chemicals in fish. PBTK models have been developed for many different fish species with promising results. In this study, we developed PBTK models for zebrafish (Danio rerio) and roach (Rutilus rutilus) and combined them with existing models for rainbow trout (Onchorhynchus mykiss), lake trout (Salvelinus namaycush), and fathead minnow (Pimephales promelas). The resulting multispecies model framework allows for cross-species extrapolation of the bioaccumulative potential of neutral organic compounds. Predictions were compared with experimental data and were accurate for most substances. Our model can be used for probabilistic risk assessment of chemical bioaccumulation, with particular emphasis on cross-species evaluations.

Ecological Recovery Potential of Freshwater Organisms: Consequences for Environmental Risk Assessment of Chemicals.

Chemical contaminants released into the in the environment may have adverse effects on (non-target) species, populations and communities. The return of a stressed system to its pre-disturbance or other reference state, i.e. the ecological recovery, may depend on various factors related to the affected taxon, the ecosystem of concern and the type of stressor with consequences for the assessment and management of risks associated with chemical contaminants. Whereas the effects caused by short-term exposure might be acceptable to some extent, the conditions under which ecological recovery can serve as a decision criterion in the environmental risk assessment of chemical stressors remains to be evaluated. For a generic consideration of recovery in the risk assessment of chemicals, we reviewed case studies of natural and artificial aquatic systems and evaluate five aspects that might cause variability in population recovery time: (1) taxonomic differences and life-history variability, (2) factors related to ecosystem type and community processes, (3) type of disturbance, (4) comparison of field and semi-field studies, and (5) effect magnitude, i.e., the decline in population size following disturbance. We discuss our findings with regard to both retrospective assessments and prospective risk assessment.

Family-portraits for daphnids: scanning living individuals and populations to measure body length.

A method has been developed and tested to determine the body length of living daphnids. The purpose of the method was the simple, accurate, repeatable, quick, and to the living organism, harmless measurement of body length of all individuals in a population to enhance the capability of observing population development over time. Generally, organisms are transferred to a petri dish and temporarily fixed by removing access medium. A picture of the petri dish is taken using an ordinary flatbed scanner. Pictures are manually analysed with purposely developed software. We found no significant impact of the method on either individual performance (growth and reproduction) or population development (abundance and structure) of daphnids in comparison to the previously used method for data gathering (sieving, counting and length measurement of a subsample via microscopy). The disadvantage of our method, an increased demand in time for picture analysis, is negligible compared to the advantages this method has. Data generated with the new method do represent the population structure more accurately than those data generated with the previously used method. Scanning organisms does also allow a retrospective quality control for generated data as pictures can securely be stored. The quality of the pictures is furthermore sufficient to include additional endpoints to the analysis (e.g., number and size of aborts, number and size of eggs in the brood pouch, spine length). Here, we present, test and discuss an alternative approach to automated image analysis for data gathering in single and multiple individual and species experiments.

Population level effects of multiwalled carbon nanotubes in Daphnia magna exposed to pulses of triclocarban.

Due to the rapid increase of carbon nanotubes (CNT) applications and their inevitable release into the aquatic environment, CNT may interact with and further influence the fate and transport of other pollutants such as triclocarban (TCC). TCC is a high-production-volume chemical that is widely used as an antimicrobial agent, is continually released into the aquatic environment, and is biologically active and persistent. In the present study, the population test with Daphnia magna was performed over 93 days. Different treatments were examined: (a) control, (b) solvent control, (c) 1 mg CNT/L from the beginning, (d) 1 mg CNT/L as of day 14, (e) control with a 2-day pulse of 25 µg TCC/L on day 14, 41 µg TCC/L (day 54), and 61 µg TCC/L (day 68) and (f) same pulses of TCC with co-exposure to 1 mg CNT/L. Significant changes in all three size classes were observed as a result of the long-term exposure to 1 mg CNT/L. Increasing in number of neonates, and decreasing in number of juveniles and adults were observed. Moreover, daphnids were significantly smaller when they were exposed to MWCNT. The exposure with TCC led to size-dependent mortality in Daphnia magna populations and a subsequent recovery. Lower toxicity of TCC was observed, with the presence of MWCNT in the medium. The reported effects of TCC on population level were compared to the output of an individual-based Daphnia magna population model, in order to verify the model predictions with laboratory data.

Limitations of extrapolating toxic effects on reproduction to the population level.

For the ecological risk assessment of toxic chemicals, standardized tests on individuals are often used as proxies for population-level effects. Here, we address the utility of one commonly used metric, reproductive output, as a proxy for population-level effects. Because reproduction integrates the outcome of many interacting processes (e.g., feeding, growth, allocation of energy to reproduction), the observed toxic effects in a reproduction test could be due to stress on one of many processes. Although this makes reproduction a robust endpoint for detecting stress, it may mask important population-level consequences if the different physiological processes stress affects are associated with different feedback mechanisms at the population level. We therefore evaluated how an observed reduction in reproduction found in a standard reproduction test translates to effects at the population level if it is caused by hypothetical toxicants affecting different physiological processes (physiological modes of action; PMoA). For this we used two consumer­resource models: the Yodzis-Innes (YI) model, which is mathematically tractable, but requires strong assumptions of energetic equivalence among individuals as they progress through ontogeny, and an individual-based implementation of dynamic energy budget theory (DEB-IBM), which relaxes these assumptions at the expense of tractability. We identified two important feedback mechanisms controlling the link between individual- and population-level stress in the YI model. These mechanisms turned out to also be important for interpreting some of the individual-based model results; for two PMoAs, they determined the population response to stress in both models. In contrast, others stress types involved more complex feedbacks, because they asymmetrically stressed the production efficiency of reproduction and somatic growth. The feedbacks associated with different PMoAs drastically altered the link between individual- and population-level effects. For example, hypothetical stressors with different PMoAs that had equal effects on reproduction had effects ranging from a negligible decline in biomass to population extinction. Thus, reproduction tests alone are of little use for extrapolating toxicity to the population level, but we showed that the ecological relevance of standard tests could easily be improved if growth is measured along with reproduction.

Understanding receptor-mediated effects in rainbow trout: in vitro-in vivo extrapolation using physiologically based toxicokinetic models.

The European REACH regulation requires the use of animal experimentation to assess the risk of industrial chemicals. However, the 3R principle (reduction, replacement, refinement) demands the use of suitable alternative test methods. Many dossiers submitted for the authorization of chemicals have attempted to provide the required data without performing new experiments, relying heavily on in silico methods; in vitro assays were scarcely used. We propose a methodology that uses physiologically based toxicokinetic (PBTK) models to extrapolate in vitro data to the in vivo level. We collected experimental results for in vitro and in vivo ethoxyresorufin-O-deethylase and vitellogenin induction following chemical exposure and compared those results with model predictions. We found that the predictive power of aqueous chemical concentrations was limited; median effect concentrations (EC50s) based on internal concentrations in fish correlated better with in vitro EC50s. Our data show that in vitro assays could offer a substitute for fish studies when combined with PBTK models.

A spatiotemporally explicit modeling approach for more realistic exposure and risk assessment of off-field soil organisms.

Natural and seminatural habitats of soil living organisms in cultivated landscapes can be subject to unintended exposure by active substances of plant protection products (PPPs) used in adjacent fields. Spray-drift deposition and runoff are considered major exposure routes into such off-field areas. In this work, we develop a model (xOffFieldSoil) and associated scenarios to estimate exposure of off-field soil habitats. The modular model approach consists of components, each addressing a specific aspect of exposure processes, for example, PPP use, drift deposition, runoff generation and filtering, estimation of soil concentrations. The approach is spatiotemporally explicit and operates at scales ranging from local edge-of-field to large landscapes. The outcome can be aggregated and presented to the risk assessor in a way that addresses the dimensions and scales defined in specific protection goals (SPGs). The approach can be used to assess the effect of mitigation options, for example, field margins, in-field buffers, or drift-reducing technology. The presented provisional scenarios start with a schematic edge-of-field situation and extend to real-world landscapes of up to 5 km × 5 km. A case study was conducted for two active substances of different environmental fate characteristics. Results are presented as a collection of percentiles over time and space, as contour plots, and as maps. The results show that exposure patterns of off-field soil organisms are of a complex nature due to spatial and temporal variabilities combined with landscape structure and event-based processes. Our concepts and analysis demonstrate that more realistic exposure data can be meaningfully consolidated to serve in standard-tier risk assessments. The real-world landscape-scale scenarios indicate risk hot-spots that support the identification of efficient risk mitigation. As a next step, the spatiotemporally explicit exposure data can be directly coupled to ecological effect models (e.g., for earthworms or collembola) to conduct risk assessments at biological entity levels as required by SPGs. Integr Environ Assess Manag 2024;20:263-278. © 2023 Applied Analysis Solutions LLC and WSC Scientific GmbH and Bayer AG and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Comparing Sensitivity of Different Bee Species to Pesticides: A TKTD modeling approach.

Risk assessment for bees is mainly based on data for honey bees; however, risk assessment is intended to protect all bee species. This raises the question of whether data for honey bees are a good proxy for other bee species. This issue is not new and has resulted in several publications in which the sensitivity of bee species is compared based on the values of the 48-h median lethal dose (LD50) from acute test results. When this approach is used, observed differences in sensitivity may result both from differences in kinetics and from inherent differences in species sensitivity. In addition, the physiology of the bee, like its overall size, the size of the honey stomach (for acute oral tests), and the physical appearance (for acute contact tests) also influences the sensitivity of the bee. The recently introduced Toxicokinetic-Toxicodynamic (TKTD) model that was developed for the interpretation of honey bee tests (Bee General Uniform Threshold Model for Survival [BeeGUTS]) could integrate the results of acute oral tests, acute contact tests, and chronic tests within one consistent framework. We show that the BeeGUTS model can be calibrated and validated for other bee species and also that the honey bee is among the more sensitive bee species. In addition, we found that differences in sensitivity between species are smaller than previously published comparisons based on 48-h LD50 values. The time-dependency of the LD50 and the specifics of the bee physiology are the main causes of the wider variation found in the published literature. Environ Toxicol Chem 2024;43:1431-1441. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Predicting population dynamics from the properties of individuals: a cross-level test of dynamic energy budget theory.

Individual-based models (IBMs) are increasingly used to link the dynamics of individuals to higher levels of biological organization. Still, many IBMs are data hungry, species specific, and time-consuming to develop and analyze. Many of these issues would be resolved by using general theories of individual dynamics as the basis for IBMs. While such theories have frequently been examined at the individual level, few cross-level tests exist that also try to predict population dynamics. Here we performed a cross-level test of dynamic energy budget (DEB) theory by parameterizing an individual-based model using individual-level data of the water flea, Daphnia magna, and comparing the emerging population dynamics to independent data from population experiments. We found that DEB theory successfully predicted population growth rates and peak densities but failed to capture the decline phase. Further assumptions on food-dependent mortality of juveniles were needed to capture the population dynamics after the initial population peak. The resulting model then predicted, without further calibration, characteristic switches between small- and large-amplitude cycles, which have been observed for Daphnia. We conclude that cross-level tests help detect gaps in current individual-level theories and ultimately will lead to theory development and the establishment of a generic basis for individual-based models and ecology.

Feeding inhibition explains effects of imidacloprid on the growth, maturation, reproduction, and survival of Daphnia magna.

Effects of some xenobiotics on aquatic organisms might not be caused directly by the compound but rather arise from acclimation of the organism to stress invoked by feeding inhibition during exposure. Experiments were conducted to identify effects of imidacloprid on individual performance (feeding, growth, maturation, reproduction, and survival) of Daphnia magna under surplus and reduced food availability. Concentrations inhibiting feeding by 5, 50, and 95% after one day of exposure were 0.19, 1.83, and 8.70 mg/L, respectively. Exposure with imidacloprid at ≥ 3.7 mg/L reduced growth by up to 53 ± 11% within one week. Surplus food availability after inhibition allowed recovery from this growth inhibition, whereas limited food supply eliminated the potential for recovery in growth even for exposure at 0.15 mg/L. A shift in the distribution of individual energy reserves toward reproduction rather than growth resulted in increased reproduction after exposure to concentrations ≤ 0.4 mg/L. Exposure to imidacloprid at ≥ 4.0 mg/L overwhelmed this adaptive response and reduced reproduction by up to 57%. We used the individual based Daphnia magna population model IDamP as a virtual laboratory to demonstrate that only feeding was affected by imidacloprid, and that in turn this caused the other impacts on individual performance. Consideration of end points individually would have led to a different interpretation of the effects. Thus, we demonstrate how multiple lines of evidence linked by understanding the ecology of the organism are necessary to elucidate xenobiotic impacts along the effect cascade.

Extrapolating ecotoxicological effects from individuals to populations: a generic approach based on Dynamic Energy Budget theory and individual-based modeling.

Individual-based models (IBMs) predict how dynamics at higher levels of biological organization emerge from individual-level processes. This makes them a particularly useful tool for ecotoxicology, where the effects of toxicants are measured at the individual level but protection goals are often aimed at the population level or higher. However, one drawback of IBMs is that they require significant effort and data to design for each species. A solution would be to develop IBMs for chemical risk assessment that are based on generic individual-level models and theory. Here we show how one generic theory, Dynamic Energy Budget (DEB) theory, can be used to extrapolate the effect of toxicants measured at the individual level to effects on population dynamics. DEB is based on first principles in bioenergetics and uses a common model structure to model all species. Parameterization for a certain species is done at the individual level and allows to predict population-level effects of toxicants for a wide range of environmental conditions and toxicant concentrations. We present the general approach, which in principle can be used for all animal species, and give an example using Daphnia magna exposed to 3,4-dichloroaniline. We conclude that our generic approach holds great potential for standardized ecological risk assessment based on ecological models. Currently, available data from standard tests can directly be used for parameterization under certain circumstances, but with limited extra effort standard tests at the individual would deliver data that could considerably improve the applicability and precision of extrapolation to the population level. Specifically, the measurement of a toxicant's effect on growth in addition to reproduction, and presenting data over time as opposed to reporting a single EC50 or dose response curve at one time point.

A Framework for Algae Modeling in Regulatory Risk Assessment.

The use of toxicokinetic-toxicodynamic (TKTD) modeling in regulatory risk assessment of plant protection products is increasingly popular, especially since the 2018 European Food Safety Authority (EFSA) opinion on TKTD modeling announced that several established models are ready for use in risk assessment. With careful adherence to the guidelines laid out by EFSA, we present a stepwise approach to validation and use of the Simple Algae Model Extended (SAM-X) for regulatory submission in Tier 2C. We demonstrate how the use of moving time windows across time-variable exposure profiles can generate thousands of virtual laboratory mimic simulations that seamlessly predict the effects of time-variable exposures across a full exposure profile while maintaining the laboratory conditions of the standard Organisation for Economic Co-operation and Development (OECD) growth inhibition test. Thus, every virtual laboratory test has a duration of 72 h, with OECD medium and constant light and temperature conditions. The only deviation from the standard test setup is the replacement of constant exposure conditions for time-variable concentrations. The present study demonstrates that for simulation of 72-h toxicity tests, the nutrient dynamics in the SAM-X model are not required, and we propose the alternative use of a simplified model version. For risk assessment, in accordance with the EFSA guidelines we use a median exposure profile of 10 as a threshold, meaning that if a time window within the exposure profile causes 50% growth inhibition when magnified by a factor of 10, the threshold will have been exceeded. We present a simplified example for chlorotoluron and isoproturon. The present case study brings to life our proposed framework for TKTD modeling of algae to establish whether a given exposure can be considered to be of low risk. Environ Toxicol Chem 2023;42:1823-1838. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Aquatic macrophyte growth season in Central and Northern European Union and the implications for aquatic macrophyte risk assessments for herbicides.

Under current European Union regulation, the risks to aquatic organisms must be assessed for uses of plant protection products (PPPs) that may result in exposure to the environment. For herbicidal PPPs, aquatic macrophytes are often the most sensitive taxa. For some herbicidal modes of action, macrophytes may be affected only while they are actively growing. For the risk assessment, it is therefore useful to know whether application timings would result in surface water exposure during periods when aquatic macrophytes are actively growing (therefore potentially resulting in effects). Toxicity endpoints, which are based on studies with active growth, may be overconservative in cases where exposure of PPPs will not co-occur with active macrophyte growth. A comprehensive literature search was performed, using systematic and manual approaches, with the aim of identifying the main active growth period for macrophytes in natural freshwater bodies in climates relevant to the Central and Northern zones of the European Union. The results of the searches were screened initially to identify all potentially relevant references, for which a full evaluation was then performed. Reliability was assessed using the principles of the Klimisch scoring system. As part of the full evaluation, growth periods were identified for each macrophyte species studied. Finally, the extracted growth periods were considered together to determine an overall active growth period for aquatic macrophytes representative of the Central and Northern EU zones. Based on this literature review, the active growth period identified for most aquatic macrophyte species representative of the Central and Northern EU zones is April to September. Relating to the regulatory implication of these results, it may be possible to conclude a low risk for aquatic macrophytes if the predicted surface water exposure period for certain PPPs is demonstrated to be outside the periods of active growth. Integr Environ Assess Manag 2023;00:1-15. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Evaluating and Explaining the Variability of Honey Bee Field Studies across Europe Using BEEHAVE.

To assess the effect of plant protection products on pollinator colonies, the higher tier of environmental risk assessment (ERA), for managed honey bee colonies and other pollinators, is in need of a mechanistic effect model. Such models are seen as a promising solution to the shortcomings, which empirical risk assessment can only overcome to a certain degree. A recent assessment of 40 models conducted by the European Food Safety Authority (EFSA) revealed that BEEHAVE is currently the only publicly available mechanistic honey bee model that has the potential to be accepted for ERA purposes. A concern in the use of this model is a lack of model validation against empirical data, spanning field studies conducted in different regions of Europe and covering the variability in colony and environmental conditions. We filled this gap with a BEEHAVE validation study against 66 control colonies of field studies conducted across Germany, Hungary, and the United Kingdom. Our study implements realistic initial colony size and landscape structure to consider foraging options. Overall, the temporal pattern of colony strength is predicted well. Some discrepancies between experimental data and prediction outcomes are explained by assumptions made for model parameterization. Complementary to the recent EFSA study using BEEHAVE, our validation covers a large variability in colony conditions and environmental impacts representing the Northern and Central European Regulatory Zones. Thus we believe that BEEHAVE can be used to serve the development of specific protection goals as well as the development of simulation scenarios for the European Regulatory Zone. Subsequently, the model can be applied as a standard tool for higher tier ERA of managed honey bees using the mechanistic ecotoxicological module for BEEHAVE, BEEHAVEecotox . Environ Toxicol Chem 2023;42:1839-1850. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

BeeGUTS-A Toxicokinetic-Toxicodynamic Model for the Interpretation and Integration of Acute and Chronic Honey Bee Tests.

Understanding the survival of honey bees after pesticide exposure is key for environmental risk assessment. Currently, effects on adult honey bees are assessed by Organisation for Economic Co-operation and Development standardized guidelines, such as the acute and chronic oral exposure and acute contact exposure tests. The three different tests are interpreted individually, without consideration that the same compound is investigated in the same species, which should allow for an integrative assessment. In the present study we developed, calibrated, and validated a toxicokinetic-toxicodynamic model with 17 existing data sets on acute and chronic effects for honey bees. The model is based on the generalized unified threshold model for survival (GUTS), which is able to integrate the different exposure regimes, taking into account the physiology of the honey bee: the BeeGUTS model. The model is able to accurately describe the effects over time for all three exposure routes combined within one consistent framework. The model can also be used as a validity check for toxicity values used in honey bee risk assessment and to conduct effect assessments for real-life exposure scenarios. This new integrative approach, moving from single-point estimates of toxicity and exposure to a holistic link between exposure and effect, will allow for a higher confidence of honey bee toxicity assessment in the future. Environ Toxicol Chem 2022;41:2193-2201. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

TWAc-Check: A New Approach to Determine the Appropriate Use of Time-Weighted Average Concentration in Aquatic Risk Assessment.

In pesticide risk assessment, regulatory acceptable concentrations for surface water bodies (RACsw,ch) are used that are derived from standard studies with continuous exposure of organisms to a test compound for days or months. These RACsw,ch are compared with the maximum tested concentration of more realistic exposure scenarios. However, the actual exposure duration could be notably shorter (e.g., hours) than the standard study, which intentionally leads to an overly conservative Tier 1 risk assessment. This discrepancy can be addressed in a risk assessment using the time-weighted average concentration (TWAc). In Europe, the applicability of TWAc for a particular risk assessment is evaluated using a complex decision scheme, which has been controversial; thus we propose an alternative approach: We used TWAc-check (which is based on the idea that the TWAc concept is just a model for aquatic risk assessment) to test whether the use of a TWAc is appropriate for such assessment. The TWAc-check method works by using predicted-measured diagrams to test how well the TWAc model predicts experimental data from peak exposure experiments. Overestimated effects are accepted because the conservatism of the TWAc model is prioritized over the goodness of fit. We illustrate the applicability of TWAc-check by applying it to various data sets for different species and substances. We demonstrate that the applicability is case dependent. Specifically, TWAc-check correctly identifies that the use of TWAc is not appropriate for early onset of effects or delayed effects. The proposed concept shows that the time window is a decisive factor as to whether or not the model is acceptable and that this concept can be used as a potential refinement option prior to the use of toxicokinetic-toxicodynamic models. Environ Toxicol Chem 2022;41:1778-1787. © 2022 Bayer AG. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The BEEHAVEecotox Model-Integrating a Mechanistic Effect Module into the Honeybee Colony Model.

Mechanistic effect models are powerful tools for extrapolating from laboratory studies to field conditions. For bees, several good models are available that can simulate colony dynamics. Controlled and reliable experimental systems are also available to estimate the inherent toxicity of pesticides to individuals. However, there is currently no systematic and mechanistic way of linking the output of experimental ecotoxicological testing to bee models for bee risk assessment. We introduce an ecotoxicological module that mechanistically links exposure with the hazard profile of a pesticide for individual honeybees so that colony effects emerge. This mechanistic link allows the translation of results from standard laboratory studies to relevant parameters and processes for simulating bee colony dynamics. The module was integrated into the state-of-the-art honeybee model BEEHAVE. For the integration, BEEHAVE was adapted to mechanistically link the exposure and effects on different cohorts to colony dynamics. The BEEHAVEecotox model was tested against semifield (tunnel) studies, which were deemed the best study type to test whether BEEHAVEecotox predicted realistic effect sizes under controlled conditions. Two pesticides used as toxic standards were chosen for this validation to represent two different modes of action: acute mortality of foragers and chronic brood effects. The ecotoxicological module was able to predict effect sizes in the tunnel studies based on information from standard laboratory tests. In conclusion, the BEEHAVEecotox model is an excellent tool to be used for honeybee risk assessment, interpretation of field and semifield studies, and exploring the efficiency of different mitigation measures. The principles for exposure and effect modules are portable and could be used for any well-constructed honeybee model. Environ Toxicol Chem 2022;41:2870-2882. © 2022 Bayer AG & Sygenta, et al. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Quantitative Morphometric, Physiological, and Metabolic Characteristics of Chickens and Mallards for Physiologically Based Kinetic Model Development.

Physiologically based kinetic (PBK) models are a promising tool for xenobiotic environmental risk assessment that could reduce animal testing by predicting in vivo exposure. PBK models for birds could further our understanding of species-specific sensitivities to xenobiotics, but would require species-specific parameterization. To this end, we summarize multiple major morphometric and physiological characteristics in chickens, particularly laying hens (Gallus gallus) and mallards (Anas platyrhynchos) in a meta-analysis of published data. Where such data did not exist, data are substituted from domesticated ducks (Anas platyrhynchos) and, in their absence, from chickens. The distribution of water between intracellular, extracellular, and plasma is similar in laying hens and mallards. Similarly, the lengths of the components of the small intestine (duodenum, jejunum, and ileum) are similar in chickens and mallards. Moreover, not only are the gastrointestinal absorptive areas similar in mallard and chickens but also they are similar to those in mammals when expressed on a log basis and compared to log body weight. In contrast, the following are much lower in laying hens than mallards: cardiac output (CO), hematocrit (Hct), and blood hemoglobin. There are shifts in ovary weight (increased), oviduct weight (increased), and plasma/serum concentrations of vitellogenin and triglyceride between laying hens and sexually immature females. In contrast, reproductive state does not affect the relative weights of the liver, kidneys, spleen, and gizzard.