Identifying and Predicting Delayed Mortality with Toxicokinetic-Toxicodynamic Models

The prevalence of standardized toxicity testing in ecotoxicology has largely obscured the notion that toxicity is a function of time as well. The necessity of considering time is vividly demonstrated by observations of delayed mortality, that is, deaths continue to occur even when animals are no longer exposed to a toxicant. In this contribution, I explore to what extent toxicokinetic-toxicodynamic (TKTD) models from the framework of the General Unified Threshold model for Survival (GUTS) can capture delayed mortality, and to what extent this phenomenon can be predicted from short-term standard tests. I use a previously published data set for fluoroquinolones in Daphnia magna that shows strongly delayed mortality (using immobilization as a proxy for death). The model analysis shows that the GUTS stochastic death models can capture delayed mortality in the complete data set with a long recovery phase, but that the delayed effects would not have been predicted from a 2-day standard test. The study underlines the limited information content of standard acute test designs. Toxicokinetic-toxicodynamic modeling offers a handle on the time aspects of toxicity but cannot always be relied on to provide accurate extrapolations based on severely limited standard tests. The phenomenon of delayed toxicity requires more structured study to clarify its prevalence and impact; I discuss several avenues for further investigation. Environ Toxicol Chem 2024;00:1-6. © 2024 SETAC.

Changes in Temperature Alter the Toxicity of Polycyclic Aromatic Compounds to American Lobster (Homarus americanus) Larvae

Polycyclic aromatic compounds (PACs) present in the water column are considered to be one of the primary contaminant groups contributing to the toxicity of a crude oil spill. Because crude oil is a complex mixture composed of thousands of different compounds, oil spill models rely on quantitative structure-activity relationships like the target lipid model to predict the effects of crude oil exposure on aquatic life. These models rely on input provided by single species toxicity studies, which remain insufficient. Although the toxicity of select PACs has been well studied, there is little data available for many, including transformation products such as oxidized hydrocarbons. In addition, the effect of environmental influencing factors such as temperature on PAC toxicity is a wide data gap. In response to these needs, in the present study, Stage I lobster larvae were exposed to six different understudied PACs (naphthalene, fluorenone, methylnaphthalene, phenanthrene, dibenzothiophene, and fluoranthene) at three different relevant temperatures (10, 15, and 20 °C) all within the biological norms for the species during summer when larval releases occur. Lobster larvae were assessed for immobilization as a sublethal effect and mortality following 3, 6, 12, 24, and 48 h of exposure. Higher temperatures increased the rate at which immobilization and mortality were observed for each of the compounds tested and also altered the predicted critical target lipid body burden, incipient median lethal concentration, and elimination rate. Our results demonstrate that temperature has an important influence on PAC toxicity for this species and provides critical data for oil spill modeling. More studies are needed so oil spill models can be appropriately calibrated and to improve their predictive ability. Environ Toxicol Chem 2023;42:2389-2399. © 2023 SETAC.

Pathway network-based quantitative modeling of the time-dependent and dose-response anti-inflammatory effect of Reduning Injection

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional Chinese medicine (TCM) has extensive healing effects on inflammatory diseases with few side effects. Reduning injection (RDNI), a TCM prescription composed of Lonicera japonica Thunb., Gardenia jasminoides Ellis. and Artemisia annua L., is wildly used for treating inflammatory diseases. However, the mechanism of action of RDNI, like most TCM prescriptions, is unclear due to the complexity of relationships between components and their curative effects.

AIM OF THE STUDY

To develop a universal systems pharmacology protocol for mechanism modeling of TCM and apply it to reveal the real-time anti-inflammatory effect of Reduning Injection.

MATERIALS AND METHODS

Combined with database mining and references, a regulatory mechanism network of inflammation was constructed. A quantitative model was established afterwards by integrating pharmacokinetic data and two network parameters, namely Network Efficiency and Network Flux. The time-dependent and dose-response relationship of RDNI on the regulation of inflammation was then quantitatively evaluated. ELISA tests were performed to verify the reliability of the model.

RESULTS

Three cytokines, namely IL-1β, IL-6 and TNF-α were screened out to be markers for evaluation of the anti-inflammatory effect of RDNI. An HPLC method for the simultaneous determination of 10 RDNI compounds in SD rat plasma was established and then applied to pharmacokinetic studies. Based on compound activity and pharmacokinetic data, the time-dependent effect of RDNI were quantitatively predicted by the pathway network-based modeling procedure.

CONCLUSIONS

The quantitative model established in this work was successfully applied to predict a TCM prescription's real-time dynamic healing effect after administration. It is qualified to provide novel insights into the time-dependent and dose-effect relationship of drugs in an intricate biological system and new strategies for investigating the detailed molecular mechanisms of TCM.

Investigating Population-Level Toxicity of the Antidepressant Citalopram in Harpacticoid Copepods Using In Vivo Methods and Bioenergetics-Based Population Modeling

