Cross-species extrapolation of chemical sensitivity

Environmentally relevant concentrations of the tricyclic antidepressant, amitriptyline, affect feeding and reproduction in a freshwater mollusc

Antidepressant drugs (ADDs) are one of the most extensively used pharmaceuticals globally. They act at particularly low therapeutic concentrations to modulate monoamine neurotransmission, which is one of the most evolutionary conserved pathways in both humans and animal species including invertebrates. As ADDs are widely detected in the aquatic environment at low concentrations (ng/L to low µg/L), their potential to exert drug-target mediated effects in aquatic species has raised serious concerns. Amitriptyline (AMI) is the most widely used tricyclic ADD, while monoamines, the target of ADDs, are major bioregulators of multiple key physiological processes including feeding, reproduction and behaviour in molluscs. However, the effects of AMI on feeding, reproduction and mating behaviour are unknown in molluscs despite their ecological importance, diversity and reported sensitivity to ADDs. To address this knowledge gap, we investigated the effects of environmentally relevant concentrations of AMI (0, 10, 100, 500 and 1000 ng/L) on feeding, reproduction and key locomotor behaviours, including mating, in the freshwater gastropod, Biomphalaria glabrata over a period of 28 days. To further provide insight into the sensitivity of molluscs to ADDs, AMI concentrations (exposure water and hemolymph) were determined using a novel extraction method. The Fish Plasma Model (FPM), a critical tool for prioritization assessment of pharmaceuticals with potential to cause drug target-mediated effects in fish, was then evaluated for its applicability to molluscs for the first time. Disruption of food intake (1000 ng/L) and reproductive output (500 and 1000 ng/L) were observed at particularly low hemolymph levels of AMI, whereas locomotor behaviours were unaffected. Importantly, the predicted hemolymph levels of AMI using the FPM agreed closely with the measured levels. The findings suggest that hemolymph levels of AMI may be a useful indicator of feeding and reproductive disruptions in wild population of freshwater gastropods, and confirm the applicability of the FPM to molluscs for comparative pharmaceutical hazard identification.

A systematic evidence map and bibliometric analysis of the behavioural impacts of pesticide exposure on zebrafish

Pesticides are indispensable in agriculture and have become ubiquitous in aquatic environments. Pesticides in natural environments can cause many negative impacts on aquatic species, ranging from mortality to sub-lethal physiological and behavioural changes. The complex sub-lethal impacts of pesticides are routinely tested on model species, with zebrafish (Danio rerio) being regularly used as a behavioural model. Although behavioural ecotoxicology research using zebrafish is increasing rapidly, we lack quantitative evidence to support which pesticides have been tested and how study designs are carried out. This shortcoming not only limits the deliberate planning for future primary studies to fill the knowledge gaps but also hinders evidence synthesis. To provide quantitative evidence of what pesticides are currently studied and what study designs are used, we combined a systematic evidence map approach and bibliometric analysis. This novel method has been coined research weaving and allows us to elicit gaps and clusters in our evidence base, whilst showing connections between authors and institutions. The methodology can be summarised in five primary steps: literature searching, screening, extraction, data analysis and bibliometric analysis. We identified four areas where research on the sub-lethal effects of pesticide exposure on zebrafish is lacking. First, some widely used pesticides, such as neonicotinoids, are understudied. Second, most studies do not report important elements of the study design, namely the sex and the life-stage of the zebrafish. Third, some behaviours, such as impacts of pesticide exposure on zebrafish cognition, are underexplored. And last, we revealed through the bibliometric analysis that most of the research is conducted in developed countries and there is limited cross country co-authorships. Upon identifying these gaps, we offer solutions for each limitation, emphasizing the importance of diverse global research output and cross-country co-authorships. Our systematic evidence map and bibliometric analysis provide valuable insights for helping to guide future research, which can be used to help support evidence-based policy decisions.

Climate warming shifts riverine macroinvertebrate communities to be more sensitive to chemical pollutants

