Benchmark dose analysis framework for developing wildlife toxicity reference values

Methylmercury Effects on Birds: A Review, Meta-Analysis, and Development of Toxicity Reference Values for Injury Assessment Based on Tissue Residues and Diet

Birds are used as bioindicators of environmental mercury (Hg) contamination, and toxicity reference values are needed for injury assessments. We conducted a comprehensive review, summarized data from 168 studies, performed a series of Bayesian hierarchical meta-analyses, and developed new toxicity reference values for the effects of methylmercury (MeHg) on birds using a benchmark dose analysis framework. Lethal and sublethal effects of MeHg on birds were categorized into nine biologically relevant endpoint categories and three age classes. Effective Hg concentrations where there was a 10% reduction (EC10) in the production of juvenile offspring (0.55 µg/g wet wt adult blood-equivalent Hg concentrations, 80% credible interval: [0.33, 0.85]), histology endpoints (0.49 [0.15, 0.96] and 0.61 [0.09, 2.48]), and biochemical markers (0.77 [<0.25, 2.12] and 0.57 [0.35, 0.92]) were substantially lower than those for survival (2.97 [2.10, 4.73] and 5.24 [3.30, 9.55]) and behavior (6.23 [1.84, >13.42] and 3.11 [2.10, 4.64]) of juveniles and adults, respectively. Within the egg age class, survival was the most sensitive endpoint (EC10 = 2.02 µg/g wet wt adult blood-equivalent Hg concentrations [1.39, 2.94] or 1.17 µg/g fresh wet wt egg-equivalent Hg concentrations [0.80, 1.70]). Body morphology was not particularly sensitive to Hg. We developed toxicity reference values using a combined survival and reproduction endpoints category for juveniles, because juveniles were more sensitive to Hg toxicity than eggs or adults. Adult blood-equivalent Hg concentrations (µg/g wet wt) and egg-equivalent Hg concentrations (µg/g fresh wet wt) caused low injury to birds (EC1) at 0.09 [0.04, 0.17] and 0.04 [0.01, 0.08], moderate injury (EC5) at 0.6 [0.37, 0.84] and 0.3 [0.17, 0.44], high injury (EC10) at 1.3 [0.94, 1.89] and 0.7 [0.49, 1.02], and severe injury (EC20) at 3.2 [2.24, 4.78] and 1.8 [1.28, 2.79], respectively. Maternal dietary Hg (µg/g dry wt) caused low injury to juveniles at 0.16 [0.05, 0.38], moderate injury at 0.6 [0.29, 1.03], high injury at 1.1 [0.63, 1.87], and severe injury at 2.4 [1.42, 4.13]. We found few substantial differences in Hg toxicity among avian taxonomic orders, including for controlled laboratory studies that injected Hg into eggs. Our results can be used to quantify injury to birds caused by Hg pollution. Environ Toxicol Chem 2024;43:1195-1241. Published 2024. This article is a U.S. Government work and is in the public domain in the USA. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Meta-metabolomic responses of river biofilms to cobalt exposure and use of dose-response model trends as an indicator of effects

The response of the meta-metabolome is rarely used to characterize the effects of contaminants on a whole community. Here, the meta-metabolomic fingerprints of biofilms were examined after 1, 3 and 7 days of exposure to five concentrations of cobalt (from background concentration to 1 × 10-5 M) in aquatic microcosms. The untargeted metabolomic data were processed using the DRomics tool to build dose-response models and to calculate benchmark-doses. This approach made it possible to use 100% of the chemical signal instead of being limited to the very few annotated metabolites (7%). These benchmark-doses were further aggregated into an empirical cumulative density function. A trend analysis of the untargeted meta-metabolomic feature dose-response curves after 7 days of exposure suggested the presence of a concentration range inducing defense responses between 1.7 × 10-9 and 2.7 × 10-6 M, and of a concentration range inducing damage responses from 2.7 × 10-6 M and above. This distinction was in good agreement with changes in the other biological parameters studied (biomass and chlorophyll content). This study demonstrated that the molecular defense and damage responses can be related to contaminant concentrations and represents a promising approach for environmental risk assessment of metals.

Wildlife ecological risk assessment in the 21st century: Promising technologies to assess toxicological effects