Recent research has revealed various lethal and sublethal effects of the selective serotonin reuptake inhibitor citalopram hydrobromide on the harpacticoid copepod Nitocra spinipes. In the present study, an individual-based model (IBM) grounded in the dynamic energy budget (DEB) theory was developed to extrapolate said effects to the population level. Using a generic DEB-IBM as a template, the model was designed to be as simple as possible, keeping model components that are outside the scope of the core DEB theory to a minimum. To test the model, a 56-day population experiment was performed at 0, 100, and 1000 μg citalopram hydrobromide L^-1 . In the experiment, the populations quickly reached a plateau in the control and at 100 μg L^-1 , which was correctly reproduced by the model and could be explained by food limitations hindering further population growth. At 1000 μg L^-1 , a clear mismatch occurred: Whereas in the experiment the population size increased beyond the supposed (food competition-induced) capacity, the model predicted a suppression of the population size. It is assumed that the IBM still misses important components addressing population density-regulating processes. Particularly crowding effects may have played an important role in the population experiment and should be further investigated to improve the model. Overall, the current DEB IBM for N. spinipes should be seen as a promising starting point for bioenergetics-based copepod population modeling, which-with further improvements-may become a valuable individual-to-population extrapolation tool in the future. Environ Toxicol Chem 2023;42:1094-1108. © 2023 SETAC.

Calibration of an acute toxicity model for the marine crustacean, Artemia franciscana, nauplii to support oil spill effect assessments

Oil spill risk and impact assessments rely on time-dependent toxicity models to predict the hazard of the constituents that comprise crude oils and petroleum substances. Dissolved aromatic compounds (ACs) are recognized as a primary driver of aquatic toxicity in surface spill exposure scenarios. However, limited time-dependent toxicity data are available for different classes of ACs to calibrate such models. This study examined the acute toxicity of 14 ACs and 3 binary AC mixtures on Artemia franciscana nauplii at 25 °C. Toxicity tests for 3 ACs were also conducted at 15 °C to evaluate the role of temperature on toxicity. The ACs investigated represented parent and alkylated homocyclic and nitrogen-, sulfur- and oxygen-containing heterocyclic structures with octanol-water partition coefficients (log K_ow) ranging from 3.2 to 6.6. Passive dosing was used to expose and maintain concentrations in toxicity tests which were confirmed using fluorometry, and independently validated for 6 ACs using GC-MS analysis. Mortality was assessed at 6, 24, and 48 h to characterize the time course of toxicity. No mortality was observed for the most hydrophobic AC tested, 7,12-dimethylbenz[a]anthracene, due to apparent water solubility constraints. Empirical log LC₅₀ s for the remaining ACs were fit to a linear regression with log K_ow to derive a critical target lipid body burden (CTLBB) based on the target lipid model. The calculated 48 h CTLBB of 47.1 ± 8.1 μmol/g octanol indicates that Artemia nauplii exhibited comparable sensitivity to other crustaceans. A steep concentration-response was found across all compounds as evidenced by a narrow range (1.0-3.1) in the observed LC₅₀ /LC₁₀ ratio. Differences in toxicokinetics were noted, and no impacts of temperature-dependence of AC toxicity were found. Toxicity data obtained for individual ACs yielded acceptable predictions of observed binary AC mixture toxicity. Results from this study advance toxicity models used in oil spill assessments.

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.

Transcriptomic Points of Departure Calculated from Rainbow Trout Gill, Liver, and Gut Cell Lines Exposed to Methylmercury and Fluoxetine

Ethical and resource limitation concerns are pushing chemicals management to develop alternatives to animal testing strategies. The objective of our study was to determine whether transcriptomic point of departure (tPOD) values could be derived from studies that followed Organisation for Economic Co-operation and Development (OECD) Test No. 249 (rainbow trout gill cell line), as well as from studies on trout liver and gut cells. Gill, liver, and gut cell lines were exposed to methylmercury and fluoxetine. Concentrations causing 50% cytotoxicity (LC50) were derived, the whole transcriptome was sequenced, and gene tPOD and pathway benchmark dose (BMD) values were derived from transcriptomic dose-response analysis. Differences in LC50 and transcriptomic responses across the cell lines were noted. For methylmercury, the tPOD_mode values were 14.5, 20.5, and 17.8 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway (pathway BMDs in parentheses) was ferroptosis in the gill (3.1 ppb) and liver (3.5 ppb), and glutathione metabolism in the gut (6.6 ppb). For fluoxetine, the tPOD_mode values were 109.4, 108.4, and 97.4 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway was neurotrophin signaling in the gill (147 ppb) and dopaminergic signaling in the gut (86.3 ppb). For both chemicals, the gene tPOD and pathway BMD values were lower than cytotoxic concentrations in vitro, and within 10-fold below the in vivo LC50s. By bringing together transcriptomics and dose-response analysis with an OECD test method in three cell lines, the results help to establish an in vitro method yielding tPOD values that are hypothesized to be protective of in vivo concentrations associated with adverse outcomes, and also give insights into mechanisms of action. Environ Toxicol Chem 2022;41:1982-1992. © 2022 SETAC.

Making Sense of Life-History Effects of the Antidepressant Citalopram in the Copepod Nitocra spinipes Using a Bioenergetics Model