Freshwaters are highly threatened ecosystems that are vulnerable to chemical pollution and climate change. Freshwater taxa vary in their sensitivity to chemicals and changes in species composition can potentially affect the sensitivity of assemblages to chemical exposure. Here we explore the potential consequences of future climate change on the composition and sensitivity of freshwater macroinvertebrate assemblages to chemical stressors using the UK as a case study. Macroinvertebrate assemblages under end of century (2080-2100) and baseline (1980-2000) climate conditions were predicted for 608 UK sites for four climate scenarios corresponding to mean temperature changes of 1.28 to 3.78°C. Freshwater macroinvertebrate toxicity data were collated for 19 chemicals and the hierarchical species sensitivity distribution model was used to predict the sensitivity of untested taxa using relatedness within a Bayesian approach. All four future climate scenarios shifted assemblage compositions, increasing the prevalence of Mollusca, Crustacea and Oligochaeta species, and the insect taxa of Odonata, Chironomidae, and Baetidae species. Contrastingly, decreases were projected for Plecoptera, Ephemeroptera (except for Baetidae) and Coleoptera species. Shifts in taxonomic composition were associated with changes in the percentage of species at risk from chemical exposure. For the 3.78°C climate scenario, 76% of all assemblages became more sensitive to chemicals and for 18 of the 19 chemicals, the percentage of species at risk increased. Climate warming-induced increases in sensitivity were greatest for assemblages exposed to metals and were dependent on baseline assemblage composition, which varied spatially. Climate warming is predicted to result in changes in the use, environmental exposure and toxicity of chemicals. Here we show that, even in the absence of these climate-chemical interactions, shifts in species composition due to climate warming will increase chemical risk and that the impact of chemical pollution on freshwater macroinvertebrate biodiversity may double or quadruple by the end of the 21st century.

Towards Precision Ecotoxicology: Leveraging Evolutionary Conservation of Pharmaceutical and Personal Care Product Targets to Understand Adverse Outcomes Across Species and Life Stages

Translation of environmental science to the practice aims to protect biodiversity and ecosystem services, and our future ability to do so relies on the development of a precision ecotoxicology approach wherein we leverage the genetics and informatics of species to better understand and manage the risks of global pollution. A little over a decade ago, a workshop focusing on the risks of pharmaceuticals and personal care products (PPCPs) in the environment identified a priority research question, "What can be learned about the evolutionary conservation of PPCP targets across species and life stages in the context of potential adverse outcomes and effects?" We review the activities in this area over the past decade, consider prospects of more recent developments, and identify future research needs to develop next-generation approaches for PPCPs and other global chemicals and waste challenges. Environ Toxicol Chem 2024;43:526-536. © 2023 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Spatial variation in the recovery potential of freshwater macroinvertebrate assemblages: Moving towards spatially defined assemblage vulnerability to chemicals

The vulnerability of freshwater biodiversity to chemical stressors is dependent on its ability to resist chemical stress and recover from any stress-induced effects. Spatial variation in recovery has the potential to exacerbate or mitigate assemblage vulnerability but this has not been explored in detail. By combining information on assemblage-specific recovery potential with information on assemblage-specific chemical sensitivity, we have demonstrated that the vulnerability of 3307 macroinvertebrate assemblages to 18 different chemicals is spatially dependent and that recovery potential may reduce chemical risk. The recovery potential of each assemblage was quantified based on trait information and landscape factors using a weighted sum method, but it did not consider succession processes. Recovery potential varied by river type with assemblages in mid-altitude siliceous rivers with small catchments in the west of England having the lowest recovery potential. For 17 or the 18 chemicals investigated, there was a positive correlation between the recovery potential and sensitivity and this was strongest for assemblages exposed to metals. More sensitive assemblages had a higher recovery potential and were therefore potentially less vulnerable than would be expected based on sensitivity alone. Assemblages in rivers with small catchments were the most vulnerable to chemical exposure. Furthermore, assemblages with high vulnerability to insecticide exposure were more prevalent in mid-altitude rivers with siliceous geology in the west of England, whereas assemblages with high vulnerability to metals were more prevalent in lowland rivers with calcareous or mixed geology in the midlands. This study: (i) highlights the importance of spatial context in determining the risk of chemical pollution to freshwater biodiversity; (ii) demonstrates how spatial variation in taxonomic composition influences both the internal and external recovery of assemblages and how landscape factors modify trait-based recovery capabilities; (iii) provides the foundations for spatially-defined vulnerability assessment by identifying ecological scenarios for assessing chemical risk.

Spatial variation in the sensitivity of freshwater macroinvertebrate assemblages to chemical stressors

