Use of effect-directed analysis for the identification of organic toxicants in surface flow constructed wetland sediments

High-efficiency effect-directed analysis (EDA) advancing toxicant identification in aquatic environments: Latest progress and application status

Facing the great threats to ecosystems and human health posed by the continuous release of chemicals into aquatic environments, effect-directed analysis (EDA) has emerged as a powerful tool for identifying causative toxicants. However, traditional EDA shows problems of low-coverage, labor-intensive and low-efficiency. Currently, a number of high-efficiency techniques have been integrated into EDA to improve toxicant identification. In this review, the latest progress and current limitations of high-efficiency EDA, comprising high-coverage effect evaluation, high-resolution fractionation, high-coverage chemical analysis, high-automation causative peak extraction and high-efficiency structure elucidation, are summarized. Specifically, high-resolution fractionation, high-automation data processing algorithms and in silico structure elucidation techniques have been well developed to enhance EDA. While high-coverage effect evaluation and chemical analysis should be further emphasized, especially omics tools and data-independent mass acquisition. For the application status in aquatic environments, high-efficiency EDA is widely applied in surface water and wastewater. Estrogenic, androgenic and aryl hydrocarbon receptor-mediated activities are the most concerning, with causative toxicants showing the typical structural features of steroids and benzenoids. A better understanding of the latest progress and application status of EDA would be beneficial to further advance in the field and greatly support aquatic environment monitoring.

Novel polar AhR-active chemicals detected in sediments of an industrial area using effect-directed analysis based on in vitro bioassays with full-scan high resolution mass spectrometric screening

Studies investigating aryl hydrocarbon receptor (AhR)-active compounds in the environment typically focus on non- and mid-polar substances, such as PAHs; while, information on polar AhR agonists remains limited. Here, we identified polar AhR agonists in sediments collected from the inland creeks of an industrialized area (Lake Sihwa, Korea) using effect-directed analysis combined with full-scan screening analysis (FSA; using LC-QTOFMS). Strong AhR-mediated potencies were observed for the polar and latter fractions of RP-HPLC (F3.5-F3.8) from sediment organic extracts in the H4IIE-luc in vitro bioassays. FSA was performed on the corresponding fractions. Twenty-eight tentative AhR agonists were chosen using a five-step process. Toxicological confirmation using bioassay revealed that canrenone, rutaecarpine, ciprofloxacin, mepanipyrim, genistein, protopine, hydrocortisone, and medroxyprogesterone were significantly active. The relative potencies of these AhR-active compounds compared to that of benzo[a]pyrene ranged from 0.00002 to 2.0. Potency balance analysis showed that polar AhR agonists explained, on average, ~6% of total AhR-mediated potencies in samples. Some novel polar AhR agonists also exhibited endocrine-disrupting potentials capable of binding to estrogen and glucocorticoid receptors, as identified by QSAR modeling. In conclusion, the focused studies on distributions, sources, fate, and ecotoxicological effects of novel polar AhR agonists in the environment are necessary.

Major AhR-active chemicals in sediments of Lake Sihwa, South Korea: Application of effect-directed analysis combined with full-scan screening analysis

This study utilized effect-directed analysis (EDA) combined with full-scan screening analysis (FSA) to identify aryl hydrocarbon receptor (AhR)-active compounds in sediments of inland creeks flowing into Lake Sihwa, South Korea. The specific objectives were to (i) investigate the major AhR-active fractions of organic extracts of sediments by using H4IIE-luc in vitro bioassay (4 h and 72 h exposures), (ii) quantify known AhR agonists, such as polycyclic aromatic hydrocarbons (PAHs) and styrene oligomers (SOs), (iii) identify unknown AhR agonists by use of gas chromatography-quadrupole time-of-flight mass spectrometry (GC-QTOFMS), and (iv) determine contributions of AhR agonists to total potencies measured by use of the bioassay. FSA was conducted on fractions F2.6 and F2.7 (aromatics with log K_ow 5-7) in extracts of sediment from Siheung Creek (industrial area). Those fractions exhibited significant AhR-mediated potency as well as relatively great concentrations of PAHs and SOs. FSA detected 461 and 449 compounds in F2.6 and F2.7, respectively. Of these, five tentative candidates of AhR agonist were selected based on NIST library matching, aromatic structures and numbers of rings, and available standards. Benz[b]anthracene, 11H-benzo[a]fluorene, and 4,5-methanochrysene exhibited significant AhR-mediated potency in the H4IIE-luc bioassay, and relative potencies of these compounds were determined. Potency balance analysis demonstrated that these three newly identified AhR agonists explained 1.1% to 67% of total induced AhR-mediated potencies of samples, which were particularly great for industrial sediments. Follow-up studies on sources and ecotoxicological effects of these compounds in coastal environments would be required.

