Effect-Directed Analysis of Toxicants in Sediment with Combined Passive Dosing and in Vivo Toxicity Testing

Identifying and prioritizing organic toxicants in treated flowback and produced water from shale gas exploitation sites using an integrative effect-directed analysis and nontarget screening method

The use of hydraulic fracturing in shale gas exploitation has generated substantial amount of flowback and produced water (FPW), and ecological risk of these highly complex chemical mixtures has raised worldwide concern. Herein, an integrative effect-directed analysis (EDA) and nontarget screening (NTS) workflow was developed to identify and prioritize main toxicants in the treated FPW (T-FPW). The workflow included sample extraction and fractionation, zebrafish embryo toxicity tests, target and nontarget chemical analyses, and toxicity prioritization and confirmation using toxicological priority index (ToxPi). Results showed that less hydrophobic compounds (log Kow < 3.7) which were used in fracturing fluid and their degradation products were the potentially high-risk toxicants in T-FPW. Thirty-nine target compounds identified in toxic fraction explained 4.82% of the mortality. Additional 584 nontarget contaminants were annotated by NTS. Risk prioritization was achieved for 470 identified contaminants with ecotoxicity data available using a ToxPi method. Six nontarget toxicants were identified with higher ecological risks than all target contaminants, and their presence in FPW were confirmed using reference standards. A principal component analysis of NTS features revealed that EDA fractionation reduced mixture complexity and focused toxicant screening, which significantly improved NTS efficiency, highlighting advantages of integrative EDA and NTS for mixture risk assessment.

Progress, applications, and challenges in high-throughput effect-directed analysis for toxicity driver identification - is it time for HT-EDA?

The rapid increase in the production and global use of chemicals and their mixtures has raised concerns about their potential impact on human and environmental health. With advances in analytical techniques, in particular, high-resolution mass spectrometry (HRMS), thousands of compounds and transformation products with potential adverse effects can now be detected in environmental samples. However, identifying and prioritizing the toxicity drivers among these compounds remain a significant challenge. Effect-directed analysis (EDA) emerged as an important tool to address this challenge, combining biotesting, sample fractionation, and chemical analysis to unravel toxicity drivers in complex mixtures. Traditional EDA workflows are labor-intensive and time-consuming, hindering large-scale applications. The concept of high-throughput (HT) EDA has recently gained traction as a means of accelerating these workflows. Key features of HT-EDA include the combination of microfractionation and downscaled bioassays, automation of sample preparation and biotesting, and efficient data processing workflows supported by novel computational tools. In addition to microplate-based fractionation, high-performance thin-layer chromatography (HPTLC) offers an interesting alternative to HPLC in HT-EDA. This review provides an updated perspective on the state-of-the-art in HT-EDA, and novel methods/tools that can be incorporated into HT-EDA workflows. It also discusses recent studies on HT-EDA, HT bioassays, and computational prioritization tools, along with considerations regarding HPTLC. By identifying current gaps in HT-EDA and proposing new approaches to overcome them, this review aims to bring HT-EDA a step closer to monitoring applications.

Deep Learning Bridged Bioactivity, Structure, and GC-HRMS-Readable Evidence to Decipher Nontarget Toxicants in Sediments

Identifying causative toxicants in mixtures is critical, but this task is challenging when mixtures contain multiple chemical classes. Effect-based methods are used to complement chemical analyses to identify toxicants, yet conventional bioassays typically rely on an apical and/or single endpoint, providing limited diagnostic potential to guide chemical prioritization. We proposed an event-driven taxonomy framework for mixture risk assessment that relied on high-throughput screening bioassays and toxicant identification integrated by deep learning. In this work, the framework was evaluated using chemical mixtures in sediments eliciting aryl-hydrocarbon receptor activation and oxidative stress response. Mixture prediction using target analysis explained <10% of observed sediment bioactivity. To identify additional contaminants, two deep learning models were developed to predict fingerprints of a pool of bioactive substances (event driver fingerprint, EDFP) and convert these candidates to MS-readable information (event driver ion, EDION) for nontarget analysis. Two libraries with 121 and 118 fingerprints were established, and 247 bioactive compounds were identified at confidence level 2 or 3 in sediment extract using GC-qToF-MS. Among them, 12 toxicants were analytically confirmed using reference standards. Collectively, we present a "bioactivity-signature-toxicant" strategy to deconvolute mixtures and to connect patchy data sets and guide nontarget analysis for diverse chemicals that elicit the same bioactivity.

High-Efficiency Effect-Directed Analysis Leveraging Five High Level Advancements: A Critical Review

Organic contaminants are ubiquitous in the environment, with mounting evidence unequivocally connecting them to aquatic toxicity, illness, and increased mortality, underscoring their substantial impacts on ecological security and environmental health. The intricate composition of sample mixtures and uncertain physicochemical features of potential toxic substances pose challenges to identify key toxicants in environmental samples. Effect-directed analysis (EDA), establishing a connection between key toxicants found in environmental samples and associated hazards, enables the identification of toxicants that can streamline research efforts and inform management action. Nevertheless, the advancement of EDA is constrained by the following factors: inadequate extraction and fractionation of environmental samples, limited bioassay endpoints and unknown linkage to higher order impacts, limited coverage of chemical analysis (i.e., high-resolution mass spectrometry, HRMS), and lacking effective linkage between bioassays and chemical analysis. This review proposes five key advancements to enhance the efficiency of EDA in addressing these challenges: (1) multiple adsorbents for comprehensive coverage of chemical extraction, (2) high-resolution microfractionation and multidimensional fractionation for refined fractionation, (3) robust in vivo/vitro bioassays and omics, (4) high-performance configurations for HRMS analysis, and (5) chemical-, data-, and knowledge-driven approaches for streamlined toxicant identification and validation. We envision that future EDA will integrate big data and artificial intelligence based on the development of quantitative omics, cutting-edge multidimensional microfractionation, and ultraperformance MS to identify environmental hazard factors, serving for broader environmental governance.