Assessing spatial variation in the chemical sensitivity of natural assemblages will enhance ecological relevance and reduce uncertainty in ecological risk assessments and the derivation of environmental quality standards (EQSs). However, the majority of species in natural communities have not undergone toxicity testing for any chemical, which poses a major challenge when assessing their sensitivity. We investigated spatial variation and patterns in the sensitivity of 4084 freshwater macroinvertebrate assemblages across England to 5 general-acting chemicals (heavy metals) and 13 specifically acting chemicals (insecticides) using a novel hierarchical species sensitivity distribution method based on taxonomic relatedness. Furthermore, we explored how river typology relates to spatial variation in assemblage sensitivity to chemicals and the potential impacts of such variation on current EQSs. Our findings revealed that, whereas assemblages with similar taxonomic compositions exhibit comparable sensitivity distributions, assemblages with different taxonomic compositions could have very similar or very different sensitivity distributions. The variation in assemblage sensitivity was greater for specifically acting chemicals than for general-acting chemicals and exhibited spatial clustering patterns. These spatial clustering patterns varied depending on the chemical, and the regions where assemblages were most sensitive to metals were generally not the same as the regions where assemblages were most sensitive to insecticides. Spatial variation in assemblage sensitivity was related to river typology with sensitive assemblages being more common than expected in lowland calcareous (or mixed geology) rivers within very small to small catchments. Comparing spatial variation in assemblage-specific chemical sensitivity to EQSs, we found that the operational EQSs in England would protect most study assemblages (i.e., > 99.5 %), although a small proportion of assemblages may face potential risks associated with azinphos-methyl, copper, and malathion. In many cases the EQSs were very precautionary, potentially requiring expensive control measures or restricting beneficial chemical use with no additional environmental benefit. The development of spatially defined EQSs, possibly based on river types, could be developed to target areas that require the highest level of protection and thus strike a balance between the benefits of chemical use and environmental protection.

Uncertainty analysis and updated user guidance for interspecies correlation estimation models and low toxicity compounds

Interspecies correlation estimation (ICE) models are log-linear relationships of acute sensitivity between two species that estimate the sensitivity of an untested species from the known sensitivity of a surrogate. As ICE model use increases globally, additional user guidance is required to ensure consistent use across chemicals and applications. The present study expands ICE uncertainty analyses and user guidance with a focus on low toxicity compounds whose acute values (i.e., reported as mg/L) can be greater than those used to develop a model. In these cases, surrogate values may be outside the ICE model domain and require additional extrapolations to predict acute toxicity. We use the extensive, standardized acute toxicity database underlying ICE models to broadly summarize inter-test variability of acute toxicity data as a measure by which model prediction accuracy can be evaluated. Using the data and models found on the USEPA Web-ICE (www3.epa.gov/webice), we created a set of "truncated" models from data corresponding to the lower 75th percentile of surrogate toxicity. We predicted toxicity for chemicals in the upper 25th percentile as both μg/L beyond the model domain and converted to mg/L (i.e., "scaled" value) and compared these predictions with those from cross-validation of whole ICE models and to the measured value. For ICE models with slopes in the range 0.66-1.33, prediction accuracy of scaled values did not differ from the accuracy of the models when data were entered as μg/L within or beyond the model domain. An uncertainty analysis of ICE confidence intervals was conducted and an interval range of two orders of magnitude was determined to minimize type I and II errors when accepting or rejecting ICE predictions. We updated the ICE user guidance based on these analyses to advance the state of the science for ICE model application and interpretation. Integr Environ Assess Manag 2024;00:1-12. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Hepatic Transcriptomic Responses to Ethinylestradiol in Embryonic Japanese Quail and Double-Crested Cormorant

Understanding species differences in sensitivity to toxicants is a critical issue in ecotoxicology. We recently established that double-crested cormorant (DCCO) embryos are more sensitive than Japanese quail (JQ) to the developmental effects of ethinylestradiol (EE2). We explored how this difference in sensitivity between species is reflected at a transcriptomic level. The EE2 was dissolved in dimethyl sulfoxide and injected into the air cell of eggs prior to incubation at nominal concentrations of 0, 3.33, and 33.3 µg/g egg weight. At midincubation (JQ 9 days; DCCO 16 days), livers were collected from five embryos/treatment group for RNA sequencing. Data were processed and analyzed using EcoOmicsAnalyst and ExpressAnalyst. The EE2 exposure dysregulated 238 and 1,987 genes in JQ and DCCO, respectively, with 78 genes in common between the two species. These included classic biomarkers of estrogen exposure such as vitellogenin and apovitellenin. We also report DCCO-specific dysregulation of Phase I/II enzyme-coding genes and species-specific transcriptional ontogeny of vitellogenin-2. Twelve Kyoto Encyclopedia of Genes and Genomes pathways and two EcoToxModules were dysregulated in common in both species including the peroxisome proliferator-activated receptor (PPAR) signaling pathway and fatty acid metabolism. Similar to previously reported differences at the organismal level, DCCO were more responsive to EE2 exposure than JQ at the gene expression level. Our description of differences in transcriptional responses to EE2 in early life stage birds may contribute to a better understanding of the molecular basis for species differences. Environ Toxicol Chem 2024;00:1-12. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Physiological variables in machine learning QSARs allow for both cross-chemical and cross-species predictions