The global consumption of human antidepressants has steadily increased over the last years. The most widely prescribed antidepressants are the selective serotonin reuptake inhibitors (SSRIs), which have been linked to various life-history effects in nontarget organisms. We investigated the effects of the SSRI citalopram hydrobromide on the life history of the copepod Nitocra spinipes. Slight but significant developmental delay effects were observed at nominal concentrations of 0.1 and 1 µg/L, with stronger effects occurring at measured concentrations of 178 µg/L and above. At 77 µg/L and above, a significant increase in adult body length and offspring production/brood was found, although the time between brood releases remained unaffected. The pre-adult surviving fraction was significantly reduced (by 44%) at 765 µg/L. For a mechanistic evaluation of these observations, we used a bioenergetics model for N. spinipes based on the dynamic energy budget theory. Toxicokinetic and toxicodynamic submodels were used to dynamically simulate the chemical uptake and elimination, as well as dose-response relationships for hypothetical physiological modes of action and survival over time. Although none of the commonly invoked physiological modes of action, acting on assimilation, maintenance, growth, or offspring production, could explain the observed combination of effects, a newly proposed physiological mode of action acting on the process of maturation delivered correct predictions in terms of each effect's direction. The model fits could be further improved by allowing for a gentler concentration-effect slope and by adding an auxiliary physiological mode of action acting on the reproduction efficiency. The quantitative explanations provided in the present study offer a starting point for exploratory simulation studies investigating the effects of SSRIs at higher ecological levels. Environ Toxicol Chem 2021;40:1928-1939. © 2021 SETAC.

Assessing the Toxicity of Individual Aromatic Compounds and Mixtures to American Lobster (Homarus americanus) Larvae Using a Passive Dosing System

Aquatic exposures to aromatic compounds (ACs) may be important contributors to biological effects of oil spills. The present study examined the acute toxicity of 11 ACs and 3 binary AC mixtures on stage 1 American lobster larvae using a passive dosing test design. The ACs investigated covered a range of classes and log octanol-water partition coefficient values (K_OW ; 2.5-5.5). Silicone O-rings were used to partition ACs into seawater and maintain stable exposures. Exposed lobster larvae were assessed for mobility and survival at 3, 6, 12, 24, 36, and 48 h. Fluorometry and gas chromatography-mass spectrometry measurements confirmed well-defined substance exposures. Expressing lethality in terms of chemical activities yielded values between 0.01 and 0.1, consistent with a baseline mode of action. Analysis of time-dependent median lethal/effect concentration (L/EC50) values were used to determine incipient values. An expected linear relationship between the incipient log L/EC50 and log K_OW was fit to the empirical toxicity data to derive critical target lipid body burdens for immobilization and lethality endpoints. These values indicate that American lobster larvae fall on the sensitive end of the acute species sensitivity distribution. We used AC toxicity data to successfully predict toxicity of binary mixtures assuming additive toxicity. The observed time-dependent toxicity was inversely related to log K_OW and occurred more quickly than reported previously. The results contribute to improving models for predicting oil spill impacts on American lobster larvae populations. Environ Toxicol Chem 2021;40:1379-1388. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Impurities in technical mixtures of chlorinated paraffins show AhR agonist properties as determined by the DR-CALUX bioassay

Chlorinated paraffins (CPs) are produced at more than one million tons per year. Technical CPs mixtures may contain impurities, which end up in consumer products. In the present study, 17 technical CPs mixtures were investigated for the potential occurrence of potential impurities. By applying the DR-CALUX bioassay, 3 out of 17 technical mixtures were shown to elicit responses at 4 h exposure time, but much lower at 48 h. Constitutional defined CPs materials did not show responses. Subsequently different groups of known AhR-agonists and compounds suspected to be present in technical CPs mixtures were investigated. Benzene, (poly)chlorobenzene, non-dioxin like polychlorinated naphthalenes (PCNs), and three-ringed polyaromatic hydrocarbons (PAHs) did not result in a significant response at 4 h or 48 h. TCDD, non-ortho PCBs, dioxin-like PCNs, four or five ringed PAHs and their chlorinated analogues resulted in a significant response. TCDD and the non-ortho PCBs showed the highest potency and stability, while dioxin-like PCNs, PAHs, and the chlorinated PAHs were clearly inactivated (metabolized) at longer incubation. Altogether, the present findings substantiate that AhR-mediated responses of CPs technical mixtures in the DR-CALUX bioassay are caused by impurities, most likely some intermediate stable AhR-agonists such as dioxin-like PCNs or (chlorinated) PAHs. The current study shows that impurities in CPs technical mixtures need to be investigated for assessing the safety of technical CPs mixtures.

Estrogenic activity of food contact materials-evaluation of 20 chemicals using a yeast estrogen screen on HPTLC or 96-well plates

Food contact materials (FCM) may contain complex mixtures of estrogenic chemicals. A yeast estrogen screen performed on high performance thin-layer chromatography plates (planar-YES, P-YES) is promising for analysis of such mixtures, as it could allow for better elucidation of effects compared with established methods in microtiter plates. However, the P-YES has not been directly compared with established methods. We compared the performance of a microtiter plate YES (lyticase-YES, L-YES) to P-YES on silica gel HPTLC plates using 17β-estradiol (E2), 20 chemicals representative of migrants from plastic FCM, and three migrates of coated metal food cans. Effective doses (ED10, ED50) and estradiol equivalencies were calculated for each chemical. Thirteen chemicals had calculable EDs in the L-YES or P-YES, with average EDs 13-fold (range 0.63-36) more potent in P-YES than in the L-YES. Normalized to E2, the median estrogenicity was within 1.5-fold (0.43-8.8) between the assays. Therefore, P-YES was as or more sensitive than L-YES but potencies relative to E2 were comparable between assays. With chromatography, the P-YES detected estrogenicity in coated metal cans, effects that were unmeasurable in L-YES. With the sample preparation methods used in this study, both YES assays are sufficiently sensitive to detect bisphenol A below the specific migration limit for plastic packaging (0.05 mg/kg food). This study demonstrates that P-YES outperforms L-YES because it is more sensitive, provides comparable estradiol equivalents, and circumvents confounding mixture effects. The P-YES will be useful for routine monitoring of FCM and toxicant identification in problematic materials. Graphical abstract.