Rapid screening of acetylcholinesterase active contaminants in water: A solid phase microextraction-based ligand fishing approach

Effect-directed analysis (EDA) has been increasingly used for screening toxic contaminants in the environment, but conventional EDA procedures are often time-consuming and labor-extensive. This challenges the use of EDA for toxicant identification in the scenarios when quick answers are demanded. Herein, a solid phase microextraction ligand fishing (SPME-LF) strategy has been proposed as a rapid EDA approach for identifying acetylcholinesterase (AChE) active compounds in water. The feasibility of ligand fishing techniques for screening AChE active chemicals from environmental mixtures was first verified by a membrane separation method. Then, SPME fibers were prepared through self-assembly of boronic acid groups with AChE via co-bonding and applied for SPME-LF. As AChE coated SPME fibers selectively enriched AChE-active compounds from water, comparing sorbing compounds by the SPME fibers with and without AChE coating can quickly distinguish AChE toxicants in mixtures. Compared with conventional EDA, SPME-LF does not require repeating sample separations and bioassays, endowing SPME-LF with the merits of low-cost, labor-saving, and user-friendly. It is believed that cost-efficient and easy-to-use SPME-LF strategy can potentially be a rapid EDA method for screening receptor-specific toxicants in aquatic environment, especially applicable in time-sensitive screening.

Determining buffering capacity of polydimethylsiloxane-based passive dosing for hydrophobic organic compounds in large-volume bioassays

Polydimethylsiloxane (PDMS) is the most widely used material for passive dosing. However, the ability of PDMS to maintain constant water concentrations of chemicals in large-volume bioassays was insufficiently investigated. In this study, we proposed a kinetic-based method to determine the buffering capacity of PDMS for maintaining constant water concentrations of hydrophobic organic contaminants (HOCs) in large-volume bioassays. A good correlation between log Kow and PDMS-water partitioning coefficients (log KPW) was observed for HOCs with log Kow values ranging from 3.30 to 7.42. For low-molecular-weight HOCs, volatile loss was identified as the primary cause of unstable water concentrations in passive dosing systems. Slow desorption from PDMS resulted in a reduction of water concentrations for high-molecular-weight HOCs. The volume ratio of PDMS to water (RV) was the key factor controlling buffering capacity. As such, buffering capacity was defined as the minimum RV required to maintain 90% of the initial water concentration and was determined to be 0.0076-0.032 for six representative HOCs. Finally, passive dosing with an RV of 0.014 was validated to effectively maintain water concentrations of phenanthrene in 2-L and 96-h toxicity tests with adult mosquitofish. By determining buffering capacity of PDMS, this study recommended specific RV values for cost-efficient implementation of passive dosing approaches in aquatic toxicology, particularly in large-volume bioassays.

Effect-Directed Analysis Combined with Nontarget Screening to Identify Unmonitored Toxic Substances in the Environment

Effect-directed analysis (EDA) combined with nontarget screening (NTS) has established a valuable tool for the identification of unmonitored toxic substances in environmental samples. It consists of three main steps: (1) highly potent fraction identification, (2) toxicant candidate selection, and (3) major toxicant identification. Here, we discuss the methodology, current status, limitations, and future challenges of EDA combined with NTS. This method has been applied successfully to various environmental samples, such as sediments, wastewater treatment plant effluents, and biota. We present several case studies and highlight key results. EDA has undergone significant technological advancements in the past 20 years, with the establishment of its key components: target chemical analysis, bioassays, fractionation, NTS, and data processing. However, it has not been incorporated widely into environmental monitoring programs. We provide suggestions for the application of EDA combined with NTS in environmental monitoring programs and management, with the identification of further research needs.

Noncovalent Tagging for Identifying Unknown Contaminants of Specific Bioactivity in Environmental Water

Identifying contaminants of specific bioactivities from complicated environmental matrices remains costly and time-consuming, as it requires us to not only resolve their structures but also determine their bioactivities. Herein, a novel noncovalent tagging method is integrated in mass spectrometry for identifying unknown contaminants that target dopamine (DA) receptors. Via proteolysis of bovine serum albumin, a stereoselective hexapeptide (ACFAVE) is selected for noncovalently tagging the contaminants that possess the stereostructural characteristics of binding to DA receptors. The tagged contaminants can be readily distinguished from the coexisting species for subsequent structural analysis based on the tagging-induced shifts of the mass-to-charge ratios. Thus, both bioactivity evaluation and structure analysis are accomplished via mass spectrometry. By using this method, 1,3-diphenylguanidine (DPG), a widely used additive in rubber and plastics, is successfully identified out of 2495 features detected in the Pearl River water, with its concentration determined as only 9.8 μg L-1. Furthermore, DPG is confirmed as a potential disrupter to the DA receptors via a simulated docking experiment, which has not been reported before. The present noncovalent tagging method provides a cost-effective and time-efficient way of identifying bioactive molecules in complicated matrices. And proteolysis of proteins is promising for developing more taggants with other desired stereoselectivities in the future.

Identification and Prioritization of Environmental Organic Pollutants: From an Analytical and Toxicological Perspective

Exposure to environmental organic pollutants has triggered significant ecological impacts and adverse health outcomes, which have been received substantial and increasing attention. The contribution of unidentified chemical components is considered as the most significant knowledge gap in understanding the combined effects of pollutant mixtures. To address this issue, remarkable analytical breakthroughs have recently been made. In this review, the basic principles on recognition of environmental organic pollutants are overviewed. Complementary analytical methodologies (i.e., quantitative structure-activity relationship prediction, mass spectrometric nontarget screening, and effect-directed analysis) and experimental platforms are briefly described. The stages of technique development and/or essential parts of the analytical workflow for each of the methodologies are then reviewed. Finally, plausible technique paths and applications of the future nontarget screening methods, interdisciplinary techniques for achieving toxicant identification, and burgeoning strategies on risk assessment of chemical cocktails are discussed.

