Modeling the Sensitivity of Aquatic Macroinvertebrates to Chemicals Using Traits

Building species trait-specific nano-QSARs: Model stacking, navigating model uncertainties and limitations, and the effect of dataset size

A strong need exists for broadly applicable nano-QSARs, capable of predicting toxicological outcomes towards untested species and nanomaterials, under different environmental conditions. Existing nano-QSARs are generally limited to only a few species but the inclusion of species characteristics into models can aid in making them applicable to multiple species, even when toxicity data is not available for biological species. Species traits were used to create classification- and regression machine learning models to predict acute toxicity towards aquatic species for metallic nanomaterials. Afterwards, the individual classification- and regression models were stacked into a meta-model to improve performance. Additionally, the uncertainty and limitations of the models were assessed in detail (beyond the OECD principles) and it was investigated whether models would benefit from the addition of more data. Results showed a significant improvement in model performance following model stacking. Investigation of model uncertainties and limitations highlighted the discrepancy between the applicability domain and accuracy of predictions. Data points outside of the assessed chemical space did not have higher likelihoods of generating inadequate predictions or vice versa. It is therefore concluded that the applicability domain does not give complete insight into the uncertainty of predictions and instead the generation of prediction intervals can help in this regard. Furthermore, results indicated that an increase of the dataset size did not improve model performance. This implies that larger dataset sizes may not necessarily improve model performance while in turn also meaning that large datasets are not necessarily required for prediction of acute toxicity with nano-QSARs.

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.

Assessment of the impacts of GABA and AChE targeting pesticides on freshwater invertebrate family richness in English Rivers

Globally, riverine system biodiversity is threatened by a range of stressors, spanning pollution, sedimentation, alterations to water flow, and climate change. Pesticides have been associated with population level impacts on freshwater invertebrates for acute high-level exposures, but far less is known about the chronic impact of episodic exposure to specific classes of pesticides or their mixtures. Here we employed the use of the UK Environment Agency's monitoring datasets over 40 years (covering years 1980 to 2019) to assess the impacts of AChE (acetylcholinesterase) and GABA (gamma-aminobutyric acid) receptor targeting pesticides on invertebrate family richness at English river sites. Concentrations of AChE and GABA pesticides toxic to freshwater invertebrates occurred (measured) across 18 of the 66 river sites assessed. For one of the three river sites (all found in the Midlands region of England) where data recorded over the past 40 years were sufficient for robust modelling studies, both AChE and GABA pesticides associated with invertebrate family richness. Here, where AChE total pesticide concentrations were classified as high, 46 of 64 invertebrate families were absent, and where GABA total pesticide concentration were classified as high, 16 of 64 invertebrate families were absent. Using a combination of field evidence and laboratory toxicity thresholds for population relevant endpoints we identify families of invertebrates most at risk in the selected English rivers to AChE and GABA pesticides. We, furthermore, provide strong evidence that the absence of the invertebrate family Polycentropodidae (caddisfly) from one field site is due to exposure effects to AChE pesticides.

A broad-taxa approach as an important concept in ecotoxicological studies and pollution monitoring

Aquatic invertebrates play a pivotal role in (eco)toxicological assessments because they offer ethical, cost-effective and repeatable testing options. Additionally, their significance in the food chain and their ability to represent diverse aquatic ecosystems make them valuable subjects for (eco)toxicological studies. To ensure consistency and comparability across studies, international (eco)toxicology guidelines have been used to establish standardised methods and protocols for data collection, analysis and interpretation. However, the current standardised protocols primarily focus on a limited number of aquatic invertebrate species, mainly from Arthropoda, Mollusca and Annelida. These protocols are suitable for basic toxicity screening, effectively assessing the immediate and severe effects of toxic substances on organisms. For more comprehensive and ecologically relevant assessments, particularly those addressing long-term effects and ecosystem-wide impacts, we recommended the use of a broader diversity of species, since the present choice of taxa exacerbates the limited scope of basic ecotoxicological studies. This review provides a comprehensive overview of (eco)toxicological studies, focusing on major aquatic invertebrate taxa and how they are used to assess the impact of chemicals in diverse aquatic environments. The present work supports the use of a broad-taxa approach in basic environmental assessments, as it better represents the natural populations inhabiting various ecosystems. Advances in omics and other biochemical and computational techniques make the broad-taxa approach more feasible, enabling mechanistic studies on non-model organisms. By combining these approaches with in vitro techniques together with the broad-taxa approach, researchers can gain insights into less-explored impacts of pollution, such as changes in population diversity, the development of tolerance and transgenerational inheritance of pollution responses, the impact on organism phenotypic plasticity, biological invasion outcomes, social behaviour changes, metabolome changes, regeneration phenomena, disease susceptibility and tissue pathologies. This review also emphasises the need for harmonised data-reporting standards and minimum annotation checklists to ensure that research results are findable, accessible, interoperable and reusable (FAIR), maximising the use and reusability of data. The ultimate goal is to encourage integrated and holistic problem-focused collaboration between diverse scientific disciplines, international standardisation organisations and decision-making bodies, with a focus on transdisciplinary knowledge co-production for the One-Health approach.

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.

Proposing Effective Ecotoxicity Test Species for Chemical Safety Assessment in East Asia: A Review