Risk-based principles and incompleteness theorems for linear dose-response extrapolation for carcinogenic chemicals

To conduct better health risk assessments, this study introduced two risk-based principles and a series of line-lognormal-intersection theorems that helped derive the safe ranges of the cancer slope factors (CSFs) for 708 carcinogenic chemicals. The extrapolated linear dose-response relationships presented in this study can ensure safety with respect to both static and dose-based instantaneous risks compared to the lognormal dose-response model. The theorems proved that the maximum static and dose-based hazard risk ratios of a lognormal curve and a linear model are independent of a chemical's toxicity (the effect dose that corresponds to a 50% response, or ED₅₀), where the selected linear extrapolation (m value) and the individual variability (σ) of the responses to carcinogenic chemicals are two determining factors. The theorems also indicated that individual variability determines the range of m if the acceptable risk ratios were regulated. When σ was 1.36 (i.e., the 50th percentile of the individual variability's lognormal distribution), the safe range of m was derived as [11.22, 21.46] (i.e., from ED_11.22 to ED_21.46); if the 95th percentile of the σ lognormal distribution was used, the safe range of m was [1.13, 4.57] (i.e., from ED_1.13 to ED_4.57). This study also showed that for a relatively homogenous population (i.e., σ is relatively small) that has similar characteristics, the linear dose-response extrapolation method might not be completely effective due to the shape shift of the lognormal curve that draws the static risk of the extrapolated linear model away from the lognormal model.

Influence of Soil Phosphate on the Accumulation and Toxicity of Arsenic and Antimony in Choy Sum Cultivated in Individually and Co-contaminated Soils

Fertilizers containing phosphate (PO₄^3- ) are commonly used within the agricultural industry and are known to increase the bioavailability and mobility of metalloids like arsenic (As). This may increase plant uptake of As and hence pose a risk to human health. Arsenic and antimony (Sb) often co-occur in contaminated soils; however, little is known about the interactions between As and Sb with PO₄^3- on their bioavailability, accumulation, and toxicity in plants. The present study investigated individual and combined As and Sb-contaminated soils across 2 soil PO₄^3- concentrations using a commonly consumed leafy vegetable, choy sum (Brassica chinensis var. parachinensis). Increased soil PO₄^3- had no clear influence on the bioavailability of As or Sb (derived from a sequential extraction procedure). At high PO₄^3- concentration, B. chinensis accumulated higher amounts of As in the shoots and roots in both individual and co-contaminated soil, whereas Sb accumulation increased only when Sb was the only contaminant. When As was the only contaminant, the translocation of As from roots to shoots decreased as soil PO₄^3- increased. Increased soil PO₄^3- had no influence on Sb translocation from root to shoot. Although As was toxic (impaired growth) at low PO₄^3- soil concentration, no toxicity was observed in the high-PO₄^3- soil. No toxicity was observed for Sb in either low- or high-PO₄^3- soils. Increased soil PO₄^3- concentration ameliorated or masked As toxicity to plant growth and led to higher As concentration in the plant's edible parts. The addition of high soil PO₄^3- concentrations ameliorated or masked As toxicity to plant growth in both individually and As + Sb co-contaminated soil; however, the plant's edible parts accumulated higher As and Sb concentrations. Environ Toxicol Chem 2020;39:1233-1243. © 2020 SETAC.

Soil exposed to silver nanoparticles reveals significant changes in community structure and altered microbial transcriptional profiles

Anthropogenic activities can disrupt soil ecosystems, normally resulting in reduced soil microbial health. Regulatory agencies need to determine the effects of uncharacterized substances on soil microbial health to establish the safety of these chemicals if they end up in the environment. Previous work has focused on measuring traditional ecotoxicologial endpoints within the categories of microbial biomass, activity, and community structure/diversity. Because these tests can be labor intensive, lengthy to conduct, and cannot measure changes in individual gene functions, we wanted to establish whether metatranscriptomics could be used as a more sensitive endpoint and provide a perspective on community function that is more informative than taxonomic identification of microbes alone. We spiked a freshly collected sandy loam soil (Vulcan, Alberta, Canada) with 0, 60, 145, 347, 833, and 2000 mg kg^-1 of silver nanoparticles (AgNPs), a known antagonist of microorganisms due to its propensity for dissolution of toxic silver ions. Assessments performed in our previous work using traditional tests demonstrated the toxicity of AgNPs on soil microbial processes. We expanded this analysis with genomics-based tests by measuring changes in community taxonomic structure and function using 16S rDNA profiling and metatranscriptomics. In addition to identifying bacterial taxa affected by AgNPs, we found that genes involved in heavy metal resistance (e.g., the CzcA efflux pump) and other toxicity response pathways were highly upregulated in the presence of silver. Dose-response analysis using BMDExpress2 software successfully modeled many physiologically relevant genes responding to low concentrations of AgNPs. We found that the transcriptomic point of departure (BMD₅₀) was lower than the IC50s calculated using the traditional tests in our previous work. These results suggest that dose-response modeling of metatranscriptomic gene expression is a useful tool in soil microbial health assessment. SUMMARY: Genomics-based endpoints for the assessment of soil microbial health can be used to perform quantitative dose-response modeling, and soil-based RNAseq adds functional insights.