Methodology for Effect-Based Identification of Bioconcentratable Endocrine Disrupting Chemicals (EDCs) in Water: Establishment, Validation, and Application

Since the wide occurrence of endocrine disrupting chemicals (EDCs) in water is associated with various adverse effects in aquatic organisms, it is urgent to identify key bioconcentratable EDCs. Currently, bioconcentration is generally ignored during the identification of key EDCs. Thus, a methodology for effect-based identification of bioconcentratable EDCs was established in Microcosm, validated in the field, and applied to typical surface water in Taihu Lake. In Microcosm, an inverted U-shaped relationship between logBCFs and logKows was observed for typical EDCs, with medium hydrophobic EDCs (3 ≤ logKow ≤ 7) exhibiting the greatest bioconcentration potentials. On this basis, enrichment methods for bioconcentratable EDCs were established using POM and LDPE, which better fitted the bioconcentration characteristics and enabled the enrichment of 71 ± 8% and 69 ± 6% bioconcentratable compounds. The enrichment methods were validated in the field, where LDPE exhibited a more significant correlation with the bioconcentration characteristics than POM, with mean correlation coefficients of 0.36 and 0.15, respectively, which was selected for further application. By application of the new methodology in Taihu Lake, 7 EDCs were prioritized from 79 identified EDCs as key bioconcentratable EDCs on consideration of their great abundance, bioconcentration potentials, and anti-androgenic potencies. The established methodology could support the evaluation and identification of bioconcentratable contaminants.

Bioassay-based identification and removal of target and suspect toxicants in municipal wastewater: Impacts of chemical properties and transformation

Municipal wastewater contains numerous chemicals and transformation products with highly diverse physiochemical properties and intrinsic toxicity; thus, it is imperative but challenging to identify major toxicants. Herein, toxicity identification evaluation (TIE) was applied to identify major toxicants in a typical municipal wastewater treatment plant (WWTP). Impacts of chemical properties on the removal of contaminants and toxicity at individual treatment stages were also examined. The WWTP influent caused 100% death of Daphnia magna and zebrafish embryos, and toxicity characterization suggested that organics, metals, and volatiles all contributed to the toxicity. Toxicity identification based on 189 target and approximately one-thousand suspect chemicals showed that toxicity contributions of organic contaminants, metals, and ammonia to D. magna were 77%, 4%, and 19%, respectively. Galaxolide, pyrene, phenanthrene, benzo[a]anthracene, fluoranthene, octinoxate, silver, and ammonia were identified as potential toxicants. Comparatively, the detected transformation products elicited lower toxicity than their respective parent contaminants. In contrast, the analyzed contaminants showed negligible contributions to the toxicity of zebrafish embryos. Removal efficiencies of these toxicants in WWTP were highly related to their hydrophobicity. Diverse transformation and removal efficiencies of contaminants in WWTPs may influence the chemical compositions in effluent and ultimately the risk to aquatic organisms in the receiving waterways.

Effect-directed analysis for revealing aryl hydrocarbon receptor agonists in sediment samples from an electronic waste recycling town in China

Exposure to electronic and electrical waste (e-waste) has been related to a few adverse health effects. In this study, sediment samples from an e-waste recycling town in China were collected, and aryl hydrocarbon receptor (AhR) agonists in the samples were identified using an effect-directed analysis (EDA) strategy. The CBG2.8D cell line reporter gene bioassay was used as a toxicity test, while suspect screening against chemical databases was performed for potential AhR agonist identification where both gas chromatography- and liquid chromatography-high resolution mass spectrometry analyses were run. When the original sample extract showed high AhR-mediated activity, sample fractionation was performed, and fractions exhibiting high bioactivity were chemically analyzed again to reveal the corresponding AhR agonists. In total, 23 AhR agonists were identified, including 14 commonly known ones and 9 new ones. Benzo [k]fluoranthene and 6-nitrochrysene were the dominant AhR agonists, covering 16-71% and 2.7-12%, respectively, of the AhR activation effects measured in the parent extracts. The newly identified AhR-active chemicals combined explained 0.13-0.20% of the parent extracts' effects, with 7,12-dimethylbenz [a]anthracene and 8,9,11-trimethylbenz [a]anthracene being the major contributors. A diagnostic isomer ratio analysis of polycyclic aromatic hydrocarbons suggested that the major source of AhR agonists identified in these e-waste related sediment samples were probably petroleum product combustion and biomass combustion. In the future, for a more comprehensive AhR agonist investigation, in-house chemical synthesis and purification, and, when necessary, a secondary sample fractionation, would be beneficial.

Effect-directed analysis of toxic organics in PM2.5 exposure to the cellular bioassays in vitro: Application in Shanxi of China

To optimize the effect-directed analysis (EDA) approach to identify the fine particulate matter (PM_2.5) bound organic toxicants, Jinzhong city, in the Shanxi Province of China, was selected as the object of our study. First, PM_2.5 samples were collected and their organic extracts were separated out in 9 fractions (F1-F9) using reversed-phase high performance liquid chromatography after purification using gel permeation chromatography. Second, the toxicity effects of each fraction were measured by human bronchial epithelial cells (BEAS-2B) in vitro. And toxicity effects included antioxidant stress (ROS, LDH, and CAT) and an inflammatory response (IL-6, IL-1β, and TNF-α). The results showed that the scores of the toxicity effects on multiple lines of evidence were the highest in the F3 and F4 fractions compared with those of the control. Subsequently, the main poisons, o-cymene, p-cymene, benzene, ethylbenzene, xylene, and styrene, were identified using GC×GC-TOF/MS. Finally, to confirm the above possible candidates, (1) the levels of o-cymene, p-cymene and BTEXS in daily PM_2.5 were measured using GC-MS in November 2020, and the rates of detection of these pollutants were 100% in PM_2.5. Among them, o-cymene and p-cymene were first reported as the key toxic substances of PM_2.5, and their average concentration values were 0.16 ± 0.11 and 0.18 ± 0.15 ng‧m^-3, respectively. (2) the toxicity of p-cymene may be no less than that of other benzene derivatives according to their LC50 in Daphnia magna. (3) based on canonical correlation analysis, the exposure to p-cymene, benzene, and styrene in PM_2.5 was most likely associated with the toxicity effects (CAT, IL-6, and TNF-α), which in turn caused the observed toxicity. In conclusion, p-cymene, benzene, and styrene were found to be the key toxic organics in PM_2.5 for cells in vitro. EDA technology avoids the limitations of chemical analysis and uncertainty of the biological testing and adds new toxicants to the control list of PM_2.5, contributing to this study field. However, the application of EDA to PM_2.5 still faces challenges such as the selection of biological effects, loss of toxicity with the separation process, influence of the dosing method, and identification of the unknown effects of pollutants.

