Making the Biotic Ligand Model kinetic, easier to develop, and more flexible for deriving water quality criteria

Quantitative differentiation of toxicity contributions and predicted global risk of fipronil and its transformation products to aquatic invertebrates

Biotransformation often alters chemical toxicity, yet its impacts on risk assessment are hardly quantified due to the challenges in acquiring internal exposure-based thresholds for chemicals that are readily metabolizable. Here, we integrated toxic unit and toxicokinetics to quantitatively assess toxicity contributions and potential risk of both parent compound and transformation products (TPs) to aquatic organisms, using fipronil (FIP) as a representative toxicant. In aquatic invertebrates Chironomus dilutus and Hyalella azteca, approximately 90 % of FIP was transformed to fipronil sulfone (SUL). FIP and SUL exhibited similar intrinsic toxicity to these organisms, which was contrary to conventional perception that SUL was more toxic than FIP. However, biotransformation was still important in risk assessment because the TP had 10-fold slower depuration rate than FIP. The amphipod H. azteca was found to be as sensitive to FIPs as the insect C. dilutus, which was previously considered ten times more sensitive based on external thresholds. This discrepancy has led to overlooking the toxicity of FIP to H. azteca in regional risk assessments. Lastly, we predicted the lethal risk of FIPs in global surface water. When using external thresholds for prediction, FIPs in 3.4 % of the water samples were lethally toxic to H. azteca, and the percentage of water samples at risk increased to 14 % when internal thresholds were used and SUF dominated the risk. This study presents an improved method for quantifying aquatic risk of readily metabolized toxicants. Our findings underscore the urgency of considering TPs in water quality assessments, especially for sensitive species that are at risk in the environment.

Turbid Waters and Clearer Standards: Refining Water Quality Criteria for Coastal Environments by Encompassing Metal Bioavailability from Suspended Particles

Suspended particulate matter (SPM) carries a major fraction of metals in turbid coastal waters, markedly influencing metal bioaccumulation and posing risks to marine life. However, its effects are often overlooked in current water quality criteria for metals, primarily due to challenges in quantifying SPM's contribution. This contribution depends on the SPM concentration, metal distribution coefficients (Kd), and the bioavailability of SPM-bound metals (assimilation efficiency, AE), which can collectively be integrated as a modifying factor (MF). Accordingly, we developed a new stable isotope method to measure metal AE by individual organisms from SPM, employing the widely distributed filter-feeding clam Ruditapes philippinarum as a representative species. Assessing SPM from 23 coastal sites in China, we found average AEs of 42% for Zn, 26% for Cd, 20% for Cu, 8% for Ni, and 6% for Pb. Moreover, using stable isotope methods, we determined metal Kd of SPM from these sites, which can be well predicted by the total organic carbon and iron content (R2 = 0.977). We calculated MFs using a Monte Carlo method. The calculated MFs are in the range 9.9-43 for Pb, 8.5-37 for Zn, 2.9-9.7 for Cu, 1.4-2.7 for Ni, and 1.1-1.6 for Cd, suggesting that dissolved-metal-based criteria values should be divided by MFs to provide adequate protection to aquatic life. This study provides foundational guidelines to refine water quality criteria in turbid waters and protect coastal ecosystems.

Heavy metals and metalloid in aquatic invertebrates: A review of single/mixed forms, combination with other pollutants, and environmental factors

Heavy metals (HMs) and metalloid occur naturally and are found throughout the Earth's crust but they are discharged into aquatic environments at high concentrations by human activities, increasing heavy metal pollution. HMs can bioaccumulate in higher organisms through the food web and consequently affect humans. In an aquatic environment, various HMs mixtures can be present. Furthermore, HMs adsorb on other environmental pollutants, such as microplastics and persistent organic pollutants, causing a synergistic or antagonistic effect on aquatic organisms. Therefore, to understand the biological and physiological effects of HMs on aquatic organisms, it is important to evaluate the effects of exposure to combinations of complex HM mixtures and/or pollutants and other environmental factors. Aquatic invertebrates occupy an important niche in the aquatic food chain as the main energy link between higher and lower organisms. The distribution of heavy metals and the resulting toxic effects in aquatic invertebrates have been extensively studied, but few reports have dealt with the relationship between HMs, pollutants, and environmental factors in biological systems with regard to biological availability and toxicity. This review describes the overall properties of individual HM and their effects on aquatic invertebrates and comprehensively reviews physiological and biochemical endpoints in aquatic invertebrates depending on interactions among HMs, other pollutants, and environmental factors.