Importance of Surface Coating to Accumulation Dynamics and Acute Toxicity of Copper Nanomaterials and Dissolved Copper in Daphnia magna

We evaluated the effect of copper oxide nanomaterials (CuO NMs), uncoated and with 3 different surface coatings (carboxylated, pegylated, and ammonia groups), on acute toxicity and accumulation dynamics in Daphnia magna. With the use of biodynamic modelling, biosorption and elimination rate constants were determined for D. magna following waterborne exposure to dissolved Cu and CuO NMs. The relationship between modeled parameters and acute toxicity endpoints was evaluated to investigate whether accumulation dynamics parameters could be used as a predictor of acute toxicity. The Langmuir equation was used to characterize the biosorption dynamics of Cu NMs and Cu chloride, used as dissolved Cu control. Uptake rates showed the following NM rankings: pristine-CuO > NH₃ -CuO > aqueous Cu > polyethylene glycol (PEG)-CuO > COOH-CuO. To determine Cu elimination by D. magna, a one-compartment model was used. Different elimination rate constants were estimated for each chemical substance tested. Those that were easily biosorbed were also easily removed from organisms. Biosorption and depuration properties of NMs were correlated with zeta potential values and diameters of NM agglomerates in the suspensions. No link was found between biosorption and toxicity. Waterborne exposures to more difficult-to-biosorb CuO NMs were more likely to induce adverse effects than those that biosorbed easily. It is proposed that some physicochemical properties of NMs in media, including zeta potential and agglomerate diameter, can lead to higher biosorption but do not necessarily affect toxicity. The mode of interaction of the NMs with the organism seems to be complex and to depend on chemical speciation and physicochemical properties of the NMs inside an organism. Moreover, our findings highlight that coating type affects the biosorption dynamics, depuration kinetics, and dissolution rate of NMs in media. Environ Toxicol Chem 2020;39:287-299. © 2019 SETAC.

Ordinal dose-response modeling approach for the phthalate syndrome

BACKGROUND

The phthalate syndrome (PS) is a collection of related male reproductive developmental effects, ranging in severity, that have been observed in rats after gestational exposure to developmentally-toxic phthalates. For statistical purposes, the PS is defined as a single endpoint and one dose-response analysis is conducted, rather than conducting multiple analyses on each individual endpoint.

OBJECTIVE

To improve dose-response modeling approaches for the PS and other syndromes of effects by accounting for differing severity levels among the endpoints.

METHODS

Ordinal dose-response modeling was performed on PS data from a published study of diisobutyl phthalate (DIBP) gestational exposure to male Sprague-Dawley rats. To incorporate PS endpoint severity, the endpoints were categorized into ordinal levels based on the expected impact of male developmental endpoint's on fertility. Then, a benchmark dose was estimated for each ordinal level. A bootstrap procedure was used to account for the nested nature of the data, and a sensitivity analysis was performed to assess the bootstrap results. A comparison of the estimates between the ordinal and the dichotomous model was performed.

RESULTS

The ordinal version of the log-logistic model applied to the data categorized by PS endpoint severity level provided benchmark dose estimates that were closer to each other in value and had lower variability than the traditional dichotomous application. The sensitivity analysis confirmed the validity of the bootstrap results.

CONCLUSION

The ordinal dose-response modeling method accounts for severity differences among dichotomous PS endpoints, can be expanded in the future to include more severity levels, and can be used in both single and cumulative phthalate risk assessments.

Long-term effects and recovery of streams from acid mine drainage and evaluation of toxic metal threshold ranges for macroinvertebrate community reassembly

Monitoring of benthic invertebrates in streams receiving acidic metal-contaminated water over an 18-yr period revealed both degraded conditions and recovery along a network of downstream locations. Compared with reference streams, and over the course of clean-up remediation efforts below an abandoned open-pit sulfur mine in the central Sierra Nevada of California, improving water quality was accompanied by recovery of benthic communities at some sites. Years of high flow resulted in degraded biological status when acid mine drainage capture was incomplete and metal loading had increased with runoff. Seasonal patterns of recovery evident in the fall after the summer treatment season reverted in the next spring after overwinter periods when sources were not captured. As the metal load has been reduced, phased recovery of community structure, function, and similarity progressed toward that of reference assemblage taxonomic composition. From impacted communities dominated by relatively tolerant midges, reassembly involved an increase in density, return of long-lived taxa, an increased ratio of sensitive-to-tolerant forms, then overall diversity and community composition, and eventually large predators and grazers reappearing along with mayfly, stonefly, and caddisfly richness. Threshold effect levels defined using several analysis methods showed that the response range of biological indicators corresponds to US Environmental Protection Agency guidelines of predicted effects utilizing cumulative criterion units (CCUs) of metal toxicity (i.e., CCU ∼ 1). All sites have shown improved function with increased density of some or all trophic groups over time. Although recovery is progressing, year-around treatment may be necessary to fully restore biological integrity in streams nearest the mine. Environ Toxicol Chem 2018;37:2575-2592. © 2018 SETAC.