A major challenge in ecological risk assessment is estimating chemical-induced effects across taxa without species-specific testing. Where ecotoxicological data may be more challenging to gather, information on species physiology is more available for a broad range of taxa. Physiology is known to drive species sensitivity but understanding about the relative contribution of specific underlying processes is still elusive. Consequently, there remains a need to understand which physiological processes lead to differences in species sensitivity. The objective of our study was to utilize existing knowledge about organismal physiology to both understand and predict differences in species sensitivity. Machine learning models were trained to predict chemical- and species-specific endpoints as a function of both chemical fingerprints/descriptors and physiological properties represented by dynamic energy budget (DEB) parameters. We found that random forest models were able to predict chemical- and species-specific endpoints, and that DEB parameters were relatively important in the models, particularly for invertebrates. Our approach illuminates how physiological properties may drive species sensitivity, which will allow more realistic predictions of effects across species without the need for additional animal testing.

Environmental risks of a commonly used pyrethroid: Insights from temporary pond species of the Lake Manyara Basin, Tanzania

Environmental risks posed by widespread pesticide application have attracted global attention. Currently, chemical risk assessments in aquatic environments rely on extrapolation of toxicity data from classic model species. However, similar assessments based on local species could be complementary, particularly for unusual living environments such as temporary ponds. Here, we carried out an environmental risk assessment (ERA) of a pyrethroid model compound, cypermethrin, based on local temporary pond species. First, we measured cypermethrin residue concentrations in rivers, irrigation canals and temporary ponds in the Lake Manyara Basin (LMB). Then, we estimated the environmental risks of cypermethrin by combining these data with acute toxicity data of three resident species across three trophic levels: primary producers (Arthrospira platensis), invertebrate grazers (Streptocephalus lamellifer) and fish (Nothobranchius neumanni). Furthermore, we compared the derived ERA to that obtained using toxicity data from literature of classic model species. Cypermethrin residue concentrations in contaminated systems of the LMB ranged from 0.01 to 57.9 ng/L. For temporary pond species, S. lamellifer was the most sensitive one with a 96 h-LC₅₀ of 0.14 ng/L. Regardless of the assumed exposure concentration (0.01 and 57.9 ng/L), the estimated risks were low for primary producers and high for invertebrate grazers, both for local species as well as for classic model species. The highest detected cypermethrin concentration resulted in a moderate risk estimation for local fish species, while the estimated risk was high when considering classic fish models. Our results confirm that, at least for pyrethroids, ERAs with classic model species are useful to estimate chemical risks in temporary pond ecosystems, and suggest that complementary ERAs based on local species could help to fine-tune environmental regulations to specific local conditions and conservation targets.

From Protein Sequence to Structure: The Next Frontier in Cross-Species Extrapolation for Chemical Safety Evaluations

Computational screening for potentially bioactive molecules using advanced molecular modeling approaches including molecular docking and molecular dynamic simulation is mainstream in certain fields like drug discovery. Significant advances in computationally predicting protein structures from sequence information have also expanded the availability of structures for nonmodel species. Therefore, the objective of the present study was to develop an analysis pipeline to harness the power of these bioinformatics approaches for cross-species extrapolation for evaluating chemical safety. The Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool compares protein-sequence similarity across species for conservation of known chemical targets, providing an initial line of evidence for extrapolation of toxicity knowledge. However, with the development of structural models from tools like the Iterative Threading ASSEmbly Refinement (ITASSER), analyses of protein structural conservation can be included to add further lines of evidence and generate protein models across species. Models generated through such a pipeline could then be used for advanced molecular modeling approaches in the context of species extrapolation. Two case examples illustrating this pipeline from SeqAPASS sequences to I-TASSER-generated protein structures were created for human liver fatty acid-binding protein (LFABP) and androgen receptor (AR). Ninety-nine LFABP and 268 AR protein models representing diverse species were generated and analyzed for conservation using template modeling (TM)-align. The results from the structural comparisons were in line with the sequence-based SeqAPASS workflow, adding further evidence of LFABL and AR conservation across vertebrate species. The present study lays the foundation for expanding the capabilities of the web-based SeqAPASS tool to include structural comparisons for species extrapolation, facilitating more rapid and efficient toxicological assessments among species with limited or no existing toxicity data. Environ Toxicol Chem 2023;42:463-474. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Joint survival modelling for multiple species exposed to toxicants