Multiple evaluation of the potential toxic effects of sediments and biota collected from an oil-polluted area around Abu Ali Island, Saudi Arabia, Arabian Gulf

After the Gulf War Oil Spill, there have been many investigations about distributions of oil-derived pollutants nearby areas, but lacking in ecotoxicological assessment. We evaluated the potential toxicity of asphalt mats, sediments, and biota (polychaetes, chitons, snapping shrimps, and crabs) by combining two bioassays (H4IIE-luc and Vibrio fischeri) and in situ microbial community (eDNA). Samples were collected from Abu Ali Island, and organic extracts were bioassayed and further fractionated according to the chemical polarity using silica gel column. Great aryl hydrocarbon receptor (AhR)-mediated potencies and inhibition of bioluminescence were mainly found in aromatics (F2) and saturates (F1) fractions of asphalt mat and sediments, respectively, while great toxicological responses in biota samples were found in resins and polar (F3) fraction. We also confirmed that potential toxicities of biota were species-specific; great AhR-mediated potencies were found in polychaetes and great bioluminescence inhibitions were found in crabs. In microbial communities, most genera (up to 90%) were associated with polycyclic aromatic hydrocarbons (PAHs)-degrading bacteria, supporting that PAHs are the primary stressors of the benthic community around Abu Ali Island. The present study provides useful information on the contamination status, risk assessment of environmental matrices and benthic organisms in Abu Ali Island.

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.

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.

Characterization of endocrine disruption potentials of coastal sediments of Taean, Korea employing H295R and MVLN assays-Reconnaissance at 5years after Hebei Spirit oil spill

Endocrine disrupting potentials were assessed for sediment samples collected near Hebei Spirit oil spill (HSOS) site, between December 2007 and January 2012. For comparison, major crude oil (CO) of HSOS, or its weathered form were assessed. Both raw extracts (REs) and their fractionated samples were tested using H295R and MVLNluc bioassays. In H295R cells, REs of crude and weathered oil (WO), and nine of 14 sediments significantly increased E2 levels, which were correlated with the concentrations of PAHs. Steroidogenic disruption potentials of the sediments generally decreased over time. Among silica fractions of all REs, aromatic hydrocarbons (F2) and polar compounds (F3) caused greater E2 levels. While, in MVLN cell bioassay, only three of 14 sediment REs showed estrogen receptor binding potencies, and no temporal trend was observed. In conclusion, oil spill can cause endocrine disruption in the affected ecosystem through steroidogenic alteration for years, and such potencies attenuate over time.

Normal-phase liquid chromatography retention behavior of polycyclic aromatic sulfur heterocycles and alkyl-substituted polycyclic aromatic sulfur heterocycle isomers on an aminopropyl stationary phase

Retention indices for 67 polycyclic aromatic sulfur heterocycles (PASHs) and 80 alkyl-substituted PASHs were determined using normal-phase liquid chromatography (NPLC) on an aminopropyl (NH₂) stationary phase. The retention behavior of PASH on the NH₂ phase is correlated with the number of aromatic carbon atoms and two structural characteristics have a significant influence on their retention: non-planarity (thickness, T) and the position of the sulfur atom in the bay-region of the structure. Correlations between solute retention on the NH₂ phase and T of PASHs were investigated for three cata-condensed (cata-) PASH isomer groups: (a) 13 four-ring molecular mass (MM) 234 Da cata-PASHs, (b) 20 five-ring MM 284 Da cata-PASHs, and (c) 12 six-ring MM 334 Da cata-PASHs. Correlation coefficients ranged from r = -0.49 (MM 234 Da) to r = -0.65 (MM 334 Da), which were significantly lower than structurally similar PAH isomer groups (r = -0.70 to r = -0.99). The NPLC retention behavior of the PASHs are compared to similar results for PAHs.