The advantages of linear concentration-response curves for in vitro bioassays with environmental samples

In vitro assays and high-throughput screening (HTS) tools are increasingly being employed as replacements for animal testing, but most concentration-response curves are still evaluated with models developed for animal testing. We argue that application of in vitro assays, particularly reporter gene assays, to environmental samples can benefit from a different approach to concentration-response modeling. First, cytotoxicity often occurs at higher concentrations, especially for weakly acting compounds and in complex environmental mixtures with many components. In these cases, specific effects can be masked by cytotoxicity. Second, for many HTS assays, low effect levels can be precisely quantified because of the low variability of controls in cell-based assays and the opportunity to run many concentrations and replicates when using high-density well-plate formats (e.g., 384 or more wells per plate). Hence, we recommend focusing concentration-response modeling on the lower portion of the concentration-response curve, which is approximately linear. Effect concentrations derived from low-effect level linear concentration-response models facilitate simple derivation of relative effect potencies and the correct application of mixture toxicity models in the calculation of bioanalytical equivalent concentrations. Environ Toxicol Chem 2018;37:2273-2280. © 2018 SETAC.

Parameterizing the binding properties of dissolved organic matter with default values skews the prediction of copper solution speciation and ecotoxicity in soil

Parameterizing speciation models by setting the percentage of dissolved organic matter (DOM) that is reactive (% r-DOM) toward metal cations at a single 65% default value is very common in predictive ecotoxicology. The authors tested this practice by comparing the free copper activity (pCu²⁺  = -log₁₀ [Cu²⁺ ]) measured in 55 soil sample solutions with pCu²⁺ predicted with the Windermere humic aqueous model (WHAM) parameterized by default. Predictions of Cu toxicity to soil organisms based on measured or predicted pCu²⁺ were also compared. Default WHAM parameterization substantially skewed the prediction of measured pCu²⁺ by up to 2.7 pCu²⁺ units (root mean square residual = 0.75-1.3) and subsequently the prediction of Cu toxicity for microbial functions, invertebrates, and plants by up to 36%, 45%, and 59% (root mean square residuals ≤9 %, 11%, and 17%), respectively. Reparametrizing WHAM by optimizing the 2 DOM binding properties (i.e., % r-DOM and the Cu complexation constant) within a physically realistic value range much improved the prediction of measured pCu²⁺ (root mean square residual = 0.14-0.25). Accordingly, this WHAM parameterization successfully predicted Cu toxicity for microbial functions, invertebrates, and plants (root mean square residual ≤3.4%, 4.4%, and 5.8%, respectively). Thus, it is essential to account for the real heterogeneity in DOM binding properties for relatively accurate prediction of Cu speciation in soil solution and Cu toxic effects on soil organisms. Environ Toxicol Chem 2017;36:898-905. © 2016 SETAC.

An automated approach to Litchfield and Wilcoxon's evaluation of dose-effect experiments using the R package LW1949

The authors developed a package, LW1949, for use with the statistical software R to automatically carry out the manual steps of Litchfield and Wilcoxon's method of evaluating dose-effect experiments. The LW1949 package consistently finds the best fitting dose-effect relation by minimizing the chi-squared statistic of the observed and expected number of affected individuals and substantially speeds up the line-fitting process and other calculations that Litchfield and Wilcoxon originally carried out by hand. Environ Toxicol Chem 2016;35:3058-3061. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US Government work and, as such, is in the public domain in the United States of America.

The acute toxicity of major ion salts to Ceriodaphnia dubia: I. influence of background water chemistry

The ions Na⁺ , K⁺ , Ca²⁺ , Mg²⁺ , Cl^- , SO₄^2- , and HCO₃^- /CO₃^2- (referred to in the present study as "major ions") are present in all freshwaters and physiologically required by aquatic organisms but can increase to harmful levels from a variety of anthropogenic activities. It is also known that the toxicities of major ion salts can vary depending on the concentrations of other ions, and understanding these relationships is key to establishing appropriate environmental limits. The authors present a series of experiments with Ceriodaphnia dubia to evaluate the acute toxicity of 12 major ion salts and to determine how toxicity of these salts varies as a function of background water chemistry. All salts except CaSO₄ and CaCO₃ were acutely toxic below saturation, with the lowest median lethal concentrations found for K salts. All 10 salts that showed toxicity also showed some degree of reduced toxicity as the ionic content of the background water increased. Experiments that independently varied Ca:Mg ratio, Na:K ratio, Cl:SO₄ ratio, and alkalinity/pH demonstrated that Ca concentration was the primary factor influencing the toxicities of Na and Mg salts, whereas the toxicities of K salts were primarily influenced by the concentration of Na. These experiments also indicated multiple mechanisms of toxicity and suggested important aspects of dosimetry; the toxicities of K, Mg, and Ca salts were best related to the chemical activity of the cation, whereas the toxicities of Na salts also reflected an influence of the anions and were well correlated with osmolarity. Understanding these relationships between major ion toxicity and background water chemistry should aid in the development of sensible risk-assessments and regulatory standards. Environ Toxicol Chem 2016;35:3039-3057. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Comparison and evaluation of pesticide monitoring programs using a process-based mixture model