In environmental risk assessment (ERA), the multitude of compounds and taxa demands cross-species extrapolation to cover the variability in sensitivity to toxicants. However, only the impact of a single compound to a single species is addressed by the general unified threshold model of survival (GUTS). The reduced GUTS is the recommended model to analyse lethal toxic effects in regulatory aquatic ERA. GUTS considers toxicokinetics and toxicodynamics. Two toxicodynamic approaches are considered: Stochastic death (SD) assumes that survival decreases with an increasing internalized amount of the toxicant. Individual tolerance (IT) assumes that individuals vary in their tolerance to toxic exposure. Existing theory suggests that the product of the threshold z_w and killing rate b_w (both SD toxicodynamic parameters) are constant across species or compounds if receptors and target sites are shared. We extend that theory and show that the shape parameter β of the loglogistic threshold distribution in IT is also constant. To verify the predicted relationships, we conducted three tests using toxicity studies for eight arthropods exposed to the insecticide flupyradifurone. We confirmed previous verifications of the relation- between SD parameters, and the newly established relation for the IT parameter β. We enhanced GUTS to jointly model survival for multiple species with shared receptors and pathways by incorporating the relations among toxicodynamic parameters described above. The joint GUTS exploits the shared parameter relations and therefore constrains parameter uncertainty for each of the separate species. Particularly for IT, the joint GUTS more precisely predicted risk to the separate species than the standard single species GUTS under environmentally realistic exposure. We suggest that joint GUTS modelling can improve cross-species extrapolation in regulatory ERA by increasing the reliability of risk estimates and reducing animal testing. Furthermore, the shared toxicodynamic response provides potential to reduce complexity of ecosystem models.

Three perspectives on the prediction of chemical effects in ecosystems

The increasing production, use and emission of synthetic chemicals into the environment represents a major driver of global change. The large number of synthetic chemicals, limited knowledge on exposure patterns and effects in organisms and their interaction with other global change drivers hamper the prediction of effects in ecosystems. However, recent advances in biomolecular and computational methods are promising to improve our capacity for prediction. We delineate three idealised perspectives for the prediction of chemical effects: the suborganismal, organismal and ecological perspective, which are currently largely separated. Each of the outlined perspectives includes essential and complementary theories and tools for prediction but captures only part of the phenomenon of chemical effects. Links between the perspectives may foster predictive modelling of chemical effects in ecosystems and extrapolation between species. A major challenge for the linkage is the lack of data sets simultaneously covering different levels of biological organisation (here referred to as biological levels) as well as varying temporal and spatial scales. Synthesising the three perspectives, some central aspects and associated types of data seem particularly necessary to improve prediction. First, suborganism- and organism-level responses to chemicals need to be recorded and tested for relationships with chemical groups and organism traits. Second, metrics that are measurable at many biological levels, such as energy, need to be scrutinised for their potential to integrate across levels. Third, experimental data on the simultaneous response over multiple biological levels and spatiotemporal scales are required. These could be collected in nested and interconnected micro- and mesocosm experiments. Lastly, prioritisation of processes involved in the prediction framework needs to find a balance between simplification and capturing the essential complexity of a system. For example, in some cases, eco-evolutionary dynamics and interactions may need stronger consideration. Prediction needs to move from a static to a real-world eco-evolutionary view.

Ecotoxicological Assessment of "Glitter" Leachates in Aquatic Ecosystems: An Integrated Approach

The most worrisome fraction within plastic pollution is that of microplastics (MP). A category of MP almost completely ignored is that of glitter. The objective of this study is to test the toxicity of nine types of glitter leachate (3 soak times: 3, 90 and 180 days) on model organisms in freshwater (Allivibrio fischeri, Raphidocelis subcapitata, Daphnia magna) and saltwater (Allivibrio fischeri, Phaeodactylum tricornutum, Paracentrotus lividus). An integrated approach was applied to obtain the percentage of ecotoxicological risk. The results show that (i) photosynthesizing primary producers are the most sensitive trophic level; (ii) algae transitioned from growth inhibition to biostimulation; (iii) D. magna showed higher sensitivity after 48 h compared to 24 h; (iv) A. fischeri responded more strongly in saltwater than in freshwater. The integrated data show a greater risk associated with the marine environment, with the highest risk for glitters that are hexagonal and composed of poly-methyl-methacrylate. Our multivariate analysis shows that the toxicity of plastic leaching is a complex phenomenon that depends on the sensitivity of the species, in some cases on the soaking time and on the medium, and is not clearly linked to the polymer type, the contact area or the colors of the particles.

Heterogeneity in biological assemblages and exposure in chemical risk assessment: Exploring capabilities and challenges in methodology with two landscape-scale case studies