Point or non-point source: Toxicity evaluation using m-POCIS and zebrafish embryos in municipal sewage treatment plants and urban waterways

Municipal sewage treatment plants (STPs) have been regarded as an important source of organic contaminants in aquatic environment. To assess the impact of STPs on occurrence and toxicity of STP-associated contaminants in receiving waterways, a novel passive sampler modified from polar organic chemical integrative sampler (m-POCIS) was deployed at the inlet and outlet of a STP and several upstream and downstream sites along a river receiving STP effluent in Guangzhou, China. Eighty-seven contaminants were analyzed in m-POCIS extracts, along with toxicity evaluation using zebrafish embryos. Polycyclic musks were the predominant contaminants in both STP and urban waterways, and antibiotics and current-use pesticides (e.g., neonicotinoids, fiproles) were also ubiquitous. The m-POCIS extracts from downstream sites caused significant deformity in embryos, yet the toxicity could not be explained by the measured contaminants, implying the presence of nontarget stressors. Sewage treatment process substantially reduced embryo deformity, chemical oxygen demand, and contamination levels of some contaminants; however, concentrations of neonicotinoids and fiproles increased after STP treatment, possibly due to the release of chemicals from perturbed sludge. Source identification showed that most of the contaminants found in urban waterways were originated from nonpoint runoff, while cosmetics factories and hospitals were likely point sources for musks and antibiotics, respectively. Although the observed embryo toxicity could not be well explained by target contaminants, the present study showed a promising future of using passive samplers to evaluate chemical occurrence and aquatic toxicity concurrently. Zebrafish embryo toxicity significantly decreased after sewage treatment, but higher toxicity was observed for downstream samples, demonstrating that urban runoff may produce detrimental effects to aquatic life, particularly in rainy season. These results highlight the relevance of monitoring nonpoint source pollution along with boosting municipal sewage treatment infrastructure.

Toxicity assessment of typical polycyclic aromatic hydrocarbons to Daphnia magna and Hyalella azteca in water-only and sediment-water exposure systems

While the equilibrium partitioning (EqP) method has been demonstrated to effectively predict the adverse effects of nonionic organic chemicals in sediment on benthic organisms by sediment toxicity tests, only a limited number of studies have been performed both in water-only and whole-sediment toxicity tests using the same species and verified the validity of EqP-based toxicity assessment. To further examine the validity of the EqP method for application in a wide range of hydrophobicity, we conducted sorption/desorption experiments and both water-only and sediment toxicity tests using a popular aquatic crustacean species, Daphnia magna (48 h), and benthic species Hyalella azteca (96 h) for six typical polycyclic aromatic hydrocarbons (PAHs) with three to five rings and an amine derivative: anthracene, phenanthrene, fluoranthene, pyrene, benzo[a]pyrene, dibenzo[a,h]anthracene, and 1-aminopyrene. The linear sorption coefficient was determined and ranged from 2.7 × 10² (phenanthrene) and 1.2 × 10⁴ L/kg (benzo[a]pyrene) highly depending on the hydrophobicity while the aqueous concentrations were stable after 24 h in the desorption test. As result of acute toxicity tests in the water-only exposure system, anthracene and dibenz[a,h]anthracene were found to be nontoxic to both species, while median effect/lethal concentrations (EC₅₀/LC₅₀) were determined as ranging from 0.66 (benzo[a]pyrene) to 330 μg/L (phenanthrene), and from 11 (1-aminopyrene) to 180 μg/L (phenanthrene) for D. magna and H. azteca, respectively. Among these compounds, three PAHs with three, four, and five rings each, and 1-aminopyrene were subjected to sediment-water toxicity tests. In the sediment-water tests, the LC₅₀ of phenanthrene and pyrene was three to six times higher than that of the water-only tests for H. azteca while the EC₅₀ was 1.1 to 2.0 times higher for D. magna. In contrast, the EC₅₀/LC₅₀ of benzo[a]pyrene (BaP) in the sediment-water toxicity test was more than 5 times higher than that in the water-only test for both H. azteca and D. magna. The EC₅₀/LC₅₀ values of 1-aminopyrene were similar in both the sediment-water and the water-only toxicity tests, ranging narrowly from 21 to 28 μg/L and 8.8 to 11 μg/L for D. magna and H. azteca, respectively. The EC₅₀/LC₅₀ based on the body residue (ER₅₀/LR₅₀) was investigated for two of the representative PAHs, pyrene, and BaP. The ER₅₀/LR₅₀ of pyrene in both species was 2.3 and 11 times higher in the water-only toxicity test for D. magna and H. azteca, respectively, while those of BaP in the sediment-water toxicity test were not calculated for the sediment-water toxicity tests, and the highest body concentration in the sediment-water tests was lower than the ER₅₀/LR₅₀ in the water-only toxicity test. Although the experimental results were comparable with the predicted sediment toxicity values based on the EqP method for the selected PAHs in this study, there is a risk of phenanthrene and pyrene being slightly underestimated (1.4-1.9 fold for phenanthrene and 3.7-6.1 fold for pyrene) by the EqP method for H. azteca. These results reaffirm that the bioavailability of poorly water-soluble chemicals is important for sediment toxicity and that the exposure pathway should be further investigated to avoid under- and overestimation via the EqP method.