East Asia leads the global chemical industry, but environmental chemical risk in these countries is an emerging concern. Despite this, only a few native species that are representative of East Asian environments are listed as test species in international guidelines compared with those native to Europe and America. This review suggests that Zacco platypus, Misgurnus anguillicaudatus, Hydrilla verticillata, Neocaridina denticulata spp., and Scenedesmus obliquus, all resident to East Asia, are promising test species for ecotoxicity tests. The utility of these five species in environmental risk assessment (ERA) varies depending on their individual traits and the state of ecotoxicity research, indicating a need for different applications of each species according to ERA objectives. Furthermore, the traits of these five species can complement each other when assessing chemical effects under diverse exposure scenarios, suggesting they can form a versatile battery for ERA. This review also analyzes recent trends in ecotoxicity studies and proposes emerging research issues, such as the application of alternative test methods, comparative studies using model species, the identification of specific markers for test species, and performance of toxicity tests under environmentally relevant conditions. The information provided on the utility of the five species and alternative issues in toxicity tests could assist in selecting test species suited to study objectives for more effective ERA.

Integrative ecotoxicity evaluation of Cd, Cu, Zn and Ni in aquatic animals reveals high tolerance of Artemia franciscana

Heavy metals pose a significant threat to animals in aquatic environments due to the adverse effects they exert. Species of the genus Artemia have been described as heavy metal tolerant, but the sensitivity/tolerance range for these species has not been established. In the present study, the toxicity of Cd, Cu, Zn and Ni as reported in the ECOTOX and Web of Science databases was examined for Artemia franciscana and compared with other species and taxonomic groups using an integrative ecotoxicity evaluation. The hazard concentration for 5% of the species (HC5) of acute toxicity tests (24-96 h), obtained through a species sensitivity distribution (SSD) indicated that Cu (0.02 mg/L) and Cd (0.03 mg/L) were the metals with the highest toxicity to aquatic animals followed by Zn (0.15 mg/L) and Ni (0.23 mg/L). In addition to the higher hazard of Cu and Cd to aquatic animals, the comparison of acute LC50 values for A. franciscana indicates lower toxicity of Cd followed by Cu, Zn, and Ni (200.0, 14.5, 9.5, and 0.6 mg/L, respectively). Using the SSD and physiological sensitivity (S) approaches, it was demonstrated that A. franciscana is relatively tolerant to Cd (SSD= HC99; S = 2.21), Cu (SSD= HC97; S = 2.00), Zn (SSD= HC90; S = 1.29) and Ni (SSD= HC83; S = 0.96) compared with other species and taxonomic groups. It appears that tolerance to the metals Cd, Cu, Zn and Ni is unique to the family Artemiidae within the order Anostraca, as the families Streptocephalidae and Thamnocephalidae are not tolerant (have negative S values). Our study confirmed that as expected, A. franciscana presents higher tolerance to Cd, Cu, Zn, and Ni than other aquatic animals. Our findings confirm that A. franciscana can be used as a model organism to understand mechanisms involved in tolerance to heavy metals, mainly Cd and Cu, which are considered highly toxic to other animals.

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.

Effects of heat and pesticide stress on life history, physiology and the gut microbiome of two congeneric damselflies that differ in stressor tolerance

The combined impact of toxicants and warming on organisms is getting increased attention in ecotoxicology, but is still hard to predict, especially with regard to heat waves. Recent studies suggested that the gut microbiome may provide mechanistic insights into the single and combined stressor effects on their host. We therefore investigated effects of sequential exposure to a heat spike and a pesticide on both the phenotype (life history and physiology) and the gut microbiome composition of damselfly larvae. We compared the fast-paced Ischnura pumilio, which is more tolerant to both stressors, with the slow-paced I. elegans, to obtain mechanistic insights into species-specific stressor effects. The two species differed in gut microbiome composition, potentially contributing to their pace-of-life differences. Intriguingly, there was a general resemblance between the stressor response patterns in the phenotype and in the gut microbiome, whereby both species responded broadly similar to the single and combined stressors. The heat spike negatively affected the life history of both species (increased mortality, reduced growth rate), which could be explained not only by shared negative effects on physiology (inhibition of acetylcholinesterase, increase of malondialdehyde), but also by shared effects on gut bacterial species' abundances. The pesticide only had negative effects (reduced growth rate, reduced net energy budget) in I. elegans. The pesticide generated shifts in the bacterial community composition (e.g. increased abundance of Sphaerotilus and Enterobacteriaceae in the gut microbiome of I. pumilio), which potentially contributed to the relatively higher pesticide tolerance of I. pumilio. Moreover, in line with the response patterns in the host phenotype, the effects of the heat spike and the pesticide on the gut microbiome were mainly additive. By contrasting two species differing in stress tolerance, our results suggest that response patterns in the gut microbiome may improve our mechanistic understanding of single and combined stressor effects.

Genes-to-Pathways Species Conservation Analysis: Enabling the Exploration of Conservation of Biological Pathways and Processes Across Species

The last two decades have witnessed a strong momentum toward integration of cell-based and computational approaches in safety assessments. This is fueling a global regulatory paradigm shift toward reduction and replacement of the use of animals in toxicity tests while promoting the use of new approach methodologies. The understanding of conservation of molecular targets and pathways provides an opportunity to extrapolate effects across species and ultimately to determine the taxonomic applicability domain of assays and biological effects. Despite the wealth of genome-linked data available, there is a compelling need for improved accessibility, while ensuring that it reflects the underpinning biology. We present the novel pipeline Genes-to-Pathways Species Conservation Analysis (G2P-SCAN) to further support understanding on cross-species extrapolation of biological processes. This R package extracts, synthetizes, and structures the data available from different databases, that is, gene orthologs, protein families, entities, and reactions, linked to human genes and respective pathways across six relevant model species. The use of G2P-SCAN enables the overall analysis of orthology and functional families to substantiate the identification of conservation and susceptibility at the pathway level. In the present study we discuss five case studies, demonstrating the validity of the developed pipeline and its potential use as species extrapolation support. We foresee this pipeline will provide valuable biological insights and create space for the use of mechanistically based data to inform potential species susceptibility for research and safety decision purposes. Environ Toxicol Chem 2023;42:1152-1166. © 2023 UNILEVER GLOBAL IP LTD. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Species sensitivity distributions for two widely used anti-sea lice chemotherapeutants in the salmon aquaculture industry