WHAM-FTOXβ - An aquatic toxicity model based on intrinsic metal toxic potency and intrinsic species sensitivity

We developed a model that quantifies aquatic cationic toxicity by a combination of the intrinsic toxicities of metals and protons and the intrinsic sensitivities of the test species. It is based on the WHAM-F_TOX model, which combines the calculated binding of cations by the organism with toxicity coefficients (α_H, α_M) to estimate the variable F_TOX, a measure of toxic effect; the key parameter α_M,max (applying at infinite time) depends upon both the metal and the test species. In our new model, WHAM-F_TOXβ, values of α_M,max are given by the product α_M* × β, where α_M* has a single value for each metal, and β a single value for each species. To parameterise WHAM-F_TOXβ, we assembled a set of 2182 estimates of α_M,max obtained by applying the basic model to laboratory toxicity data for 76 different test species, covering 15 different metals, and including results for metal mixtures. Then we fitted the log₁₀α_M,max values with α_M* and β values (a total of 91 parameters). The resulting model accounted for 72% of the variance in log₁₀α_M,max. The values of α_M* increased markedly as the chemical character of the metal changed from hard (average α_M* = 4.4) to intermediate (average α_M* = 25) to soft (average α_M* = 560). The values of log₁₀β were normally distributed, with a 5-95 percentile range of -0.73 to +0.56, corresponding to β values of 0.18 to 3.62. The WHAM-F_TOXβ model entails the assumption that test species exhibit common relative sensitivity, i.e. the ratio α_M,max / α_M* is constant across all metals. This was tested with data from studies in which the toxic responses of a single organism towards two or more metals had been measured (179 examples for the most-tested metals Ni, Cu, Zn, Ag, Cd, Pb), and statistically-significant (p < 0.003) results were obtained.

Salinity-dependent aquatic life criteria of inorganic mercury in coastal water and its ecological risk assessment

Mercury (Hg) is one of the most toxic pollutants to aquatic organisms. The influence of salinity on Hg toxicity, an important factor restricting the development of global marine aquatic life criteria (ALC), is unclear. Therefore, mercury toxicity data were corrected based on salinity using the aggregate slope method, and the ALC values were derived. Short-term aquatic life criteria (SALC) and long-term aquatic life criteria (LALC) were derived using the species sensitivity distribution method based on Log-logistic, Log-normal, Burr III, Gumbel, and Weibull models. The hazard quotient (HQ) and joint probability curve (JPC) methods were used to evaluate the ecological risk of Hg in the coastal waters of China. The results showed that the SALC and LALC of Hg in the coastal waters of China were 2.21 and 0.54 μg/L. The toxicity data and salinity were positively correlated for Chordate and Arthropoda and negatively correlated for Mollusca. The SALC values increased by approximately 75%, with salinities ranging from 10 to 20 ppt. A slight peak in the SALC at mid-salinities was also observed. The ecological risk assessment of Hg in China's coastal waters showed that attention should be paid to Hg pollution in the Bohai Sea and East China Sea, especially the ecological risk of Hg to crustacean organisms. This study could promote the development of water quality criteria for coastal waters and provide a technical reference for mercury management in the coastal waters of China.

Effect of dissolved organic matter on the bioavailability and toxicity of cadmium in zebrafish larvae: Determination based on toxicokinetic-toxicodynamic processes

The presence of dissolved organic matter (DOM) strongly influences the bioavailability of metals in aquatic environments; however, the association between the binding activities and the concentrations of DOM compositions is not well documented, leading to uncertainties in metal toxicity assessment. We creatively quantify the mitigation and acceleration effects of DOM compositions on cadmium (Cd) bioaccumulation and toxicity in zebrafish larvae using abiotic ligand (ABLs) and biotic ligand (BLs) in a toxicokinetic-toxicodynamic (TK-TD) model. The BL-TK-TD model could accurately predict the protective effect of fulvic acid while overestimating the complexing capacity of citric acid. The model also could successfully simulate the protective effects of native DOM in most cases from 32 natural water bodies in China. The observed LC50 values of Cd showed a peak effect for the native DOM fraction comprising hydrophilic acidic contents (3.55 ± 0.44 mg L ^- 1) in natural water from 32 sites. The BL-TK-TD model provides practically useful information to identify the effect of different DOM compositions on metal bioavailability and toxicity in aquatic environments and guides future water management policies aimed at controlling aquatic heavy metal pollution.