Can the hydrophobic organic contaminants in the filtrate passing through 0.45 μm filter membranes reflect the water quality?

In the traditional water quality assessment, the concentration of total dissolved hydrophobic organic compounds (HOCs) passing through 0.45 μm filter membranes is usually used to evaluate the influence of HOCs on water quality. However, the bioavailability of dissolved organic matter (DOM)-associated and particle-associated HOCs is not considered. In the present work, pyrene, fulvic acid, and natural suspended particles (SPS) were used to simulate natural water (raw water). The immobilization and pyrene content in the tissues of D. magna caused by total pyrene in the raw water and those caused by freely dissolved pyrene with the concentration equal to the total dissolved pyrene in the filtrate of raw water were compared to determine whether the total dissolved pyrene concentration can reflect the water quality. The results indicated that when the DOM concentration was 5 mg C L^-1 and the SPS concentration was higher than 0.2-0.4 g L^-1, the bioavailability of pyrene was underestimated by the traditional water quality assessment because of the SPS-associated pyrene, and it was underestimated by 23.6-63.9% when SPS concentration was higher than 0.6 g L^-1 due to the neglection of SPS-associated pyrene. Furthermore, the threshold value of SPS concentration was related to the SPS size and composition, and the effects of SPS and DOM on water quality were influenced by the concentration, size, and composition of SPS as well as the molecular weight of DOM. This study suggests that the traditional water quality assessment should be improved by comprehensively considering concentrations and characteristics of SPS and DOM.

The distributions of three fungicides in vegetables and their potential health risks in Zhejiang, China: A 3-year study (2015-2017)

Fungicides have been extensively used around the world in agriculture due to their effectiveness of sterilization. Recent evidences have shown that fungicides would exert a negative effect on gut microbiota and result in gut microbiota dysbiosis and metabolism disorder on non-target organisms and even humans. However, research on residues and potential health risks of fungicides in daily consumed vegetables has received less attention compared to insecticides. In this study, we studied three widely applied fungicides, procymidone, dimethomorph, and azoxystrobin, in China. We collected 551 samples of 10 different vegetables in 11 cities from Zhejiang province during 2015-2017. Three fungicides were analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The average apparent recoveries of three fungicides ranged from 84.2% to 110% with the relative standard deviations lower than 10%. The LOD values of procymidone, dimethomorph, and azoxystrobin was 2, 0.09, and 1 μg/kg, respectively. The levels of procymidone, dimethomorph, and azoxystrobin in those vegetables ranged from ND-875, ND-238, and ND-76 μg/kg, respectively. The highest mean concentrations of procymidone, dimethomorph, and azoxystrobin were found in eggplant (68 μg/kg), spinach (16.4 μg/kg), and kidney bean (4 μg/kg), respectively. Tomato (62.6% of samples), eggplant (44.3% of samples), and cucumber (41.6% of samples) were most frequently detected with fungicides. Solanaceous fruit vegetables have the highest detection rate than other vegetables, and fungicides were most frequently detected in winter. The mean concentrations of three fungicides in different vegetables were all below the maximum residue limits for the national food safety standards of China, and the health risks resulting from consuming those vegetables in adults and children were all within the safe ranges. The data provided here clarify the distributions of fungicides in commonly consumed vegetables and their potential health risks.

Identifying bioaccessible suspect toxicants in sediment using adverse outcome pathway directed analysis

Chemical mixtures are a common occurrence in contaminated sediment and determining causal relationship between sediment contamination and adverse outcomes is challenging. The bioavailability and choice of bioassay endpoints played important roles in elucidating causality. As such, bioaccessibility-based XAD extraction and adverse outcome pathway (AOP) guided bioassays were incorporated into an effect-directed analysis to more effectively determine sediment causality. XAD extracts of sediments from urban waterways in Guangzhou, China were examined using cell viability bioassays with four human tumor cells from lung, liver, breast, and bone marrow. Pronounced effects to SH-SY5Y cells were noted, thus neurotoxicity was subsequently focused in the AOP-guided bioassays. Intracellular calcium influx, mitochondrial membrane potential inhibition, reactive oxygen species generation, and cell viability were utilized as evidence for neurotoxicity AOP-guided analysis. Suspect toxicants were identified in active fractions using GC-MS. Toxicity confirmation was performed by evaluating toxicity contributions of the candidates to the pathway. Cypermethrin, bisphenol A, galaxolide, tonalide, and versalide were found as the major stressors across key events of the studied pathway. Moreover, good correlations among key events validated the feasibility of method to predict in vivo response, suggesting that considering bioavailability and AOP improved environmental relevance for toxicant identification in a complex mixture.

Accelerated Passive Dosing of Hydrophobic Complex Mixtures-Controlling the Level and Composition in Aquatic Tests

Petroleum products and essential oils are complex mixtures of hydrophobic and volatile chemicals and are categorized as substances of unknown or variable composition, complex reaction products, or biological materials (UVCBs). In aquatic testing and research of such mixtures, it is challenging to establish initial concentrations without the addition of cosolvents, to maintain constant concentrations during the test, and to keep a constant mixture composition in dilution series and throughout test duration. Passive dosing was here designed to meet these challenges by maximizing the surface area (A_donor/V_medium = 3.8 cm²/mL) and volume (V_donor/V_medium > 0.1 L/L) of the passive dosing donor in order to ensure rapid mass transfer and avoid donor depletion for all mixture constituents. Cracked gas oil, cedarwood Virginia oil, and lavender oil served as model mixtures. This study advances the field by (i) showing accelerated passive dosing kinetics for 68 cracked gas oil constituents with typical equilibration times of 5-10 min and for 21 cederwood Virginia oil constituents with typical equilibration times < 1 h, (ii) demonstrating how to control mixture concentration and composition in aquatic tests, and (iii) discussing the fundamental differences between solvent spiking, water-accommodated fractions, and passive dosing.