The main objective of the present study is to construct acute aquatic species sensitivity distributions (SSD) and generate proposed HC₅ values (i.e. the hazardous concentration for which 5 % of species are affected or potentially affected) for two aquaculture anti-sea lice bath pesticides, azamethiphos, and hydrogen peroxide. These values could be used as the basis for the establishment of environmental quality standards (EQS). We have generated SSDs and inferred HC₅ values for mortality and sublethal endpoints using LC50, EC50, and NOEC/LOEC data points separately and for each bath pesticide. Through the examination of literature data on the toxicity of both compounds, we opted to use tests with limited exposure times to ensure environmental relevance for bath pesticides. We also separated life stages for some of the sensitive taxa to account for differences in sensitivities and risk of exposure. The resulting threshold concentrations in environmental seawater are 0.10 μg/L for azamethiphos and 0.15 mg/L for hydrogen peroxide. These suggested azamethiphos and hydrogen peroxide thresholds are comparable to some previously reported EQS values. Further considerations need to be included in how to better use these thresholds in a regulatory context in relation to dispersion patterns. It is also clear that delayed mortality and sublethal effects documented in the literature require further study to fully anticipate the environmental risks of using these two bath pesticides.

Using life-history trait variation to inform ecological risk assessments for threatened and endangered plant species

Developing population models for assessing risks to terrestrial plant species listed as threatened or endangered under the Endangered Species Act (ESA) is challenging given a paucity of data on their life histories. The purpose of this study was to develop a novel approach for identifying relatively data-rich nonlisted species that could serve as representatives for species listed under the ESA in the development of population models to inform risk assessments. We used the USDA PLANTS Database, which provides data on plants present in the US territories, to create a list of herbaceous plants. A total of 8742 species was obtained, of which 344 were listed under the ESA. Using the most up-to-date phylogeny for vascular plants in combination with a database of matrix population models for plants (COMPADRE) and cluster analyses, we investigated how listed species were distributed across the plant phylogeny, grouped listed and nonlisted species according to their life history, and identified the traits distinguishing the clusters. We performed elasticity analyses to determine the relative sensitivity of population growth rate to perturbations of species' survival, growth, and reproduction and compared these across clusters and between listed and nonlisted species. We found that listed species were distributed widely across the plant phylogeny as well as clusters, suggesting that listed species do not share a common evolution or life-history characteristics that would make them uniquely vulnerable. Lifespan and age at maturity were more important for distinguishing clusters than were reproductive traits. For clusters that were intermediate in their lifespan, listed and nonlisted species responded similarly to perturbations of their life histories. However, for clusters at either extreme of lifespan, the response to survival perturbations varied depending on conservation status. These results can be used to guide the choice of representative species for population model development in the context of ecological risk assessment. Integr Environ Assess Manag 2023;19:213-223. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

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.

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.

Contribution of sediment contamination to multi-stress in lowland waters

Water bodies in densely populated lowland areas are often impacted by multiple stressors. At these multi-stressed sites, it remains challenging to quantify the contribution of contaminated sediments. This study, therefore, aimed to elucidate the contribution of sediment contamination in 16 multi-stressed drainage ditches throughout the Netherlands. To this end an adjusted TRIAD framework was applied, where 1) contaminants and other variables in the sediment and the overlying water were measured, 2) whole-sediment laboratory bioassays were performed using larvae of the non-biting midge Chironomus riparius, and 3) the in situ benthic macroinvertebrate community composition was determined. It was hypothesized that the benthic macroinvertebrate community composition would respond to all jointly present stressors in both water and sediment, whereas the whole-sediment bioassays would only respond to the stressors present in the sediment. The benthic macroinvertebrate community composition was indeed related to multiple stressors in both water and sediment. Taxa richness was positively correlated with the presence of PO₄-P in the water, macrophyte cover and some pesticides. Evenness, the number of Trichoptera families and the SPEAR_pesticides were positively correlated to the C:P ratios in the sediment, whilst negative correlations were observed with various contaminants in both the water and sediment. The whole-sediment bioassays with C. riparius positively related to the nutrient content of the sediment, whereas no negative relations to the sediment-associated contaminants were observed, even though the lowered SPEAR_pesticides index indicated contaminant effects in the field. Therefore, it was concluded that sediment contamination was identified as one of the various stressors that potentially drove the benthic macroinvertebrate community composition in the multi-stressed drainage ditches, but that nutrients may have masked the adverse effects caused by low and diverse sediment contaminants on C. riparius in the bioassays.

Functional measures as potential indicators of down-the-drain chemical stress in freshwater ecological risk assessment