Despite advances in toxicity testing and the development of new approach methodologies (NAMs) for hazard assessment, the ecological risk assessment (ERA) framework for terrestrial wildlife (i.e., air-breathing amphibians, reptiles, birds, and mammals) has remained unchanged for decades. While survival, growth, and reproductive endpoints derived from whole-animal toxicity tests are central to hazard assessment, nonstandard measures of biological effects at multiple levels of biological organization (e.g., molecular, cellular, tissue, organ, organism, population, community, ecosystem) have the potential to enhance the relevance of prospective and retrospective wildlife ERAs. Other factors (e.g., indirect effects of contaminants on food supplies and infectious disease processes) are influenced by toxicants at individual, population, and community levels, and need to be factored into chemically based risk assessments to enhance the "eco" component of ERAs. Regulatory and logistical challenges often relegate such nonstandard endpoints and indirect effects to postregistration evaluations of pesticides and industrial chemicals and contaminated site evaluations. While NAMs are being developed, to date, their applications in ERAs focused on wildlife have been limited. No single magic tool or model will address all uncertainties in hazard assessment. Modernizing wildlife ERAs will likely entail combinations of laboratory- and field-derived data at multiple levels of biological organization, knowledge collection solutions (e.g., systematic review, adverse outcome pathway frameworks), and inferential methods that facilitate integrations and risk estimations focused on species, populations, interspecific extrapolations, and ecosystem services modeling, with less dependence on whole-animal data and simple hazard ratios. Integr Environ Assess Manag 2024;20:725-748. © 2023 His Majesty the King in Right of Canada and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). Reproduced with the permission of the Minister of Environment and Climate Change Canada. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Integrating emerging science to improve estimates of risk to wildlife from chemical exposure: What are the challenges?

Many jurisdictions require ecological risk assessments for terrestrial wildlife (i.e., terrestrial vertebrates) to assess potential adverse effects from exposure to anthropogenic chemicals. This occurs, for example, at contaminated sites and when new pesticides are proposed, and it occurs for chemicals that are in production and/or proposed for wide-scale use. However, guidance to evaluate such risks has not changed markedly in decades, despite the availability of new scientific tools to do so. In 2019, the Wildlife Toxicology World Interest Group of the Society of Environmental Toxicology and Chemistry (SETAC) initiated a virtual workshop that included a special session coincident with the annual SETAC North America meeting and which focused on the prospect of improving risk assessments for wildlife and improving their use in implementing chemical regulations. Work groups continued the work and investigated the utility of integrating emerging science and novel methods for improving problem formulation (WG1), exposure (WG2), toxicology (WG3), and risk characterization (WG4). Here we provide a summary of that workshop and the follow-up work, the regulations that drive risk assessment, and the key focus areas identified to advance the ability to predict risks of chemicals to wildlife. Integr Environ Assess Manag 2024;20:645-657. © 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Using radish (Raphanus lativus L.) germination to establish a benchmark dose for the toxicity of ozonated-petroleum byproducts in soil

The concentration-response relationship between the germination outcome of radish (Raphanus lativus L.) and ozonated petroleum residuals was determined experimentally. The outcomes were used to produce an ecological risk assessment model to predict the extra risk of adverse outcomes based on the concentration of ozonated residuals. A test soil with low organic matter (0.5% w/w) was mixed with raw crude oil, artificially weathered, and treated at three doses of ozone (O₃) gas (5 g, 10 g, and 40 g O₃ per 600 g of soil). Total petroleum hydrocarbons (TPH) and produced dissolved organic carbon (DOC) were measured. TREATMENT categories (control, petroleum, petroleum + 5 g O₃, petroleum + 10 g O₃, and petroleum + 40 g O₃) were then used to create a dilution series using different proportions of the test soil and a commercially available potting mix (∼75% w/w organic matter) to evaluate the effects of background organic matter (b-ORGANIC) in conjunction with TPH and DOC. Multivariable logistic regression was performed on the adverse germination outcome as a function of TPH, DOC, TREATMENT, and b-ORGANIC. The parameters controlling germination were the continuous variable DOC and the categorical variables TREATMENT and b-ORGANIC. Radish germination was strongly harmed by DOC from ozonation, but DOC's ecotoxicity decreased with increasing O₃ dose and the presence of b-ORGANIC beyond 10% (w/w). We used the germination outcome of radish to produce a logistic regression model that computes margins of DOC (± std. error) that create 10%, 25%, and 50% extra risk of adverse germination effects.

Reproductive and developmental toxicity of perfluorooctane sulfonate (PFOS) in the white-footed mouse (Peromyscus leucopus)

Concerns about per- and polyfluoroalkyl substances (PFAS) stem from their ubiquitous presence in the environment, bioaccumulation, resistance to degradation, and toxicity. Previously, toxicity data relevant to ecological risk assessment has largely been aquatic, terrestrial invertebrates, or avian in origin. In this study, repeated oral exposures of perfluorooctane sulfonate (PFOS) were administered to white-footed mice (Peromyscus leucopus) to evaluate effects on reproduction and development. Prenatal exposure to high doses of PFOS caused neonatal mortality, though growth and development were unaffected by low doses. Additionally, parental (P) generation animals exhibited increased liver:body weight, increased hepatocyte cytoplasmic vacuolization, and decreased serum thyroxine (T4) levels. Total litter loss was selected as the protective critical effect in this study resulting in a benchmark dose low (BMDL) of 0.12 mg/kg-d PFOS. Importantly, PFOS exposure has been linked to reduced adult recruitment in myriad species and at similar thresholds to this study. Similarities in critical/toxicologic effects across taxa may add confidence in risk assessments at sites with multiple taxa or environments.

Benchmark dose approach in investigating the relationship between blood metal levels and reproductive hormones: Data set from human study