Using the HPTLC-bioluminescence bacteria assay for the determination of acute toxicities in marine sediments and its eligibility as a monitoring assessment tool

For an integrated ecological risk assessment of marine sediment contamination, the determination of target-compound concentrations by e.g. mass spectrometric methods is not sufficient to explain sediment toxicity. Due to the presence of a multitude of environmental contaminants in this complex matrix causing a mixed toxicity, the identification and assessment of main toxicants is a challenge. One approach in identifying main toxicants is the application of effect-directed analysis (EDA). In this study, an EDA approach was developed using high performance thin layer chromatography (HPTLC) coupled to bioluminescence bacteria detection with Aliivibrio fischeri for the determination of marine sediments acute toxicity. In a first attempt, the HPTLC separation was optimized with a fast, two-step gradient to separate main hydrophobic organic contaminant (HOC) classes found in marine sediments. An easy-to-use evaluation script for the resulting bioluminescence inhibition images was programed using R. The developed method was applied to sediment extracts of two different sample sets: (i) Fourteen marine sediment samples from the German Bight representing a wide range of contaminant loads and sediment properties and (ii) sediment samples from a core representing temporal trends of contamination. Results from the HPTLC-bioluminescence bacteria assay were compared to HOC concentrations determined by GC-MS/MS. A correlation of the determined inhibition Γ-values for the PAH inhibition zone to PAH concentrations showed a very good agreement (R² = 0.91). The results of this study were used to evaluate the suitability of the EDA approach to be used as an assessment tool for marine sediments.

Understanding the bioavailability of pyrethroids in the aquatic environment using chemical approaches

Pyrethroids are a class of commonly used insecticides and are ubiquitous in the aquatic environment in various regions. Aquatic toxicity of pyrethroids was often overestimated when using conventional bulk chemical concentrations because of their strong hydrophobicity. Over the last two decades, bioavailability has been recognized and applied to refine the assessment of ecotoxicological effects of pyrethroids. This review focuses on recent advances in the bioavailability of pyrethroids, specifically in the aquatic environment. We summarize the development of passive sampling and Tenax extraction methods for assessing the bioavailability of pyrethroids. Factors affecting the bioavailability of pyrethroids, including physicochemical properties of pyrethroids, and quality and quantity of organic matter, were overviewed. Various applications of bioavailability on the assessment of bioaccumulation and acute toxicity of pyrethroids were also discussed. The final section of this review highlights future directions of research, including development of standardized protocols for measurement of bioavailability, establishment of bioavailability-based toxicity benchmarks and water/sediment quality criteria, and incorporation of bioavailability into future risk assessment and management actions.

Expanding phytoremediation to the realms of known and unknown organic chemicals of concern

Recent advancements in analytical chemistry and data analyses via high-resolution mass spectrometry (HRMS) are evolving scientific understanding of the potential totality of organic chemical exposure and pollutant risk. This review addresses the importance of HRMS approaches, namely suspect screening and nontarget chemical analyses, to the realm of phytoremediation. These analytical approaches are not without caveats and constraints, but they provide an opportunity to understand in greater totality how plant-based technologies contribute, mitigate, and reduce organic chemical exposure across scales of experimental and system-level studies. These analytical tools can enlighten the complexity and efficacy of plant-contaminant system design and expand our understanding of biogenic and anthropogenic chemicals at work in phytoremediation systems. Advances in data analytics from biological sciences, such as metabolomics, are crucial to HRMS analysis. This review provides an overview of targeted, suspect screening, and nontarget HRMS approaches, summarizes the expanding knowledge of regulated and unregulated organic chemicals in the environment, addresses requisite HRMS instrumentation, analysis cost, uncertainty, and data processing techniques, and offers potential bridges of HRMS analyses to phytoremediation research and application.

Occurrence and risk of neonicotinoid insecticides in surface water in a rapidly developing region: Application of polar organic chemical integrative samplers

Extensive use of neonicotinoid insecticides (NNIs) worldwide calls for further knowledge on their environmental occurrence and risk. The present study highlighted the need for more research on long-term exposure of NNIs in aquatic environment, which is important but remains elusive. Time weighted average concentrations of five commonly used NNIs in urban waterways of Guangzhou, China were measured using newly developed polar organic chemical integrative samplers (POCIS). Acetamiprid (from 18.8 ± 1.9 to 157 ± 31 ng/L; mean ± standard deviation), clothianidin (from 14.8 ± 3.7 to 47.6 ± 10.0 ng/L) and imidacloprid (from 32.9 ± 11.6 to 249 ± 19 ng/L) were detectable in all samples. Thiamethoxam was found at 71.4% of the 21 sampling sites (from not detected to 52.4 ± 9.4 ng/L), while thiacloprid was not detected at any site. Vegetable planting and sewage effluent were the main sources of NNIs in surface water in Guangzhou. Probabilistic environmental exposure distributions were subsequently constructed using the measured concentrations and the exceedances of predicted environmental concentrations of NNI to ecological thresholds were assessed. In Guangzhou, 63.5%, 16.2%, 87.8% and 17.2% of acetamiprid, clothianidin, imidacloprid and thiamethoxam, respectively, exceeded an interim chronic threshold of 35 ng/L for NNIs. Further risk assessment and control measures for the use of NNIs are advocated for protecting the integrity of aquatic ecosystems.

Identifying Organic Toxicants in Sediment Using Effect-Directed Analysis: A Combination of Bioaccessibility-Based Extraction and High-Throughput Midge Toxicity Testing