Conventional ecological risk assessment (ERA) predominately evaluates the impact of individual chemical stressors on a limited range of taxa, which are assumed to act as proxies to predict impacts on freshwater ecosystem function. However, it is recognized that this approach has limited ecological relevance. We reviewed the published literature to identify measures that are potential functional indicators of down-the-drain chemical stress, as an approach to building more ecological relevance into ERA. We found wide variation in the use of the term "ecosystem function," and concluded it is important to distinguish between measures of processes and measures of the capacity for processes (i.e., species' functional traits). Here, we present a classification of potential functional indicators and suggest that including indicators more directly connected with processes will improve the detection of impacts on ecosystem functioning. The rate of leaf litter breakdown, oxygen production, carbon dioxide consumption, and biomass production have great potential to be used as functional indicators. However, the limited supporting evidence means that further study is needed before these measures can be fully implemented and interpreted within an ERA and regulatory context. Sensitivity to chemical stress is likely to vary among functional indicators depending on the stressor and ecosystem context. Therefore, we recommend that ERA incorporates a variety of indicators relevant to each aspect of the function of interest, such as a direct measure of a process (e.g., rate of leaf litter breakdown) and a capacity for a process (e.g., functional composition of macroinvertebrates), alongside structural indicators (e.g., taxonomic diversity of macroinvertebrates). Overall, we believe that the consideration of functional indicators can add value to ERA by providing greater ecological relevance, particularly in relation to indirect effects, functional compensation (Box 1), interactions of multiple stressors, and the importance of ecosystem context. Environ Assess Manag 2022;18:1135-1147. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Size- and sex-related sensitivity differences of aquatic crustaceans to imidacloprid

Field collected aquatic invertebrates are often used as test organisms in the refinement of the standard Tier 1 risk assessment of various pollutants. This approach can provide insights into the effects of pollutants on the natural environment. However, researchers often pragmatically select test organisms of a specific sex and/or size, which may not represent the sensitivity of the whole population. To investigate such intraspecies sensitivity differences, we performed standard acute toxicity and toxicokinetic tests with different size classes and sex of Gammarus pulex and Asellus aquaticus. Furthermore, toxicokinetics and toxicodynamics models were used to understand the mechanism of the intraspecies sensitivity differences. We used neonates, juveniles and male and female adults in separate dedicated experiments, in which we exposed the animals to imidacloprid and its bioactive metabolite, imidacloprid-olefin. For both species, we found that neonates were the most sensitive group. For G. pulex, the sensitivity decreased linearly with size, which can be explained by the size-related uptake rate constant in the toxicokinetic process and size-related threshold value in the toxicodynamic process. For A. aquaticus, female adults were least sensitive to imidacloprid, which could be explained by a low internal biotransformation of imidacloprid to imidacloprid-olefin. Besides, imidacloprid-olefin was more toxic than imidacloprid to A. aquaticus, with differences being 8.4 times for females and 2.7 times for males. In conclusion, we established size-related sensitivity differences for G. pulex and sex-related sensitivity for A. aquaticus, and intraspecies differences can be explained by both toxicokinetic and toxicodynamic processes. Our findings suggest that to protect populations in the field, we should consider the size and sex of focal organisms and that a pragmatic selection of test organisms of equal size and/or sex can underestimate the sensitivities of populations in the field.

A QSAR-ICE-SSD model prediction of the PNECs for alkylphenol substances and application in ecological risk assessment for rivers of a megacity

Alkylphenols (APs) are ubiquitous and generally present in higher residue levels in the environment. The present work focuses on the development of a set of in silico models to predict the aquatic toxicity of APs with incomplete/unknown toxicity data in aquatic environments. To achieve this, a QSAR-ICE-SSD model was constructed for aquatic organisms by combining quantitative structure-activity relationship (QSAR), interspecies correlation estimation (ICE), and species sensitivity distribution (SSD) models in order to obtain the hazardous concentrations (HCs) of selected APs. The research indicated that the keywords "alkylphenol" and "nonylphenol" were most commonly studied. The selected ICE models were robust (R2: 0.70-0.99; p-value < 0.01). All models had a high reliability cross- validation success rates (>75%), and the HC5 predicted with the QSAR-ICE-SSD model was 2-fold than that derived with measured experimental data. The HC5 values demonstrated nearly linear decreasing trend from 2-MP to 4-HTP, while the decreasing trend from 4-HTP to 4-DP became shallower, indicates that the toxicity of APs to aquatic organisms increases with the addition of alkyl carbon chain lengths. The ecological risks assessment (ERA) of APs revealed that aquatic organisms were at risk from exposure to 4-NP at most river stations (the highest risk quotient (RQ) = 1.51), with the highest relative risk associated with 2.9% of 4-NP detected in 82.9% of the sampling sites. The targeted APs posed potential ecological risks in the Yongding and Beiyun River according to the mixture ERA. The potential application of QSAR-ICE-SSD models could satisfy the immediate needs for HC5 derivations without the need for additional in vivo testing.

ECORISK2050: An Innovative Training Network for predicting the effects of global change on the emission, fate, effects, and risks of chemicals in aquatic ecosystems

By 2050, the global population is predicted to reach nine billion, with almost three quarters living in cities. The road to 2050 will be marked by changes in land use, climate, and the management of water and food across the world. These global changes (GCs) will likely affect the emissions, transport, and fate of chemicals, and thus the exposure of the natural environment to chemicals. ECORISK2050 is a Marie Skłodowska-Curie Innovative Training Network that brings together an interdisciplinary consortium of academic, industry and governmental partners to deliver a new generation of scientists, with the skills required to study and manage the effects of GCs on chemical risks to the aquatic environment. The research and training goals are to: (1) assess how inputs and behaviour of chemicals from agriculture and urban environments are affected by different environmental conditions, and how different GC scenarios will drive changes in chemical risks to human and ecosystem health; (2) identify short-to-medium term adaptation and mitigation strategies, to abate unacceptable increases to risks, and (3) develop tools for use by industry and policymakers for the assessment and management of the impacts of GC-related drivers on chemical risks. This project will deliver the next generation of scientists, consultants, and industry and governmental decision-makers who have the knowledge and skillsets required to address the changing pressures associated with chemicals emitted by agricultural and urban activities, on aquatic systems on the path to 2050 and beyond.

(Eco)toxicological tests for assessing impacts of chemical stress to aquatic ecosystems: Facts, challenges, and future

Monitoring of chemicals in the aquatic environment by chemical analysis alone cannot completely assess and predict the effects of chemicals on aquatic species and ecosystems. This is primarily because of the increasing number of (unknown) chemical stressors and mixture effects present in the environment. In addition, the ability of ecological indices to identify underlying stressors causing negative ecological effects is limited. Therefore, additional complementary methods are needed that can address the biological effects in a direct manner and provide a link to chemical exposure, i.e. (eco)toxicological tests. (Eco)toxicological tests are defined as test systems that expose biological components (cells, individuals, populations, communities) to (environmental mixtures of) chemicals to register biological effects. These tests measure responses at the sub-organismal (biomarkers and in vitro bioassays), whole-organismal, population, or community level. We performed a literature search to obtain a state-of-the-art overview of ecotoxicological tests available for assessing impacts of chemicals to aquatic biota and to reveal datagaps. In total, we included 509 biomarkers, 207 in vitro bioassays, 422 tests measuring biological effects at the whole-organismal level, and 78 tests at the population- community- and ecosystem-level. Tests at the whole-organismal level and biomarkers were most abundant for invertebrates and fish, whilst in vitro bioassays are mostly based on mammalian cell lines. Tests at the community- and ecosystem-level were almost missing for organisms other than microorganisms and algae. In addition, we provide an overview of the various extrapolation challenges faced in using data from these tests and suggest some forward looking perspectives. Although extrapolating the measured responses to relevant protection goals remains challenging, the combination of ecotoxicological experiments and models is key for a more comprehensive assessment of the effects of chemical stressors to aquatic ecosystems.

Pharmaceuticals and other urban contaminants threaten Amazonian freshwater ecosystems

Urban areas in the Brazilian Amazon have grown at an unprecedented rate during the last years. About 90% of the wastewater produced by these urban areas are discharged untreated into Amazonian freshwater ecosystems, constituting a potential environmental pathway for pharmaceuticals and other chemicals consumed by modern societies (e.g. psychostimulants, personal-care products, hormones). The distribution of these chemicals into the Amazon River and their potential risks for freshwater biodiversity have not been evaluated so far. Here, we show the results of the largest chemical monitoring campaign conducted in the Amazon region. We assessed exposure patterns for 43 pharmaceuticals and other urban contaminants in 40 sampling sites distributed along the Amazon River, three major tributaries (Negro, Tapajós and Tocantins Rivers), and four large cities of the Brazilian Amazon (Manaus, Santarém, Macapá, Belém). We assessed risks for freshwater biodiversity using species sensitivity distributions and mixture toxicity approaches. We found that urban areas constitute important hot-spots for chemical contamination, with mixtures containing up to 40 different compounds and exposure concentrations reaching the world's maxima for some of them. We show that chemical pollution can result in long-term effects for up to 50-80% of aquatic species next to urban areas. Moreover, we identified several ubiquitous compounds which can be used as tracers of anthropogenic pressure in the Amazon basin. We conclude that the chemical burden created by urbanization significantly contributes to a biodiversity loss in the region and should be further controlled.

Litter decomposition can be reduced by pesticide effects on detritivores and decomposers: Implications for tropical stream functioning

Understanding which factors affect the process of leaf litter decomposition is crucial if we are to predict changes in the functioning of stream ecosystems as a result of human activities. One major activity with known consequences on streams is agriculture, which is of particular concern in tropical regions, where forests are being rapidly replaced by crops. While pesticides are potential drivers of reduced decomposition rates observed in agricultural tropical streams, their specific effects on the performance of decomposers and detritivores are mostly unknown. We used a microcosm experiment to examine the individual and joint effects of an insecticide (chlorpyrifos) and a fungicide (chlorothalonil) on survival and growth of detritivores (Anchytarsus, Hyalella and Lepidostoma), aquatic hyphomycetes (AH) sporulation rate, taxon richness, assemblage structure, and leaf litter decomposition rates. Our results revealed detrimental effects on detritivore survival (which were mostly due to the insecticide and strongest for Hyalella), changes in AH assemblage structure, and reduced sporulation rate, taxon richness and microbial decomposition (mostly in response to the fungicide). Total decomposition was reduced especially when the pesticides were combined, suggesting that they operated differently and their effects were additive. Importantly, effects on decomposition were greater for single-species detritivore treatments than for the 3-species mixture, indicating that detritivore species loss may exacerbate the consequences of pesticides of stream ecosystem functioning.

Role of Toxicokinetic and Toxicodynamic Parameters in Explaining the Sensitivity of Zebrafish Larvae to Four Metals

Given the persistence and toxic potencies of metal contaminants in ecosystems, animals, and human beings, they are considered to be hazardous global pollutants. While the lethality of metal toxicities (e.g., LC₅₀) can significantly vary, even within the same species, the underlying mechanisms are less well-understood. In this study, we developed a subcellular two-compartment toxicokinetic-toxicodynamic (TK-TD) model for zebrafish larvae when exposed to four metals (cadmium, lead, copper, and zinc) to reveal whether differences in metal toxicity (LC₅₀ values) were dominated by the TK or TD processes. Results showed that the subcellular TK and TD parameters of the four metals were significantly different, and the bioconcentration factor (BCF) value of copper was higher than those of the other metals. We also found that the TD parameter internal threshold concentration (C_IT) was significantly positively correlated to the LC₅₀ values (R² = 0.7), suggesting a dominant role of TD processes in metal toxicity. Furthermore, the combined parameter C_IT/BCF for a metal-sensitive fraction (BCF_MSF), which linked exposure to effects through the TK-TD approach, explained up to 89% of the variation in toxicity to the four metals. The present study suggests that the observed variation in toxicity of these four metals was mainly determined by TD processes but that TK processes should not be ignored, especially for copper.

Effects of anthropogenic pollution and hydrological variation on macroinvertebrates in Mediterranean rivers: A case-study in the upper Tagus river basin (Spain)

Seasonal hydrological variation and chemical pollution represent two main drivers of freshwater biodiversity change in Mediterranean rivers. We investigated to what extent low flow conditions can modify the effects of chemical pollution on macroinvertebrate communities. To that purpose, we selected twelve sampling sites in the upper Tagus river basin (central Spain) having different sources of chemical pollution and levels of seasonal hydrological variation. The sites were classified as natural (high flow variation, low chemical impact), agricultural (high flow variation, high agricultural chemical inputs) and urban (limited flow variation, high urban chemical inputs). In these sites, we measured daily water discharge, nutrients, and contaminant concentrations, and we sampled benthic macroinvertebrates, in spring, summer and autumn. Significant differences related to toxic pressure and nutrient concentrations were observed between the three groups of sites. Seasonal patterns were found for some water quality parameters (e.g. nitrites, ammonia, suspended solids, metal toxicity), particularly in agricultural sites. Taxonomic and functional richness were slightly lower in the polluted sites (agricultural and urban), particularly during low flow periods (summer and early autumn). Functional diversity was significantly lower in sites with seasonal flow variation (agricultural sites) as compared to the more constant ones (urban sites). The frequency of traits such as large size, asexual reproduction, aquatic passive dispersion and the production of cocoons increased in response to pollution during low flow periods. This study shows that the impacts of anthropogenic chemical pollution on taxonomic and functional characteristics of macroinvertebrate communities seem to be larger during low flow periods. Therefore, further studies and monitoring campaigns assessing the effects of chemical pollution within these periods are recommended.

High sensitivity of invertebrate detritivores from tropical streams to different pesticides

Freshwater organisms are often sensitive to pesticides, but their sensitivity varies across different taxa and with pesticide type and action mode, as shown by multiple acute toxicity tests. Such variability hampers predictions about how freshwater ecosystems may be altered by pesticide toxicity, which is especially critical for understudied areas of the world such as the tropics. Furthermore, there is little information about the sensitivity of some organisms that are key components of stream food webs; this is the case of litter-feeding detritivorous invertebrates, which contribute to the fundamental process of litter decomposition. Here, we examined the sensitivity of three common detritivores [Anchytarsus sp. (Coleoptera: Ptilodactylidae), Hyalella sp. (Amphipoda: Hyalellidae) and Lepidostoma sp. (Trichoptera: Lepidostomatidae)] to three pesticides commonly used (the insecticides bifenthrin and chlorpyrifos and the fungicide chlorothalonil) using acute (48 or 96 h) toxicity tests. Our study demonstrates that common-use pesticides provoke the mortality of half their populations at concentrations of 0.04-2.7 μg L-1. We found that all species were sensitive to the three pesticides, with the highest sensitivity found for chlorpyrifos. Additionally, we used the approach of species sensitivity distributions (SSD) to compare our study species with Daphnia magna and other temperate and tropical invertebrates. We found that the study species were among the most sensitive species to chlorpyrifos and chlorothalonil. Our results suggest that tropical detritivores merit special attention in ecological risk assessment of pesticides and highlight the need for accurate ecotoxicological information from ecologically relevant species in the tropics.

Off-Target Stoichiometric Binding Identified from Toxicogenomics Explains Why Some Species Are More Sensitive than Others to a Widely Used Neonicotinoid

Neonicotinoids are currently licensed for use in 120 countries, making accurate nontarget species sensitivity predictions critical. Unfortunately, such predictions are fraught with uncertainty, as sensitivity is extrapolated from only a few test species and neonicotinoid sensitivities can differ greatly between closely related taxa. Combining classical toxicology with de novo toxicogenomics could greatly improve sensitivity predictions and identify unexpectedly susceptible species. We show that there is a >30-fold differential species sensitivity (DSS) for the neonicotinoid imidacloprid between five earthworm species, a critical nontarget taxon. This variation could not be explained by differential toxicokinetics. Furthermore, comparing key motif expression in subunit genes of the classical nicotinic acetylcholine receptor (nAChR) target predicts only minor differences in the ligand binding domains (LBDs). In contrast, predicted dissimilarities in LBDs do occur in the highly expressed but nonclassical targets, acetylcholine binding proteins (AChBPs). Critically, the predicted AChBP divergence is capable of explaining DSS. We propose that high expression levels of putative nonsynaptic AChBPs with high imidacloprid affinities reduce imidacloprid binding to critical nAChRs involved in vital synaptic neurotransmission. This study provides a clear example of how pragmatic interrogation of key motif expression in complex multisubunit receptors can predict observed DSS, thereby informing sensitivity predictions for essential nontarget species.

Double constrained ordination for assessing biological trait responses to multiple stressors: A case study with benthic macroinvertebrate communities

Benthic macroinvertebrate communities are used as indicators for anthropogenic stress in freshwater ecosystems. To better understand the relationship between anthropogenic stress and changes in macroinvertebrate community composition, it is important to understand how different stressors and species traits are associated, and how these associations influence variation in species occurrence and abundances. Here, we show the capacity of the multivariate technique of double constrained correspondence analysis (dc-CA) to analyse trait-environment relationships, and we compare it with the redundancy analysis method on community weighted mean values of traits (CWM-RDA), which is frequently used for this type of analysis. The analyses were based on available biomonitoring data for macroinvertebrate communities from the Danube River. Results from forward selection of traits and environmental variables using dc-CA analyses showed that aquatic stages, reproduction techniques, dispersal tactics, locomotion and substrate relations, altitude, longitudinal and transversal distribution, and substrate preferendum were significantly related to habitat characteristics, hydromorphological alterations and water quality measurements such as physico-chemical parameters, heavy metals, pesticides and pharmaceuticals. Environmental variables significantly associated with traits using the CWM-RDA method were generally consistent with those found in dc-CA analysis. However, the CWM-RDA does neither test nor explicitly select traits, while dc-CA tests and selects both traits and environmental variables. Moreover, the dc-CA analysis revealed that the set of environmental variables was much better in explaining the community data than the available trait set, a kind of information that can neither be obtained from CWM-RDA nor from RLQ (Environment, Link and Trait data), which is a close cousin of dc-CA but not regression-based. Our results suggest that trait-based analysis based on dc-CA may be useful to assess mechanistic links between multiple anthropogenic stressors and ecosystem health, but more data sets should be analysed in the same manner.

Cross-species extrapolation of chemical sensitivity

Ecosystems are usually populated by many species. Each of these species carries the potential to show a different sensitivity towards all of the numerous chemical compounds that can be present in their environment. Since experimentally testing all possible species-chemical combinations is impossible, the ecological risk assessment of chemicals largely depends on cross-species extrapolation approaches. This review overviews currently existing cross-species extrapolation methodologies, and discusses i) how species sensitivity could be described, ii) which predictors might be useful for explaining differences in species sensitivity, and iii) which statistical considerations are important. We argue that risk assessment can benefit most from modelling approaches when sensitivity is described based on ecologically relevant and robust effects. Additionally, specific attention should be paid to heterogeneity of the training data (e.g. exposure duration, pH, temperature), since this strongly influences the reliability of the resulting models. Regarding which predictors are useful for explaining differences in species sensitivity, we review interspecies-correlation, relatedness-based, traits-based, and genomic-based extrapolation methods, describing the amount of mechanistic information the predictors contain, the amount of input data the models require, and the extent to which the different methods provide protection for ecological entities. We develop a conceptual framework, incorporating the strengths of each of the methods described. Finally, the discussion of statistical considerations reveals that regardless of the method used, statistically significant models can be found, although the usefulness, applicability, and understanding of these models varies considerably. We therefore recommend publication of scientific code along with scientific studies to simultaneously clarify modelling choices and enable elaboration on existing work. In general, this review specifies the data requirements of different cross-species extrapolation methods, aiming to make regulators and publishers more aware that access to raw- and meta-data needs to be improved to make future cross-species extrapolation efforts successful, enabling their integration into the regulatory environment.