Predicting Metal Bioavailability and Risk of Toxicity in Nigerian Surface Waters: Are the Existing User-Friendly Bioavailability Tools Applicable?

In the present study, we assessed the use of existing user-friendly bioavailability tools to predict metal bioavailability and the risk of toxicity in a typical Nigerian surface water. The effect of prevailing water chemistry on toxicity of metals was also assessed in laboratory studies. Surface water samples were collected at four locations downstream in the Ogun River in Lagos State, and water chemistry analysis was carried out using standard methods. Relevant parameters were inputted into the Bio-Met software and Metal Bioavailability Assessment Tool (M-Bat), which are run on a Microsoft Excel spreadsheet, to obtain the site-specific hazardous concentrations affecting 5% of the population (HC5)/predicted-no-effect concentrations (PNECs) for selected metals and risk characterization ratios (RCRs) at locations downstream in the river. Assessments using the bioavailability tools showed that two locations, Owode-Onirin and Kara, were more sensitive to Zn and Pb inputs compared to the other locations. The RCRs for Cu, Zn, and Pb, respectively, were >1 at all locations, indicating a potential risk of toxicity to aquatic life. Results from laboratory studies conformed with predictions from the bioavailability tools because Zn and Pb were more toxic to aquatic organisms in surface water collected from Owode-Onirin and Kara locations, respectively, compared to toxicity in surface water from other locations. The issue of bioavailability in metal pollution control is widely accepted and has been incorporated into water quality guidelines (WQGs) in the United States, the United Kingdom, and Europe. In the absence of regulations incorporating bioavailability in developing countries including Nigeria, it is essential to assess the applicability of existing methods incorporating bioavailability to surface waters in this region. Findings from such studies will facilitate the development of region-specific WQGs for metals which will be protective of locally relevant aquatic life. Environ Toxicol Chem 2022;41:2537-2547. © 2022 SETAC.

Investigating and modeling the toxicity of arsenate on wheat root elongation: Assessing the effects of pH, sulfate and phosphate

Excessive arsenic in soil and groundwater will not only seriously affect the growth of plants, but also endanger human health through the food chain. However, there are few studies on the effects of metalloid speciation and anion competition on the toxicity of arsenate [As(Ⅴ)]. To investigate the effects of accompanying anions and pH on the toxicity of As(Ⅴ) on wheat root elongation, wheat roots were exposed to the concentrations of As(Ⅴ) in the solution ranged from 0 to 500 mM and different levels of pH (4.5-8.0) and different accompanying anions (H₂PO₄^-, SO₄^2-, NO₃^- and Cl^-) for five days. The root length of wheat was measured and the biotic ligand model (BLM) was developed to predict the potential toxicity of As(V) speciation to wheat roots. The results illustrated that EC50 of total As(V) (EC50{As(Ⅴ)_T}) values increased from 6.88 to 33.9 μM with increasing pH values from 4.5 to 8.0, suggesting that increasing pH alleviated As(Ⅴ) toxicity. The EC50{AsO₄^3-} and EC50{HAsO₄^2-} values increased from 0.001 to 4342 μM and from 0.0214 to 27.4 μM, respectively, while the EC50{H₂AsO₄^-} and EC50{H₃AsO₄} values sharply decreased from 6.62 to 2.68 μM and from 41.8 μM to 5.34 nm, respectively, when pH increased from 4.5 to 8.0. The toxicity of As(Ⅴ) decreased as the H₂PO₄^- and SO₄^2- activities increased but not when the activities of NO₃^- and Cl^- increased, indicating that SO₄^2- and H₂PO₄^- showed competitive effects with As(Ⅴ) on the binding sites. Based on BLM theory, the stability constants were obtained: [Formula: see text] = 3.70; [Formula: see text] = 4.08; [Formula: see text] = 4.77; [Formula: see text] = 6.50; [Formula: see text] = 2.09 and [Formula: see text] = 1.86, with f_AsBL^50%= 0.30 and β = 1.73. Results implied that BLM performed well in As(Ⅴ) toxicity prediction when coupling toxic species AsO₄^3-, HAsO₄^2-, H₂AsO₄^-, and H₃AsO₄, and the competition of SO₄^2- and H₂PO₄^- for binding sites. The current study provides a useful tool to accurately predict As(V) toxicity to wheat roots.

The spectral characteristics and cadmium complexation of soil dissolved organic matter in a wide range of forest lands