Normal-phase liquid chromatography retention behavior of polycyclic aromatic hydrocarbon and their methyl-substituted derivatives on an aminopropyl stationary phase

Retention indices for 124 polycyclic aromatic hydrocarbons (PAHs) and 62 methyl-substituted (Me-) PAHs were determined using normal-phase liquid chromatography (NPLC) on a aminopropyl (NH2) stationary phase. PAH retention behavior on the NH2 phase is correlated to the total number of aromatic carbons in the PAH structure. Within an isomer group, non-planar isomers generally elute earlier than planar isomers. MePAHs generally elute slightly later but in the same region as the parent PAHs. Correlations between PAH retention behavior on the NH2 phase and PAH thickness (T) values were investigated to determine the influence of non-planarity for isomeric PAHs with four to seven aromatic rings. Correlation coefficients ranged from r = 0.19 (five-ring peri-condensed molecular mass (MM) 252 Da) to r = -0.99 (five-ring cata-condensed MM 278 Da). In the case of the smaller PAHs (MM ≤ 252 Da), most of the PAHs had a planar structure and provided a low correlation. In the case of larger PAHs (MM ≥ 278 Da), nonplanarity had a significant influence on the retention behavior and good correlation between retention and T was obtained for the MM 278 Da, MM 302 Da, MM 328 Da, and MM 378 Da isomer sets. Graphical abstract NPLC separation of the three-, four-, five-, and six-ring PAH isomers with different number of aromatic carbon atoms and degrees of non-planarity (Thickness, T). The inserted figure plots the number of aromatic carbon atoms vs. the log I value for the 124 parent PAHs.

Qualitative characterization of SRM 1597a coal tar for polycyclic aromatic hydrocarbons and methyl-substituted derivatives via normal-phase liquid chromatography and gas chromatography/mass spectrometry

A normal-phase liquid chromatography (NPLC) fractionation procedure was developed for the characterization of a complex mixture of polycyclic aromatic hydrocarbons (PAHs) from a coal tar sample (Standard Reference Material (SRM) 1597a). Using a semi-preparative aminopropyl (NH₂) LC column, the coal tar sample was separated using NPLC based on the number of aromatic carbons; a total of 14 NPLC fractions were collected. SRM 1597a was analyzed before and after NPLC fractionation by using gas chromatography/mass spectrometry (GC/MS) with a 50% phenyl stationary phase. The NPLC-GC/MS method presented in this study allowed for the identification of 72 PAHs and 56 MePAHs. These identifications were based on the NPLC retention times for authentic reference standards, GC retention times for authentic reference standards, and the predominant molecular ion peak in the mass spectrum. Most noteworthy was the determination of dibenzo[a,l]pyrene, which could not be measured directly by GC/MS because of low concentration and co-elution with dibenzo[j,l]fluoranthene. The NPLC-GC/MS procedure also allowed for the tentative identification of 74 PAHs and 117 MePAHs based on the molecular ion peak only. This study represents the most comprehensive qualitative characterization of SRM 1597a to date. Graphical abstract NPLC-GC/MS analysis for the six-ring MM 302 Da PAH isomers in SRM 1597a.

Assessment of potential biological activities and distributions of endocrine-disrupting chemicals in sediments of the west coast of South Korea

The west coast of Korea has experienced environmental deterioration for more than half a century. In the present study, we specifically aimed to: i) evaluate potential toxicities of contaminants in sediments that cause effects mediated through the aryl hydrocarbon receptor (AhR) and estrogen receptor (ER); ii) determine spatio-temporal distributions of polycyclic aromatic hydrocarbons (PAHs) and alkylphenols (APs); and iii) identify causes of greater potencies of samples. From 2010 to 2014, sediments were collected from 12 major estuarine and coastal regions along the west coast of South Korea. In vitro cell bioassays were performed to determine AhR- and ER-mediated potencies using H4IIE-luc and MVLN cells, respectively. Fifteen PAHs and six APs in sediments were identified by GC/MSD. Results of bioassays generally showed a low-to-moderate degree of contamination, however, greater AhR- and ER-mediated potencies were measured at some locations. Concentrations of PAHs and APs varied among locations, which indicated that sources were independently affected by the surrounding environment (e.g., industrial complex and cities). Results of bioassays were generally well correlated with concentrations of putative causative chemicals. Benzo[k]fluoranthene, dibenz[a,h]anthracene, and benzo[b]fluoranthene were the major AhR agonists, explaining approximately 30% of the bioassay-derived benzo[a]pyrene equivalent concentration (BaP-EQ). Unknown AhR and ER agonists and potential mixture effects remain in question. Overall, the present study provides baseline information on chemical contaminations and potential toxicity of sediments in a fairly wide geographical region of the west coast of South Korea.