Effect Modeling Quantifies the Difference Between the Toxicity of Average Pesticide Concentrations and Time-Variable Exposures from Water Quality Monitoring

Synthetic chemicals are frequently detected in water bodies, and their concentrations vary over time. Water monitoring programs typically employ either a sequence of grab samples or continuous sampling, followed by chemical analysis. Continuous time-proportional sampling yields the time-weighted average concentration, which is taken as proxy for the real, time-variable exposure. However, we do not know how much the toxicity of the average concentration differs from the toxicity of the corresponding fluctuating exposure profile. We used toxicokinetic-toxicodynamic models (invertebrates, fish) and population growth models (algae, duckweed) to calculate the margin of safety in moving time windows across measured aquatic concentration time series (7 pesticides) in 5 streams. A longer sampling period (14 d) for time-proportional sampling leads to more deviations from the real chemical stress than shorter sampling durations (3 d). The associated error is a factor of 4 or less in the margin of safety value toward underestimating and an error of factor 9 toward overestimating chemical stress in the most toxic time windows. Under- and overestimations occur with approximate equal frequency and are very small compared with the overall variation, which ranged from 0.027 to 2.4 × 10¹⁰ (margin of safety values). We conclude that continuous, time-proportional sampling for a period of 3 and 14 d for acute and chronic assessment, respectively, yields sufficiently accurate average concentrations to assess ecotoxicological effects. Environ Toxicol Chem 2020;39:2158-2168. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Potential impact of chemical stress on freshwater invertebrates: A sensitivity assessment on continental and national scale based on distribution patterns, biological traits, and relatedness