The quality and quantity of dissolved organic matter (DOM) greatly controls the fate of heavy metals. The characteristics of DOM and its interaction with metals are essential for the metal ecological risk assessment of soils. In this study, the DOM spectral characteristics of representative forest soils and the complex capacities between fluorescent DOM components and cadmium (Cd) were analyzed. Functional groups, such as carboxylic acids, alcohols and phenols, were determined by FT-IR analysis. Chromophoric DOM, fluorescent DOM and dissolved organic carbon (DOC) concentrations exhibited strong correlations with each other, indicating that variations of DOC could be well explained by Chromophoric DOM or fluorescent DOM due to high correlation coefficients. The spectral slope ratio was in the range of 0.85-5.90, implying an abundance of heavy macromolecular humic acids, peptides, and polycondensates. The absorbance spectral at 254 nm (SUVA₂₅₄) strongly correlated with SUVA₂₆₀ (r = 0.992, P < 0.01), indicating that hydrophobicity closely related with aromatic structure, and aromatic groups could be broadly hydrophobic. Fluorescence indices were from 1.62 to 2.21 and biological index values ranged from 0.54 to 1.14, where the DOM was mainly sourced from mixed terrestrial and autogenous inputs in most sites. Four universal fluorescence components were identified and characterized by fluorescence EEM-PARAFAC, including two humic-like (components 1 and 2), one tyrosine-like (components 3) and one fulvic-like (components 4) component. Both components 3 and 4 showed fluorescence quenching with increasing Cd concentrations, while components 1 and 2 had no evident change in fluorescence intensity. The logK₃ and logK₄ values ranged from 4.41 to 5.29 and 4.71 to 5.54, respectively, with most logK values of component 3 for Cd binding being smaller than that of component 4, thus, indicating that the fulvic acid substances exhibited stronger and more stable interactions with Cd than protein-like components.

Modelling chronic toxicokinetics and toxicodynamics of copper in mussels considering ionoregulatory homeostasis and oxidative stress

Chronic toxicity of copper (Cu) at sublethal levels is associated with ionoregulatory disturbance and oxidative stress. These factors were considered in a toxicokinetic-toxicodynamic model in the present study. The ionoregulatory disturbance was evaluated by the activity of the Na⁺/K⁺-ATPase enzyme (NKA), while oxidative stress was presented by lipid peroxidation (LPO) and glutathione-S-transferase (GST) activity. NKA activity was related to the binding of Cu²⁺ and Na ⁺ to NKA. LPO and GST activity were linked with the simulated concentration of unbound Cu. The model was calibrated using previously reported data and empirical data generated when zebra mussels were exposed to Cu. The model clearly demonstrated that Cu might inhibit NKA activity by reducing the number of functional pump sites and the limited Cu-bound NKA turnover rate. An ordinary differential equation was used to describe the relationship between the simulated concentration of unbound Cu and LPO/GST activity. Although this method could not explain the fluctuations in these biomarkers during the experiment, the measurements were within the confidence interval of estimations. Model simulation consistently shows non-significant differences in LPO and GST activity at two exposure levels, similar to the empirical observation.

Dynamic interaction processes of rare earth metal mixtures in terrestrial organisms interpreted by toxicokinetic and toxicodynamic model

Despite the progress in explanation of mixture toxicity of rare earth elements (REEs), a large knowledge gap still exists in interpreting their mixed effects from a dynamic perspective. Here, we investigated the effects of La-Ce mixtures in Enchytraeus crypticus at different exposure times. The single and mixture toxicity of La and Ce increased with time, as reflected by the reduced LC50/MT50 values. With concentration addition as the reference model, the interactions between La and Ce were quantified by MIXTOX modelling tool, showing a time-dependent pattern with antagonistic effect after 1 and 2 d but additive effects afterwards. The dynamic accumulation and toxicity of La/Ce in organisms exposed to REE mixtures was fitted using a process-based toxicokinetic and toxicodynamic (TK-TD) model to unravel how the elements interacted. Generally, the estimated uptake, elimination, and damage rate constants of La/Ce declined with increasing level of each other, suggesting inhibited uptake and subsequently reduced toxicity of La/Ce due to competition effect. The interplay of La and Ce in TK and TD processes seemed responsible for the observed antagonism. Our study showed that mixture toxicity and interaction of REEs are time-dependent processes and application of TK-TD model may provide more insight into this dynamic effect.