Toxicity identification evaluation (TIE) and effect-directed analysis (EDA) were integrated to diagnose toxicity drivers in a complex system, such as sediment. In TIE manipulation, XAD resin was utilized as an amending agent for characterizing organic toxicants, which also facilitate a large-volume bioaccessibility-based extraction for EDA purposes. Both raw sediments in TIE and extract fractions in EDA were tested with Chironomus dilutus for toxicity using whole-sediment testing and a high-throughput microplate assay. This allowed for a direct link between whole-sediment TIE and EDA, which strongly strengthened the characterization and identification of toxicants. Sediments amended with XAD resin, as part of the TIE, significantly reduced midge mortality compared with unamended sediments, suggesting that organics were one class of main toxicants. On the basis of bioaccessible concentrations in sediment measured by XAD extraction, a group of previously unidentified contaminants, synthetic polycyclic musks (versalide, tonalide, and galaxolide), were found to explain 32-73% of the observed toxicity in test sediments. Meanwhile, three pyrethroids contributed to an additional 17-35% of toxicity. Surprisingly, the toxicity contribution of musks and pyrethroids reached 58-442 and 56-1625%, respectively, based on total sediment concentrations measured by exhaustive extraction. This suggested that total sediment concentrations significantly overestimated toxicity and that bioavailability should be considered in toxicity identification. Identifying nontarget toxicants sheds a light on application of the integrated TIE and EDA method in defining causality in a complex environment.

Emerging investigator series: effect-based characterization of mixtures of environmental pollutants in diverse sediments

This study investigated whether cell-based bioassays were suitable to characterize profiles of mixture effects of hydrophobic pollutants in multiple sediments covering remote Arctic and tropical sites to highly populated sites in Europe and Australia. The total contamination was determined after total solvent extraction and the bioavailable contamination after silicone-based passive equilibrium sampling. In addition to cytotoxicity, we observed specific responses in cell-based reporter gene bioassays: activation of metabolic enzymes (arylhydrocarbon receptor: AhR, peroxisome proliferator activated receptor gamma: PPARγ) and adaptive stress responses (oxidative stress response: AREc32). No mixture effects were found for effects on the estrogen, androgen, progesterone and glucocorticoid receptors, or they were masked by cytotoxicity. The bioanalytical equivalent concentrations (BEQ) spanned several orders of magnitude for each bioassay. The bioavailable BEQs (passive equilibrium sampling) typically were 10-100 times and up to 420 times lower than the total BEQ (solvent extraction) for the AhR and AREc32 assays, indicating that the readily desorbing fraction of the bioactive chemicals was substantially lower than the fraction bound strongly to the sediment sorptive phases. Contrarily, the bioavailable BEQ in the PPARγ assay was within a factor of five of the total BEQ. We identified several hotspots of contamination in Europe and established background contamination levels in the Arctic and Australia.

Sublethal and chronic effects of reclaimed water on aquatic organisms. Looking for relationships between physico-chemical characterisation and toxic effects

The use of reclaimed water for irrigation and aquaculture purposes is generally considered a reliable alternative for sustainable water management in regions with water scarcity. Many organic compounds, generally called compounds of emerging concern (CECs), have been detected in reclaimed water, which implies continuous exposure for aquatic organisms. To date no quality criteria have been proposed for this group of compounds. This work aims to assess the acute, sublethal and chronic effects of reclaimed water using two representative organisms of the aquatic compartment; the green alga Chlorella vulgaris and the microcrustacean Daphnia magna. The study comprises the 72 h-algal growth inhibition test, the D. magna feeding bioassay and the D. magna reproduction test. The results highlighted, for the selected characterised compounds, no differences in the concentrations between the different tertiary WWTP treatments, except for the particular case of carbamazepine. Considering seasonality, no differences were observed between the two different sample collection campaigns. The sublethal and chronic effects observed for these samples could not be explained by the lower concentrations found in the chemical characterisation. However, in the majority of cases, dilution of raw reclaimed water reduced the toxic effects of these samples. Several interactions among compounds can affect the mixture's toxicity. Canonical correlation analyses (CCA) were included to explore the potential relationships between the physico-chemical characterisation of reclaimed water and effects on aquatic organisms. The results corroborated the toxic effect of some pharmaceuticals, in particular beta-blockers and antibiotics, on the growth and yield of green algae, as well as inhibition of daphnia reproduction. Thus the CCA methods could help to elucidate the potential relationships between the physico-chemical characterisation and toxic effects by considering all the potential interactions.

Dissolved organic matter affects both bioconcentration kinetics and steady-state concentrations of polycyclic aromatic hydrocarbons in zebrafish (Danio rerio)

Dissolved organic matter (DOM) is ubiquitous in natural aquatic ecosystems. The association of hydrophobic organic compounds (HOCs), such as polycyclic aromatic hydrocarbons (PAHs), with DOM may have a large impact on HOC fractions in water and their bioconcentration in fish. However, the effects of DOM on HOC bioconcentration in fish are not well understood, especially whether DOM will affect the bioconcentration steady-state concentrations of HOCs in fish is still confusing. Thus, this study investigated the effects of three DOM including gallic acid (GA), tannic acid (TA), and humic acid (HA) with molecular weights ranging from 170 Da to about 10 kDa at different concentrations (1, 5, and 15 mgOC L^-1) on the bioconcentration of PAHs including phenanthrene, anthracene, fluoranthene, and pyrene in zebrafish (Danio rerio), with the PAH freely dissolved concentrations maintained constant by passive dosing systems. The results revealed that the presence of DOM generally increased the bioconcentration steady-state concentrations of the PAHs in zebrafish (C_b-ss), with the increase ranging from 28.1% to 204.0%, and the increase of C_b-ss promoted by TA with middle molecular weight (1700 Da) was the highest among the studied DOM. Moreover, the C_b-ss increased with the concentrations of GA with low molecular weight and TA with middle molecular weight in water, whereas decreased with increasing concentrations of HA with high molecular weight. The uptake rate constants of the PAHs in zebrafish with DOM were higher than that without DOM. Ingestion of DOM and direct accumulation of PAHs associated with DOM might be primary influencing mechanisms of DOM on the C_b-ss, and whether the facilitated diffusive mass transfer by DOM will affect the C_b-ss needs to be further studied. This study suggested that DOM-associated HOCs should be considered in future HOC risk assessment in addition to the freely dissolved HOCs.