Current chemical risk assessment approaches rely on a standard suite of test species to assess toxicity to environmental species. Assessment factors are used to extrapolate from single species to communities and ecosystem effects. This approach is pragmatic, but lacks resolution in biological and environmental parameters. Novel modelling approaches can help improve the biological resolution of assessments by using mechanistic information to identify priority species and priority regions that are potentially most impacted by chemical stressors. In this study we developed predictive sensitivity models by combining species-specific information on acute chemical sensitivity (LC50 and EC50), traits, and taxonomic relatedness. These models were applied at two spatial scales to reveal spatial differences in the sensitivity of species assemblages towards two chemical modes of action (MOA): narcosis and acetylcholinesterase (AChE) inhibition. We found that on a relative scale, 46% and 33% of European species were ranked as more sensitive towards narcosis and AChE inhibition, respectively. These more sensitive species were distributed with higher occurrences in the south and north-eastern regions, reflecting known continental patterns of endemic macroinvertebrate biodiversity. We found contradicting sensitivity patterns depending on the MOA for UK scenarios, with more species displaying relative sensitivity to narcotic MOA in north and north-western regions, and more species with relative sensitivity to AChE inhibition MOA in south and south-western regions. Overall, we identified hotspots of species sensitive to chemical stressors at two spatial scales, and discuss data gaps and crucial technological advances required for the successful application of the proposed methodology to invertebrate scenarios, which remain underrepresented in global conservation priorities.