Chemical exposure concentrations and the composition of ecological receptors (e.g., species) vary in space and time, resulting in landscape-scale (e.g. catchment) heterogeneity. Current regulatory, prospective chemical risk assessment frameworks do not directly address this heterogeneity because they assume that reasonably worst-case chemical exposure concentrations co-occur (spatially and temporally) with biological species that are the most sensitive to the chemical's toxicity. Whilst current approaches may parameterise fate models with site-specific data and aim to be protective, a more precise understanding of when and where chemical exposure and species sensitivity co-occur enables risk assessments to be better tailored and applied mitigation more efficient. We use two aquatic case studies covering different spatial and temporal resolution to explore how geo-referenced data and spatial tools might be used to account for landscape heterogeneity of chemical exposure and ecological assemblages in prospective risk assessment. Each case study followed a stepwise approach: i) estimate and establish spatial chemical exposure distributions using local environmental information and environmental fate models; ii) derive toxicity thresholds for different taxonomic groups and determine geo-referenced distributions of exposure-toxicity ratios (i.e., potential risk); iii) overlay risk data with the ecological status of biomonitoring sites to determine if relationships exist. We focus on demonstrating whether the integration of relevant data and potential approaches is feasible rather than making comprehensive and refined risk assessments of specific chemicals. The case studies indicate that geo-referenced predicted environmental concentration estimations can be achieved with available data, models and tools but establishing the distribution of species assemblages is reliant on the availability of a few sources of biomonitoring data and tools. Linking large sets of geo-referenced exposure and biomonitoring data is feasible but assessment of risk will often be limited by the availability of ecotoxicity data. The studies highlight the important influence that choices for aggregating data and for the selection of statistical metrics have on assessing and interpreting risk at different spatial scales and patterns of distribution within the landscape. Finally, we discuss approaches and development needs that could help to address environmental heterogeneity in chemical risk assessment.

Are resource recovery insects safe for feed and food? A screening approach for bioaccumulative trace organic contaminants

Most bioaccumulation assessments select one or several compound classes a priori for analysis performed by either liquid or gas chromatography coupled with mass spectrometry (LC-MS or GC-MS). When organisms are exposed to complex mixtures of trace organic contaminants (TOrCs), targeted chemical assays limit understanding of contaminant profiles in biological tissues and associated risks. We used a semi-quantitative suspect-screening approach to assess the bioaccumulation potential of diverse TOrCs in black soldier fly larvae (BSFL) using almond hulls (by-products of the booming almond industry in California) as test substrates. BSFL digestion is gaining traction as a resource recovery strategy to generate animal feed from low-value organic wastes. We screened almond hulls from six California farms for the presence of 5728 TOrCs using high resolution mass spectrometry. We then categorized the risk potential of 46 TOrCs detected in the hulls based on their predicted bioaccumulation, persistence, and toxicity in order to select two hulls for an in situ BSFL bioaccumulation screening study. We analyzed larvae tissues and feeding substrate initially and after 14 days of growth using targeted, suspect-screening, and nontarget-screening methods. The survival rate of BSFL in all rearing reactors was greater than 90%, indicating low toxicity of the substrates to BSFL. Esfenvalerate, cyhalothrin, and bifenthrin were the most abundant pyrethroids quantified (81.7 to 381.6 ng/g-dw) in the hulls. Bifenthrin bioaccumulated in BSFL tissues (14-day bioaccumulation factor, BAF, of 2.17 ± 0.24). For nontarget analysis, kendrick mass defect (KMD) analysis of PFAS homologous series revealed hydrogen-substituted perfluoroalkyl carboxylic acids (H-PFCAs) in the hulls and BSFL tissues after growth. Our approach demonstrates the utility of suspect-screening in chemical safety assessments when organic wastes with highly diverse and variable contaminant profiles are used in resource recovery pipelines.

Pesticides alter ecosystem respiration via phytoplankton abundance and community structure: Effects on the carbon cycle?