Evaluation of estrogenic activity in the Pearl River by using effect-directed analysis

This study investigated estrogenic activity of water, sediment, and fish bile of the Pearl River in southern China by effect-directed analysis based on in vitro yeast screen assay and chemical analysis. Results showed higher estradiol equivalents (EEQ) for surface water in dry season than in wet season. Simple risk assessment suggested that high estrogenic risk would be expected in Shima River and Danshui River receiving discharge of effluents from cities in the region. Fractionation and effect-directed analysis showed that estrogenic activity mainly occurred in relatively polar fractions of surface water. Seven target estrogenic compounds (bisphenol A, 4-nonylphenol, 4-tert-octylphenol, 17α-ethynyl estradiol, estrone, diethylstilbestrol, and 17β-estradiol) only accounted for part of the measured estrogenic activity, with the rest contributions from other potential estrogenic chemicals such as phenols. Findings from this study suggest that fish in the river could be affected by those estrogenic chemicals. Proper measures should be taken to reduce the estrogenic activity in wastewaters before they are discharged into the riverine system in order to protect aquatic organisms.

Effect-directed analysis and mixture effects of AhR-active PAHs in crude oil and coastal sediments contaminated by the Hebei Spirit oil spill

The major AhR-active PAHs were identified in crude oil and oil-contaminated sediments by use of effect-directed analysis. As part of the study, an enhanced potency balance analysis was conducted by establishing the novel relative potency values of (alkyl)-PAHs from the H4IIE-luc bioassay. Silica gel column fractionation of crude oil resulted in greater AhR-mediated potencies in fractions of aromatics (F2) and resins (F3), and such trend was also observed for field collected sediment samples. AhR-mediated potencies of six F2 sub-fractions from HPLC indicated that the majority of F2 responses were attributable to 3-4 ring aromatics. Target PAHs including C4-phenanthrene, C1-chrysene, and C3-chrysene in sediments explained ∼ 18% of the bioassay-derived TCDD-EQs, however, the unknown AhR agonists and potential mixture effects remain in question. Overall, the AhR-potency and antagonistic potential of residual oil in sediment tended to decrease over time, thus monitoring of weathering process would be key for the post management of oil-contaminated sites.

Cytostatic drugs and metabolites in municipal and hospital wastewaters in Spain: filtration, occurrence, and environmental risk

Concerns about cytostatic anticancer drugs in the environment are increasing, mainly due to the lack of knowledge about the fate and impact of these cytotoxic compounds in the water cycle. In this context, the present work investigated the occurrence of 13 cytostatics and 4 metabolites in wastewater samples from various wastewater treatment plants (WWTPs) and from a large hospital from Spain. The target compounds belong to five different classes according to the Anatomical Therapeutic Classification (ATC), namely, alkylating agents, antimetabolites, plant alkaloids and other natural products, cytotoxic antibiotics and related substances, and other antineoplastic agents. Some of them have been classified as carcinogens in humans by the International Agency for Research on Cancer (IARC). These compounds were determined by an automated on line solid-phase extraction-liquid chromatography-tandem mass spectrometry (SPE-LC-MS/MS) method. Results showed the presence of methotrexate (MET), ifosfamide (IF), cyclophosphamide (CP), irinotecan (IRI), doxorubicin (DOX), capecitabine (CAP), tamoxifen (TAM) and the metabolites endoxifen (OH-D-TAM), hydroxytamoxifen (OH-TAM) and hydroxypaclitaxel (OH-PAC) at levels ranging from 2 ng L(-1) (for MET) to 180 ng L(-1) (for TAM). Some of these compounds were found to be efficiently removed after wastewater treatment, e.g. MET, DOX and IRI, whereas other compounds, such as TAM, CP and IF remained largely unaltered. The behaviour of the target compounds during the common filtration step of the water samples was also investigated with the finding that some compounds are strongly adsorbed to nylon filters, while cellulose acetate appears as the best choice for the filter material. The aquatic environmental risk associated to the detected compounds was also assessed. To the best of the authors' knowledge, this is the first report of the presence of the metabolites OH-D-TAM and OH-TAM in the water cycle.