Standartox: Standardizing Toxicity Data

An increasing number of chemicals such as pharmaceuticals, pesticides and synthetic hormones are in daily use all over the world. In the environment, chemicals can adversely affect populations and communities and in turn related ecosystem functions. To evaluate the risks from chemicals for ecosystems, data on their toxicity, which are typically produced in standardized ecotoxicological laboratory tests, is required. The results from ecotoxicological tests are compiled in (meta-)databases such as the United States Environmental Protection Agency (EPA) ECOTOXicology Knowledgebase (ECOTOX). However, for many chemicals, multiple ecotoxicity data are available for the same test organism. These can vary strongly, thereby causing uncertainty of related analyses. Given that most current databases lack aggregation steps or are confined to specific chemicals, we developed Standartox, a tool and database that continuously incorporates the ever-growing number of test results in an automated process workflow that ultimately leads to a single aggregated data point for a specific chemical-organism test combination, representing the toxicity of a chemical. Standartox can be accessed through a web application and an R package.

Reliable and representative in silico predictions of freshwater ecotoxicological hazardous concentrations

A reliable quantification of the potential effects of chemicals on freshwater ecosystems requires ecotoxicological response data for a large set of species which is typically not available in practice. In this study, we propose a method to estimate hazardous concentrations (HCs) of chemicals on freshwater ecosystems by combining two in silico approaches: quantitative structure activity relationships (QSARs) and interspecies correlation estimation (ICE) models. We illustrate the principle of our QSAR-ICE method by quantifying the HCs of 51 chemicals at which 50% and 5% of all species are exposed above the concentration causing acute effects. We assessed the bias of the HCs, defined as the ratio of the HC based on measured ecotoxicity data and the HC based on in silico data, as well as the statistical uncertainty, defined as the ratio of the 95th and 5th percentile of the HC. Our QSAR-ICE method resulted in a bias that was comparable to the use of measured data for three species, as commonly used in effect assessments: the average bias of the QSAR-ICE HC50 was 1.2 and of the HC5 2.3 compared to 1.2 when measured data for three species were used for both HCs. We also found that extreme statistical uncertainties (>10⁵) are commonly avoided in the HCs derived with the QSAR-ICE method compared to the use of three measurements with statistical uncertainties up to 10¹². We demonstrated the applicability of our QSAR-ICE approach by deriving HC50s for 1,223 out of the 3,077 organic chemicals of the USEtox database. We conclude that our QSAR-ICE method can be used to determine HCs without the need for additional in vivo testing to help prioritise which chemicals with no or few ecotoxicity data require more thorough assessment.