Freshwater systems are critical to life on earth, yet they are threatened by the increasing rate of synthetic chemical pollution. Current predictions of the effects of synthetic chemicals on freshwater ecosystems are hampered by the sheer number of chemical contaminants entering aquatic systems, the diversity of organisms inhabiting these systems, the myriad possible direct and indirect effects resulting from these combinations, and uncertainties concerning how contaminants might alter ecosystem metabolism via changes in biodiversity. To address these knowledge gaps, we conducted a mesocosm experiment that elucidated the responses of ponds composed of phytoplankton and zooplankton to standardized concentrations of 12 pesticides, nested within four pesticide classes, and two pesticide types. We show that the effects of the pesticides on algae were consistent within herbicides and insecticides and that responses of over 70 phytoplankton species and genera were consistent within broad taxonomic groups. Insecticides generated top-down effects on phytoplankton community composition and abundance, which were associated with persistent increases in ecosystem respiration. Insecticides had direct toxic effects on cladocerans, which led to competitive release of copepods. These changes in the zooplankton community led to a decrease in green algae and a modest increase in diatoms. Herbicides did not change phytoplankton composition but reduced total phytoplankton abundance. This reduction in phytoplankton led to short-term decreases in ecosystem respiration. Given that ponds release atmospheric carbon and that worldwide pesticide pollution continues to increase exponentially, scientists and policy makers should pay more attention to the ways pesticides alter the carbon cycle in ponds via changes in communities, as demonstrated by our results. Our results show that these predictions can be simplified by grouping pesticides into types and species into functional groups. Adopting this approach provides an opportunity to improve the efficiency of risk assessment and mitigation responses to global change.

Using Transcriptomics and Metabolomics to Understand Species Differences in Sensitivity to Chlorpyrifos in Japanese Quail and Double-Crested Cormorant Embryos

Modern 21st-century toxicity testing makes use of omics technologies to address critical questions in toxicology and chemical management. Of interest are questions relating to chemical mechanisms of toxicity, differences in species sensitivity, and translation of molecular effects to observable apical endpoints. Our study addressed these questions by comparing apical outcomes and multiple omics responses in early-life stage exposure studies with Japanese quail (Coturnix japonica) and double-crested cormorant (Phalacrocorax auritus), representing a model and ecological species, respectively. Specifically, we investigated the dose-dependent response of apical outcomes as well as transcriptomics and metabolomics in the liver of each species exposed to chlorpyrifos, a widely used organophosphate pesticide. Our results revealed a clear pattern of dose-dependent disruption of gene expression and metabolic profiles in Japanese quail but not double-crested cormorant at similar chlorpyrifos exposure concentrations. The difference in sensitivity between species was likely due to higher metabolic transformation of chlorpyrifos in Japanese quail compared to double-crested cormorant. The most impacted biological pathways after chlorpyrifos exposure in Japanese quail included hepatic metabolism, oxidative stress, endocrine disruption (steroid and nonsteroid hormones), and metabolic disease (lipid and fatty acid metabolism). Importantly, we show consistent responses across biological scales, suggesting that significant disruption at the level of gene expression and metabolite profiles leads to observable apical responses at the organism level. Our study demonstrates the utility of evaluating effects at multiple biological levels of organization to understand how modern toxicity testing relates to outcomes of regulatory relevance, while also highlighting important, yet poorly understood, species differences in sensitivity to chemical exposure. Environ Toxicol Chem 2021;40:3019-3033. © 2021 SETAC.

Assessing chemical risk within an ecosystem services framework: Implementation and added value

An ecosystem services (ES) approach to chemical risk assessment has many potential advantages, but there are also substantial challenges regarding its implementation. We report the findings of a multi-stakeholder workshop that evaluated the feasibility of adopting an ES approach to chemical risk assessment using currently available tools and data. Also evaluated is the added value such an approach would bring to environmental decision making. The aim was to build consensus across disparate stakeholders and to co-produce a common understanding of the regulatory benefits and feasibility of implementing an ES approach in European chemicals regulation. Workshop discussions were informed by proof of concept studies and resulted in the development of a novel tiered framework for assessing chemical risk to ES delivery. There was consensus on the substantial added value of adopting an ES-based approach for regulatory decision making. Ecosystem services provide a common currency and a 'unifying approach' across environmental compartments, stressors and regulatory frameworks. The ES approach informs prioritisation of risk and remedial action and aids risk communication and risk management. It facilitates a more holistic assessment, enables ES trade-offs to be compared across alternative interventions, and supports comparative risk assessments and a socio-economic analysis of management options and decisions. Key to realising this added value is a shift away from using a single threshold value to categorise risk, towards a consideration of the exposure-effect distribution for individual ES of interest. Also required is the development of an integrated systems-level approach across regulatory frameworks and agreement on specific protection goals and scenarios for framing environmental risk assessments. The need to further develop tools for extrapolating toxicity data to service providers and ES delivery, including logic chains and ecological production functions, was highlighted. Also agreed was the need for methods and metrics for ES valuation to be used in assessing trade-offs.

A Method to Screen for Consistency of Effect in Laboratory Toxicity Tests: A Case Study with Anurans and the Herbicide Atrazine