Identification of photosynthesis inhibitors of pelagic marine algae using 96-well plate microfractionation for enhanced throughput in effect-directed analysis

Because of large-scale production and use of an increasing diversity of chemicals in modern society, estuarine and coastal waters may be contaminated with numerous substances. Some of these compounds have the potential to affect microalgae at the base of the pelagic food chain. Therefore, we identified the main chemical stressors that negatively affect the effective photosystem II efficiency (ϕPSII) in marine microalgae of the Dutch estuarine and coastal waters. An enhanced effect-directed analysis (EDA) was carried out by combining reversed-phase ultra performance liquid chromatography fractionation of extracts from passive samplers, followed by effect assessment using the pulse amplitude modulation fluorometry assay and chemical analysis of biologically active fractions using high-resolution mass spectrometry. This study focuses on a novel microfractionation technique using 96-well plates to enhance throughput in EDA, structure elucidation, and the analytical and effect confirmation of the compounds that are identified. Although there are numerous unknown compounds present in estuarine and coastal waters, our EDA study shows that atrazine, diuron, irgarol, isoproturon, terbutryn, and terbutylazine are the main contributors to the observed effect on the ϕPSII of marine microalgae.

Effect-directed analysis to explore the polar bear exposome: identification of thyroid hormone disrupting compounds in plasma

Compounds with transthyretin (TTR)-binding potency in the blood plasma of polar bear cubs were identified with effect-directed analysis (EDA). This approach contributes to the understanding of the thyroid disrupting exposome of polar bears. The selection of these samples for in-depth EDA was based on the difference between the observed TTR-binding potency on the one hand and the calculated potency (based on known concentrations of TTR-binding compounds and their relative potencies) on the other. A library-based identification was applied to the liquid chromatography-time-of-flight-mass spectrometry (LC-ToF-MS) data by screening for matches between compound lists and the LC-ToF-MS data regarding accurate mass and isotope pattern. Then, isotope cluster analysis (ICA) was applied to the LC-ToF-MS data allowing specific screening for halogen isotope patterns. The presence of linear and branched nonylphenol (NP) was observed for the first time in polar bears. Furthermore, the presence of one di- and two monohydroxylated octachlorinated biphenyls (octaCBs) was revealed in the extracts. Linear and branched NP, 4'-OH-CB201 and 4,4'-OH-CB202 could be successfully confirmed with respect to their retention time in the analytical system. In addition, branched NP, mono- and dihydroxylated-octaCBs showed TTR-binding potencies and could explain another 32 ± 2% of the total measured activities in the extracts.

Occurrence of PCDD/PCDFs and PCBs in soil and comparison with CYP1A response in PLHC-1 cell line

The responsiveness of CYP1A (gene transcription and EROD enzyme activity) in the cell line Poeciliopsis lucida hepatoma (PLHC-1) upon exposure to extracts of contaminated soil samples was investigated and compared to levels of PCDD/PCDFs and PCBs including non-ortho obtained by GC/MS analysis. Soil samples A and B were collected in sites A and B. Two fractions, not purified (np) and purified (p), were obtained from each sample and analyzed for PCDD/PCDF and PCB content by GC/MS; in parallel they were tested for 24 h with PLHC-1. CYP1A response was investigated at gene (RT-qPCR) level and as 7-ethoxyresorufin-O-deethylase (EROD) enzyme activity. Chem-TEQs and Bio-TEQs were then calculated. ∑TEQ calculated for PCDD/Fs and PCBs was 0.081 pg/g and 20.32 pg/g for samples A and B, respectively. PLHC-1 showed less up-regulation of cyp1a gene on exposure to the two purified fractions (Ap 2.1-fold and Bp 1.8-fold) than to non-purified fractions (up to 15-fold for Anp and 13-fold for Bnp). EROD was also induced 2.38- and 9.44-fold in the two purified fractions (Ap and Bp) compared to model inducer 2,3,7,8-TCDD, and up to 16.03-fold for non-purified Anp and 33.79-fold for Bnp. The combination of CYP1A response, obtained in a PLHC-1 cell-based bioassay, with contaminant residue analysis provided a better description of the presence and toxicity of dioxin-like compounds in an environmental matrix.