Bioavailability of Pyrene Associated with Different Types of Protein Compounds: Direct Evidence for Its Uptake by Daphnia magna

The protein-like dissolved organic matter (DOM) is ubiquitous in aquatic environments. However, the bioavailability of protein-like DOM-associated hydrophobic organic compounds (HOCs) is not well-understood, and in particular, the direct evidence of their uptake by organisms is scarce. In the present work, tryptone (2000 Da), bovine serum albumin (BSA; 66 000 Da), and phycocyanin (120 000 Da) were chosen as model protein-like DOM, which were labeled by commercial fluorescein (cy5) to investigate the uptake mechanisms of protein compound-associated pyrene (a typical HOC) by Daphnia magna. The pyrene concentration in the tissues except the gut and immobilization of D. magna were detected to calculate the bioavailable fraction of protein compound-associated pyrene when the freely dissolved pyrene concentration was controlled through passive dosing devices. The results demonstrated that the tryptone could permeate cellular membrane and directly enter the tissues of  D. magna from the exposure solutions, whereas BSA and phycocyanin might indirectly enter the tissues from the gut. A part of pyrene associated with protein compounds was bioavailable to D. magna; the order of their bioavailable fractions was trypone (54.6-58.1%) > phycocyanin (21.6-32.8%) > BSA (17.7-26.8%). The difference was principally related to the uptake mechanisms of pyrene associated with different types of protein. This work suggests that the protein compound-associated HOCs should be considered to evaluate the bioavailability and eco-environmental hazard of HOCs in natural waters.

Diagnosis of complex mixture toxicity in sediments: Application of toxicity identification evaluation (TIE) and effect-directed analysis (EDA)

Determining causality of sediment toxicity is of great importance in aquatic risk assessment, but there are tremendous challenges due to joint toxicity of trace pollutants in complex sediment matrices. Two approaches, namely toxicity identification evaluation (TIE) and effect-directed analysis (EDA) have been developed. Conventional sediment TIEs take the advantage of environmental relevance by using whole organism bioassays; however, they suffer from lack of effective methods for specifically identifying major contributors as it typically only evaluates contaminant class rather than specific contaminants. Alternatively, EDA is a powerful tool in identifying causes of sediment toxicity with sophisticated fractionation and chemical analysis of targeted and non-targeted non-polar organic toxicants, but it is not always environmentally relevant due to the use of in-vitro bioassays and exhaustive solvent extraction. An integrated TIE and EDA method would provide an environmentally relevant and toxicant specific approach to effectively determine causality of sediment toxicity by combining the merits of the two methods. Bioavailability-based extraction and dosing techniques are recommended to be incorporated into the integrated method to improve the accuracy of toxicity diagnosis. Besides considering bioavailability in the integrated TIE and EDA approach, the premise of adverse outcome pathways should also be considered. Generally speaking, both TIE and EDA have focused on adverse effects at cellular and organism levels. The addition of trait-based approaches in screening multiple toxicological endpoints helps to extend effects on cellular and organism levels to population level, and provides a better understanding of potential impacts to the community and ecosystem. The outcome pathway underlies the critical role of determining causality in interpreting impacts of complex mixtures to benthic community and aquatic ecosystem.

Improving the accuracy of effect-directed analysis: the role of bioavailability

Aquatic ecosystems have been suffering from contamination by multiple stressors. Traditional chemical-based risk assessment usually fails to explain the toxicity contributions from contaminants that are not regularly monitored or that have an unknown identity. Diagnosing the causes of noted adverse outcomes in the environment is of great importance in ecological risk assessment and in this regard effect-directed analysis (EDA) has been designed to fulfill this purpose. The EDA approach is now increasingly used in aquatic risk assessment owing to its specialty in achieving effect-directed nontarget analysis; however, a lack of environmental relevance makes conventional EDA less favorable. In particular, ignoring the bioavailability in EDA may cause a biased and even erroneous identification of causative toxicants in a mixture. Taking bioavailability into consideration is therefore of great importance to improve the accuracy of EDA diagnosis. The present article reviews the current status and applications of EDA practices that incorporate bioavailability. The use of biological samples is the most obvious way to include bioavailability into EDA applications, but its development is limited due to the small sample size and lack of evidence for metabolizable compounds. Bioavailability/bioaccessibility-based extraction (bioaccessibility-directed and partitioning-based extraction) and passive-dosing techniques are recommended to be used to integrate bioavailability into EDA diagnosis in abiotic samples. Lastly, the future perspectives of expanding and standardizing the use of biological samples and bioavailability-based techniques in EDA are discussed.

Extended Virtual Screening Strategies To Link Antiandrogenic Activities and Detected Organic Contaminants in Soils

A tiered screening strategy based on extensive virtual fractionation and elucidation was developed to simplify identification of toxicants in complex environments. In tier1-virtual fractionation, multivariate analysis (MVA) was set up as an alternative of physical fractionation. In tier2-virtual structure elucidation, in-house quantitative structure-retention relationship (QSRR) models and toxicity simulation methods were developed to simplify nontarget identification. The efficiency of the tiered virtual strategy was tentatively verified by soil samples from a chemical park contaminated by antiandrogenic substances. Eight out of 18 sites were detected as antiandrogenic, while none of them exhibited androgenic agonist potencies. Sixty-seven peaks were selected for further identification by MVA, among which over 90% were verified in androgenic fractions in traditional effect-directed analysis (EDA). With 579 tentative structures generated by in silico fragmentation, 74% were elucidated by QSRR and 65% were elucidated by in silico toxicity prediction. All prior peaks were identified at different confidence levels with over 40% of the identified peaks above confidence level 2b, which has been increased over 40% with less than half of the time spent compared to traditional EDA. Such a combination of tiered virtual screening methods provides more efficient and rapid identifications of key toxicants at contaminated sites.