A number of European countries run large-scale pesticide monitoring schemes in watersheds aimed at identifying and evaluating the presence of pesticide residues in the environment. These schemes provide national and regional scale assessments of pesticide concentrations within the context of environmental quality assessment, aiming to ensure some degree of ecological protection. The present study is aimed at evaluating the joint effects of the pesticide mixtures detected in monitoring programs, using a process-based mixture model that was parameterized for Daphnia magna. In total, over 15 000 samples containing over 1 million individual measurements were evaluated for effects. It was found that there are only a small number of places where one can expect to have effects on daphnids, based on measured concentrations. The most polluted samples would cause extinction of a daphnid population within only 30 h. The results show that effects are mostly triggered by a limited number of pesticide residues at locations with high emissions. It was also shown that the analytical detection limits are basically too high to exclude mixture effects. So, despite all the effort that is put into chemical monitoring programs, it remains a challenge to make statements on whether or not the environment is protected. Recommendations are offered for a different setup of monitoring programs to improve this situation. Environ Toxicol Chem 2016;35:3113-3123. © 2016 SETAC.

An analysis of lethal and sublethal interactions among type I and type II pyrethroid pesticide mixtures using standard Hyalella azteca water column toxicity tests

A novel 2-tiered analytical approach was used to characterize and quantify interactions between type I and type II pyrethroids in Hyalella azteca using standardized water column toxicity tests. Bifenthrin, permethrin, cyfluthrin, and lambda-cyhalothrin were tested in all possible binary combinations across 6 experiments. All mixtures were analyzed for 4-d lethality, and 2 of the 6 mixtures (permethrin-bifenthrin and permethrin-cyfluthrin) were tested for subchronic 10-d lethality and sublethal effects on swimming motility and growth. Mixtures were initially analyzed for interactions using regression analyses, and subsequently compared with the additive models of concentration addition and independent action to further characterize mixture responses. Negative interactions (antagonistic) were significant in 2 of the 6 mixtures tested, including cyfluthrin-bifenthrin and cyfluthrin-permethrin, but only on the acute 4-d lethality endpoint. In both cases mixture responses fell between the additive models of concentration addition and independent action. All other mixtures were additive across 4-d lethality, and bifenthrin-permethrin and cyfluthrin-permethrin were also additive in terms of subchronic 10-d lethality and sublethal responses. Environ Toxicol Chem 2016;35:2542-2549. © 2016 SETAC.

Challenges in radiobiological modeling: can we decide between LQ and LQ-L models based on reviewed clinical NSCLC treatment outcome data?

Aim

To study the dose-response of stage I non-small-cell lung cancer (NSCLC) in terms of long-term local tumor control (LC) after conventional and hypofractionated photon radiotherapy, modeled with the linear-quadratic (LQ) and linear-quadratic-linear (LQ-L) approaches and to estimate the clinical α/β ratio within the LQ frame.

Material and methods

We identified studies of curative radiotherapy as single treatment through MedLine search reporting 3-year LC as primary outcome of interest. Logistic models coupled with the biologically effective dose (BED) at isocenter and PTV edge according to both the LQ and LQ-L models with α/β = 10 Gy were fitted. Additionally, α/β was estimated from direct LQ fits.

Results

Thirty one studies were included reporting outcome of 2319 patients. The LQ-L fit yielded a significant value of 11.0 ± 5.2 Gy for the dose threshold (D_t) for BED₁₀ at the isocenter. The LQ and LQ-L fits did not differ substantially. Concerning the estimation of α/β, the value obtained from the direct LQ fit for the complete fractionation range was 3.9 [68 % CI: 2.2–9.0] Gy (p > 0.05).

Conclusion

Both LQ and LQ-L fits can model local tumor control after conventionally and hypofractionated irradiation and are robust methods for predicting clinical effects. The observed dose-effect for local control in NSCLC is weaker at high doses due to data dispersion. For BED₁₀ values of 100–150 Gy in ≥3 fractions, the differences in isoeffects predicted by both models can be neglected.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-016-0643-5) contains supplementary material, which is available to authorized users.

Critical perspectives on mercury toxicity reference values for protection of fish

Environmental management decisions at mercury-contaminated sediment sites are predicated on the understanding of risks to various receptors, including fish. Toxicity reference values (TRVs) for interpreting risks to fish have been developed to assess mercury concentrations in fish or fish prey. These TRVs were systematically evaluated based on several lines of evidence. First, their conceptual basis and specific derivation were evaluated, including a close review of underlying toxicity studies. Second, case studies were reviewed to investigate whether TRVs are predictive of effects on fish populations in the field. Third, TRVs were compared with available information regarding preindustrial and present-day background concentrations of mercury in fish. The findings show that existing TRVs are highly uncertain, because they were developed using limited data from studies not designed for TRV derivation. Although field studies also entail uncertainty, several case studies indicate no evidence of adverse effects despite mercury exposures that exceed the available TRVs. Some TRVs also fall within the range of background mercury concentrations in predatory or prey fish. Lack of information on the selenium status of mercury-exposed fish is a critical confounding factor, and the form of methylmercury used in toxicity testing may also contribute to differences between TRV-based predictions and field observations of mercury effects on fish. On balance, the available information indicates that several of the TRVs reviewed are lower than necessary to protect fish populations. The 20% effect concentration from a previously published dose-response analysis appears closer to an effect threshold, based on available laboratory data. Additional research is needed to provide a stronger basis to establish dose-response relationships for mercury effects on fish.