The main objective of this research was to conduct a dose-response modeling between the internal dose of measured blood Cd, As, Hg, Ni, and Cr and hormonal response of serum testosterone, luteinizing hormone (LH), and follicle-stimulating hormone (FSH). The study included 207 male participants from subjects of 5 different cohorts (patients with prostate, testicular, and pancreatic cancer, patients suffering from various thyroid and metabolic disorders, as well as healthy volunteers), enrolled from January 2019 to May 2021 at the Clinical Centre of Serbia in Belgrade, Serbia. Benchmark dose-response modeling analysis was performed with the PROAST software version 70.1, showing the hormone levels as quantal data. The averaging technique was applied to compute the Benchmark dose (BMD) interval (BMDI), with benchmark response set at 10%. Dose-response relationships between metal/metalloid blood concentration and serum hormone levels were confirmed for all the investigated metals/metalloid and hormones. The narrowest BMDI was found for Cd-testosterone and Hg-LH pairs, indicative of high confidence in these estimates. Although further research is needed, the observed findings demonstrate that the BMD approach may prove to be significant in the dose-response modeling of human data.

A Simplified Population-Level Landscape Model Identifying Ecological Risk Drivers of Pesticide Applications, Part One: Case Study for Large Herbivorous Mammals

Environmental risk assessment is a key process for the authorization of pesticides, and is subjected to continuous challenges and updates. Current approaches are based on standard scenarios and independent substance-crop assessments. This arrangement does not address the complexity of agricultural ecosystems with mammals feeding on different crops. This work presents a simplified model for regulatory use addressing landscape variability, co-exposure to several pesticides, and predicting the effect on population abundance. The focus is on terrestrial vertebrates and the aim is the identification of the key risk drivers impacting on mid-term population dynamics. The model is parameterized for EU assessments according to the European Food Safety Authority (EFSA) Guidance Document, but can be adapted to other regulatory schemes. The conceptual approach includes two modules: (a) the species population dynamics, and (b) the population impact of pesticide exposure. Population dynamics is modelled through daily survival and seasonal reproductions rates; which are modified in case of pesticide exposure. All variables, parameters, and functions can be modified. The model has been calibrated with ecological data for wild rabbits and brown hares and tested for two herbicides, glyphosate and bromoxynil, using validated toxicity data extracted from EFSA assessments. Results demonstrate that the information available for a regulatory assessment, according to current EU information requirements, is sufficient for predicting the impact and possible consequences at population dynamic levels. The model confirms that agroecological parameters play a key role when assessing the effect of pesticide exposure on population abundance. The integration of laboratory toxicity studies with this simplified landscape model allows for the identification of conditions leading to population vulnerability or resilience. An Annex includes a detailed assessment of the model characteristics according to the EFSA scheme on Good Modelling Practice.

Cadmium levels in human breast tissue and estradiol serum levels: Is there a connection?

Cadmium (Cd), one of the most abundant environmental pollutants, is considered to have endocrine disrupting properties. However, data on the dose-response relationship between Cd dose and levels of hormones have been insufficiently studied, especially in human data sets. Thus, the aim of this study was to determine the possibility of analyzing data obtained from a case-control study in female patients with benign/malignant breast tumors, using the Benchmark dose (BMD) concept. The collected data on Cd levels in breast tissue and estrogen serum levels were processed in PROAST software using different variables. The dose-response relationship between the internal dose of Cd and estradiol levels in the serum was investigated and BMD intervals were calculated. The dose-response relationship between the Cd concentration in breast tissue and the estradiol serum level was shown, indicating lower estradiol serum levels as a consequence of higher Cd concentrations in breast tissue. As one of the few studies analyzing human data using the BMD approach, these findings could have a pivotal role in dose response analysis of data collected from human studies.

Transcriptomic effects of tributyltin (TBT) in zebrafish eleutheroembryos. A functional benchmark dose analysis

Exposure to the antifouling tributyltin (TBT) has been related to imposex in mollusks and to obesogenicity, adipogenesis and masculinization in fish. To understand the underlying molecular mechanisms, we evaluated dose-response effects of TBT (1.7-56 nM) in zebrafish eleutheroembryos transcriptome exposed from 2 to 5 days post-fertilization. RNA-sequencing analysis identified 3238 differentially expressed transcripts in eleutheroembryos exposed to TBT. Benchmark dose analyses (BMD) showed that the point of departure (PoD) for transcriptomic effects (9.28 nM) was similar to the metabolomic PoD (11.5 nM) and about one order of magnitude lower than the morphometric PoD (67.9 nM) or the median lethal concentration (LC₅₀: 93.6 nM). Functional analysis of BMD transcriptomic data identified steroid metabolism and cholesterol and vitamin D₃ biosynthesis as the most sensitive pathways to TBT (<50% PoD). Conversely, transcripts related to general stress and DNA damage became affected only at doses above the PoD. Therefore, our results indicate that transcriptomes can act as early molecular indicators of pollutant exposure, and illustrates their usefulness for the mechanistic identification of the initial toxic events. As the estimated molecular PoDs are close to environmental levels, we concluded that TBT may represent a substantial risk in some natural environments.