This paper presents a semiquantitative method to help ecotoxicologists evaluate the consistency of data within the available peer-reviewed literature. In this case study, we queried whether there is consistent evidence of direct toxicity in Anurans exposed to atrazine at concentrations ≤ 100 μg/L under laboratory conditions. Atrazine was selected because of the relatively large repository of Anuran toxicity data. To accomplish this, we interrogated available data found in recent quantitative weight-of-evidence risk assessments for atrazine with a series of yes or no questions developed a priori. The questions examined consistency of reported effects within and between studies, within and between species, and across a wide range of endpoints categories (e.g., survivorship, growth and development, reproduction). The analysis found no compelling evidence of a consistent direct effect in Anurans around growth and development, reproduction, or survivorship at concentrations of up to at least 100 μg/L atrazine in laboratory studies. Further work is needed to refine the approach, including accounting for the magnitude of the reported effects. However, we recommend that ecotoxicologists employ some method of formal consistency of effects assessment method routinely before performing toxicity tests, in the contextualizing of new data, and in reviews of contaminants.

Mechanistic Effect Modeling of Earthworms in the Context of Pesticide Risk Assessment: Synthesis of the FORESEE Workshop

Earthworms are important ecosystem engineers, and assessment of the risk of plant protection products toward them is part of the European environmental risk assessment (ERA). In the current ERA scheme, exposure and effects are represented simplistically and are not well integrated, resulting in uncertainty when the results are applied to ecosystems. Modeling offers a powerful tool to integrate the effects observed in lower tier laboratory studies with the environmental conditions under which exposure is expected in the field. This paper provides a summary of the (In)Field Organism Risk modEling by coupling Soil Exposure and Effect (FORESEE) Workshop held 28-30 January 2020 in Düsseldorf, Germany. This workshop focused on toxicokinetic-toxicodynamic (TKTD) and population modeling of earthworms in the context of ERA. The goal was to bring together scientists from different stakeholder groups to discuss the current state of soil invertebrate modeling and to explore how earthworm modeling could be applied to risk assessments, in particular how the different model outputs can be used in the tiered ERA approach. In support of these goals, the workshop aimed at addressing the requirements and concerns of the different stakeholder groups to support further model development. The modeling approach included 4 submodules to cover the most relevant processes for earthworm risk assessment: environment, behavior (feeding, vertical movement), TKTD, and population. Four workgroups examined different aspects of the model with relevance for risk assessment, earthworm ecology, uptake routes, and cross-species extrapolation and model testing. Here, we present the perspectives of each workgroup and highlight how the collaborative effort of participants from multidisciplinary backgrounds helped to establish common ground. In addition, we provide a list of recommendations for how earthworm TKTD modeling could address some of the uncertainties in current risk assessments for plant protection products. Integr Environ Assess Manag 2021;17:352-363. © 2020 SETAC.

Pesticide effects on fish cholinesterase variability and mean activity: A meta-analytic review

Fish cholinesterases (ChEs) - like acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) - are common biomarkers of environmental contamination due to their sensitivity to a variety of toxicants. To understand pesticide effects on fish ChEs mean activity and variability, we conducted a systematic review and meta-analyses. Our goal was to verify (i) if brain and muscle ChEs responded differently to pesticide exposure; (ii) how fish size and life stage (i.e., juvenile and adult) influence ChEs variability and mean activity; (iii) what type of pesticides (i.e., herbicide, insecticide, and fungicide) has the strongest effect, and if the analytical-grade compounds differ from commercial formulations; (iv) if increasing concentrations combined with prolonged exposure leads to stronger ChEs inhibition; and (v) how each class of pesticide affects these enzymes. We validated ChEs reliability as biomarkers and identified factors influencing their response. Regardless of tissue, BChE response was more variable than AChE, and no difference between their average activity was detected. The size of juvenile fish is an important factor affecting ChEs mean activity and variability, whereas pesticide had no significant effect on adult fish ChEs. Insecticides were stronger inhibitors compared to herbicides and fungicides. Analytical-grade compounds decreased ChEs mean activity to a higher degree than commercial formulations. The combined effect of concentration and time was only significant for fungicides and insecticides. Among classes, organophosphorus insecticides had the strongest effect on ChEs, followed by carbamates, organochlorines, and pyrethroids. Organophosphorus herbicides and oxazolidinones were the only herbicides to decrease ChEs mean activity significantly, and their effects were similar from those of pyrethroids and organochlorines. Additionally, our results identified research gaps, such as the small number of studies on fungicides, neonicotinoids and other relevant pesticides. These findings suggest future directions, which might help researchers identify robust cause-effect relationships between fish ChEs and pesticides.