High-throughput concentration-response analysis for omics datasets

Omics-based methods are increasingly used in current ecotoxicology. Therefore, a large number of observations for various toxic substances and organisms are available and may be used for identifying modes of action, adverse outcome pathways, or novel biomarkers. For these purposes, good statistical analysis of toxicogenomic data is vital. In contrast to established ecotoxicological techniques, concentration-response modeling is rarely used for large datasets. Instead, statistical hypothesis testing is prevalent, which provides only a limited scope for inference. The present study therefore applied automated concentration-response modeling for 3 different ecotoxicotranscriptomic and ecotoxicometabolomic datasets. The modeling process was performed by simultaneously applying 9 different regression models, representing distinct mechanistic, toxicological, and statistical ideas that result in different curve shapes. The best-fitting models were selected by using Akaike's information criterion. The linear and exponential models represented the best data description for more than 50% of responses. Models generating U-shaped curves were frequently selected for transcriptomic signals (30%), and sigmoid models were identified as best fit for many metabolomic signals (21%). Thus, selecting the models from an array of different types seems appropriate, because concentration-response functions may vary because of the observed response type, and they also depend on the compound, the organism, and the investigated concentration and exposure duration range. The application of concentration-response models can help to further tap the potential of omics data and is a necessary step for quantitative mixture effect assessment at the molecular response level.

Sensitivity of shovelnose sturgeon (Scaphirhynchus platorynchus) and pallid sturgeon (S. albus) early life stages to 3,3',4,4',5-pentachlorobiphenyl and 2,3,7,8-tetrachlorodibenzo-P-dioxin exposure

Concern exists that polychlorinated biphenyls (PCBs) may be contributing to the current decline of shovelnose sturgeon (Scaphirhynchus platorynchus) and the US federally endangered pallid sturgeon (Scaphirhynchus albus). Waterborne exposures with newly fertilized eggs were used to assess developmental and morphological effects of 2 of the most potent aryl hydrocarbon receptor (AhR) agonists, 3,3',4,4',5-pentachlorobiphenyl (PCB-126) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), on early life stage shovelnose and pallid sturgeon. No dose-related effects of PCB-126 were observed on percent development or hatch in either species at concentrations as high as 1711 ng/g egg. Effects of TCDD on percent development were not assessed in shovelnose sturgeon. However, percent development was not affected by TCDD in pallid sturgeon, and percent hatch was unaffected by TCDD doses as high as 60 ng/g egg to 81 ng/g egg in either species. Morphological pathologies such as yolk sac edema and craniofacial deformities were typical of AhR agonist exposure and were similar in both species. Calculated PCB-126 50% lethal dose (LD50, 95% fiducial limits) values were 196 ng/g egg (188-203 ng/g) for shovelnose and 159 ng/g egg (122-199 ng/g) for pallid sturgeon. Likewise, calculated TCDD LD50 values were 13 ng/g egg (11-15 ng/g) for shovelnose and 12 ng/g egg (10-14 ng/g) for pallid sturgeon. These LD50 values are among the highest recorded in early life stage fish, suggesting that early life stage Scaphirhynchus sturgeon may be comparatively insensitive to AhR agonists.

A practical approach to determine dose metrics for nanomaterials

Traditionally, administered mass is used to describe doses of conventional chemical substances in toxicity studies. For deriving toxic doses of nanomaterials, mass and chemical composition alone may not adequately describe the dose, because particles with the same chemical composition can have completely different toxic mass doses depending on properties such as particle size. Other dose metrics such as particle number, volume, or surface area have been suggested, but consensus is lacking. The discussion regarding the most adequate dose metric for nanomaterials clearly needs a systematic, unbiased approach to determine the most appropriate dose metric for nanomaterials. In the present study, the authors propose such an approach and apply it to results from in vitro and in vivo experiments with silver and silica nanomaterials. The proposed approach is shown to provide a convenient tool to systematically investigate and interpret dose metrics of nanomaterials. Recommendations for study designs aimed at investigating dose metrics are provided.

Selecting the best design for nonstandard toxicology experiments

Although many experiments in environmental toxicology use standard statistical experimental designs, there are situations that arise where no such standard design is natural or applicable because of logistical constraints. For example, the layout of a laboratory may suggest that each shelf serve as a block, with the number of experimental units per shelf either greater than or less than the number of treatments in a way that precludes the use of a typical block design. In such cases, an effective and powerful alternative is to employ optimal experimental design principles, a strategy that produces designs with precise statistical estimates. Here, a D-optimal design was generated for an experiment in environmental toxicology that has 2 factors, 16 treatments, and constraints similar to those described above. After initial consideration of a randomized complete block design and an intuitive cyclic design, it was decided to compare a D-optimal design and a slightly more complicated version of the cyclic design. Simulations were conducted generating random responses under a variety of scenarios that reflect conditions motivated by a similar toxicology study, and the designs were evaluated via D-efficiency as well as by a power analysis. The cyclic design performed well compared to the D-optimal design.