Smarter Sediment Screening: Effect-Based Quality Assessment, Chemical Profiling, and Risk Identification

Sediments play an essential role in the functioning of aquatic ecosystems but simultaneously retain harmful compounds. However, sediment quality assessment methods that consider the risks caused by the combined action of all sediment-associated contaminants to benthic biota are still underrepresented in water quality assessment strategies. Significant advancements have been made in the application of effect-based methods, but methodological improvements can still advance sediment risk assessment. The present study aimed to explore such improvements by integrating effect-monitoring and chemical profiling of sediment contamination. To this end, 28 day life cycle bioassays with Chironomus riparius using intact whole sediment cores from contaminated sites were performed in tandem with explorative chemical profiling of bioavailable concentrations of groups of legacy and emerging sediment contaminants to investigate ecotoxicological risks to benthic biota. All contaminated sediments caused effects on the resilient midge C. riparius, stressing that sediment contamination is ubiquitous and potentially harmful to aquatic ecosystems. However, bioassay responses were not in line with any of the calculated toxicity indices, suggesting that toxicity was caused by unmeasured compounds. Hence, this study underlines the relevance of effect-based sediment quality assessment and provides smarter ways to do so.

Shrinking Body Size and Physiology Contribute to Geographic Variation and the Higher Toxicity of Pesticides in a Warming World

To improve current and future risk assessment of pesticides under global warming, mechanistic insights and consideration of daily temperature fluctuations (DTFs) are needed. One overlooked mechanism how both higher mean temperatures and DTFs may increase toxicity is by reducing body size (temperature-size-rule). We studied whether a higher mean temperature and DTF magnified chlorpyrifos toxicity in Ischnura elegans damselfly larvae, and whether this was mediated by temperature-induced reductions in body size and/or physiological changes. The lethal effects of chlorpyrifos were magnified at the high mean temperature (up to ∼15%) and under DTF (up to ∼33%), and especially at their combination (up to ∼46%) indicating synergisms. This highlights that not only considering DTFs, but also their interaction with higher mean temperatures is pivotal for realistic predictions of pesticide toxicity. Both higher temperatures and DTFs resulted in smaller larvae, which were more sensitive to chlorpyrifos. Notably, the DTF-induced smaller body sizes, as well as the higher oxidative damage to lipids, contributed to the higher chlorpyrifos toxicity under DTF. By integrating the temperature-size rule and size-pesticide sensitivity pattern we provide proof-of-principle for a novel, likely general mechanism contributing to geographic variation and the higher toxicity of pesticides in a warming world.

Effects of multiple stressors associated with agriculture on stream macroinvertebrate communities in a tropical catchment

Tropical forests are declining at unprecedented rates in favour of agriculture, and streams can be severely impacted due to effects of multiple stressors that have rarely been considered together in tropical studies. We studied the effects of multiple stressors associated with agricultural practices (pesticide toxicity, nutrient enrichment and habitat alteration-quantified as TUmax, soluble reactive phosphorus concentration and sedimentation, respectively) on macroinvertebrate communities in a tropical catchment in Panama (13 stream sites sampled in 20 occasions from 2015 to 2017, with 260 samples in total). We examined how macroinvertebrate abundance, taxonomic richness, community composition and biotic indices (SPEAR and BMWP/PAN, which were specifically designed to detect pesticide toxicity and nutrient enrichment, respectively) varied depending on the studied stressors, considering their single and combined effects. Our analyses revealed significant effects of the studied stressors on macroinvertebrate communities, with two particular results that merit further attention: (1) the fact that pesticide toxicity affected BMWP/PAN, but not SPEAR, possibly because the former had been adapted for local fauna; and (2) that most stressors showed antagonistic interactions (i.e., lower combined effects than expected from their individual effects). These results highlight the need for toxicity bioassays with tropical species that allow adaptations of biotic indices, and of observational and manipulative studies exploring the combined effects of multiple stressors on tropical macroinvertebrate communities and ecosystems, in order to predict and manage future anthropogenic impacts on tropical streams.